Two- and three-dimensional isogeometric cohesive elements for composite delamination analysis

We propose an automatic numerical method requiring minimal user intervention to simulate delamination in composite structures. We develop isogeometric cohesive elements for two- and three-dimensional delamination by exploiting the knot insertion algorithm directly from CAD data to generate cohesive elements along delamination. A complete computational framework is presented including pre-processing, processing and post-processing. They are explained in detail and implemented in MIGFEM – an open source Matlab Isogemetric Analysis code developed by the authors. The composite laminates are modeled using both NURBS solid and rotation-free shell elements. Several two and three dimensional examples ranging from standard delamination tests (the mixed mode bending test) to the L-shaped specimen with a fillet, three dimensional (3D) double cantilever beam and a 3D singly curved thick-walled laminate are provided. The method proposed provides a bi-directional system in which one can go forward from CAD to analysis and backwards from analysis to CAD. This is believed to facilitate the design of composite structures.

[1]  Cv Clemens Verhoosel,et al.  Modelling inter- and transgranular fracture in piezoelectric polycrystals , 2009 .

[2]  F. Stan.  Discontinuous Galerkin Method for Interface Crack Propagation , 2008 .

[3]  De Xie,et al.  Discrete cohesive zone model for mixed-mode fracture using finite element analysis , 2006 .

[4]  Ronald Krueger,et al.  The Virtual Crack Closure Technique : History , Approach and Applications , 2002 .

[5]  S. Bordas,et al.  A robust preconditioning technique for the extended finite element method , 2011 .

[6]  T. Hughes,et al.  Isogeometric analysis : CAD, finite elements, NURBS, exact geometry and mesh refinement , 2005 .

[7]  M. A. Crisfield,et al.  Progressive Delamination Using Interface Elements , 1998 .

[8]  Jens Gravesen,et al.  Isogeometric Shape Optimization of Vibrating Membranes , 2011 .

[9]  de R René Borst,et al.  Delamination buckling of fibre–metal laminates , 2001 .

[10]  T. Belytschko,et al.  A generalized finite element formulation for arbitrary basis functions: From isogeometric analysis to XFEM , 2010 .

[11]  L. J. Sluys,et al.  Simulation of fracture of cementitious composites with explicit modeling of microstructural features , 2001 .

[12]  T. Hughes,et al.  Isogeometric analysis of the Cahn–Hilliard phase-field model , 2008 .

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

[14]  G. LaPlante,et al.  Modeling delamination growth in composites under fatigue loadings of varying amplitudes , 2012 .

[15]  Michael R Wisnom,et al.  The open hole tensile test: a challenge for virtual testing of composites , 2009 .

[16]  Stéphane Bordas,et al.  An extended finite element library , 2007 .

[17]  Vinh Phu Nguyen,et al.  Isogeometric analysis suitable trivariate NURBS representation of composite panels with a new offset algorithm , 2014, Comput. Aided Des..

[18]  Anath Fischer,et al.  New B‐Spline Finite Element approach for geometrical design and mechanical analysis , 1998 .

[19]  Arif Masud,et al.  B-splines and NURBS based finite element methods for Kohn–Sham equations , 2012 .

[20]  G.A.O. Davies,et al.  Element-Free Galerkin modelling of composite damage , 2009 .

[21]  T. Belytschko,et al.  X‐FEM in isogeometric analysis for linear fracture mechanics , 2011 .

[22]  Thomas J. R. Hughes,et al.  A large deformation, rotation-free, isogeometric shell , 2011 .

[23]  I. Akkerman,et al.  Large eddy simulation of turbulent Taylor-Couette flow using isogeometric analysis and the residual-based variational multiscale method , 2010, J. Comput. Phys..

[24]  Thomas J. R. Hughes,et al.  Patient-specific isogeometric fluid–structure interaction analysis of thoracic aortic blood flow due to implantation of the Jarvik 2000 left ventricular assist device , 2009 .

[25]  Thomas J. R. Hughes,et al.  An isogeometric approach to cohesive zone modeling , 2011 .

[26]  R. H. Dodds,et al.  Simulation of ductile crack growth in thin aluminum panels using 3-D surface cohesive elements , 2001 .

[27]  Roland Wüchner,et al.  Isogeometric shell analysis with Kirchhoff–Love elements , 2009 .

[28]  Ole Sigmund,et al.  Isogeometric shape optimization of photonic crystals via Coons patches , 2011 .

[29]  Utkarsh Ayachit,et al.  The ParaView Guide: A Parallel Visualization Application , 2015 .

[30]  Alessandro Reali,et al.  GeoPDEs: A research tool for Isogeometric Analysis of PDEs , 2011, Adv. Eng. Softw..

[31]  W. Wall,et al.  Isogeometric structural shape optimization , 2008 .

[32]  Jia Lu,et al.  Isogeometric contact analysis: Geometric basis and formulation for frictionless contact , 2011 .

