DAMAGE MECHANICS MODEL FOR OFF-AXIS FATIGUE BEHAVIOR OF UNIDIRECTIONAL CARBON FIBER-REINFORCED COMPOSITES AT ROOM AND HIGH TEMPERATURES

Similar documents
RATE-DEPENDENT OFF-AXIS COMPRESSIVE STRENGTH OF A UNIDIRECTIONAL CARBON/EPOXY LAMINATE AT HIGH TEMPERATURE

EFFECT OF ALTERNATE CHANGE IN STRESS RATIO ON FATIGUE STRENGTH OF WOVEN FABRIC CFRP LAMINATE AND LIFE PREDICTION USING THE ANISOMORPHIC CFL DIAGRAM

Nonlinear Constant Fatigue Life Diagrams for CFRP Laminates

IJSER 1. INTRODUCTION. M.Elhadary

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

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

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

American Society for Testing and Materials (ASTM) Standards. Mechanical Testing of Composites and their Constituents

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

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

Open-hole compressive strength prediction of CFRP composite laminates

THE MUTUAL EFFECTS OF SHEAR AND TRANSVERSE DAMAGE IN POLYMERIC COMPOSITES

Materials and Structures. Indian Institute of Technology Kanpur

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

THE ROLE OF DELAMINATION IN NOTCHED AND UNNOTCHED TENSILE STRENGTH

MODELING OF THE BEHAVIOR OF WOVEN LAMINATED COMPOSITES UNTIL RUPTURE

BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS

EFFECTS OF MODELING ASSUMPTIONS ON THE ACCURACY OF SPECTRUM FATIGUE LIFETIME PREDICTIONS FOR A FIBERGLASS LAMINATE

QUESTION BANK Composite Materials

ME 582 Advanced Materials Science. Chapter 2 Macromechanical Analysis of a Lamina (Part 2)

ME 243. Mechanics of Solids

Composite Structures. Indian Institute of Technology Kanpur

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

DESIGN FOR FATIGUE STRENGTH

Hygrothermal stresses in laminates

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

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

Stress-Strain Behavior

SPECTRUM FATIGUE LIFETIME AND RESIDUAL STRENGTH FOR FIBERGLASS LAMINATES IN TENSION

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

Effects of Resin and Fabric Structure

Chapter 7. Highlights:

MECE 3321 MECHANICS OF SOLIDS CHAPTER 3

MECHANICS OF MATERIALS

Composite Materials Research Laboratory, State University of New York at Buffalo, Buffalo, NY , USA

Module 7: Micromechanics Lecture 25: Strength of Materials Approach. Introduction. The Lecture Contains. Effective Transverse Modulus

MICROMECHANICAL ANALYSIS OF FRP COMPOSITES SUBJECTED TO LONGITUDINAL LOADING

Module 4: Behaviour of a Laminae-II. Learning Unit 1: M1. M4.1 Mechanics of Composites. M4.1.1 Introduction to Mechanics of Composites

ME 2570 MECHANICS OF MATERIALS

Outline. Tensile-Test Specimen and Machine. Stress-Strain Curve. Review of Mechanical Properties. Mechanical Behaviour

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

BioMechanics and BioMaterials Lab (BME 541) Experiment #5 Mechanical Prosperities of Biomaterials Tensile Test

Fatigue lifetime of glass fabric/epoxy composites

Most of the material in this package is based on a recently published book. This is:

EFFECTS OF MOLECULAR STRUCTURE ON MACROSCOPIC MECHANICAL PROPERTIES OF AN ADVANCED POLYMER (LARC -SI)

Accelerated Testing Methodology for Long Term Durability of CFRP

Failure analysis of serial pinned joints in composite materials

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

Part 7. Nonlinearity

Materials and Structures. Indian Institute of Technology Kanpur

NUMERICAL MODELLING OF COMPOSITE PIN- JOINTS AND EXPERIMENTAL VALIDATION

