Numerical simulation of fatigue crack growth

A numerical simulation of fatigue crack growth which uses currently available crack tip stress and strain fields is described. The essential features of the numerical model are the concepts of damage accumulation cycle by cycle and repeated re-initiation at the tip of the growing crack. The failure criteria employed are a combination of a failure condition and a critical distance over which this condition must be achieved. This critical distance, the material size parameter, has a magnitude which depends on the failure mechanism.

The use of the model to illustrate the effects of stress ratio and environmental effects is described and the ability of the model to predict the onset of bursts of crack growth due to static failure mechanisms is demonstrated. The phenomenon of self-arresting cracks is also displayed.

Material characteristics are included in the model and comparisons with experimental data are presented for a C-Mn steel used in the fabrication of offshore structures.     

Stage II fatigue crack growth

The linear part of the fatigue crack growth diagram is found to be divided into Stages IIa and IIb by the point O whose coordinates K^* and A are dependent on the physical and structural characteristics of the material. In Stage IIa K_eff remains constant as the microcrack advances in increments corresponding to the dislocation cell structure size, λ, pausing for (dN−1) cycles to accumulate the elastic energy required for the crack opening. During Stage IIb K_op remains constant and the microcrack opens during each cycle and advances irrespective of the substructure but in accordance with an increasing value of K_eff. The effects of temperature and vacuum on K^* are considered; the A values correspond to those of λ and are independent of the above effects.     

Time-derivative equations for fatigue crack growth in metals

Predicting fatigue crack growth in metals remains a difficult task because available models are based on cycle-derivative equations, such as the Paris law, while service loads are often far from being cyclic. The main objective of this paper is therefore to propose a set of time-derivative equations for fatigue crack growth. The model is based on the thermodynamics of dissipative processes. For this purpose, three global state variables are introduced in order to characterize the state of the crackthe crack length a, the plastic blunting at crack tip and the intensity of crack opening C. Thermodynamics counterparts are introduced for each variable. Special attention is paid to the elastic energy stored inside the crack tip plastic zone, because, in practice, residual stresses at crack tip are known to considerably influence fatigue crack growth. The stored energy is included in the energy balance equation, and this leads to the appearance of a kinematics hardening term in the yield criterion for the cracked structure. No dissipation is associated with crack opening, but to crack growth and to crack tip blunting. Finally, the model consists in two laws: a crack propagation law, which is a relationship between d dt and da/dt and which observes the inequality stemmed from the second principle, and an elastic-plastic constitutive behaviour for the cracked structure, which provides d dt versus applied-load. The model was implemented and tested. It reproduces successfully the main features of fatigue crack growth as reported in the literature, such as the Paris law, the stress-ratio effect and the overload retardation effect.     

Micromechanisms of fatigue crack growth in aluminium alloys

The fatigue crack growth characteristics of high-strength aluminium alloys are discussed in terms of behaviour during mechanical testing and fracture surface appearance. For a wide range of crack growth rates, the crack extends both by the formation of ductile striations and by the coalescence of micro-voids. Dimples are observed at stress intensities very much less than the plane strain fracture toughness, and this is explained in terms of the probability of inclusions lying close to the crack tip. The striation formation process is described as a combination of environmentally-enhanced cleavage processes and plastic blunting of the crack tip.     

Mixed mode fatigue crack growth: A literature survey

The applications of fracture mechanics have traditionally concentrated on crack growth problems under an opening or mode I mechanism. However, many service failures occur from growth of cracks subjected to mixed mode loadings. This paper reviews the various criteria and parameters proposed in the literature for predictions of mixed mode crack growth directions and rates. The physical basis and limitations for each criterion are briefly reviewed, and the corresponding experimental supports are discussed. Results from experimental studies using different specimen geometries and loading conditions are presented and discussed. The loading conditions discussed consist of crack growth under mode II, mode III, mixed mode I and II, and mixed mode I and III loads. The effects of important variables such as load magnitudes, material strength, initial crack tip condition, mean stress, load non-proportionality, overloads and crack closure on mixed mode crack growth directions and/or rates are also discussed.     

