Two-dimensional analysis of anisotropic crack problems using coupled meshless and fractal finite element method

This paper presents a coupling technique for integrating the element-free Galerkin method (EFGM) with fractal the finite element method (FFEM) for analyzing homogeneous, anisotropic, and two dimensional linear-elastic cracked structures subjected to mixed-mode (modes I and II) loading conditions. FFEM is adopted for discretization of domain close to the crack tip and EFGM is adopted in the rest of the domain. In the transition region interface elements are employed. The shape functions within interface elements which comprises both the element-free Galerkin and the finite element shape functions, satisfies the consistency condition thus ensuring convergence of the proposed coupled EFGM-FFEM. The proposed method combines the best features of EFGM and FFEM, in the sense that no structured mesh or special enriched basis functions are necessary and no post-processing (employing any path independent integrals) is needed to determine fracture parameters such as stress-intensity factors (SIFs) and T − stress. The numerical results based on all four orthotropic cases show that SIFs and T − stress obtained using the proposed method are in excellent agreement with the reference solutions for the structural and crack geometries considered in the present study. Also a parametric study is carried out to examine the effects of the integration order, the similarity ratio, the number of transformation terms, and the crack length to width ratio on the quality of the numerical solutions.     

Computation of the stress intensity factors of sharp notched plates by the fractal‐like finite element method

The fractal‐like finite element method (FFEM) is an accurate and efficient method to compute the stress intensity factors (SIFs) of different crack configurations. In the FFEM, the cracked/notched body is divided into singular and regular regions; both regions are modelled using conventional finite elements. A self‐similar fractal mesh of an ‘infinite’ number of conventional finite elements is used to model the singular region. The corresponding large number of local variables in the singular region around the crack tip is transformed to a small set of global co‐ordinates after performing a global transformation by using global interpolation functions. In this paper, we extend this method to analyse the singularity problems of sharp notched plates. The exact stress and displacement fields of a plate with a notch of general angle are derived for plane‐stress/strain conditions. These exact analytical solutions which are eigenfunction expansion series are used to perform the global transformation and to determine the SIFs. The use of the global interpolation functions reduces the computational cost significantly and neither post‐processing technique to extract SIFs nor special singular elements to model the singular region are needed. The numerical examples demonstrate the accuracy and efficiency of the FFEM for sharp notched problems. Copyright © 2008 John Wiley & Sons, Ltd.     

An element‐free method for in‐plane notch problems with anisotropic materials

This study developed an element‐free Galerkin method (EFGM) to simulate notched anisotropic plates containing stress singularities at the notch tip. Two‐dimensional theoretical complex displacement functions are first deduced into the moving least‐squares interpolation. The interpolation functions and their derivatives are then determined to calculate the nodal stiffness using the Galerkin method. In the numerical validation, an interface layer of the EFGM is used to combine the mesh between the traditional finite elements and the proposed singular notch EFGM. The H‐integral determined from finite element analyses with a very fine mesh is used to validate the numerical results of the proposed method. The comparisons indicate that the proposed method obtains more accurate results for the displacement, stress, and energy fields than those determined from the standard finite element method. Copyright © 2013 John Wiley & Sons, Ltd.     

A brief note on elastic T-stress for centred crack in anisotropic plate

The stress intensity factors (SIFs) and the T-stress for a planar crack with anisotropic materials are evaluated by the fractal finite element method (FFEM). The FFEM combines an exterior finite element model and a localized inner model near the crack tip. The mesh geometry of the latter is self-similar in radial layers around the tip. A higher order displacement series derived from Laurent series and Goursat functions is used to condense the large numbers of nodal displacements at the inner model near the crack tip into a small set of unknown coefficients. In this study, the variations of the SIFs and the T-stress with material properties and orientations of a crack are presented. The separation of the analytical displacement series into four fundamental cases has shown to be necessary in order to cover all the material variations and the orientations of a crack in the plate with general rectilinear anisotropic materials.     

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.     

XFEM for direct evaluation of mixed mode SIFs in homogeneous and bi‐materials

