求解双材料界面裂纹应力强度因子的扩展有限元法

基于双材料界面裂纹尖端的基本解,构造扩展有限元法(eXtended Finite Element Methods, XFEM)裂尖单元结点的改进函数。有限元网格剖分不遵从材料界面,考虑3种类型的结点改进函数:弱不连续改进函数、Heaviside改进函数和裂尖改进函数,建立XFEM的位移模式,给出计算双材料界面裂纹应力强度因子(Stress Intensity Factors, SIFs)的相互作用积分方法。数值结果表明:XFEM无需遵从材料界面剖分网格,该文的方法能够准确评价双材料界面裂纹尖端的SIFs。     

Material‐dependent crack‐tip enrichment functions in XFEM for modeling interfacial cracks in bimaterials

A novel set of enrichment functions within the framework of the extended finite element method is proposed for linear elastic fracture analysis of interface cracks in bimaterials. The motivation for the new enrichment set stems from the revelation that the accuracy and conditioning of the widely accepted 12‐fold bimaterial enrichment functions significantly deteriorates with the increase in material mismatch. To this end, we propose an 8‐fold material‐dependent enrichment set, derived from the analytical asymptotic displacement field, that well captures the near‐tip oscillating singular fields of interface cracks, including the transition to weak discontinuities of bimaterials. The performance of the proposed material‐dependent enrichment functions is studied on 2 benchmark examples. Comparisons are made with the 12‐fold bimaterial enrichment as well as the classical 4‐fold homogeneous branch functions, which have also been used for bimaterials. The numerical studies clearly demonstrate the superiority of the new enrichment functions, which yield the most accurate results but with less number of degrees of freedom and significantly improved conditioning than the 12‐fold functions.     

Partition of unity enrichment for bimaterial interface cracks

Partition of unity enrichment techniques are developed for bimaterial interface cracks. A discontinuous function and the two‐dimensional near‐tip asymptotic displacement functions are added to the finite element approximation using the framework of partition of unity. This enables the domain to be modelled by finite elements without explicitly meshing the crack surfaces. The crack‐tip enrichment functions are chosen as those that span the asymptotic displacement fields for an interfacial crack. The concept of partition of unity facilitates the incorporation of the oscillatory nature of the singularity within a conforming finite element approximation. The mixed‐mode (complex) stress intensity factors for bimaterial interfacial cracks are numerically evaluated using the domain form of the interaction integral. Good agreement between the numerical results and the reference solutions for benchmark interfacial crack problems is realized. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

The extended finite element method for fracture in composite materials

Methods for treating fracture in composite material by the extended finite element method with meshes that are independent of matrix/fiber interfaces and crack morphology are described. All discontinuities and near‐tip enrichments are modeled using the framework of local partition of unity. Level sets are used to describe the geometry of the interfaces and cracks so that no explicit representation of either the cracks or the material interfaces are needed. Both full 12 function enrichments and approximate enrichments for bimaterial crack tips are employed. A technique to correct the approximation in blending elements is used to improve the accuracy. Several numerical results for both two‐dimensional and three‐dimensional examples illustrate the versatility of the technique. The results clearly demonstrate that interface enrichment is sufficient to model the correct mechanics of an interface crack. Copyright © 2008 John Wiley & Sons, Ltd.     

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.     

Stable 3D XFEM/vector level sets for non‐planar 3D crack propagation and comparison of enrichment schemes

We present a three‐dimensional vector level set method coupled to a recently developed stable extended finite element method (XFEM). We further investigate a new enrichment approach for XFEM adopting discontinuous linear enrichment functions in place of the asymptotic near‐tip functions. Through the vector level set method, level set values for propagating cracks are obtained via simple geometrical operations, eliminating the need for solution of differential evolution equations. The first XFEM variant ensures optimal convergence rates by means of geometrical enrichment, ie, the use of enriched elements in a fixed volume around the crack front, without giving rise to conditioning problems. The linear enrichment approach, significantly simplifies implementation and reduces the computational cost associated with numerical integration, while providing nonoptimal convergence rates similar to standard finite elements. The 2 dicretization schemes are tested for different benchmark problems, and their combination to the vector level set method is verified for nonplanar crack propagation problems.     

