Tensile behaviour of anti-symmetric CFRP composite

Similar documents
THE ROLE OF DELAMINATION IN NOTCHED AND UNNOTCHED TENSILE STRENGTH

FRACTURE TOUGHNESS OF ADHESIVE BONDED COMPOSITE JOINTS UNDER MIXED MODE LOADING.

QUESTION BANK Composite Materials

Capability Assessment of Finite Element Software in Predicting the Last Ply Failure of Composite Laminates

University of Bristol - Explore Bristol Research. Early version, also known as pre-print

Prediction of Delamination Growth Behavior in a Carbon Fiber Composite Laminate Subjected to Constant Amplitude Compression-Compression Fatigue Loads

NUMERICAL INVESTIGATION OF DELAMINATION IN L-SHAPED CROSS-PLY COMPOSITE BRACKET

COMPARISON OF COHESIVE ZONE MODELS USED TO PREDICT DELAMINATION INITIATED FROM FREE-EDGES : VALIDATION AGAINST EXPERIMENTAL RESULTS

BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS

Open-hole compressive strength prediction of CFRP composite laminates

Strength of GRP-laminates with multiple fragment damages

Calculation of Energy Release Rate in Mode I Delamination of Angle Ply Laminated Composites

Finite element modelling of infinitely wide Angle-ply FRP. laminates

INITIATION AND PROPAGATION OF FIBER FAILURE IN COMPOSITE LAMINATES

THE MUTUAL EFFECTS OF SHEAR AND TRANSVERSE DAMAGE IN POLYMERIC COMPOSITES

INTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 1, 2011

EXPERIMENTAL CHARACTERIZATION AND COHESIVE LAWS FOR DELAMINATION OF OFF-AXIS GFRP LAMINATES

Crack Length Dependence of Mode III Delamination Using Edge Crack Torsion Test

FASTENER PULL-THROUGH FAILURE IN GFRP LAMINATES

Numerical simulation of delamination onset and growth in laminated composites

Modelling the nonlinear shear stress-strain response of glass fibrereinforced composites. Part II: Model development and finite element simulations

NUMERICAL AND EXPERIMENTAL ANALYSES OF MULTIPLE DELAMINATIONS IN CURVED COMPOSITE LAMINATES

Discrete ply modelling of open hole tensile tests

Interlaminar fracture characterization in composite materials by using acoustic emission

FLOATING NODE METHOD AND VIRTUAL CRACK CLOSURE TECHNIQUE FOR MODELING MATRIX CRACKING- DELAMINATION MIGRATION

Autodesk Helius PFA. Guidelines for Determining Finite Element Cohesive Material Parameters

Modeling of Interfacial Debonding Induced by IC Crack for Concrete Beam-bonded with CFRP

Impact and Crash Modeling of Composite Structures: A Challenge for Damage Mechanics

PREDICTION OF OUT-OF-PLANE FAILURE MODES IN CFRP

PERFORMANCE OF COMPOSITE PANELS SUBJECTED TO UNDERWATER IMPULSIVE LOADING

Failure analysis of serial pinned joints in composite materials

Published in: Composites Part B: Engineering. Document Version: Peer reviewed version

Module III - Macro-mechanics of Lamina. Lecture 23. Macro-Mechanics of Lamina

DAMAGE SIMULATION OF CFRP LAMINATES UNDER HIGH VELOCITY PROJECTILE IMPACT

PRELIMINARY PREDICTION OF SPECIMEN PROPERTIES CLT and 1 st order FEM analyses

Passive Damping Characteristics of Carbon Epoxy Composite Plates

ID-1160 REAL-TIME DETECTION AND EXPLICIT FINITE ELEMENT SIMULATION OF DELAMINATION IN COMPOSITE LAMINATES UNDER IMPACT LOADING

Fracture Mechanics, Damage and Fatigue: Composites

This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International licence

BEARING STRENGTH ASSESSMENT OF COMPOSITE MATERIAL THROUGH NUMERICAL MODELS

Crash and Impact Simulation of Composite Structures by Using CAE Process Chain

SKIN-STRINGER DEBONDING AND DELAMINATION ANALYSIS IN COMPOSITE STIFFENED SHELLS

IMPACT ON LAMINATED COMPOSITE PLATES: COMPARISON OF TEST AND SIMULATION RESULTS OBTAINED WITH LMS SAMTECH SAMCEF

KINK BAND FORMATION OF FIBER REINFORCED POLYMER (FRP)

