Extended finite element method and fast marching method for three-dimensional fatigue crack propagation

A numerical technique for planar three-dimensional fatigue crack growth simulations is proposed. The new technique couples the extended finite element method (X-FEM) to the fast marching method (FMM). In the X-FEM, a discontinuous function and the two-dimensional asymptotic crack-tip displacement fields are added to the finite element approximation to account for the crack using the notion of partition of unity. This enables the domain to be modeled by finite elements with no explicit meshing of the crack surfaces. The initial crack geometry is represented by level set functions, and subsequently signed distance functions are used to compute the enrichment functions that appear in the displacement-based finite element approximation. The FMM in conjunction with the Paris crack growth law is used to advance the crack front. Stress intensity factors for planar three-dimensional cracks are computed, and fatigue crack growth simulations for planar cracks are presented. Good agreement between the numerical results and theory is realized.     

Extended finite element method for three‐dimensional crack modelling

An extended finite element method (X‐FEM) for three‐dimensional crack modelling is described. A discontinuous function and the two‐dimensional asymptotic crack‐tip displacement fields are added to the finite element approximation to account for the crack using the notion of partition of unity. This enables the domain to be modelled by finite elements with no explicit meshing of the crack surfaces. Computational geometry issues associated with the representation of the crack and the enrichment of the finite element approximation are discussed. Stress intensity factors (SIFs) for planar three‐dimensional cracks are presented, which are found to be in good agreement with benchmark solutions. Copyright © 2000 John Wiley & Sons, Ltd.     

Non‐planar 3D crack growth by the extended finite element and level sets—Part II: Level set update

We present a level set method for treating the growth of non‐planar three‐dimensional cracks.The crack is defined by two almost‐orthogonal level sets (signed distance functions). One of them describes the crack as a two‐dimensional surface in a three‐dimensional space, and the second is used to describe the one‐dimensional crack front, which is the intersection of the two level sets. A Hamilton–Jacobi equation is used to update the level sets. A velocity extension is developed that preserves the old crack surface and can accurately generate the growing surface. The technique is coupled with the extended finite element method which approximates the displacement field with a discontinuous partition of unity. This displacement field is constructed directly in terms of the level sets, so the discretization by finite elements requires no explicit representation of the crack surface. Numerical experiments show the robustness of the method, both in accuracy and in treating cracks with significant changes in topology. Copyright © 2002 John Wiley & Sons, Ltd.     

hp‐Generalized FEM and crack surface representation for non‐planar 3‐D cracks

A high‐order generalized finite element method (GFEM) for non‐planar three‐dimensional crack surfaces is presented. Discontinuous p‐hierarchical enrichment functions are applied to strongly graded tetrahedral meshes automatically created around crack fronts. The GFEM is able to model a crack arbitrarily located within a finite element (FE) mesh and thus the proposed method allows fully automated fracture analysis using an existing FE discretization without cracks. We also propose a crack surface representation that is independent of the underlying GFEM discretization and controlled only by the physics of the problem. The representation preserves continuity of the crack surface while being able to represent non‐planar, non‐smooth, crack surfaces inside of elements of any size. The proposed representation also provides support for the implementation of accurate, robust, and computationally efficient numerical integration of the weak form over elements cut by the crack surface. Numerical simulations using the proposed GFEM show high convergence rates of extracted stress intensity factors along non‐planar curved crack fronts and the robustness of the method. Copyright © 2008 John Wiley & Sons, Ltd.     

Fatigue crack propagation of multiple coplanar cracks with the coupled extended finite element/fast marching method

A numerical technique for modeling fatigue crack propagation of multiple coplanar cracks is presented. The proposed method couples the extended finite element method (X-FEM) [Int. J. Numer. Meth. Engng. 48 (11) (2000) 1549] to the fast marching method (FMM) [Level Set Methods & Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science, Cambridge University Press, Cambridge, UK, 1999]. The entire crack geometry, including one or more cracks, is represented by a single signed distance (level set) function. Merging of distinct cracks is handled naturally by the FMM with no collision detection or mesh reconstruction required. The FMM in conjunction with the Paris crack growth law is used to advance the crack front. In the X-FEM, a discontinuous function and the two-dimensional asymptotic crack-tip displacement fields are added to the finite element approximation to account for the crack using the notion of partition of unity [Comput. Meth. Appl. Mech. Engng. 139 (1996) 289]. This enables the domain to be modeled by a single fixed finite element mesh with no explicit meshing of the crack surfaces. In an earlier study [Engng. Fract. Mech. 70 (1) (2003) 29], the methodology, algorithm, and implementation for three-dimensional crack propagation of single cracks was introduced. In this paper, simulations for multiple planar cracks are presented, with crack merging and fatigue growth carried out without any user-intervention or remeshing.     

