A two-dimensional augmented finite element for dynamic crack initiation and propagation

In this study, an implicit formulation of a 2D finite element based on a recently developed augmented finite element method is proposed for stable and efficient simulation of dynamic fracture in elastic solids. The 2D A-FE ensures smooth transition from a continuous state to a discontinuous state with an arbitrary intra-element cohesive crack, without the need of additional degree of freedoms (DoFs). Internal nodal DoFs are introduced for sub-domain integration and cohesive stress integration and they are then condensed at elemental level by a consistency-check based algorithm. The numerical performance of the proposed A-FE has been assessed through simulations of several benchmark dynamic fracture problems and in all cases the numerical results are in good agreement with the respective experimental results and other simulation results in literature. It has further been demonstrated that, (i) the dynamic A-FE is rather insensitive to mesh sizes and mesh structures; (ii) with similar solution accuracy it allows for the use of time steps 1–2 orders of magnitude larger than those used in other similar studies; and (iii) the implicit nature of the proposed A-FE allows for the use of a Courant number as large as 3.0–3.5 while maintaining solution stability.     

Incremental‐secant modulus iteration scheme and stress recovery for simulating cracking process in quasi‐brittle materials using XFEM

In this paper, an incremental‐secant modulus iteration scheme using the extended/generalized finite element method (XFEM) is proposed for the simulation of cracking process in quasi‐brittle materials described by cohesive crack models whose softening law is composed of linear segments. The leading term of the displacement asymptotic field at the tip of a cohesive crack (which ensures a displacement discontinuity normal to the cohesive crack face) is used as the enrichment function in the XFEM. The opening component of the same field is also used as the initial guess opening profile of a newly extended cohesive segment in the simulation of cohesive crack propagation. A statically admissible stress recovery (SAR) technique is extended to cohesive cracks with special treatment of non‐homogeneous boundary tractions. The application of locally normalized co‐ordinates to eliminate possible ill‐conditioning of SAR, and the influence of different weight functions on SAR are also studied. Several mode I cracking problems in quasi‐brittle materials with linear and bilinear softening laws are analysed to demonstrate the usefulness of the proposed scheme, as well as the characteristics of global responses and local fields obtained numerically by the XFEM. Copyright © 2006 John Wiley & Sons, Ltd.     

Thermoelectroelastic analysis of cracks in piezoelectric half-plane by BEM

The Green's function and the boundary element method for analysing fracture behaviour of cracks in piezoelectric half-plane are presented in this paper. By combining Stroh formalism and the concept of perturbation, a general thermoelectroelastic solution for half-plane solid subjected to point heat source and/or temperature discontinuity has been derived. Using the proposed solution and the potential variational principle, a boundary element model (BEM) for 2-D half-plane solid with multiple cracks has been developed and used to calculate the stress intensity factors of the multiple crack problem. The method is available for multiple crack problems in both finite and infinite solids. Numerical results for a two-crack system are presented and compared with those from finite element method (FEM).     

Combined extended and superimposed finite element method for cracks

A combination of the extended finite element method (XFEM) and the mesh superposition method (s‐version FEM) for modelling of stationary and growing cracks is presented. The near‐tip field is modelled by superimposed quarter point elements on an overlaid mesh and the rest of the discontinuity is implicitly described by a step function on partition of unity. The two displacement fields are matched through a transition region. The method can robustly deal with stationary crack and crack growth. It simplifies the numerical integration of the weak form in the Galerkin method as compared to the s‐version FEM. Numerical experiments are provided to demonstrate the effectiveness and robustness of the proposed method. Copyright © 2004 John Wiley & Sons, Ltd.     

Mesh‐independent matrix cracking and delamination modeling in laminated composites

The initiation and evolution of transverse matrix cracks and delaminations are predicted within a mesh‐independent cracking (MIC) framework. MIC is a regularized extended finite element method (x‐FEM) that allows the insertion of cracks in directions that are independent of the mesh orientation. The Heaviside step function that is typically used to introduce a displacement discontinuity across a crack surface is replaced by a continuous function approximated by using the original displacement shape functions. Such regularization allows the preservation of the Gaussian integration schema regardless of the enrichment required to model cracking in an arbitrary direction. The interaction between plies is anchored on the integration point distribution, which remains constant through the entire simulation. Initiation and propagation of delaminations between plies as well as intra‐ply MIC opening is implemented by using a mixed‐mode cohesive formulation in a fully three‐dimensional model that includes residual thermal stresses. The validity of the proposed methodology was tested against a variety of problems ranging from simple evolution of delamination from existing transverse cracks to strength predictions of complex laminates withouttextita priori knowledge of damage location or initiation. Good agreement with conventional numerical solutions and/or experimental data was observed in all the problems considered. Published 2011. This article is a US Government work and is in the public domain in the USA.     

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.     