The extended finite element method (XFEM) is improved to directly evaluate mixed mode stress intensity factors (SIFs) without extra post‐processing, for homogeneous materials as well as for bimaterials. This is achieved by enriching the finite element (FE) approximation of the nodes surrounding the crack tip with not only the first term but also the higher order terms of the crack tip asymptotic field using a partition of unity method (PUM). The crack faces behind the tip(s) are modelled independently of the mesh by displacement jump functions. The additional coefficients corresponding to the enrichments at the nodes of the elements surrounding the crack tip are forced to be equal by a penalty function method, thus ensuring that the displacement approximations reduce to the actual asymptotic fields adjacent to the crack tip. The numerical results so obtained are in excellent agreement with analytical and numerical results available in the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

A quasi‐static crack growth simulation based on the singular ES‐FEM

In the edge‐based smoothed finite element method (ES‐FEM), one needs only the assumed displacement values (not the derivatives) on the boundary of the edge‐based smoothing domains to compute the stiffness matrix of the system. Adopting this important feature, a five‐node crack‐tip element is employed in this paper to produce a proper stress singularity near the crack tip based on a basic mesh of linear triangular elements that can be generated automatically for problems with complicated geometries. The singular ES‐FEM is then formulated and used to simulate the crack propagation in various settings, using a largely coarse mesh with a few layers of fine mesh near the crack tip. The results demonstrate that the singular ES‐FEM is much more accurate than X‐FEM and the existing FEM. Moreover, the excellent agreement between numerical results and the reference observations shows that the singular ES‐FEM offers an efficient and high‐quality solution for crack propagation problems. Copyright © 2011 John Wiley & Sons, Ltd.     

Direct evaluation of stress intensity factors for curved cracks using Irwin's integral and XFEM with high‐order enrichment functions

This paper presents a comprehensive study on the use of Irwin's crack closure integral for direct evaluation of mixed‐mode stress intensity factors (SIFs) in curved crack problems, within the extended finite element method. The approach employs high‐order enrichment functions derived from the standard Williams asymptotic solution, and SIFs are computed in closed form without any special post‐processing requirements. Linear triangular elements are used to discretize the domain, and the crack curvature within an element is represented explicitly. An improved quadrature scheme using high‐order isoparametric mapping together with a generalized Duffy transformation is proposed to integrate singular fields in tip elements with curved cracks. Furthermore, because the Williams asymptotic solution is derived for straight cracks, an appropriate definition of the angle in the enrichment functions is presented and discussed. This contribution is an important extension of our previous work on straight cracks and illustrates the applicability of the SIF extraction method to curved cracks. The performance of the method is studied on several circular and parabolic arc crack benchmark examples. With two layers of elements enriched in the vicinity of the crack tip, striking accuracy, even on relatively coarse meshes, is obtained, and the method converges to the reference SIFs for the circular arc crack problem with mesh refinement. Furthermore, while the popular interaction integral (a variant of the J‐integral method) requires special auxiliary fields for curved cracks and also needs cracks to be sufficiently apart from each other in multicracks systems, the proposed approach shows none of those limitations. Copyright © 2017 John Wiley & Sons, Ltd.     

Cracks emanating from a square hole in rectangular plate in tension

A numerical analysis of cracks emanating from a square hole in a rectangular plate in tension is performed using a hybrid displacement discontinuity method (a boundary element method). Detailed solutions of the stress intensity factors (SIFs) of the plane elastic crack problem are given, which can reveal the effect of geometric parameters of the cracked body on the SIFs. By comparing the calculated SIFs of the plane elastic crack problem with those of the centre crack in a rectangular plate in tension, in addition, an amplifying effect of the square hole on the SIFs is found. The numerical results reported here also prove that the boundary element method is simple, yet accurate, for calculating the SIFs of complex crack problems in finite plate.     

An extended finite element method with analytical enrichment for cohesive crack modeling

