Weight functions for the determination of stress intensity factor and T‐stress for semi‐elliptical cracks in finite thickness plate

This paper presents the application of weight function method for the calculation of stress intensity factors (K) and T‐stress for surface semi‐elliptical crack in finite thickness plates subjected to arbitrary two‐dimensional stress fields. New general mathematical forms of point load weight functions for K and T have been formulated by taking advantage of the knowledge of a few specific weight functions for two‐dimensional planar cracks available in the literature and certain properties of weight function in general. The existence of the generalised forms of the weight functions simplifies the determination of specific weight functions for specific crack configurations. The determination of a specific weight function is reduced to the determination of the parameters of the generalised weight function expression. These unknown parameters can be determined from reference stress intensity factor and T‐stress solutions. This method is used to derive the weight functions for both K and T for semi‐elliptical surface cracks in finite thickness plates, covering a wide range of crack aspect ratio (a/c) and relative depth (a/t) at any point along the crack front. The derived weight functions are then validated against stress intensity factor and T‐stress solutions for several linear and nonlinear two‐dimensional stress distributions. These derived weight functions are particularly useful for the development of two‐parameter fracture and fatigue models for surface cracks subjected to fluctuating nonlinear stress fields, such as these resulting from surface treatment (shot peening), stress concentration or welding (residual stress).     

Weight functions and Stress Intensity Magnification factors for elliptical and semi-elliptical cracks under variable normal stress – part I

