Local weight functions for semi-elliptical surface cracks in finite thickness plates

Weight functions for any local point, 0 < Φ < π/2 along a semi-elliptical surface crack in finite thickness plates were derived from an assumed approximate general weight function and two reference stress intensity factors. The resulting weight functions were verified using available finite element results for two nonlinear stress fields and good agreement was achieved. When used together with weight functions for Φ = 0 and Φ = π/2 the weight functions are suitable for the calculation of stress intensity factors anywhere along the crack front for semi-elliptical surface cracks in complex stress fields with aspect ratios in the range 0 ≤ a/c ≤ 1 and relative depths 0 ≤ a/t ≤ 0.8.     

Effect of stress ratio and specimen thickness on fatigue crack growth of CK45 steel

Presented are the effect of stress ratio and thickness on the fatigue crack growth rate of CK45 steel according to DIN 17200. Test results are obtained for constant amplitude load in tension with three stress ratios of R=0, 0.2 and 0.4 and three specimen thicknesses of B=6, 12 and 24 mm. Microgauge crack opening values were used to calculate ΔK_eff values from which the da/dN − ΔK_eff curves are obtained. Crack closure can be applied to explain the influence of mean stress and specimen thickness on the fatigue crack growth rate in the second regime of the two-parameter crack growth rate relation. An empirical model is chosen for calculating the normalized load ratio parameter U as a function of R, B and ΔK and, for correlating the test data.     

The effect of contact stress intensity factors on fatigue crack propagation

Using the technique of Dimensional Analysis the phenomenon of crack closure is modelled using the concept of a contact stress intensity factor Kc. For constant amplitude loading, a simple expression, Kc_max = g(R) ΔK, is obtained without making idealized assumptions concerning crack tip behaviour. Further, by assuming that crack closure arises from the interaction of residual plasticity in the wake of the crack and crack tip compressive stresses, the function g(R) is shown to be constant for non-workhardening materials. This implies that any dependency of Kc_max on R must be attributed to the workhardening characteristic of the material. With Kc known, an “effective” stress intensity factor Ke may be calculated and incorporated into a crack growth law of the form da/dn = f(ΔKe). From analysis, it can be deduced that for a workhardening material, Kc_max will decrease as R increases and the effective stress intensity factor will increase. This means that the fatigue crack propagation rate will increase with R, in accordance with experimental observations.     

Fatigue crack growth rate in ferrite-martensite dual-phase steel

An empirical study is made on the fatigue crack growth rate in ferrite-martensite dual-phase (FMDP) steel. Particular attention is given to the effect of ferrite content in the range of 24.2% to 41.5% where good fatigue resistance was found at 33.8%. Variations in ferrite content did not affect the crack growth rate da/dN when plotted against the effective stress intensity factor range ΔK_eff which was assumed to follow a linear relation with the crack tip stress intensity factor range ΔK. A high ΔK_eff corresponds to uniformly distributed small size ferrite and martensite. No other appreciable correlation could be ralated to the microstructure morphology of the FMDP steel. The closure stress intensity factor K_cl, however, is affected by the ferrite content with K_cl/Kmax reaching a maximum value of 0.7. In general, crack growth followed the interphase between the martensite and ferrite.Dividing the fatigue crack growth process into Stage I and II where the former would be highly sensitive to changes in ΔK and the latter would increase with ΔK depending on the R = σ_min/σ_max ratio. The same data when correlated with the strain energy density factor range ΔS showed negligible dependence on mean stress or R ratio for Stage I crack growth. A parameter α involving the ratio of ultimate stress to yield stress, percent reduction of area and R is introduced for Stage II crack growth so that the da/dN data for different R would collapse onto a single curve with a narrow scatter band when plotted against αΔS.     

A centrally cracked thin circular disk, Part I: 3D Elastic–plastic finite element analysis

