Experimental setup for and numerical modelling of bending fatigue experiments on plain woven glass/epoxy composites

Similar documents
Poisson s ratio as a sensitive indicator of (fatigue) damage in fibre-reinforced plastics

Nonlinear correction to the bending stiffness of a damaged composite beam

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

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

International Journal of Fatigue

BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS

IJSER 1. INTRODUCTION. M.Elhadary

A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE

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

DEVELOPMENT OF THERMOELASTIC STRESS ANALYSIS AS A NON-DESTRUCTIVE EVALUATION TOOL

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

MESO-SCALE MODELLING IN THERMOPLASTIC 5-HARNESS SATIN WEAVE COMPOSITE

INVESTIGATION OF THE PROCESSING PARAMETERS OF A 3D WOVEN REINFORCEMENT

IDENTIFICATION OF THE ELASTIC PROPERTIES ON COMPOSITE MATERIALS AS A FUNCTION OF TEMPERATURE

In Situ Ultrasonic NDT of Fracture and Fatigue in Composites

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

Multi Disciplinary Delamination Studies In Frp Composites Using 3d Finite Element Analysis Mohan Rentala

EFFECT OF THERMAL FATIGUE ON INTRALAMINAR CRACKING IN LAMINATES LOADED IN TENSION

TENSILE FATIGUE BEHAVIOR OF SINGLE FIBRES AND FIBRE BUNDLES

2 Experiment of GFRP bolt

Laboratory 4 Bending Test of Materials

Open-hole compressive strength prediction of CFRP composite laminates

SCALING EFFECTS IN THE LOW VELOCITY IMPACT RESPONSE OF FIBRE METAL

QUESTION BANK Composite Materials

Computational Analysis for Composites

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

Influence of residual stresses in the structural behavior of. tubular columns and arches. Nuno Rocha Cima Gomes

EXPERIMENTAL AND THEORETICAL STUDY OF THE DAMAGE ONSET IN BIAXIAL CRUCIFORM SPECIMENS UNDER STATIC AND HYSTERESIS LOADING

Flexural properties of polymers

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

NUMERICAL MODELLING OF COMPOSITE PIN- JOINTS AND EXPERIMENTAL VALIDATION

THREE DIMENSIONAL STRESS ANALYSIS OF THE T BOLT JOINT

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

A fatigue design methodology for GRP composites in offshore underwater applications

Probabilistic fatigue life prediction of multidirectional composite laminates

Strength of GRP-laminates with multiple fragment damages

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

Fig. 1. Different locus of failure and crack trajectories observed in mode I testing of adhesively bonded double cantilever beam (DCB) specimens.

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

Fatigue Damage Modelling of Fibre-reinforced Composite Materials: Review

Meso-scale Modeling of Static and Fatigue Damage in Woven Composite Materials with Finite Element Method

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

MODELING OF THE BEHAVIOR OF WOVEN LAMINATED COMPOSITES UNTIL RUPTURE

Composites Science and Technology

Calibration and Experimental Validation of LS-DYNA Composite Material Models by Multi Objective Optimization Techniques

INTERNAL STRAIN MEASUREMENTS IN CFRP PLATES SUBJECTED TO IMPACT LOAD USING FBG SENSORS

NUMERICAL FEM ANALYSIS FOR THE PART OF COMPOSITE HELICOPTER ROTOR BLADE

DESIGNING A FLEXIBLE BELLOWS COUPLING MADE FROM COMPOSITE MATERIALS USING NUMERICAL SIMULATIONS SVOČ FST 2018

Estimation of the Residual Stiffness of Fire-Damaged Concrete Members

SPECTRUM FATIGUE LIFETIME AND RESIDUAL STRENGTH FOR FIBERGLASS LAMINATES IN TENSION

Sabah Shawkat Cabinet of Structural Engineering Walls carrying vertical loads should be designed as columns. Basically walls are designed in

Stress-strain response and fracture behaviour of plain weave ceramic matrix composites under uni-axial tension, compression or shear

Implementation of fatigue model for unidirectional laminate based on finite element analysis: theory and practice

