Stress computation with peridynamics

LAMMPS LAMMPS Mailing List Mirror
1

Hello,

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.

Hugo

2

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.

Steve