Thus in the Hertzian case Kn and Kt can be set to values that corresponds to properties of the material being modeled. This is also true in the Hookean case, except that a spring constant must be chosen that is appropriate for the absolute size of particles in the model. diff --git a/doc/pair_gran.txt b/doc/pair_gran.txt index e15a6aca5f..1027bb77f1 100644 --- a/doc/pair_gran.txt +++ b/doc/pair_gran.txt @@ -117,7 +117,10 @@ constant with units of force/distance. In the Hertzian case, Kn is like a non-linear spring constant with units of force/area or pressure, and as shown in the "(Zhang)"_#Zhang paper, Kn = 4G / (3(1-nu)) where nu = the Poisson ratio, G = shear modulus = E / -(2(1+nu)), and E = Young's modulus. Similarly, Kt = 8G / (2-nu). +(2(1+nu)), and E = Young's modulus. Similarly, Kt = 4G / (2-nu). +(NOTE: in an earlier version of the manual, we incorrectly stated that +Kt = 8G / (2-nu).) + Thus in the Hertzian case Kn and Kt can be set to values that corresponds to properties of the material being modeled. This is also true in the Hookean case, except that a spring constant must be chosen