Dear liggghts users,
I want to break multisphere through liggghts. I want to break multispheres when the z coordinate of the multsphere exceeds certain height(here .4 m).
i considered gravity in the +ve z coordinate.
here is my input code.
#Multisphere
atom_style sphere
atom_modify map array sort 0 0
boundary m m m
newton off
communicate single vel yes
#PUBLIC version only supports serial execution
processors * * *
units si
region reg block -0.3 0.3 -0.3 0.3 0.0 0.5 units box
create_box 1 reg
neighbor 0.004 bin
neigh_modify delay 0
region bc1 block -.21 .21 -.21 .21 0.0 0.35 units box
group nve_group region bc1
#Material properties required for new pair styles
fix m1 all property/global youngsModulus peratomtype 1.e7
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
fix m5 all property/global characteristicVelocity scalar 2.
fix z all property/atom z scalar yes no no 0.1
#New pair style
pair_style gran model hertz tangential history #Hertzian without cohesion
pair_coeff * *
timestep 0.00001
fix gravi all gravity 9.81 vector 0.0 0.0 1.0
fix zwalls all wall/gran model hertz tangential history primitive type 1 zplane 0.0
#distributions for insertion
fix pts1 all particletemplate/multisphere 15485863 atom_type 1 density constant 2500 nspheres 50 ntry 1000000 spheres file data/stone1.multisphere scale 0.001 type 1
fix pdd1 all particledistribution/discrete 15485867 1 pts1 1.0
#region and insertion
region bc cylinder z 0.0 0.0 0.25 0.0 0.2 units box
fix ins all insert/pack seed 32452843 distributiontemplate pdd1 vel constant 0. 0. 0. &
insert_every once overlapcheck yes region bc ntry_mc 10000 volumefraction_region 0.002
#integrator for multisphere rigid bodies
#fix integr all multisphere
#fix integr1 all multisphere/break allow_group_and_set yes trigger_threshold .001 trigger_timeStep 0 trigger_name v_b
#fix integr all multisphere/break allow_group_and_set yes trigger_threshold 10 trigger_timeStep 0 trigger_name v_npartilces
fix integr all multisphere/break allow_group_and_set yes trigger_threshold .4 trigger_timeStep 0 trigger_name f_z
fix 2 all nve/sphere
#output settings, include total thermal energy
compute 1 all property/atom x y z
compute 2 all erotate/sphere
fix ts all check/timestep/gran 1000 0.1 0.1
thermo_style custom step atoms #ke c_1 f_ts[1] f_ts[2] vol
thermo 1000
thermo_modify lost ignore norm no
run 1
#insert the first particles so that dump is not empty
dump dmp all custom/vtk 200 post/multi_*.vtk id type mol x y z vx vy vz fx fy fz omegax omegay omegaz radius
#insert particles
run 40000 upto
but its not working .
Any suggestion would be deeply appreciated. Thank you very much for your time.
regrads,
Tarun