Numerical modeling of two-dimensional delamination growth in composite laminates with in-plane isotropy

Published on May 13, 2021in Engineering Fracture Mechanics3.426
· DOI :10.1016/J.ENGFRACMECH.2021.107787
Congzhe Wang (EPFL: École Polytechnique Fédérale de Lausanne), Anastasios P. Vassilopoulos28
Estimated H-index: 28
(EPFL: École Polytechnique Fédérale de Lausanne),
Thomas Keller43
Estimated H-index: 43
(EPFL: École Polytechnique Fédérale de Lausanne)
Sources
Abstract
Abstract The two-dimensional (2D) delamination growth in fiber-reinforced polymer (FRP) laminates with in-plane isotropy under Mode I loading condition was numerically investigated using finite element analyses. Two sizes of plate models were developed, focusing on different fracture stages. Cohesive elements were employed to simulate the fracture behavior in the presence of large-scale bridging (LSB). The influences of the pre-crack shape/area, loading zone shape/area and fracture resistance were parametrically studied. It was found that either a flatter pre-crack shape or a flatter loading zone shape could result in higher initial structural stiffness and less uniform distribution of the strain energy release rate (SERR) along the pre-crack perimeter during crack initiation and early propagation. However, they had only a minor effect on the stiffness after full fiber bridging development in all directions. The plates finally achieved constant stiffness, which increased linearly with the fracture resistance. The final crack shape was dependent on the loading zone shape and area, but the effects were relatively weak.
📖 Papers frequently viewed together
6 Citations
2 Citations
References33
Newest
#1G. Pappas (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 8
#2John Botsis (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 37
Abstract Traction-separation relations due to large scale bridging in composites are very important in modeling their fracture response. Several works demonstrated that such relations are dependent upon the specimen's stiffness and loading conditions. Experimental data and micromechanics of bridging suggested that, in addition to the crack opening displacement (COD), the local curvature can be used as an additional kinematic parameter to incorporate the stiffness dependence on traction-separatio...
2 CitationsSource
#1Zhengwen Jiang (Hunan University)H-Index: 2
#2Shui Wan (SEU: Southeast University)H-Index: 2
Last. Anastasios P. Vassilopoulos (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 28
view all 5 authors...
Abstract This work investigated the evolution of the non-uniform distribution of Mode-I strain energy release rate (SERR) along a straight pre-crack front of double cantilever beams (DCB), caused by the anticlastic deformation and the longitudinal free-edge effects. A modified DCB specimen with a curved pre-crack front was introduced in order to alleviate this problem and produce an uniform SERR distribution at the initiation and propagation stages. The first order optimization method was used t...
2 CitationsSource
#1Aida Cameselle-Molares (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 4
#2Anastasios P. Vassilopoulos (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 28
Last. Thomas Keller (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 43
view all 5 authors...
Abstract A new numerically-based method suitable for determining the total strain energy release rate (SERR) involved in two-dimensional (2D) Mode I-dominated delamination under opening loads in FRP laminates is presented. The method is based on the mutual dependence of the load vs opening displacement curves slope exhibited after the full development of the fracture process zone (FPZ) and the total SERR involved in the delamination process. The equation relating these parameters is derived from...
1 CitationsSource
#1F. Naya (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 1
#1F. Naya (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 9
Last. John Botsis (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 37
view all 3 authors...
Abstract Fiber reinforced polymers (FRPs) subjected to mode I fracture show important toughening due to the development of large scale bridging (LSB). Experimental studies of this phenomenon in unidirectional carbon/epoxy laminates using double cantilever beam specimens, demonstrate important differences in R-curve response for inter- and intralaminar fracture. Post fracture observation of composite’s cross-section pointed out dissimilar fiber bundle size and shape, as the main origin of their d...
15 CitationsSource
#1Aida Cameselle-Molares (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 4
#2Anastasios P. Vassilopoulos (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 28
Last. Thomas Keller (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 43
view all 3 authors...
Abstract The 2D delamination behavior of composite laminates under quasi-static out-of-plane opening loading has been experimentally investigated. A circular embedded pre-crack was introduced in the center and at the midplane of the laminates. To investigate the effect of the fiber architecture on the fracture behavior, six different types of fabrics were selected. Increasing load–displacement curves were obtained due to the increasing crack front length during propagation. During the loading pr...
6 CitationsSource
#1Aida Cameselle-Molares (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 4
#2Anastasios P. Vassilopoulos (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 28
Last. Thomas Keller (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 43
view all 5 authors...
Abstract A numerical investigation was carried out in order to simulate experimental results previously obtained concerning two-dimensional (2D) in-plane crack propagation in laminated glass fiber-reinforced polymer (GFRP) plates. The laminated plates were designed with an embedded circular pre-crack and subjected to quasi-static out-of-plane loads. In order to study the transition from standard fracture mechanics tests, where the crack propagates only in one dimension (1D), to 2D scenarios, add...
6 CitationsSource
#1K. Senthil (VSSC: Vikram Sarabhai Space Centre)H-Index: 4
#2A. Arockiarajan (Indian Institute of Technology Madras)H-Index: 16
Last. R. Palaninathan (Indian Institute of Technology Madras)H-Index: 10
view all 3 authors...
Abstract Adhesively bonded composite joints plays, an important role in the modern composite structures assembly to realize light weight structures than the mechanical fastening. However, in the process of adhesive bonding of composite components can result in manufacturing flaws namely debonds. Debonds can change the structural response of adhesively bonded joints significantly in compressively loaded aerospace structures. In general, the prediction of the maximum load that composite structures...
5 CitationsSource
#1Andrew Rhead (University of Bath)H-Index: 10
#2Richard Butler (University of Bath)H-Index: 21
Last. Giles W Hunt (University of Bath)H-Index: 32
view all 3 authors...
Compression After Impact (CAI) strength is critical to the safety and weight of carbon fibre aircraft. In this paper, the standard aerospace industry practice of using separate analyses and tests for panel buckling and CAI strength is challenged. Composite panels with a range of stacking sequences were artificially delaminated and subject to compression testing in a fixture that allowed local sublaminate and global panel buckling modes to interact. Compared to panels without delamination, intera...
24 CitationsSource
#1Luis P. Canal (École Polytechnique)H-Index: 11
#2Marco Alfano (University of Calabria)H-Index: 22
Last. John Botsis (École Polytechnique)H-Index: 37
view all 3 authors...
Abstract A cohesive zone model (CZM) is proposed to assess the thickness scaling effect associated with fiber bridging during fracture. The CZM was developed through a multi-scale simulation approach and utilizes an embedded cell model of the Double Cantilever Beam (DCB) that explicitly accounts for the bridging bundles on the fracture plane. In particular, micromechanical simulations of failure were carried out, for varying arms thickness, in order to determine the homogenized fracture behavior...
26 CitationsSource
#1Guillaume Frossard (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 3
#2Joël Cugnoni (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 26
Last. John Botsis (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 37
view all 4 authors...
Abstract In this work, the influence of ply thickness on strain energy release rate (ERR) in delamination of carbon-epoxy laminates is addressed. Specimens with three ply thicknesses: 0.030, 0.075 and 0.150 mm are tested. While the ERR at onset of crack propagation is independent of ply thickness, the plateau level is much lower in thin-ply laminates than in thicker ply ones. This effect is attributed to changes in microstructure, caused by the tow spreading process involved in the fabrication o...
27 CitationsSource
Cited By0
Newest