A recent approach to fracture modeling has combined the extended finite element method (XFEM) with cohesive zone models. Most studies have used simplified enrichment functions to represent the strong discontinuity but have lacked an analytical basis to represent the displacement gradients in the vicinity of the cohesive crack. In this study enrichment functions based upon an existing analytical investigation of the cohesive crack problem are proposed. These functions have the potential of representing displacement gradients in the vicinity of the cohesive crack and allow the crack to incrementally advance across each element. Key aspects of the corresponding numerical formulation and enrichment functions are discussed. A parameter study for a simple mode I model problem is presented to evaluate if quasi‐static crack propagation can be accurately followed with the proposed formulation. The effects of mesh refinement and mesh orientation are considered. Propagation of the cohesive zone tip and crack tip, time variation of the cohesive zone length, and crack profiles are examined. The analysis results indicate that the analytically based enrichment functions can accurately track the cohesive crack propagation of a mode I crack independent of mesh orientation. A mixed mode example further demonstrates the potential of the formulation. Copyright © 2008 John Wiley & Sons, Ltd.     

The extended finite element method with new crack‐tip enrichment functions for an interface crack between two dissimilar piezoelectric materials

This paper studies the static fracture problems of an interface crack in linear piezoelectric bimaterial by means of the extended finite element method (X‐FEM) with new crack‐tip enrichment functions. In the X‐FEM, crack modeling is facilitated by adding a discontinuous function and crack‐tip asymptotic functions to the classical finite element approximation within the framework of the partition of unity. In this work, the coupled effects of an elastic field and an electric field in piezoelectricity are considered. Corresponding to the two classes of singularities of the aforementioned interface crack problem, namely, ? class and κ class, two classes of crack‐tip enrichment functions are newly derived, and the former that exhibits oscillating feature at the crack tip is numerically investigated. Computation of the fracture parameter, i.e., the J‐integral, using the domain form of the contour integral, is presented. Excellent accuracy of the proposed formulation is demonstrated on benchmark interface crack problems through comparisons with analytical solutions and numerical results obtained by the classical FEM. Moreover, it is shown that the geometrical enrichment combining the mesh with local refinement is substantially better in terms of accuracy and efficiency. Copyright © 2015 John Wiley & Sons, Ltd.     

Boundary element analysis of three‐dimensional cracks in anisotropic solids

This paper presents a boundary element analysis of linear elastic fracture mechanics in three‐dimensional cracks of anisotropic solids. The method is a single‐domain based, thus it can model the solids with multiple interacting cracks or damage. In addition, the method can apply the fracture analysis in both bounded and unbounded anisotropic media and the stress intensity factors (SIFs) can be deduced directly from the boundary element solutions. The present boundary element formulation is based on a pair of boundary integral equations, namely, the displacement and traction boundary integral equations. While the former is collocated exclusively on the uncracked boundary, the latter is discretized only on one side of the crack surface. The displacement and/or traction are used as unknown variables on the uncracked boundary and the relative crack opening displacement (COD) (i.e. displacement discontinuity, or dislocation) is treated as a unknown quantity on the crack surface. This formulation possesses the advantages of both the traditional displacement boundary element method (BEM) and the displacement discontinuity (or dislocation) method, and thus eliminates the deficiency associated with the BEMs in modelling fracture behaviour of the solids. Special crack‐front elements are introduced to capture the crack‐tip behaviour. Numerical examples of stress intensity factors (SIFs) calculation are given for transversely isotropic orthotropic and anisotropic solids. For a penny‐shaped or a square‐shaped crack located in the plane of isotropy, the SIFs obtained with the present formulation are in very good agreement with existing closed‐form solutions and numerical results. For the crack not aligned with the plane of isotropy or in an anisotropic solid under remote pure tension, mixed mode fracture behavior occurs due to the material anisotropy and SIFs strongly depend on material anisotropy. Copyright © 2000 John Wiley & Sons, Ltd.     

Evaluating MTS criterion in predicting mixed-mode crack extension under different loading conditions

