Influence of singularity dominated zone for propagating cracks in finite size specimens

The singularity dominated zones for straight as well as curved cracks propagating in finite size specimens were determined experimentally by using the optical method of dynamic photoelasticity using the near-field stress equations. Experimental data was carefully analyzed using improved numerical schemes to get the complete stress field around the propagating crack. This stress field was critically examined to evaluate the size of singularity dominated zones for cracks propagating in straight as well as curved paths. For this purpose, the exact solution was compared with the singular solution using stress components σ_x, σ_y, τ_xy and the maximum shear stress τ_max as a criterion respectively. For straight cracks where the stress field is symmetric about the crack path, the singularity dominated zones can be determined by using any one of the stresses. However, for a curved crack, the zones were unsymmetric. This study shows that σ_y, the crack opening stress, yields the best result for characterizing the singularity dominated zone around a running crack tip.     

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.     

Fatigue crack opening stress based on the strip-yield model

The modified strip-yield model based on the Dugdale model and two-dimensional approximate weight function method were utilized to evaluate the effect of in-plane constraint, transverse stress, on the fatigue crack closure. The plastic zone sizes and the crack opening stresses considering transverse stress were calculated for four specimens: single edge-notched tension (SENT) specimen, single edge-notched bend (SENB) specimen, center-cracked tension (CCT) specimen, double edge-notched tension (DENT) specimen under uniaxial loading. And the crack opening behavior of the center-cracked specimen under biaxial loading was also evaluated. Normalized crack opening stresses σ_op/σ_max for four specimens were successfully described by the normalized plastic zone parameter Δω^′_rev/ω^′ considering transverse stress, where Δω^′_rev and ω^′ are the size of the reversed plastic zone at the moment of first crack tip closure and the size of the forward plastic zone for maximum stress, respectively. The normalized plastic zone parameter with transverse stress also was satisfactorily correlated with the behavior of crack closure for CCT specimen under biaxial loading.     

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.     

A field model interpretation of crack initiation and growth behavior in ferroelectric ceramics: change of poling direction and boundary condition

Strain energy density expressions are obtained from a field model that can qualitatively exhibit how the electrical and mechanical disturbances would affect the crack growth behavior in ferroelectric ceramics. Simplification is achieved by considering only three material constants to account for elastic, piezoelectric and dielectric effects. Cross interaction of electric field (or displacement) with mechanical stress (or strain) is identified with the piezoelectric effect; it occurs only when the pole is aligned normal to the crack. Switching of the pole axis by 90° and 180° is examined for possible connection with domain switching. Opposing crack growth behavior can be obtained when the specification of mechanical stress σ^∞ and electric field E^∞ or (σ^∞,E^∞) is replaced by strain ε^∞ and electric displacement D^∞ or (ε^∞,D^∞). Mixed conditions (σ^∞,D^∞) and (ε^∞,E^∞) are also considered. In general, crack growth is found to be larger when compared to that without the application of electric disturbances. This includes both the electric field and displacement. For the eight possible boundary conditions, crack growth retardation is identified only with (E_y^∞,σ_y^∞) for negative E_y^∞ and (D_y^∞,ε_y^∞) for positive D_y^∞ while the mechanical conditions σ_y^∞ or ε_y^∞ are not changed. Suitable combinations of the elastic, piezoelectric and dielectric material constants could also be made to suppress crack growth.     

Crack reinforcement by distributed springs

An analysis is presented for a centre crack which is reinforced by linear springs over a distance l from both tips. The reinforcement introduces a characteristic length k⁻¹ where k denotes a suitably defined spring constant. Two functions F, V are defined corresponding respectively to the stress intensity factor and the maximum spring-stretch, normalized relative to their values in the absence of reinforcement, for modes I, II or III. Asymptotic expansions for these functions are derived for the limiting cases of soft springs (kl ⪡ 1) and hard springs (kl ⪢ 1). Interpolating functions constructed from these asymptotic expansions are shown to agree with the numerical solution of the governing integral equation over the intermediate range of kl, to within an error of less than 1%. It is noted that the leading term of the expansions for hard springs can be derived by physical reasoning, requiring relatively simple calculations, which can also be used for nonlinear springs.     

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.     

