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