Determination of approximate point load weight functions for embedded elliptical cracks

This paper presents the application of weight function method for the calculation of stress intensity factors in embedded elliptical cracks under complex two-dimensional loading conditions. A new general mathematical form of point load weight function is proposed based on the properties of weight functions and the available weight functions for two-dimensional cracks. The existence of this general weight function form has simplified the determination of point load weight functions significantly. For an embedded elliptical crack of any aspect ratio, the unknown parameters in the general form can be determined from one reference stress intensity factor solution. This method was used to derive the weight functions for embedded elliptical cracks in an infinite body and in a semi-infinite body. The derived weight functions are then validated against available stress intensity factor solutions for several linear and non-linear stress distributions. The derived weight functions are particularly useful for the fatigue crack growth analysis of planer embedded cracks subjected to fluctuating non-linear stress fields resulting from surface treatment (shot peening), stress concentration or welding (residual stress).     

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.     

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.     

Calculation of mixed mode stress intensity factors for an elliptical subsurface crack under arbitrary normal loading

Normal loading causes mixed fracture modes in an elliptical subsurface crack because of the nonsymmetrical geometry with respect to the crack face. In this paper, mixed mode weight functions (MMWFs) for elliptical subsurface cracks in an elastic semi‐infinite space under normal loading are derived. Reference mixed mode stress intensity factors (MMSIFs), calculated by finite element analysis, under uniform normal loading are used to derive MMWFs. The cracks have aspect ratios and crack depth to crack length ratios of 0.2–1.0 and 0.05 to infinity, respectively. MMWFs are used to calculate MMSIFs for any point of the crack front under linear and nonlinear two‐dimensional (2D) loadings. So, in order to evaluate the fatigue crack growth phenomenon under complicated 2D stress distributions, MMWFs can be easily used. The comparison between the MMSIFs obtained from the MMWFs and finite element analysis indicates high accuracy.     

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.     

Stress intensity factor solutions for a cracked bolt loaded by a nut

This paper presents the calculation of stress intensity factor (K) solutions for surface cracks in the thread ground of bolts subjected to axial loading directly applied by the nut. The stress-strain computations have been done by means of the finite element method with quarter-point singular isoparametric elements along the crack front. The stress intensity factor is calculated through the stiffness derivative method, by using a virtual crack extension technique to compute the energy release rate. Two modifications are made to improve the accuracy of the results: the displacement not only of the main node, but also of the quarter-point nodes located in the normal plane and the adjacent nodes in the crack line, avoiding both the change of the singularity and the crack curving. The results show that direct loading on the thread flank by a nut increases the stress intensity factor. This effect decreases with the crack length. For the deepest circular cracks, however, nut loading relaxes the K-value, mainly at the crack surface.     

Stress intensity factors of parallel cracks in a finite width sheet

Two and three parallel cracks in a finite sheet subjected to remote tensile loading have been studied. This paper presents empirical stress intensity factor formulae for these crack configurations. The stress intensity factors used to develop these formulae were obtained from finite element analysis. For central cracks and edge cracks, the formulae were within 1 and 3% of the finite element results, respectively.     

Weight function,stress intensity factor and crack opening displacement solutions to periodic collinear edge hole cracks

Periodic collinear edge hole cracks and arbitrary small cracks emanating from collinear holes, which are two typical multiple site damages occurred in the aircraft structures, are studied by using the weigh function method. An explicit closed form weight function for periodic edge hole cracks in an infinite sheet is obtained and further used to calculate the stress intensity factor and crack opening displacement for various loading cases. Compared to finite element method, the present weight function is accurate and highly efficient. The interactions of the holes and cracks on the stress intensity factor and crack opening displacement are quantitatively determined by using the present weight function. An approximate weight function method is also proposed for arbitrary small cracks emanating from multiple collinear holes. This method is very useful for calculating the stress intensity factor for arbitrary small cracks.     

轮轨摩擦接触下钢轨多裂纹相互作用研究

