Hello everybody!
I was wondering if the order of the material properties fix and the contact model in the LIGGGHTS script matters? In the doc it is not written that the material-fix has to be specified before the pair style command, but I wanted to make sure I can change the script like this and the result will be the same:
=== 1 ===
fix m1 all property/global youngsModulus peratomtype 5.e6
fix m2 all property/global poissonsRatio peratomtype 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.3
fix m4 all property/global coefficientFriction peratomtypepair 1 0.5
pair_style gran model hertz tangential history
pair_coeff * *
========
=== 2 ===
pair_style gran model hertz tangential history
pair_coeff * *
fix m1 all property/global youngsModulus peratomtype 5.e6
fix m2 all property/global poissonsRatio peratomtype 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.3
fix m4 all property/global coefficientFriction peratomtypepair 1 0.5
========
Best regards,
Harald
ckloss | Thu, 09/24/2015 - 09:58
Hi Harald,
Hi Harald,
the material properties must be specified then the first "run" command is issued, the results of the two snippets you showed will be identical (however, I can not assure that there might be different round off-errors, but this should probably not be the case)
If LIGGGHTS requires something to be specified beforehand and it's not there, it will exit with an error which will tell you the source of the problem
Christoph
hunger | Thu, 09/24/2015 - 17:08
Thanks for your answer,
Thanks for your answer, Christoph!!