New development in extended finite element modeling of large elasto‐plastic deformations

This paper presents new achievements in the extended finite element modeling of large elasto‐plastic deformation in solid problems. The computational technique is presented based on the extended finite element method (X‐FEM) coupled with the Lagrangian formulation in order to model arbitrary interfaces in large deformations. In X‐FEM, the material interfaces are represented independently of element boundaries, and the process is accomplished by partitioning the domain with some triangular sub‐elements whose Gauss points are used for integration of the domain of elements. The large elasto‐plastic deformation formulation is employed within the X‐FEM framework to simulate the non‐linear behavior of materials. The interface between two bodies is modeled by using the X‐FEM technique and applying the Heaviside‐ and level‐set‐based enrichment functions. Finally, several numerical examples are analyzed, including arbitrary material interfaces, to demonstrate the efficiency of the X‐FEM technique in large plasticity deformations. Copyright © 2008 John Wiley & Sons, Ltd.     

Non‐planar mixed‐mode growth of initially straight‐fronted surface cracks,in cylindrical bars under tension,torsion and bending,using the symmetric Galerkin boundary element method‐finite element method alternating method

In this paper, the stress intensity factor (SIF) variations along an arbitrarily developing crack front, the non‐planar fatigue‐crack growth patterns, and the fatigue life of a round bar with an initially straight‐fronted surface crack, are studied by employing the 3D symmetric Galerkin boundary element method‐finite element method (SGBEM‐FEM) alternating method. Different loading cases, involving tension, bending and torsion of the bar, with different initial crack depths and different stress ratios in fatigue, are considered. By using the SGBEM‐FEM alternating method, the SIF variations along the evolving crack front are computed; the fatigue growth rates and directions of the non‐planar growths of the crack surface are predicted; the evolving fatigue‐crack growth patterns are simulated, and thus, the fatigue life estimations of the cracked round bar are made. The accuracy and reliability of the SGBEM‐FEM alternating method are verified by comparing the presently computed results to the empirical solutions of SIFs, as well as experimental data of fatigue crack growth, available in the open literature. It is shown that the current approach gives very accurate solutions of SIFs and simulations of fatigue crack growth during the entire crack propagation, with very little computational burden and human–labour cost. The characteristics of fatigue growth patterns of initially simple‐shaped cracks in the cylindrical bar under different Modes I, III and mixed‐mode types of loads are also discussed in detail.     

Fracture in magnetoelectroelastic materials using the extended finite element method

Static fracture analyses in two‐dimensional linear magnetoelectroelastic (MEE) solids is studied by means of the extended finite element method (X‐FEM). In the X‐FEM, crack modeling is facilitated by adding a discontinuous function and the crack‐tip asymptotic functions to the standard finite element approximation using the framework of partition of unity. In this study, media possessing fully coupled piezoelectric, piezomagnetic and magnetoelectric effects are considered. New enrichment functions for cracks in transversely isotropic MEE materials are derived, and the computation of fracture parameters using the domain form of the contour interaction integral is presented. The convergence rates in energy for topological and geometric enrichments are studied. Excellent accuracy of the proposed formulation is demonstrated on benchmark crack problems through comparisons with both analytical solutions and numerical results obtained by the dual boundary element method. Copyright © 2011 John Wiley & Sons, Ltd.     

A marching cubes based failure surface propagation concept for three‐dimensional finite elements with non‐planar embedded strong discontinuities of higher‐order kinematics

