  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)
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.
