Mixed-mode stress intensity factors for arbitrarily oriented cracks in cylindrical shells under axial pull and torsion

Results of photoelastic investigations conducted on cylindrical tubes (made of Araldite material) containing cracks oriented at 0°, 30°, 45°, 60° and 90° to the axis of the tube and subjected to axial and torsional loads are reported. The stress-intensity factors (SIFs) were determined by analysing the crack-tip stress fields. Smith and Smith's method [Engng Fracture Mech.4, 357–366 (1972)] and a new method developed by the authors by modifying Rakesh et al.'s method [Proc. 26th Congress of ISTAM, India (1981)] were employed to evaluate the mixed-mode SIFs.

Effect of interaction of macrocracks on the stress intensity factor in a beam

Beams with a central edge crack, as well as other noncentral vertical and inclined edge cracks distributed symmetrically, subjected to three-point as well as four-point bending, are analysed using the finite element technique. Values of stress intensity factor K1 at the central crack tip for a crack-to-beam depth ratio Image equal to 0.5, are calculated for various cracked-beam configurations, using the compliance calibration technique as well as method of strain energy release rate. These are compared with the value of K1 for the case of a beam with central edge crack alone. Results of the present parametric study are used to specify the range of values pertaining to basic parameters such as crack-to-beam depth ratios, geometry and position with respect to central edge crack, of other macrocracks for which interaction exists. Accordingly, the macrocracks are classified as either interacting or noninteracting types. Hence for noninteracting types of cracks, analytical expressions available for the determination of K1 in the case of beam with a central edge crack alone, are applicable.

Stress intensity factor determination of radially cracked circular rings subjected to tension using photoelastic technique

A parametric study was carried out to determine the Stress Intensity Factor (SIF) in a cracked circular ring by using the photoelastic technique. The stress intensity factors for mode I deformation were determined by subjecting the specimens to the tensile loading from inner boundary and through the holes. The results of Non-Dimensional Stress Intensity Factor (NDSIF) variation with non-dimensional crack length for both methods of loading are compared with each other and with published results.

New a posteriori error estimator and adaptive mesh refinement strategy for 2-D crack problems

A posteriori error estimation and adaptive refinement technique for fracture analysis of 2-D/3-D crack problems is the state-of-the-art. The objective of the present paper is to propose a new a posteriori error estimator based on strain energy release rate (SERR) or stress intensity factor (SIF) at the crack tip region and to use this along with the stress based error estimator available in the literature for the region away from the crack tip. The proposed a posteriori error estimator is called the K-S error estimator. Further, an adaptive mesh refinement (h-) strategy which can be used with K-S error estimator has been proposed for fracture analysis of 2-D crack problems. The performance of the proposed a posteriori error estimator and the h-adaptive refinement strategy have been demonstrated by employing the 4-noded, 8-noded and 9-noded plane stress finite elements. The proposed error estimator together with the h-adaptive refinement strategy will facilitate automation of fracture analysis process to provide reliable solutions.

A two-dimensional model for crack closure effect in plates under bending

A two-dimensional model is proposed for taking into account the establishment of contact on the compression side of crack faces in plates under bending. An approximate but simple method is developed for evaluating reduction of stress intensity factor due to such ‘crack closure’. Analysis is first carried out permitting interference of the crack faces. Contact forces are then introduced on the crack faces and their magnitudes determined from the consideration that the interference is just eliminated. The method is based partly on finite element analysis and partly on a continuum analysis using Irwin's solution for point loads on the crack line.

Effect of prestrain on fatigue crack propagation in mild steel containing non-metallic inclusions

The effect of tensile prestrain on fatigue crack propagation behaviour of commercial mild steel with significant amount of stringer inclusions has been studied. In prestrained materials the usual stable stage II crack growth region is preceded by a phase wherein a retardation in crack growth rate occurs. No such behaviour is observed in annealed material. The amount of retardation is found to increase with increase in prestrain. A mechanism for the observed retardation in crack growth rate is also presented.

Fatigue crack propagation model for plain concrete - An analogy with population growth

