Developing new enrichment functions for crack simulation in orthotropic media by the extended finite element method

New enrichment functions are proposed for crack modelling in orthotropic media using the extended finite element method (XFEM). In this method, Heaviside and near‐tip functions are utilized in the framework of the partition of unity method for modelling discontinuities in the classical finite element method. In this procedure, by using meshless based ideas, elements containing a crack are not required to conform to crack edges. Therefore, mesh generation is directly performed ignoring the existence of any crack while the method remains capable of extending the crack without any remeshing requirement. Furthermore, the type of elements around the crack‐tip remains the same as other parts of the finite element model and the number of nodes and consequently degrees of freedom are reduced considerably in comparison to the classical finite element method. Mixed‐mode stress intensity factors (SIFs) are evaluated to determine the fracture properties of domain and to compare the proposed approach with other available methods. In this paper, the interaction integral (M‐integral) is adopted, which is considered as one of the most accurate numerical methods for calculating stress intensity factors. Copyright © 2006 John Wiley & Sons, Ltd.     

Experimental and numerical investigation of the effect of temperature on mixed‐mode fracture behaviour of AM60 Mg alloy

The aim of present paper is to experimentally investigate mixed‐mode fracture behaviour of AM60 Mg alloy at low and elevated temperatures. For this purpose, mode I, 45° mixed‐mode, and mode II tests were conducted using a modified version of Arcan device at three different temperatures. An elastic‐plastic finite element model was used to extract necessary geometric parameters. Crack resistance curves (J‐R) and critical J‐integral of the material were extracted. The results indicated that, for all loading modes, maximum critical J‐integral value was observed at ambient temperature and decreased by either increasing or decreasing the temperature. It was observed that effect of temperature on fracture behaviour is much larger at temperatures above 0°C rather than sub‐zero temperatures. By changing the loading angle to go from mode I to mode II, a decreasing trend was observed in the values of critical fracture parameters at all temperatures. Finally, the surfaces were examined using scanning electron microscopy (SEM).     

Stress-intensity factors for circular hollow section V-joints with a rack-plate chord

This study investigates the fatigue crack‐driving force, measured by the linear‐elastic stress‐intensity factors (SIFs), for a surface crack at the root of the welds in a thick‐walled, circular hollow section (CHS) V‐shape joint, typically installed in modern offshore jack‐up platforms. The primary (chord) member of the V‐joint consists of two half CHSs welded to both sides of a thick rack plate, while the secondary (brace) member adopts thick‐walled CHSs. The surface‐breaking crack considered in this study locates at the interface between the rack plate and the weld metal joining the half CHS, and represents an initial flaw introduced by lack of penetration in the welding procedure. The finite‐element model incorporates a very detailed, local crack‐front mesh in a global continuous mesh through a mesh‐tying procedure, which ensures displacement continuity between the independent master surface and the dependent slave surface. A simple plate model verifies the mesh‐tying procedure in computing the linear‐elastic stress‐intensity factors for two remote loading conditions. The computation of the stress‐intensity factors employs a linear‐elastic interaction integral approach. The comparison of the computed SIF values with a previous experimental measurement for a CHS T‐joint verifies the accuracy and feasibility of the interaction integral approach in computing SIF values for surface cracks in welded tubular connections. Subsequent numerical analysis on the gapped V‐joints examines the mixed‐mode SIF values for different loading conditions and includes an array of practical joint geometric parameters and crack sizes. The nondimensional mode I stress‐intensity factors generally increases with the following variations in the joint geometric parameters: an increase in the chord radius to the wall thickness ratio (γ=d₀/2t₀) , an increase in the brace diameter to the chord diameter ratio (β=d₁/d₀) , a decrease in the crack depth ratio (a/t) or an increase in the crack length c. The current study identifies a practical group of V‐joints that requires detailed treatment in the fatigue assessment procedure. These V‐joints adopt a large β ratio and demonstrate high mode‐mixity angles [ψ= tan^?1(K_II/K_I)] with correspondingly high mode I and mode II stress‐intensity factors.     

Failure mode of laser welds in lap‐shear specimens of high strength low alloy (HSLA) steel sheets

