Determination of fiber/matrix interface debond growth parameters from cyclic loading of single fiber composites

Similar documents
Micro-Crack Initiation and Propagation in Fiber Reinforced Composites. Andrejs Pupurs

Effect of fiber volume fraction on fracture mechanics in continuously reinforced fiber composite materials

Long-term Life Prediction of CFRP Structures Based on MMF/ATM Method

Buckling Behavior of 3D Randomly Oriented CNT Reinforced Nanocomposite Plate

Fatigue Life Prediction in Composite Materials

Modelling of damage in composite materials using interface elements

Xiaoming Mao. Department of Physics and Astronomy, University of Pennsylvania. Collaborators: Tom Lubensky, Ning Xu, Anton Souslov, Andrea Liu

Journal of Solid Mechanics and Materials Engineering

FATIGUE DURABILITY OF CONCRETE EXTERNALLY STRENGTHENED WITH FRP SHEETS

A method to determine relative stroke detection efficiencies from multiplicity distributions

Block failure in connections - including effets of eccentric loads

EVALUATION OF THERMAL CONDUCTIVITY IN PITCH- BASED CARBON FIBER REINFORCED PLASTICS

Strain Rate and Temperature Effects on the Nonlinear Behavior of Woven Composites

Predictions of Storage Modulus of Glass Bead-Filled Low-Density-Polyethylene Composites

Nonmonotonic Networks. a. IRST, I Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I Povo (Trento) Italy

CHAPTER 5 THEORETICAL MODELING OF ELASTIC PROPERTIES

POST-PERFORATION LENGTH AND VELOCITY OF KE PROJECTILES WITH SINGLE OBLIQUE TARGETS

ANALYTICAL INVESTIGATION AND PARAMETRIC STUDY OF LATERAL IMPACT BEHAVIOR OF PRESSURIZED PIPELINES AND INFLUENCE OF INTERNAL PRESSURE

Density. Density 22/10/2015. Rules of Mixture for Elastic Properties. Materi ke-6

Critical Analysis of FRP-NSM double-shear pull-out tests. Abstract

DIRECT NUMERICAL SIMULATION OF DAMAGE PROGRESSION IN LAMINATED COMPOSITE PLATES USING MULTI-SCALE MODELLING

Dynamic buckling of laminated composite stringer stiffened CFRP panels under axial compression

Damage Detection using Stochastic Subspace Identification

Micromechanical models for stiffness prediction of alumina trihydrate (ATH)

Revealing spatially heterogeneous relaxation in a model. nanocomposite

ESTIMATION OF THE VISCOELASTIC PARAMETERS OF LAMINATED COMPOSITES. PART I. ANALYTICAL CONSIDERATIONS

A NEW APPROACH TO DYNAMIC BUCKLING LOAD ESTIMATION FOR PLATE STRUCTURES

Estimating flow properties of porous media with a model for dynamic diffusion

Materials and Design

Benchmarking of lamina failure tests from WWFE-I and WWFE-II with a three parameter micromechanics based matrix failure theory

Fiber / Toughened Epoxy Composites under

Supplement of Detailed characterizations of the new Mines Douai comparative reactivity method instrument via laboratory experiments and modeling

2nd Workshop on Joints Modelling Dartington April 2009 Identification of Nonlinear Bolted Lap Joint Parameters using Force State Mapping

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects

Debond crack growth in fatigue along fiber in UD composite with broken fibers

Dual porosity DRM formulation for flow and transport through fractured porous media

Random Vibration Fatigue Analysis with LS-DYNA

THERMAL ENDURANCE OF UNREINFORCED UNSATURATED POLYESTERS AND VINYL ESTER RESINS = (1) ln = COMPOSITES & POLYCON 2009

Optimization of Flywheel Weight using Genetic Algorithm

Supplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators

Supplementary Materials for

Fatigue Prediction for Composite Materials and Structures

EFFECT OF SHEAR STRESS ON FATIGUE LIFE OF COMPOSITE LAMINATES

An Inverse Interpolation Method Utilizing In-Flight Strain Measurements for Determining Loads and Structural Response of Aerospace Vehicles