The maximum tangential stress (MTS) criterion is one of the most widely used criteria for predicting the direction of crack extension. The suitability of this criterion is examined under different loading conditions using the extended finite element method (XFEM). Experimental and numerical results reported in the literature are considered to evaluate the validity and accuracy of the criterion. The results demonstrate that the MTS criterion evaluated by stress intensity factors (SIFs) can accurately predict the direction of crack propagation in specimens under direct tensile loading. This criterion overestimates the angle of crack initiation in the specimen under indirect tensile loading but underestimates the angle in the specimen subjected to 3-point bending. It is concluded that the MTS criterion based on SIFs could not accurately predict the direction of the crack initiation, which could, however, be determined properly based on the stress distribution around the crack tip obtained by the XFEM numerical models.     

A unified and integrated numerical‐analytical approach for evaluation of Jk‐integrals of an interfacial crack in orthotropic and isotropic bimaterials

A new unified and integrated method for the numerical‐analytical calculation of J_k‐integrals of an in‐plane traction free interfacial crack in homogeneous orthotropic and isotropic bimaterials is presented. The numerical algorithm, based on the boundary element crack shape sensitivities, is generic and flexible. It applies to both straight and curved interfacial cracks in anisotropic and/or isotropic bimaterials. The shape functions of semidiscontinuous quadratic quarter point crack tip elements are correctly scaled to adapt the singular oscillatory near tip field of tractions. The length of crack is designated as the design variable to compute the strain energy release rate precisely. Although an analytical equation relating J₁ and stress intensity factors (SIFs) exists, a similar relation for J₂ in debonded anisotropic solids for decoupling SIFs is not available. An analytical expression was recently derived by this author for J₂ in aligned orthotropic/orthotropic bimaterials with a straight interface crack. Using this new relation and the present computed J_k values, the SIFs can be decoupled without the need for an auxiliary equation. Here, the aforementioned analytical relation is reconstructed for cubic symmetry/isotropic bimaterials and used with the present numerical algorithm. An example with known analytical SIFs is presented. The numerical and analytical magnitudes of J_k for an interface crack in orthotropic/orthotropic and cubic symmetry/isotropic bimaterials show an excellent agreement.     

Analyse the role of the non‐singular stress in brittle fracture by BEM coupled with eigen‐analysis

A coupled model resulting from the boundary element method and eigen‐analysis is proposed in this paper to analyse the stress field at crack tip. This new combine method can yield several terms of the non‐singular stress in the Williams asymptotic expansion. Then the maximum circumferential stress (MCS) criterion taken the non‐singular stress into account is introduced to predict the brittle fracture of cracked structures. Two earlier experiments are re‐examined by the present numerical method and the role of the non‐singular stress in the brittle fracture is investigated. Results show that if more terms of non‐singular stress are taken into account, the predicted crack propagation direction and the critical loading by MCS criterion are much closer to the existing experimental results, especially for dominating mode II loading conditions. Moreover, numerical results manifest that Williams series expansion can describe the stress field further from the crack tip if more non‐singular stress terms are adopted.     

Polygon scaled boundary finite elements for crack propagation modelling

An automatic crack propagation modelling technique using polygon elements is presented. A simple algorithm to generate a polygon mesh from a Delaunay triangulated mesh is implemented. The polygon element formulation is constructed from the scaled boundary finite element method (SBFEM), treating each polygon as a SBFEM subdomain and is very efficient in modelling singular stress fields in the vicinity of cracks. Stress intensity factors are computed directly from their definitions without any nodal enrichment functions. An automatic remeshing algorithm capable of handling any n‐sided polygon is developed to accommodate crack propagation. The algorithm is simple yet flexible because remeshing involves minimal changes to the global mesh and is limited to only polygons on the crack paths. The efficiency of the polygon SBFEM in computing accurate stress intensity factors is first demonstrated for a problem with a stationary crack. Four crack propagation benchmarks are then modelled to validate the developed technique and demonstrate its salient features. The predicted crack paths show good agreement with experimental observations and numerical simulations reported in the literature. Copyright © 2012 John Wiley & Sons, Ltd.     

A computational strategy for damage‐tolerant design of hollow shafts under mixed‐mode loading condition