A full field solution, based on small deformation, three-dimensional elastic–plastic finite element analysis of the centrally cracked thin disk under mode I loading has been performed. The solution for the stresses under small-scale yielding and lo!cally fully plastic state has been compared with the HRR plane stress solution. At the outside of the 3D zone, within a distance of rσ_o/J=18, HRR dominance is maintained in the presence of a significant amount of compressive stress along the crack flanks. Ahead of this region, the HRR field overestimate the stresses. These results demonstrate a completely reversed state of stress in the near crack front compared to that in the plane strain case. The combined effect of geometry and finite thickness of the specimen on elastic–plastic crack tip stress field has been explored. To the best of our knowledge, such an attempt in the published literature has not been made yet. For the qualitative assessment of the results some of the field parameters have been compared to the available experimental results of K, gives a fair estimate of the crack opening stress near the crack front at a distance of order 10⁻² in. On the basis of this analysis, the Linear Elastic Fracture Mechanics approach has been adopted in analyzing the fatigue crack extension experiments performed in the disk (Part II).     

Rail web fracture in the presence of residual stresses

A recent accident involving roller-straighened alloy rail has raised the question of the safety of such rails. This work shows that the residual stresses in roller-straightened rail can indeed self-drive a long web crack. Specifically, the stress intensity K_I due to release of the key component, longitudinal stress, if of the order of the critical stress intensity for initiation K_Ic for both plain carbon and alloy rail. At cut ends, the resulting vertical residual stresses can give rise to K_Ic if there are 0.1–1 in (3–25 mm) cracks. In this work, checks of the existing residual stress data for self-consistency suggest that the data are only accurate within a factor of two. Therefore, a more direct method is proposed for measuring K_I on a web crack by saw-cutting the web.     

Residual stresses near a rail end

Stress fields near a cut end of a rail containing longitudinal residual stress typical of roller-straightened rail were studied using analysis and a finite element model. For a self-equilibrating residual stress distribution with equal maximum and minimum stresses, the distance to reach 95% of the mid-rail residual stress field is from 0.7 to 1.8 times the rail height, with the finite element model predicting a length of 1.1 times the rail height. This gives a measure of the accuracy of the simpler analytical models. At the rail end, the longitudinal residual stress goes to zero, and the vertical residual stress near mid-web reaches a maximum of approximately 27 ksi (186 MPa) (1.35 times the maximum mid-rail longitudinal residual stress of 20 ksi, or 138 MPa). The maximum shear stresses are 6 ksi (41 MPa) and −8 ksi (−55 MPa) near the head-web and web-base intersections, respectively, approximately 2 in. (51 mm) from the end of a 7.3 in. (185 mm) high rail. The shear stress is zero at the cut end and in mid-rail. The worst possible end-crack is a horizontal web crack in the vertical residual stress field at the rail end. The stress intensity K_I on such a crack is estimated to reach 20 ksi√in. (22 MPa√m) for cracks 0.5 in. (13mm) long. This is already 0.4 to 0.8 times K_I for carbon and alloy rails, and about 0.5 times K_Ic for a long crack.     

Influence of mean stress on the fatigue crack growth characreristics of CrlMo steel tube at elevated temperature