MICROMECHANICAL FAILURE ANALYSIS OF UNIDIRECTIONAL FIBER-REINFORCED COMPOSITES UNDER IN-PLANE AND TRANSVERSE SHEAR

Fracture Mechanics of soft Materials: crack tip field of a material that exhibits finite extensibility

Characterization of Internal State Variable for fiber fracture in UD Composite

Extension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels

Spine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area

Numerical predictions of the anisotropic. viscoelastic response of uni-directional fibre composites

Engineering Solid Mechanics

Uniaxial compressive stress strain model for clay brick masonry

Research and Experiments on Electromagnetic Field Induced by Two. Coaxial Solenoid Coils of Axially Mag-lev Driving Needle

IDENTIFICATION OF STABILITY OF CONCRETE TUNNEL LINING USING COUPLED MODELING

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

Excitability of guided waves in composites with PWAS transducers

Successive Model-Updating of the dynamic behaviour of casing bodies on a practical example of an axial piston pump

Interfacial Stress Transfer in a Graphene Monolayer Nanocomposite

REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION

2. A crack which is oblique (Swedish sned ) with respect to the xy coordinate system is to be analysed. TMHL

Physical Chemistry I for Biochemists Chem340. Lecture 26 (3/14/11)

A model study of the effective Young s modulus for randomly distributed short-fiber composites

CONTRIBUTIONS TO ULTRASONIC WELDING OF INTELLIGENT COMPOSITE MATERIALS USED IN VEHICLE MANUFACTURING

A Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair

Chapter 8 Deflection. Structural Mechanics 2 Dept of Architecture

Chapter 8. Lecture Notes Dr. Rakhmad Arief Siregar Kolej Universiti Kejuruteraan Utara Malaysia

Nonlinear Viscoelastic Viscoplastic Material Model Including Stiffness Degradation for Hemp/Lignin Composites

DECONVOLUTION VERSUS CONVOLUTION A COMPARISON FOR MATERIALS WITH CONCENTRATION GRADIENT

Bond strength model for interfaces between nearsurface mounted (NSM) CFRP strips and concrete

Page 1. Physics 131: Lecture 22. Today s Agenda. SHM and Circles. Position

ACTIVE VIBRATION CONTROL FOR STRUCTURE HAVING NON- LINEAR BEHAVIOR UNDER EARTHQUAKE EXCITATION

Assessment of wind-induced structural fatigue based on joint probability density function of wind speed and direction

Topic 5a Introduction to Curve Fitting & Linear Regression

Mathematical Models to Determine Stable Behavior of Complex Systems

Approach to Operational Experimental Estimation of Static Stresses of Elements of Mechanical Structures

Nonlinear Analysis of Reinforced Masonry Shear Walls with ASCE 41

Modeling and Analysis of Thermal Bimorph Using COMSOL

Department of Electronic and Optical Engineering, Ordnance Engineering College, Shijiazhuang, , China

Chapter 3 Micromechanical Analysis of a Lamina Ultimate Strengths of a Unidirectional Lamina

Feature Extraction Techniques

Life Prediction Under Multiaxial Fatigue

I. Understand get a conceptual grasp of the problem

A Finite Element Propagation Model For Extracting Normal Incidence Impedance In Nonprogressive Acoustic Wave Fields

Delamination Prediction of Uni-Directional Composite Laminates using Shell Elements and a Strain Rate Dependent Micro-mechanical Material Model

TEST OF HOMOGENEITY OF PARALLEL SAMPLES FROM LOGNORMAL POPULATIONS WITH UNEQUAL VARIANCES

NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT

Application of Differential Geometry in Magnetic Levitation Control

A DESIGN GUIDE OF DOUBLE-LAYER CELLULAR CLADDINGS FOR BLAST ALLEVIATION

SPECTRUM sensing is a core concept of cognitive radio

S. Srinivasan, Technip Offshore, Inc., Houston, TX

COS 424: Interacting with Data. Written Exercises

Determination of Relaxation and Retardation Spectrum from Modulus of Complex Frequency-Domain Material functions

