Microsoft Word - Khelil-ce_r3.doc ETASR - Engineering, Technology & Applied Science Research Vol. 3, No. 4, 2013, 488-496 488 www.etasr.com Khelil et al.: Modeling of Fatigue Crack Propagation in Aluminum Alloys Using an Energy… Modeling of Fatigue Crack Propagation in Aluminum Alloys Using an Energy Based Approach F. Khelil Laboratoire de Mécanique de Lille (LML), University of Lille1, France foudil.khelil@univ-lille1.fr B. Aour Laboratory of Environmental Technology Research, ENSET d’Oran, Algeria Ben_aour@yahoo.fr M. Belhouari Department of Mechanical Engineering, University of Sidi Bel Abbes, Algeria m_belhouari@yahoo.com N. Benseddiq LML, University of Lille1, France noureddine.benseddiq@ univ-lille1.fr Abstract—Materials fatigue is a particularly serious and unsafe kind of material destruction. Investigations of the fatigue crack growth rate and fatigue life constitute very important and complex problems in mechanics. The understanding of the cracking mechanisms, taking into account various factors such as the load pattern, the strain rate, the stress ratio, etc., is of a first need. In this work an energy approach of the Fatigue Crack Growth (FCG) was proposed. This approach is based on the numerical determination of the plastic zone by introducing a novel form of plastic radius. The experimental results conducted on two aluminum alloys of types 2024-T351 and 7075-T7351 were exploited to validate the developed numerical model. A good agreement has been found between the two types of results. Keywords- fatigue crack growth; energetic approach; Plastic zone; aluminum alloys. I. INTRODUCTION In these last years, the concepts of fracture mechanics allowed a better definition of the stresses and strains fields at the vicinity of crack tips under static and dynamic loadings. Cracking laws, empirical or formal, were developed in order to explore with an acceptable approximation the fatigue crack growth. Indeed, cracks’ growth is related to the existence of a plastic zone (PZ) at the crack tip, the formation and the intensification of which is accompanied by energy dissipation. Thus, the amount of the cyclic plastic strain energy may represent with precision the rate of damage at the crack tip. The use of a cyclic plastic dissipation criterion for fatigue crack growth was first proposed by Rice [1]. From this date, plastic energy approaches to fatigue crack extension prediction have been the subject of several experimental, analytical and numerical investigations [2-10]. Weertman [11] proposed that the crack advances when the accumulated plastic energy at the crack tip reaches a critical value. Then, Shozo et al [2] measured the cyclic work to produce a unit area of fatigue crack for a steel of low carbon content and for high resistance aluminum alloys, using micro strain gages stuck in the plastic zone associated with a fatigue crack. Subsequently, different techniques have been developed to evaluate the plastic energy, such as sub-grain size measurements [3], infrared thermography [4], micro-calorimetry [6] and by direct measurement of hysteresis energy under the loading line of a compact tensile (CT) specimen [12, 20]. Following the work of Bodner et al [5], Klingbeil [7] has proposed a crack growth law, in which the fatigue crack growth rate was related to the total plastic energy dissipated ahead of a crack tip under cyclic loading. This model has been further extended to mixed mode fatigue delamination of layered materials across the interface [10, 14]. Recently, Mazari et al [15] proposed an empirical correction factor which takes into account the over evaluations obtained by hysteresis loops and shows the different effects of plasticity, crack closure and opening mode. In this paper, a new approach, for the evaluation of the cyclic plastic strain energy at the crack tip in mode I, has been proposed. This approach is based on the numerical determination of the plastic zone by introducing a novel form of plastic radius. The theoretical basis related to the surface energy creation and the evolution of the energy parameters will be discussed in section 2. Section 3 is devoted to the presentation of the experimental data exploited for the validation. Then, section 4 describes the development of the numerical algorithm used for the evaluation of the cyclic plastic strain energy. The obtained results are presented and discussed in the last section. II. THEORETICAL BACKGROUND A. Energetic description of fatigue crack growth The description of kinetics of fatigue failure is very important for estimating fatigue lifetime of a component. The knowledge of crack propagation direction and crack growth rate makes it possible to predict the lifetime by means of Kinetics Fatigue Failure Diagrams (KFFD). However, for an estimation of the cyclic plastic strain energy Wp, the area of the hysteresis loop (Figure 1), which characterizes the energy corresponding to one loading cycle, can be used. To this end, a power law relationship between stress and strain has been proposed by Morrow [13] as follows: 1 1 N " W . N "       (1) where "N is the exponent linking the stress amplitude  and the plastic strain amplitude p . ETASR - Engineering, Technology & Applied Science Research Vol. 3, No. 4, 2013, 488-496 489 www.etasr.com Khelil et al.: Modeling of Fatigue Crack Propagation in Aluminum Alloys Using an Energy… On the other hand, a specific energy, U, is defined as, 2 pW U B da / dN    (2) where B is the specimen thickness and /da dN is the fatigue crack growth rate. Fig. 1. Schematic drawing showing hysteresis loop. The evolution of U as a function of /da dN can be subdivided according to three stages of the KFFD. It should be noted that this diagram can be easily obtained experimentally by measuring crack propagation as a function of the stress intensity factor (see Figure 2). The obtained curve is characterized by three stages which are commonly referred to as stage I, II and III respectively [16, 17]. Fig. 2. A typical fatigue crack growth rate curve (kinetic fatigue failure diagram: KFFD). The relation between R and log K is linear [11, 18, 19], whereas that between /da dN and U can be written as: 4 2 c da A K dN G U   (3) B. The cyclic plastic strain enery– Proposed model Assuming that the energy is primarily dissipated in the plastic zone, a comparison can be made between the measured values and those predicted theoretically by assuming propagation in mode I. Rice [1] and Tracey [21] gave an expression for the equivalent shear strain 3  near the crack tip defined in terms of the amplitude function  NR  as follows: 1 1 ' 0 / ( ( ) / ) N NG R r    (4) where r and  are the polar coordinates at the crack tip, G the shear modulus, 0 0 / 3  the yield stress under pure shear and  NR  can be considered the dominant singular term approximation to the elastic-plastic boundary which depends upon the hardening exponent N’ and is given as a function of  in the normalized form [21]:      20/ /N Nf R K   (5) where  Nf  is a dimensionless function which defines the profile of iso-deformation as a function of polar coordinates at the crack tip. According to Rice [1], the singularity of equivalent cyclic strain can be described by applying the tensile form of (4-5) by simply replacing K by ΔK and σ0 by Δσ0 related to the cyclic stress-strain law. 3p / G    (6) where p and  are the plastic strain and stress amplitude. Consequently the cyclic stress strain curve gives the hardening law for the material at the crack tip [22]: 0 0 N'              (7) where  and  are the equivalent strain and stress amplitude, 0 and 0 0 / 3G    represent the cyclic yield strength of the material. Using the expressions 3    and 0 0 / 3  , equation (4) can be rewritten as follows :   1/1 0 3 N NR G r                  (8) Therefore within the hypothesis given by Rice [1], the equivalent strain amplitude near the crack tip is given by:   1/1 2 0 0 . N Nf K r                   (9)   =2a   = 2 a ap ae a Wp ETASR - Engineering, Technology & Applied Science Research Vol. 3, No. 4, 2013, 488-496 490 www.etasr.com Khelil et al.: Modeling of Fatigue Crack Propagation in Aluminum Alloys Using an Energy… where K is the stress intensity range which can be given as a function of the maximum and minimum stresses and the crack length a as follows:  max minK a     . In order to evaluate the amplitude of the equivalent average strain  Chalant [23] considered an element with a rectangular or circular form located at the crack tip so that: 0 1 S dS S     (10) where S is the surface of the element at the crack tip. In the case of a rectangular element with dimensions D1 and D2 (see Figure 3a), we get [16, 23]: 2 1 1 1 1 1 0 1 1 2 2 0 1 1 2 N ' N ' N 'N ' K D I D I N '                          (1) The terms 1I and 2I are given by:         0 0 1 1 1 1 20 1 1 2 2 1 2 2 N ' N N ' N ' / N N ' f I .d co s f I .d sin                                (12) with 1 2 0 12 D tan D   (13)  (a) (b) Fig. 3. Elements at crack tip: (a) rectangular, (b) circular. The drawing shows the hysteresis loop. In the case of a circular crack tip element with a radius R1 (Figure 3b), the expression of the amplitude of the equivalent average strain is given by [12]: 2 1 1 1 0 1 0 2N ' N ' K R I                       (14) with   12 1 0 / N 'NI f d         (15) The plastic energy throughout the plastic zone is obtained by integrating the plastic energy per surface element given by (1), i.e. 0 1 " ( ) 4 . . 1 " pS p N W PZ dS N       (16) where pS is the surface of a quarter of the plastic zone. Hence 2 0 0 1 4 1 P/ r p N " W ( PZ ) . .r dr d N "          (17) with pr indicates the limit of cyclic plastic zone. If we assume that pr is defined by the distance for which the total equivalent strain is equal to 0 , we have, from (9): 2 0 p N K r f ( )     (18) The integral (17) can be evaluated by substituting the expressions (5-8) for  and p , which gives after simplification [12] :       4 22 0 0 0 0 1 " 2 1 ' . 1 " I N KN W PZ N f d N                   (19) This expression can be rewritten as follows:     0 011 1 pz N" W PZ N ' . .S N"      (20) with   4 22 0 0 2 Ipz N K S f d                (21) indicating simply the surface of the plastic zone. On the other hand, Engerand [24] proosed the Tresca or Von Mises criteria to compute the limit of the plastic zone. It is interesting to note that (19) enables us to express the energy dissipated throughout the plastic zone per unit thickness as a function of 4K , which conforms to the theoretical models given by Klingbeil [7], Mazari et al [15]. and Ranganathan et al [8]. Hence, in order to obtain a similar variation as a function of 4K , we propose to compute the plastic zone surface using the following expressions for the plastic rays:   4 2 2 2 0 1 1 2 3 2 2 2p K r ( PD ) cos sin B                       (22) for plane strain, and: 4 2 2 0 1 ( ) cos 1 3sin 2 2 2p K r PS B                      (23) ETASR - Engineering, Technology & Applied Science Research Vol. 3, No. 4, 2013, 488-496 491 www.etasr.com Khelil et al.: Modeling of Fatigue Crack Propagation in Aluminum Alloys Using an Energy… for plane stress. By using (22) and (23), the new model proposed for the calculation of plastic energy per unit thickness can be written as follows:   0 0. . .M pzMMW PZ S      (24) where M is a constant which depends on the material and the criterion used and MS indicates the surface of the plastic zone determined by expressions (22) or (23). Hence, the total dissipated energy Q in the specimen is given by:  Q W PZ B   (25) III. EXPERIMENTAL DETAILS A. Material and specimen configurations The tests were conducted on two aluminum alloys, the 2024 alloy in the T351 condition and the 7075 alloy in the T7351 condition. The nominal composition and mechanical properties of these alloys are given in Tables I and II respectively. In Table II, the parameters K and n are computed from the relationship of Ludwik [25]: .( )npK  (26) TABLE I. NOMINAL COMPOSITION (IN %) OF THE STUDIED ALLOYS. Alloy Si Fe Cu Mn Mg Cr Zn Ti Al 2024 0.10 0.22 4.46 0.66 1.50 0.01 0.04 0.02 Remain 7075 0.07 0.16 1.52 0.04 2.55 0.20 6.00 0.04 Remain TABLE II. NOMINAL MECHANICAL PROPERTIES OF THE STUDIED ALLOYS Material 2024-T351 7075-T7351 Conventional yield stress at 0.2% of plastic strain 0.2  (MPa) 318 470 Stress at fracture R  (MPa) 524 539 Elongation A% 12.8 11.7 Strength coefficient K(MPa) 652 960.5 Hardening coefficient n 0.104 0.051 The tests were carried out using compact tension specimens (CT) with thicknesses of 10 mm for 2024 and 6 mm for 7075 (Figure 4). The direction of cracking has been taken in the rolling direction. All mechanical tests were conducted using an Instron servohydraulic machine at a typical test frequency of 20 Hz in temperature room with a stress ratio of R=0.5. maxmin /R P P (27) where Pmin and Pmax are the minimum and maximum load in the cycle. Fig. 4. Cyclic test configuration on a compact tension specimen. The crack propagation rate was measured by optical techniques on the polished side of the specimen using a travelling microscope with a precision of 0.01 mm. The total dissipated energy was evaluated by the area locked up inside the recorded hysteresis cycles. The stress intensity range for this geometry is given by Newman [26]:   with . P a K f z z wB w     (28) where w and B are respectively the width and the thickness of the specimen, a, is the crack length, and max minP P P   , is the amplitude applied load. In order to obtain more precision, two functions of compliance  f z were used:  for 0.2 0.3z  [26]:   2 3 4 5 6 4.55 40.32 414.7 1698 3781 4287 2017 f z z z z z z z        (29)  for 0.3 0.7z  [27]:   0 5 1 5 2 5 3 5 4 5 29 6 185 5 655 7 1017 638 9 . . . . . f z . z . z . z z . z      (30) B. Identification of the cyclic plastic strain energy parameters The results obtained after the dentification of the cyclic plastic strain energy parameters are summarized in Table III. TABLE III. CYCLIC PLASTIC STRAIN ENERGY PARAMETERS For 2024-T351 0 0 M 'N "N Integral 4W K  914 0.0111 6.67E-4 0.148 0.078 0.0138 2.92E-13 For 7075-T7351 0 0 M 'N "N Integral 4W K  705 0.00849 1.26E-5 0.166 0.0765 0.0140 4.85E-13 For the particular materials under study, we finally obtain:  For 2024-T351: 13 42 92 10 Joule mW , . K /  (31)  For 7075-T7351: P Pmax Pmin P Time P P B a w ETASR - Engineering, Technology & Applied Science Research Vol. 3, No. 4, 2013, 488-496 492 www.etasr.com Khelil et al.: Modeling of Fatigue Crack Propagation in Aluminum Alloys Using an Energy… 13 44 85 10 Joule mW , . K /   (32) C. Identification of Paris law Fatigue crack growth in the CT specimens was modeled according to the Paris Law, where the FCG rate, da/dN, is described in terms of the stress intensity range K according to the following relationship [28]:   mda C K dN  (33) where C and m are the fatigue crack growth coefficient and exponent, respectively.   y = 6E-09x4.5849 y = 3E-08x3.1987 1.0E-07 1.0E-06 1.0E-05 1.0E-04 1.0E-03 1.0E-02 1 10 100 K (MP a.m^1) d a/ d N ( m m /c y cl e) 2024T351 7075T7351 Fig. 5. Fitting of experimental data with a Paris law Figure 5 shows typical results of the experiments used in constructing plots of log(da/dN) in terms of log(K) for each aluminum alloy. A power law model was fit to the steady state region (stage II) of fatigue crack growth and the Paris Law coefficient (C) and exponent (m) were determined for each specimen that underwent stable fatigue crack growth. A logarithmic scale was used to represent the curve as a straight line. Thus, a linear regression returns the material parameters C and m. By identifying the Paris law parameters according to experimental data's of the two alloys, we found that C=6.0E-9 and m=4.5849 with a slope of 4.58 for 2024-T351. However, for 7075-T7351: C=3.0E-8 and m=3.1987 with a slope of 3.20. It is worth noting that the 2024 aluminum alloy remains as an important aircraft structural material due to its extremely good damage tolerance and high resistance to fatigue crack propagation [29]. However, the 7075-T7351 offers good stress- corrosion cracking resistance [30]. IV. IMPLEMENTATION OF THE PROPOSED APPROACH A computer program has been written in Matlab language. The program was developed to calculate the stress intensity factor, the size of the plastic zone, the crack growth rate, the number of cycles and the cyclic plastic strain energy. Conditions of plane strain and Von Mises criterion were considered. First, we specify CT specimen dimensions (B: thickness, w: width, a0: initial crack length, aEnd: crack length for forced termination), quantities relevant to material properties (E: Young’s modulus, ν: Poisson’s ratio, y: yield stress, C and m: material constants of Paris law and the amplitude of the applied load). Then, the program computes the stress intensity factor range K, the plastic radius rp, the fatigue crack growth rate ,da dN and the energetic parameters as shown by the flaw-chart in Figure 6. Fig. 6. Flow chart of the numerical process. ETASR - Engineering, Technology & Applied Science Research Vol. 3, No. 4, 2013, 488-496 493 www.etasr.com Khelil et al.: Modeling of Fatigue Crack Propagation in Aluminum Alloys Using an Energy… V. RESULTS AND DISCUSSION It is worth noting that the evolution of the cyclic plastic strain energy is directly related to the change of the plastic zone size at the crack tip as described in section (II.B). In what follows, we will first deal with the evolution of the plastic zone. Then, we will analyze in detail the evolution of the cyclic plastic strain energy and the speed of crack propagation. A. Evolution of the plastic zone The purpose of this section is to study in detail the evolution of the plastic zone at the crack tip during the crack growth on the two aluminum alloys. It should be noted that, at the crack tip of a ductile material, the strain fields are significant and leading to a considerable extent of the plastic zone. In this case, the mechanical energy at the crack tip is absorbed by the material in the form of linear defects (dislocations) [20,31].On the other hand, the size of the plastic zone depends not only on the nature of the materials, but also on the intensity of the mechanical energy at the crack tip, the geometry and the size of the crack. Noting that the hardening (maximum consolidation of material) characterized by its rate, can significantly slow the extension of the plastic zone. Figure 7 shows the evolution of the plastic zone for both aluminum alloys (2024 and 7075) using Tracey model (Figure 7a) and the proposed model (Figure 7b). It is found that the size of the plastic zone calculated by the proposed model is about 16 times (for 2024) and 18 times (for 7075) greater than that calculated by the Tracey model. Furthermore, we observe that the 2024 alloy, whose mechanical properties of ductility are higher than 7075, presents the largest size of the plastic zone. We have obtained for:  Tracey model: S(PZ) of 2024 = 2,37  S(PZ) of 7075.  Proposed model: S(PZ) of 2024 = 2,12  S(PZ) of 7075. B. Evolution of the total dissipated energy in the specimen It should be noted that the crack extension which leads to the fracture occurs when the provided energy is sufficient to overcome the material strength. Assuming that this energy is mainly dissipated in the plastic zone, a comparison can be made between the measured values and those predicted theoretically as shown in Figure 8 by assuming propagation in mode I. C. Evolution of the total dissipated energy in the specimen It should be noted that the crack extension which leads to the fracture occurs when the provided energy is sufficient to overcome the material strength. Assuming that this energy is mainly dissipated in the plastic zone, a comparison can be made between the measured values and those predicted theoretically as shown in Figure 8 by assuming propagation in mode I. (a) ‐3,00 ‐2,00 ‐1,00 0,00 1,00 2,00 3,00 ‐1,50 ‐0,50 0,50 1,50 y / ( K /   0 )^ 4 x/(K/0)̂ 4 2024T351 7075T7351 Crack (b) ‐0,30 ‐0,20 ‐0,10 0,00 0,10 0,20 0,30 ‐0,30 ‐0,20 ‐0,10 0,00 0,10 y / ( K /   0 )^ 4 x/(K/0)^4 2024T351 7075T7351 Crack Fig. 7. Evolution of the normalized plastic zone in the case of : (a) the Tracey model and (b) the proposed model for both aluminum alloys. Figure 8 illustrates the evolution of the total dissipated energy according to the amplitude of the stress intensity factor K for both types of aluminum alloys in the case of a constant stress ratio of R=0.5. It can be seen that the theoretical estimates given by the Tracey model are much lower than the experimental measurements. However, a good prediction has been obtained by the proposed model. Noting that for 2024 ETASR - Engineering, Technology & Applied Science Research Vol. 3, No. 4, 2013, 488-496 494 www.etasr.com Khelil et al.: Modeling of Fatigue Crack Propagation in Aluminum Alloys Using an Energy… (Figure 8a), Q calculated by the Tracey model is about 11 to 76 times weaker than the measured Q; the difference is less for high values of ∆K. The same trends have been also found for 7075 (Figure 8b). The gap between the theoretical estimates of Tracey and the measured values ranges from 3, for high values of ∆K, to 10 for the low values of ∆K. This difference can be allotted in a great part to the size of the estimated plastic zone which is much lower than that measured especially for materials that exhibit a high ductility as that of 2024-T351 in comparison with that of 7075-T7531. On the other hand, it may be noted from Figure 8, that the results obtained by the proposed model are in good agreement with experimental data for both types of aluminum alloys. It is also observed that for the same value of ∆K, the total dissipated energy of 2024-T351 is higher than that of 7075-T7531. (a)   1,0E-07 1,0E-06 1,0E-05 1,0E-04 1,0E-03 1,0E-02 1 10 100 K (MP a.m^1/2) Q ( J/ cy cl e) Exp. Result s T racey Model P roposed Model (b)   1,0E-06 1,0E-05 1,0E-04 1,0E-03 1,0E-02 1,0E-01 1 10 100 K (MPa.m^1/2) Q ( J/ cy cl e) Exp. Res ults Tracey Mo del P ro po s ed Mo del Fig. 8. Comparison of measured and estimated dissipated energy per cycle for (a) 2024-T351 and (b) 7075-T7351. VI. RELATIONSHIP BETWEEN Q AND da dN The evolution of the crack growth rate is studied as a function of energetic parameters in order to interpret the crack behavior in various elucidated regimes. Figure 9 shows the evolution of da/dN in terms of the total dissipated energy Q for R=0.5. For both aluminum alloys, the experimental results can be subdivided into two distinct stages as shown in Figure 9. In the case of 2024, the stage I is defined by da/dN210-5 mm/cycle. In this stage we note a strong decrease in crack propagation speed with the total dissipated energy. A logarithmic approximation can be used in this stage. The relationship obtained is given by: y = 0.7346x 1.272 y = 2E-05Ln(x) + 0.0001 1.0E-07 1.0E-06 1.0E-05 1.0E-04 1.0E-03 1.0E-02 0.00001 0.0001 0.001 0.01 0.1 Q (J/cycle) d a/ d N ( m m /c y cl e) Exp . Results Prop osed M odel Fig. 9. Comparison of measured Evolution of da/dN with Q for 2024- T351 y = 0.1123x 0.8382 y = 0.0001Ln(x) + 0.0011 1.0E-05 1.0E-04 1.0E-03 1.0E-02 0.0001 0.001 0.01 0.1 Q (J/cy cle) d a/ d N ( m m /c y cl e) Exp . Results Prop osed M odel Fig. 10. Evolution of da/dN with Q for 7075-T7351. The black lines represent the regression curves of experimental results using power and logarithmic functions ETASR - Engineering, Technology & Applied Science Research Vol. 3, No. 4, 2013, 488-496 495 www.etasr.com Khelil et al.: Modeling of Fatigue Crack Propagation in Aluminum Alloys Using an Energy… Stage II is defined by da/dN10-5 mm/cycle. This region exhibits a stable evolution of the crack growth and a power law can be used as an approximation in this stage: 1 2720 7346 . da . Q dN  (35) In the case of 7075, stage I is defined by da/dN2.10-4 mm/cycle. The evolution of the crack growth in this regime is unstable. If we use a logarithmic function, we can determine a relationship of the form: 0 0001 0 0011 da . Ln( Q ) . dN   (36) Stage II is defined by da/dN 2.10-4 mm/cycle. The obtained law has the following form: 0 83820 1123 . da . Q dN  (37) One can observe that there is a good agreement in the stage II between experimental data and the results obtained by the proposed model for both types of aluminum alloys. Furthermore, the straight lines of slope 0.316 for 2024 and 0.086 for 7075 can be deduced from the proposed model. Indeed, we obtain the following expressions: 1 14620 316 . da . Q dN  for 2024T351 (38) 0.79970.086 da Q dN  for 7075T7531 (39) However a significant difference was found in Stage I between experimental data and the results of the proposed model. This limits, therefore, the application of the proposed model in the stage where the crack propagation is stable. VII. CONCLUSIONS In this study an energy based approach of fatigue crack growth has been proposed. This approach provides a direct link between the cyclic plastic strain energy and the plastic zone at the crack tip. In order to validate the model, experimental data conducted on CT specimens of aluminum alloys (2024-T351 and 7075-T7351), under constant amplitude with a stress ratio of R=0.5 and mode I loading, have been exploited. The following conclusions have been drawn:  The measured values of the cyclic plastic strain energy for both types of aluminum alloys are substantially higher than those calculated by the Tracey model. However, a relatively good agreement has been found between the experimental data and the results obtained by the proposed model.  A correct modeling of the plastic zone is necessary to accurately determine the total cyclic plastic dissipation at the crack tip. The size of the plastic zone calculated by the proposed model is about 16 to 18 times greater than that calculated by the Tracey model. Finally, it should be noted that the simplicity of the current modeling approach limits its ability to account for the crack closure, the environment and the variable amplitude loadings effects, which are topics of ongoing research. ACKNOWLEDGEMENTS It should be noted that the experimental tests were carried out and provided within the framework of a cooperation project "CMEP 389 MDU 97" between the Laboratory of Mechanics and Energetic (University Djillali Liabes of Sidi Bel Abbes (Algeria)) and the Laboratory of Mechanics and Rheology (University François Rabelais of Tours (France)). The authors especially thank Professor N. Ranganathan and his team for their support and assistance. REFERENCES [1] J. R. Rice, “The mechanics of crack tip deformation and extension by fatigue”, Fatigue Crack Propagation Special Technical Publication 415, ASTM, pp. 247-311, Philadelphia, 1967 [2] I. Shozo, I. Yoshito, M. E. Fine, “Plastic work during fatigue crack propagation in a high strength low alloy steel and in 7050 AL-Alloy”, Engineering Fracture Mechanics, Vol. 9, No. 1, pp. 123-136, 1977 [3] P. K. Liaw, S. I. Kwun, M. E. Fine, “Plastic work of fatigue crack propagation in steels and aluminum alloys”, Metallurgical Transactions A, Vol. 12, No. 1, pp. 49-55, 1981 [4] C. Saix, P. Jouanna, “Analyse de la dissipation plastique dans des pièces métalliques minces”, Journal de Mécanique Appliquée, Vol. 5, No. 1, pp. 65-93, 1981 [5] S. R. Bodner, D. L. Davidson, J. Lankford “A description of fatigue crack growth in terms of plastic work”, Engineering Fracture Mechanics, Vol. 17. No. 2, pp. 189-191, 1983 [6] A. D. Joseph, T. S. Gross, “Comparison of techniques for the measurement of plastic work of fatigue crack growth in low carbon steel”, Engineering Fracture Mechanics, Vol. 21, No. 1, pp. 63-74, 1985 [7] N. W. Klingbeil, “A total dissipated energy theory of fatigue crack growth in ductile solids”, International Journal of Fatigue, Vol. 25, No. 2, pp. 117-128, 2003 [8] N. Ranganathan, F. Chalon, S. Meo, “Some aspects of the energy based approach to fatigue crack propagation”, International Journal of Fatigue, Vol. 30, No. 10-11, pp. 1921-1929, 2008 [9] R. Jones, M. Krishnapillai, K. Cairns, N. Matthews, “Application of infrared thermography to study crack growth and fatigue life extension procedures”, Fatigue & Fracture of Engineering Materials & Structures, Vol. 33, No. 12, pp. 871-884, 2010 [10] J. S. Daily, N. W. Klingbeil, “Plastic dissipation energy at a bimaterial crack tip under cyclic loading”, International Journal of Fatigue, Vol. 32, No. 10, pp. 1710-1723, 2010 [11] J. Weertman, “Theory of fatigue crack growth based on a BCS crack theory with work hardening”, International Journal of Fracture.; Vol. 9, No. 2. 125-131, 1973 [12] M .Mazari, “Contribution à l’étude d’une approche énergétique de la propagation des fissures de fatigue”, Thèse de doctorat, Université de Sidi Bel Abbès, Algérie, 2003 [13] J. Morrow, “Cyclic Plastic Strain Energy and Fatigue of Metals” in Internal friction, damping, and cyclic plasticity, ASTM STP 378, American Society for Testing and Material, 1965 [14] J. S. Daily, N. W. Klingbeil. “Plastic dissipation in fatigue crack growth under mixed mode loading”, International Journal of Fatigue, Vol. 26, No. 7, pp.727-738, 2004 [15] M. Mazari, B. Bouchouicha, M. Zemri, M. Benguediab, N. Ranganathan, “Fatigue crack propagation analyses based on plastic energy approach”, Computational Materials Science, Vol. 41, No. 3, pp. 344-349, 2008 [16] N. Ranganathan, “Contribution au développement d’une approche énergétique à la propagation d’une fissure de fatigue”, Thèse de doctorat, Université de Poitiers. France, 1985 ETASR - Engineering, Technology & Applied Science Research Vol. 3, No. 4, 2013, 488-496 496 www.etasr.com Khelil et al.: Modeling of Fatigue Crack Propagation in Aluminum Alloys Using an Energy… [17] W. Wang, H. Hsu, “Fatigue crack growth rate of metal by plastic energy damage accumulation theory”, Journal of Engineering Mechanics, Vol. 120, No. 4, pp. 776-795, 1994 [18] S. M. Beden, S. Abdullah, A. K. Ariffin. “Review of Fatigue Crack Propagation Models for Metallic Components”. European Journal of Scientific Research, Vol. 28. No. 3, pp. 364-397, 2009 [19] R. O. Ritchie, “Mechanisms of fatigue-crack propagation in ductile and brittle solids”. International Journal of Fracture., Vol. 100, No. 1, pp. 55- 83, 1999 [20] N. Ranganathan, K. Jendoubi, M. Benguediab, J. Petit, “Effect of R ratio and K level on the hysteretic energy dissipated during fatigue crack propagation”, Scripta Metallurgica, Vol. 21, No. 8, pp. 1045-1049, 1987 [21] D. M. Tracey, “Finite element solution for crack-tip behavior in small- scale yielding”, Journal of Engineering Materials and Technology, Vol. 98, No. 2, pp. 146-151, 1976 [22] G. Chalant, L. Remy, “Plastic strain distribution at the tip of a fatigue crack. Application to fatigue crack closure in the threshold regime”, Engineering Fracture Mechanics, Vol. 16, No. 5, pp. 707-720, 1982 [23] G. Chalant, “Fissuration par fatigue d’alliages cobalt-Nickel : Discussion d’un modèle mécanique de propagation”. Thèse de doctorat de l’Ecole des Mines de Paris, 1981 [24] J. L. Engerand, Mécanique de la rupture, Ed. Techniques Ingénieur, 1990 [25] P. Ludwik, Elemente der Technologischen Mechanik, Springer-Verlag OHG, Berlin, 1909 [26] J. C. Newman, “Stress analysis of the compact specimen including the effects of pin loading fracture analysis”, ASTM STP 560, pp. 105-121, 1974 [27] J. E. Srawley, B. Gross, “Stress intensity factors for bend and compact specimens”, Engineering Fracture Mechanics, Vol. 4, No. 3, pp. 587- 589, 1972 [28] P. C. Paris, F. A. Erdogan, “A critical analysis of crack propagation laws”, Journal of Basic Engineering, Vol. 85, No. 4, pp. 528-533, 1963 [29] D. Altenpohl, Aluminium: Technology, application and environment. A profile of a modern metal: Aluminum from Within, 6th Edition, Wiley, 2010 [30] P. S. Pao, S. J. Gill , C. R. Feng, “On fatigue crack initiation from corrosion pits in 7075-T7351 aluminum alloy”, Scripta Materialia, Vol. 43, No. 5, pp. 391-396, 2000 [31] J. W. Kysar, “Energy dissipation mechanisms in ductile fracture”, Journal of the Mechanics and Physics of Solids, Vol. 51, No. 5, pp. 795- 824, 2003