how to calculate the total volume of particles after deformation?

Hi Qiang,

Like most DEM codes, ESyS-Particle uses a soft contact model. Particles are allowed to overlap and the repulsive force between particles is a function of the overlap distance. The particles themselves remain spherical so the total volume of the particles is always the same regardless of the deformation.

I presume what you wish to calculate is the total volume of the space containing the particles once deformed. This can be obtained via a Runnable during simulations using the LsmMpi.getCurMinMaxPt(..) subroutine. A simple usage example would be:

minPt = Vec3(0,0,0)
maxPt = Vec3(0,0,0)

sim.getCurMinMaxPt(minPt, maxPt)

print minPt, maxPt

This method will only be accurate if the deformed volume of the particle assembly is a rectangular prism.

I hope this helps.

Cheers,

Dion