In this paper, the failure mode of laser welds in lap‐shear specimens of non‐galvanized SAE J2340 300Y high strength low alloy steel sheets under quasi‐static loading conditions is examined based on experimental observations and finite element analyses. Laser welded lap‐shear specimens with reduced cross sections were made. Optical micrographs of the cross sections of the welds in the specimens before and after tests are examined to understand the microstructure and failure mode of the welds. Micro‐hardness tests were also conducted to provide an assessment of the mechanical properties in the base metal, heat‐affected and fusion zones. The micrographs indicate that the weld failure appears to be initiated from the base metal near the boundary of the base metal and the heat‐affected zone at a distance away from the pre‐existing crack tip, and the specimens fail due to the necking/shear of the lower left load carrying sheets. Finite element analyses based on non‐homogenous multi‐zone material models were conducted to model the ductile necking/shear failure and to obtain the J integral solutions for the pre‐existing cracks. The results of the finite element analyses are used to explain the ductile failure initiation sites and the necking/shear of the lower left load carrying sheets. The J integral solutions obtained from the finite element analyses based on the 3‐zone finite element model indicate that the J integral for the pre‐existing cracks at the failure loads are low compared to the fracture toughness and the specimens should fail in a plastic collapse or necking/shear mode. The effects of the sheet thickness on the failure mode were then investigated for laser welds with a fixed ratio of the weld width to the thickness. For the given non‐homogenous material model, the J integral solutions appear to be scaled by the sheet thickness. With consideration of the plastic collapse failure mode and fracture initiation failure mode, a critical thickness can be obtained for the transition of the plastic collapse or necking/shear failure mode to the fracture initiation failure mode. Finally, the failure load is expressed as a function of the sheet thickness according to the governing equations based on the two failure modes. The results demonstrate that the failure mode of welds of thin sheets depends on the sheet thickness, ductility of the base metal and fracture toughness of the heat‐affected zone. Therefore, failure criteria based on either the plastic collapse failure mode or the fracture initiation failure mode should be used cautiously for welds of thin sheets.     

Constraint‐corrected J–R curve based on three‐dimensional finite element analyses

Three‐dimensional (3D) finite element analyses are carried out on single‐edge bend [SE(B)] specimens for which the J‐integral resistance curves (J–R curves) have been experimentally determined to develop the constraint‐corrected J–R curves for the X80 grade pipe steel. The constraint parameters considered in this study include Q_HRR, Q_SSY, Q_{SSY_m}, Q_LM, Q_BM1, Q_BM2, A₂, h and T_z. The constraint‐corrected J–R curves were developed on the basis of the constraint parameters obtained from finite element analysis and experimentally determined J–R curves associated with deeply cracked and medium‐cracked SE(B) specimens and validated against shallow‐cracked SE(B) specimens. The analysis results indicate that all the constraint parameters considered in this study except Q_HRR, Q_SSY, Q_{SSY_m} and Q_LM lead to reasonably accurate constraint‐corrected J–R curves if the crack extensions are relatively small (≤0.7 mm). For larger crack extensions (≤1.5 mm), the Q_BM1‐based constraint‐corrected J–R curve leads to the most accurate predictions of J among all the constraint parameters considered.     

Mode III interlaminar fracture of carbon/epoxy laminates using a four-point bending plate test

A new four-point bending plate (4PBP) test was used for characterising the mode III interlaminar fracture of carbon/epoxy laminates. The specimen has a cross-ply lay-up and two edge delaminations whose propagation becomes visible at the edges. Although the test setup is very simple, determination of the mode III critical strain energy release rate G_IIIc requires finite element analyses (FEA). The virtual crack closure technique with an assumed initiation region was first proposed for computing G_IIIc. This scheme was subsequently validated by crack growth simulations with a cohesive zone model. The results showed an average G_IIIc = 1550 J/m², which is significantly higher than the G_IIIc = 850–1100 J/m² and G_IIc = 800 J/m² measured in previous studies.     

Fatigue crack propagation behavior of RC beams strengthened with CFRP under cyclic bending loads

A fatigue crack propagation equation of reinforced concrete (RC) beams strengthened with a new type carbon fiber reinforced polymer was proposed in this paper on the basis of experimental and numerical methods. Fatigue crack propagation tests were performed to obtain the crack propagation rate of the strengthened RC beams. Digital image correlation method was used to capture the fatigue crack pattern. Finite element model of RC beam strengthened with carbon fiber reinforced polymer was established to determinate J‐integral of a main crack considering material nonlinearities and degradation of material properties under cyclic loading. Paris law with a parameter of J‐integral was developed on the basis of the fatigue tests and finite element analysis. This law was preliminarily verified, which can be applied for prediction of fatigue lives of the strengthened RC beams.     

3D J‐integral evaluation using the computation of line and surface integrals

In the analysis of fracture mechanics of structures using three‐dimensional (3D) J‐integral, an integral evaluation of line and surface is required. However, because surface integral evaluation requires the calculation of the second derivative of displacement field and commercial finite element codes cannot calculate it, then this portion of the integral is neglected in some research. In this paper, a method for computing 3D J‐integral is presented using finite element analysis. In the analysis, the second derivative evaluation of displacement field is employed. The method is implemented in calculating the J‐integral of some 3D cracks and results are compared to well‐known reference values. The results show that the method is reliable and is suitable for applications in engineering. The portion of 3D J‐integral, namely the surface integral value is investigated and it is shown that neglecting this portion can introduce considerable error in the final results.     

