Im trying to insert particle into a cylindrical container and apply axial load on it over the top plate. But it is showing error
ERROR: Fix wall/gran (id geometry): More than one wall of type 'mesh' is not supported (/build/buildd/liggghts-3.0.3+repack/src/fix_wall_gran.cpp:505)
# Shear Cell example
### Initialization
# Preliminaries
units si
atom_style sphere
boundary f f f
newton off
communicate single vel yes
#processors 2 2 2
# Declare domain
region reg block -0.021 0.021 -0.021 0.021 -0.001 0.051 units box
create_box 2 reg
neighbor 0.0002 bin
neigh_modify delay 0
### Setup
# Material and interaction properties
fix m1 all property/global youngsModulus peratomtype 1.0e5 1.0e5
fix m2 all property/global poissonsRatio peratomtype 0.31 0.31
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 characteristicVelocity scalar 2.
# Physics
pair_style gran model hertz tangential history
pair_coeff * *
# Particle setup – insertions to be handled below
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 8200 radius constant 0.0001 volume_limit 4.18879e-15
fix pdd1 all particledistribution/discrete 63243 1 pts1 1.0
# Geometry
fix cad1 all mesh/surface/stress/servo file without_dimples.stl type 1 scale 1. stress on com 0. 0. 0. ctrlPV force axis 0. 0. 0. vel_max 1.0
fix cad2 all mesh/surface/stress file bottom_plate.stl type 1 scale 1. stress on
fix ins_mesh all mesh/surface/planar file factory.stl type 1 scale 1.
fix geometry2 all wall/gran model hertz tangential history mesh n_meshes 2 meshes cad1 cad2
### Detailed settings
# Integration
fix integrator all nve/sphere
# Gravity
fix grav all gravity 9.81 vector 0.0 0.0 -1.0
# Timestep
timestep 0.00002
# Thermodynamic output
thermo_style custom step atoms ke cpu
thermo 10000
thermo_modify lost ignore norm no
### Execution and further settings
# Particle insertion – large particles
fix ins1 all insert/stream seed 123481 distributiontemplate pdd1 nparticles 2000 particlerate 2000 overlapcheck yes vel constant 0. 0. -0.4005 all_in yes insertion_face ins_mesh extrude_length 0.01
run 50000
# Axial load
fix cad3 all mesh/surface/stress/servo file top_plate.stl type 1 scale 1. stress on com 0. 0. 0. ctrlPV force axis 0. 0. -117.75 vel_max 1.0
fix geometry2 all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad3
#parameters for gradually growing particle diameter
variable alphastart equal 0.25
variable alphatarget equal 0.50
variable growts equal 50000
variable growevery equal 40
variable relaxts equal 50000
# Set the dumps
dump dumpmesh all mesh/vtk 10000 dump*.vtk stress
dump dmp all custom 10000 dump.1 id type type x y z ix iy iz vx vy vz fx fy fz &
omegax omegay omegaz radius
#dump mov all movie 10000 Ullas/shearcell/movie.m4v type type size 1280 720
#calculate grow rate
variable Rgrowrate equal (${alphatarget}/${alphastart})^(${growevery}/(3.*${growts}))
print "The radius grow rate is ${Rgrowrate}"
#do the diameter grow
compute rad all property/atom radius
variable dgrown atom ${Rgrowrate}*2.*c_rad
fix grow all adapt ${growevery} atom diameter v_dgrown
#run
run ${growts}
#let the packing relax
unfix grow
run ${relaxts}