A Lagrangian DG-Method for Wave Propagation in a Cracked Solid with Frictional Contact Interfaces

We developed a discontinuous Galerkin (DG) numerical scheme for wave propagation in elastic solids with frictional contact interfaces. This type of numerical scheme is useful in investigations of wave propagation in elastic solids with micro-cracks (cracked solid) that involve modeling the damage in brittle materials or architected meta-materials. Only processes with mild loading that do not trigger crack fracture extension or the nucleation of new fractures are considered. The main focus lies on the contact conditions at crack surfaces, including provisions for crack opening and closure and stick-and-slip with Coulomb friction. The proposed numerical algorithm uses the leapfrog scheme for the time discretization and an augmented Lagrangian algorithm to solve the associated non-linear problems. For efficient parallelization, a Galerkin discontinuous method was chosen for the space discretization. The frictional interfaces (micro-cracks), where the numerical flux is obtained by solving non-linear and non-smooth variational problems, concerns only a limited number the degrees of freedom, implying a small additional computational cost compared to classical bulk DG schemes. The numerical method was tested through two model problems with analytical solutions. The proposed Lagrangian approach of the nonlinear interfaces had excellent results (stability and high accuracy) and only required a reasonable additional amount of computational effort. To illustrate the method, we conclude with some numerical simulations on the blast propagation in a cracked material and in a meta-material designed for shock dissipation.