Limitation on crack speed for a strip-craze model applied to PMMA

A strip-craze model is proposed to study crack propagation in polymers. A nonlinear differential equation is derived to govern the dynamic process of crack propagation. The viscous feature of the material in the craze zone is taken into account by means of an experimentally determined relationship between the craze stress and crack speed. By fitting experimental data of PMMA into the model, some parameters including the strip-craze length are deduced. A non-singular stress is introduced to control the crack propagation with a strip craze at its tip. Variations of the crack length and the crack speed with time are computed and their dependence on the non-singular stress is investigated. For PMMA, three stages of crack propagation are identified in terms of initial non-singular stress σ_ns0. When σ_ns0<60 MPa, the crack speed mm/s and the crack is basically stationary; when 60 <σ_ns0<95 MPa, then mm/s the crack is in slow propagation; when σ_ns0>95 MPa, then mm/s and the crack is in rapid propagation. The proposed model is applicable only in slow crack propagation.     

Crack front curvature effect on stress intensity factor of through and part-through thickness cracks

The character of the local stresses and displacements are determined for a through crack with finite radius of curvature in a finite thickness plate. Numerical results obtained from the boundary element method show that the solutions are sufficiently accurate for /a ≤ 0.03 and 0.03 ≤ /a ≤ 0.1, where and a represent, respectively, the crack front radius of curvature and crack dimension such that a is the width of a through thickness crack and the depth of a part-through crack. For /a ≤ 0.1, the asymptotic singular stress field dominates such that the Mode I stress intensity factor K₁ can be evaluated. As the crack border radius of curvature is increased for /a ≥ 0.1, the non-singular terms become significant such that K_I would no longer dominate. Other failure criteria would have to be invoked to address fracture initiation.     

Fatigue life prediction of out-of-plane gusset welded joints using strain energy density factor approach

Fatigue growth behavior of out-of-plane gusset welded joints is studied using the strain energy density factor approach. Fatigue tests on two types of specimens with curvatures of ρ = 0 and ρ = 30 were performed in order to estimate fatigue strength under tension. Fatigue crack growth analysis is carried out to show the effects of initial crack shape, initial crack length and stress ratio. Fatigue crack growth parameters were obtained from crack growth curves assuming constant crack shapes. The results of analysis for the assumed crack shapes agreed well with the experimental data. Fatigue propagation life of the ρ = 30 specimen was larger than that of the ρ = 0 specimen.     

Near tip field for pseudo Mode II crack growth: Power law hardening material

The asymptotic behavior of stress and strain near the tip of a Mode II crack growing in power law hardening material is analyzed by assuming that the crack grows straight ahead even though tests show otherwise. The results show that the stress and strain possess the singularities of (ln r)^2/(n−1) and (ln r)^2n/(n−1) respectively. The distance from the crack tip is r, and n is the hardening exponent, i.e. σⁿ. The amplitudes of the stress and strain near the crack tip are determined by the asymptotic analysis.     

Isoparametric shear spring element applied to crack patching and instability

In this work the isoparametric shear spring element is applied to the stress and energy analysis of a center-crack panel reinforced by a rectangular patch. In this model, only transverse shears are assumed to prevail in the adhesive layer. The stresses and crack-tip stress intensity factors are obtained for reinforcement on both sides and one side of the panel, and are found to be in agreement with those obtained by previous authors using the triangular shear spring element.Crack stability that tends to vary with patch thickness is determined from the local and global maximum of the minimum strain energy density function denoted, respectively, as [(dW/dV)_min^max]_L at point L and [(dW/dV)_min^max]_G at point G. The distance l between L and G gives the prospective path of subcritical crack growth and its magnitude provides a measure of the degree of crack stability. A patched panel with small l tends to be more stable than that with large l. By increasing the patch thickness beyond a certain value, l can be contained within the patch such that failure, if initiated, will be highly localized. Such a behavior is exhibited. Numerical results are provided to support the foregoing conclusion.     