The fatigue crack growth characteristics of CrlMo steel have been investigated at 861 K over the R-ratio range 0.1–0.7 utilising a dwell time of 10 min. at maximum load. All tests were conducted under load control in a laboratory air environment. It was established that the R-ratio significantly affected the fatigue crack extension behaviour inasmuch that with increasing R-ratio, the critical ΔK level for the onset of creep fatigue interactive growth, ΔK^IG, decreased from 20 to 7 MPa√m and the threshold stress intensity, ΔK_th, decreased from 9 to about 3 MPa√m. At intermediate ΔK levels, i.e. between ΔK_th and ΔK^IG, the fatigue crack extension rates, for all R-ratio values, resided on or slightly below the CTOD line, which represents the upper bound for contrnuum controlled fatigue crack growth. Creep fatigue interactive growth was typified by crack extension rates that reside above the CTOD line with a ΔK^IG dependence; the attainment of some critical creep condition or crack linkage condition which causes the abrupt change in crack extension behaviour at ΔK^IG; and crack extension occurs almost exclusively in an intergranular manner. The R-ratio and ΔK^IG followed a linear relation. A literature review concerning the effect of temperature on the threshold fatigue crack growth characteristics of low alloy ferritic steels demonstrated powerful effects of temperature; the magnitude of these effects, however, were dependent upon the testing temperature regime and R-ratio level. The effect of R-ratio on ΔK_th was greatest at temperatures >400°C, significant at ambient temperatures and least in the temperature range 90°C to <300°C. The relationship between temperature and ΔK_th, at a given R-ratio, exhibited a through and a minimum ΔK_th value was observed in the temperature range 200–250°C. The magnitude of the temperature effects on ΔK_th decreased with increasing R-ratio. Such effects of temperature and R-ratio on ΔK_th was reasonably explained in terms of crack closure effects. Finally, the present elevated temperature fatigue crack growth data exhibited massive crack extension enhancement values when compared to ambient near-threshold fatigue crack growth data for CrlMo steel. Such large enhancement values were the combined effects of temperature (environment) and frequency.     

Kinetics of microcrack blunting ahead of macrocrack approaching shear wave speed

The fracturing of glass and tearing of rubber both involve the separation of material but their crack growth behavior can be quite different, particularly with reference to the distance of separation of the adjacent planes of material and the speed at which they separate. Relatively speaking, the former and the latter are recognized, respectively, to be fast and slow under normal conditions. Moreover, the crack tip radius of curvature in glass can be very sharp while that in the rubber can be very blunt. These changes in the geometric features of the crack or defect, however, have not been incorporated into the modeling of running cracks because the mathematical treatment makes use of the Galilean transformation where the crack opening distance or the change in the radius of curvature of the crack does not enter into the solution. Change in crack speed is accounted for only via the modulus of elasticity and mass density. For this simple reason, many of the dynamic features of the running crack have remained unexplained although speculations are not lacking. To begin with, the process of energy dissipation due to separation is affected by the microstructure of the material that distinguishes polycrystalline from amorphous form. Energy extracted from macroscopic reaches of a solid will travel to the atomic or smaller regions at different speeds at a given instance. It is not clear how many of the succeeding size scales should be included within a given time interval for an accurate prediction of the macroscopic dynamic crack characteristics. The minimum requirement would therefore necessitate the simultaneous treatment of two scales at the same time. This means that the analysis should capture the change in the macroscopic and microscopic features of a defect as it propagates. The discussion for a dual scale model has been invoked only very recently for a stationary crack. The objective of this work is to extend this effort to a crack running at constant speed beyond that of Rayleigh wave. Developed is a dual scale moving crack model containing microscopic damage ahead of a macroscopic crack with a gradual transition. This transitory region is referred to as the mesoscopic zone where the tractions prevail on the damaged portion of the material ahead of the original crack known as the restraining stresses, the magnitude of which depends on the geometry, material and loading. This damaged or restraining zone is not assumed arbitrarily nor assumed to be intrinsically a constant in the cohesive stress approach; it is determined for each step of crack advancement. For the range of micronotch bluntness with 0 < β < 30° and 0.2 σ_∞/σ₀ 0.5, there prevails a nearly constant restraining zone size as the crack approaches the shear wave speed. Note that β is the half micronotch angle and the applied stress ratio is σ_∞/σ₀ with σ₀ being the maximum of the restraining stress. For σ_∞/σ₀ equal to or less than 0.5, the macrocrack opening displacement COD is nearly constant and starts to decrease more quickly as the crack approaches the shear wave speed. For the present dual scale model where the normalized crack speed v/c_s increases with decreasing with the one-half microcrack tip angle β. There prevails a limit of crack tip bluntness that corresponds to β 36° and v/c_s 0.15. That is a crack cannot be maintained at a constant speed if the bluntness is increased beyond this limiting value. Such a feature is manifestation of the dependency of the restraining stress on crack velocity and the applied stress or the energy pumped into the system to maintain the crack at a constant velocity. More specifically, the transitory character from macro to micro is being determined as part of the unknown solution. Using the energy density function dW/dV as the indicator, plots are made in terms of the macrodistance ahead of the original crack while the microdefect bluntness can vary depending on the tip geometry. Such a generality has not been considered previously. The macro-dW/dV behavior with distance remains as the inverse r relation yielding a perfect hyperbola for the homogeneous material. This behavior is the same as the stationary crack. The micro-dW/dV relations are expressed in terms of a single undetermined parameter. Its evaluation is beyond the scope of this investigation although the qualitative behavior is expected to be similar to that for the stationary crack. To reiterate, what has been achieved as an objective is a model that accounts for the thickness of a running crack since the surface of separation representing damage at the macroscopic and microscopic scale is different. The transitory behavior from micro to macro is described by the state of affairs in the mesoscopic zone.     

