Element‐local level set method for three‐dimensional dynamic crack growth

An approximate level set method for three‐dimensional crack propagation is presented. In this method, the discontinuity surface in each cracked element is defined by element‐local level sets (ELLSs). The local level sets are generated by a fitting procedure that meets the fracture directionality and its continuity with the adjacent element crack surfaces in a least‐square sense. A simple iterative procedure is introduced to improve the consistency of the generated element crack surface with those of the adjacent cracked elements. The discrete discontinuity is treated by the phantom node method which is a simplified version of the extended finite element method (XFEM). The ELLS method and the phantom node technology are combined for the solution of dynamic fracture problems. Numerical examples for three‐dimensional dynamic crack propagation are provided to demonstrate the effectiveness and robustness of the proposed method. Copyright © 2009 John Wiley & Sons, Ltd.     

A novel numerical model of debonding of FRP-plated concrete beam

Accurate modeling is required to estimate the debonding in a plated fiber-reinforced polymer (FRP) concrete beam. In the present investigation, a numerical method is developed to model a crack in the FRP–concrete interface. An initial notch is located at the mid-span of the concrete beam. A modified crack closure integral method is implemented to model Mode-I fracture in the concrete. In the present research, a special interface element is formulated to simulate and to predict the distribution of interfacial shear stresses by using drilling degrees of freedom in the nodes of interface elements. Cohesive forces in the nodes of interface elements are formulated by finite element methods. A crack propagation criterion is presented to evaluate when the crack grows in FRP–concrete interface. If the principal stress in the node at the tip of an interface element reaches the maximum shear stress along the FRP–concrete interface, debonding happens. The model is robust, accurate, independent of mesh size, and it is able to model the crack growth in the concrete and debonding of the FRP–concrete interface, simultaneously. The model presented in this study showed acceptable similarity to previous research data.     

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.     

Development of a new semianalytical approach for 2D analysis of crack propagation problems

This study offers novel progresses for a semianalytical method to simulate 2D propagation of cracks based on linear elastic fracture mechanics. For this purpose, a new algorithm is proposed on the basis of the linear elastic fracture mechanics for crack propagation in single‐mode and mixed mode. Herein, discretization is only performed on the boundaries of the problem by using specific subparametric elements and higher‐order Chebyshev polynomials as mapping functions. Implementing the weighted residual method and by taking Clenshaw‐Curtis numerical quadrature, diagonal Euler's differential equations are obtained. Consequently, once the local coordinate origin is assumed at the tip of the crack, the stress intensity factors can be derived directly. In accordance with the rate of maximum energy release as a propagation criterion and by proposing a new quasi‐automatic remeshing procedure based on domain division, the crack propagation is applied here. Based on the new presented algorithm, application of the new semianalytical method to crack propagation analysis is more flexible and efficient. By taking this advantage, relatively coarse and simple discretization compared with other computational methods may be used. By modelling 4 benchmark problems with a few numbers of degrees of freedom, the validity and accuracy of the current method is illustrated. Results show that the presented algorithm is applicable for efficient and precise prediction of crack trajectories.     

A new crack tip element for the phantom‐node method with arbitrary cohesive cracks

We have developed a new crack tip element for the phantom‐node method. In this method, a crack tip can be placed inside an element. Therefore, cracks can propagate almost independent of the finite element mesh. We developed two different formulations for the three‐node triangular element and four‐node quadrilateral element, respectively. Although this method is well suited for the one‐point quadrature scheme, it can be used with other general quadrature schemes. We provide some numerical examples for some static and dynamic problems. Copyright © 2007 John Wiley & Sons, Ltd.     

A quasi‐static discontinuous Galerkin configurational force crack propagation method for brittle materials