Crack size and speed interaction characteristics at micro-, meso- and macro-scale

The motivation to examine physical events at even smaller size scale arises from the development of use-specific materials where information transfer from one micro- or macro-element to another could be pre-assigned. There is the growing belief that the cumulated macroscopic experiences could be related to those at the lower size scales. Otherwise, there serves little purpose to examine material behavior at the different scale levels. Size scale, however, is intimately associated with time, not to mention temperature. As the size and time scales are shifted, different physical events may be identified. Dislocations with the movements of atoms, shear and rotation of clusters of molecules with inhomogeneity of polycrystals; and yielding/fracture with bulk properties of continuum specimens. Piecemeal results at the different scale levels are vulnerable to the possibility that they may be incompatible. The attention should therefore be focused on a single formulation that has the characteristics of multiscaling in size and time. The fact that the task may be overwhelmingly difficult cannot be used as an excuse for ignoring the fundamental aspects of the problem.Local nonlinearity is smeared into a small zone ahead of the crack. A “restrain stress” is introduced to also account for cracking at the meso-scale.The major emphasis is placed on developing a model that could exhibit the evolution characteristics of change in cracking behavior due to size and speed. Material inhomogeneity is assumed to favor self-similar crack growth although this may not always be the case. For relatively high restrain stress, the possible nucleation of micro-, meso- and macro-crack can be distinguished near the crack tip region. This distinction quickly disappears after a small distance after which scaling is no longer possible. This character prevails for Mode I and II cracking at different speeds. Special efforts are made to confine discussions within the framework of assumed conditions. To be kept in mind are the words of Isaac Newton in the Fourth Regula Philosophandi:

“Men are often led into error by the love of simplicity which disposes us to reduce things to few principles, and to conceive a greater simplicity in nature than there really is … We may learn something of the way in which nature operates from fact and observation; but if we conclude that it operates in such a manner, only because to our understanding that operates to be the best and simplest manner, we shall always go wrong.”––Isaac Newton

Effect of grain size on slow fatigue crack propagation and plastic deformation near crack tip

Fatigue crack growth and its threshold are investigated at a stress ratio of 0.5 for the three-point bend specimen made of Austenitic stainless steel. The effect of grain size on the crack tip plastic deformation is investigated. The results show that the threshold value Δk_th increases linearly with the square root of grain size d and the growth rate is slower for materials with larger grain size. The plastic zone size and ratio for different grain sizes are different at the threshold. The maximum stress intensity factor is k_max and σ_ys is the yield strength. At the same time, the characteristics of the plastic deformation development is discontinuous and anti-symmetric as the growth rate is increased from 2·10^—8 to 10⁻⁷ mm/cycle.A dimensionless relation of the form for collating fatigue crack starting growth data is proposed in which Δk_th represents the stress intensity factor range at the threshold. Based on experimental results, this relation attains the value of 0.6 for a fatigue crack to start growth in the Austenitic stainless steel investigated in this work. Metallurgical examinations were also carried out to show a transgranular shear mode of cyclic cleavage and plastic shear.     

Crack repair using an elastic filler

The effect of repairing a crack in an elastic body using an elastic filler is examined in terms of the stress intensity levels generated at the crack tip. The effect of the filler is to change the stress field singularity from order 1/r^1/2 to 1/r^(1-λ) where r is the distance from the crack tip, and λ is the solution to a simple transcendental equation. The singularity power (1-λ) varies from (the unfilled crack limit) to 1 (the fully repaired crack), depending primarily on the scaled shear modulus ratio γ_r defined by G₂/G₁=γ_rε, where 2πε is the (small) crack angle, and the indices (1, 2) refer to base and filler material properties, respectively. The fully repaired limit is effectively reached for γ_r≈10, so that fillers with surprisingly small shear modulus ratios can be effectively used to repair cracks. This fits in with observations in the mining industry, where materials with G₂/G₁ of the order of 10^-3 have been found to be effective for stabilizing the walls of tunnels. The results are also relevant for the repair of cracks in thin elastic sheets.     

Micro/macro-crack growth due to creep–fatigue dependency on time–temperature material behavior