NON-LINEAR FRACTURE BEHAVIOR OF DOUBLE CANTILEVER BEAM

An investigation of the mechanical behaviour of carbon epoxy cross ply cruciform specimens under biaxial loading

Prediction of failure in notched CFRP laminates under multi-axial loading. J L Y Tan, V S Deshpande and N A Fleck*,

Finite element analysis of drilled holes in uni-directional composite laminates using failure theories

DYNAMIC DELAMINATION OF AERONAUTIC STRUCTURAL COMPOSITES BY USING COHESIVE FINITE ELEMENTS

Dynamic analysis of Composite Micro Air Vehicles

Failure Mode Predictions in the Compressive Response of Laminated Composites

ISSN: ISO 9001:2008 Certified International Journal of Engineering Science and Innovative Technology (IJESIT) Volume 2, Issue 4, July 2013

A RESEARCH ON NONLINEAR STABILITY AND FAILURE OF THIN- WALLED COMPOSITE COLUMNS WITH OPEN CROSS-SECTION

Effects of Resin and Fabric Structure

Non-conventional Glass fiber NCF composites with thermoset and thermoplastic matrices. F Talence, France Le Cheylard, France

TABLE OF CONTENTS. Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA

Keywords: CFRP, compressive failure, kink-band, cohesive zone model. * Corresponding author

Progressive Failure Analysis of Laminated Composite Structures

NUMERICAL MODELLING OF COMPOSITE PIN- JOINTS AND EXPERIMENTAL VALIDATION

CHARACTERIZATION, ANALYSIS AND PREDICTION OF DELAMINATION IN COMPOSITES USING FRACTURE MECHANICS

ADVANCES IN THE PROGRESSIVE DAMAGE ANALYSIS OF COMPOSITES

Probabilistic fatigue life prediction of multidirectional composite laminates

SSRG International Journal of Mechanical Engineering (SSRG-IJME) volume1 issue5 September 2014

Effect of Specimen Dimensions on Flexural Modulus in a 3-Point Bending Test

PROGRESSIVE DAMAGE ANALYSES OF SKIN/STRINGER DEBONDING. C. G. Dávila, P. P. Camanho, and M. F. de Moura

A FINITE ELEMENT MODEL FOR THE ANALYSIS OF DELAMINATIONS IN FRP SHELLS

Prediction of Elastic Constants on 3D Four-directional Braided

Progressive Damage of GFRP Composite Plate Under Ballistic Impact: Experimental and Numerical Study

FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS

Fracture mechanics fundamentals. Stress at a notch Stress at a crack Stress intensity factors Fracture mechanics based design

Assessment Methods of Mechanical Properties of Composite Materials

Numerical Analysis of Delamination Behavior in Laminated Composite with Double Delaminations Embedded in Different Depth Positions

NUMERICAL FEM ANALYSIS FOR THE PART OF COMPOSITE HELICOPTER ROTOR BLADE

*Corresponding author: Keywords: Finite-element analysis; Multiscale modelling; Onset theory; Dilatational strain invariant.

Failure Analysis of Unidirectional Composite Pinned- Joints

Analysis of nonlinear shear deformations in CFRP and GFRP textile laminates

ARTICLE IN PRESS Materials Science and Engineering A xxx (2009) xxx xxx

Project MMS13 Task 5 Report No 3 (M6/D3)

Finite Element Analysis of FRP Debonding Failure at the Tip of Flexural/Shear Crack in Concrete Beam

Modeling Hailstone Impact onto Composite Material Panel Under a Multi-axial State of Stress

IJSER 1. INTRODUCTION. M.Elhadary

Calculation of Damage-dependent Directional Failure Indices from the Tsai-Wu Static Failure Criterion

Experimental characterization of interlaminar fracture toughness of composite laminates assembled with three different carbon fiber lamina

EXPERIMENTAL AND NUMERICAL INVESTIGATION ON THE FAILURE MODES OF THICK COMPOSITE LAMINATES

Crashworthy Design of Composite Structures Using CAE Process Chain

Fracture Behaviour of FRP Cross-Ply Laminate With Embedded Delamination Subjected To Transverse Load

Computational Analysis for Composites

Material Characterization of Carbon Fiber Reinforced Polymer Laminate Using Virtual Fields Method

FRACTURE MECHANICS OF COMPOSITES WITH RESIDUAL STRESSES, TRACTION-LOADED CRACKS, AND IMPERFECT INTERFACES

