Dear all
I’m trying to perform the 3D crack simulation of a metal using peridynamic package.
While doing so, only atom is being created and hence I cannot apply velocity to the system.
I am sharing the input script here, kindly take a look.
Thank you
2d peridynamic crack simulation
units si
dimension 3
boundary s s p
atom_style peri
neighbor 0.0001 bin
neigh_modify delay 5
create geometry
lattice bcc 0.00287
#lattice structure and constant will change with metal
region box block 0 0.1 0 0.1 -0.01 0.01
create_box 5 box
create_atoms 1 box
#mass 1 1.0
#mass 2 1.0
#mass 3 1.0
#mass 4 1.0
#mass 5 1.0
peridynamic potentials
pair_style peri/eps
pair_coeff * * 1.6e11 8.4e10 0.0015001 0.0005 0.25 25e7
define groups
region 1 block INF INF INF 0.001 INF INF
group lower region 1
region 2 block INF INF 0.099 INF INF INF
group upper region 2
group boundary union lower upper
group mobile subtract all boundary
region leftupper block INF 0.05 0.05 INF INF INF
region leftlower block INF 0.05 INF 0.05 INF INF
group leftupper region leftupper
group leftlower region leftlower
set group leftupper type 2
set group leftlower type 3
set group lower type 4
set group upper type 5
initial velocities
compute new mobile temp
velocity mobile create 0.01 887723 temp new
velocity upper set 0.0 0.3 0.0
velocity mobile ramp vy 0.0 0.3 y 0.001 0.099 sum yes
compute stress all stress/atom NULL
fixes
fix 1 all nve
fix 2 boundary setforce NULL 0.0 0.0
run
timestep 0.003
thermo 200
thermo_modify temp new
neigh_modify exclude type 2 3
#dump 1 all atom 500 dump-v3.crack
dump 1 all custom 100 dump-v3.crack.* id type x y z vx vy vz c_stress[1] c_stress[2] c_stress[3]
#dump 2 all image 250 image.*.jpg type type &
zoom 1.6 adiam 1.5
#dump_modify 2 pad 4
dump 3 all movie 100 dump-v3.mp4 type type &
zoom 1.6 adiam 1.5
dump_modify 3 pad 4
run 10000