Peridynamics mesoscale modeling of impact fracture

This is work of Mike Parks at Sandia. We implemented a Peridynamics (PD) model in LAMMPS to enable meso- and continuum-scale simulations of materials response. PD is a particle-based meshless continuum model that is well-suited for hi-deformation problems, such as impact and fracture scenarios. The theory underlying PD was developed by Stewart Silling at Sandia.

In these images a 1 cm ball (not shown) moving at 100 m/sec impacts a circular disk of brittle material composed of about 100K particles. The timestep is 1 nsec and the simulation runs for 200K timesteps.