This paper presents a framework for r‐adaptive quasi‐static configurational force (CF) brittle crack propagation, cast within a discontinuous Galerkin (DG) symmetric interior penalty (SIPG) finite element scheme. Cracks are propagated in discrete steps, with a staggered algorithm, along element interfaces, which align themselves with the predicted crack propagation direction. The key novelty of the work is the exploitation of the DG face stiffness terms existing along element interfaces to propagate a crack in a mesh‐independent r‐adaptive quasi‐static fashion, driven by the CF at the crack tip. This adds no new degrees of freedom to the data structure. Additionally, as DG methods have element‐specific degrees of freedom, a geometry‐driven p‐adaptive algorithm is also easily included allowing for more accurate solutions of the CF on a moving crack front. Further, for nondeterminant systems, we introduce an average boundary condition that restrains rigid body motion leading to a determinant system. To the authors' knowledge, this is the first time that such a boundary condition has been described. The proposed formulation is validated against single and multiple crack problems with single‐ and mixed‐mode cracks, demonstrating the predictive capabilities of the method.     

Adaptive crack propagation analysis with the element‐free Galerkin method

In this paper, an adaptive analysis of crack propagation based on the error estimation by the element‐free Galerkin (EFG) method is presented. The adaptivity analysis in quasi‐static crack propagation is achieved by adding and/or removing the nodes along the background integration cells, those are refined or recovered according to the estimated errors. These errors are obtained basically by calculating the difference between the values of the projected stresses and original EFG stresses. To evaluate the performance of the proposed adaptive procedure, the crack propagation behaviour is investigated for several examples. The results of these examples show the efficiency and accuracy of the proposed scheme in crack propagation analysis. Copyright © 2002 John Wiley & Sons, Ltd.     

Enhanced mixed interpolation XFEM formulations for discontinuous Timoshenko beam and Mindlin‐Reissner plate

Shear locking is a major issue emerging in the computational formulation of beam and plate finite elements of minimal number of degrees of freedom as it leads to artificial overstiffening. In this paper, discontinuous Timoshenko beam and Mindlin‐Reissner plate elements are developed by adopting the Hellinger‐Reissner functional with the displacements and through‐thickness shear strains as degrees of freedom. Heterogeneous beams and plates with weak discontinuity are considered, and the mixed formulation has been combined with the extended finite element method (FEM); thus, mixed enrichment functions are used. Both the displacement and the shear strain fields are enriched as opposed to the traditional extended FEM where only the displacement functions are enriched. The enrichment type is restricted to extrinsic mesh‐based topological local enrichment. The results from the proposed formulation correlate well with analytical solution in the case of the beam and in the case of the Mindlin‐Reissner plate with those of a finite element package (ABAQUS) and classical FEM and show higher rates of convergence. In all cases, the proposed method captures strain discontinuity accurately. Thus, the proposed method provides an accurate and a computationally more efficient way for the formulation of beam and plate finite elements of minimal number of degrees of freedom.     

An enriched element‐failure method (REFM) for delamination analysis of composite structures

This paper develops an enriched element‐failure method for delamination analysis of composite structures. This method combines discontinuous enrichments in the extended finite element method and element‐failure concepts in the element‐failure method within the finite element framework. An improved discontinuous enrichment function is presented to effectively model the kinked discontinuities; and, based on fracture mechanics, a general near‐tip enrichment function is also derived from the asymptotic displacement fields to represent the discontinuity and local stress intensification around the crack‐tip. The delamination is treated as a crack problem that is represented by the discontinuous enrichment functions and then the enrichments are transformed to external nodal forces applied to nodes around the crack. The crack and its propagation are modeled by the ‘failed elements’ that are applied to the external nodal forces. Delamination and crack kinking problems can be solved simultaneously without remeshing the model or re‐assembling the stiffness matrix with this method. Examples are used to demonstrate the application of the proposed method to delamination analysis. The validity of the proposed method is verified and the simulation results show that both interlaminar delamination and crack kinking (intralaminar crack) occur in the cross‐ply laminated plate, which is observed in the experiment. Copyright © 2009 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.     

无网格自动加密方法及其在土体裂缝分析中的应用

将无网格法和弥散裂缝模型结合用于裂缝分析时,需要较细密的节点分布。为了准确描述土体裂缝开展,提出了一种无网格自动加密方法。通过在等参单元上布点后变换到实际的背景积分单元中,实现了裂缝开展区域的局部加密。在背景积分单元中使用有限元插值来计算新增节点和高斯积分点的状态变量。将提出的方法与弥散裂缝模型相结合,可用于土体裂缝的分析。对土梁弯曲试验进行模拟计算,验证了方法的有效性。     

