Revisit to the determination of stress-intensity factors and J-integrals using the caustics method

Applications of the optical (shadow) method of reflective caustics to measurement of the stress-intensity factor andJ integral in various specimens are investigated. The necessary experimental requirements to help determine accurate stress-intensity factor andJ integral are described. The ratios ofr _o (radius of initial curve)/r _p (plastic-zone size) andr _o/t(thickness of specimen) are found to be very important experimental parameters to obtain meaningful stress and/or strain intensities surrounding crack tips. The appropriate ranges to determine accurate values of stress-intensity factor andJ integral for polycarbonate (compact tension) and aluminum (c-shaped tension) specimens are presented. Paper was presented at the 1992 SEM Spring Conference on Experimental Mechanics held in Las Vegas, NV on June 8–11.     

Elastic-plastic analysis of biaxially loaded center-cracked plates

Center-cracked panels loaded in biaxial tension are examined in this paper. Calibration relations for the J integral and the Q constraint factor are presented for a Ramberg–Osgood power law hardening material under plane stress and plane strain loadings. Two cases are examined: an isolated crack and a periodic array of cracks both under biaxial loading conditions. The latter has previously been studied for plane stress conditions. A number of different J estimation schemes are proposed based on the remote load and displacement and their dependence on geometry, biaxiality, and material properties is discussed. The variation of constraint, as characterised by Q, is also presented for plane stress and plane strain conditions. Simple slip line field solutions are derived for perfectly plastic conditions and the resulting limit load solutions are compared with numerically determined values. Implications for failure of cracked plates under biaxial loading are discussed.     

Dynamic response for functionally graded materials with penny-shaped cracks

This paper provides a method for studying the penny-shaped cracks configuration in functionally graded material(FGM) structures subjected to dynamic or steady loading. It is assumed that the FGMs are transversely isotropic and all the material properties only depend on the axial coordinatez. In the analysis, the elastic region is treated as a number of layers. The material properties are taken to be constants for each layer. By utilizing the Laplace transform and Hankel transform technique, the general solutions for the layers are derived. The dual integral equations are then obtained by introducing the mechanical boundary and layer interface conditions via the flexibility/stiffness matrix approach. The stress intensity factors are computed by solving dual integral equations numerically in Laplace transform domain. The solution in time domain is obtained by utilizing numerical Laplace inverse. The main advantage of the present model is its ability for treating multiple crack configurations in FGMs with arbitrarily distributed and continuously varied material properties by dividing the FGMs into a number of layers with the properties of each layer slightly different from one another. This work was supported by Failure Mechanics Laboratory of State Education Commission and the Post-doctor Research Fund of China.     

A domain independent integral expression for the crack extension force of a curved crack in three dimensions

An integral expression that is domain independent in curvilinear coordinates and compatible with zero divergence of Eshelby's (Phil. Trans. Roy. Soc. (London) 244 (1951) 87.) energy momentum tensor was obtained from the principle of virtual work. By applying Eshelby's definition of the force on a material defect a general expression of the crack extension force for a curved crack in three dimensions, here called the F-integral, was derived from the domain independent integral expression. The F-integral is given explicitly for a number of curved cracks and found to be in agreement with previously known solutions, when available. The influence of crack surface and crack front curvature upon the various forms of the F-integral is discussed. The F-integral presented in this work is a generalisation of the J-integral (Rice, J. Appl. Mech. 35 (1968) 379.) to curved cracks in orthogonal curvilinear coordinates.     

Computation of thermal fracture parameters for orthotropic functionally graded materials using Jk-integral

This article introduces a computational method based on the J_k-integral for mixed-mode fracture analysis of orthotropic functionally graded materials (FGMs) that are subjected to thermal stresses. The generalized definition of the J_k-integral is recast into a domain independent form composed of line and area integrals by utilizing the constitutive relations of plane orthotropic thermoelasticity. Implementation of the domain independent J_k-integral is realized through a numerical procedure developed by means of the finite element method. The outlined computational approach enables the evaluation of the modes I and II stress intensity factors, the energy release rate, and the T-stress. The developed technique is validated numerically by considering two different problems, the first of which is the problem of an embedded crack in an orthotropic FGM layer subjected to steady-state thermal stresses; and the second one is that of periodic cracks under transient thermal loading. Comparisons of the mixed-mode stress intensity factors evaluated by the J_k-integral based method to those calculated through the displacement correlation technique (DCT) and to those available in the literature point out that, the proposed form of the J_k-integral possesses the required domain independence and leads to numerical results of high accuracy. Further results are presented to illustrate the influences of the geometric and material constants on the thermal fracture parameters.     

动态载荷下功能梯度复合材料的圆币形裂纹问题