Analysis and Implementation of a Hardware-in-the-Loop Simulation

A 3D Discrete Damage Modeling Methodology for Abaqus for Fatigue Damage Evaluation in Bolted Composite Joints

Measuring Temperature with a Silicon Diode

Description of the nonlinear viscoelastic and thixotropic properties of human blood by the modified model of Maxwell-Gurevich-Rabinovich

Envelope frequency Response Function Analysis of Mechanical Structures with Uncertain Modal Damping Characteristics

Transfer from Arbitrary Pure State to Target Mixed State for Quantum Systems

Determining a Function for the Damping Coefficient of a laminated Stack

Numerical Modeling of Self-Compacting Mortar Flow Using Discrete Element Method

Transcription:

Deterination o iber/atrix interace debond growth paraeters ro cyclic loading o single iber coposites Andrejs Pupurs 1,, Janis Varna 1, Povl Brøndsted 3, Stergios Goutianos 3 1 Luleå University o Technology, S-97187 Luleå, Sweden Swerea SICOMP, SE-9416 Piteå, Sweden 3 Technical University o Denark, Risø Capus, DK-4000 Roskilde, Denark 6 th International Conerence on Coposites Testing and Model Identiication (CopTest 013) -4 April 013, Aalborg, Denark.

c Fatigue daage in UD coposites c FL FL Failure o coposite siilar to static ailure. Rando iber breaks and ailure in ew load cycles. N=1 Rando iber breaks during irst cycle, N=1. Further daage developent in or o growing iber/atrix interace debond cracks. Ater certain aount o load cycles, aligned atrix cracks appear and propagate as debond cracks. Bridged cracks ollowed by new iber cracks Applied strain level is too low (below atigue liit). N=

3 Objective To ocus on iber/atrix interace debond growth in tensiontension atigue To use racture echanics (energy release rate) or the debond growth analysis To calculate the energy release rate and to deterine its dependence on debond length and to describe it by siple itting unction Using experiental data and the calculated energy release rate to evaluate whether a power law can describe the debond length growth rate in cyclic loading To deterine aterial constants in this power law, i it is applicable

Debond growth in tensiontension atigue l d (N) 4 ax Paris law: da dn B K l d in ax 0 in 0 da increase o the crack surace area dn increase o nuber o applied load cycles K stress intensity actor range B, aterial constants da dn B G II dldn dn B * G II where l dn l d /r B * B r G II G( ax ) G( in )

5 Deterination o power law paraeters dl dn dn B * Proper values o Paris law paraeters and B can only be ound i experiental data o l d =(N) are available. Experiental easureents o debond growth in atigue is coplicated or UD coposites. Single iber (SF) coposite with transparent atrix, is a convenient tool to easure debond crack growth in cyclic loading. Analytical and FEM calculations o G II or SF coposites are necessary to account or the stress state. G II

6 Concentric cylinder odel T F M L L l d iber break debonded iber region

7 Long and short debonds F M Stress state regions d) Plateau stress-state in bonded iber region. Long debonds Debond growth in selsiilar anner Analytical odel or G II calculation l d c) Debond tip region with stress singularity. b) Plateau stress-state in debonded iber region. a) Very coplex stress state at the iber break. Short debonds Strong interaction between stress states in regions a) and c) Nuerical ethods or G II calculation

8 Analytical solution. Energy release rate or long debonds J.A. Nairn, Y.C. Liu, On the use o energy ethods or interpretation o results o single-iber ragentation experients, Coposite Interaces, vol.4, pp. 41-67, 1996. 4 T k k r E G th ech II 1 d Q d E k r r r th d Q d E k 1 r r E E d 1 1 1 1 d E d r 1 1 d E Q r E = 3 [GPa] E = 3.5 [GPa] E = 4 [GPa] = 0.3 0.99747 0.997011 0.996598 1.014151 1.016443 1.018716 k k th

