An interface debonding law subject to viscous regularization for avoiding instability: Application to the delamination problems

Abstract The prediction of ply delamination in laminated composites is modeled by using interface elements. The numerical approach is based on the cohesive zone model, which is shown to provide an efficient description of the delamination growth. The theoretical study is performed in quasi-static regime and an implicit finite element scheme is used. A comprehensive 1D example shows that cohesive elements may induce numerical instability and that the use of a viscous regularization is relevant to suppress the instability and to obtain a reasonable numerical convergence. This paper describes a technique that can be used to introduce damping into cohesive zone finite element simulations of crack nucleation and growth, with a view to avoiding convergence difficulties in quasi-static finite element simulations. A new kind of viscous regularization is proposed, which applies to the bilinear debonding model and has a limited rate dependency. An experimental study is performed by using DCB (double cantilever beam) samples for the mode I interlaminar fracture toughness test; a good agreement between numerical predictions and experimental results is obtained.

[1]  P. Suppakul,et al.  The effect of weave pattern on the mode-I interlaminar fracture energy of E-glass/vinyl ester composites , 2002 .

[2]  F. Feyel,et al.  Interface debonding models: a viscous regularization with a limited rate dependency , 2001 .

[3]  Pedro P. Camanho,et al.  A damage model for the simulation of delamination in advanced composites under variable-mode loading , 2006 .

[4]  Pedro P. Camanho,et al.  An Interface Damage Model for the Simulation of Delamination Under Variable-Mode Ratio in Composite Materials , 2004 .

[5]  René de Borst,et al.  Computational modelling of delamination , 2006 .

[6]  P. Davies Fracture of Marine Composites , 1996 .

[7]  Stephen R Reid,et al.  A continuum damage model for delaminations in laminated composites , 2003 .

[8]  G. Sih,et al.  Mathematical theories of brittle fracture. , 1968 .

[9]  Alberto Corigliano,et al.  Damage analysis of interlaminar fracture specimens , 1995 .

[10]  Xiaopeng Xu,et al.  Numerical simulations of fast crack growth in brittle solids , 1994 .

[11]  A. A. Griffith The Phenomena of Rupture and Flow in Solids , 1921 .

[12]  Larsgunnar Nilsson,et al.  Simulation of delamination in fiber composites with a discrete cohesive failure model , 2001 .

[13]  P. Camanho,et al.  Mixed-Mode Decohesion Finite Elements for the Simulation of Delamination in Composite Materials , 2002 .

[14]  J. Reeder,et al.  An Evaluation of Mixed-Mode Delamination Failure Criteria , 1992 .

[15]  G. D. Sims,et al.  Comparison of test configurations for determination of mode II interlaminar fracture toughness results from international collaborative test programme , 1999 .

[16]  M. Crisfield,et al.  Finite element interface models for the delamination analysis of laminated composites: mechanical and computational issues , 2001 .

[17]  Rhj Ron Peerlings,et al.  Gradient enhanced damage for quasi-brittle materials , 1996 .

[18]  Alberto Corigliano,et al.  Modeling and simulation of crack propagation in mixed-modes interlaminar fracture specimens , 1996 .

[19]  V. Tvergaard Effect of fibre debonding in a whisker-reinforced metal , 1990 .

[20]  Yanfei Gao,et al.  A simple technique for avoiding convergence problems in finite element simulations of crack nucleation and growth on cohesive interfaces , 2004 .

[21]  O. Allix,et al.  Interlaminar interface modelling for the prediction of delamination , 1992 .

[22]  Pierre Ladevèze,et al.  A mesomodel for localisation and damage computation in laminates , 2000 .

[23]  M. Kanninen,et al.  A finite element calculation of stress intensity factors by a modified crack closure integral , 1977 .

[24]  Laurent Blanchard,et al.  Mesomodeling of delamination: towards industrial applications , 2006 .

[25]  A. Hillerborg,et al.  Analysis of crack formation and crack growth in concrete by means of fracture mechanics and finite elements , 1976 .

[26]  P. Geubelle,et al.  Impact-induced delamination of composites: A 2D simulation , 1998 .

[27]  Olivier Allix,et al.  A delay damage mesomodel of laminates under dynamic loading: basic aspects and identification issues , 2003 .

[28]  M. D. Moura,et al.  Modelling single and double-lap repairs on composite materials , 2005 .

[29]  P.M.S.T. de Castro,et al.  Prediction of compressive strength of carbon–epoxy laminates containing delamination by using a mixed-mode damage model , 2000 .

[30]  J. Mazars APPLICATION DE LA MECANIQUE DE L'ENDOMMAGEMENT AU COMPORTEMENT NON LINEAIRE ET A LA RUPTURE DU BETON DE STRUCTURE , 1984 .

[31]  J. C. Simo,et al.  Strain- and stress-based continuum damage models—I. Formulation , 1987 .

[32]  A. Needleman Material rate dependence and mesh sensitivity in localization problems , 1988 .

[33]  A. Needleman An analysis of tensile decohesion along an interface , 1990 .

[34]  G. I. Barenblatt THE MATHEMATICAL THEORY OF EQUILIBRIUM CRACKS IN BRITTLE FRACTURE , 1962 .