Non-self-similar crack growth in elastic-plastic finite thickness plate

This work is concerned with non-self-similar crack growth in medium strength metal plates while the loading step, plate thickness and material properties are altered. The three-dimensional elastic-plastic finite element stress analysis is combined with the strain energy density criterion for modeling the material damage process from crack initiation to final global instability including the intervening stage of slow crack growth. Both inelastic deformation and crack growth are accounted for each increment of loading such that the redistribution of stresses and strains are made for each new crack profile. Numerical results are obtained for the center cracked plate configuration under uniform extension with twenty-seven (27) different combinations of specimen thickness, loading step and material type. The fracture toughness S_c being related to K_1c for three different materials are predicted analytically from the corresponding uniaxial tensile test data. Effective strain energy density factor and half crack length are defined so that the results can be compared with their two-dimensional counterparts. Crack growth resistance curves (R-curves) are constructed by plotting as a function of . The condition is found to prevail during slow crack growth. Translation and/or rotation of the lines can yield results other than those calculated and serve a useful purpose for scaling component size and test time. The minimum thickness requirement for the ASTM valid K_1c test is also discussed in connection with predictions based on the strain energy density criterion. The corresponding K_1c for smaller specimens that exhibit moderate ductility and nonlinearity can also be obtained analytically. In such cases, the influence of loading step can be significant and should not be neglected. Notwithstanding the shortcomings of the theory of plasticity, the qualitative features of non-self-similar crack growth are predicted by the strain energy density criterion. Any refinements on the analytical modeling of the material damage process would only affect the results qualitatively, a subject that is left for future investigation.     

Stress intensity factor range and propagation mode of surface cracks under rolling–sliding contact

Finite element analyses were conducted in order to evaluate the mode I and mode II stress intensity factors for inclined edge cracks under cyclic contact load under rolling and rolling–sliding condition. The SIF range depends on crack orientation, crack length to Hertzian contact zone half-width ratio, friction between the crack faces and friction on the contact surface. The results were combined in two compact functions that determine the ΔK_I and ΔK_II values. The crack propagation mode and direction were investigated using both the maximum stress criterion and the minimum strain energy density criterion. The results are displayed in graph form, which allows a fast evaluation of the crack growth condition.     

A two-strain-gage technique for determining mode I stress-intensity factor

A finite element analysis and a experimental test were performed to show that the terms with r^−1/2 and ^1/2 in the eigenfunction expansion of the strains can describe the crack tip strain distribution with sufficient accuracy. A set of two linear equations can be obtained to determine the stress intensity factor K_I using only two strain-gages. Errors within 5% can be achieved provided that the two strain-gages are placed at the appropriate locations. The technique can be developed to treat crack bodies with irregular geometry and complex loading.     

Numerical simulation of stress intensity factors via the thermoelastic technique