New crack‐tip elements for XFEM and applications to cohesive cracks

An extended finite element method scheme for a static cohesive crack is developed with a new formulation for elements containing crack tips. This method can treat arbitrary cracks independent of the mesh and crack growth without remeshing. All cracked elements are enriched by the sign function so that no blending of the local partition of unity is required. This method is able to treat the entire crack with only one type of enrichment function, including the elements containing the crack tip. This scheme is applied to linear 3‐node triangular elements and quadratic 6‐node triangular elements. To ensure smooth crack closing of the cohesive crack, the stress projection normal to the crack tip is imposed to be equal to the material strength. The equilibrium equation and the traction condition are solved by the Newton–Raphson method to obtain the nodal displacements and the external load simultaneously. The results obtained by the new extended finite element method are compared to reference solutions and show excellent agreement. Copyright © 2003 John Wiley & Sons, Ltd.     

A framework for fracture modelling based on the material forces concept with XFEM kinematics

A theoretical and computational framework which covers both linear and non‐linear fracture behaviour is presented. As a basis for the formulation, we use the material forces concept due to the close relation between on one hand the Eshelby energy–momentum tensor and on the other hand material defects like cracks and material inhomogeneities. By separating the discontinuous displacement from the continuous counterpart in line with the eXtended finite element method (XFEM), we are able to formulate the weak equilibrium in two coupled problems representing the total deformation. However, in contrast to standard XFEM, where the direct motion discontinuity is used to model the crack, we rather formulate an inverse motion discontinuity to model crack development. The resulting formulation thus couples the continuous direct motion to the inverse discontinuous motion, which may be used to simulate linear as well as non‐linear fracture in one and the same formulation. In fact, the linear fracture formulation can be retrieved from the non‐linear cohesive zone formulation simply by confining the cohesive zone to the crack tip. These features are clarified in the two numerical examples which conclude the paper. Copyright © 2005 John Wiley & Sons, Ltd.     

Finite elements with embedded strong discontinuities for the modeling of failure in solids

This paper presents new finite elements that incorporate strong discontinuities with linear interpolations of the displacement jumps for the modeling of failure in solids. The cases of interest are characterized by a localized cohesive law along a propagating discontinuity (e.g. a crack), with this propagation occurring in a general finite element mesh without remeshing. Plane problems are considered in the infinitesimal deformation range. The new elements are constructed by enhancing the strains of existing finite elements (including general displacement based, mixed, assumed and enhanced strain elements) with a series of strain modes that depend on the proper enhanced parameters local to the element. These strain modes are designed by identifying the strain fields to be captured exactly, including the rigid body motions of the two parts of a splitting element for a fully softened discontinuity, and the relative stretching of these parts for a linear tangential sliding of the discontinuity. This procedure accounts for the discrete kinematics of the underlying finite element and assures the lack of stress locking in general quadrilateral elements for linearly separating discontinuities, that is, spurious transfers of stresses through the discontinuity are avoided. The equations for the enhanced parameters are constructed by imposing the local equilibrium between the stresses in the bulk of the element and the tractions driving the aforementioned cohesive law, with the proper equilibrium operators to account for the linear kinematics of the discontinuity. Given the locality of all these considerations, the enhanced parameters can be eliminated by their static condensation at the element level, resulting in an efficient implementation of the resulting methods and involving minor modifications of an existing finite element code. A series of numerical tests and more general representative numerical simulations are presented to illustrate the performance of the new elements. Copyright © 2007 John Wiley & Sons, Ltd.     

Embedded interfaces by polytope FEM

