Hello,
I want to insert particles in cad files. after that compress it and rotate it. When I'm trying to execute the program. particles are falling down to cad files. .Plate is not rotating. Can anyone help me . Please find the code below:
### Initialization
# Preliminaries
soft_particles yes
units cgs
atom_style granular
boundary f f f
atom_modify map array sort 10000 2.0
neigh_modify delay 0 every 1 check yes page 500000 one 50000
newton off
communicate single vel yes
processors * * *
# Domain
region reg block -15 15 -3.6 3.6 -15 15 units box
create_box 3 reg
### Setup
# Material and interaction properties required
fix m1 all property/global youngsModulus peratomtype 250 250 250
fix m2 all property/global poissonsRatio peratomtype 0.25 0.25 0.25
fix m3 all property/global coefficientRestitution peratomtypepair 3 1 0.5 1 0.5 1 1 1 1 1
fix m4 all property/global coefficientFriction peratomtypepair 3 1 0.5 1 0.5 1 1 1 1 1
fix m5 all property/global coefficientRollingFriction peratomtypepair 3 1 0.5 1 0.5 1 1 1 1 1
fix m6 all property/global characteristicVelocity scalar 0.1
# Define the physics
pair_style gran model hooke tangential history rolling_friction cdt
pair_coeff * *
###Detailed Settings
#Integrator
fix integrate all nve/sphere
#TimeStep
timestep 0.000001
fix ts all check/timestep/gran 1000 0.1 0.1
#Define the Gravity
fix grav all gravity 980.665 vector 0.0 -1.0 0.0
###Particle Insertion
region bc1 cylinder y 0 0 12.24 -3.60 0 units box
fix pts1 all particletemplate/sphere 12347 atom_type 1 density constant 1.8 radius constant 0.0047
fix pts2 all particletemplate/sphere 12373 atom_type 1 density constant 1.8 radius constant 0.00345
fix pts3 all particletemplate/sphere 12377 atom_type 1 density constant 1.8 radius constant 0.00255
fix pts4 all particletemplate/sphere 12379 atom_type 1 density constant 1.8 radius constant 0.0019
fix pts5 all particletemplate/sphere 12391 atom_type 1 density constant 1.8 radius constant 0.0014
fix pts6 all particletemplate/sphere 12401 atom_type 1 density constant 1.8 radius constant 0.001
fix pts7 all particletemplate/sphere 12409 atom_type 1 density constant 1.8 radius constant 0.00075
fix pts8 all particletemplate/sphere 12413 atom_type 1 density constant 1.8 radius constant 0.00055
fix pts9 all particletemplate/sphere 12421 atom_type 1 density constant 1.8 radius constant 0.0004
fix pts10 all particletemplate/sphere 12433 atom_type 1 density constant 1.8 radius constant 0.0003
fix pts11 all particletemplate/sphere 12437 atom_type 1 density constant 1.8 radius constant 0.0002
fix pts12 all particletemplate/sphere 12451 atom_type 1 density constant 1.8 radius constant 0.00015
fix pts13 all particletemplate/sphere 12457 atom_type 1 density constant 1.8 radius constant 0.0001
fix pdd1 all particledistribution/discrete/numberbased 15485959 13 pts1 0.007 pts2 0.045 pts3 0.1004 pts4 0.1322 pts5 0.1244 pts6 0.1007 pts7 0.0794 pts8 0.0655 pts9 0.0608 pts10 0.0601 pts11 0.0556 pts12 0.0449 pts13 0.124
fix ins all insert/rate/region seed 32452843 distributiontemplate pdd1 nparticles 100000000 massrate 1.66666667 insert_every 10 overlapcheck yes vel constant 0 -1 0 region bc1 ntry_mc 10000
fix cad3 all mesh/surface file Simu/plate2.stl type 1 scale 0.1 curvature 1e-5
fix meshwalls1 all wall/gran model hooke tangential history mesh n_meshes 1 meshes cad3
run 1
dump dumpstl all stl 10000 post1/dump*.stl
dump dmp all custom 10000 post1/dump*.flow id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
run 50000
unfix ins
unfix meshwalls1
fix cad4 all mesh/surface/stress/servo file Simulation/top2.stl type 1 scale 0.01 curvature 1e-5 com 0 0 0 ctrlPV force axis 0 -1 0 target_val 100e6 vel_max 1
fix meshwalls1 all wall/gran model hooke tangential history mesh n_meshes 2 meshes cad4 cad3
run 1
unfix meshwalls1
unfix cad4
fix cad2 all mesh/surface/stress/servo file Simulation/rubber_2.stl type 1 scale 0.1 curvature 1e-5 com 0 0 0 ctrlPV force axis 0 -1 0 target_val 3e6 vel_max 1
fix meshwalls2 all wall/gran model hooke tangential history mesh n_meshes 2 meshes cad2 cad3
dump dumpstl2 all stl 10000 post1/dump*.stl
fix momemti all ave/time 1 1 10000 f_cad2[3] f_cad2[6] file post/force_of_cad2.profile
run 50000
###Rotate the Soil Plate
fix movecad3 all move/mesh mesh cad3 rotate origin 0 -3.6 0 axis 0 1 0 period 38.4846
##Execution and further settings
#run remainder
run 20000000 upto