利用热机耦合有限元法,建立了轮轨摩擦接触时钢轨表面多裂纹的热弹性平面应变有限元模型。数值模型中,考虑轮轨摩擦温升对轮轨材料参数的影响,通过移动载荷和热源来模拟运动车轮对钢轨的作用。分析了轮轨滑动接触时多裂纹相互作用和表面裂纹数量对钢轨疲劳裂纹扩展特性的影响。计算结果表明:与单个裂纹相比,多裂纹有降低钢轨疲劳裂纹扩展的作用;钢轨裂纹尖端应力强度因子K1和应力强度因子范围?K2均随裂纹数的增多而减小;钢轨表面裂纹数为5条时可以反映更多裂纹时的裂纹扩展特性。     

Stress intensity factors and weight functions for deep semi-elliptical surface cracks in finite-thickness plates

ABSTRACT Three-dimensional finite element analyses have been conducted to calculate the stress intensity factors for deep semi-elliptical cracks in flat plates. The stress intensity factors are presented for the deepest and surface points on semi-elliptic cracks with a/t -values of 0.9 and 0.95 and aspect ratios ( a/c ) from 0.05 to 2. Uniform, linear, parabolic or cubic stress distributions were applied to the crack face. The results for uniform and linear stress distributions were combined with corresponding results for surface cracks with a/t = 0.6 and 0.8 to derive weight functions over the range 0.05 ≤  a/c  ≤ 2.0 and 0.6 ≤  a/t  ≤ 0.95. The weight functions were then verified against finite element data for parabolic or cubic stress distributions. Excellent agreements are achieved for both the deepest and surface points. The present results complement stress intensity factors and weight functions for surface cracks in finite thickness plate developed previously.     

Calculation of approximate weight functions in fracture mechanics by FEM

A method is presented for the calculation of weight functions used in fracture mechanics to determine stress intensity factors of cracks loaded by stress gradients. The reference solution for the stress intensity factor and for the reference crack opening displacement field is computed numerically by use of finite elements. The accuracy of the method is checked by comparison with well-known solutions from the literature.     

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

Growth of three-dimensional cracks in finite-thickness plates

The stress field in a finite-thickness plate weakened by a three-dimensional crack and subjected to tension in a direction perpendicular to the crack plane is studied. The cases of an embedded elliptical, semi-elliptical and quarter-elliptical surface crack are considered. The stress analysis takes place by a finite element computer program which uses twenty-node isoparametric and fifteen-node enriched elements. The stress intensity factor variations along the periphery of the elliptical cracks are given for various plate thicknesses. The results of the stress analysis are used in conjunction with the strain energy density theory to study the growth characteristics of the cracks. The history of non-self-similar crack growth from initiation to final instability through the intermediate stage of stable growth is analyzed. The increments of crack growth from each point of its front are determined on the basis that the critical element in the direction of crack growth absorbs a critical amount of strain energy density. Crack growth becomes unstable when the last increment from the critical point of the crack front takes a limiting value. Results for the crack growth characteristics are given for the three types of cracks considered and various plate thicknesses.     

Weight function calculation for surface cracks using the line-spring model

A numerical method for calculating weight functions for surface cracks in plates and shells is proposed. Thick-shell finite elements are used to create the discrete model of a body with a through-wall flaw. Line-spring elements transform the through-wall flaw into a surface crack. A quadratic line-spring element is presented. Weight functions for some semielliptical surface cracks in a plate have been calculated. The weight functions obtained may be used for computing stress intensity factors related to two-dimensional stress fields at the crack surface.     

Fatigue crack growth analyses of offshore structural tubular joints