研究了动态载荷下功能梯度材料中的圆币形裂纹问题.假设材料为横观各向同性,并且含有多个垂直于厚度方向的裂纹,材料参数沿轴向(与裂纹面垂直的方向)为变化的,沿该方向将材料划分为许多单层,各单层材料参数为常数,利用Hankel变换法,在Laplace域内推导出了控制问题的对偶积分方程组.利用Laplace数值反演,得出了裂纹尖端的动态应力强度因子和能量释放率.研究了含两个裂纹的功能梯度接头结构,分析了材料非均匀性参数对应力强度因子和能量释放率的影响.     

Experimental determination ofK I for surface cracks in a square tube under torsion

This research used the fatigue-crack-growth-rate experimental method to obtainK _I for initial semi-elliptical 45-deg inclined, longitudinal and transverse surface cracks in a hollow square tube subjected to pure torsion. Both side and corner surface cracks were investigated. Reference crack-growth-rate data using CT and SEN specimens taken from the mild-steel tubes were used in calculations. All cracks propagated in mode I on the surface at 45 deg to the longitudinal axis in the plane of the maximum tensile stress. Mixed-mode crack growth occurred subsurface for the longitudinal and transverse cracks. The experimental results were compared with other semi-elliptical surface crackK _I solutions.     

Crack separation energy-rates for the DBCS model under biaxial modes of loading

A Modified version of the Dugdale-Bilby-Cottrell-Swinden (DBCS) model simulating the effect of plasticity at the tip of a crack in an infinite region was used by kfouri and rice (1978) to calculate the crack separation energy-rate G^Δ corresponding to a finite crack growth step Δa during plane strain mode I crack extension. The loading consisted of a remote uniaxial tension σp applied normally to the plane of the crack. Using Rice's path-independent integral J to characterize the applied load in the crack tip region, and assuming the length R of the crack tip plastic zone to be small compared with the length of the crack, an analytical expression was derived relating the ratios (J/G^Δ) and (2a/R) for small values of (2a/R). The analytical solution was incomplete in itself in that the value assumed in the plastic region of the DBCS model for the normal stress Y acting on the extending crack surfaces was the product of the yield stress in uniaxial tension σ_Y and an unknown parameter C, the value of which depends on the effect of the local hydrostatic stresses in the case of plane strain conditions. The analytical solution was compared with a numerical solution obtained from a plane strain elastic-plastic finite element analysis on a centre-cracked plate (CCP) of material obeying the von Mises yield criterion. The value used for the yield stress was 310 MN/m² and moderate isotropic linear hardening was applied with a tangent modulus of 4830 MN/m². A uniaxial tension σp was applied on the two appropriate sides of the plate. The comparisons showed that the analytical and finite element solutions were mutually consistent and they enabled the value of C to be established at 1.91. In the present paper similar comparisons are made between the analytical solution and the finite element solutions for the CCP of the same material under different biaxial modes of loading. By assuming further that the form of the analytical solution does not depend on the details of the geometry and of the loading at remote boundaries, a comparison has also been made with the results of a finite element analysis on a compact tension specimen (CTS) made of the same material as the CCP. The different values of C obtained in each case are consistent with investigations by other authors on the effect of load biaxiality on crack tip plasticity.     

Dynamic stress-intensity-factor determination from isopachics

The governing equations for determination of dynamic stress-intensity factor at the tip of a running crack are developed from a dynamic analysis of dynamic isopachicfringe patterns. The equations are applied to investigate dynamic crack propagation in Homalite 100 by means of dynamic holographic interferometry. A simple method based on simultaneous measurement of the widths of two isopachics allows determination of Irwin's additional stress field, and a dynamic correction function for the stress-intensity factor is derived. It was found that dynamicK-values obtained from dynamic isopachic-fringe-pattern analysis are lower than their corresponding static values. This implies a modification of the crack velocity vs. stress-intensity-factor relationship towards lower values ofK for dynamic crack propagation.     

Matched asymptotic expansions in nonlinear fracture mechanics—II. Longitudinal shear of an elastic work-hardening plastic specimen

Earlier analysis given by T.M. Edmunds and J.R. Willis (1976) is extended to deal with cracks in elastic work-hardening plastic specimens subjected to longitudinal shear loads. Solutions are expressed in terms of a set of parameters that are determined from linear elastic solutions alone. It is proved, for any specimen geometry and any loading symmetric about the plane of the crack, that a ‘plastic-zone correction’, obtained by solving a linear elastic problem for a crack which is a length r_y longer than the actual crack, provides a two-term asymptotic expansion for the J-integral, if r_y is defined suitably in terms of the linear elastic stress concentration factor and the initial slope of the work-hardening curve. The general method is applied in detail to a strip of finite width containing an edge crack, for which the effect of the work-hardening on the maximum extent of the plastic zone and on the J-integral is summarized graphically.     

