Preprint
Article

This version is not peer-reviewed.

Effect of Pre-Strain on Fatigue Crack Growth in AA2024-T351

Submitted:

31 January 2025

Posted:

03 February 2025

You are already at the latest version

Abstract
The objective here is to study the effect of pre-strain on fatigue crack growth (FCG) in 2024-T351 aluminium alloy. Three pre-strain conditions were considered: without pre-strain, compressive and tensile permanent pre-strains of 4% . A numerical approach based on cumulative plastic strain at the crack tip was followed to predict FCG rate. The compressive pre-strain increased FCG rate, while the tensile pre-strain reduced the da/dN relatively to the conditions without pre-strain. A strong link was established between the influence of pre-strain and plasticity induced crack closure. This relevant conclusion was based on the linear trend obtained between da/dN and Keff for three crack lengths (a=16.184; a=15.048 mm and a=15.152 mm) and three the pre-strain conditions. The increase of stress ratio from R=0.1 to R=0.5 and the elimination of the contact of crack flanks reduced significantly the effect of pre-strain, also pointing for the huge relevance of crack closure in this context. The effect of pre-strain on da/dN after an overload was also explained by crack closure variations.
Keywords: 
;  ;  ;  ;  

1. Introduction

Real structural components are submitted to plastic deformation during fabrication. Operations like blanking, trimming, stretching, bending and deep drawing are involved in structure fabrication. Sheet metal parts for automotive industry are produced by stamping [1], involving large plastic deformation. Material deformation can also be done to improve the mechanical properties, like tensile strength and yield stress [2,3]. The plastic deformation resulting from fabrication affects material behaviour [4,5], therefore must be included in the design of components and structures.
In laboratory testing or numerical simulation, the effect of pre-strain is studied loading a material up to the desired amount, so that it can experience deformation in the plastic strain regime, and then unload the material. The pre-strain increases yield stress and tensile strength, while decreases ductility. In other words, the material becomes less ductile with pre-strain. This was observed by Ghosal et al. in dual-phase (DP590) steel [6,7], by Le et al. [8] in DP600 steel, by Fredriksson et al. [9] in DP400, DP600, HSLA steel (HSLA500) and a deep-drawing quality steel. Similar trends were also found by Wang et al. [10] for TWIP steel, Walker et al. [11] for DP and FB steels, Ji et al. [12] for different grades of DP steels, Wu et al. [13] for stainless steel (SUH660). These effects of pre-strain are explained by a substantial increase in dislocations [6], leading to localized alterations in lattice orientations. The increased dislocation density of the pre-strain material leads to significant work hardening.
In the high cycle fatigue regime (HCF), the fatigue strength can be defined as the particular stress amplitude below that a material can endure an infinite number of cycles without failure. Ghosal et al. [6] and Le et al. [7] observed an increase of endurance limit with pre-strain in dual-phase steel. Wang et al. [10] studied AA2024-T3 and found a decrease in fatigue strength. The influence of tensile pre-strain on high-cycle fatigue life is therefore material-dependent.
The effect of pre-strain on low-cycle fatigue (LCF) has also been widely studied. Le Roux et al. [14] found a decrease in the uniaxial pre-strain 304L stainless steel’s fatigue life compared to the as-received condition under strain-controlled cycling. The decrease in the specimen’s fatigue life was attributed to the reduction in ductility after the pre-straining operation. Gustavsson et al. [15] also observed a reduction in fatigue life of uniaxial pre-strain DP steel, due to the presence of tensile mean stress under fully reversed strain-controlled fatigue loading. Parker et al. [16] investigated the effect of 40% balanced biaxial stretching on hot rolled low carbon steel SAE1008. They have reported that fatigue resistance starts to degrade in the LCF region compared to the as-received condition due to the cyclic softening nature of the 40% balanced biaxial stretched low carbon steel. The pre-straining decreases ductility, which can cause poor fatigue performance in the LCF region. In contrast, it was found that fatigue resistance improves in the HCF region. The same trends were observed in DP600 steel submitted to 12.5% uniaxial pre-strain [17]. Froustey et al. [18] studied AA5454 and AA2017. The AA2017 showed a considerable reduction in fatigue life, but AA5454 revealed unaffected fatigue properties as increasing the tensile pre-straining magnitudes. Rex et al. [19] studied the 2024-T4 aluminium alloy. A considerable increase in fatigue life is found for uniaxial tensile pre-strained than the 0% prestrained condition at higher strain amplitude. Furthermore, equi-biaxial tensile pre-strained specimens showed a quite small increment in fatigue life compared to 0% pre-strained material. As-received specimens exhibit noticeable cyclic hardening at all the strain amplitude values. The pre-strained specimen shows an almost stable cyclic stress–strain response for each strain amplitude. Branco et al. [20] studied the effect of 4 and 8% tensile pre-strains on LCF in AA7050-T6. They observed a reduction of fatigue life with pre-strain. Without pre-strain the dislocations tended to be distributed in a more homogeneous manner, while under higher tensile pre-strains, the dislocations were mainly observed at the slip bands.
The effect on fatigue crack growth (FCG) has also been studied. Schijve [21] studied aluminium alloys and observed a substantial rise in FCG rate of 2024-T3 alloy after applying a 3% tensile pre-strain. He suggested that the elevated yield strength due to pre-straining leads to a smaller plastic zone size, intensifying tensile stresses and diminishing crack closure within the plastic zone near the crack tip, thereby accelerating crack propagation. Al-Rubaie et al. [22] studied AA7475-T7351 and found a limited effect of pre-straining in regimes I and II, but an effect in regime III associated with a loss of fracture toughness. Anandavijayan et al. [23] also found a limited effect in S355 steel submitted to pre-strains of 0, 5 and 10%. On the contrary, a great influence was found on high-cycle fatigue life. Tai Shan and Liu [24] observed an increase of FCG rate in AA2024-T351. Wasen et al. [25] found that pre-straining reduces FCPR in dual-phase steel. Arora et al. [26] studied mild steel and observed a decrease of FCG rate in regime II. Radhakrishnan et al. [27] investigated the influence of pre-straining and stress ratio on FCG in mild and stainless steel. They observed that the FCG rate decreased by factors of 2.4 and 3.7 for 6% and 9% pre-strain and 1.4 and 2.4 for 9% and 14% pre-strain in mild and stainless steel, respectively. In contrast to earlier findings, Kim et al. [28] showed that FCG rate increases with higher pre-strain levels in the lower and intermediate ΔK regimes in PH steel. Leitner et al. [29] also found that increased plastic pre-deformation increases the FCG rate and ΔKth in pearlitic steel. The conclusions of the literature mentioned above demonstrate that the plastic pre-straining's influence on FCG response is primarily material-dependent. There is a lack of understanding of the mechanisms explaining the effect of pre-strain. Besides, the effect of compressive pre-strain on FCG was not reported.
The objective of this work is to study the effect of pre-strain on FCG in 2024-T351 aluminium alloy. A numerical was built to simulate the application of pre-strain and the subsequent propagation of a fatigue crack by node release. Tensile and compressive pre-strains of 4% were studied.

