Hi, I am trying to run a triaxial test. I have defined properties as:
# interaction property [NOTE: Young's modulus value is same for both particle and wall]
variable y1 equal 65e9 # Young's Modulus particle-particle [Pa]
variable y2 equal 6.5e9 # Young's Modulus particle-wall [Pa]
variable us_pp equal 0.15 # Coefficient of static friction particle-particle
variable us_pw equal 0.05 # Coefficient of static friction particle-wall
variable ur_pp equal 0.01 # Coefficient of rolling friction particle-particle
variable ur_pw equal 0.01 # Coefficient of rolling friction particle-wall
variable p_ratio equal 0.25 # Poisson's ratio
variable damping equal 1.5e-6 # Damping coefficient
variable cor equal 0.8 # Coefficient of Restitution
# Material properties required for new pair styles
fix p1 all property/global youngsModulus peratomtypepair ${y1} ${y2}
fix p3 all property/global gamman peratomtypepair 2 ${damping} ${damping} ${damping} ${damping}
fix p4 all property/global gammat peratomtypepair 2 0 0 0 0
fix p5 all property/global coefficientFriction peratomtypepair 2 ${us_pp} ${us_pw} ${us_pw} ${us_pw}
fix p6 all property/global coefficientRollingFriction peratomtypepair 2 ${ur_pp} ${ur_pw} ${ur_pw} ${ur_pw}
fix p7 all property/global coefficientRestitution peratomtypepair 2 ${cor} ${cor} ${cor} ${cor}
fix p8 all property/global poissonsRatio peratometype ${p_ratio} ${p_ratio}
But this piece of code is giving error during execution.
ERROR: Fix property/global (id p1): the number of default values must be a multiple of nCols. (/build/liggghts-YO7u74/liggghts-3.8.0+repack1/src/fix_property_global.cpp:129)
I have tried different combinations, all giving the same error. Please let me know if you have any suggestion.
mschramm | Tue, 09/29/2020 - 15:18
Wrong type
You should use
peratomtype
for the Young's modulus.