Strong tangential discontinuity modeling of shear bands using the extended finite element method

A method is developed for modeling of shear band with strong tangential discontinuity by means of cohesive surfaces within the extended finite element method (XFEM). A rate-independent non-associated plasticity model is incorporated along the strong discontinuity to consider the highly localized regions. Once the localization is occurred, tangential enrichment degrees of freedom are added to the localized elements, and the discontinuity is captured regardless of mesh resolution and alignment. By introducing the tangential enrichment function, the discontinuity is only imposed in the tangential direction, while the continuity across the shear band is automatically fulfilled. Adopting bilinear quadrilateral elements within the context of XFEM allows for the plastic deformation of shear band to be obtained with quadratic distribution within an enriched element. Since the strong discontinuity approach is employed, the singularity of strain field at the position of displacement jump is attained through a Dirac delta distribution. By means of this singularity, the cohesive shear traction is derived for the J2 plasticity model and is applied onto the band interfaces in order to reproduce the dissipative mechanism of the band. Several numerical examples are analyzed to assess the accuracy and robustness of the proposed approach.     

An XFEM method for modeling geometrically elaborate crack propagation in brittle materials

We present a method for simulating quasistatic crack propagation in 2‐D which combines the extended finite element method (XFEM) with a general algorithm for cutting triangulated domains, and introduce a simple yet general and flexible quadrature rule based on the same geometric algorithm. The combination of these methods gives several advantages. First, the cutting algorithm provides a flexible and systematic way of determining material connectivity, which is required by the XFEM enrichment functions. Also, our integration scheme is straightforward to implement and accurate, without requiring a triangulation that incorporates the new crack edges or the addition of new degrees of freedom to the system. The use of this cutting algorithm and integration rule allows for geometrically complicated domains and complex crack patterns. Copyright © 2011 John Wiley & Sons, Ltd.     

Modeling of dynamic crack branching by enhanced extended finite element method

The conventional extended finite element method (XFEM) is enhanced in this paper to simulate dynamic crack branching, which is a top challenge issue in fracture mechanics and finite element method. XFEM uses the enriched shape functions with special characteristics to represent the discontinuity in computation field. In order to describe branched cracks, it is necessary to set up the additional enrichment. Here we have developed two kinds of branched elements, namely the “element crossed by two separated cracks” and “element embedded by a junction”. Another series of enriched degrees of freedom are introduced to seize the additional discontinuity in the elements. A shifted enrichment scheme is used to avoid the treatment of blending element. Correspondingly a new mass lumping method is developed for the branched elements based on the kinetic conservation. The derivation of the mass matrix of a four-node quadrilateral element which contains two strong discontinuities is specially presented. Then by choosing crack speed as the branching criterion, the branching process of a single mode I crack is simulated. The results including the branching angle and propagation routes are compared with that obtained by the conventionally used element deletion method.     

A new method for modelling cohesive cracks using finite elements

A model which allows the introduction of displacements jumps to conventional finite elements is developed. The path of the discontinuity is completely independent of the mesh structure. Unlike so‐called ‘embedded discontinuity’ models, which are based on incompatible strain modes, there is no restriction on the type of underlying solid finite element that can be used and displacement jumps are continuous across element boundaries. Using finite element shape functions as partitions of unity, the displacement jump across a crack is represented by extra degrees of freedom at existing nodes. To model fracture in quasi‐brittle heterogeneous materials, a cohesive crack model is used. Numerical simulations illustrate the ability of the method to objectively simulate fracture with unstructured meshes. Copyright © 2001 John Wiley & Sons, Ltd.     

A general mass lumping scheme for the variants of the extended finite element method

