Hello Users,

I am attempting to set up a peridynamics tensile model. I would like to be able to calculate the modulus at various strain rates. The strain is easy to calculate but I would like to get some input on the best way to calculate the stress. The tensile strain is along the y-axis.

Would I want to use “compute property/atom” to get the forces on the atoms in my “clamp” region, sum those forces in the y-direction and divide by the total area to get stress?

Or would I use “compute stress/atom” in the clamp region and sum those stresses for the total stress?

Or could I simply use the negative pressure in the y-direction?

For clarity, to conduct the simulation, I have a tensile bar aligned in the y-axis. The first 3 particle layers on the top and bottom of the bar are being pulled using “fix move” with a constant velocity.

Thanks for the help,

Abby