Numerical computation and analytical estimation of stress-intensity factors for strips with multiple edge cracks

In this paper, stress-intensity factors for a two-dimensional problem are determined. Strips with multiple symmetrical edge cracks in tension are investigated. A simple analytical estimation is compared to numerical results. The influence of penetration of the crack faces and mixed-mode loading on the numerical results is investigated. A simple method to estimate stress-intensity factors for strips with multiple edge cracks is proposed.     

Dynamic curving cracks in functionally graded materials under thermo-mechanical loading

Mixed-mode dynamic crack growth along an arbitrarily smoothly varying path in functionally graded materials (FGMs) under thermo-mechanical loading is studied. The property gradation in FGMs is considered by varying shear-modulus, mass density, thermal conductivity and coefficient of thermal expansion exponentially along the gradation direction. Asymptotic analysis in conjunction with displacement potentials is used to develop the stress fields around propagating cracks in FGMs. Asymptotic temperature fields are developed first for the exponential variation of thermal conductivity and later these temperature fields are used to derive thermo-mechanical stress fields for a curving crack in FGMs. Using these thermo-mechanical stress fields, various components of the stresses are developed and the effect of curvature parameters, temperature and gradation on these stresses are discussed. Finally, using the minimum strain energy density criterion, the effect of curvature parameters, crack-tip speeds, non-homogeneity values and temperature gradients on crack growth directions are determined and discussed.     

Evaluation of stress intensity factors for bi-material interface cracks using displacement jump methods

The aim of the present work is to investigate the numerical modeling of interfacial cracks that may appear at the interface between two isotropic elastic materials. The extended finite element method is employed to analyze brittle and bi-material interfacial fatigue crack growth by computing the mixed mode stress intensity factors (SIF). Three different approaches are introduced to compute the SIFs. In the first one, mixed mode SIF is deduced from the computation of the contour integral as per the classical J-integral method, whereas a displacement method is used to evaluate the SIF by using either one or two displacement jumps located along the crack path in the second and third approaches. The displacement jump method is rather classical for mono-materials, but has to our knowledge not been used up to now for a bi-material. Hence, use of displacement jump for characterizing bi-material cracks constitutes the main contribution of the present study. Several benchmark tests including parametric studies are performed to show the effectiveness of these computational methodologies for SIF considering static and fatigue problems of bi-material structures. It is found that results based on the displacement jump methods are in a very good agreement with those of exact solutions, such as for the J-integral method, but with a larger domain of applicability and a better numerical efficiency (less time consuming and less spurious boundary effect).     

The nonlocal solution of two parallel cracks in functionally graded materials subjected to harmonic anti-plane shear waves

In this paper, the dynamic interaction of two parallel cracks in functionally graded materials (FGMs) is investigated by means of the non-local theory. To make the analysis tractable, the shear modulus and the material density are assumed to vary exponentially with the coordinate vertical to the crack. To reduce mathematical difficulties, a one-dimensional non-local kernel is used instead of a two-dimensional one for the dynamic problem to obtain stress fields near the crack tips. By use of the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations, in which the unknown variables are the jumps of displacements across the crack surfaces. To solve the dual integral equations, the jumps of displacements across the crack surfaces are expanded in a series of Jacobi polynomials. Unlike the classical elasticity solutions, it is found that no stress singularity is present at the crack tips. The non-local elastic solutions yield a finite hoop stress at the crack tips. The present result provides theoretical references helpful for evaluating relevant strength and preventing material failure of FGMs with initial cracks. The magnitude of the finite stress field depends on relevant parameters, such as the crack length, the distance between two parallel cracks, the parameter describing the FGMs, the frequency of the incident waves and the lattice parameter of materials. The project supported by the National Natural Science Foundation of China (90405016, 10572044) and the Specialized Research Fund for the Doctoral Program of Higher Education (20040213034). The English text was polished by Yunming Chen.     

A critical evaluation of numerical solutions to the ‘Benchmark’ surface-flaw problem

Recent precise, elastic numerical solutions to a surface flaw in a plate under remote tension and bending have been compared to arrive at a ‘best estimate’ of the stress-intensity-factor variation along the flaw border. The geometry of the semielliptical surface flaw examined had a depth to length ratio of 0.25 (a/2c=0.25) and 25- and 75-percent relative crack depths (a/t=0.25, 0.75). The analysis methods used to determine the solutions included: Schwarz alternating technique, finite-element method and boundary-integral-equation method. The derived best-estimate curve for the stress-intensity factor is believed within 3 percent of the actual value along the crack front. The best-estimate curve compared well with scarce experimental data (±10 percent). The difference between the best-estimate curve and experiment is thought largely due to differences in geometry and Poisson's ratio.     

