MIXED FINITE ELEMENT FOR CRACKED INTERFACE

Several disorders observed in an existing work on civil engineering may have their origin in local phenomena which reveal the weak points of this work. These critical zones are located, on the one hand in the links between materials or interfaces, on the other hand in singularly formed areas such as cavities, angles and cracks, seats of strong stress concentrations. In this paper, the finite element method (mixed formulation) is used for the study of the interfacial cracks in bimaterials. A special finite element based on the mixed formulation, able to take into account the continuity of the interface on the coherent part (mechanical and geometrical continuity), and the discontinuity of this one on the cracked part (edge effect), is used to model at best this type of interface. This element was developed by H. Bouzerd using a direct formulation: the shape functions of the displacement and stress fields are built directly starting from the real configuration of the element in a reference (x, y). In this purpose, this element was reformulated starting from a reference element in a natural plan (ζ,η). This formulation presents, in addition to the simplification of calculations, the enormous advantage of modelling the types of cracks and their orientations. This interface element was associated with the virtual extension—crack to evaluate the energy release rates using only one meshing by finite elements. Several numerical examples concerning the interfacial cracks are analysed to assess the validity of this element by comparing with an available analytical solution or numerical ones obtained from others finite elements.

[1]  R. Mahapatra,et al.  The mixed finite element method in elastic and elasto-plastic axisymmetric problems , 1988 .

[2]  A. H. Shah,et al.  Mixed Finite-Element Model for Thick Composite Laminated Plates , 2002 .

[3]  J. Rice,et al.  Plane Problems of Cracks in Dissimilar Media , 1965 .

[4]  Erasmo Carrera,et al.  Transverse Normal Stress Effects in Multilayered Plates , 1999 .

[5]  D. M. Parks A stiffness derivative finite element technique for determination of crack tip stress intensity factors , 1974 .

[6]  Y. M. Desai,et al.  Hybrid-interface finite element for laminated composite and sandwich beams , 2007 .

[7]  F. A. Mirza,et al.  The mixed finite element method in plane elasticity , 1980 .

[8]  N. J. Pagano,et al.  Exact solutions for rectangular bidirectional composites and sandwich plates : Journal of Composite Materials, Vol 4, pp 20 – 34 (January 1970) , 1970 .

[9]  A. H. Shah,et al.  Application of a three-dimensional mixed finite element model to the flexure of sandwich plate , 2003 .

[10]  Erasmo Carrera,et al.  Mixed layer-wise models for multilayered plates analysis , 1998 .

[11]  E. Carrera C0 REISSNER–MINDLIN MULTILAYERED PLATE ELEMENTS INCLUDING ZIG-ZAG AND INTERLAMINAR STRESS CONTINUITY , 1996 .

[12]  K. Lin,et al.  Finite element analysis of stress intensity factors for cracks at a bi-material interface , 1976 .

[13]  Shohei Nakazawa Hybrid and mixed finite element methods: S.N. Atluri, R.H. Gallagher and O.C. Zienkiewicz, eds. (Wiley, Chichester, 1983) , 1984 .

[14]  Y. Desai,et al.  Mixed finite element model for laminated composite beams , 2002 .