The purpose of this study is to investigate the accuracy of the least squares method for finding the in-plane stress intensity factorsK _I andK _II using thermoelastic data from isotropic materials. To fully understand the idealized condition ofK _I andK _II calculated from thermoelastic experiments, the total stress field calculated from finite element analysis is used to take the place of data obtained from real thermoelastic experiments. In the finite element analysis, theJ-integral is also calculated to compare with (K _I ² +K _II ² )/E evaluated by the least squares method. The stress fields near the crack tip are dominated by the two stress intensity factors; however, the edge effect will cause inaccuracy of the thermoelastic data near the crack tip. Furthermore, the scan area of thermoelastic experiments cannot be too small. Therefore, we suggest that three or four terms of stress function be included in the least squares method for evaluating stress intensity factors via the thermoelastic technique. In the idealized condition, the error can be smaller than 3 percent from our numerical simulations. If only ther ^–1/2 term (K _I andK _II ) is included in the least squares method, even in the idealized case the error can be up to 20 percent.     

Resolved shear stress intensity coefficient and fatigue crack growth in large crystals

Fatigue crack growth is caused primarily by shear decohesion due to dislocation motion in the crack tip region. The resolved shear stress, which drives dislocation in a crystal, is strongly orientation dependent, and therefore, the cyclic plastic deformation of the shear decohesion process is highly anisotropic.The crack planes are often inclined to the loading axis both in the inplane orientation and in the thickness direction. This inclination induces all three modes of the crack tip stress field, K_I, K_II, and K_III.Fatigue crack growth in large-grain Al 7029 aluminum alloy was studied. The crack tip stress fields of the test specimens are calculated with the finite element method. The values of K_I, K_II, and K_III are evaluated. The orientation of the crystal at a crack tip was determined with the Laue X-ray method. The crystal orientation and the calculated crack tip stress fields are used to compute the resolved shear stress intensity of each of the twelve slip systems of the crystal at the crack tip. The resolved shear stress field of a slip system is linearly proportional to the resolved shear stress intensity coefficient, RSSIC.The values of RSSIC thus evaluated are used to analyze the orientations of the crack plane and to correlate with the shear fatigue crack growth rate.     

Stress intensification in cracked shank of tightened bolt

Defects or cracks in the shank of bolts can degrade their load carrying capacity. The ways with which loading and residual stress intensify the crack border stress field can be reflected through the stress intensity factor quantity as defined in the linear elastic fracture mechanics theory. Use is made of the stiffness derivative method where quarter-point singular finite elements are used in the numerical calculation. Improved accuracy is achieved by considering the displacements not only of the main nodes but also of those quarter-point nodes in plane normal and adjacent to the crack.Numerical results are obtained for a semi-elliptical shaped crack in the bolt shank owing to tension, bending, residual stress and stress caused by tightening of the bolt. Maximum value of the Mode I stress intensity factor K_a due to tension or bending could prevail either at the deepest point on the crack border or at the root of the shank where the crack border terminates depending on the aspect ratio of the ellipse. In general, K₁ at the deepest point of crack penetration is larger than that at the free surface for tension and bending for a fixed crack depth with reference to the bolt diameter. Tightening of the bolt tends to increase K₁ at the free surface if the crack depth is small. The opposite is obtained for deeper cracks. Assumed residual stress effect obtained from experimental data is found to have negligible influence on the stress intensity factor when compared with that arising from tensile load.     

Edge fracture in cone-plate and parallel plate flows