[33]  Manfred Bischoff,et al.  A point to segment contact formulation for isogeometric, NURBS based finite elements , 2013 .

[34]  Peter Wriggers,et al.  Contact treatment in isogeometric analysis with NURBS , 2011 .

[35]  I. Babuska,et al.  The Partition of Unity Method , 1997 .

[36]  Paolo Ermanni,et al.  Enhanced model for interlaminar normal stress in singly curved laminates , 2007 .

[37]  M. A. Gutiérrez.  Energy release control for numerical simulations of failure in quasi‐brittle solids , 2004 .

[38]  L. Beirao da Veiga,et al.  An isogeometric method for the Reissner-Mindlin plate bending problem , 2011, 1106.4436.

[39]  Thomas J. R. Hughes,et al.  Isogeometric shell analysis: The Reissner-Mindlin shell , 2010 .

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

[41]  Sundararajan Natarajan,et al.  Integrating strong and weak discontinuities without integration subcells and example applications in an XFEM/GFEM framework , 2010, 1107.4732.

[42]  John A. Evans,et al.  Isogeometric finite element data structures based on Bézier extraction of NURBS , 2011 .

[43]  Michael Ortiz,et al.  An Efficient Adaptive Procedure for Three-Dimensional Fragmentation Simulations , 2001, Engineering with Computers.

[44]  Heinz E. Pettermann,et al.  A semi-analytical model for the simulation of delamination in laminated composites , 2008 .

[45]  Ludovic Noels,et al.  A micro-meso-model of intra-laminar fracture in fiber-reinforced composites based on a Discontinuous Galerkin/Cohesive Zone Method , 2013 .

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

[47]  T. Hughes,et al.  Efficient quadrature for NURBS-based isogeometric analysis , 2010 .

[48]  R. Echter,et al.  A hierarchic family of isogeometric shell finite elements , 2013 .

[49]  van Jaw Hans Dommelen,et al.  A self-adaptive finite element approach for simulation of mixed-mode delamination using cohesive zone models , 2011 .

[50]  N. Valizadeh,et al.  Extended isogeometric analysis for simulation of stationary and propagating cracks , 2012 .

[51]  D. Borst,et al.  A NON-LINEAR FINITE ELEMENT APPROACH FOR THE ANALYSIS OF MODE-I FREE EDGE DELAMINATION IN COMPOSITES , 1993 .

[52]  Olivier Allix,et al.  A three-scale domain decomposition method for the 3D analysis of debonding in laminates , 2009, 1109.6111.

[53]  G. Subbarayan,et al.  Isogeometric enriched field approximations , 2012 .

[54]  J. Prévost,et al.  Modeling quasi-static crack growth with the extended finite element method Part I: Computer implementation , 2003 .

[55]  John A. Evans,et al.  An Isogeometric design-through-analysis methodology based on adaptive hierarchical refinement of NURBS, immersed boundary methods, and T-spline CAD surfaces , 2012 .

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

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

[58]  M. Ortiz,et al.  Computational modelling of impact damage in brittle materials , 1996 .

[59]  Thomas J. R. Hughes,et al.  Blended isogeometric shells , 2013 .

[60]  Alessandro Reali,et al.  Isogeometric Analysis of Structural Vibrations , 2006 .

[61]  Garth N. Wells,et al.  A solid‐like shell element allowing for arbitrary delaminations , 2003 .

[62]  Peter Wriggers,et al.  A large deformation frictional contact formulation using NURBS‐based isogeometric analysis , 2011 .

[63]  Niels Leergaard Pedersen,et al.  Discretizations in isogeometric analysis of Navier-Stokes flow , 2011 .

[64]  John A. Evans,et al.  Isogeometric analysis using T-splines , 2010 .

[65]  Hiroshi Suemasu,et al.  STRESS ANALYSES OF COMPOSITE LAMINATE WITH DELAMINATION USING X-FEM , 2006 .

[66]  Thomas J. R. Hughes,et al.  An isogeometric analysis approach to gradient damage models , 2011 .

[67]  Giulio Alfano,et al.  Solution strategies for the delamination analysis based on a combination of local‐control arc‐length and line searches , 2003 .

[68]  G.A.O. Davies,et al.  Decohesion finite element with enriched basis functions for delamination , 2009 .

[69]  Alessandro Reali,et al.  Duality and unified analysis of discrete approximations in structural dynamics and wave propagation : Comparison of p-method finite elements with k-method NURBS , 2008 .

[70]  Zafer Gürdal,et al.  A layerwise theory for laminated composites in the framework of isogeometric analysis , 2014 .

[71]  B. Simeon,et al.  A hierarchical approach to adaptive local refinement in isogeometric analysis , 2011 .

[72]  Krueger Ronald,et al.  A Shell/3D Modeling Technique for the Analysis of Delaminated Composite Laminates , 2000 .

[73]  Xiaoping Qian.  Full analytical sensitivities in NURBS based isogeometric shape optimization , 2010 .