Three‐dimensional numerical analyses, using the finite element method (FEM), have been adopted to simulate fatigue crack propagation in a hollow cylindrical specimen, under pure axial or combined axial‐torsion loading conditions. Specimens, made of Al alloys B95AT and D16T, have been experimentally tested under pure axial load and combined in‐phase constant amplitude axial and torsional loadings. The stress intensity factors (SIFs) have been calculated, according to the J‐integral approach, along the front of a part through crack, initiated in correspondence of the outer surface of a hollow cylindrical specimen. The crack path is evaluated by using the maximum energy release rate (MERR) criterion, whereas the Paris law is used to calculate crack growth rates. A numerical and experimental comparison of the results is presented, showing a good agreement in terms of crack growth rates and paths.     

An enriched meshless method for non‐linear fracture mechanics

This paper presents an enriched meshless method for fracture analysis of cracks in homogeneous, isotropic, non‐linear‐elastic, two‐dimensional solids, subject to mode‐I loading conditions. The method involves an element‐free Galerkin formulation and two new enriched basis functions (Types I and II) to capture the Hutchinson–Rice–Rosengren singularity field in non‐linear fracture mechanics. The Type I enriched basis function can be viewed as a generalized enriched basis function, which degenerates to the linear‐elastic basis function when the material hardening exponent is unity. The Type II enriched basis function entails further improvements of the Type I basis function by adding trigonometric functions. Four numerical examples are presented to illustrate the proposed method. The boundary layer analysis indicates that the crack‐tip field predicted by using the proposed basis functions matches with the theoretical solution very well in the whole region considered, whether for the near‐tip asymptotic field or for the far‐tip elastic field. Numerical analyses of standard fracture specimens by the proposed meshless method also yield accurate estimates of the J‐integral for the applied load intensities and material properties considered. Also, the crack‐mouth opening displacement evaluated by the proposed meshless method is in good agreement with finite element results. Furthermore, the meshless results show excellent agreement with the experimental measurements, indicating that the new basis functions are also capable of capturing elastic–plastic deformations at a stress concentration effectively. Copyright © 2003 John Wiley & Sons, Ltd.     

An over‐deterministic method for calculation of coefficients of crack tip asymptotic field from finite element analysis

An over‐deterministic method has been employed for calculating the stress intensity factors (SIFs) as well as the coefficients of the higher‐order terms in the Williams series expansions in cracked bodies, using the conventional finite element analysis. For a large number of nodes around the crack tip, an over‐determined set of simultaneous linear equations is obtained, and using the fundamental concepts of the least‐squares method, the coefficients of the Williams expansion can be calculated for pure mode I, pure mode II and mixed mode I/II conditions. A convergence study has been conducted to examine the effects of the number of nodes used, the number of terms in Williams expansion and the distance of the selected nodes from the crack tip, on the accuracy of the results. It is shown that the simple method presented in this paper, yields accurate results even for coarse finite element meshes or in the absence of singular elements. The accuracy of SIFs and the coefficients of higher‐order terms are validated by using the available results in the literature.     

Three‐dimensional non‐planar crack growth by a coupled extended finite element and fast marching method

A numerical technique for non‐planar three‐dimensional linear elastic crack growth simulations is proposed. This technique couples the extended finite element method (X‐FEM) and the fast marching method (FMM). In crack modeling using X‐FEM, the framework of partition of unity is used to enrich the standard finite element approximation by a discontinuous function and the two‐dimensional asymptotic crack‐tip displacement fields. The initial crack geometry is represented by two level set functions, and subsequently signed distance functions are used to maintain the location of the crack and to compute the enrichment functions that appear in the displacement approximation. Crack modeling is performed without the need to mesh the crack, and crack propagation is simulated without remeshing. Crack growth is conducted using FMM; unlike a level set formulation for interface capturing, no iterations nor any time step restrictions are imposed in the FMM. Planar and non‐planar quasi‐static crack growth simulations are presented to demonstrate the robustness and versatility of the proposed technique. Copyright © 2008 John Wiley & Sons, Ltd.