Surface cracks under peak stresses are investigated. The calculational procedure is based on the general form of the weight function for an elliptical crack embedded in an infinite solid. Two points on the contour of the ellipse are investigated. The superposition method is used for transfer from the embedded crack to surface crack configurations. Weight functions for both points have been found with the crack aspect ratio a/c as parameter. For the point at the end of the minor axis all weight functions are describable by one equation only (Heuman's lambda function). For various a/c ratios of the surface crack under different stress distributions the stress intensity magnification factors are given.     

Weight functions and Stress Intensity Magnification factors for elliptical and semi-elliptical cracks under variable normal stress – Part II

Surface cracks under peak stresses are investigated. The calculational procedure is based on the general form of the weight function for an elliptical crack embedded in an infinite solid. Two points on the contour of the ellipse are investigated. A new correction procedure for transfer from the embedded crack to surface crack configurations is presented, which is valid for all a/t-values. Weight functions for both points have been found with the crack aspect ratio a/c as parameter. For the point at the end of the minor axis all weight functions for embedded cracks are describable by one equation only (using Heuman's lambda function). For various a/c-ratios of the surface crack under different stress distributions the stress intensity magnification factors are given.     

Weight functions and stress intensity factors for quarter-elliptical corner cracks in fastener holes

Approximate weight functions for a quarter‐elliptical crack in a fastener hole were derived from a general weight function form and two reference stress intensity factors. Closed‐form expressions were obtained for the coefficients of the weight functions. The derived weight functions were validated against numerical data by comparison of stress intensity factors calculated for several nonlinear stress fields. Good agreements were achieved. These derived weight functions are valid for the geometric range of 0.5 ≤a/c≤ 1.5 and 0 ≤a/t≤ 0.8 and R/t= 0.5; and are given in forms suitable for computer numerical integration. The weight functions appear to be particularly suitable for fatigue crack growth prediction of corner cracks in fastener holes and fracture analysis of such cracks in complex stress fields.     

STRESS INTENSITY FACTORS AND WEIGHT FUNCTIONS FOR SEMI-ELLIPTICAL SURFACE CRACKS IN FINITE-THICKNESS PLATES UNDER TWO-DIMENSIONAL STRESS DISTRIBUTION

Abstract— A Fourier series approach is proposed to calculate stress intensity factors using weight functions for semi-elliptical surface cracks in flat plates subjected to two-dimensional stress distributions. The weight functions were derived from reference stress intensity factors obtained by three-dimensional finite element analyses. The close form weight functions derived are suitable for the calculation of stress intensity factors for semi-elliptical surface cracks in flat plates under two-dimensional stress distributions with the crack aspect ratio in the range of 0.1 ≤ a/c ≤ 1 and relative depth in the range of 0 ≤ a/t ≤ 0.8. Solutions were verified using several two-dimensional non-linear stress distributions; the maximum difference being 6%.     

Approximate weight function method for elliptical arc through cracks in mechanical joints

Mechanical joints such as bolted, riveted or pinned joints are widely used to join the constituent parts of structural components. Reliable stress intensity factor analysis of arbitrary cracks in mechanical joints is required for the safety evaluation or fracture mechanics design. It has been reported that cracks in mechanical joints usually nucleate as the corner crack and grow as the elliptical arc through crack. The weight function method is a useful technique to calculate the stress intensity factor using the appropriate weight function for a cracked body and the stress field of an uncracked body. In this paper, the weight function method for the two surface points of elliptical arc through cracks in mechanical joints is developed to analyze the mixed-mode stress intensity factors. Unknown coefficients included in the weight function are determined using the reference stress intensity factors obtained from finite element analysis.     

Mode-I Elliptical Crack in an Infinite Elastic Body. Part 1. Crack-Face Displacement for the Polynomial Law of Loading

For a mode-I embedded elliptical crack in an infinite elastic body, we propose a modified procedure for the calculation of the weight function and stress intensity factor based on our earlier development. Analytical and numerical values of the stress intensity factor along the crack front have been obtained for different cases of the polynomial law of loading. We propose an approach to the determination of the crack-face displacements from the stress intensity factor values, in which the Rice energy-balance equation, Dyson's theorem, and the theory of crack translation in a nonuniform stress field are used. An expression of closed form for the elliptical crack-face displacement for a polynomial law of loading of any degree has been derived, which can be employed in solving three-dimensional problems of the elasticity theory for cracked bodies.     

Weight Function Method for Stress Intensity Factor of Internal Surface Semi-elliptical Crack in Elliptical Holes

The hatches for inspecting are usually designed with elliptical holes in airplane structures, so computation of the stress intensity factor of three dimensional crack at elliptical holes is pivotal for damage tolerance analysis of these structures. In this paper, weight function is derived for a two dimensional through cracks at elliptical holes by applying a compounding method. Stress intensity factor formulas for an internal surface semi-elliptical crack in elliptical holes are obtained wing the three dimensional weight function method. Stress intensity factors for an internal surface semi-elliptical crack in elliptical holes under remote tension are computed. At the same time, research on how radius of curvature for elliptical holes affect stress intensity factors was conducted. Stress intensity factors decrease when radius of curvature increases. Some results and conclusions which are of practical value are given.     

Analytical solution for embedded elliptical cracks,and finite element alternating method for elliptical surface cracks,subjected to arbitrary loadings

The complete solution for an embedded elliptical crack in an infinite solid and subjected to arbitrary tractions on the crack surface is rederived from Vijayakumar and Atluri's general solution procedure. The general procedure for evaluating the necessary elliptic integrals in the generalized solution for elliptical crack is also derived in this paper. The generalized solution is employed in the Schwartz alternating technique in conjunction with the finite element method. This finite element-alternating method gives an inexpensive way to evaluate accurate stress intensity factors for embedded or elliptical cracks in engineering structural components.     

WEIGHT FUNCTIONS AND STRESS INTENSITY FACTORS FOR SEMI-ELLIPTICAL CRACKS IN T-PLATE WELDED JOINTS

Weight functions were derived for the deepest point and surface point of a semi-elliptical surface crack in T-plate joints with weld angles between 0 and 45°. These weight functions were derived from reference stress intensity factor solutions obtained from three-dimensional finite element calculations, and verified using stress intensity factors for different non-linear stress fields and for far-field tension and bending cases. The differences between the weight function predictions and the finite element data were less than 10%. They are suitable for semi-elliptical surface cracks with aspect ratios in the range 0.05 ≤ a/c ≤ 1, together with relative depths 0 ≤ a/t ≤ 0.6 and weld angles 0 ≤ φ ≤ 45°.     

Interacting dissimilar semi-elliptical surface flaws under tension and bending

Stress intensity factors for two dissimilar interacting semi-elliptical coplanar surface flaws (cracks) in a semi-infinite elastic body are obtained under overall tension and bending. First the basic equations for a general planar crack normal to the free surface are established, using the method of equivalent eigen- or transformation strains (the body force method). Then the results are specialized for application to elliptical cracks. Numerical values are obtained for various configurations and crack shapes. Results are compared with those of two-dimensional collinear cracks. Finally, an approximate procedure for estimating the stress intensity factors for a general three-dimensional crack is suggested.     

A generalised solution for the crack surface displacement of mode I two-dimensional part-elliptical crack

A generalised approximate crack surface displacement solution for the two-dimensional part-elliptical mode I crack was developed. This solution includes the surface crack, corner crack and embedded crack, which is subjected to the arbitrary crack surface pressure. The crack surface displacement is derived from stress intensity factor solution and corresponding crack surface pressure distribution. Comparisons of the solution with accurate solutions showed that rather high accuracy has been achieved with the developed solution for various surface, embedded and corner crack problems. This solution can be used to derive three-dimensional weight functions as long as the stress intensity factor and the corresponding crack surface pressure are available for arbitrary mode I problems.     

Computation of shear mode weight functions for circular and elliptical cracks

Three-dimensional shear mode fundamental fields in finite bodies with mixed boundary conditions are analyzed by a special finite element method for circular and elliptical cracks. A procedure for determining the Fourier coefficients of the stress intensity factor for circular cracks is presented. A special series is proposed to represent the computed crack face weight functions for elliptical cracks.     

The effect of an elliptical inclusion on a crack

The interaction between an elliptical inclusion and a crack is analyzed by body force method. The investigated stress field is simulated by superposing the fundamental solutions for a point force applied at a point in an infinite plate containing an elliptical inclusion. Based on numerical results, effects of the inclusion shape on the crack tip stress intensity factor are discussed. It is found that for small cracks emanating from a stress-higher point on the inclusion interface the stress intensity factors are mainly determined by the stresses, occurring at the crack starting point before the crack initiation, and the inclusion root radius, besides the crack length. However, for the cracks occurring in a stress-lower region around the inclusion, it is difficult to characterize the effect of the inclusion geometry on the stress intensity factors of small cracks by the inclusion root radius alone. This revised version was published online in July 2006 with corrections to the Cover Date.     

A weight function method for mixed modes hole‐edge cracks

A weight function approach is proposed to calculate the stress intensity factor and crack opening displacement for cracks emanating from a circular hole in an infinite sheet subjected to mixed modes load. The weight function for a pure mode II hole‐edge crack is given in this paper. The stress intensity factors for a mixed modes hole‐edge crack are obtained by using the present mode II weight function and existing mode I Green (weight) function for a hole‐edge crack. Without complex derivation, the weight functions for a single hole‐edge crack and a centre crack in infinite sheets are used to study 2 unequal‐length hole‐edge cracks. The stress intensity factor and crack opening displacement obtained from the present weight function method are compared well with available results from literature and finite element analysis. Compared with the alternative methods, the present weight function approach is simple, accurate, efficient, and versatile in calculating the stress intensity factor and crack opening displacement.     

Variations of stress intensity factor of a semi-elliptical surface crack subjected to mixed mode loading

Maximum stress intensity factors of a surface crack usually appear at the deepest point of the crack, or a certain point along crack front near the free surface depending on the aspect ratio of the crack. However, generally it has been difficult to obtain smooth distributions of stress intensity factors along the crack front accurately due to the effect of corner point singularity. It is known that the stress singularity at a corner point where the front of 3 D cracks intersect free surface is depend on Poisson's ratio and different from the one of ordinary crack. In this paper, a singular integral equation method is applied to calculate the stress intensity factor along crack front of a 3-D semi-elliptical surface crack in a semi-infinite body under mixed mode loading. The body force method is used to formulate the problem as a system of singular integral equations with singularities of the form r ⁻³ using the stress field induced by a force doublet in a semi-infinite body as fundamental solution. In the numerical calculation, unknown body force densities are approximated by using fundamental density functions and polynomials. The results show that the present method yields smooth variations of mixed modes stress intensity factors along the crack front accurately. Distributions of stress intensity factors are indicated in tables and figures with varying the elliptical shape and Poisson's ratio.     

含界面裂纹圆形夹杂的十二次对称二维准晶热应力分析

针对点热源作用下,无限大十二次对称二维准晶基体和圆形弹性夹杂界面之间含多条裂纹的问题进行了研究。基于复变函数分区全纯理论、留数定理、广义 Liouville 定理、Riemann-Schwarz 解析延拓定理及复应力函数奇性主部分析方法,获得了集中热源作用于准晶基体内任意一点时,准晶基体和圆形弹性夹杂内外温度场、声子场热应力的一般复势解。由此获得了含一条界面裂纹和两条界面裂纹时温度场以及声子场热应力的封闭形式解答,将所得结果与已有结果进行了对比,验证了该方法的有效性。最后通过数值算例分析了夹杂半径、点热源强度及裂纹角度对热应力和裂纹尖端热应力强度因子的影响规律。结果表明：随着热源强度的增大,裂纹尖端的声子场热应力也逐渐增大;随着裂纹角度的增大,裂纹尖端的声子场热应力强度因子变大;随着半径的增大,热应力强度因子的变化趋势越来越明显,并且取得的峰值越高,即裂纹角度和夹杂半径的增加,促进了裂纹的扩展。这些结论为准晶材料的结构设计和使用提供了科学依据。     

Calculation of stress intensity factors for functionally graded materials by using the weight functions derived by the virtual crack extension technique

The weight function method provides a powerful approach for calculating the stress intensity factors for a homogeneous cracked body subjected to mechanical loadings. In this paper, the basic equations of weight function method for mode I and mixed mode crack problems in a two-dimensional functionally graded crack system are derived based on the Betti’s reciprocal theorem. The weight functions derived by the virtual crack extension technique are further used to calculate the stress intensity factors of functionally graded materials (FGMs). The practicability and accuracy of this proposed method has been confirmed by the comparison with theoretical or numerical solutions available in the literatures. On account that the repeated extractions of the distributions of normal stress and shear stress in the uncracked component along the prospective crack line under different loadings can be avoided using the method presented in this paper, this method can be potentially used for optimal design for FGMs under multiple-load cases.     

Stress intensity factors for an interface crack along an elliptical inclusion

The problem of a crack along the interface of an elliptical elastic inclusion embedded in an infinite plate subjected to uniform stresses at infinity is analyzed by the body force method. The crack tip stress intensity factors are calculated for various inclusion geometries and material combinations. Based on numerical results, the effect of the inclusion geometry on the stress intensity factors is investigated. It is found that for small interface cracks the stress intensity factors are mainly determined by the stresses, occurring at the crack center point before the crack initiation, and interface curvature radius alone.     

Weight functions and strip yield solution for two equal-length collinear cracks in an infinite sheet

The problem of two equal-length collinear cracks in an infinite sheet is treated using the weight function method. Exact weight functions for the inner and outer crack tips are derived based on the crack opening displacement solution for a reference load case. These weight functions are used to calculate stress intensity factors for different load cases, plastic zone sizes and crack tip opening displacements of the strip yield model. The approach is validated by the perfect agreement between the present strip yield model solutions and Collins and Cartwright’s analytical results based on the direct complex stress function formulation.