A cumulative model of fatigue crack growth and the crack closure effect

A model of fatigue crack growth based on an analysis of elastic/plastic stress and strain at the crack tip is presented. It is shown that the fatigue crack growth rate can be calculated using the local stress/strain at the crack tip by assuming that a small highly strained area $\text{x}{}^1$, existing at the crack tip, is responsible for the fatigue crack growth, and that the fatigue crack growth may be regarded as the cumulation of successive crack re-initiations over a distance $\text{x}{}^1$. It is shown that crack closure can be modelled using the effective contact zone g behind the crack tip. The model allows the fatigue crack growth rate over the near threshold and linear ranges of the general da/dN versus ΔK curve to be calculated. The fatigue crack growth retardation due to overload and fatigue crack arrest can also be analysed in terms of g and $\text{x}{}^1$.Calculated fatigue crack growth rates are compared with experimental ones for low and high strength steel.     

Examination of fatigue crack driving force parameter

Most of the previous parameters that utilized as a crack driving force were established in modifying the parameter K_op in Elber's effective SIF range ΔK_eff(=K_max?K_op). However, the parameters that replaced the traditional parameter K_op were based on different measurements or theoretical calculations, so it is difficult to distinguish their differences. This paper focuses on the physical meaning of compliance changes caused by plastic deformation at the crack tip; the tests were carried out under different amplitude loading for structural steel. Based on these test results, differences of several parameter ΔK_eff in literature are analysed and an improved two‐parameter driving force ΔK_drive(=(K_max)ⁿ(ΔK^∧)^1‐n) has been proposed. Experimental data for several different types of materials taken from literature were used in the analyses. Presented results indicate that the ΔK_drive parameter was equally effective or better than ΔK(=K_max?K_min), ΔK_eff(=K_max?K_op) and ΔK*(= (K_max)^α(ΔK⁺)^1?α) in correlating and predicting the R‐ratio effects on fatigue crack growth rate.     

Bifurcation analysis of the Hobson short fatigue crack growth law

A bifurcation analysis of the Hobson short fatigue crack growth law is presented. The analysis reveals that, although the growth law is non‐linear, it contains no bifurcation points.     

Time derivative equations for mode I fatigue crack growth in metals

Predicting fatigue crack growth in metals remains a difficult task since the available models based on the Paris law are cycle-derivative equations (da/dN), while service loads are often far from being cyclic. This imposes a cycle-reconstruction of the load sequence, which significantly modifies the load history in the signal. The main objective of this paper is therefore to propose a set of time-derivative equations for fatigue crack growth in order to avoid any cycle reconstruction. The model is based on the thermodynamics of dissipative processes. Its main originality lies in the introduction of a supplementary state variable for the crack, which allows describing continuously the state of the crack throughout any complex load sequence. The state of the crack is considered to be fully characterized at the global scale by its length a, its plastic blunting ρ, and its elastic opening. In the equations, special attention is paid to the elastic energy stored inside the crack tip plastic zone, since, in practice, residual stresses at the crack tip are known to considerably influence fatigue crack growth. The model consists finally in two laws: a crack propagation law, which is a relationship between dρ/dt and da/dt and which observes the inequality stemming from the inequality of Clausius Duhem, and an elastic–plastic constitutive behaviour for the cracked structure, which provides dρ/dt versus load and which stems from the energy balance equation. The model was implemented and tested. It successfully reproduces the main features of fatigue crack growth as reported in the literature, such as the Paris law, the stress ratio effect, and the overload retardation effect.     

The influence of cross-sectional thickness on fatigue crack growth