2. Numerical Model

Figure 1a shows the CT specimens considered in the numerical study. These specimens had a width W=36 mm and an initial crack length a0=15 mm. Only 1/4 of CT specimen was modelled, as indicated in Figure 1b and Figure 1c, considering adequate boundary conditions. The contact of crack flanks was simulated using a rigid plane placed at the crack symmetry plane. A small thickness of 0.1 mm was considered for the specimen, in order to model a plane stress state. Simulations without contact between crack flanks were also run, to understand the effect of crack closure.
The material studied was the 2024-T351 aluminium alloy. The 2XXX series of aluminium alloys are widely used in aerospace applications owing to their high strength, good fracture toughness, and low density [30]. As a result of its high strength-to-weight ratio, the Al–Cu–Mg–Mn based AA2024 aluminium alloy is widely used for various engineering applications, including aircraft and transportation applications [31,32]. Typically, this material is used in the design of fuselage, wings, and lightweight components in the aerospace industry. Table 1 shows the chemical composition of the 2024-T351 aluminium alloy used in this work. Table 2 presents the material properties of the 2024-T351 aluminium alloy without pre-strain. Young modulus, Yield stress, Ultimate Tensile Strength (UTS) and elongation at break were determined with a tensile test.
The accurate modelling of material elastic-plastic behaviour is fundamental to obtain good quality numerical predictions of crack tip plastic deformation. The elastic behaviour is assumed isotropic and described by Hooke’s law. On the other hand, the plastic behaviour is characterised by the von Mises yield criterion and the Swift isotropic hardening law coupled with Lemaître-Chaboche kinematic hardening law under an associated flow rule. Swift hardening law [33] is described by:
Y = C Y 0 C 1 n + ε - p n
where Y0, C and n are the material parameters and ε - p denotes the equivalent plastic strain. The Lemaître-Chaboche kinematic hardening law [34] is:
X ˙ = C x X Sat σ ' - X σ ̄ - X   ε ̄ ˙ p
where CX and XSat are the material parameters of Lemaître-Chaboche law, σ ̄ is the equivalent stress and ε ̄ ˙ p is the equivalent plastic strain rate. The material parameters were identified using LCF fatigue curves obtained for Rε =-1 and Δε/2=1.5% [35] and are presented in Table 3.
The numerical model was implemented using the in-house software DD3IMP (Three-Dimensional Elasto-Plastic Finite Element Code), originally developed for the simulation of deep-drawing [36]. An updated Lagrangian approach is used to describe the evolution of the deformation process, assuming a hypoelastic-plastic model. The numerical model considers large elastoplastic strains and rotations, while the elastic strains are assumed negligibly small. Figure 2 shows the finite element mesh considered, composed of linear isoparametric elements. It was refined in the crack growth region, having elements with 8×8 μm2. Only one layer of elements was considered along the thickness direction.
Crack propagation was numerically modelled by successive debonding of both current crack front nodes. The release occurs at minimum load to avoid convergence problems that could arise by propagating the crack at maximum load. The size of each crack increment corresponded to the finite element size in the refined region. The two crack front nodes were released when the accumulation of plastic strain reached the critical value of strain. This critical value was calibrated using one experimental value of da/dN, and a value of 110% was obtained in a previous work of the authors [37].
The pre-strain was applied by displacing a range of nodes along the entire width of the specimen. This displacement simulated what would experimentally be a distributed tensile or compressive load on the specimen. The region where plastic deformation was induced is indicated by the red rectangle in Figure 3. The region above it is free from pre-deformation. The displacements were adjusted to obtain plastic deformations of 4% in tension and compression (εps=-4, +4%). After the application of pre-deformation, a crack with a length a0=15 mm is simulated and submitted to cyclic loading and crack propagation by node release.
The CT specimen was submitted to constant amplitude loading with stress ratios R=0.1 and R=0.5. In the first case, the maximum and minimum forces were 41.67 N e 4.17 N, being 75 N and 37.5 N for a stress ratio R=0.5. Simulations of overloads were also considered. In the overload cycle the maximum and minimum loads were 79.17 N and 4.17, respectively, giving an overload ratio OLR=1.5, calculated according:
O L R = K O L K m i n K m a x K m i n
The baseline loading had Fmax=41.67 N and Fmin= 4.17 N, respectively. Table 4 presents all the load cases considered. Simulations without contact of crack flanks were also considered, as already mentioned.

3. Results

3.1. Effect of Pre-Strain on da/dN

