1# small Peridynamic cylinder hit by projectile 2 3units si 4boundary s s s 5atom_style peri 6atom_modify map array 7neighbor 0.0010 bin 8 9# small target 10 11lattice sc 0.0005 12region target cylinder y 0.0 0.0 0.0050 -0.0050 0.0 units box 13create_box 1 target 14create_atoms 1 region target 15 16pair_style peri/eps 17pair_coeff * * 14.9e9 14.9e9 0.0015001 0.0005 0.25 10.0e8 18set group all density 2200 19set group all volume 1.25e-10 20velocity all set 0.0 0.0 0.0 sum no units box 21fix 1 all nve 22 23# spherical indenter to shatter target 24 25variable y0 equal 0.00155 26variable vy equal -100 27variable y equal "v_y0 + step*dt*v_vy" 28 29fix 2 all indent 1e17 sphere 0.0000 v_y 0.0000 0.0015 units box 30 31compute 1 all damage/atom 32timestep 1.0e-7 33thermo 100 34 35#dump 1 all custom 100 dump.peri id type x y z c_1 36 37#dump 2 all image 50 image.*.jpg type type & 38# axes yes 0.8 0.02 view 80 -30 adiam 0.0006 39#dump_modify 2 pad 4 40 41#dump 3 all movie 50 movie.mpg type type & 42# axes yes 0.8 0.02 view 80 -30 adiam 0.0006 43#dump_modify 3 pad 4 44 45run 1000 46