Hi all,
I am using pair_style gran model hertz tangential history in my input script.
I am getting the following error, when i am running the simulation.
ERROR: Could not locate a fix/property storing value(s) for poissonsRatio as requested by model hertz. (../modify_liggghts.cpp:313)
Can someone please suggest, what can be the reason for this.
Thank you,
Narendra
richti83 | Thu, 11/13/2014 - 14:38
you did not define the
you did not define the material properties for liggghts-granular, have a look in the examples/LIGGGHTS/Tutorials_public/contactModels/ folder
you need at least
fix m1 all property/global youngsModulus peratomtype 1.e7
fix m2 all property/global poissonsRatio peratomtype 0.3
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.2
fix m4 all property/global coefficientFriction peratomtypepair 1 0.3
fix m5 all property/global coefficientRollingFriction peratomtypepair 1 0.3
sometimes the problem is, when the fix-id (fix mN all) is used twice that the last fix overwrites the first fix and liggghts can not find a needed property.
naren | Fri, 11/14/2014 - 06:17
Hi richti83
I have found the error. It was the spelling mistake in my input script.
Thanks,
Narendra
naren | Fri, 11/14/2014 - 05:42
Hi rishti83,
Hi rishti83,
I have defined the material properties, and i used it only once. Can you please look into the input script attached below.
This is one of the example tutorial, i want to run.
units si
atom_style sphere
boundary f f p
newton off
communicate single vel yes
region reg block -0.078 0.078 -0.078 0.078 -0.1524 0.306 units box
create_box 2 reg
#neigh_modify delay 0
fix m1 all property/global youngsModulus peratomtype 2.5e7 2.5e7
fix m2 all property/global poissonRatio peratomtype 0.25 0.25
fix m3 all property/global coefficientRestitution peratomtypepair 2 0.5 0.5 0.5 0.5
fix m4 all property/global coefficientFriction peratomtypepair 2 0.5 0.5 0.5 0.5
fix m5 all property/global coefficientRollingFriction peratomtypepair 2 0.1 0.1 0.1 0.1
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 1000 radius constant 0.00125
fix pts2 all particletemplate/sphere 1 atom_type 2 density constant 1000 radius constant 0.00100
fix pdd all particledistribution/discrete 12345 2 pts1 0.5 pts2 0.5
fix ins_mesh all mesh/surface file meshes/factory.stl type 1 scale 0.001
fix ins all insert/stream seed 123481 distributiontemplate pdd nparticles 450000 particlerate 900000 overlapcheck yes &
vel constant 0. 0. -3.0 insertion_face ins_mesh extrude_length 0.02
fix cad1 all mesh/surface file meshes/outer_cylinder.stl type 1 scale 0.001 curvature 1e-5
fix cad2 all mesh/surface file meshes/inner_cylinder.stl type 1 scale 0.001 curvature 1e-5
fix cad3 all mesh/surface file meshes/funnel.stl type 1 scale 0.001 curvature 1e-5
fix cad4 all mesh/surface file meshes/plate.stl type 1 scale 0.001 curvature 1e-5
fix geometry all wall/gran model hertz tangential history mesh n_meshes 4 meshes cad1 cad2 cad3 cad4 #rolling_friction cdt
pair_style gran model hertz tangential history #rolling_friction cdt
pair_coeff * *
fix integrate all nve/sphere
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0
timestep 0.00000625
thermo_style custom step atoms ke cpu
thermo 8000
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes
fix timecheck all check/timestep/gran 1 0.01 0.01
run 1
unfix timecheck
#run 1
dump dmp all custom 16000 post/dump*.tardos id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
#dump dumpstress all mesh/gran/VTK 200 post/dump*.vtk stress wear cad
run 160000 upto
unfix geometry
fix geometry all wall/gran model hertz tangential history mesh n_meshes 3 meshes cad1 cad2 cad3 rolling_friction cdt
fix movecad all move/mesh mesh cad2 rotate origin 0. 0. 0. axis 0. 0. 1. period 1.
run 4800000