Bending Load & Calibration Module

VALIDATION of CoDA SOFTWARE for COMPOSITES SYNTHESIS AND PRELIMINARY DESIGN (or GETTING COMPOSITES USED - PART 2 )

Mechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002

EVALUATION OF DAMAGE DEVELOPMENT FOR NCF COMPOSITES WITH A CIRCULAR HOLE BASED ON MULTI-SCALE ANALYSIS

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

Failure analysis of serial pinned joints in composite materials

Interlaminar fracture characterization in composite materials by using acoustic emission

DRAPING SIMULATION. Recent achievements and future trends. Dr. Sylvain Bel LGCIE University Lyon 1

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

5. STRESS CONCENTRATIONS. and strains in shafts apply only to solid and hollow circular shafts while they are in the

Effects of Resin and Fabric Structure

ANALYSIS OF YARN BENDING BEHAVIOUR

High Tech High Top Hat Technicians. An Introduction to Solid Mechanics. Is that supposed to bend there?

An orthotropic damage model for crash simulation of composites

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

NORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts.

SERVICEABILITY LIMIT STATE DESIGN

Tensile behaviour of anti-symmetric CFRP composite

SUB-SURFACE DAMAGE LOCATION AND IDENTIFICATION USING INFRA-RED TECHNIQUES

Discrete Element Modelling of a Reinforced Concrete Structure

ASPECTS CONCERNING TO THE MECHANICAL PROPERTIES OF THE GLASS / FLAX / EPOXY COMPOSITE MATERIAL

A Constitutive Model for DYNEEMA UD composites

In-situ local strain measurement in textile composites with embedded optical fibre sensors

Analysis of nonlinear shear deformations in CFRP and GFRP textile laminates

Advanced Strength of Materials Prof S. K. Maiti Mechanical Engineering Indian Institute of Technology, Bombay. Lecture 27

Chapter 8 Structural Design and Analysis. Strength and stiffness 5 types of load: Tension Compression Shear Bending Torsion

Static and Time Dependent Failure of Fibre Reinforced Elastomeric Components. Salim Mirza Element Materials Technology Hitchin, UK

Evaluation Axisymmetric Analysis of Thermal Stress Residual Near Fiber/Epoxy Interface

Hardened Concrete. Lecture No. 16

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

The objective of this experiment is to investigate the behavior of steel specimen under a tensile test and to determine it's properties.

THE ROLE OF DELAMINATION IN NOTCHED AND UNNOTCHED TENSILE STRENGTH

SOME RESEARCH ON FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS

On characterising fracture resistance in mode-i delamination

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

Chapter 7. Highlights:

A Study on the Tube of Integral Propeller Shaft for the Rear-wheel Drive Automobile Using Carbon Composite Fiber

Feasibility of dynamic test methods in classification of damaged bridges

FINITE ELEMENT AND EXPERIMENTAL STUDY OF NOVEL CONCEPT OF 3D FIBRE CELL STRUCTURE

Fatigue Analysis of Wind Turbine Composites using Multi-Continuum Theory and the Kinetic Theory of Fracture

MATERIALS AND METHODS

Analysis on propulsion shafting coupled torsional-longitudinal vibration under different applied loads

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

Fatigue lifetime of glass fabric/epoxy composites

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

Mechanical Behavior of Circular Composite Springs with Extended Flat Contact Surfaces

SKIN-STRINGER DEBONDING AND DELAMINATION ANALYSIS IN COMPOSITE STIFFENED SHELLS

Transcription:

Experimental setup for and numerical modelling of bending fatigue experiments on plain woven glass/epoxy composites W. Van Paepegem a* and J. Degrieck b a Research Assistant of the Fund for Scientific Research Flanders (F.W.O.), Dpt. of Mechanical Construction and Production, Sint-Pietersnieuwstraat 41, 9000 Gent, Belgium b Professor, Dpt. of Mechanical Construction and Production, Sint-Pietersnieuwstraat 41, 9000 Gent, Belgium Abstract In general fatigue of fibre-reinforced composite materials is a quite complex phenomenon, and a large research effort is being spent on it today. Due to deficiencies in the current life prediction methodologies for these materials, composite structures are often overdesigned: large factors of safety are adopted and extensive prototype-testing is required to allow for an acceptable life time prediction. This paper presents an investigation of the fatigue performance of plain woven glass/epoxy composite materials and of the numerical modelling of these composites behaviour under fatigue. First the experimental setup which has been developed for bending fatigue experiments, is discussed. The materials used are plain woven glass/epoxy specimens in two configurations: [#0º] 8 and [#45º] 8. Experiments show that these two specimen types, although being made of the same material, have a quite different damage behaviour and that the stiffness degradation follows a different path. Next a numerical model is presented which allows one to describe the degradation behaviour of the composite specimen during its fatigue life. This model has been implemented in a mathematical software package (Mathcad TM ) and proves to be a useful tool to study the fatigue degradation behaviour of composite materials. Keywords: fatigue, composites, bending experiments, damage model, residual stiffness 1 Introduction Fibre-reinforced composite materials are finding increasing application in aerospace structures, naval industry and other high-tech designs, because of their high specific stiffness and strength. Although the fatigue behaviour of fibre-reinforced composites has been studied for a long time, it is still not possible to make adequate predictions about the fatigue life and degradation of stiffness and strength without extensive ad hoc investigation. Research on fatigue of composite materials is conducted by performing numerous fatigue experiments, whereby tension-tension fatigue and tension-compression fatigue are the most favourable ways of working, since damage is developing more or less equally in all layers of the composite specimen (Fuji et al. [1], Schulte et al. [] and Hansen []). Bending fatigue experiments on the other hand have been reported by only a few authors (Ferry et al. [3], Herrington and Doucet [5] and Chen and Matthews [6]). Although fatigue experiments in tension and compression remain indispensable for accurately modelling the fatigue behaviour of composites, bending fatigue experiments yield a lot of useful information as well: - stresses and strains can vary along the gauge length of the specimen. Damage is thus varying gradually along the specimen length, providing much more information from one experiment, - the different behaviour and damage development of the specimen at the tension- and compression side can be monitored for any tension-compression ratio, * Author to whom correspondence should be addressed (E-mail: Wim.VanPaepegem@rug.ac.be).

- and permanent deformations of the composite specimen can be easily measured. Although stresses nor strains are constant along the length of the specimen (only a constant displacement is imposed at the lower end of the specimen (see infra)), the geometry and loading conditions are very simple, such that numerical validation and geometrical modelling are relatively straightforward. Material and testing procedures The material used in this study is a plain woven glass/epoxy composite. The plain woven glass fabric was stacked in eight layers, with two different stacking sequences. With the first type, the warp direction of all eight layers was aligned with the loading direction (denoted as [#0º] 8, where 0 means that the warp direction of each of the eight layers has been aligned with the loading direction and where the symbol # refers to the fabric reinforcement type). With the second stacking sequence, the angle between the warp direction of all layers and the loading direction was 45º (denoted as [#45º] 8 ). All composite specimens were manufactured using the resin-transfer-moulding technique. After curing they had a thickness of.7 mm. The samples were cut to dimensions of 145 mm long by 30 mm wide on a water-cooled diamond saw. The in-plane elastic properties of the composite laminates were determined using dynamic modal analyis (Sol et al. [7]). They are listed in Table 1. Table 1 Measured in-plane elastic properties of the composite laminates E 11 (GPa) 4.57 E (GPa) 3.94 ν 1 (-) 0.153 G 1 (GPa) 4.83 These values compare very well with a micromechanical analysis, where the properties of the E-glass fibre and the epoxy matrix are filled in separately. Each woven layer is modelled as two unidirectional layers, stacked in a [0º/90º] sequence. With the fibre volume fraction ν f being 0.48, the estimated homogenized properties are listed in Table. Table Estimated in-plane elastic properties of the composite laminates E 11 (GPa) 3.51 E (GPa) 3.51 ν 1 (-) 0.143 G 1 (GPa) 3.95 3 Experimental setup for bending fatigue The experimental setup was developed specially for flexural fatigue tests on cantilever beam specimens and is shown in Figure 1. Although more sophisticated commercial servohydraulic testing machines are capable of performing the same fatigue experiments, the maintenance and service costs are very high, compared to the costs of use of this testing machine.

Figure 1 Experimental setup for bending fatigue experiments. The outcoming shaft of the motor has a rotational speed of 185 rpm. The power is transmitted by a V-belt to a second shaft, which provides a fatigue testing frequency of.3 Hz. The influence of frequency can be assumed to be small in this range of values (Manger et al. [8]). The power transmission through a V-belt ensures that the earthing of the motor and the earthing of the measurement system are electrically isolated. Figure Schematic drawing of the crank-linkage mechanism. The second shaft bears a crank-linkage mechanism, which is shown in Figure. This mechanism imposes an alternating displacement on the hinge (point B) that connects the linkage with the lower clamp of the composite specimen. At the upper end the specimen is clamped. Hence the samples are loaded as a composite cantilever beam. The amplitude of the imposed displacement is a controllable parameter and the adjustable crank allows to choose

between single-faced and double-faced bending, i.e. the deflection can vary from zero to the maximum deflection in one direction, respectively in two opposite directions. A full Wheatstone bridge on the linkage is used to measure the force acting on the composite specimen. The selection of the strain gauges had to be done carefully, since the linkage is fatigued as well and internal energy dissipation leads to small temperature rises at the site of the active strain gauges, which are not compensated by the dummy strain gauges in the neighbourhood. Therefore two self-temperature compensated strain gauges are placed on the two opposite surfaces of the linkage and they measure both the longitudinal and transverse strains. In that way temperature effects, as well as small bending moments, are compensated. 4 Experimental results As mentioned above, two stacking sequences were used: [#0º] 8 and [#45º] 8 specimens. The stacking sequence [#0º] 8 is schematically illustrated in Figure 3. Figure 3 Definition of the [#0º] 8 stacking sequence for plain woven glass/epoxy specimens. Fatigue experiments were performed with different values of the imposed displacement, as well as with single- and double-faced bending. To characterize each experiment, the u min displacement ratio R d = (analogous to the stress ratio R) is defined, whereby the u max minimum deflection is not necessarily zero. When the displacement u max and the length L between the both clamps (Figure ) are given too, all parameters of the fatigue experiment are determined. Now some remarkable observations from the bending fatigue experiments are discussed. When the maximum displacement u max is the same for both specimens, the stiffness degradation path is quite different for [#0º] 8 and [#45º] 8 specimens. Figure 4 shows typical force-cycle histories for both stacking sequences. The abscis contains the number of cycles;

the ordinate axis shows the force (Newton), which is measured by the strain gauge bridge during the fatigue tests at constant bending displacement. Force-cycle histories of [#0 ] 8 and [#45 ] 8 specimens 10 110 100 90 80 u max = 3.3 mm R d = 0.0 L = 54.0 mm [#0 ] 8 glass/epoxy [#45 ] 8 glass/epoxy Force [N] 70 60 50 40 30 0 10 0 0 100000 00000 300000 400000 No. of cycles [-] Figure 4 Bending fatigue performance of [#0 ] 8 and [#45 ] 8 specimens at constant bending displacement. The [#0 ] 8 specimens degrade gradually and their stiffness is reduced significantly after 400,000 cycles. The initial force on the [#45 ] 8 specimens is smaller, because their stiffness is lower. However, after nearly 300,000 cycles, their remaining stiffness has become larger than that of the [#0 ] 8 specimens. Not only the force-cycle histories are different for both specimen types; Figure 5 shows that their fatigue damage pattern differs as well. Figure 5 Degradation of a [#0 ] 8 (left) and a [#45 ] 8 (right) specimen. In the [#0 ] 8 specimens, several fibres are broken at the clamped end of the specimen and a sort of hinge is formed, while in the [#45 ] 8 specimens a more regular pattern is found. Further away from the fixation a regular pattern of transverse cracking in the weft tows at the fabric cross-over points is shown in the [#0 ] 8 specimens (see also zoomed-in drawing in Figure 3). This characteristic damage state is very similar to the steady state matrix cracks in tension fatigue (Fujii et al. [1] and Manger et al. [8]).

With both types of specimens, the outermost layers which have been subjected to the largest tensile stresses, are severely damaged, while the outermost layers at the compression side show no observable sign of damage. This brings about an important result: due to the growth of damage and the degradation of the bending stiffness, there is a continuous redistribution of the stresses in the cross-section, especially near the fixation where damage growth is dominant. The position of the neutral fibre (according to its definition in the classical beam theory) does not remain in the middle of the cross-section, but tends to move towards the compression side and transfers the load to that zone. Further a considerable permanent deformation did remain after removing the grips from the [#45 ] 8 specimens, while this is not the case for the [#0 ] 8 specimens. Such deformations are difficult to observe when specimens are loaded in fatigue tension and/or compression because of the very small displacements, but are easy to measure in bending, because the displacements are much larger. The bending fatigue experiments thus show that a more generalized fatigue damage model should take into account the accumulation of permanent strain for certain stacking sequences. Another interesting observation is that there exists a fatigue threshold stress in bending. Figure 6 illustrates a typical experiment on a [#0 ] 8 specimen. The imposed displacement ranges from zero to a maximum value and the specimens always bends in the same direction (single-faced bending). First the composite specimen has been fatigued during 770,000 cycles with a small constant-amplitude displacement (u max = 7.1 mm). The force, necessary to bend the composite specimen and measured by the Wheatstone bridge, is constant during these 770,000 cycles and equals approximately 0 N. Next a larger displacement (u max = 9.5 mm) is imposed, which is held constant during 530,000 cycles more. Then indeed degradation is initiating and growing fast, until complete damage has occurred. Degradation of [#0 ] 8 specimen with small and large displacement 110 100 R d = 0.0 L = 54.0 mm 90 80 70 u max = 7.1 mm u max = 9.5 mm Force [N] 60 50 40 30 0 10 0 00000 400000 600000 800000 1000000 100000 No. of cycles [-] Figure 6 Degradation of a [#0 ] 8 specimen at a small and a large constant displacement. Figure 4 and Figure 5 illustrate that the degradation of the [#45 ] 8 specimens is mainly due to matrix cracks and shows a very gradual path with a sharp initial stiffness reduction. Therefore in the next paragraph a damage accumulation model is proposed to describe the experimental observations through numerical simulations for the [#45 ] 8 specimens. The simulation for the

[#0 ] 8 specimens requires a more complex model, because several damage mechanisms are acting together (matrix cracks, fibre-matrix interface failure, fibre fracture). 5 Local damage model In the opinion of the authors, fatigue models can be classified in three major categories: - fatigue life models, which do not take into account the actual degradation mechanisms but use S-N curves (influence of stress amplitude on occurrence of final failure) or Goodmantype diagrams (influence of mean stress level) and introduce some sort of fatigue failure criterion, - residual strength models, which describe the degradation of the initial static strength, equally without taking into account stiffness reduction, - and finally damage accumulation models which use one or more damage variables related to measurable manifestations of damage (number of transverse matrix cracks, delamination size) or to the residual stiffness. Although the fatigue behaviour of fibre-reinforced composites is fundamentally different from the behaviour exposed by metals, many models have been established which are based on the well-known S-N curves. This approach requires extensive experimental work and does not take into account the actual damage mechanisms, such as matrix cracks and fibre fracture. With regard to residual strength models, measurement of strength during damage development in a material is not feasible, because it is a destructive method and only one measurement can be made for one specimen. Moreover it is very difficult to compare damage states between two specimens, because the scatter on the strength is large compared to the scatter on the stiffness. On the other hand stiffness can be measured frequently or even continuously during damage development and does not affect the damage development as such. Therefore stiffness is a potential non-destructive parameter which can be used to monitor the in-service fatigue damage in the composite specimen. Further, in a fibre-reinforced composite the extent of the damage zones grows, while the damage-type in these zones can change (e.g. small matrix cracks leading to large size delaminations). The gradual deterioration of a fibre-reinforced composite with a loss of stiffness in the damaged zones leads to a redistribution of stress and a reduction of stress concentrations in a structural component. As a consequence an appraisal of the actual state or a prediction of the final state (when and where final failure is to be expected) requires the simulation to follow the complete path of successive damage states. This is confirmed by the experiments which show that the normal stresses acting on the crosssectional area of the clamped end, are redistributed during fatigue life and that the neutral fibre (according to the definition used in beam theory) is moving across the height of the specimen, because damage in the outermost layers at the tensile side is growing much faster than in the neighbouring zones. Therefore numerical calculations at decisive stages in fatigue life, which take into account the stiffness degradation, are indispensable to make an estimate of the redistribution of the stresses over the structure. Here a unidirectional model for fatigue damage is proposed. One local damage variable D is associated with the longitudinal stiffness loss. The stresses and strains are then related by the commonly used equation: σ = E 0.(1 D). ε (1) Strains are calculated with the equations of classical beam theory for a cantilever beam that is loaded with a force and a moment in the point A (see Figure ).

The form for the fatigue evolution law is given by: where: - D : local damage variable - N : number of cycles - σ : amplitude of the applied cyclic loading - A, b and c : three material constants c A.( σ) dd in tension = b (1 D) () dn 0 in compression This damage law is similar to the one proposed by Sidoroff and Subagio [9], but uses the stress amplitude instead of the strain amplitude. When the axial force is supposed to remain zero and when only a bending moment exists, the position of the neutral fibre y 0 at each moment of time is calculated as: h + [ 1 D(y) ] h h + y = (3) 0 h.y dy 1 D(y) dy where: - y : thickness-coordinate, with y = 0 in the middle of the specimen s thickness - h : total thickness of the specimen The degraded bending stiffness EI becomes (with b the specimen s width): h + [ 1 D(y) ] EI = b.e 0..y dy (4) h The remaining equation is formed by using the principle of virtual work and where F is solved from this transcendental equation: u ( F.a. α + F α ) max F.u F + F.a.u M + a.sin M. = (5) F where: - F : force measured by the strain gauges and acting on the hinge (point B in Figure ) - a : length of the lower clamp (Figure ) - u F and α F : displacement, respectively rotation angle for a unit force F = 1 N in the point A (Figure ) - u M and α M : displacement, respectively rotation angle for a unit moment M = 1 Nm in the point A (Figure )

This numerical model can be easily implemented in a mathematical software package such as Mathcad. The numerical integration formulas must be chosen such that the second degree polynomials are exactly integrated. This is the case for the Simpson s rule, which is a Newton-Cotes quadrature formula. Because the increase of the damage variable D during one cycle is so small, the integrations must be exact indeed, otherwise the relative error on the calculation of the bending stiffness EI may be larger than the increase of the damage variable itself. Therefore the conventional first-order trapezium method is not suited for this purpose. Figure 7 shows the experimental and simulated results for the [#45 ] 8 specimen for which the force-cycle history was shown in Figure 4. The imposed displacement u max is 3.3 mm, the displacement ratio R d is zero and the length L (Figure ) is 54.0 mm. The measured force initially equals 74.0 N, the total number of cycles in the experiment is 400,000 and the measured remaining force then equals 53.3 N. With regard to the simulation, the composite specimen is modelled with a mesh of 400 integration points and the calculation uses 40 integration steps. The values of the constants A, b and c (see equation ()) are respectively: 10-18, 0.45 and 6.5. The remaining force at 400,000 cycles is 51.7 N, which results in an error of only 3% after 400,000 cycles. There is a very good agreement between the experimental data and the theoretical simulation, which proves that the proposed damage law is capable of simulating the kind of fatigue degradation behaviour that has been observed with [#45 ] 8 specimens. 80 70 Force versus number of cycles for [#45] 8 specimen u max = 3.3 mm R d = 0.0 L = 54.0 mm 60 Force [N] 50 40 30 experiment simulation 0 10 0 0 100000 00000 300000 400000 No. of cycles [-] Figure 7 Degradation of the force versus the number of cycles for a [#45 ] 8 specimen. Figure 8 shows the results from another experiment for a [#45 ] 8 specimen. Again the imposed displacement u max is 3.3 mm, the displacement ratio R d is zero and the length L is 54.0 mm. This specimen has a slightly lower modulus of elasticity, which results in a lower initial force of 65.3 N. After 399,000 cycles, the remaining force is 50.1 N. For the simulation the values of the constants A, b and c remain the same (respectively: 10-18, 0.45 and 6.5). The force after 399,000 cycles equals 48. N for the simulation, which results in an error of 3.8%. The simulated results from this experiment can be used now for a more detailed discussion of the strain-, stress- and damage distributions in the composite specimen during fatigue life.

70 60 Force versus number of cycles for [#45] 8 specimen u max = 3.3 mm R d = 0.0 L = 54.0 mm 50 Force [N] 40 30 experiment simulation 0 10 0 0 100000 00000 300000 400000 No. of cycles [-] Figure 8 Degradation of the force versus the number of cycles for a [#45 ] 8 specimen. Figure 9 shows the strain distribution over the cross-section near the clamped end at a few decisive stages in fatigue life. The ordinate axis represents the height y (mm) in the crosssection. As mentioned above the specimens have a thickness of.7 mm, so the coordinates of the integration points are in the range (-1.36 mm, +1.36 mm). At the first cycle, the tensile and compressive strain are equal and the strain is zero in the middle of the cross-section. When damage is growing, the neutral fibre is moving towards the compression side of the specimen. Strain distribution in the clamped cross-section cycle 1 cycle 18,315 cycle 100,850 cycle 399,000 1.4 1.0 0.6 height y [mm] 0. -0.05-0.015-0.005 0.005 0.015 0.05-0. -0.6-1.0-1.4 strain [-] Figure 9 Strain distribution in the clamped cross-section. The influence of damage can be observed even more clearly when the distribution of the stresses is calculated, as shown in Figure 10. At the first cycle, the stress distribution is linear and the normal stress is zero in the middle of the cross-section. When damage is initiating, the tensile stresses in the outermost layers are relaxed and load is transferred towards the inner layers. Because the damage law assumes that there is no damage growth at the compressive

side, the peak tensile stresses are moving towards the compression side and as a consequence the neutral fibre is moving down. Stress distribution in the clamped cross-section cycle 1 cycle 18,315 cycle 100,850 cycle 399,000 1.4 1.0 0.6 height y [mm] 0. -00-100 0 100 00-0. -0.6-1.0-1.4 stress [MPa] Figure 10 Stress distribution in the clamped cross-section. Finally Figure 11 shows the damage growth in the integration points in the cross-section. The value of the damage variable D is lying between zero (no damage) and one (final failure of that integration point). This damage distribution indeed reflects the experimental observations, where the outermost layers at the tension side are severely damaged ( D 1 ), while the compression side shows no observable damage. 1.4 Damage distribution in the clamped cross-section 1.0 0.6 height y [mm] 0. 0.0 0.1 0. 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0-0. -0.6-1.0 cycle 1 cycle 1,07 cycle 5,96 cycle 18,315 cycle 100,850 cycle 399,000-1.4 damage [-] Figure 11 Damage distribution in the clamped cross-section.

Further investigation is required to incorporate some other experimental observations into the damage accumulation model, such as the accumulation of permanent strain in case of the [#45 ] 8 specimens and the existence of a fatigue threshold stress. This will probably lead to a fatigue model, where damage initiation and propagation will be described by different terms in the damage evolution law. Next the fatigue damage behaviour of the [#0 ] 8 specimens will require a more complex damage accumulation model, because more than one damage mechanism is present: matrix cracking, fibre-matrix interface failures and fibre fracture. However the discussed results prove that a damage accumulation model, related to the degradation of the stiffness properties, is capable of predicting a continuous redistribution of stresses across the structure, as can be observed in the experiments. 6 Conclusions An experimental setup for bending fatigue experiments has been built. Results show that bending fatigue tests yield important additional information that cannot be recovered from the conventional tension fatigue tests. A damage accumulation model has been proposed that allows to check the experimentally observed damage initiation and growth. The fatigue model can be easily implemented in a mathematical software package and shows good agreement for simulations of the degradation behaviour of [#45 ] 8 specimens. The experiments indicate that stresses are continuously redistributed across the structure during fatigue life. As the damage accumulation model accounts for the degradation of the local stiffness properties, the stress redistribution can be theoretically simulated and the damage initiation and growth can be predicted. Acknowledgements The author W. Van Paepegem gratefully acknowledges his finance through a grant of the Fund for Scientific Research Flanders (F.W.O.). The authors also express their gratitude to Syncoglas for their support and technical collaboration. References 1. Fujii, T., Amijima, S. and Okubo, K. (1993). Microscopic fatigue processes in a plain-weave glass-fibre composite. Composites Science and Technology, 49, 37-333.. Schulte, K., Reese, E. and Chou, T.-W. (1987). Fatigue behaviour and damage development in woven fabric and hybrid fabric composites. In : Matthews, F.L., Buskell, N.C.R., Hodgkinson, J.M. and Morton, J. (eds.). Sixth International Conference on Composite Materials (ICCM-VI) & Second European Conference on Composite Materials (ECCM-II) : Volume 4. Proceedings, 0-4 July 1987, London, UK, Elsevier, pp. 4.89-4.99. 3. Hansen, U. (1997). Damage development in woven fabric composites during tension-tension fatigue. In : Andersen, S.I., Brøndsted, P., Lilholt, H., Lystrup, Aa., Rheinländer, J.T., Sørensen, B.F. and Toftegaard, H. (eds.). Polymeric Composites - Expanding the Limits. Proceedings of the 18th Risø International Symposium on Materials Science, 1-5 September 1997, Roskilde, Denmark, Risø International Laboratory, pp. 345-351.

4. Ferry, L., Gabory, D., Sicot, N., Berard, J.Y., Perreux, D. and Varchon, D. (1997). Experimental study of glass-epoxy composite bars loaded in combined bending and torsion loads. Fatigue and characterisation of the damage growth. In : Degallaix, S., Bathias, C. and Fougères, R. (eds.). International Conference on fatigue of composites. Proceedings, 3-5 June 1997, Paris, France, La Société Française de Métallurgie et de Matériaux, pp. 66-73. 5. Herrington, P.D. and Doucet, A.B. (199). Progression of bending fatigue damage around a discontinuity in glass/epoxy composites. Journal of Composite Materials, 6(14), 045-059. 6. Chen, A.S. and Matthews, F.L. (1993). Biaxial flexural fatigue of composite plates. In : Miravete, A. (ed.). ICCM/9 Composites : properties and applications. Volume VI. Proceedings of the Ninth International Conference on Composite Materials, 1-16 July 1993, Madrid, Spain, Woodhead Publishing Limited, pp. 899-906. 7. Sol, H., Hua, H., De Visscher J., Vantomme, J. and De Wilde, W.P. (1997). A mixed numerical/experimental technique for the nondestructive identification of the stiffness properties of fibre reinforced composite materials. NDT&E International, Vol. 30(), pp. 85 91. 8. Manger, C.I.C., Ogin, S.L., Smith, P.A. and Greaves, R.P. (1997). Damage development in plain weave GFRP. In : Scott, M.L. (ed.). ICCM/11 Volume V : Textile composites and characterisation. Proceedings of the Eleventh International Conference on Composite Materials (ICCM/11), 14-18 July 1997, Gold Coast, Queensland, Australia, Woodhead Publishing Limited, pp. V.58-V.66. 9. Sidoroff, F. and Subagio, B. (1987). Fatigue damage modelling of composite materials from bending tests. In : Matthews, F.L., Buskell, N.C.R., Hodgkinson, J.M. and Morton, J. (eds.). Sixth International Conference on Composite Materials (ICCM-VI) & Second European Conference on Composite Materials (ECCM-II) : Volume 4. Proceedings, 0-4 July 1987, London, UK, Elsevier, pp. 4.3-4.39.