For thin structures, fatigue crack growth rates may vary with the structure's thickness for a given stress intensity factor range. This effect is mainly due to the change in the nature of the plastic deformation when the plastic zone size becomes comparable with, or greater than, the cross-sectional thickness. Variations in the constraint affect both the crack tip plastic blunting behaviour as well as the fatigue crack closure level. Approximate expressions are constructed for the constraint factor based on asymptotic values and numerical results, which are shown to correlate well with finite element results. It is demonstrated that the present results not only permit predictions of the specimen thickness effects on fatigue crack propagation under spectrum loading, but also eliminate the need to determine the constraint factor by curve-fitting crack growth data.     

Interfacial fatigue crack growth in foam core sandwich structures

This paper deals with the experimental measurement of face/core interfacial fatigue crack growth rates in foam core sandwich beams. The so-called ‘cracked sandwich beam’ specimen is used, slightly modified, which is a sandwich beam that has a simulated face/core interface crack. The specimen is precracked so that a more realistic crack front is created prior to fatigue growth measurements. The crack is then propagated along the interface, in the core material, during fatigue loading, as is assumed to occur in a real sandwich structure. The crack growth is stable even under constant amplitude testing. Stress intensity factors are obtained from the FEM which, combined with the experimental data, result in standard da/dN versus ΔK curves for which classical Paris’ law constants can be extracted. The experiments to determine stress intensity factor threshold values are performed using a manual load-shedding technique.     

A cumulative model of fatigue crack growth

A model of fatigue crack growth based on an analysis of elastic/plastic stress and strain at the crack tip is presented. It is shown that the fatigue crack growth rate can be calculated by means of the local stress/strain at the crack tip. The local stress and strain calculations are based on the general solutions given by Hutchinson, Rice and Rosengren. It is assumed that a small highly strained area existing at the crack tip is responsible for the fatigue crack growth. It is also assumed that the fatigue crack growth rate depends mainly on the width, x¹, of the highly strained zone and on the strain range, $\text{Δ?}{}^1$, within the zone. A relationship between stress intensity factor K and the local strain and stress has been developed. It is possible to calculate the local strain for a variety of crack problems. Then, the number of cycles N¹ required for material failure inside the highly strained zone is calculated. The fatigue crack growth rate is calculated as the ratio $\text{x}{}^1\text{N}{}^1$.The calculated fatigue crack growth rates were compared to the experimental ones. Two alloys steels and two aluminium alloys were analyzed. Good agreement between experimental and theoretical results is obtained.     

A fatigue crack growth model for interacting cracks in a plate of arbitrary thickness

Fatigue and fracture assessment of structures weakened by multiple site damage, such as two or more interacting cracks, represents a very challenging problem. A proper analysis of this problem often requires advanced modelling approaches. The objective of this paper is to develop a general theoretical approach and investigate the fatigue behaviour of two interacting cracks. The developed approach is based on the classical strip yield model and plasticity induced crack closure concept. It also utilises the 3D fundamental solution for an edge dislocation. The crack advance scheme adopts the cycle‐by‐cycle calculations of the effective stress intensity factors and crack increments. The modelling results were validated against experimental data available in the literature. Further, the nonlinear effects of the crack interaction and plate thickness on the crack opening stresses and crack growth rates were studied with the new approach for the problem geometry. It was demonstrated that the both effects could have a significant influence on fatigue life and cannot be disregarded in life and integrity assessments of structural components with multiple site damage.     

Mixed-mode fatigue crack growth under biaxial loading

Studies on crack growth in a panel with an inclined crack subjected to biaxial tensile fatigue loading are presented. The strain energy density factor approach is used to characterize the fatigue crack growth. The crack growth trajectory as a function of the initial crack angle and the biaxiality ratio is also predicted. The analysis is applied to 7075-T6 aluminium alloy to predict the dependence of crack growth rate on the crack angle. The effect of crack angle on the cyclic life of the component and on the cyclic life ratio is presented and discussed.     