Figure 4 shows the variation of FCG rate versus crack length, a, with and without pre-strain. The crack propagation started at a0=15 mm and ended at a=16.2 mm.
The specimen without pre-strain presents an initial small increase of da/dN, followed by a progressive decrease to a minimum value, which occurs at about a=15.5 mm. This decrease is explained by the progressive increase of crack closure level with the formation of residual plastic wake. Above this transient regime there is a stable regime characterised by a progressive increase of da/dN, which is explained by the increase of crack tip fields with crack length.
The specimen with a tensile pre-strain of 4% presents a similar behaviour, but the variations are much smaller. There is a small increase of da/dN, followed by a smooth decrease to a minimum value which occurs at about a=15.2 mm. Finally, there is a progressive increase with crack propagation as a result of the increase of crack length.
The specimen submitted to a compressive pre-strain of 4% presents a small initial increase of da/dN, but without subsequent decrease. The slope of increase of da/dN with crack length is similar to that observed for the specimen without pre-strain. On the other hand, the specimen with tensile pre-strain has a lower rate of variation of da/dN in the stable regime. The FCG rate in the stable regime is higher for the compressive pre-strain, intermediate for the specimen without pre-strain, and lower for the tensile pre-strain. For a crack length of 16.184 mm, indicated by the vertical dashed line, da/dN has values of 0.95, 0.71 and 0.39 μm/cycle for εps=-4, 0 and +4%, respectively.
Figure 5 presents the CTOD versus load plots, which are particularly interesting to understand what is happening at the crack tip. The CTOD curves plotted were measured at the first node behind crack tip for a crack length a=16.184 mm, indicated by the vertical dashed line on the right-hand side of Figure 4. The load ranged between Fmin=4.17 N and Fmax=41.67 N, as can be seen in the horizontal axis. The crack is closed at the lower part of the load cycle, having a linear elastic behaviour after opening. At relatively high loads there is a departure from linearity, resulting from crack tip plastic deformation. The elastic and plastic CTOD ranges (δe and δp, respectively), are indicated for the compressive pre-strain.
The crack opening level is lower for the compressive pre-strain (Fopen≈10 N), and higher for the tensile pre-strain (Fopen≈20 N), having an intermediate value for the specimen without pre-strain (Fopen≈15 N). The crack closure level was quantified using:
U * = F m a x F o p e n F m a x F m i n
The compressive pre-strain has U*=18.23%, the situation without pre-strain has U*=27.15%, and the tensile pre-strain has U*=42.53%. The increase of crack opening level decreases the effective load range, which reduces maximum CTOD, the elastic and plastic deformation ranges and the area of the loop (which is proportional of the dissipated energy). The sizes of monotonic plastic and reversed plastic zones are also expected to decrease. Therefore, a relation may be established between these non-linear parameters and da/dN values presented in Figure 4.
The relation between FCG rate and the effective stress intensity factor range, ΔKeff, was studied for three crack lengths (a=16.184; a=15.048 mm and a=15.152 mm) and three pre-strain conditions (εps=-4, 0, +4%). Figure 6 shows the results obtained, which indicate a well-defined correlation between da/dN and ΔKeff, although some scatter can be observed. The increase of ΔKeff increases da/dN, and a linear trend is obtained. This correlation clearly indicates that plasticity induced crack closure (PICC) is behind the variations of da/dN observed in Figure 4 and therefore behind the effect of pre-strain on FCG.

3.2. Influence of Stress Ratio on Pre-Strain Effect

Figure 7 presents results of da/dN versus crack length for a stress ratio R=0.5 and different pre-strain values (εps=-4, 0, +4%). The dashed lines indicate results obtained for R=0.1. The increase of R increased significantly da/dN, which means that the stress ratio has a major impact. The effect of R is usually associated with crack closure phenomenon, therefore the results in Figure 7 indicate that crack closure is affecting the results. The increase of FCG rate with stress ratio has been observed by different authors, namely Mehrzadi e Taheri [38] in a magnesium alloy and Seifi and Hosseini [39] in pure copper.
On the other hand, the effect of pre-strain is less evident for R=0.5 than it was R=0.1. Note that linear scales are being used, instead of logarithmic scales. Anyway, the tensile pre-strain produces a slightly higher da/dN for the compressive pre-strain, while the pre-strain εps=0 has intermediate values. The slope, i.e., the rate of variation of da/dN with crack length, is similar for the different pre-strains. So, there is an inversion of the effects of tensile and compressive pre-strains, relatively to the results obtained for R=0.1.
Some oscillations can be observed in all curves. In fact, the node release for the simulation of crack growth occurs at minimum load, when the accumulated plastic strain reaches the critical value (110%). Normally this critical value is reached in the middle of a load cycle. So, there is some inaccuracy in the predicted value of N, i.e., the number of load cycles between node release used to calculate da/dN. The problem becomes more relevant when N is relatively small, i.e., when the crack growth rate is relatively high. In this case, although the error is just a fraction of unit, it has some impact. The results of Figure 4 also show some oscillations, but with low amplitude because N is relatively high.
Figure 8 plots the variation of CTOD with applied load for a stress ratio R=0.5 and different pre-strain values (εps=-4, 0, +4%). Both the maximum and minimum loads increased relatively to the situation with R=0.1, which is represented by the dashed lines. The curves for R=0.5 are very similar, showing a limited effect of pre-strain. The compressive pre-strain presents a slightly lower value of opening load. The comparison with the dashed lines, indicate a great decrease of crack opening level, which has values of U* in the range 6-11%, while for R=0.1 it ranged from 18 to 40%. This reduction of U* increases the effective load range, which increases significantly the total CTOD range and the elastic and plastic components, δe e δp, respectively. The variation of U* also explains the increase of da/dN produced by the increase of R and the reduced effect of pre-strain at R=0.5. In fact, since the effect of pre-strain is associated with crack closure, and there is almost no crack closure at R=0.5, a minor effect of pre-strain on FCG may be expected.

3.3. Effect of Pre-Strain on Variable Amplitude Loading

