An analysis of the mixed-type stress field near a crack-tip region

This paper proposes a complete analytical solution for stress field near a crack-tip in closed form to analyse the effects of load biaxiality and is concerned with the following problems: stress intensity factor of pure type-I and type-II fracture; patterns of maximum shear isostats (isochromatic fringe patterns); patterns of octahedral stresses σ_i and small-scale elasto-plastic stress field; the direction of initial crack extension and the biaxial effects on the crack with a blunt tip.     

A simple method for the photoelastic determination of mode I stress intensity factors

A new simple method for the photoelastic determination of Mode I stress intensity factors from isochromatics is proposed. This method takes into account the fact that a considerable part of the error committed in the photoelastic determination of Mode I stress intensity factors K_I at crack tips, based on experimentally obtained isochromatic fringe patterns, is due to ignoring the non-singular part of the stress field near the crack tips for the evaluation of these factors. This error can, in most cases, be minimized by an appropriate selection of the polar direction from the crack tip on which the experimental measurements for the subsequent evaluation of the stress intensity factors K_I are made. The suitable polar direction for determining K_I depends in general on the distance of the point where measurements on the isochromatics are made from the crack tip. The method was applied to the problem of a simple crack inside an infinite medium under uniaxial and biaxial loading. A comparison of the present method whith the employed analogous methods shows the superiority of the proposed method.     

Complex variable method of computing Jk for bi-material interface cracks

A method using functions of a complex variable is developed for evaluation of J₁ and a modified J₂ integrals for bi-material interface cracks. This method, used in conjunction with the finite element method, would be useful in the prediction of stress intensity factors for cracks lying between the interface of two dissimilar materials. Since the direct evaluation of J₂ poses difficulties in modeling the singular behavior in the near vicinity around the crack tip for bi-material crack problems, it is modified by evaluating it around a contour path of small radius from the crack tip within the singularity dominated zone. It is shown that the stress intensity factors for a bi-material interface crack can be accurately evaluated using these j_k integrals.     

Crack-tip fields for matrix cracks between dissimilar elastic and creeping materials

The stress fields near the tip of a matrix crack terminating at and perpendicular to a planar interface under symmetric in-plane loading in plane strain are investigated. The bimaterial interface is formed by a linearly elastic material and an elastic power-law creeping material in which the crack is located. Using generalized expansions at the crack tip in each region and matching the stresses and displacements across the interface in an asymptotic sense, a series asymptotic solution is constructed for the stresses and strain rates near the crack tip. It is found that the stress singularities, to the leading order, are the same in each material; the stress exponent is real. The oscillatory higher-order terms exist in both regions and stress higher-order term with the order of O(r°) appears in the elastic material. The stress exponents and the angular distributions for singular terms and higher order terms are obtained for different creep exponents and material properties in each region. A full agreement between asymptotic solutions and the full-field finite element results for a set of test examples with different times has been obtained.     

Singular behaviour at the end of a crack yielding on inclined slip-planes

Based on the Williams' method of the eigenfunction expansions, the singular stress distribution near a mode I crack tip yielding on two inclined slip-planes is obtained in the form with an undetermined constant. It is found that stresses have the singularity of the form r ^λ-1, where r is the distance from the crack tip, and the exponent λ-1 lies between-0.5 and 0, depending on the inclination angle of the slip-planes. The undetermined constant is then estimated by combining our result with an approximate solution obtained by J.R. Rice. The case of a mode III crack is also analyzed and the results are compared with a known closed-form solution and discussed.     

On the dominance of the singular dynamic crack tip stress field under high rate loading

The dynamic stress field near a running crack tip is investigated using photoelasticity and caustics. Both quasistatic and dynamic impact loads were considered. Under impact loading, it was found that the range of dominance of the singular term in the asymptotic crack tip stress field expansion was very small. The need for considering higher order transient terms in interpreting the isochromatic fringes is demonstrated. Also, the importance of considering the effect of the higher order transient terms in the characterization of failure criteria is discussed.     

Analysis of the effects of compressive stresses on fatigue crack propagation rate