This paper presents an advanced failure surface propagation concept based on the marching cubes algorithm initially proposed in the field of computer graphics and applies it to the embedded finite element method. When modeling three‐dimensional (3D) solids at failure, the propagation of the failure surface representing a crack or shear band should not exhibit a strong sensitivity to the details of the finite element discretization. This results in the need for a propagation of the discrete failure zone through the individual finite elements, which is possible for finite elements with embedded strong discontinuities. Whereas for two‐dimensional calculations the failure zone propagation location is easily predicted by the maximal principal stress direction, more advanced strategies are needed to achieve a smooth failure surface in 3D simulations. An example for such method is the global tracking algorithm, which predicts the crack path by a scalar level set function computed on the basis of the solution of a simplified heat conduction like problem. Its prediction may though lead to various scenarios on how the failure surface may propagate through the individual finite elements. In particular, for a hexahedral eight‐node finite element, 256 such cases exist. To capture all those possibilities, the marching cubes algorithm is combined with the global tracking algorithm and the finite elements with embedded strong discontinuities in this work. In addition, because many of the possible cases result in non‐planar failure surfaces within a single finite element and because the local quantities used to describe the kinematics of the embedded strong discontinuities are physically meaningful in a strict sense only for planar failure surfaces, a remedy for such scenarios is proposed. Various 3D failure propagation simulations outline the performance of the proposed concept. Copyright © 2013 John Wiley & Sons, Ltd.     

Modeling crack discontinuities without element‐partitioning in the extended finite element method

In this paper, we model crack discontinuities in two‐dimensional linear elastic continua using the extended finite element method without the need to partition an enriched element into a collection of triangles or quadrilaterals. For crack modeling in the extended finite element, the standard finite element approximation is enriched with a discontinuous function and the near‐tip crack functions. Each element that is fully cut by the crack is decomposed into two simple (convex or nonconvex) polygons, whereas the element that contains the crack tip is treated as a nonconvex polygon. On using Euler's homogeneous function theorem and Stokes's theorem to numerically integrate homogeneous functions on convex and nonconvex polygons, the exact contributions to the stiffness matrix from discontinuous enriched basis functions are computed. For contributions to the stiffness matrix from weakly singular integrals (because of enrichment with asymptotic crack‐tip functions), we only require a one‐dimensional quadrature rule along the edges of a polygon. Hence, neither element‐partitioning on either side of the crack discontinuity nor use of any cubature rule within an enriched element are needed. Structured finite element meshes consisting of rectangular elements, as well as unstructured triangular meshes, are used. We demonstrate the flexibility of the approach and its excellent accuracy in stress intensity factor computations for two‐dimensional crack problems. Copyright © 2016 John Wiley & Sons, Ltd.     

A meshfree‐enriched finite element method for compressible and near‐incompressible elasticity

In this paper, a two‐dimensional displacement‐based meshfree‐enriched FEM (ME‐FEM) is presented for the linear analysis of compressible and near‐incompressible planar elasticity. The ME‐FEM element is established by injecting a first‐order convex meshfree approximation into a low‐order finite element with an additional node. The convex meshfree approximation is constructed using the generalized meshfree approximation method and it possesses the Kronecker‐delta property on the element boundaries. The gradient matrix of ME‐FEM element satisfies the integration constraint for nodal integration and the resultant ME‐FEM formulation is shown to pass the constant stress test for the compressible media. The ME‐FEM interpolation is an element‐wise meshfree interpolation and is proven to be discrete divergence‐free in the incompressible limit. To prevent possible pressure oscillation in the near‐incompressible problems, an area‐weighted strain smoothing scheme incorporated with the divergence‐free ME‐FEM interpolation is introduced to provide the smoothing on strains and pressure. With this smoothed strain field, the discrete equations are derived based on a modified Hu–Washizu variational principle. Several numerical examples are presented to demonstrate the effectiveness of the proposed method for the compressible and near‐incompressible problems. Copyright © 2012 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.     

Improvements of explicit crack surface representation and update within the generalized finite element method with application to three‐dimensional crack coalescence

This paper presents improvements to three‐dimensional crack propagation simulation capabilities of the generalized finite element method. In particular, it presents new update algorithms suitable for explicit crack surface representations and simulations in which the initial crack surfaces grow significantly in size (one order of magnitude or more). These simulations pose problems in regard to robust crack surface/front representation throughout the propagation analysis. The proposed techniques are appropriate for propagation of highly non‐convex crack fronts and simulations involving significantly different crack front speeds. Furthermore, the algorithms are able to handle computational difficulties arising from the coalescence of non‐planar crack surfaces and their interactions with domain boundaries. An approach based on moving least squares approximations is developed to handle highly non‐convex crack fronts after crack surface coalescence. Several numerical examples are provided, which illustrate the robustness and capabilities of the proposed approaches and some of its potential engineering applications. Copyright © 2013 John Wiley & Sons, Ltd.     