This study investigated various aspects of a fatigue crack growth analysis, ranging from the stress intensity factor solutions to the simulation of a fatigue crack coalescence process of a tubular joint weld toe surface flaw. Fracture mechanics fatigue crack growth analyses for offshore structural tubular joints are not simple, because of the difficulty to calculate the stress intensity factors due to their geometric complexity. The fully mixed-mode stress intensity factors of nine weld toe surface cracks of an X-shaped tubular joint under tension loading were calculated by detailed three-dimensional finite element analyses. Using these stress intensity factor solutions, a fatigue crack growth study was performed for the X-joint until (the crack surface length grew to two times the tube thickness. Through this study, the crack shape change during the fatigue crack propagation was investigated in detail. Fatigue life calculations were also performed for a range of crack geometries using the stress intensity factor solutions of the nine flaws. These calculations indicate that the natural fatigue crack growing path for a crack is its quickest growing path. The study demonstrated that detailed fracture mechanics fatigue analyses of tubular joints can be practical using the finite element method.     

Stress intensity factor solutions for estimation of fatigue lives of laser welds in lap-shear specimens

Fatigue behavior of laser welds in lap-shear specimens of high strength low alloy (HSLA) steel is investigated based on experimental observations and two fatigue life estimation models. Fatigue experiments of laser welded lap-shear specimens are first reviewed. Analytical stress intensity factor solutions for laser welded lap-shear specimens based on the beam bending theory are derived and compared with the analytical solutions for two semi-infinite solids with connection. Finite element analyses of laser welded lap-shear specimens with different weld widths were also conducted to obtain the stress intensity factor solutions. Approximate closed-form stress intensity factor solutions based on the results of the finite element analyses in combination with the analytical solutions based on the beam bending theory and Westergaard stress function for a full range of the normalized weld widths are developed for future engineering applications. Next, finite element analyses for laser welded lap-shear specimens with three weld widths were conducted to obtain the local stress intensity factor solutions for kinked cracks as functions of the kink length. The computational results indicate that the kinked cracks are under dominant mode I loading conditions and the normalized local stress intensity factor solutions can be used in combination with the global stress intensity factor solutions to estimate fatigue lives of laser welds with the weld width as small as the sheet thickness. The global stress intensity factor solutions and the local stress intensity factor solutions for vanishing and finite kinked cracks are then adopted in a fatigue crack growth model to estimate the fatigue lives of the laser welds. Also, a structural stress model based on the beam bending theory is adopted to estimate the fatigue lives of the welds. The fatigue life estimations based on the kinked fatigue crack growth model agree well with the experimental results whereas the fatigue life estimations based on the structural stress model agree with the experimental results under larger load ranges but are higher than the experimental results under smaller load ranges.     

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

A finite element comparison between short and long cracks within a plastic zone due to a notch

This paper considers the finite element characteristics of a crack growing within an area affected by a small notch. It is a continuation of an investigation leading from the development of a finite element program which considers elastic-viscoplastic constitutive relations to the actual features present in modelling a crack growth within a specimen. A final paper actually compares given parameters developed within this paper with experimental findings. The material considered is IN 718 at a temperature of 1200°F. A compact tension with a blunt notch is the specimen eventually tested. The findings relate to the most effective technique for allowing the crack to grow within a specimen. A trace is carried out on stress and strain functions as the crack grows with a consideration of the closure phenomenon. Parameters such as the J range, COD, maximum stress and strain range, stress intensity, and effective stress intensity range are shown to be potential candidates for crack growth as short and long cracks are compared.     

CRACK-LINE AND EDGE GREEN'S FUNCTIONS FOR STRESS INTENSITY FACTORS OF INCLINED EDGE CRACKS

Abstract— Fretting loads on the surfaces of structural components can cause accelerated growth of short cracks. The rate of growth will depend on the combined stress intensity factor resulting from both remote and local loading. Many stress intensity factor solutions are available for remote loading, but solutions for arbitrary fretting loads are not readily accessible. In this paper accurate crack-line Green's functions are obtained from a boundary element analysis and then used to develop the Green's functions for loads on the edge of a half-plane containing a slant crack at various angles to the edge. These latter Green's functions can be used to obtain stress intensity factors for arbitrary stresses (normal or shear) on the edge of the half-plane without further stress analysis; simple integration procedures are all that is required.     

A method for stress intensity factor calculation of infinite plate containing multiple hole-edge cracks

Multiple site damage is the occurrence of small fatigue cracks at several sites within aging aircraft structures. Focusing on this typical structure, an analytical method for calculating the stress intensity factor of an infinite plate containing multiple hole-edge cracks was introduced in this paper. The properties of complex variable functions are used to evaluate the stress function. The approximate superposition method is applied to solve stress intensity factor problems on multiple holes. The equivalent crack is introduced to modify the method. Some numerical examples of an infinite plate containing two hole-edge cracks are examined by the method. By comparing the analytical and finite element analysis results it was realized that the analytical results are accurate and reliable. This modified analytical method is easier to apply than some traditional analytical methods and can provide stress intensity factor solutions for an infinite plate containing a random distribution of multiple hole-edge cracks.