Could not locate a fix/property storing value(s)

Submitted by naren on Thu, 11/13/2014 - 10:42

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's picture

richti83 | Thu, 11/13/2014 - 14:38

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.

I'm not an associate of DCS GmbH and not a core developer of LIGGGHTS®
ResearchGate | Contact

naren | Fri, 11/14/2014 - 06:17

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,

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