In this work, an analytical model is proposed for fatigue crack propagation in plain concrete based on population growth exponential law and in conjunction with principles of dimensional analysis and self-similarity. This model takes into account parameters such as loading history, fracture toughness, crack length, loading ratio and structural size. The predicted results are compared with experimental crack growth data for constant and variable amplitude loading and are found to capture the size effect apart from showing a good agreement. Using this model, a sensitivity analysis is carried out to study the effect of various parameters that influence fatigue failure. (C) 2010 Elsevier Ltd. All rights reserved.

A general procedure for modified crack closure integral in 3D problems with cracks

In linear elastic fracture mechanics (LEFM), Irwin's crack closure integral (CCI) is one of the signficant concepts for the estimation of strain energy release rates (SERR) G, in individual as well as mixed-mode configurations. For effective utilization of this concept in conjunction with the finite element method (FEM), Rybicki and Kanninen [Engng Fracture Mech. 9, 931 938 (1977)] have proposed simple and direct estimations of the CCI in terms of nodal forces and displacements in the elements forming the crack tip from a single finite element analysis instead of the conventional two configuration analyses. These modified CCI (MCCI) expressions are basically element dependent. A systematic derivation of these expressions using element stress and displacement distributions is required. In the present work, a general procedure is given for the derivation of MCCI expressions in 3D problems with cracks. Further, a concept of sub-area integration is proposed which facilitates evaluation of SERR at a large number of points along the crack front without refining the finite element mesh. Numerical data are presented for two standard problems, a thick centre-cracked tension specimen and a semi-elliptical surface crack in a thick slab. Estimates for the stress intensity factor based on MCCI expressions corresponding to eight-noded brick elements are obtained and compared with available results in the literature.

Cracks emanating from pin-loaded lugs

Pin-loaded lugs were analysed in the presence of cracks emanating from circular holes. The analysis presents a unified treatment of interference, push or clearance fit pins. Both metallic (isotropic) and composite (orthotropic) plates were dealt with. The finite element model used special singular six-noded quadrilateral elements at the crack tip. The non-linear load contact behaviour at the pin-hole interface was dealt with by an inverse technique. A modified crack closure integral (MCCI) technique was used to evaluate the strain energy release rates (SERRs) and stress intensity factors (SIFs) at the crack tips. Numerical results are presented showing the non-linear variation of SIF with applied stress, and the influence of the amount of interference or clearance and the interfacial friction on SIF.

Modified Crack Closure Integral (MCCI) For 3-D Problems Using 20-Noded Brick Elements

The Modified Crack Closure Integral (MCCI) technique based on Irwin's crack closure integral concept is very effective for estimation of strain energy release rates G in individual as well as mixed-mode configurations in linear elastic fracture mechanics problems. In a finite element approach, MCCI can be evaluated in the post-processing stage in terms of nodal forces and displacements near the crack tip. The MCCI expressions are however, element dependent and require a systematic derivation using stress and displacement distributions in the crack tip elements. Earlier a general procedure was proposed by the present authors for the derivation of MCCI expressions for 3-dimensional (3-d) crack problems modelled with 8-noded brick elements. A concept of sub-area integration was proposed to estimate strain energy release rates at a large number of points along the crack front. In the present paper a similar procedure is adopted for the derivation of MCCI expressions for 3-d cracks modelled with 20-noded brick elements. Numerical results are presented for centre crack tension and edge crack shear specimens in thick slabs, showing a comparison between present results and those available in the literature.

Modified crack closure integral (MCCI) for 3-d problems using 20 noded brick elememnts

The Modified Crack Closure Integral (MCCI) technique based on Irwin's crack closure integral concept is very effective for estimation of strain energy release rates G in individual as well as mixed-mode configurations in linear elastic fracture mechanics problems. In a finite element approach, MCCI can be evaluated in the post-processing stage in terms of nodal forces and displacements near the crack tip. The MCCI expressions are however, element dependent and require a systematic derivation using stress and displacement distributions in the crack tip elements. Earlier a general procedure was proposed by the present authors for the derivation of MCCI expressions for 3-dimensional (3-d) crack problems modelled with 8-noded brick elements. A concept of sub-area integration was proposed to estimate strain energy release rates at a large number of points along the crack front. In the present paper a similar procedure is adopted for the derivation of MCCI expressions for 3-d cracks modelled with 20-noded brick elements. Numerical results are presented for centre crack tension and edge crack shear specimens in thick slabs, showing a comparison between present results and those available in the literature.