[74]  Sung-Kie Youn,et al.  T‐spline finite element method for the analysis of shell structures , 2009 .

[75]  Pierre Kerfriden,et al.  Nitsche's method method for mixed dimensional analysis: conforming and non-conforming continuum-beam and continuum-plate coupling , 2013, 1308.2910.

[76]  J. L. Curiel Sosa,et al.  Delamination modelling of GLARE using the extended finite element method , 2012 .

[77]  Paul Steinmann,et al.  Isogeometric analysis of 2D gradient elasticity , 2011 .

[78]  Constantinos Soutis,et al.  Strength prediction of bolted joints in CFRP composite laminates using cohesive zone elements , 2014 .

[79]  Thomas J. R. Hughes,et al.  Isogeometric divergence-conforming B-splines for the unsteady Navier-Stokes equations , 2013, J. Comput. Phys..

[80]  T. Hughes,et al.  Isogeometric analysis of the isothermal Navier-Stokes-Korteweg equations , 2010 .

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

[82]  Hung Nguyen-Xuan,et al.  Static, free vibration, and buckling analysis of laminated composite Reissner–Mindlin plates using NURBS‐based isogeometric approach , 2012 .

[83]  Peter Wriggers,et al.  Three-dimensional mortar-based frictional contact treatment in isogeometric analysis with NURBS , 2012 .

[84]  Frans P. van der Meer,et al.  Mesolevel Modeling of Failure in Composite Laminates: Constitutive, Kinematic and Algorithmic Aspects , 2012 .

[85]  Wei Liu,et al.  Novel higher order mass matrices for isogeometric structural vibration analysis , 2013 .

[86]  M. Ortiz,et al.  FINITE-DEFORMATION IRREVERSIBLE COHESIVE ELEMENTS FOR THREE-DIMENSIONAL CRACK-PROPAGATION ANALYSIS , 1999 .

[87]  Cv Clemens Verhoosel,et al.  A phase-field description of dynamic brittle fracture , 2012 .

[88]  Thomas J. R. Hughes,et al.  Isogeometric Analysis: Toward Integration of CAD and FEA , 2009 .

[89]  Les A. Piegl,et al.  The NURBS Book , 1995, Monographs in Visual Communication.

[90]  Ganesh Subbarayan,et al.  Constructive solid analysis: a hierarchical, geometry-based meshless analysis procedure for integrated design and analysis , 2004, Comput. Aided Des..

[91]  Vinh Phu Nguyen,et al.  Isogeometric analysis: An overview and computer implementation aspects , 2012, Math. Comput. Simul..

[92]  Hung Nguyen-Xuan,et al.  Isogeometric finite element analysis of composite sandwich plates using a higher order shear deformation theory , 2013 .

[93]  Olivier Allix,et al.  A relocalization technique for the multiscale computation of delamination in composite structures , 2011, 1109.4799.

[94]  Demirkan Coker,et al.  Modeling of the dynamic delamination of L-shaped unidirectional laminated composites , 2012 .

[95]  T. Hughes,et al.  A Simple Algorithm for Obtaining Nearly Optimal Quadrature Rules for NURBS-based Isogeometric Analysis , 2012 .

[96]  Howard L. Schreyer,et al.  Modeling delamination as a strong discontinuity with the material point method , 2002 .

[97]  Rakesh K. Kapania,et al.  Geometrically nonlinear NURBS isogeometric finite element analysis of laminated composite plates , 2012 .

[98]  Marc Duflot,et al.  Meshless methods: A review and computer implementation aspects , 2008, Math. Comput. Simul..

[99]  Hung Nguyen-Xuan,et al.  High-order B-splines based finite elements for delamination analysis of laminated composites , 2013 .

[100]  Ted Belytschko,et al.  A finite element method for crack growth without remeshing , 1999 .

[101]  L. J. Sluys,et al.  Mesh-independent modeling of both distributed and discrete matrix cracking in interaction with delamination in composites , 2010 .

[102]  T Belytschko,et al.  Structured Extended Finite Element Methods of Solids Defined by Implicit Surfaces , 2002 .

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

[104]  Vinh Phu Nguyen,et al.  Nitsche’s method for two and three dimensional NURBS patch coupling , 2013, 1308.0802.

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

[106]  Cv Clemens Verhoosel,et al.  A dissipation‐based arc‐length method for robust simulation of brittle and ductile failure , 2009 .

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

[108]  M. Ortiz,et al.  Subdivision surfaces: a new paradigm for thin‐shell finite‐element analysis , 2000 .

[109]  Paolo Ermanni,et al.  Model for interlaminar normal stress in singly curved laminates , 2005 .

[110]  T. Hughes,et al.  Isogeometric fluid-structure interaction: theory, algorithms, and computations , 2008 .

[111]  Glaucio H. Paulino,et al.  Adaptive mesh refinement and coarsening for cohesive zone modeling of dynamic fracture , 2012 .