Numerical modeling of elastic waves across imperfect contacts.

A numerical method is described for studying how elastic waves interact with imperfect contacts such as fractures or glue layers existing between elastic solids. These contacts have been classicaly modeled by interfaces, using a simple rheological model consisting of a combination of normal and tangential linear springs and masses. The jump conditions satisfied by the elastic fields along the interfaces are called the "spring-mass conditions." By tuning the stiffness and mass values, it is possible to model various degrees of contact, from perfect bonding to stress-free surfaces. The conservation laws satisfied outside the interfaces are integrated using classical finite-difference schemes. The key problem arising here is how to discretize the spring-mass conditions and how to insert them into a finite-difference scheme: this was the aim of the present paper. For this purpose, we adapted an interface method previously developed for use with perfect contacts [J. Comput. Phys., 195 (2004), pp. 90-116]. This numerical method also describes closely the geometry of arbitrarily shaped interfaces on a uniform Cartesian grid, at negligible extra computational cost. Comparisons with original analytical solutions show the efficiency of this approach.