The influence of pre-strain on the effect of an overload was also studied. The overload was applied at the load cycle number 1500, at a baseline stress ratio R=0.1 and considering an overload ratio of 1.5. Figure 9 compares the results of da/dN obtained without pre-strain and with a tensile pre-strain of 4%. These tests took a long time because FCG rate drops significantly after the overload.
The variation of da/dN is according to the typical behaviour observed after an overload [40,41]. Immediately after the overload, there is a sudden increase of da/dN followed by a fast decrease to a minimum value. Some crack extension is required to reach this minimum, which is named delayed retardation. After the minimum there is a progressive increase of FCG rate to the stable value, corresponding to constant amplitude loading. The increase of da/dN is due to crack tip blunting, which eliminates previous crack closure. The variations after blunting are due to the formation of a new residual plastic wake, stronger than that corresponding to constant amplitude loading. In fact, the overload increases the monotonic plastic zone, and this has an impact on crack closure, which is felt after some crack propagation. The two main mechanisms of plasticity induced crack closure are clearly visible in the overloads, i.e., the crack tip blunting and the effect of residual plastic wake.
Comparing the curves for the two pre-strain values (εps=0, +4%), significant differences are evident. The overload is applied at a shorter crack length for εps=+4% because the crack growth rate is lower while the number of load cycles is the same (1500). Besides, the variations of da/dN are much more pronounced for the case without pre-strain. The extent of propagation affected by the overload is also much higher for εps=0.
For a better understanding of the effect of pre-strain, Figure 10 plots CTOD curves. Figure 10a presents CTOD variation exactly in the overload cycle. For relatively low loads the crack is closed, i.e., the CTOD is zero. The crack opening occurs first for εps=0, as already observed in Figure 5. After crack opening, the CTOD increases first linearly and then non-linearly with the load (region II). The linear variation is a result of elastic deformation, while the non-linear variation is due to plastic deformation. The increase of the load above the baseline values, i.e., above Fmax,BL=41.67 N, produces a sudden change of the slope of the curve (region III). The material has not experienced such high loads before, therefore there is a strong plastic deformation, significantly higher than that observed in region II. Consequently, CTOD reaches relatively high values. The subsequent decrease of the load from FOL to Fmin reduces the blunting (region IV), initially linearly and then non-linearly. The CTOD at the minimum load is the crack tip blunting produced by the overload.
Figure 10b shows the CTOD curves at the minimum values of da/dN, which occur at crack lengths of 16.336 mm and 15.448 mm for εps=0 and +4%, respectively. The crack closure is quite relevant, particularly for the tensile pre-strain, much more than observed under constant amplitude loading (Figure 5). U*=53.1% for the specimen without pre-strain and U*=71.0% for the tensile pre-strain, at the minimum da/dN after the overload. These high values of U* reduce the effective load range, and therefore the crack tip plastic deformation. This is particularly dramatic for εps=+4%, because there is no evidence of plastic CTOD, δp. This explains the very small value of da/dN in Figure 9. The slopes of the elastic regimes are very similar, which means that the elastic behaviour is not affected by the overload, as could be expected. Therefore, the minimum da/dN is linked to a peak of crack closure, as was observed by Baptista et al. [42].

4. Discussion

The first issue is the validity of the conclusions obtained using only numerical simulation. Naturally, it would be interesting to have experimental results to validate the numerical predictions. However, this is expensive, since it requires the preparation of specimens with pre-strain and their testing, and time consuming. On the other hand, the present numerical model has already proved to be competent and robust in the simulation of FCG. In fact, the numerical model based on cumulative plastic strain was able to predict the effect of ΔK [37], stress ratio [43], Superblock2020 load pattern [44], overloads [45], load blocks [46], and the effect of lateral holes [47]. These successive validations are also a clear indication that cyclic plastic deformation is the main mechanism responsible for FCG. Besides, numerical approaches are very interesting to develop parametric studies and to analyze the fundamental mechanisms. This is particularly important in the study of the effect of pre-strain, where the mechanisms are not completely understood.
The present results clearly indicate that plasticity induced crack closure is behind the effect of pre-strain. The decrease of ductility is expected to decrease the size of monotonic and reversed plastic zones. This reduces residual plastic deformations and therefore crack closure. On the other hand, there is a decrease of crack blunting which potentiates the effect of residual plastic wake. These mechanisms of opposite effect may explain the different trends observed for different materials, or, in other words, the dependence on material. Schijve [21] studied AA2024-T3 submitted to a pre-strain of 3%. He also linked the increase of da/dN with crack closure phenomenon.
The material hardening resulting from plastic deformation may also be expected to have a direct effect on crack tip damage. In fact, the cumulative damage will be lower in the hardened material. Therefore, a numerical study was developed without contact of crack flanks, which eliminates crack closure phenomenon from the results. Figure 11a presents da/dN versus crack length. The dashed lines indicate the predictions obtained with contact of crack flanks also for R=0.1. As can be seen, there is some effect of pre-strain, but significantly less than observed in the presence of contact of crack flanks. The tensile pre-strain now has a damaging effect, relatively to the situation without pre-strain. On the other hand, the compressive pre-strain reduces da/dN. Figure 11b presents the CTOD curves versus load for a crack length a=16.048 mm. Negative values of CTOD can be observed, which means that the crack flanks are overlapping. The elimination of the contact is not physically correct, but is particularly interesting to isolate crack tip phenomena, from the influence of crack closure. The applied load range is now totally effective. In other words, the crack tip feels all load cycle without any protection from crack closure phenomenon. This produces an increase of elastic, plastic and total CTOD ranges.
The inversion of the effect of tensile and compressive pre-strain relatively to the situation with contact (Figure 4) is curious and deserved a closer look. Figure 12 shows the variation of accumulated plastic strain at the crack tip with time (i.e., with load cycling). There are two aspects that deserve special attention: the slope of the curves between node releases and the initial value of plastic strain after node release. The slope represents the rate of accumulation of damage, while the initial value is a consequence of the accumulation of damage in previous positions of crack tip. Figure 12a presents the results obtained for five load increments starting at a=15.952 mm, with contact of crack flanks activated. The slopes are clearly different, being higher for the compressive pre-strain. The initial value of accumulated strain is also higher for εps=-4%, which indicates larger monotonic and reversed plastic zones. Consequently, assuming 70 load cycles, there are 5 node releases for the compressive pre-strain while the tensile pre-strain only suffers two node releases. These trends are according to the results of Figure 4, i.e., explain the higher da/dN obtained for compressive strain. Figure 12b presents similar results, but without contact of crack flanks. Now, the tensile pre-strain has slightly higher values of initial plastic strain and rate of accumulation of plastic strain. This explains the slightly higher value of da/dN obtained for εps=+4% in Figure 11a.

