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

• Published in: Engineering fracture mechanics Vol. 75; no. 10; pp. 3084 - 3100
• Main Authors: Hamitouche, L., Tarfaoui, M., Vautrin, A.
• Format: Journal Article
• Language: English
• Published: Tarrytown, NY Elsevier Ltd 01.07.2008; Oxford Elsevier
• Subjects: Cohesive zone; DCB; Delamination; Engineering Sciences; Exact sciences and technology; Fracture mechanics (crack, fatigue, damage...); Fundamental areas of phenomenology (including applications); Materials; Physics; Solid mechanics; Structural and continuum mechanics; Viscous regularization; Viscous regularization; DCB; Delamination; Cohesive zone; Stratified material; Nucleation; Cohesive end; Quasi static theory; Experimental study; Cavitation; Modeling; Crack propagation; Finite element method; Interface crack; Fracture toughness test; Regularization; Numerical convergence; Double cantilever beam; Numerical stability; Damaging; LOCALIZATION; GROWTH; MESOMODEL; COMPOSITES; SIMULATION; DAMAGE MODEL; PREDICTION
• ISSN: 0013-7944; 1873-7315
• DOI: 10.1016/j.engfracmech.2007.12.014

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.