Effects of spectrum modification on fatigue crack growth

The purpose of this study was to investigate experimentally the effects associated with modification of a loading spectrum recorded from P‐3 a maritime aircraft on fatigue crack growth behaviour. The material is 2324‐T39 Al alloy widely used in the aircraft industry. Experiments were conducted using the full spectrum and modified versions of it such as only ‘positive’ (no negative loads) or with reduced (clipped) high peaks. The results show that the compressive loads decrease fatigue life of the specimen by ∼300%. Furthermore, by running tests with clipped peaks it was found that the fatigue life was shorten significantly due to reduction of crack growth retardation caused by highest tensile peaks. Multiple tests were conducted in order to establish a scatter in the experimental data under spectrum loads.     

Effect on fatigue crack growth of interactions between overloads

ABSTRACT Various types of interactions between overloads were studied in a 0.38% C low carbon steel. The retarding effect due to consecutive overloads is found to increase with the number of overloads, until it reaches a maximum. Similarly, it is found that a critical distance between overloads ensures the highest retarding effect, while shorter or longer spacing are less efficient for retarding crack growth. These effects are successfully explained using FEM calculations of the effective stress intensity factor. The kinematic hardening of the alloy, which is very efficient in ferritic–pearlitic steels, is shown to be mostly responsible for those effects. Taking into account the amplitude of kinematic hardening allows qualitative explanation of the observed effects. The order of application of the cycles during variable amplitude fatigue is thus important and should be taken into account for predicting fatigue lives.     

Prediction of fatigue crack growth based on low cycle fatigue properties

Theoretical models of the fatigue crack growth without artificial adjustable parameters were proposed by considering the plastic strain energy and the linear damage accumulation, respectively. The crack was regarded as a sharp notch with a small curvature radius and the process zone was assumed to be the size of cyclic plastic zone. The near crack tip elastic–plastic stress and strain were evaluated in terms of modified Hutchinson, Rice and Rosengren (HRR) formulations. Predicted results from two established models have been soundly compared with open reports for frequently used materials. It is found that experimental results agree well with theoretical solutions.     

考虑结构裂纹扩展的振动疲劳寿命计算方法

工程实践中任何结构都存在不同程度的裂纹损伤,振动激励下动响应与疲劳裂纹扩展之间互相耦合,直接影响结构振动疲劳寿命.为了考虑结构振动疲劳耦合效应对疲劳寿命的影响,提出了一种考虑结构裂纹扩展的振动疲劳寿命计算方法.分析时,通过建立若干个含不同长度裂纹的结构有限元模型模拟结构裂纹扩展,采用Paris方程分段计算结构振动疲劳裂纹扩展寿命,通过试验确定的固有频率降变化规律反推结构裂纹萌生寿命,最后累计得到结构疲劳总寿命.结论表明,仿真计算结果与试验结果比较吻合.     

Prediction of threshold and ultra-low fatigue crack growth rates

A model was derived to predict the true threshold value for fatigue crack growth in the absence of crack closure. The model, based only on the tensile and cyclic properties of the material, was successfully verified against a set of experimental data on medium and high strength steels and one aluminium alloy. Good agreement with experimental results was also obtained for Region I of the $\text{d}\text{a/}\text{d}\text{N vs ΔK}$ curve using a fatigue crack growth rate equation based on the same model.Fatigue crack growth data obtained from the medium strength steel CK45 in the normalized state and two heat-treated conditions were analysed. Good data correlation was shown using a previously developed normalizing parameter, $\text{φ = (ΔK}{}^2\text{?ΔK}{}^2{}_{\mathrm{<mtext>th</mtext>}}\text{)/(K}{}^2{}_{\mathrm{<mtext>c</mtext>}}\text{?K}{}^2{}_{\mathrm{<mtext>max</mtext>}}\text{)}$, in the entire range of fatigue crack growth rates and for stress ratios ranging from 0.1 to 0.8.