Tensile Stress Acoustic Constants of Unidirectional Graphite/Epoxy Composites

Chapter 2 - Macromechanical Analysis of a Lamina. Exercise Set. 2.1 The number of independent elastic constants in three dimensions are: 2.

TENSILE FATIGUE BEHAVIOR OF SINGLE FIBRES AND FIBRE BUNDLES

Crashworthiness of composite structures: Experiment and Simulation

USING A HOMOGENIZATION PROCEDURE FOR PREDICTION OF MATERIAL PROPERTIES AND THE IMPACT RESPONSE OF UNIDIRECTIONAL COMPOSITE

POLYURETHANE SURFACE TREATMENT ON TWO KINDS OF BASALT FIBER COMPOSITE AND MECHANICAL PROPERTIES COMPARISON

Coupling of plasticity and damage in glass fibre reinforced polymer composites

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

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

9 Strength Theories of Lamina

RELATIONSHIP BETWEEN RADIAL COMPRESSIVE MODULUS OF ELASTICITY AND SHEAR MODULUS OF WOOD Jen Y. Liu Research Engineer

Fatigue life. prediction of composites. and composite. structures. Vassilopoulos WOQDHEAD PUBLISHING LIMITED. Anastasios P. Cambridge New Delhi

TESTING AND ANALYSIS OF COMPOSITE SKIN/STRINGER DEBONDING UNDER MULTI-AXIAL LOADING.

Mechanical properties 1 Elastic behaviour of materials

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

School of Materials Science and Engineering, UNSW Sydney, Australia 2. School of Mechanical and Manufacturing Engineering, UNSW Sydney, Australia

Mechanical Properties of Materials

An integrated approach to the design of high performance carbon fibre reinforced risers - from micro to macro - scale

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

STRESS ANALYSIS OF BONDED JOINTS IN PULTRUDED GRP COMPONENTS

Assessment Methods of Mechanical Properties of Composite Materials

EXPERIMENTAL CHARACTERISATION AND MODELLING OF NON- LINEAR STRESS-STRAIN BEHAVIOUR AND STRAIN RATE EFFECTS FOR UNIDIRECTIONAL CARBON-EPOXY

NUMERICAL FEM ANALYSIS FOR THE PART OF COMPOSITE HELICOPTER ROTOR BLADE

Mechanical and Thermal Properties of Coir Fiber Reinforced Epoxy Composites Using a Micromechanical Approach

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

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

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

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

Fracture Mechanics, Damage and Fatigue: Composites

Prediction of Elastic Constants on 3D Four-directional Braided

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

PREDICTION OF OPEN HOLE COMPRESSIVE FAILURE FOR QUASI-ISOTROPIC CFRP LAMINATES BY MMF/ATM METHOD

Computational Analysis for Composites

Influence of fibre proportion and position on the machinability of GFRP composites- An FEA model

Prediction of The Ultimate Strength of Composite Laminates Under In-Plane Loading Using A Probabilistic Approach

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

Evaluation of in-plane orthotropic elastic constants of paper and paperboard

Geometric and Material Property Effects on the Strength of Rubber-Toughened Adhesive Joints

APPLICATION OF A SCALAR STRAIN-BASED DAMAGE ONSET THEORY TO THE FAILURE OF A COMPLEX COMPOSITE SPECIMEN

DYNAMIC FAILURE ANALYSIS OF LAMINATED COMPOSITE PLATES

KINK BAND FORMATION OF FIBER REINFORCED POLYMER (FRP)

LAMB WAVE RESPONSE OF FATIGUED COMPOSITE SAMPLES. Reveiw of Progress in Quantitative Nondestructive Evaluation, Vol 13B (1994) pp.

SHEAR CAPACITY OF REINFORCED CONCRETE COLUMNS RETROFITTED WITH VERY FLEXIBLE FIBER REINFORCED POLYMER WITH VERY LOW YOUNG S MODULUS

