LIGGGHTS (Version LIGGGHTS-PUBLIC 3.8.0, compiled 2020-03-23-07:35:01 by root, git commit unknown) ###################################Header for General commands####################################################### atom_style granular #Simulation of particles boundary f f f #fixed boundaries -> particles will be deleted if leaving the simulation box units si communicate single vel yes #default newton off #default ###################################System variables################################################################## #Definition of boundaries variable xmin equal -0.25 variable xmax equal 0.25 variable ymin equal -0.4 variable ymax equal 0.4 variable zmin equal 0 variable zmax equal 2.3 #Definition of the timestep variable dt equal 0.5e-4 #timestep = 0.0001 second; Each iteration step represents 0.0001 seconds. ###################################Specific variables for current simulation######################################### variable natoms equal 2 #1 -> particle #2->loputa, razvod ####variable for material properties#### ####Young Modulus#### variable youngmodulus1 equal 1e8 #N/mm² variable youngmodulus2 equal 1e8 #N/mm² ####Poission ratio#### variable poission1 equal 0.3 variable poission2 equal 0.3 ####variable for contact properties#### ####coefficient of restitution#### variable CoR11 equal 0.6 variable CoR12 equal 0.6 variable CoR21 equal 0.6 variable CoR22 equal 0.6 ####sliding friction coefficient#### variable sf11 equal 0.3 variable sf12 equal 0.36 variable sf21 equal 0.36 variable sf22 equal 0.0 ####rolling friction coefficient#### variable rf11 equal 0.8 variable rf12 equal 0.8 variable rf21 equal 0.8 variable rf22 equal 0.8 ####variable for particle#### #Number of particle radius variable nradii equal 2 variable radius1 equal 0.0030 #m variable radius2 equal 0.0030 #m variable frac1 equal 0.2 #20% variable frac2 equal 0.8 #80% variable density equal 7800 #kg/m³ ####filling parameters#### variable filltime equal 3 #seconds variable fillmass equal 80 #kg variable fillmassrate equal ${fillmass}/${filltime} #kg/s variable fillmassrate equal 80/${filltime} variable fillmassrate equal 80/3 variable fillsteps equal ${filltime}/${dt} #Transform time to iteration steps variable fillsteps equal 3/${dt} variable fillsteps equal 3/5e-05 ####settle time#### variable settletime equal 2 #second variable settlesteps equal ${settletime}/${dt} #Transform time to iteration steps variable settlesteps equal 2/${dt} variable settlesteps equal 2/5e-05 ####open parameter of the loputa#### variable openvel equal -0.4 #m/s variable opentime equal 0.5 #second variable opensteps equal ${opentime}/${dt} #Transform time to iteration steps variable opensteps equal 0.5/${dt} variable opensteps equal 0.5/5e-05 ####discharge time##### variable dischargetime equal 12 #seconds variable dischargesteps equal ${dischargetime}/${dt} #Transform time to iteration steps variable dischargesteps equal 12/${dt} variable dischargesteps equal 12/5e-05 ###################################Definition of simulationbox####################################################### region reg block ${xmin} ${xmax} ${ymin} ${ymax} ${zmin} ${zmax} units box region reg block -0.25 ${xmax} ${ymin} ${ymax} ${zmin} ${zmax} units box region reg block -0.25 0.25 ${ymin} ${ymax} ${zmin} ${zmax} units box region reg block -0.25 0.25 -0.4 ${ymax} ${zmin} ${zmax} units box region reg block -0.25 0.25 -0.4 0.4 ${zmin} ${zmax} units box region reg block -0.25 0.25 -0.4 0.4 0 ${zmax} units box region reg block -0.25 0.25 -0.4 0.4 0 2.3 units box create_box 2 reg Created orthogonal box = (-0.25 -0.4 0) to (0.25 0.4 2.3) 1 by 2 by 4 MPI processor grid neighbor 0.004 bin #default neigh_modify delay 0 #default #################################Definition of the contact models#################################################### pair_style gran model hertz tangential history rolling_friction epsd2 #contact model WARNING: The contact model you specified is not located in any whitelist. Because of this the model will be run in an unoptimized version (increasing runtime by up to 20%). In order to optimize this model you have the following options: (i) Run the genAutoExamplesWhitelist.sh script in your LIGGGHTS(R) source folder to automatically parse the input script (ii) Add the model combination by hand to your style_contact_model_user.whitelist that can be found in your LIGGGHTS(R) source folder If you perform one of the steps above LIGGGHTS(R) needs to be recompiled to generate the optimized code. (/build/liggghts-YO7u74/liggghts-3.8.0+repack1/src/contact_models.h:364) pair_coeff * * #default timestep ${dt} timestep 5e-05 fix integrator all nve/sphere #default fix gravi all gravity 9.81 vector 0.0 0.0 -1.0 #gravity of 9.81 m/s² in negative z direction ###################################Definition of Material properties################################################# fix m1 all property/global youngsModulus peratomtype ${youngmodulus1} ${youngmodulus2} fix m1 all property/global youngsModulus peratomtype 100000000 ${youngmodulus2} fix m1 all property/global youngsModulus peratomtype 100000000 100000000 fix m2 all property/global poissonsRatio peratomtype ${poission1} ${poission2} fix m2 all property/global poissonsRatio peratomtype 0.3 ${poission2} fix m2 all property/global poissonsRatio peratomtype 0.3 0.3 fix m3 all property/global coefficientRestitution peratomtypepair ${natoms} 0.6 0.6 0.6 0.6 fix m3 all property/global coefficientRestitution peratomtypepair 2 0.6 0.6 0.6 0.6 fix m4 all property/global coefficientFriction peratomtypepair ${natoms} ${sf11} ${sf12} ${sf21} ${sf22} fix m4 all property/global coefficientFriction peratomtypepair 2 ${sf11} ${sf12} ${sf21} ${sf22} fix m4 all property/global coefficientFriction peratomtypepair 2 0.3 ${sf12} ${sf21} ${sf22} fix m4 all property/global coefficientFriction peratomtypepair 2 0.3 0.36 ${sf21} ${sf22} fix m4 all property/global coefficientFriction peratomtypepair 2 0.3 0.36 0.36 ${sf22} fix m4 all property/global coefficientFriction peratomtypepair 2 0.3 0.36 0.36 0 fix m5 all property/global coefficientRollingFriction peratomtypepair ${natoms} ${rf11} ${rf12} ${rf21} ${rf22} fix m5 all property/global coefficientRollingFriction peratomtypepair 2 ${rf11} ${rf12} ${rf21} ${rf22} fix m5 all property/global coefficientRollingFriction peratomtypepair 2 0.8 ${rf12} ${rf21} ${rf22} fix m5 all property/global coefficientRollingFriction peratomtypepair 2 0.8 0.8 ${rf21} ${rf22} fix m5 all property/global coefficientRollingFriction peratomtypepair 2 0.8 0.8 0.8 ${rf22} fix m5 all property/global coefficientRollingFriction peratomtypepair 2 0.8 0.8 0.8 0.8 ###################################Generation and Loading of the Geometry .stl####################################### fix razvod all mesh/surface file Razvod.stl type 2 scale 0.001 #load mesh from STL file. Type 2 for geometry. Scale down to transform mm to meters fix loputa all mesh/surface file Loputa.stl type 2 scale 0.001 fix MesureMesh all mesh/surface/planar file mesureface.stl type 1 scale 0.001 #load mesh for mesurment fix walls all wall/gran model hertz tangential history rolling_friction epsd2 mesh n_meshes 2 meshes razvod loputa WARNING: The contact model you specified is not located in any whitelist. Because of this the model will be run in an unoptimized version (increasing runtime by up to 20%). In order to optimize this model you have the following options: (i) Run the genAutoExamplesWhitelist.sh script in your LIGGGHTS(R) source folder to automatically parse the input script (ii) Add the model combination by hand to your style_contact_model_user.whitelist that can be found in your LIGGGHTS(R) source folder If you perform one of the steps above LIGGGHTS(R) needs to be recompiled to generate the optimized code. (/build/liggghts-YO7u74/liggghts-3.8.0+repack1/src/contact_models.h:364) ###################################Generation and Insertion of the particles######################################### fix pts1 all particletemplate/sphere 10487 atom_type 1 density constant ${density} radius constant ${radius1} fix pts1 all particletemplate/sphere 10487 atom_type 1 density constant 7800 radius constant ${radius1} fix pts1 all particletemplate/sphere 10487 atom_type 1 density constant 7800 radius constant 0.003 fix pts2 all particletemplate/sphere 11887 atom_type 1 density constant ${density} radius constant ${radius2} fix pts2 all particletemplate/sphere 11887 atom_type 1 density constant 7800 radius constant ${radius2} fix pts2 all particletemplate/sphere 11887 atom_type 1 density constant 7800 radius constant 0.003 fix pdd1 all particledistribution/discrete 32452867 ${nradii} pts1 ${frac1} pts2 ${frac2} fix pdd1 all particledistribution/discrete 32452867 2 pts1 ${frac1} pts2 ${frac2} fix pdd1 all particledistribution/discrete 32452867 2 pts1 0.2 pts2 ${frac2} fix pdd1 all particledistribution/discrete 32452867 2 pts1 0.2 pts2 0.8 fix ins_mesh all mesh/surface/planar file NInsertionsface.stl type 1 scale 0.001 fix ins all insert/stream seed 86028157 distributiontemplate pdd1 mass ${fillmass} massrate ${fillmassrate} overlapcheck yes all_in yes vel constant 0 0 -0.75 insertion_face ins_mesh extrude_length 0.25 fix ins all insert/stream seed 86028157 distributiontemplate pdd1 mass 80 massrate ${fillmassrate} overlapcheck yes all_in yes vel constant 0 0 -0.75 insertion_face ins_mesh extrude_length 0.25 fix ins all insert/stream seed 86028157 distributiontemplate pdd1 mass 80 massrate 26.6666666666667 overlapcheck yes all_in yes vel constant 0 0 -0.75 insertion_face ins_mesh extrude_length 0.25 ###################################Dumping of the data for post-processing to visualize############################ shell mkdir post #Definition of the dumptime variable dumptime equal 0.05 # Every 0.05 s 1 image variable dumpstep equal ${dumptime}/${dt} #Transform to iteration steps variable dumpstep equal 0.05/${dt} variable dumpstep equal 0.05/5e-05 dump dmpparticle all custom/vtk ${dumpstep} post/particles_*.vtk id radius mass vx vy vz dump dmpparticle all custom/vtk 1000 post/particles_*.vtk id radius mass vx vy vz dump dmprazvod all mesh/stl ${dumpstep} post/Razvod*.stl razvod dump dmprazvod all mesh/stl 1000 post/Razvod*.stl razvod dump dmploputa all mesh/stl ${dumpstep} post/Loputa*.stl loputa dump dmploputa all mesh/stl 1000 post/Loputa*.stl loputa ###################################Load Mesh and active mesh for the mesurment ############################ fix MesureMesh all mesh/surface/planar file mesureface.stl type 2 scale 0.001 #load mesh for mesurment fix mesurePipe all massflow/mesh mesh MesureMesh vec_side 0. 0. -1 #Activate mesurment mesh ###################################SAVE data of mesurePipe to variable ################################## variable MassTotalPipe equal f_mesurePipe[1] #field 1 ->sum total mass that has been passed [kg] variable NoPTotalPipe equal f_mesurePipe[2] #field 2 ->sum total NoP that has been passed [particles] variable MassflowRatePipe equal f_mesurePipe[3] #field 3 ->actual mass Flow rate [kg/s] variable NoPflowRatePipe equal f_mesurePipe[4] #field 4 ->actual particle flow rate [particle/s] ###################################Write Data of Mesurment############################################### variable outputtime equals 0.05 #s ERROR: Illegal variable command (/build/liggghts-YO7u74/liggghts-3.8.0+repack1/src/variable.cpp:418)