“…When the function g is linear and positive-definite, the formulation greatly simplifies, but it then allows unphysical interpenetration of the material in contact (when the interface is subject to compressive tractions). Interfacial nonlinearity may be introduced to model different situations of interest [see, for instance, applications in: fragmentation and decohesion (Camacho and Ortiz, 1996;Needleman, 1992;Rice and Wang, 1989), interactions between inclusions (Levy and Hardikar, 1999), bifurcation (Radi et al, 1999), composites (Levy and Dong, 1998;Lipton and Talbot, 2001), biomechanics (Mann et al, 1997;Gei et al, 2002)] and may avoid the interpenetration by the introduction of a suitable penalty in compression. Recently, efforts have been made to provide models of thick interfaces, where the boundary value problem consisting of a three-phase configuration is replaced by a problem which involves only two phases plus some matching condition simulating the interphase (Benveniste and Miloh, 2001;Hashin, 2002;Rubin and Benveniste, 2004).…”