# small Peridynamic cylinder hit by projectile atom_style peri atom_modify map array read_restart restart.new neighbor 0.0010 bin #read_restart restart.new # small target lattice sc 0.0005 region target block -0.0025 0.0025 -0.0025 0.0025 -0.0025 0.0025 units box group 1 id 1 group 2 id 551 312 group 3 id 179 189 190 233 234 244 245 180 181 191 192 235 236 246 247 193 237 248 200 201 211 212 255 256 266 267 202 203 213 214 257 258 268 269 204 215 259 270 288 289 299 300 343 344 354 355 399 409 410 279 290 291 301 302 345 346 356 357 400 401 411 412 292 303 347 358 413 309 310 311 321 322 365 366 376 377 420 421 431 432 313 323 324 367 368 378 379 422 423 433 434 314 315 325 369 380 424 435 477 806 871 872 807 808 873 874 875 937 938 1003 1004 1070 939 940 950 1005 1006 1071 1072 941 1007 441 442 443 452 453 454 463 464 496 497 498 507 508 509 518 519 520 552 553 562 563 564 573 574 575 554 584 #delete_bonds 2 multi #neigh_modify exclude group 2 3 set group all density 2200 set group all volume 1.25e-10 #velocity all set 0.0 0.0 0.0 sum no units box fix 1 all nve # spherical indenter to shatter target #variable y0 equal 0.00155 #variable vy equal -100 #variable y equal "v_y0 + step*dt*v_vy" #fix 2 all indent 1e17 sphere 0.0000 v_y 0.0000 0.0015 units box #unfix 2 #fix 2 1 setforce 0.0 0.0 0.0 #unfix 2 fix 2 1 addforce 1e9 0.0 0.0 compute 1 all damage/atom timestep 1.0e-7 thermo 100 dump 1 all custom 1 dump.new3 id type x y z c_1 vx vy vz fx fy fz run 100 write_restart restart.new2