Tensile behaviour of anti-symmetric CFRP composite

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

Appendix A: Laboratory Report Format

THE USE OF NEURAL NETWORKS IN THE PREDICTION OF THE FATIGUE LIFE OF DIFFERENT COMPOSITE MATERIALS

12/8/2009. Prof. A.K.M.B. Rashid Department of MME BUET, Dhaka

Laboratory 4 Bending Test of Materials

Transcription:

DAMAGE MECHANICS MODEL FOR OFF-AXIS FATIGUE BEHAVIOR OF UNIDIRECTIONAL CARBON FIBER-REINFORCED COMPOSITES AT ROOM AND HIGH TEMPERATURES M. Kawai Institute of Engineering Mechanics University of Tsukuba, Tsukuba 305-8573, Japan SUMMARY: Tension-tension fatigue tests were performed for unidirectional T800H/Epoxy at room temperature (RT) and 100 C, respectively, using seven kinds of plain coupon specimens with different fiber orientation angles. The characteristic nature of the fiber-matrix interface controlled off-axis fatigue failure was similar at RT and 100 C. Irrespective of offaxis angle, almost linear S-N relationships were obtained in the whole range of fatigue testing up to 10 6 cycles. Normalizing the maximum fatigue stress with respect to the static tensile strength at test temperature, the off-axis fatigue data eventually fell on a single S-N relationship. A non-dimensional effective stress defined on the basis of the classical static failure theory succeeded in coping with the directional nature of the off-axis fatigue behavior in a unified manner. By using the non-dimensional effective stress in the continuum damage mechanics approach, a general failure model was developed to describe the observed off-axis S-N relationships. KEYWORDS: Off-Axis Fatigue, Temperature Dependence, Fatigue Damage, S-N Curve, Non-Dimensional Effective Stress, Damage Mechanics, Carbon Fiber, Epoxy Resin INTRODUCTION For unidirectional composites consisting of high stiffness/strength fiber and brittle polymer matrix, the longitudinal fatigue failure occurs usually with some parallel splittings extending along fibers into the end tab regions and pullouts of strips fractured transversely in either of the end tab regions. In the case of the off-axis fatigue loading, on the other hand, the ultimate failure takes place in a single weakest cross section along the reinforcing fibers, regardless of off-axis angle. Therefore, it is obvious that the off-axis fatigue failure of unidirectional composites is influenced by the mechanical properties of matrix and fiber-matrix interface/interphase [1]. The ductility and strength of polymer matrix and fiber-matrix interface/interphase markedly change with temperature and moisture content. The variations in these mechanical properties result in significant effects on the inelastic deformation and damage of unidirectional laminae under off-axis fatigue loading conditions. For example, an increase in the resin ductility suppresses the extension of local damage developing in the matrix and interphase, and it enhances the nonlinear shear deformation in the fiber-matrix