In this study, the effects of compressive stresses on the crack tip parameters and its implication on fatigue crack growth have been studied. Elastic–plastic finite element analysis has been used to analyse the change of crack tip parameters with the increase of the applied compressive stress level.The near crack tip opening displacements and the reverse plastic zone size around the crack tip have been obtained. The finite element analysis shows that when unloading from peak tensile applied stress to zero applied stress, the crack tip is still kept open and the crack tip opening displacement gradually decreases further with the applied compressive stress. It has been found that for a tension–compression stress cycle these crack tip parameters are determined mainly by two loading parameters, the maximum stress intensity K_max in the tension part of the stress cycle and the maximum compressive stress σ_maxcom in the compression part of the stress cycle.Based on the two parameters, K_max, and σ_maxcom, a fatigue crack propagation model for negative R ratios only has been developed to include the compressive stress effect on the fatigue crack propagation rate.Experimental fatigue crack propagation data sets were used for the verification of this model, good agreements have been obtained.     

A multilayer hybrid-stress finite element analysis of a through-thickness edge crack in A ± 45° laminate

A solution is given for the three-dimensional stress field near a through-thickness edge crack in a thin ± 45° laminate having elastic ply moduli typical of graphite/epoxy. The stress distribution was obtained by a three-dimensional multilayer finite element analysis based on the hybrid stress model, formulated through the minimum complementary energy principle. The results indicate that the in-plane stresses of each individual ply follow the classical $\text{1}\text{√r}$ stress singularity, but that the shape of isostress contours in the crack tip region is strongly distorted from predictions based on two-dimensional anisotropic fracture mechanics theory. The interlaminar shear stresses increase rapidly as the crack tip is approached, but are restricted to a local region around the crack tip and flanks. The interlaminar normal stress is assumed to be negligible in the formulation of the analysis.     

Calculation of K-factor and T-stress for cracks in anisotropic bimaterials

The problem of a crack in a thin layer terminating perpendicular to a layer/substrate interface is analyzed for a general case of elastic anisotropy. The crack is modelled by means of continuous distribution of dislocations, which is assumed to be singular at the crack tip. A system of simultaneous functional equations is obtained that enables to find the singularity exponent λ. The reciprocal theorem (ψ-integral) is used to compute the generalized stress intensity factor (GSIF) through the remote stress and displacement field for a particular specimen geometry and boundary conditions using FEM. The results obtained are compared with the evaluation of GSIF based upon the dislocation arrays technique. Existing semi-analytical solution for singularities in anisotropic trimaterials is applied and its validity for the specimen investigated is checked by FEM. The evaluation of T-stress using the dislocation arrays technique is performed.     

A Photoelastic Study of T‐stress in Centrally Cracked Brazilian Disc Specimen Under Mode II Loading

Abstract: In the traditional formulation of the stress field near a crack tip, the presence of the T‐stress is generally considered only under mode I or mixed mode I and II conditions. In this paper its presence in almost pure Mode II is experimentally investigated by mean of photoelasticity and its effects on the isochromatic fringe patterns are discussed. The test specimens are Brazilian discs containing sharp central cracks. After crack generation, all residual stresses are removed with thermal treatment of the specimens. Then, a compressive load is applied in a specific angle to induce mode II deformation. The observed isochromatic fringes show very good consistency with theoretical predictions. Experimental results indicate that this specimen has a negative T‐stress in mode II condition. The results calculated for K_II and T from photoelastic experiments agree well with numerical results available from finite element method.     

Multiscale continuum modeling of a crack in elastic media with microstructures

Cosserat type continuum theories have been employed by many authors to study cracks in elastic solids with microstructures. Depending on which theory was used, different crack tip stress singularities have been obtained. In this paper, a microstructure continuum theory is used to model a layered elastic medium containing a crack parallel to the layers. The crack problem is solved by means of the Fourier transform. The resulting integrodifferential equations are discretized using the Chebyshev polynomial expansion method for numerical solutions. By using the present theory, the explicit internal length effects upon the crack opening displacement and stress field can be observed. It is found that the stress field near the crack tip is not singular according to the microstructure continuum solution although the level of the opening stress shows an increasing trend until it gets very close to the crack tip. The rising portion of the near tip opening stress is used to project the stress intensity factor which agrees fairly well with that obtained using the FEM to perform stress analyses of the cracked layered medium with the exact geometry. The numerical solutions also indicate that treating the layered medium as an equivalent homogeneous classical elastic solid is not adequate if cracks are present and accurate stress intensity factors in the original layered medium is desired.     