Edge fracture is an instability of cone-plate and parallel plate flows of viscoelastic liquids and suspensions, characterised by the formation of a `crack' or indentation at a critical shear rate on the free surface of the liquid. A study is undertaken of the theoretical, experimental and computational aspects of edge fracture. The Tanner-Keentok theory of edge fracture in second-order liquids is re-examined and is approximately extended to cover the Criminale-Ericksen-Filbey (CEF) model. The second-order theory shows that the stress distribution on the semi-circular crack is not constant, requiring an average to be taken of the stress; this affects the proportionality constant, K in the edge fracture equation −N _2c = KΓ/a, where N _2c is the critical second normal stress difference, Γ is the surface tension coefficient and a is the fracture diameter. When the minimum stress is used, K = 2/3 as found by Tanner and Keentok (1983). Consideration is given to the sources of experimental error, including secondary flow and slip (wall effect). The effect of inertia on edge fracture is derived. A video camera was used to record the inception and development of edge fracture in four viscoelastic liquids and two suspensions. The recorded image was then measured to obtain the fracture diameter. The edge fracture phenomenon was examined to find its dependence on the physical dimensions of the flow (i.e. parallel plate gap or cone angle), on the surface tension coefficient, on the critical shear rate and on the critical second normal stress difference. The critical second normal stress difference was found to depend on the surface tension coefficient and the fracture diameter, as shown by the theory of Tanner and Keentok (1983); however, the experimental data were best fitted by the equation −N _2c = 1.095Γ/a. It was found that edge fracture in viscoelastic liquids depends on the Reynolds number, which is in good agreement with the inertial theory of edge fracture. Edge fracture in lubricating grease and toothpaste is broadly consistent with the CEF model of edge fracture. A finite volume method program was used to simulate the flow of a viscoelastic liquid, obeying the modified Phan-Thien-Tanner model, to obtain the velocity and stress distribution in parallel plate flow in three dimensions. Stress concentrations of the second normal stress difference (N ₂) were found in the plane of the crack; the velocity distribution shows a secondary flow tending to aid crack formation if N ₂ is negative, and a secondary flow tending to suppress crack formation if N ₂ is positive. Received: 4 January 1999 Accepted: 19 May 1999     

Fatigue crack extension behaviour in porous plasma spray 80%Ni---20%Cr material: The influence of R-ratio

This paper attempts to describe the fatigue crack growth response of a plasma spray 80%---20%Cr material, utilised in the corrosion protection of engineering components, whose microstructure consisted of (i) an austenitic matrix, (ii) a secondary dispersion of a chromite non-metallic inclusion phase and (iii) regions of closed and connected porosity.It was demonstrated that little or no effect of R-ratio was observed on the threshold stress intensity range ΔK_th, which was attributed to both the materials fine to intermediate grain size and probable plain stress testing conditions which significantly decrease crack closure effects. At intermediate fatigue crack growth rates high ratio results were an order of magnitude faster than the low R-ratio data. This was the result of the high R-ratio case seeking out more regions of porosity which then increased the local ΔK on the remaining ligaments leading to accelerated crack growth rates.Porosity was shown to significantly decrease the value of ΔK_th and the extent of porosity observed on fatigue fracture surfaces increased with ΔK level and was well in excess of that of 5% recorded by metallography. Hence the growing fatigue crack preferentially sought out regions of porosity as they represented locations of low fracture energy.     

Near threshold fatigue behaviour of flake graphite cast irons microstructures

This paper describes the influence of material toughness degradation, through reversed temper embrittlement (RTE) and mean stress on the near threshold fatigue crack growth characteristics of a CrMoV turbine bolting steel at ambient and elevated temperatures. It was established at ambient temperatures that strong effects of R-ratio and material condition (toughness) were observed on near threshold fatigue crack growth characteristics. At elevated temperatures it was shown that for the non-embrittled material that only under low R-ratio conditions did increased temperature increase the level of threshold stress intensity ΔK_th, by some 20%. In the case of embrittled material, increasing the temperature increased ΔK_th levels by around 30% and decreased near threshold growth rates by an order of magnitude at low to intermediate R-ratio levels.The effects of R-ratio on ΔK_th for all material and mechanical testing conditions could be simply expressed by the difference between ΔK_th at R = O and a constant B multiplied by R.Quantitative fractographic observations indicated that, generally, the incidence of intergranular failure prevalent in embrittled and non-embrittled steels exhibited a maximum at some specific ΔK level. Also in embrittled steels large effects of environmental assisted crack (EAC) growth were observed at near threshold fatigue crack growth rates. It was suggested that this was the result of the much reduced material cohesive strength which was caused by the presence of both impurity and hydrogen atoms.