interface region. A decrease in the strength of matrix or interphase lowers the off-axis strength. These imply that the changes in ductility and strength of the matrix and in adhesive strength of the fiber-matrix interfaces due to different environmental (thermal and hygroscopic) conditions have marked effects on the off-axis fatigue behavior of unidirectional composites. Therefore, it is important to elucidate effects of the changes in ductility and strength of the matrix and fiber-matrix interface/interphase on the off-axis fatigue behavior of unidirectional laminae and to include them in the fatigue design analysis for composites. The present study aims to elucidate the effects of test temperature on the off-axis fatigue behavior of unidirectional continuous carbon fiber composite T800/Epoxy and to develop a fatigue failure model based on the continuum damage mechanics. For these purposes, tension-tension fatigue tests were performed at room temperature (RT) and 100 C, respectively, using plain coupon specimens with different off-axis angles: 0, 10, 15, 20, 30, 45, and 90. From a comparison between the relative fatigue strengths at room and high temperatures, the effects of property changes in the epoxy resin due to temperature on the offaxis fatigue behavior were elucidated. Capabilities of the conventional fatigue strength parameters to systematically describe the off-axis S-N relationships were also examined. To model the fiber-matrix interface controlled off-axis fatigue failure, a non-dimensional multiaxial effective stress was defined on the basis of the classical static failure theory [2]. The non-dimensional effective stress is interpreted as a theoretical strength ratio in the case of off-axis fatigue loading, and as an interfacial effective stress characterizing the fiber-matrix interface controlled off-axis fatigue failure of unidirectional laminae. Finally, by using this non-dimensional effective stress in the continuum damage mechanics approach [3], an attempt was made to develop a macroscopic fatigue failure model for predicting the observed off-axis fatigue behavior. Through numerical simulations, the validity of the proposed model was tested. MATERIAL AND SPECIMENS The material for the present study is a unidirectional (UD) composite laminate fabricated using the autoclave method from the prepreg (P2053-17, TORAY) which consists of carbon fiber (TORAYCA T800H, TORAY) and thermosetting epoxy resin (2500EP, the curing temperature of 130 C). The lay-up of the virgin laminate is [0] 12. From the optical micrograph of the cross-section perpendicular to the fiber direction, it was found that the fiber volume fraction V f was about 64 %. Nine kinds of plain coupon specimens with different off-axis angles (θ = 0, 5, 10, 15, 20, 30, 45, 60 and 90 ) were cut from 400 mm by 400 mm UD laminate panels. The specimen dimensions are specified as follows: the specimen length L = 200 mm; the gauge length L G = 100 mm ; and the thickness t = 2 mm, respectively. The specimen width (W) is 10 mm for θ = 0, and 25 mm for all other off-axis angles.

σ max 01020101010210 20 30 45 90 Nf 3104105106107 Fig. 1: Off-axis S-N relationships for unidirectional T800H/Epoxy laminate at RT Rectangular-shaped aluminum alloy tabs were bonded on both ends of the specimens using epoxy adhesive (ARALDITE ); the thickness of the end-tabs is 1.0 mm. EXPERIMENTAL PROCEDURES Off-axis static tension and tension-tension fatigue tests for T800H/2500EP were performed using a hydraulic servo testing machine at room temperature (RT, about 23 C) and high temperature (HT, 100 C), respectively. Specimens were clamped by the hydraulic grip fitted on the testing machine. The static tension tests were performed with stroke control at RT and HT, and the cross-head speed was specified to be 1.0 mm/min. The tension-tension fatigue tests were carried out at constant amplitudes under a load control mode. The fatigue load was applied at a constant stress ratio R = 0.1 in a sinusoidal waveform with a frequency of 10 Hz. The fatigue tests were terminated at the cut-off cycle of 10 6. Off-Axis Fatigue Behavior at RT EXPERIMENTAL RESULTS AND DISCUSSION S-N relationships using maximum fatigue stress The S-N relationships at RT for the off-axis angles θ = 10, 15, 30, 45 and 90 are presented in Fig. 1. The off-axis fatigue strength becomes smaller as the off-axis angle increases. The off-axis S-N relationships are well represented by the straight lines with negative slopes which start from the values of static tensile strength plotted on the vertical axis, although a small disorder is involved in the case θ = 30 at a low cycle fatigue range. The slopes of the

Fig. 2: Normalized S-N relationships at RT using fatigue elastic strain range σ/e x (θ) Fig. 3: Normalized S-N relationships at RT using experimental fatigue strength ratio σ max /σ B(exp) approximated straight lines tend to become more gradual as the off-axis angle increases. The static strength-dependence in the low cycle regime and the fatigue limit were not clearly observed in these off-axis fatigue tests. Effects of fatigue strength parameter on S-N relationships The S-N relationships plotted using a fatigue elastic strain range ε e instead of the maximum fatigue stress are presented in Fig. 2 for the fatigue data at RT discussed above. The fatigue elastic strain range was defined as ε e = σ/e x, where σ [= (1-R)σ max ] represents the fatigue stress range and E x is the Young's modulus in the fatigue loading direction. From Fig. 2, it is found that the fatigue data points separate into two groups associated with the fiber direction (θ = 0 ) and the other off-axis directions (0 < θ _ 90 ) and that the data band for the latter