A new method of determining the stress intensity factor K from isochromatic fringe loops

This paper describes the four parameter method of analysis for determining the stress intensity factor K. Dynamic photoelastic isochromatic fringe patterns associated with cracks propagating in centerpin loaded, eccentric-pin-loaded and crack-line-loaded SEN specimens of Homalite 100 were recorded. Data was obtained for tests over a range of crack velocities from arrest to the terminal velocity of 14,900 in. sec (378 m sec).Six measurements describing the size and shape of the experimental isochromatic loops were used to determine the stress intensity factor K by employing a comparison function to match analytical and experimental results. The analytical isochromatic loops were generated with a Westergaard stress function of the form $\text{Z(z) =}\text{K}\text{√2πz}\text{{1 + β(}\text{z}\text{a})}$ and a superimposed normal stress σ_Ox=αK|√2πz which acts parallel to the direction of crack extension. Results were obtained by the computer program (FRACTURE) for different values of the four parameters to give 8925 analytical fringe loops. Another computer program (SEARCH) was used to find a small group of solution which given very low values of the comparison function f_c. The final solution which contains the value of K was obtained from the small group by selecting the most consistent solution.The results obtained for Homalite 100 show that K_min and K_k are nearly the same and that ? increases abruptly from 0 to about 10,000 in./sec (254 msec) for modest increases in K above 400 psi √in. (4.4 × 10⁵$\text{Nm}{}^{\mathrm{<mtext>3</mtext><mtext>2</mtext>}}$). Further increases in crack velocity require significant increases in K until terminal velocity is     

Singular fields of a mode III interface crack in a power law hardening bimaterial

A high order of asymptotic solution of the singular fields near the tip of a mode III interface crack for pure power law hardening bimaterials is obtained by using the hodograph transformation. It is found that the zero order of the asymptotic solution corresponds to the assumption of a rigid substrate at the interface, and the first order of it is deduced in order to satisfy completely two continuity conditions of the stress and displacement across the interface in the asymptotic sense. The singularities of stress and strain of the zero order asymptotic solutions are –1/(n ₁+1) and –n/(n ₁+1) respectively (n=n ₁, n ₂ is the hardening exponent of the bimaterials). The applicability conditions of the asymptotic solutions are determined for both zero and first orders. It is proved that the Guo-Keer solution [23] is limited in some conditions. The angular functions of the singular fields for this interface crack problem are first expressed by closed form.     

Precise elastic-plastic analysis of crack line field for mode II plane strain crack

The near crack line analysis method has been used to investigate the exact elastic-plastic solutions of a mode II crack under plane strain condition in an elastic-perfectly plastic solid. The significance of this paper is that the assumptions of the conventional small scale yielding theory have been completely abandoned. The inappropriateness of matching conditions formerly taken at the elastic-plastic boundary ths been corrected as well. By eatching the general solution of the plastic stress (but not the special solution that was adopted) with the exact elastic stresses (but not the crack tip K-dominant field) at the elastic-plastic boundary near the crack line, the plastic stresses, the length of the plastic zone and the unit normal vector of the elastic-plastic boundary, which are sufficiently precise near the crack line region, have been given. The solutions are suitable not only under the condition that the plastic region is sufficiently small but also under the condition that the plastic region is large.     

Perturbation method for the solution of a Zener–Stroh crack with a slightly curved configuration