In this paper, the Polytope Finite Element Method is employed to model an embedded interface through the body, independent of the background FEM mesh. The elements that are crossed by the embedded interface are decomposed into new polytope elements which have some nodes on the interface line. The interface introduces discontinuity into the primary variable (strong) or into its derivatives (weak). Both strong and weak discontinuities are studied by the proposed method through different numerical examples including fracture problems with traction‐free and cohesive cracks, and heat conduction problems with Dirichlet and Dirichlet–Neumann types of boundary conditions on the embedded interface. For traction‐free cracks which have tip singularity, the nodes near the crack tip are enriched with the singular functions through the eXtended Finite Element Method. The concept of Natural Element Coordinates (NECs) is invoked to drive shape functions for the produced polytopes. A simple treatment is proposed for concave polytopes produced by a kinked interface and also for locating crack tip inside an element prior to using the singularity enrichment. The proposed method pursues some implementational details of eXtended/Generalized Finite Element Methods for interfaces. But here the additional DOFs are constructed on the interface lines in contrast to X/G‐FEM, which attach enriched DOFs to the previously existed nodes. Copyright © 2011 John Wiley & Sons, Ltd.     

Anomalous behavior of bilinear quadrilateral finite elements for modeling cohesive cracks with XFEM/GFEM

The performance of partition‐of‐unity based methods such as the generalized finite element method or the extended finite element method is studied for the simulation of cohesive cracking. The focus of investigation is on the performance of bilinear quadrilateral finite elements using these methods. In particular, the approximation of the displacement jump field, representing cohesive cracks, by extended finite element method/generalized finite element method and its effect on the overall behavior at element and structural level is investigated. A single element test is performed with two different integration schemes, namely the Newton‐Cotes/Lobatto and the Gauss integration schemes, for the cracked interface contribution. It was found that cohesive crack segments subjected to a nonuniform opening in unstructured meshes (or an inclined crack in a structured finite element mesh) result in an unrealistic crack opening. The reasons for such behavior and its effect on the response at element level are discussed. Furthermore, a mesh refinement study is performed to analyze the overall response of a cohesively cracked body in a finite element analysis. Copyright © 2013 John Wiley & Sons, Ltd.     

Effect of crack density and connectivity on the permeability of microcracked solids

In composite theory microcracks in solid are usually treated as degenerated inclusions separately embedded in matrix. For heterogeneous engineering composites like concrete and rock, the real cracking patterns are more complicate and quite different from this assumption due to the natural clustering and inter-connection of microcracks. This paper investigates the permeability of solids containing a crack network with finite connectivity following both theoretical and numerical approaches. Firstly, no connectivity is assumed for cracks and the interaction direct derivative (IDD) method is employed to obtain the crack-altered permeability of solids. Then the amplification of permeability by crack connectivity is quantified for parallel crack cases and for general crack patterns. This amplification effect is modeled by a crack length augmentation factor. In this way the IDD method is extended to evaluate the permeability of cracked solids for a finite crack connectivity before total percolation of cracks. Afterwards, by a carefully designed Monte-Carlo algorithm, the representative volume element (RVE) is built numerically for cracked solids with cracks having random spatial locations and random lengths. The permeability of 2D cracked solids is solved by finite element method (FEM). Through this numerical tool, the effect of both crack density and connectivity on the permeability is solved, and especially the relation between crack connectivity and the geometrical coefficient of crack clustering is put into evidence. From this study it is showed that the extended IDD method can be adapted to a microcracked solid with finite connectivity and can provide good estimates for the permeability.     

A stable X‐FEM in cohesive transition from closed to open crack

Ignoring crack tip effects, the stability of the X‐FEM discretizations is trivial for open cracks but remains a challenge if we constrain the crack to be closed (i.e., the bi‐material problem). Here, we develop a formulation for general cohesive interactions between crack faces within the X‐FEM framework. The stability of the new formulation is demonstrated for any cohesive crack stiffness (including the closed crack) and illustrated for a nonlinear cohesive softening law. A benchmark of the new model is carried out with simpler approaches for a closed crack (i.e., Lagrange multipliers) and for a cohesive crack (i.e., penalty approach). Due to the analogies between stable cohesive X‐FEM and Nitsche's methods, the new method simplifies the implementation and is attractive in dynamic explicit codes. Copyright © 2014 John Wiley & Sons, Ltd.     

弹塑性杆非线性断裂问题的新型增强有限元法

针对一维弹塑性材料杆件的非线性断裂,该文提出了一种新型弹塑性增强有限元。该单元采用von Mises屈服准则和线性等向硬化模型描述开裂前的材料弹塑性变形,而结合利用内聚力关系来描述随后的裂纹萌生和非线性断裂过程。引入内部节点来描述单元内由于裂纹引起的位移不连续,通过单元凝聚获得含裂纹单元的刚度矩阵,并对数值稳定性问题进行了分析。通过与基于材料力学方法推导得到的解析解的对比验证了该新型弹塑性增强有限单元法在列式上的正确性和在数值上的高效性与精确性。     