σ σmax / σb(exp) σ B(exp) group is very narrow. This result implies that the fatigue failure of unidirectional composites is dominated by different physical mechanisms for these two groups. Another kind of normalization of the maximum fatigue stress can be made using the static strength in place of the Young s modulus E x. Fig. 3 shows the off-axis S-N relationships plotted using the experimental strength ratio ; denotes the observed off-axis static tensile strength. From this figure, it is observed that the normalized fatigue data points randomly distribute in a relatively narrow band regardless of off-axis angle. The same tendency has been found by Awerbuch and Hahn [4] in their experimental results for a unidirectional graphite/epoxy composite at RT. The width of the normalized fatigue data band is not negligibly small. However, it is believed as a first approximation that the experimental strength ratio is a useful parameter to cope with the orientation dependence of the off-axis fatigue behavior of unidirectional composites. Off-Axis Fatigue Behavior at 100 C S-N relationships using maximum fatigue stress The off-axis fatigue behavior observed at 100 C for the off-axis angle range 0 < θ _ 90 is shown in Fig. 4. The S-N relationships are almost linear for every off-axis angle in the range N f _ 10 6, and the fatigue strength becomes lower as the off-axis angle increases. The static strength-dependence and the fatigue limit are not clearly observed in the range of the present fatigue testing up to 10 6. These characteristics are common to the off-axis fatigue behavior at RT and 100 C. The effect of temperature can be found in the significant reduction of the absolute fatigue strength (i.e., the level of maximum fatigue stress). Compared between the S-N data at RT and 100 C for each off-axis angle in the range 0 < θ _ 90, it is found that the fatigue strength at 100 C falls to one-half to two-thirds of that at RT. This almost corresponds to the decrease in off-axis static strength due to the increase in test temperature. The fatigue strength at 100 C in the fiber direction is almost as high as that at RT, and it cannot be shown in Fig. 4. From these observations, therefore, it is confirmed that the strength reduction of the epoxy resin, fiber-matrix interface or interphase due to high temperature have significant influences on the off-axis fatigue strength of unidirectional composites.

σij σω θ()σ σgijkl 1σij σkl x Fig. 4: Off-axis S-N relationships for unidirectional T800H/Epoxy laminate at 100 C Effects of fatigue strength parameter on S-N relationships The relationships between the fatigue elastic strain range ( ε e = σ/e x ) and the logarithm of fatigue lifetime N f for the high temperature off-axis fatigue were similar to those for the RT off-axis fatigue. This result again indicates that different fatigue failure mechanisms operate for these two groups. The macroscopic fatigue failure morphology observed at 100 C was similar to that at RT. Therefore, the microscopic fatigue failure mechanisms at 100 C which characterize these two groups are also believed to be similar. Regarding the normalized S-N relationships using the experimental strength ratio σ max /σ B(exp) at high temperature, the width of the data band is slightly wider as compared with the similar plots for the RT fatigue. No appreciable effects of the off-axis angle on the fatigue lifetime can be observed in this plot. Accepting the logic used by Awerbuch and Hahn [4], therefore, it may be concluded that the experimental strength ratio becomes a useful parameter to manipulate also the off-axis angle dependence of high-temperature fatigue. NON-DIMENSIONAL EFFECTIVE STRESS [1, 2], a non- On the basis of the Tsai-Hill quadratic interaction criterion dimensional scalar quantity σ * is defined as follows: In a particular case where a tensile stress σ x is applied to an off-axis plain coupon specimen in the longitudinal (x-axis) direction, the non-dimensional effective stress becomes (1) where θ denotes the off-axis angle. The orientation factor Ω (θ) is given by (2)

