FEM predictions of damage in continous fiber ceramic matrix composites under transverse tension using the crack band method

Abstract This paper develops a numerical micromechanics fracture model based on the crack band methodology for the prediction of damage in continuous fiber ceramic matrix composites (CMCs). The energy based crack band formulation offers mesh objective predictions in finite element simulations. Detailed finite element models of CMCs are created based on the actual microstructure and compared with randomly created representative volume elements (RVE). The influence of geometrical inhomogeneities, such as fiber clustering and variations of the fiber coating thickness, on the stress-strain response and the crack development are investigated. Furthermore, FEM models of multi-layer cross-ply laminates are created including thousands of explicitly modeled fibers and fiber coatings. Crack density predictions are compared with well established analytical models by Laws and Dvorak [1] and Nairn [2]. An eight layer FEM model is used to predict damage propagation within each layer. The Numerical prediction of the stress-strain response is compared with experimental results of smooth bar tensile coupons.

[1]  George J. Dvorak,et al.  Progressive Transverse Cracking In Composite Laminates , 1988 .

[2]  John A. Nairn,et al.  The Strain Energy Release Rate of Composite Microcracking: A Variational Approach , 1989 .

[3]  Y. D. Song,et al.  Modeling the Tensile Behavior of Unidirectional C/SiC Ceramic-Matrix Composites , 2014, Mechanics of Composite Materials.

[4]  G. Ravichandran,et al.  Computational modeling of damage evolution in unidirectional fiber reinforced ceramic matrix composites , 1997 .

[5]  S. Spearing,et al.  Matrix crack spacing in brittle matrix composites , 1992 .

[6]  Javier Segurado,et al.  Three dimensional (3D) microstructure-based modeling of interfacial decohesion in particle reinforced metal matrix composites , 2012 .

[7]  A. A. Gusev Representative volume element size for elastic composites: A numerical study , 1997 .

[8]  Samir N. Shoukry,et al.  Multi-Fiber Unit Cell for Prediction of Residual Stresses in Continuous Fiber Composites , 2007 .

[9]  Somnath Ghosh,et al.  Homogenization Based 3D Continuum Damage Mechanics Model for Composites Undergoing Microstructural Debonding , 2008 .

[10]  Mira Mitra,et al.  Modelling matrix damage and fibre-matrix interfacial decohesion in composite laminates via a multi-fibre multi-layer representative volume element (M2RVE) , 2014 .

[11]  D. Jeulin,et al.  Determination of the size of the representative volume element for random composites: statistical and numerical approach , 2003 .

[12]  Julián Bravo-Castillero,et al.  A comprehensive numerical homogenisation technique for calculating effective coefficients of uniaxial piezoelectric fibre composites , 2005 .

[13]  W. Drugan,et al.  A micromechanics-based nonlocal constitutive equation and estimates of representative volume element size for elastic composites , 1996 .

[14]  C. V. Singh,et al.  Damage and Failure of Composite Materials: Author index , 2012 .

[15]  Z. Bažant,et al.  Crack band theory for fracture of concrete , 1983 .

[16]  Residual stresses in random elastic composites: nonlocal micromechanics-based models and first estimates of the representative volume element size , 2013 .

[17]  M. Mahmoodi,et al.  Damage analysis of fiber reinforced Ti-alloy subjected to multi-axial loading—A micromechanical approach , 2011 .