This paper investigates the Zener–Stroh crack with curved configuration in plane elasticity. A singular integral equation is suggested to solve the problem. Formulae for evaluating the SIFs and T-stress at the crack tip are suggested. If the curve configuration is a product of a small parameter and a quadratic function, a perturbation method based on the singular integral equation is suggested. In the method, the singular integral equation can be expanded into a series with respect to the small parameter. Therefore, many singular integral equations can be separated from the same power order for the small parameter. These singular integral equations can be solved successively. The solution of the successive singular integral equations will provide results for stress intensity factors and T-stress at the crack tip. It is found that the behaviors for the solution of SIFs and T-stress in the Zener–Stroh crack and the Griffith crack are quite different. This can be seen from the presented comparison results.     

Dominant stress region for crack initiation at interface edge of microdot on a substrate

In order to examine the mechanics of crack initiation at the free interface edge of a microcomponent on a substrate, delamination tests are carried out for two specimen shapes of Cr microdots on a SiO₂ substrate. The microdots of the first specimen are shaped like the frustum of a round cone. The Cr microdots are successfully delaminated from the SiO₂ substrate in a brittle manner and the critical load is measured by atomic force microscopy (AFM) with a lateral loading apparatus. Stress analysis reveals that a singular stress field exists near the interface edge and the strength for the crack initiation is governed by the intensified normal stress field. The critical stress intensity parameter is evaluated as K_σC ≈ 0.24 MPa m^0.39. Similar delamination tests are conducted for microdots shaped like the frustum of an oval cone. The stress distributions at the crack initiation of this specimen shape show a higher normal stress than the first specimen shape in the region near the interface edge of about x < 40 nm, while it is lower in the region of about x > 50 nm (x: distance from the edge). This suggests a limitation of conventional fracture mechanics: namely, the crack initiation in these specimens is not uniquely governed by the intensity of the singular field. It is found that the delamination crack is initiated when the averaged stress σ_ya in the region of 90-130 nm reaches 190-270 MPa, regardless of the specimen shape. This indicates that the dominant stress region of crack initiation is roughly estimated as 90-130 nm and the criterion is given in terms of the averaged stress in the region.     

Mixed mode fracture in power law hardening materials near Mode I

Mixed mode fracture in power law hardening materials near mode I loading conditions is investigated for the case of plane strain. It is demonstrated by application of full-field finite element analysis and two-parameter asymptotic analysis that the asymptotic mixed mode solution of Shih (1974), where the effects of arbitrary loading combine to enter the leading term of the asymptotic series, does not apply close to mode I. Instead, a `mode I dominant' higher order asymptotic form applies, where the leading term is the pure mode I symmetric HRR solution and the second term is antisymmetric. The significant difference between the two asymptotic solutions is the singular nature of the antisymmetric part of the stresses and strains. Therefore, near mode I the antisymmetric part of the local fields does not contribute to the J-integral. The transition from the mixed mode asymptotic solution to this two-term solution occurs when the ratio of the radial stress on the upper crack surface to that on the lower crack surface, i.e., σ<sub>rr</sub>(r,+π)/σ<sub>rr</sub>(r,?π), for the limit as r approaches zero, switches from ?1 to +1. Budiansky and Rice (1973) predicted this change in sign somewhere between mode I and mode II. This sign change corresponds to the singular nature of the antisymmetric fields switching from the leading HRR eigenvalue to a weaker second term eigenvalue. Based on full-field finite element results, the far-field loading for which this `jump' occurs is very difficult to identify, as it originates deep within the plastic zone at unrealistically small distances from the crack tip. At the physical length scale on the order of the crack tip opening displacement, the transition from the mode II dominant asymptotic solution to the mode I dominant solution appears to be continuous, not abrupt. Identification of this transition is estimated from the full-field results by making use of the two asymptotic solutions. These results identify a breakdown in the HRR theory for the representation of the antisymmetric part of the stresses for mixed mode fracture near mode I. This breakdown is explained by the switching from a mixed mode solution to a higher order solution. Unfortunately, the practical application of this higher-order solution is limited to hardening powers of approximately 2<n<3.     

Prediction of fatigue crack growth under constant amplitude loading and a single overload based on elasto-plastic crack tip stresses and strains