σωθ( Ωθ()X2X2 + Y2 )σ ) σb(pred )1Ωθ() σmax dn KΦ Φ σ*(1r) S2 x1 Ωθ() ω max * where X and Y represent the longitudinal and transverse strengths in either tension or compression and S denotes the shear strength. The off-axis static tensile failure occurs when the following condition is satisfied: (3) The predicted tensile strength σ B(pred) in the loading direction is thus expressed as (4) The maximum fatigue stress normalized with respect to this theoretical static fracture strength is described as (5) Therefore, the maximum non-dimensional effective stress associated with the maximum fatigue stress can be interpreted as a theoretical strength ratio for the off-axis loading on plain coupon specimens. Evolution Equation for Fatigue Damage FATIGUE DAMAGE MECHANICS MODEL It is assumed that a state of the fatigue damage of unidirectional composite is described by a single scalar variable ω and its growth is prescribed by the following equation: (6) where K, n, and k are material constants. The fatigue strength parameter Φ is specified by either of the following: (7) (8)

σ Nf(k+ Nf(σ σmax 1)K Φ max *)n Fig. 5: On- and off-axis S-N relationships predicted using fatigue damage mechanics model (RT) where range σ. is the non-dimensional effective stress range corresponding to the fatigue stress Fatigue Life Equation Imposing an initial condition that ω = 0 when N = 0, Eqn 7 can be integrated to obtain the analytic expression of the fatigue damage. Assuming that a final fracture occurs when ω has reached unity, we can derive the following expression of the fatigue lifetime: (9) Setting Φ = and N f = 1 when can be represented as = 1, we find (k+1)k = 1. Then, the fatigue lifetime (10) Application to Off-Axis Fatigue In Eqn 10, the material constant to be determined is the only exponent n, except for the static strengths X, Y and S. The value of -1/n is identified as the slope of the straight line which starts from the point (logn f, log ) = (1, 1) and approximates to the normalized S-N

σmax Fig. 6: On- and off-axis S-N relationships predicted by using fatigue damage mechanics model (100 C) relationships between log and logn f. The exponents evaluated for the fatigue data at room and high temperatures are n = 22.07 and n = 22.05, respectively. The difference between these exponents is very small. Thus, the normalized master S-N relationship is independent of test temperature. Figs. 5 and 6 show the comparisons between the theoretical and experimental S-N relationships for the off-axis fatigue at RT and 100 C, respectively. It is observed that the fatigue damage mechanics model can describe the characteristic nature of the off-axis fatigue behavior at these temperatures. CONCLUSIONS Off-axis tensile fatigue behavior of unidirectional T800H/Epoxy composite was examined at room (RT) and high (100 C) temperatures under constant amplitude condition, and the influence of the temperature-dependence of the strengths of matrix and fiber-matrix interface/interphase was elucidated. A particular emphasis was placed on the identification of the mechanical parameter to describe the directional nature of the off-axis fatigue strength. It was demonstrated that a non-dimensional effective stress defined using the classical failure criterion might be used as an alternative strength parameter for a systematic description of the fatigue behavior of unidirectional composite laminae. An attempt was also made to develop a damage mechanics model for the off-axis fatigue behavior, and it was checked on the basis of the experimental results.

REFERENCES 1. Agarwal, B.D. and Broutman, L.J., Analysis and Performance of Fiber Composites, John Wiley & Sons, 1990. 2. Tsai, S.W. and Wu, E.M., «A General Theory of Strength for Anisotropic Materials», J. Composite Materials, Vol. 5, 1971, pp. 58-80. 3. Lemaitre, J., A Course on Damage Mechanics, 2nd Edition, Springer, 1996. 4. Awerbuch, J. and Hahn, H.T., «Off-Axis Fatigue of Graphite/Epoxy Composite», Fatigue of Fibrous Composite Materials, ASTM STP 723, 1981, pp. 243-273.