A new strategy for the mass matrix lumping of enriched elements for explicit transient analysis is presented. It is shown that to satisfy the kinetic energy conservation, the use of zero or negative masses for enriched degrees of freedom of lumped mass matrix may be necessary. For a completely cracked element, by lumping the mass of each side of the interface into the finite element nodes located at the same side and assigning zero masses to the enriched degrees of freedom, the kinetic energy for rigid body translations is conserved without transferring spurious energy across the interface. The time integration is performed by adopting an explicit-implicit technique, where the regular and enriched degrees of freedom are treated explicitly and implicitly, respectively. The proposed method can be viewed as a general mass lumping scheme for the variants of the extended finite element methods because it can be used irrespective of the enrichment method. It also preserves the optimal critical time step of an intact finite element by treating the enriched degrees of freedom implicitly. The accuracy and efficiency of the proposed mass matrix are validated with several benchmark examples.     

An assumed-stress hybrid 4-node shell element with drilling degrees of freedom

An assumed-stress hybrid/mixed 4-node quadrilateral shell element is introduced that alleviates most of the deficiencies associated with such elements. The formulation of the element is based on the assumed-stress hybrid/mixed method using the Hellinger-Reissner variational principle. The membrane part of the element has 12 degrees of freedom including rotational or ‘drilling’ degrees of freedom at the nodes. The bending part of the element also has 12 degrees of freedom. The bending part of the element uses the Reissner-Mindlin plate theory which takes into account the transverse shear contributions. The element formulation is derived from an 8-node isoparametric element by expressing the midside displacement degrees of freedom in terms of displacement and rotational degrees of freedom at corner nodes. The element passes the patch test, is nearly insensitive to mesh distortion, does not ‘lock’, possesses the desirable invariance properties, has no hidden spurious modes, and for the majority of test cases used in this paper produces more accurate results than the other elements employed herein for comparison.     

An efficient numerical method for modelling initiation and propagation of cracks along material interfaces

A numerical method, based on the multiregion concept in the boundary element method, is presented for modelling crack propagation. The advantage of the method is that the non-linear iterations only involve interface degrees of freedom and that the propagation of cracks is modelled by disconnecting degrees of freedom and no remeshing is necessary.     

Numerical simulation of crack deflection and penetration at an interface in a bi-material under dynamic loading by time-domain boundary element method

The hybrid time-domain boundary element method (BEM), together with the multi-region technique, is applied to simulate the dynamic process of crack deflection/ penetration at an interface in a bi-material. The whole bi-material is divided into two regions along the interface. The traditional displacement boundary integral equations (BIEs) are employed with respect to the exterior boundaries; meanwhile, the non-hypersingular traction BIEs are used with respect to the part of the crack in the matrix. Crack propagation along the interface is numerically modelled by releasing the nodes in the front of the moving crack tip and crack propagation in the matrix is modeled by adding new elements of constant length to the moving crack tip. The dynamic behaviours of the crack deflection/penetration at an interface, propagation in the matrix or along the interface and kinking out off the interface, are controlled by criteria developed from the quasi-static ones. The numerical results of the crack growth trajectory for different inclined interface and bonded strength are computed and compared with the corresponding experimental results. Agreement between numerical and experimental results implies that the present time-domain BEM can provide a simulation for the dynamic propagation and deflection of a crack in a bi-material.     

The partition of unity finite element method for elastic wave propagation in Reissner–Mindlin plates

This paper reports a numerical method for modelling the elastic wave propagation in plates. The method is based on the partition of unity approach, in which the approximate spectral properties of the infinite dimensional system are embedded within the space of a conventional finite element method through a consistent technique of waveform enrichment. The technique is general, such that it can be applied to the Lagrangian family of finite elements with specific waveform enrichment schemes, depending on the dominant modes of wave propagation in the physical system. A four‐noded element for the Reissner–Mindlin plate is derived in this paper, which is free of shear locking. Such a locking‐free property is achieved by removing the transverse displacement degrees of freedom from the element nodal variables and by recovering the same through a line integral and a weak constraint in the frequency domain. As a result, the frequency‐dependent stiffness matrix and the mass matrix are obtained, which capture the higher frequency response with even coarse meshes, accurately. The steps involved in the numerical implementation of such element are discussed in details. Numerical studies on the performance of the proposed element are reported by considering a number of cases, which show very good accuracy and low computational cost. Copyright © 2006 John Wiley & Sons, Ltd.