I am trying to run this script file but it is showing that velocity projection on face normal is less than 1e-3.
Script File
#Simple chute wear test
atom_style granular
atom_modify map array
boundary f f f
newton off
communicate single vel yes
units si
region reg block 0. 0.3 -0.05 0.05 0.0 0.2 units box # 30cm x 10cm x 20cm
create_box 1 reg
neighbor 0.003 bin
neigh_modify delay 0
#Material properties required for new pair styles
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
#New pair style
pair_style gran model hertz tangential history #Hertzian without cohesion
pair_coeff * *
timestep 0.00001
fix gravi all gravity 9.81 chute 25
#the chute
fix cad all mesh/surface file meshes/trial.stl type 1
fix inface all mesh/surface file meshes/trial2.stl type 1
fix granwalls all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad
#distributions for insertion
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2500 radius constant 0.002
fix pdd1 all particledistribution/discrete 1 1 pts1 1.0
#region and insertion
group nve_group region reg
region bc block 0 0.02 -0.02 0.02 0.05 0.15 units box
#particle insertion
fix ins nve_group insert/stream seed 5330 distributiontemplate pdd1 &
nparticles 6000 massrate 0.1 insert_every 1000 overlapcheck yes all_in no vel constant 0.0 0.0 -1.0 &
insertion_face inface
#apply nve integration to all particles that are inserted as single particles
fix integr nve_group nve/sphere
#output settings, include total thermal energy
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 vol
thermo 1000
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes
#insert the first particles so that dump is not empty
run 1
dump dmp all custom 1000 post/dump*.txt x y z vx vy vz
#dump_modify dmp sort 1
#dump_modify dmp sort 1 element S1 every 1000
dump pic all image 1000 dump.*.ppm type type size 1024 1024 zoom 1.5 adiam 0.0015 adiam 0.0025
#dump stress all mesh/gran/VTK 1000 post/dump*.vtk stress wear cad
#dump mov all movie 1000 movie.m4v type type size 1280 720
#insert particles
run 10000 upto
unfix ins
and I attached the mesh file which is in txt extension just change the extension to stl.
Attachment | Size |
---|---|
![]() | 282.02 KB |
![]() | 16.44 KB |