Stress computation with peridynamics


I’m trying to compute stresses in peridynamic simulations.
I would like to know if this command (stress/atom) is a good way to do it :

compute 1 middle_atoms stress/atom
compute 2 middle_atoms reduce sum c_1[1]
variable stress equal c_2/(0.94*vol)

This command is ok for molecular dynamics but I’m not sure it is well suited for peridynamics.


It computes the same thing no matter what model you
are running. It’s up to you to figure out if that is a good

measure of what you want for PD.