Photoelastic determination of stress-intensity factors

The increasing number of analytical and numerical solutions for the crack-tip stress-intensity factor has greatly widened the scope of application of linear elastic fracture-mechanics technology. Experimental verification of a particular solution by elastic stress analysis is often a necessary supplement to provide the criteria for proper application to actual design problems. In this paper, it is shown that the photoelastic technique can be used to obtain rather good estimates of the stress-intensity factor for various specimen geometries and loading conditions. Treated are the following cases: wedge-opening load specimen, several notched rotating-disk configurations, and a notched pressure vessel. A sharp crack is simulated by a relatively narrow notch terminating in a root radius of 0.010 in or less. Stress distributions along the section of symmetry ahead of the notch tip are obtained using three-dimensional frozen-stress photoelasticity. The results are used to determine the stress-intensity factor, cK _I , by three methods. Two of these are based on Irwin's expressions for the elastic stress field at the tip cf a crack, and the other is a result of Neuber's hyperbolic-notch analysis. Agreement, with available analytical solutions is good.     

Interaction between an interface crack and subinterface microcracks in metal/piezoelectric bimaterials

The J-integral analysis is presented for the interaction problem between a semi-infinite interface crack and subinterface matrix microcracks in dissimilar anisotropic materials. After deriving the fundamental solutions for an interface crack subjected to different loads and the fundamental solutions for an edge dislocation beneath the interface, the interaction problem is deduced to a system of singular integral equations with the aid of a superimposing technique. The integral equations are then solved numerically and a conservation law among three values of the J-integral is presented, which are induced from the interface crack tip, the microcracks and the remote field, respectively. The conservation law not only provides a necessary condition to confirm the numerical results derived, but also reveals that the microcrack shielding effect in such materials could be considered as a redistribution of the remote J-integral. It is this redistribution that does lead to the phenomenological shielding effect.     

Variation in stress-intensity factor and back-surface displacement for surface cracks

Knowledge of the magnitude and variation in the stress-intensity factor (SIF) around the perimeter of a surface crack is essential for an accurate analysis of a flawed component. SIFs for surface flaws of various semi-elliptical geometries were analytically determined. Three-dimensional linear-elastic finite-element analysis was performed to determine the maximum SIF for bending and tension for each of 12 crack geometries which represent deep surface flaws in finite-thickness plates. Experimental verification of one of the crack geometries was performed. Interferometry techniques were used to determine the actual variation in SIF along the curve crack front due to bending. In addition to the SIF calculations, physical characteristics are noted as observed in the analytical and experimental evaluations.     

Variation of the stress intensity factor along the crack front of interacting semi-elliptical surface cracks

Summary  In this study, the interaction between two semi-elliptical co-planar surface cracks is considered when Poisson's ratio ν = 0.3. The problem is formulated as a system of singular integral equations, based on the idea of the body force method. In the numerical calculation, the unknown density of body force density is approximated by the product of a fundamental density function and a polynomial. The results show that the present method yields smooth variations of stress intensity factors along the crack front very accurately, for various geometrical conditions. When the size of crack 1 is larger than the size of crack 2, the maximum stress intensity factor appears at a certain point, β₁=177^∘, of crack 1. Along the outside of crack 1, that is at β₁=0∼90^∘, the interaction can be negligible even if the two cracks are very close. The interaction can be negligible when the two cracks are spaced in such a manner that their two closest points are separated by a distance exceeding the small crack's major diameter. The variations of stress intensity factor of a semi-elliptical crack are tabulated and charted. Received 30 August 1999; accepted for publication 22 February 2000     

The sliding interface crack with friction between elastic and rigid bodies

The paper analyzes the frictional sliding crack at the interface between a semi-infinite elastic body and a rigid one. It gives solutions in complex form for non-homogeneous loading at infinity and explicit solutions for polynomial loading at the interface. It is found that the singularities at the crack tips are different and that they are related to distinct kinematics at the crack tips. Firstly, we postulate that the geometry of the equilibrium crack with crack-tip positions b and a is determined by the conditions of square integrable stresses and continuous displacement at both crack tips. The crack geometry solution is not unique and is defined by any compatible pair (b,a) belonging to a quasi-elliptical curve. Then we prove that, for an equilibrium crack under given applied load, the “energy release rate” G_tip, defined at each crack tip by the J_ε-integral along a semi-circular path, centered at the crack tip, with vanishing radius ε, vanishes. For arbitrarily shaped paths embracing the whole crack, with end points on the unbroken zone, the J-integral is path-independent and has the significance of the rate, with respect to the crack length, of energy dissipated by friction on the crack.