The implicit character of micro-structural degradation is determined by specifying the time history of crack growth caused by creep–fatigue interaction at high temperature. A dual scale micro/macro-equivalent crack growth model is used to illustrate the underlying principle of multiscaling which can be applied equally well to nano/micro. A series of dual scale models can be connected to formulate triple or quadruple scale models. Temperature and time-dependent thermo-mechanical material properties are developed to dictate the design time history of creep–fatigue cracking that can serve as the master curve for health monitoring.In contrast to the conventional procedure of problem/solution approach by specifying the time- and temperature-dependent material properties as a priori, the desired solution is then defined for a class of anticipated loadings. A scheme for matching the loading history with the damage evolution is then obtained. The results depend on the initial crack size and the extent of creep in proportion to fatigue damage. The path dependent nature of damage is demonstrated by showing the range of the pertinent parameters that control the final destruction of the material. A possible scenario of 20 yr of life span for the 38Cr2Mo2VA ultra-high strength steel is used to develop the evolution of the micro-structural degradation. Three micro/macro-parameters μ^*, d^* and σ^* are used to exhibit the time-dependent variation of the material, geometry and load effects. They are necessary to reflect the scale transitory behavior of creep–fatigue damage. Once the algorithm is developed, the material can be tailor made to match the behavior. That is a different life span of the same material would alter the time behavior of μ^*, d^* and σ^* and hence the micro-structural degradation history. The one-to-one correspondence of the material micro-structure degradation history with that of damage by cracking is the essence of path dependency. Numerical results and graphs are obtained to demonstrate how the inherently implicit material micro-structure parameters can be evaluated from the uniaxial bulk material properties at the macroscopic scale.The combined behavior of creep and fatigue can be exhibited by specifying the parameter ξ with reference to the initial defect size a₀. Large ξ (0.90 and 0.85) gives critical crack size a_cr = 11–14 mm (at t < 20 yr) for a₀ about 1.3 mm. For small ξ (0.05 and 0.15), there results critical a_cr = 6–7 mm (at t < 20 yr) for a₀ about 0.7–0.8 mm. The initial crack is estimated to increase its length by an order of magnitude before triggering global to the instability. This also applies ξ ≈ 0.5 where creep interacts severely with fatigue. Fine tuning of a_cr and a₀ can be made to meet the condition oft = 20 yr.Trade off among load, material and geometric parameters are quantified such that the optimum conditions can be determined for the desired life qualified by the initial–final defect sizes. The scenario assumed in this work is indicative of the capability of the methodology. The initial–final defect sizes can be varied by re-designing the time–temperature material specifications. To reiterate, the uniqueness of solution requires the end result to match with the initial conditions for a given problem. This basic requirement has been accomplished by the dual scale micro/macro-crack growth model for creep and fatigue.     

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).     

Plane-strain crack tip field and a dual-parameter fracture-criterion for non-linear fracture mechanics

Hancock and Cowling measured the critical crack tip opening displacements, δ_f, at fracture initiation in HY-80 steel specimens of six different configurations. δ_f varied from 90 μm in a deeply double-edge-cracked tensile panel to 900 μm in a single-edge-cracked tensile panel.McMeeking and Parks, and Shih and German have shown by their finite element calculations that the characteristics of the plane strain crack tip fields in both large scale yielding and general yielding are strongly dependent on specimen geometry and load level.In this study, the plane strain crack tip fields in the specimens tested by Hancock and Cowling were calculated using the finite element method. The crack tip triaxial tensile stress field is strongly affected by specimen geometric constraint, and the state of the triaxial tensile stress in a crack tip region is monitored by the ratio between the local tensile stress and the effective stress, i.e., ( ), at a distance x=2δ from the crack tip. The values of ( ) vary from 3.1 for the double-edge-cracked tensile panel to 1.7 for the single-edge-cracked tensile panel. The δ_f measured by Hancock and Cowling correlates very well with the ratio ( ). δ_f is a measure of the fracture ductility of the material ahead of the crack tip, and the ductility decreases with an increase in the triaxial tensile stress, i.e., the ratio ( ).