On the performance of strain smoothing for quadratic and enriched finite element approximations (XFEM/GFEM/PUFEM)

By using the strain smoothing technique proposed by Chen et al. (Comput. Mech. 2000; 25: 137-156) for meshless methods in the context of the finite element method (FEM), Liu et al. (Comput. Mech. 2007; 39(6): 859-877) developed the Smoothed FEM (SFEM). Although the SFEM is not yet well understood mathematically, numerical experiments point to potentially useful features of this particularly simple modification of the FEM. To date, the SFEM has only been investigated for bilinear and Wachspress approximations and is limited to linear reproducing conditions. The goal of this paper is to extend the strain smoothing to higher order elements and to investigate numerically in which condition strain smoothing is beneficial to accuracy and convergence of enriched finite element approximations. We focus on three widely used enrichment schemes, namely: (a) weak discontinuities; (b) strong discontinuities; (c) near-tip linear elastic fracture mechanics functions. The main conclusion is that strain smoothing in enriched approximation is only beneficial when the enrichment functions are polynomial (cases (a) and (b)), but that non-polynomial enrichment of type (c) lead to inferior methods compared to the standard enriched FEM (e.g. XFEM). Copyright (C) 2011 John Wiley & Sons, Ltd.

Role of elasto-plastic analysis under cyclic loading in fatigue crack growth studies

Linear Elastic Fracture Mechanics (LEFM) has been widely used in the past for fatigue crack growth studies, but this is acceptable only in situations which are within small scale yielding (SSY). In many practical structural components, conditions of SSY could be violated and one has to look for fracture criteria based on elasto-plastic analysis. Crack closure phenomenon, one of the most striking discoveries based on inelastic deformations during crack growth, has significant effect on fatigue crack growth rate. Numerical simulation of this phenomenon is computationally intensive and involved but has been successfully implemented. Stress intensity factors and strain energy release rates lose their meaning, J-integral (or its incremental) values are applicable only in specific situations, whereas alternate path independent integrals have been proposed in the literature for use with elasto-plastic fracture mechanics (EPFM) based criteria. This paper presents certain salient features of two independent finite element (numerical) studies of relevance to fatigue crack growth, where elasto-plastic analysis becomes significant. These problems can only be handled in the current day computational environment, and would have been only a dream just a few years ago.

The fatigue crack growth rate in L-72 Al-alloy plate specimens with cold worked holes

A fatigue crack growth rate study has been carried out on L-72 aluminium alloy plate specimens with and without cold worked holes. The cold worked specimens showed significantly increased fatigue life compared to unworked specimens. Computer software is developed to evaluate the stress intensity factor for non-uniform stress distributions using Green's function approach. The exponents for the Paris equation in the stable crack growth region for cold worked and unworked specimens are 1.26 and 3.15 respectively. The reduction in exponent value indicates the retardation in crack growth rate. An SEM study indicates more plastic deformation at the edge of the hole for unworked samples as compared to the worked samples during the crack initiation period.

Modelling of a patch repair to a thin cracked sheet

Integrity enhancement of damaged or design deficient structures through repairs is attracting considerable engineering attention. Bonded composite patch repairs to cracked metallic sheets offer various advantages over riveted doubler type, particularly for airframe applications. This paper first reviews the R&D activity in the area of structural repairs. It then approaches the problem of a composite patch repair to a cracked aluminium sheet with different finite element modelling strategies and compares their outcome. The efficient finite element modelling approach thus established is used to study the effect of patch material, patch size, patch symmetry and adhesive thickness on repair performance as the crack grows in the repair configuration. (C) 1999 Elsevier Science Ltd. All rights reserved.