5. Conclusions

The main objective here is to study the effect of pre-strain on FCG in AA2024-T351. A numerical model was created to predict FCG assuming that cyclic plastic deformation is the main damage mechanism and using the cumulative plastic strain at the crack tip as driving parameter. Two values of pre-strain were considered: a compressive pre-strain of 4%, a tensile pre-strain of 4%, which were compared with the situation without pre-strain (i.e., εps=-4, 0, +4%).
After stabilization of FCG rate, da/dN is higher for the compressive pre-strain, intermediate for the specimen without pre-strain, and lower for the tensile pre-strain. Values of U* of 18.23, 27.15 and 42.53% were obtained for εps=-4, 0 and +4%, respectively. The increase of crack opening level decreases the effective load range. A well-defined correlation was obtained between da/dN and ΔKeff for three crack lengths (a=16.184; a=15.048 mm and a=15.152 mm) and the three pre-strain conditions (εps=-4, 0, +4%). This clearly indicates that crack closure is explaining the variations of da/dN produced by the pre-strain. The increase of stress ratio from R=0.1 to R=0.5 reduced significantly the relevance of pre-strain, which is also explained by the reduction of crack closure phenomenon. The effect of pre-strain on da/dN after an overload was also explained by crack closure variations.

Author Contributions

Conceptualization, D.N. and F.A.; methodology, D.N.; formal analysis, A.A. and E.S. ; investigation, A.A.; data curation, A.A.; writing—original draft preparation, F.A.; writing—review and editing, D.N. and E.S.. All authors have read and agreed to the published version of the manuscript.

Funding

The authors gratefully acknowledge the financial support from the Portuguese Foundation of Science and Technology (FCT) under the projects UIDB/00285/2020 and LA/P/0112/2020. Edmundo Sérgio is also grateful to FCT for the PhD grant with reference 2022.11438.BD.

Data Availability Statement

Dataset available on request from the authors.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AA Aluminium alloy
CA Constant amplitude
CT Compact tension
CTOD Crack tip opening displacement
FCG Fatigue crack growth
OL Overload
OLR Overload ratio

