how to compute Elastic potential energy due to particle and wall contact

Submitted by Adithya Ramgopal on Sat, 12/15/2018 - 16:59

Hello,
Can anyone tell me how to get the elastic potential energy due to the particle-wall contact. pe/atom or pe gives only gravitational potential energy. (RESTRICTION: This compute does not include the potential energy due to the overlap of granular particles.).
There is an option 'computeElasticPotential' as a keyword for the hertz contact which saves the elastic potential energy, but I don't know how to use the saved value.
My wall is not a primitive wall it is a mesh.
Can anyone guide me through the procedure?
Thanks in advance.

Forums:

arnom's picture

arnom | Mon, 12/17/2018 - 10:15

There is a fix to compute the elastic energy of particle-particle and particle-wall contacts that uses the computeElasticPotential keyword. However, it is only part of LIGGGHT premium at the moment.

DCS team member & LIGGGHTS(R) core developer

Adithya Ramgopal | Wed, 12/19/2018 - 12:20

Thanks for the reply arnom.
I have one more doubt in fix multisphere/break command.
Can I give overlap (contact overlap between two particles) as a trigger_ name to break the multisphere. As overlap is not a peratom property , can I use this as a property to initiate the breaking of multisphere.
Thanks in advance.