An edge-based smoothed XFEM for fracture in composite materials

In this work, an edge-based smoothed extended finite element method (ES-XFEM) is extended to fracture analysis in composite materials. This method, in which the edge-based smoothing technique is married with enrichment in XFEM, shows advantages of both the extended finite element method (XFEM) and the edge-based smoothed finite element method (ES-FEM). The crack tip enrichment functions are specially derived to represent the characteristic of the displacement field around the crack tip in composite materials. Due to the strain smoothing, the necessity of integrating the singular derivatives of the crack tip enrichment functions is eliminated by transforming area integration into path integration, which is an obvious advantage compared with XFEM. Two examples are presented to testify the accuracy and convergence rate of the ES-XFEM.     

二维压电材料动强度因子的扩展有限元计算

采用扩展有限元求解二维弹性压电材料动断裂问题。扩展有限元的网格独立于裂纹,因此网格生成可大大地简化,且裂纹扩展时不需重构网格。采用相互作用积分技术计算动强度因子。比较了标准的力裂尖加强函数和力-电裂尖加强函数对动强度因子的影响,结果表明标准的力裂尖加强函数能有效地分析压电材料动断裂问题。分析了极化方向对动强度因子的影响。数值分析表明采用扩展有限元获得的动强度因子与其他数值方法解吻合得很好。     

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.     

XFEM analysis of the effects of voids,inclusions and other cracks on the dynamic stress intensity factor of a major crack

This paper is devoted to the extraction of the dynamic stress intensity factor (DSIF) for structures containing multiple discontinuities (cracks, voids and inclusions) by developing the extended finite element method (XFEM). In this method, four types of enrichment functions are used in the framework of the partition of unity to model interface discontinuity within the classical finite element method. In this procedure, elements that include a crack segment, the boundary of a void or the boundary of an inclusion are not required to conform to discontinuous edges. The DSIF is evaluated by the interaction integral. After the effectiveness of the implemented XFEM program is verified, the effects of voids, inclusions and other cracks on the DSIF of a stationary major crack are investigated by using XFEM. The results show that the dynamic effects have an influence on the path independence of the interaction integral, and these voids, inclusions and other cracks have a significant effect on the DSIF of the major crack.     

Numerical simulation of bi-material interfacial cracks using EFGM and XFEM

In this paper, bi-material interfacial cracks have been simulated using element free Galerkin method (EFGM) and extended finite element method (XFEM) under mode-I and mixed mode loading conditions. Few crack interaction problems of dissimilar layered materials are also simulated using extrinsic partition of unity enriched approach. Material discontinuity has been modeled by a signed distance function whereas strong discontinuity has been modeled by two functions i.e. Heaviside and asymptotic crack tip enrichment functions. The stress intensity factors for bi-material interface cracks are numerically evaluated using the modified domain form of interaction integral. The results obtained by EFGM and XFEM for bi-material edge and center cracks are compared with those available in literature. In order to check the validity of simulations, the results have been obtained for two different ratio of Young’s modulus.     

A consistent partly cracked XFEM element for cohesive crack growth

Present extended finite element method (XFEM) elements for cohesive crack growth may often not be able to model equal stresses on both sides of the discontinuity when acting as a crack‐tip element. The authors have developed a new partly cracked XFEM element for cohesive crack growth with extra enrichments to the cracked elements. The extra enrichments are element side local and were developed by superposition of the standard nodal shape functions for the element and standard nodal shape functions for a sub‐triangle of the cracked element. With the extra enrichments, the crack‐tip element becomes capable of modelling variations in the discontinuous displacement field on both sides of the crack and hence also capable of modelling the case where equal stresses are present on each side of the crack. The enrichment was implemented for the 3‐node constant strain triangle (CST) and a standard algorithm was used to solve the non‐linear equations. The performance of the element is illustrated by modelling fracture mechanical benchmark tests. Investigations were carried out on the performance of the element for different crack lengths within one element. The results are compared with previously obtained XFEM results applying fully cracked XFEM elements, with computational results achieved using standard cohesive interface elements in a commercial code, and with experimental results. The suggested element performed well in the tests. Copyright © 2007 John Wiley & Sons, Ltd.     