Robust and direct evaluation of J2 in linear elastic fracture mechanics with the X‐FEM

The aim of the present paper is to study the accuracy and the robustness of the evaluation of J_k‐integrals in linear elastic fracture mechanics using the extended finite element method (X‐FEM) approach. X‐FEM is a numerical method based on the partition of unity framework that allows the representation of discontinuity surfaces such as cracks, material inclusions or holes without meshing them explicitly. The main focus in this contribution is to compare various approaches for the numerical evaluation of the J₂‐integral. These approaches have been proposed in the context of both classical and enriched finite elements. However, their convergence and the robustness have not yet been studied, which are the goals of this contribution. It is shown that the approaches that were used previously within the enriched finite element context do not converge numerically and that this convergence can be recovered with an improved strategy that is proposed in this paper. Copyright © 2008 John Wiley & Sons, Ltd.     

Application of the small punch test to determine the fracture toughness of metallic materials

A new methodology to determine the elasto‐plastic fracture toughness, J_Ic, by means of notched small punch tests (SPT) samples is reported. Standard SPT samples were used after being longitudinally notched machined from the centre of one side of the sample to the centre of the opposite side, producing a notch depth‐to‐thickness ratio a/t= 0.4. The onset of crack initiation was experimentally determined directly from the experimental load‐displacement plot of the test and also with the aid of scanning electron microscope observations performed on different samples, with tests being interrupted at different percentages of the maximum registered load. The test was also modelled using finite element analysis and the J‐integral was evaluated as a contour integral in ABAQUS. The obtained results were compared with the critical J values of the steel determined using standard tests (J–R curves) and the differences found were duly justified.     

A direct analytical method to extract mixed‐mode components of strain energy release rates from Irwin's integral using extended finite element method

A new analytical approach, within the extended finite element framework, is proposed to compute mixed‐mode components of strain energy release rates directly from Irwin's integral. Crack tip enrichment functions in extended FEM allow for evaluation of integral quantities in closed form (for some crack configurations studied) and therefore resulting in a simple and accurate method. Several benchmark examples on pure and mixed‐mode problems are studied. In particular, we analyze the effects of high‐order enrichments, mesh refinement, and the integration limits of Irwin's integral. The results indicate that high‐order enrichment functions have significant effect on the convergence, in particular when the integral limits are finite. When the integral limits tend to zero, simpler strain energy release rate expressions are obtained, and high‐order terms vanish. Nonetheless, these terms contribute indirectly via coefficients of first‐order terms. The numerical results show that high accuracy can be achieved with high‐order enrichment terms and mesh refinement. However, the effect of the integral limits remains an open question, with finite integration intervals chosen as h ∕ 2 tending to give more accurate results. Copyright © 2013 John Wiley & Sons, Ltd.     

Studies of elastic and elastic–plastic J integral for mixed mode cracked plate under biaxial loading

Analyses of I–II mixed mode central cracked plate by finite element method are performed in this paper, and some different phenomena are found. First for I–II mixed mode crack, the distribution of J integral along crack tip thickness depends on biaxiality factors because of the existence of vertex (corner) singularity, which is unlike that for mode I or mode II crack. Then J integrals at middle layer keep constant for any cracked plates with different inclined angles β when the biaxiality ratio is equal to 1 or ?1, which implies that the inclined angle or the extent of I–II mixed mode has no effect on the J integral for positive or negative equal axial loading conditions. And the decreasing trend of J integral with the inclined angle β for biaxiality ratio λ being between?1 and 1 is just opposite with that for biaxiality ratio λ being larger than 1 and smaller than ?1. Finally, proposed h₁ (a/W, n, λ, β) of cracked plate with different inclined angles under different biaxial loading are calculated.     

An ablation problem in a finite medium with time‐dependent heat flux

Abstract

The objective of this work is to develop an approximate analytical solution for the transient ablation problem in a finite medium. The medium is subjected to time‐dependent boundary heat flux, i. e., q₀=at^P , and with this refined integral technique, the complicated nonlinear problem is reduced to an initial value problem, which is then solved by the Runge‐Kutta method. Results are more accurate than with the classical heat balance integral method and also indicate that the dependence of the solution on the assumed temperature profile is much weaker than is the case with the classical integral method.     

The effect of T-stress on plane strain dynamic crack growth in elastic–plastic materials