An analysis of competing toughening mechanisms in layered and particulate solids

The relative potency of common toughening mechanisms is explored for layered solids and particulate solids, with an emphasis on crack multiplication and plasticity. First, the enhancement in toughness due to a parallel array of cracks in an elastic solid is explored, and the stability of co-operative cracking is quantified. Second, the degree of synergistic toughening is determined for combined crack penetration and crack kinking at the tip of a macroscopic, mode I crack; specifically, the asymptotic problem of self-similar crack advance (penetration mode) versus $90^{\circ }$ symmetric kinking is considered for an isotropic, homogeneous solid with weak interfaces. Each interface is treated as a cohesive zone of finite strength and toughness. Third, the degree of toughening associated with crack multiplication is assessed for a particulate solid comprising isotropic elastic grains of hexagonal shape, bonded by cohesive zones of finite strength and toughness. The study concludes with the prediction of R-curves for a mode I crack in a multi-layer stack of elastic and elastic–plastic solids. A detailed comparison of the potency of the above mechanisms and their practical application are given. In broad terms, crack tip kinking can be highly potent, whereas multiple cracking is difficult to activate under quasi-static conditions. Plastic dissipation can give a significant toughening in multi-layers especially at the nanoscale.     

The effects of inter-surface cohesive tractions on linear and penny-shaped cracks

A mesoscopic fracture model of equilibrium slit cracks in brittle solids, including inter-surface cohesive tractions acting near the crack tip, is analyzed and the effects of the cohesive tractions on the in-plane stress fields, crack-opening displacement profiles, and crack driving forces examined quantitatively for linear and penny-shaped cracks. The (numerical) analysis method is described in detail, along with results for four different cohesive forces. The assumed distribution of cohesive tractions were found to suppress the in-plane stress field adjacent to cracks in a homogeneous, isotropic medium when uniformly loaded in mode-I, and the suppression was a function of crack length. The crack-opening displacement profile was also perturbed and a new regime identified between the near-field Barenblatt zone and the far-field continuum zone. The extent of this `cohesive zone' was quantified by use of an interpolating function fit to the calculated profiles and found to be independent of crack size for a given cohesive tractions distribution. Furthermore, the crack-opening displacement at the edge of the cohesive zone was found to be independent of crack size, implying that despite significant perturbations to the stress field, the crack driving force at unstable equilibrium remains unchanged with crack size.     

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.     

Theory and numerics for finite deformation fracture modelling using strong discontinuities

A general finite element approach for the modelling of fracture is presented for the geometrically non‐linear case. The kinematical representation is based on a strong discontinuity formulation in line with the concept of partition of unity for finite elements. Thus, the deformation map is defined in terms of one continuous and one discontinuous portion, considered as mutually independent, giving rise to a weak formulation of the equilibrium consisting of two coupled equations. In addition, two different fracture criteria are considered. Firstly, a principle stress criterion in terms of the material Mandel stress in conjunction with a material cohesive zone law, relating the cohesive Mandel traction to a material displacement ‘jump’ associated with the direct discontinuity. Secondly, a criterion of Griffith type is formulated in terms of the material‐crack‐driving force (MCDF) with the crack propagation direction determined by the direction of the force, corresponding to the direction of maximum energy release. Apart from the material modelling, the numerical treatment and aspects of computational implementation of the proposed approach is also thoroughly discussed and the paper is concluded with a few numerical examples illustrating the capabilities of the proposed approach and the connection between the two fracture criteria. Copyright © 2005 John Wiley & Sons, Ltd.     

An augmented cohesive zone element for arbitrary crack coalescence and bifurcation in heterogeneous materials

We demonstrate that traditional cohesive zone (CZ) elements cannot be accurate when used in conjunction with solid elements with arbitrary intra‐element cracking capability, because they cannot capture the load transfer between cohesive interfaces and the solid elements when crack bifurcation or coalescence occurs. An augmented cohesive zone (ACZ) element based on the augmented finite element method formulation is therefore proposed. The new element allows for arbitrary separation of the cohesive element in accordance with the crack configuration of the abutting solid elements, thus correctly maintaining the non‐linear coupling between merging or bifurcating cracks. Numerical accuracy and effectiveness of the proposed ACZ element are demonstrated through several examples. Copyright © 2011 John Wiley & Sons, Ltd.