The Accuracy of Characteristic Length Method on Failure Load Prediction of Composite Pinned Joints

MODELING DYNAMIC FRACTURE AND DAMAGE IN A FIBER-REINFORCED COMPOSITE LAMINA WITH PERIDYNAMICS

REPRESENTING MATRIX CRACKS THROUGH DECOMPOSITION OF THE DEFORMATION GRADIENT TENSOR IN CONTINUUM DAMAGE MECHANICS METHODS

Characterization of Fiber Bridging in Mode II Fracture Growth of Laminated Composite Materials

SIMULATION OF PROGRESSIVE FAILURE PREDICTION OF FILAMENT WOUND COMPOSITE TUBES SUBJECTED TO MULTIPLE LOADING WITH MEMBRANE-FLEXION COUPLING EFFECTS

Experimental methods to determine model parameters for failure modes of CFRP

TENSILE FATIGUE BEHAVIOR OF SINGLE FIBRES AND FIBRE BUNDLES

PLY LEVEL UNCERTAINTY EFFECTS ON FAILURE OF COMPOSITE

SOME RESEARCH ON FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS

Laboratory 4 Bending Test of Materials

Micromechanical analysis of FRP hybrid composite lamina for in-plane transverse loading

Transcription:

Available online at www.sciencedirect.com Procedia Engineering 1 (211) 1865 187 ICM11 Tensile behaviour of anti-symmetric CFRP composite K. J. Wong a,b, *, X. J. Gong a, S. Aivazzadeh a, M. N. Tamin b a Département de Recherche en Ingénierie des Véhicules pour l Environnement, Université de Bourgogne, 58 Nevers, France b Centre for Composites, Faculty of Mechanical Engineering, Universiti Teknologi Malaysia, 8131 Johor, Malaysia Elsevier use only: Received date here; revised date here; accepted date here Abstract This paper addresses the response of a 17-ply anti-symmetric carbon/epoxy composite subjected to uniaxial tensile loading. Hashin ply damage model is adopted to describe the damage behaviour of the plies, whereas damage initiation and progression of the interfaces are characterised by mixed-mode cohesive damage model. Forcedisplacement curves obtained numerically and experimentally show good agreement. Results show that all laminae and interfaces experience the damage except laminae with o fibre. In addition, damage is concentrated at the tab and central regions of the tensile specimen. Edge delamination is observed in all interfaces. 211 Published by Elsevier Ltd. Open access under CC BY-NC-ND license. Selection and peer-review under responsibility of ICM11 Keywords: CFRP composite; Tensile behaviour; Lamina damage; Interface failure; Finite Element Method (FEM) 1. Introduction Carbon fibre reinforced polymer (CFRP) composites are widely employed in advanced structural applications. In order to be used safely as structural materials, thorough understanding on the response of the composites under different loading conditions is inevitable. Studies on glass fibre reinforced polymer (GFRP) composites with central cut-out were reported by Rakesh and his coworkers [1]. Hashin model [2,3] was adopted to predict the damage initiation of the lamina. Besides, Lapczyk and Hurtado studied blunt notch fibre metal laminate (FML) behaviour numerically by considering both lamina and interface failure [4]. Through both experimental and numerical works, it was reported that matrix cracking and delamination are the two dominant failure modes in quasi-isotropic (QI) carbon/epoxy composite laminates [5]. From the works mentioned above, it * Corresponding author. Tel.: +33 3 86 71 5 36; fax: +33 3 86 71 5 18. E-mail address: jye1189@gmail.com. 1877-758 211 Published by Elsevier Ltd. Open access under CC BY-NC-ND license. Selection and peer-review under responsibility of ICM11 doi:1.116/j.proeng.211.4.31