References

  1. Hariharan, K.; Prakash, R.V. Integrating effect of forming in fatigue life prediction: review of present scenario and challenges. Proc. Inst. Mech. Eng. Part B: J. Eng. Manuf. 2012, 226, 967–979. [Google Scholar] [CrossRef]
  2. Wang, B.; Zhang, P.; Duan, Q.; Yang, H.; Pang, J.; Tian, Y.; Li, X.; Zhang, Z. High-cycle fatigue properties and damage mechanisms of pre-strained Fe-30Mn-0.9C twinning-induced plasticity steel. Mater. Sci. Eng. A 2017, 679, 258–271. [Google Scholar] [CrossRef]
  3. Zhou, Y.; Lin, M.; Liu, C.; Wang, L.; Chen, H.; Dan, C.; Ma, S.; Chen, Z.; Wang, H. Enhancing mechanical properties of uniformly distributed nano TiB2/2024 Al composite rolling sheet by pre-stretch aging. J. Alloy. Compd. 2022, 913. [Google Scholar] [CrossRef]
  4. Baek, J.-H.; Kim, Y.-P.; Kim, C.-M.; Kim, W.-S.; Seok, C.-S. Effects of pre-strain on the mechanical properties of API 5L X65 pipe. Mater. Sci. Eng. A 2010, 527, 1473–1479. [Google Scholar] [CrossRef]
  5. Mehmanparast, A.; Davies, C.M.; Dean, D.W.; Nikbin, K. Effects of plastic pre-straining level on the creep deformation, crack initiation and growth behaviour of 316H stainless steel. Int. J. Press. Vessel. Pip. 2016, 141, 1–10. [Google Scholar] [CrossRef]
  6. Ghosal, P.; Raj, A.; Paul, S.K. Influence of uniaxial and biaxial pre-straining on the low cycle fatigue performance of DP590 steel. Int. J. Fatigue 2021, 149, 106260. [Google Scholar] [CrossRef]
  7. Ghosal, P.; Paul, S.K.; Raj, A. Influence of uniaxial and biaxial pre-straining on the notch fatigue performance of DP590 steel. Theor. Appl. Fract. Mech. 2021, 115, 103072. [Google Scholar] [CrossRef]
  8. Le, Q.; Kang, H.T.; Kridli, G.; Khosrovaneh, A.; Yan, B. Effect of prestrain paths on mechanical behavior of dual phase sheet steel. Int. J. Fatigue 2009, 31, 607–615. [Google Scholar] [CrossRef]
  9. Fredriksson, K.; Melander, A.; Hedman, M. Influence of prestraining and ageing on fatigue properties of high-strength sheet steels. Int. J. Fatigue 1988, 10, 139–151. [Google Scholar] [CrossRef]
  10. Wang, B.; Zhang, P.; Duan, Q.; Yang, H.; Pang, J.; Tian, Y.; Li, X.; Zhang, Z. High-cycle fatigue properties and damage mechanisms of pre-strained Fe-30Mn-0.9C twinning-induced plasticity steel. Mater. Sci. Eng. A 2017, 679, 258–271. [Google Scholar] [CrossRef]
  11. Walker, J.; Thomas, D.J.; Gao, Y. Effects of shot peening and pre-strain on the fatigue life of dual phase Mar-tensitic and Bainitic steels. J Manuf Processes 2017, 26, 419–424. [Google Scholar] [CrossRef]
  12. Ji, D.; Zhang, M.; Zhu, D.; Luo, S.; Li, L. Influence of microstructure and pre-straining on the bake hardening response for ferrite-martensite dual-phase steels of different grades. Mater. Sci. Eng. A 2017, 708, 129–141. [Google Scholar] [CrossRef]
  13. Wu, H.; Hamada, S.; Noguchi, H. Pre-strain effect on fatigue strength characteristics of SUH660 plain specimens. Int. J. Fatigue 2013, 55, 291–298. [Google Scholar] [CrossRef]
  14. Le Roux, J.-C.; Taheri, S.; Sermage, J.-P.; Colin, J.; Fatemi, A. Cyclic Deformation and Fatigue Behaviors of Stainless Steel 304L Including Mean Stress and Pre-Straining Effects. ASME 2008 Pressure Vessels and Piping Conference. LOCATION OF CONFERENCE, United StatesDATE OF CONFERENCE; pp. 411–420.
  15. Gustavsson, A.; Melander, A. Variable-amplitude fatigue of a dual-phase sheet steel subjected to prestrain. Int. J. Fatigue 1994, 16, 503–509. [Google Scholar] [CrossRef]
  16. Parker, T.E.; Montgomery, G.L. Effect of Balanced Biaxial Stretching on the Low Cycle Fatigue Behavior of SAE 1008 Hot Rolled Low Carbon Steel. 1975 Automotive Engineering Congress and Exposition. LOCATION OF CONFERENCE, COUNTRYDATE OF CONFERENCE;
  17. Das, B.; Singh, A.; Paul, S.K. Low cycle fatigue performance of DP600 steel under various pre-straining paths. Int. J. Fatigue 2019, 132, 105331. [Google Scholar] [CrossRef]
  18. Froustey, C.; Lataillade, J.L. Influence of large pre-straining of aluminium alloys on their residual fatigue resistance. Int. J. Fatigue 2008, 30, 908–916. [Google Scholar] [CrossRef]
  19. Rex, A.V.; Paul, S.K.; Singh, A. The influence of equi-biaxial and uniaxial tensile pre-strain on the low cycle fatigue performance of the AA2024-T4 aluminium alloy. Int. J. Fatigue 2023, 173. [Google Scholar] [CrossRef]
  20. Branco, R.; Costa, J.; Borrego, L.; Wu, S.; Long, X.; Antunes, F. Effect of tensile pre-strain on low-cycle fatigue behaviour of 7050-T6 aluminium alloy. Eng. Fail. Anal. 2020, 114, 104592. [Google Scholar] [CrossRef]
  21. Schjjve, J. THE EFFECT OF PRE-STRAIN ON FATIGUE CRACK GROWTH AND CRACK CLOSURE. Engineering Fracture Mechanics 1976, 8, 575–581. [Google Scholar] [CrossRef]
  22. Alrubaie, K.; Barroso, E.K.; Godefroid, L.B. Fatigue crack growth analysis of pre-strained 7475–T7351 aluminum alloy. Int. J. Fatigue 2006, 28, 934–942. [Google Scholar] [CrossRef]
  23. Anandavijayan, S.; Mehmanparast, A.; Braithwaite, J.; Brennan, F.; Chahardehi, A. Material pre-straining effects on fatigue behaviour of S355 structural steel. J. Constr. Steel Res. 2021, 183, 106707. [Google Scholar] [CrossRef]
  24. Kang, T.S.; Liu, H. The effect of pre-stress cycles on fatigue crack growth—an analysis of crack growth mechanism. Eng. Fract. Mech. 1974, 6, 631–638. [Google Scholar] [CrossRef]
  25. Wasén, J.; Hamberg, K.; Karlsson, B. The influence of prestrain and ageing on fatigue crack growth in a dual-phase steel. Scr. Met. 1984, 18, 621–624. [Google Scholar] [CrossRef]
  26. Arora, P.; Raghavan, M.; Prasad, Y. Effect of prestrain on fatigue crack propagation in mild steel containing non-metallic inclusions. Eng. Fract. Mech. 1988, 29, 67–79. [Google Scholar] [CrossRef]
  27. Radhakrishnan, V.; Baburamani, P. An investigation of the effect of pre-straining on fatigue crack growth. Mater. Sci. Eng. 1975, 17, 283–288. [Google Scholar] [CrossRef]
  28. Kim, S.; Song, T.; Sung, H.; Kim, S. Effect of Pre-Straining on High Cycle Fatigue and Fatigue Crack Propagation Behaviors of Precipitation Hardened Steel. Met. Mater. Int. 2019, 27, 1383–1391. [Google Scholar] [CrossRef]
  29. Leitner, T.; Trummer, G.; Pippan, R.; Hohenwarter, A. Influence of severe plastic deformation and specimen orientation on the fatigue crack propagation behavior of a pearlitic steel. Mater. Sci. Eng. A 2018, 710, 260–270. [Google Scholar] [CrossRef]
  30. Ashkenazi, D. How aluminum changed the world: A metallurgical revolution through technological and cultural perspectives. Technol. Forecast. Soc. Chang. 2019, 143, 101–113. [Google Scholar] [CrossRef]
  31. Zhang, X.; Chen, Y.; Hu, J. Recent advances in the development of aerospace materials. Prog. Aerosp. Sci. 2018, 97, 22–34. [Google Scholar] [CrossRef]
  32. Dursun, T.; Soutis, C. Recent developments in advanced aircraft aluminium alloys. Mater. Des. 2014, 56, 862–871. [Google Scholar] [CrossRef]
  33. Swift, H. Plastic instability under plane stress. J. Mech. Phys. Solids 1952, 1, 1–18. [Google Scholar] [CrossRef]
  34. Lemaitre, J.; Chaboche, J.-L. Mechanics of Solid Materials; Cambridge University Press: 1990. [CrossRef]
  35. Borges, M.; Lopez-Crespo, P.; Antunes, F.; Moreno, B.; Prates, P.; Camas, D.; Neto, D. Fatigue crack propagation analysis in 2024-T351 aluminium alloy using nonlinear parameters. Int. J. Fatigue 2021, 153, 106478. [Google Scholar] [CrossRef]
  36. Menezes, L.; Teodosiu, C. Three-dimensional numerical simulation of the deep-drawing process using solid finite elements. J. Mech. Work. Technol. 2000, 97, 100–106. [Google Scholar] [CrossRef]
  37. Borges, M.; Neto, D.; Antunes, F. Numerical simulation of fatigue crack growth based on accumulated plastic strain. Theor. Appl. Fract. Mech. 2020, 108, 102676. [Google Scholar] [CrossRef]
  38. Mehrzadi, M.; Taheri, F. The influence of negative and positive stress ratios on crack growth rate in AM60B magnesium alloy. Mater. Sci. Eng. A 2012, 545, 68–77. [Google Scholar] [CrossRef]
  39. Seifi, R.; Hosseini, R. Experimental study of fatigue crack growth in raw and annealed pure copper with considering cyclic plastic effects. Theor. Appl. Fract. Mech. 2018, 94, 1–9. [Google Scholar] [CrossRef]
  40. Borrego, L.F.; Ferreira, J.; Costa, J. Fatigue crack growth and crack closure in an AlMgSi alloy. Fatigue Fract. Eng. Mater. Struct. 2001, 24, 255–265. [Google Scholar] [CrossRef]
  41. Zheng, S.; Yu, Q.; Gao, Z.; Jiang, Y. Loading history effect on fatigue crack growth of extruded AZ31B magnesium alloy. Eng. Fract. Mech. 2013, 114, 42–54. [Google Scholar] [CrossRef]
  42. Baptista, J.; Antunes, F.; Correia, L.; Branco, R. A numerical study of the effect of single overloads on plasticity induced crack closure. Theor. Appl. Fract. Mech. 2017, 88, 51–63. [Google Scholar] [CrossRef]
  43. Neto, D.; Borges, M.; Antunes, F.; Jesus, J. Mechanisms of fatigue crack growth in Ti-6Al-4V alloy subjected to single overloads. Theor. Appl. Fract. Mech. 2021, 114. [Google Scholar] [CrossRef]
  44. Neto, D.; Borges, M.; Antunes, F.; Sunder, R. Numerical analysis of fatigue crack growth under SuperBlock2020 load sequence. Eng. Fract. Mech. 2022, 260. [Google Scholar] [CrossRef]
  45. Neto, D.M.; Pedro, J.; Borges, M.F.; Borrego, L.F.P.; Sérgio, E.R.; Antunes, F.V. Numerical prediction of fatigue crack growth based on cumulative plastic strain versus experimental results for AA6082-T6. Int. J. Fract. 2022, 240, 167–181. [Google Scholar] [CrossRef]
  46. Neto, D.; Sérgio, E.; Borges, M.; Borrego, L.; Antunes, F. Effect of load blocks on fatigue crack growth. Int. J. Fatigue 2022, 162. [Google Scholar] [CrossRef]
  47. Neto, D.; Cavaleiro, N.; Sérgio, E.; Jesus, J.; Camacho-Reyes, A.; Antunes, F. Effect of crack flank holes on fatigue crack growth. Int. J. Fatigue 2023, 170. [Google Scholar] [CrossRef]
  48. Rana, M.S.; Makabe, C.; Fujiwara, G. The effect of hole shape on the extent of fatigue life improvement by cold expansions. Eng. Fail. Anal. 2009, 16, 2081–2090. [Google Scholar] [CrossRef]
  49. Chakherlou, T.; Shakouri, M.; Aghdam, A.; Akbari, A. Effect of cold expansion on the fatigue life of Al 2024-T3 in double shear lap joints: Experimental and numerical investigations. Mater. Des. 2012, 33, 185–196. [Google Scholar] [CrossRef]
  50. Semari, Z.; Aid, A.; Benhamena, A.; Amrouche, A.; Benguediab, M.; Sadok, A.; Benseddiq, N. Effect of residual stresses induced by cold expansion on the crack growth in 6082 aluminum alloy. Eng. Fract. Mech. 2013, 99, 159–168. [Google Scholar] [CrossRef]
  51. Gopalakrishna, H.; Murthy, H.N.; Krishna, M.; Vinod, M.; Suresh, A. Cold expansion of holes and resulting fatigue life enhancement and residual stresses in Al 2024 T3 alloy – An experimental study. Eng. Fail. Anal. 2010, 17, 361–368. [Google Scholar] [CrossRef]
  52. Radhakrishnan, V.M.; Baburamani, P.S. Initiation and propagation of fatigue crack in pre-strained material. International Journal of Fracture 1976. [Google Scholar] [CrossRef]
  53. Li, H.F.; Wang, Z.B.; Sun, X.J.; Qian, C.F. Experimental investigation of the strain-strengthen effect on the existed I+II mixed mode fatigue crack transformation propagation in an austenite stainless steel [C]. In: Proceed-ings of the ASME 2013 Pressure Vessels & Piping Division Conference. 2013, 7: 14-18.
  54. Moreno, B.; Zapatero, J.; Domínguez, J. An experimental analysis of fatigue crack growth under random loading. Int. J. Fatigue 2003, 25, 597–608. [Google Scholar] [CrossRef]
