Damage initialization techniques for non-sequential FE propagation analysis of delaminations in composite aerospace structures

The experimental effort required to develop, damage tolerant, aerospace composite structures could be significantly reduced if reliable numerical simulations were used to perform engineering studies of complex damaged structures. Finite element (FE) simulations of impact damaged structures typically follow a sequential approach that require large computational resources to reproduce complex damage scenarios. A numerical tool capable to reconstruct such scenarios using data from previous impact simulations or NDI could noticeably improve the simulation workflow for damaged composite structures. The paper proposes a method to inizialize the damage variables in numerical analyses aimed at assessing damage propagation, and that are potentially able to evaluate the residual strength of damaged structures. The approach is developed within FE software ABAQUS, and uses SDVINI subroutine to initialize damage variables defined by a user-material-subroutine (UMAT), that provides the constitutive models of the lamina and of the interlaminar layers. Albeit the proposed technique might deal with both inter-laminar and intra-laminar damage, the paper is focused on delaminations. A user defined traction-separation law is coded in an UMAT that endows ABAQUS cohesive elements with damage initialization capabilities. Then, results of test cases, of increasing complexity, are presented in order to assess the damage initialization procedure and verify the performances of its different operating modes. Two test-cases are based on plate-like specimens for which literature data exist: the first is relevant to a circular artificial delamination while the second presents multiple delaminations caused by an impact and measured via NDI techniques. The last test-case is a stiffened panel which incorporates the typical complexities of aerospace structures, but is still tractable with the sequential simulation approach whose results are used as a term of comparison.

[1]  P. Camanho,et al.  Numerical Simulation of Mixed-Mode Progressive Delamination in Composite Materials , 2003 .

[2]  M. Benzeggagh,et al.  Measurement of mixed-mode delamination fracture toughness of unidirectional glass/epoxy composites with mixed-mode bending apparatus , 1996 .

[3]  Lorenzo Iannucci,et al.  A progressive failure model for mesh-size-independent FE analysis of composite laminates subject to low-velocity impact damage , 2012 .

[4]  Pedro P. Camanho,et al.  Accurate simulation of delamination growth under mixed-mode loading using cohesive elements: Definition of interlaminar strengths and elastic stiffness , 2010 .

[5]  Dirk Vandepitte,et al.  Failure analysis of low velocity impact on thin composite laminates : Experimental and numerical approaches , 2008 .

[6]  Francesco Aymerich,et al.  Finite element modelling of damage induced by low-velocity impact on composite laminates , 2014 .

[7]  S. Hallett,et al.  A study on the influence of cohesive zone interface element strength parameters on mixed mode behaviour , 2012 .

[8]  Daniele Fanteria,et al.  A non-linear shear damage model to reproduce permanent indentation caused by impacts in composite laminates , 2014 .

[9]  J. Ankersen,et al.  Virtual testing of realistic aerospace composite structures , 2008 .

[10]  D. S. Dugdale.  Yielding of steel sheets containing slits , 1960 .

[11]  de R René Borst,et al.  On the numerical integration of interface elements , 1993 .

[12]  Brian Falzon,et al.  A progressive failure model for composite laminates subjected to low velocity impact damage , 2008 .

[13]  Zafer Gürdal,et al.  Low-velocity impact damage on dispersed stacking sequence laminates. Part II: Numerical simulations , 2009 .

[14]  Prasad B. Chunchu,et al.  POSTBUCKLING AND GROWTH OF DELAMINATIONS IN COMPOSITE PLATES SUBJECTED TO AXIAL COMPRESSION , 2002 .

[15]  Luis Reis,et al.  Failure mechanisms on composite specimens subjected to compression after impact , 1998 .

[16]  Zafer Gürdal,et al.  LOW-VELOCITY IMPACT DAMAGE ON DISPERSED STACKING SEQUENCE LAMINATES. PART I: EXPERIMENTS , 2009 .

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

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

[19]  Brian Falzon,et al.  Predicting low-velocity impact damage on a stiffened composite panel , 2010 .

[20]  Pedro P. Camanho,et al.  An engineering solution for mesh size effects in the simulation of delamination using cohesive zone models , 2007 .

[21]  Constantinos Soutis,et al.  Modelling damage evolution in composite laminates subjected to low velocity impact , 2012 .

[22]  Michael R Wisnom,et al.  Modelling discrete failures in composites with interface elements , 2010 .

[23]  Cheryl A. Rose,et al.  Guidelines and Parameter Selection for the Simulation of Progressive Delamination , 2008 .

[24]  Pedro P. Camanho,et al.  Effective simulation of delamination in aeronautical structures using shells and cohesive elements , 2008 .