In this paper, the effects of T‐stress on steady, dynamic crack growth in an elastic–plastic material are examined using a modified boundary layer formulation. The analyses are carried out under mode I, plane strain conditions by employing a special finite element procedure based on moving crack tip coordinates. The material is assumed to obey the J₂ flow theory of plasticity with isotropic power law hardening. The results show that the crack opening profile as well as the opening stress at a finite distance from the tip are strongly affected by the magnitude and sign of the T‐stress at any given crack speed. Further, it is found that the fracture toughness predicted by the analyses enhances significantly with negative T‐stress for both ductile and cleavage mode of crack growth.     

On the use of an anti‐symmetric four‐point bend specimen for mode II fracture experiments

The edge‐cracked beam specimen subjected to anti‐symmetric four‐point bend (ASFPB) loading has been conventionally used in the past for investigating the pure mode II fracture experiments in many engineering materials. However, it is shown through finite element analysis that the ASFPB specimen sometimes fails to produce pure mode II conditions. For anti‐symmetric loads applied close to the crack line, there are considerable effects from K_I and T‐stress in the ASFPB specimen. Pure mode II is provided only when the applied loads are sufficiently far from the crack plane.     

Centroidal Voronoi tessellation‐based finite element superconvergence

In this article, a finding on finite element superconvergence is reported. The Laplacian operator with Dirichlet boundary condition is considered. The linear finite element solutions have an O(h^2+α)(α≈0.5)‐superconvergence in l₂ norm at nodes on an almost equilateral triangular mesh generated based on centroidal Voronoi tessellation, for an arbitrary 2D bounded domain. Extensive numerical examples are presented to demonstrate the superconvergence property. Copyright © 2008 John Wiley & Sons, Ltd.     

Hierarchic finite element bases on unstructured tetrahedral meshes

The problem of constructing hierarchic bases for finite element discretization of the spaces H¹, H ( curl ), H ( div ) and L₂ on tetrahedral elements is addressed. A simple and efficient approach to ensuring conformity of the approximations across element interfaces is described. Hierarchic bases of arbitrary polynomial order are presented. It is shown how these may be used to construct finite element approximations of arbitrary, non‐uniform, local order approximation on unstructured meshes of curvilinear tetrahedral elements. Copyright © 2003 John Wiley & Sons, Ltd.     

Finite elements for materials with strain gradient effects

A finite element implementation is reported of the Fleck–Hutchinson phenomenological strain gradient theory. This theory fits within the Toupin–Mindlin framework and deals with first‐order strain gradients and the associated work‐conjugate higher‐order stresses. In conventional displacement‐based approaches, the interpolation of displacement requires C¹‐continuity in order to ensure convergence of the finite element procedure for higher‐order theories. Mixed‐type finite elements are developed herein for the Fleck–Hutchinson theory; these elements use standard C⁰‐continuous shape functions and can achieve the same convergence as C¹ elements. These C⁰ elements use displacements and displacement gradients as nodal degrees of freedom. Kinematic constraints between displacement gradients are enforced via the Lagrange multiplier method. The elements developed all pass a patch test. The resulting finite element scheme is used to solve some representative linear elastic boundary value problems and the comparative accuracy of various types of element is evaluated. Copyright © 1999 John Wiley & Sons, Ltd.     

Engineering C*‐integral and COD estimates for nonidealized circumferential through‐wall cracked pipes at elevated temperature

In this study, creep fracture mechanics parameters, C^*‐integral and crack opening displacement (COD) rate, are estimated for a nonidealized circumferential through‐wall crack (TWC) in pipes. The GE/EPRI and enhanced reference stress (ERS) methods are employed. As for creep condition, the Norton and RCC‐MRx creep models are considered for secondary and primary‐secondary creep strain, respectively. The bending moment, axial tension, and internal pressure are applied to a pipe with a nonidealized circumferential TWC, as individual loads. Three‐dimensional elastic‐creep finite element (FE) analyses are performed, and the predictions from the GE/EPRI and ERS methods are compared with FE results. For the Norton creep model, both methods show good agreement with the FE results. For the RCC‐MRx creep model, only the ERS method can be used, and it provides accurate predictions comparing with FE results. Based on the comparison results, the use of the present engineering C^*‐integral and COD estimation methods for nonidealized circumferential TWC is validated.     

An adaptive finite element method for second‐order plate theory

We present a discontinuous finite element method for the Kirchhoff plate model with membrane stresses. The method is based on P₂‐approximations on simplices for the out‐of‐plane deformations, using C⁰‐continuous approximations. We derive a posteriori error estimates for linear functionals of the error and give some numerical examples. Copyright © 2009 John Wiley & Sons, Ltd.