Figure 1. (a) Geometry of CT specimen (dimensions in mm). (b) (c) Physical model.
Figure 1. (a) Geometry of CT specimen (dimensions in mm). (b) (c) Physical model.
Preprints 147895 g001
Figure 2. Finite element mesh.
Figure 2. Finite element mesh.
Preprints 147895 g002
Figure 3. Geometry of CT specimen (dimensions in mm). (b) Physical model.
Figure 3. Geometry of CT specimen (dimensions in mm). (b) Physical model.
Preprints 147895 g003
Figure 4. Effect of pre-strain on FCG (Plane stress; R=0.1; Fmin=4.17 N; Fmax=41.67 N; with contact).
Figure 4. Effect of pre-strain on FCG (Plane stress; R=0.1; Fmin=4.17 N; Fmax=41.67 N; with contact).
Preprints 147895 g004
Figure 5. Variation of Crack tip Opening Displacement (CTOD) with applied load (a=16.184 mm; Plane stress; R=0.1; Fmin=4.17 N; Fmax=41.67 N; with contact).
Figure 5. Variation of Crack tip Opening Displacement (CTOD) with applied load (a=16.184 mm; Plane stress; R=0.1; Fmin=4.17 N; Fmax=41.67 N; with contact).
Preprints 147895 g005
Figure 6. FCG rate versus effective stress intensity factor range for three crack lengths (a=16.184; a=15.048 mm and a=15.152 mm) and three pre-strain conditions (εps=-4, 0, +4%). (Plane stress; R=0.1; Fmin=4.17 N; Fmax=41.67 N; with contact).
Figure 6. FCG rate versus effective stress intensity factor range for three crack lengths (a=16.184; a=15.048 mm and a=15.152 mm) and three pre-strain conditions (εps=-4, 0, +4%). (Plane stress; R=0.1; Fmin=4.17 N; Fmax=41.67 N; with contact).
Preprints 147895 g006
Figure 7. Effect of stress ratio on FCG for different pre-strains (Plane stress; R=0.5, Fmin=4.17 N; Fmax=41.67 N; with contact).
Figure 7. Effect of stress ratio on FCG for different pre-strains (Plane stress; R=0.5, Fmin=4.17 N; Fmax=41.67 N; with contact).
Preprints 147895 g007
Figure 8. CTOD versus load (a=16.184 mm; Plane stress; Fmin=37.5 N; Fmax=75 N; R=0.5; with contact).
Figure 8. CTOD versus load (a=16.184 mm; Plane stress; Fmin=37.5 N; Fmax=75 N; R=0.5; with contact).
Preprints 147895 g008
Figure 9. Influence of pre-strain on the effect of overload (OLR=1.5; FOL=79.17 N; Plane stress; Fmin,BL=4.17 N; Fmax,BL=41.67 N; R=0.1; with contact).
Figure 9. Influence of pre-strain on the effect of overload (OLR=1.5; FOL=79.17 N; Plane stress; Fmin,BL=4.17 N; Fmax,BL=41.67 N; R=0.1; with contact).
Preprints 147895 g009
Figure 10. CTOD versus applied load (OLR=1.5; FOL=79.17 N; Plane stress; Fmin,BL=4.17 N; Fmax,BL=41.67 N; R=0.1; with contact). (a) Overload cycle. (b) At minimum da/dN (a=16.336 mm and a=15.448 mm for εps=0 and 4%, respectively.
Figure 10. CTOD versus applied load (OLR=1.5; FOL=79.17 N; Plane stress; Fmin,BL=4.17 N; Fmax,BL=41.67 N; R=0.1; with contact). (a) Overload cycle. (b) At minimum da/dN (a=16.336 mm and a=15.448 mm for εps=0 and 4%, respectively.
Preprints 147895 g010aPreprints 147895 g010b
Figure 11. (a) da/dN versus crack length. (b) CTOD versus load plots at a=16.048 mm (Plane stress; Fmin=4.17 N; Fmax=41.67 N; R=0.1; without contact).
Figure 11. (a) da/dN versus crack length. (b) CTOD versus load plots at a=16.048 mm (Plane stress; Fmin=4.17 N; Fmax=41.67 N; R=0.1; without contact).
Preprints 147895 g011
Figure 12. Accumulation of plastic strain at the crack tip with cyclic loading. (a) With contact of crack flanks. (b) Without contact of crack flanks (Plane stress; Fmin=4.17 N; Fmax=41.67 N; R=0.1).
Figure 12. Accumulation of plastic strain at the crack tip with cyclic loading. (a) With contact of crack flanks. (b) Without contact of crack flanks (Plane stress; Fmin=4.17 N; Fmax=41.67 N; R=0.1).
Preprints 147895 g012
Table 1. Chemical composition in weight % of 2024-T351 aluminium alloy. The balance is Al.
Table 1. Chemical composition in weight % of 2024-T351 aluminium alloy. The balance is Al.
Si Fe Cu Mn Mg Cr Zn Ti
0.50 0.50 3.8-4.9 0.3-0.9 1.2-1.8 0.10 0.25 0.15
Table 2. Mechanical properties of 2024-T351 aluminium alloy without pre-strain [24].
Table 2. Mechanical properties of 2024-T351 aluminium alloy without pre-strain [24].
Young modulus Yield stress UTS Elongation at break Brinell Hardness
73 GPa 325 MPa 470 MPa 20% 137
Table 3. List of material parameters involved in the Swift and Lemaître-Chaboche laws.
Table 3. List of material parameters involved in the Swift and Lemaître-Chaboche laws.
Material Y0 [MPa] C [N] n Cx [mm] XSat [MPa]
2024-T351 288.96 389.00 0.056 138.80 111.84
Table 4. Loading conditions (a0=15 mm; plane stress).
Table 4. Loading conditions (a0=15 mm; plane stress).
Test Loading R εps [%] Contact
1 CA 0.1 0 Yes
2 CA 0.1 +4 Yes
3 CA 0.1 -4 Yes
4 CA 0.1 0 No
5 CA 0.1 +4 No
6 CA 0.1 -4 No
7 CA 0.5 0 Yes
8 CA 0.5 +4 Yes
9 CA 0.5 -4 Yes
10 OL 0.1 (BL) 1 0 Yes
11 OL 0.1 (BL) +4 Yes
1 BL- Baseline loading.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2025 MDPI (Basel, Switzerland) unless otherwise stated