9 Nuerical calculation o energy release rate or short debonds u L syetry axis outer boundary T Energy release rate calculation using Virtual Crack Closure Technique (VCCT) Calculations perored using FEM sotware ANSYS G II (l d ) li dl d 0 1 dl d l d dl d l d u l d dl d l d r d r iber break Syetry b.c.

10 Magniication o energy release rate or short debonds.5 CF, k 1.8 FEM.0 CF, k th BEM k j * GF, k k * 1.4 1.5 GF, k th 1.0 0 10 0 30 1.0 0 5 10 l dn l dn BEM results taken ro: E. Graciani, V. Mantič, F. París, J. Varna, Nuerical analysis o debond propagation in the Single Fibre Fragentation Test, Coposites Science and Technology, vol.69, pp. 514-50, 009.

11 General expression or energy release rate Long debonds: G k k T II E 4 r ech th Short debonds: k k * (l dn )k k th k th * (l dn )k th Fitting the whole range: k 0 k * * ( ldn) kth( ldn) k0 a1 1 a a3l ldn 1; dn ; ( l dn ) 1 l Energy release rate change in one load cycle: l dn dn 30 30 G II (l dn ) r E k 0 (l dn ) 4 ( ech ax ) ( in ech ) ( )T( ax ech in ech )

Experiental easureents 1 1 st iber break Glass iber/epoxy E r [GPa] [-] [1/ C] E [GPa] [-] [1/ C] 70 0.0 4.70 10-6 3.17 0.33 91.00 10-6 Fiber was pre-stressed with 7g o weight suspended at each end o iber First a static load was applied until irst break was initiated in iber ( 1st ) Tension-tension cyclic loading with requency =H, R=0.1 was applied under load control. Measureents o debond length l d were perored ater selected nuber o load cycles. Two dierent values o axial atigue strain level were studied: Saples A, B: ax = 1.76%; Saple C: ax = 1.3%. The strain levels approxiately correspond to 80% and 60% o the stress 1st respectively.

Experiental results. Debond growth in cyclic loading 13 N=1 N=1000 N=5000 80 60 l d () 40 ld1 Saple A ld SapleA ld3 Saple A ld4 Saple A ld1 Saple B ld Saple B ld3 Saple B ld4 Saple B N=10000 0 0 1 10 100 1000 10000 100000 N

dldn dn B * G Deterination o the power law paraeters II r E k0 ( ldn) ax in ax in GII ( ldn) ( ech ) ( ech ) ( ) T ( ech ech ) 4 r 0 ( ) 0 E k ldn ax in ax in GII g ( ech) ( ech) ( ) T ( ech ech) 4 14 ln( l dn /N ) - -4-6 -8-10 Saple A Saple B y = 8.799x - 353.69 y =.938x - 86.96-1 11.9 1.0 1.1 1. 1.3 dl ln dn dn dldn dn ln B Saple B * * G 0 II g ln G 0 II ln g ln B * A.94-99.7 B 8.80-118.60 lng II 0

Validation o power law paraeters 15 10 Saple A experiental Saple A odeling 8 Saple B experiantal Saple B odeling Saple C experiental 6 Saple C odeling l dn 4 0 0 0000 40000 60000 80000 100000 N

Conclusions 16 The strain energy release rate in Mode II or debond growth was analyed cobining analytical solution or long debonds and FEM solution or short debonds. Using the quantiied debond length versus nuber o cycles data it was shown that the power law with respect to the strain energy release rate change is applicable or debond growth characteriation in tension-tension atigue. Siulations showed that the obtained paraeters give acceptable predictions or cases, when the debond grows as well as when it does not grow. Being aterial properties the power law paraeters deterined in this study using single iber coposites can also be applied or the case o UD coposites ade o the sae aterial syste

17 Related papers A. Pupurs, S. Goutianos, P. Brøndsted, J. Varna, Interace debond crack growth in tension-tension cyclic loading o single iber polyer coposites, Coposites Part A, vol.44, pp. 86-94, 013. A. Pupurs, J. Varna, FEM odeling o iber/atrix debond growth in tension-tension cyclic loading o unidirectional coposites, International Journal o Daage Mechanics, In Press, 013. (available online)

Thank you or your attention!