It is generally accepted that the fatigue crack growth (FCG) depends mainly on the stress intensity factor range (ΔK) and the maximum stress intensity factor (K_max). The two parameters are usually combined into one expression called often as the driving force and many various driving forces have been proposed up to date. The driving force can be successful as long as the stress intensity factors are appropriately correlated with the actual elasto-plastic crack tip stress-strain field. However, the correlation between the stress intensity factors and the crack tip stress-strain field is often influenced by residual stresses induced in due course.A two-parameter (ΔK_tot, K_max,tot) driving force based on the elasto-plastic crack tip stress-strain history has been proposed. The applied stress intensity factors (ΔK_appl, K_max,appl) were modified to the total stress intensity factors (ΔK_tot, K_max,tot) in order to account for the effect of the local crack tip stresses and strains on fatigue crack growth. The FCG was predicted by simulating the stress-strain response in the material volume adjacent to the crack tip and estimating the accumulated fatigue damage. The fatigue crack growth was regarded as a process of successive crack re-initiations in the crack tip region. The model was developed to predict the effect of the mean and residual stresses induced by the cyclic loading. The effect of variable amplitude loadings on FCG can be also quantified on the basis of the proposed model. A two-parameter driving force in the form of: was derived based on the local stresses and strains at the crack tip and the Smith-Watson-Topper (SWT) fatigue damage parameter: D = σ_maxΔε/2. The effect of the internal (residual) stress induced by the reversed cyclic plasticity manifested itself in the change of the resultant (total) stress intensity factors controlling the fatigue crack growth.The model was verified using experimental fatigue crack growth data for aluminum alloy 7075-T6 obtained under constant amplitude loading and a single overload.     

Dynamic fields near a crack tip growing in an elastic-perfectly-plastic solid

A complete asymptotic solution is given for the fields in the neighborhood of the tip of a steadily advancing crack in an incompressible elastic-perfectly-plastic solid.For Mode I crack growth in the plane strain condition, the following noteworthy results are revealed: (1) The entire crack tip in steady crack growth is surrounded by a plastic region, and no elastic unloading is predicted by the complete dynamic asymptotic solution. Thus, the elastic unloading region predicted by the result of neglecting the important influence of the inertia terms in the equations of motion. (2) Unlike the quasi-static solution, the dynamic solution yields strain fields with a logarithmic singularity everywhere near the crack tip. (3) The stress field varies throughout the entire crack tip neighborhood, but does display behavior which can be approximated by a constant field followed by an essentially centered-fan field and then by another constant field, especially for small crack growth speeds. Indeed, the stress field reduces to that for the stationary crack, as the crack tip velocity - measured by the Mach number, M - reduces to zero; the strain field, however, does not reduce to that for the static solution, as M vanishes. (4) There are two shock fronts emanating from the crack tip across which certain stress and strain components undergo jump discontinuities. The location of the shock fronts and the magnitude of the jumps depend on the crack growth speed. The stress jump vanishes while the strain jump becomes unbounded, as the crack tip speed goes to zero.Finally, the Mode III steady-state crack growth is reviewed and, on the basis of Mode I and Mode III results, it is concluded that ductile fracture criteria for nonstationary cracks must be based on solutions which include the inertia effects, and that for this purpose, quasi-static solutions may be inadequate. Then, a possible ductile fracture criterion is suggested and discussed.One interesting feature of the complete dynamic asymptotic solution is that, unlike the quasi-static solution, it yields the same strain singularity for all three fracture modes.     

Determination of higher order coefficients and zones of dominance using a singular integral equation approach

A method to determine higher order coefficients from the solution of a singular integral equation is presented. The coefficients are defined by , which gives the radial stress at a distance, r, in front of the crack tip. In this asymptotic series the stress intensity factor, k₀, is the first coefficient, and the T-stress, T₀, is the second coefficient. For the example of an edge crack in a half space, converged values of the first 12 mode I coefficients (k_n and T_n, n = 0, … , 5) have been determined, and for an edge crack in a finite width strip, the first six coefficients are presented. Coefficients for an internal crack in a half space are also presented. Results for an edge crack in a finite width strip are used to quantify the size of the k-dominant zone, the kT-dominant zone and the zones associated with three and four terms, taking into account the entire region around the crack tip.