1866 K. J. Wong et al. / Procedia Engineering 1 (211) 1865 187 is recommended that both lamina and interface failures should be modelled in order to accurately characterise the damage behaviour of the composites. In this study, the response of a 17-layer anti-symmetric CFRP composite arranged in the ply sequence of [-6 2 //6 5 //-6 5 //6 2 ] subjected to tensile load is examined. This particular layup could eliminate coupling stiffness matrix [B] and the entries A 16, A 26, D 16 and D 26 in extensional [A] and bending [D] stiffness matrices respectively. Besides, for other entries, A ij is proportional to D ij. Nomenclature Lamina G ft,c, G fc,c G mt,c, G mc,c S L, X T, Y T ˆij Interface D G C G j, G jc t n, t s, t t t u,n, t u,s, t u,t o, m, f fibre critical fracture energy, t and c refer to tensile and compressive respectively matrix critical fracture energy, t and c refer to tensile and compressive respectively longitudinal shear, longitudinal and transverse strength coefficient of shear stress contribution effective stress components on ijth-plane damage parameter effective strain energy release rate mode j strain energy release rate, C refers to critical strain energy release rate nominal traction stress vector in normal, shear and tangent directions interface strength in normal, shear and tangent directions interface relative displacement at damage initiation, intermediate and total failure material parameter 2. Numerical models 2.1. Lamina damage model The different failure modes of Hashin damage failure criterion [2,3] are treated separately and the tensile failure conditions are described by equations (1)-(2) as follow: 2 2 ˆ ˆ 11 12 Fibre tension, F ft T L 1 X S 2 2 ˆ ˆ 22 12 Matrix tension, F mt T L 1 Y S (1) (2) Damage is initiated for any of the failure modes whenever the criterion given by the equations is satisfied. Damage evolution is described by linear softening law, which is to be discussed in the next section.

K. J. Wong et al. / Procedia Engineering 1 (211) 1865 187 1867 2.2. Interface damage model A quadratic nominal stress criterion is used to indicate the initiation of interface damage, as shown in equation (3). 2 2 2 t n t s t t 1 tun, tus, tut, (3) Following interface damage initiation event, subsequent damage propagation upon further loading is predicted using a mixed-mode energy-based criterion given by equation (4) [6,7]. GII G III G ( G G ) G GI GII GIII IC IIC IC C (4) In the current work, linear softening law is used and the damage variable is described as follows: f ( m o) D ( ) m f o (5) In the lamina damage model, D refers to fibre and matrix damage variables respectively. 2.3. Finite element model The finite element model of the composite specimen consists of seventeen laminae and six interfaces, as shown in Fig. 1. Besides, four cross-ply GFRP tabs are modelled. Both laminae and tabs are modelled using continuum shell elements and the interfaces are modelled using cohesive elements. Fig. 2 shows the finite element model used in this study along with the loading and boundary conditions used. L17 I6 I1 L1 6 2-6 5 6 5-6 2 lamina interface Fig. 1. Layup and designation of the composite used in this study Fig. 2. The finite element model used in this study The properties of the lamina, interface and GFRP tab are shown in Tables 1 thru 6 below. The parameter, is chosen to be unity, as proposed by Hashin [3]. Table 1. Material properties of the lamina E 1 (GPa) E 2 (GPa) E 3 (GPa) G 12 (GPa) G 13 (GPa) G 23 (GPa) 12 13 23 131.9 9.51 9.43 5.27 7.3 3.39.326.341.485

1868 K. J. Wong et al. / Procedia Engineering 1 (211) 1865 187 Table 2. Strengths of the lamina X T (MPa) X C (MPa) Y T (MPa) X C (MPa) S L (MPa) S T (MPa) 1328 164 7.9 221 71.2 94.5 Table 3. Fracture energies of the lamina G mt,c (N/mm) G mc,c (N/mm) G ft,c (N/mm) G fc,c (N/mm).33.33 2 2 Table 4. Material properties of the interface K nn (MPa) K ss (MPa) K tt (MPa) 4 14 14 Table 5. Damage properties of the interface t u,n (MPa) t u,s (MPa) t u,t (MPa) G IC (N/mm) G IIC =G IIIC (N/mm) 25 13.5 13.5.33.79 1.45 Table 6. Material properties of the tab E 1 (GPa) E 2 (GPa) G 12 = G 13 (GPa) G 23 (GPa) 12 55 9.5 5.5 3.33 3. Results and discussion 3.1. Load-displacement curves The measured and predicted load-displacement responses of the composite specimen under tensile loading are compared in Fig. 3. The strain gauges readout is used to compute the relative displacement over the gauge length of 3 mm. The linear plot suggests that no apparent damage is observed during loading until the final catastrophic fracture occurs. The deviation in the measured and predicted stiffness could be attributed to imperfect composite plate fabrication. The peak load at fracture is compared well to predicted load magnitude. The difference could be due to inability of the numerical computation to account for the load carrying ability of the fibre after complete failure of the matrix. Force, P (kn) 1 8 6 4 P = 398.5 P = 366.82 R² =.9985 P = 351.7 R² =.9988 2 Gauge-1 Gauge-2 Num.5.1.15.2.25.3 Gauge relative displacement, (mm) Fig. 3. Experimental and predicted load-displacement curves

K. J. Wong et al. / Procedia Engineering 1 (211) 1865 187 1869 3.2 Damage behaviour of laminae and interfaces Evolution of internal states and damage in the laminae and interfaces throughout the loading, as predicted by the finite element model is now presented and discussed. Fig. 4 describes the corresponding time at damage initiation of each lamina and interface. Damage is initiated in all laminae and interfaces except laminae with o fibre. It is predicted that only tensile matrix failure occurs in the lamina during the tension test. This matrix damage occurs due to the combined effect of 22 and 12, that leads to F mt = 1 (see equation (2)). Interlaminar failure is a result of combined stresses in 33 and 13-direction. Time (s) 8 7 6 5 4 3 2 1 Time at peak load L1 L2 I1 L3 I2 L4 L5 L6 L7 L8 I3 L9 I4 L1 L11 L12 L13 L14 I5 L15 I6 L16 L17 Lamina/interface Fig. 4. Damage initiation time of damaged laminae and interfaces Next, the progression of damage in L1 is described. Fig. 5 illustrates the sequence of matrix damage propagation with continuously increasing tensile loading of the specimen. Results show damage is initiated at near tab region (Fig. 5(a)). Then, the stresses and thus matrix damage variable are increased uniformly across the whole lamina (Fig. 5. (b)). This is followed by matrix damage initiation at isolated region nearby (Fig. 5. (c)), that propagated to the edges (Fig. 5. (d)). Fig. 5. (e) illustrates the damage contour at peak load, where matrix cracking at the central region is observed. Besides, anti-symmetric feature of the lamina is obvious as reflected in anti-symmetric distribution of matrix cracking damage. Fig. 5. Progression of damage initiation in L1 Fig. 6 describes the initiation and progression of delamination damage in I1. Since damage is initiated at the edges near the tabs (circled in red in Fig. 6. (a)), delamination propagates across the width and connects after that (Fig. 6. (b)). It is accompanied by edge delamination along the edges. Then, delamination is predicted to propagate slightly further from the tabs (Fig. 6. (c)). When peak load is attained, delamination is found at the central edges of the specimen (Fig. 6. (d)).

187 K. J. Wong et al. / Procedia Engineering 1 (211) 1865 187 Edge delamination Fig. 6. Interface damage initiation and propagation in I1 4. Conclusions Failure process of anti-symmetric CFRP composite laminate specimen under tensile loading has been examined both experimentally and using finite element method. Both lamina damage and interface delamination were modelled. Results show that: 1. Numerical prediction and experimental measurement of load-displacement curves are in good agreement. 2. Damage occurs in two forms, which are matrix cracking (in laminae) and delamination (in interfaces). 3. For both lamina and interface, damage is first initiated at the tab region, which eventually propagated slightly and accompanied by damage at the central region. As for interfaces, delamination occurs along the edges as well. References [1] Rakesh PK, Singh I, Kumar D. Failure prediction in glass fiber reinforced plastics laminates with drilled hole under uniaxial loading Mater Des 21;31(6):32-7. [2] Hashin Z, Rotem A. A fatigue failure criterion for fiber reinforced materials. J Compos Mater 1973;7:448 64. [3] Hashin Z. Failure criteria for unidirectional fiber composites. J Appl Mech 198:47(2):329-34. [4] Lapczyk I, Hurtado JA. Progressive damage modeling in fiber-reinforced materials. Compos Part A-Appl S 27;38(11):2333-41. [5] Hallett SR, Jiang WG, Khan B, Wisnom MR. Modelling the interaction between matrix cracks and delamination damage in scaled quasi-isotropic specimens. Compos Sci Technol 28;68(1):8-9. [6] Gong XJ, Benzeggagh ML. Mixed mode interlaminar fracture toughness of unidirectional glass/epoxy composite. In: Martin RH, editor. Composite materials: Fatigue and Fracture - Fifth Volume, ASTM STP 123, Philadelphia: American Society for Testing and Materials; 1995, p. 1-23. [7] Benzeggagh ML, Kenane M. Measurement of mixed-mode delamination fracture toughness of unidirectional glass/epoxy composites with mixed-mode bending apparatus. Compos Sci Technol 1996;66:1479-96.