Moving up a mesh with a constant velocity

Submitted by shademani on Thu, 12/13/2018 - 18:33

Dear all,
I want to simulate a case in which a gate move up. I used move mesh command.
However, the mesh did not move.

#Multisphere

atom_style sphere
atom_modify map array sort 0 0
boundary f f f
newton off

communicate single vel yes

#PUBLIC version only supports serial execution
processors 1 1 1

units si

region reg block -0.5 0.3 -0.2 0.2 0. 0.5 units box
create_box 1 reg

neighbor 0.004 bin
neigh_modify delay 0

#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.

#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

fix zwalls all wall/gran model hertz tangential history primitive type 1 zplane 0
fix xwalls all wall/gran model hertz tangential history primitive type 1 xplane -0.5
fix xxwalls all wall/gran model hertz tangential history primitive type 1 xplane 0.3
fix ywalls all wall/gran model hertz tangential history primitive type 1 yplane 0.2
fix yywalls all wall/gran model hertz tangential history primitive type 1 yplane -0.2
#import mesh from cad:
fix cad1 all mesh/surface file meshes/gate.stl type 1 move -0.18 0 0 scale 1
#distributions for insertion
fix pts1 all particletemplate/multisphere 15485863 atom_type 1 density constant 2500 nspheres 25 ntry 1000000 spheres file data/stone1.multisphere scale 0.001 type 1
fix pdd1 all particledistribution/discrete 15485867 1 pts1 1.0
#fix cv all wall/gran model hertz tangential history primitive type 1 xplane -0.16
#region and insertion
region bc cylinder z -0.35 0 0.14 0.2 0.4 units box

fix ins all insert/pack seed 32452843 distributiontemplate pdd1 vel constant 0. 0. -1. &
insert_every once overlapcheck yes region bc ntry_mc 10000 volumefraction_region 1

#integrator for multisphere rigid bodies
fix integr all multisphere

#output settings, include total thermal energy
compute 1 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
Dear all,
I want to simulate a case in which a gate move up with constant velocity. But the mesh did not move.
could you help me.

here is my code:
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
dump dumpstl all mesh/stl 200 post/dump*.stl
#use the imported mesh as granular wall
fix gate_wall all wall/gran model hooke tangential history mesh n_meshes 1 meshes cad1

#insert particles
run 40000 upto
#upto

#run 22000
#run 30000 upto

moving mesh
fix movecad1 all move/mesh mesh cad1 linear 0 0 1
run 60000