Bubble‐enhanced smoothed finite element formulation: a variational multi‐scale approach for volume‐constrained problems in two‐dimensional linear elasticity

This paper presents a bubble‐enhanced smoothed finite element formulation for the analysis of volume‐constrained problems in two‐dimensional linear elasticity. The new formulation is derived based on the variational multi‐scale approach in which unequal order displacement‐pressure pairs are used for the mixed finite element approximation and hierarchical bubble function is selected for the fine‐scale displacement approximation. An area‐weighted averaging scheme is employed for the two‐scale smoothed strain calculation under the framework of edge‐based smoothed FEM. The smoothed fine‐scale solution is shown to naturally contain the stress field jump of the smoothed coarse‐scale solution across the boundary of edge‐based smoothing domain and thus provides the possibility to stabilize the global solution for volume‐constrained problems. A global monolithic solution strategy is employed, and the fine‐scale solution is solved without the consideration of approximating the strong form of the fine‐scale equation. Several numerical examples are analyzed to demonstrate the accuracy of the present formulation. Copyright © 2014 John Wiley & Sons, Ltd.     

几何非线性扩展有限元法及其断裂力学应用

该文将扩展有限元方法应用到几何非线性及断裂力学问题中,并研制开发了扩展有限元Fortran程序。扩展有限元法其计算网格与不连续面相互独立,因此模拟移动的不连续面时无需对网格进行重新剖分。该文推导了几何非线性扩展有限元法的公式,在常规有限元位移模式中,基于单位分解的思想加进一个阶跃函数和二维渐近裂尖位移场,反映裂纹处位移的不连续性,并用2个水平集函数表示裂纹;采用拉格朗日描述方程建立了有限变形几何非线性扩展有限元方程;采用多点位移外推法计算裂纹应力强度因子并通过最小二乘法拟合得到更精确的结果。最后给出的大变形算例表明该文提出的几何非线性的断裂力学扩展有限元方法和相应的计算机程序是合理可行的,而且对于含裂纹及裂纹扩展的问题,扩展有限元法优于传统的有限元法。     

Crack face contact in X‐FEM using a segment‐to‐segment approach

In this study, we propose a segment‐to‐segment contact formulation (mortar‐based) that uses Lagrange's multipliers to establish the contact between crack faces when modeled with the extended finite element method (X‐FEM) in 2D problems. It is shown that, in general, inaccuracies arise when the contact is formulated following a point‐to‐point approach. This is due to the non‐linear character of the X‐FEM interpolation along the crack faces that leads to crack face interpenetration. However, the segment‐to‐segment approach optimizes the fulfilment of the contact constraints along the whole crack segment, and in practice the contact is modeled precisely. Convergence studies for mesh sequences have been performed, showing the advantages of the proposed methodology. Copyright © 2009 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.     

The Discontinuity‐Enriched Finite Element Method

We introduce a new methodology for modeling problems with both weak and strong discontinuities independently of the finite element discretization. At variance with the eXtended/Generalized Finite Element Method (X/GFEM), the new method, named the Discontinuity‐Enriched Finite Element Method (DE‐FEM), adds enriched degrees of freedom only to nodes created at the intersection between a discontinuity and edges of elements in the mesh. Although general, the method is demonstrated in the context of fracture mechanics, and its versatility is illustrated with a set of traction‐free and cohesive crack examples. We show that DE‐FEM recovers the same rate of convergence as the standard FEM with matching meshes, and we also compare the new approach to X/GFEM.     

Stress‐based finite element analysis of plane plasticity problems

A stress‐based model of the finite element method is evolved for two‐dimensional quasi‐static plasticity problems. The self‐equilibrating fields of stresses are constructed by means of the Airy stress function, which is approximated by three types of elements: the Bogner–Fox–Schmit rectangle, the Hsieh–Clough–Tocher triangle and its reduced variant. Traction boundary conditions are imposed by the use of the Lagrange multiplier method which gives the possibility of calculation of displacements for boundary points. The concept of multi‐point‐constraints elements is applied in order to facilitate the application of this technique. The iterative algorithm, analogous to the closest‐point‐projection method commonly used in the displacement‐based finite element model, is proposed for solving non‐linear equations for each load increment. Two numerical examples with stress‐ and displacement‐controlled load are considered. The results are compared with those obtained by the displacement model of FEM. Bounds for limit loads are obtained. Copyright © 1999 John Wiley & Sons, Ltd.