An optimally convergent discontinuous Galerkin‐based extended finite element method for fracture mechanics

The extended finite element method (XFEM) enables the representation of cracks in arbitrary locations of a mesh. We introduce here a variant of the XFEM rendering an optimally convergent scheme. Its distinguishing features are as follows: (a) the introduction of singular asymptotic crack tip fields with support on only a small region around the crack tip (the enrichment region), (b) only one and two enrichment functions are added for anti‐plane shear and planar problems, respectively and (c) the relaxation of the continuity between the enrichment region and the rest of the domain, and the adoption of a discontinuous Galerkin (DG) method therein. The method is provably stable for any positive value of a stabilization parameter, and by weakly enforcing the continuity between the two regions it eliminates ‘blending elements’ partly responsible for the suboptimal convergence of some early XFEMs. Moreover, the particular choice of enrichment functions results in a surprisingly sparse stiffness matrix that remains reasonably conditioned as the mesh is refined. More importantly, the stress intensity factors can be extracted with a satisfactory accuracy as primary unknowns. Quadrature strategies required for the optimal convergence are also discussed. Finally, the DG method was modified to retain stability based on an inf‐sup condition. Copyright © 2009 John Wiley & Sons, Ltd.     

Error estimation for crack simulations using the XFEM

The extended finite element method (XFEM) is by now well‐established for crack calculations in linear elastic fracture mechanics. An advantage of this method is its discretization independence for crack simulations. Nevertheless, discretization errors occur when using the XFEM. In this paper, a simple recovery based error estimator for the discretization error in XFEM‐calculations for cracks is presented. The method is based on the Zienkiewicz and Zhu error estimator. Enhanced smoothed stresses incorporating the discontinuities and singularities because of the cracks are recovered to enable the error estimation for arbitrary distributed cracks. This approach also allows the consideration of materials with generally inelastic behaviour. The enhanced stresses are computed by means of a least square fit problem. To assess the quality of the error estimator, global norms and the effectivity index for the global energy norm for examples with known analytical solutions are presented. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

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.     

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.     

Cohesive and non-cohesive fracture by higher-order enrichment of XFEM

A comprehensive study is performed on the use of higher-order terms of the crack tip asymptotic fields as enriching functions for the eXtended FEM (XFEM) for both cohesive and traction-free cracks. For traction-free cracks, the Williams asymptotic field is used to obtain highly accurate stress intensity factors (SIFs), directly from the enriched degrees of freedom without any post-processing. The low accuracy of the results of the original research on this subject by Liu et al. [Int. J. Numer. Meth. Engng., 2004; 59:1103–1118] is remedied here by appropriate modifications of the enrichment scheme. The modifications are simple and can be easily included into an XFEM computer code. For cohesive cracks, the relevant asymptotic field is used, and two widely used criteria including the SIFs criterion and the stress criterion are examined for the crack growth simulation. Both linear and nonlinear cohesive laws are used. For the stress criterion, averaging is avoided due to the highly accurate crack tip approximation because of the higher-order enrichment. Then, a modified stress criterion is proposed, which is shown to be applicable to a wider class of problems. Several numerical examples, including straight and curved cracks, stationary and growing cracks, single and multiple cracks, and traction-free and cohesive cracks, are studied to investigate in detail the robustness and efficiency of the proposed enrichment scheme. Copyright © 2012 John Wiley & Sons, Ltd.