XFEM simulation of cracks, holes and inclusions in functionally graded materials

The present work aims at the numerical simulation of inhomogeneities/discontinuities (cracks, holes and inclusions) in functionally graded materials (FGMs) using extended finite element method (XFEM). A FGM with unidirectional gradation in material properties is modeled under plane strain condition. The domain contains a major crack either at the center or at the edge of the domain along with multiple minor discontinuities/flaws such as minor cracks and/or voids/inclusions distributed all over the domain. The effect of the variation in stress intensity factor (SIF) of the major crack due to the presence of the minor cracks and voids/inclusions is studied in detail. The simulations show that the presence of minor discontinuities significantly affects the values of SIFs.     

模拟三维裂纹问题的自适应多尺度扩展有限元法

为了在大型结构分析中考虑小裂纹或以小的代价提高裂纹附近求解精度,该文建立了分析三维裂纹问题的自适应多尺度扩展有限元法。基于恢复法评估三维扩展有限元后验误差,大于给定误差值的单元进行细化。所有尺度单元采用八结点六面体单元,采用六面体任意结点单元连接不同尺度单元。采用互作用积分法计算三维应力强度因子。三维I 型裂纹和I-II 复合型裂纹算例分析表明了该方法的正确性和有效性。     

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

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

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.     

Frictional contact analysis of multiple cracks by incremental displacement and resultant traction boundary integral equations

Based on the merits of the dual boundary element technique, a modified dual boundary element technique is extended to deal with the frictional contact of a finite plate with arbitrarily distributed multiple cracks. Besides establishing the incremental displacement boundary integral equation on the outer boundary, the resultant traction boundary integral equation on one of the crack surfaces is also developed. Since the resultant traction instead of incremental traction on the crack surface is introduced, the computed resultant contact tractions under sliding condition satisfy the Coulomb's friction law directly. Hence, as compared with the authors' previous work, only very few computation iterations are required by this method to accurately describe the contact situations of crack surfaces. As a result, not only the linear cracks, but also other types of multiple cracks, for example, curved and kinked cracks, can be tackled. The effects of friction and interaction among cracks on the computation of stress intensity factors are also displayed.     

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

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

Multiple holes,cracks, and inclusions in anisotropic viscoelastic solids

By using the elastic–viscoelastic correspondence principle, the problems with multiple holes, cracks, and inclusions in two-dimensional anisotropic viscoelastic solids are solved for the cases with time-invariant boundaries. Based upon this principle and the existing methods for the problems with anisotropic elastic materials, two different approaches are proposed in this paper. One is concerned with an analytical solution for certain specific cases such as two collinear cracks, collinear periodic cracks, and interaction between inclusion and crack, and the other is a boundary-based finite element method for the general cases with multiple holes, cracks, and inclusions. The former considers only specific cases in infinite domain and can be used as a reference for any other numerical methods, and the latter is applicable to any combination of holes, cracks and inclusions in finite domain, whose number, size and orientation are not restricted. Unlike the conventional finite element method or boundary element method which usually needs very fine meshes to get convergence solutions, in the proposed boundary-based finite element method no meshes are needed along the boundaries of holes, cracks and inclusions. To show the accuracy and efficiency of these two proposed approaches, several representative examples are implemented analytically and numerically, and they are compared with each other or with the solutions obtained by the finite element method.     

Solution of cracks developing from stress raisers, in finite bodies, by a hybrid dislocation density/BEM approach

A method is described for determining crack tip stress intensity factors for cracks growing from stress raisers such as voids or inclusions, but in bodies having finite external dimensions. This is based on a combination of the dislocation density approach, where the influence of the stress raiser is automatically included, and the boundary element method, which is used to account for the presence of the remote surfaces.     

A study of free surface effects on through cracks under impact loading

This paper deals with the implementation of a three-dimensional time-domain boundary integral formulation for a center-crack, finite solid under symmetrically applied step loading. The BEM displacement time domain formulations have, hitherto, been limited to analyzing two-dimensional crack problems, though hypersingular formulations have been used to analyze finite cracks in infinite domains. In this paper, variation of dynamic stress intensity factor (DSIF) along the crack front for a stationary, through-thickness straight crack is studied for a finite solid under step loading. The state of stress is evaluated at the crack vertex, where crack front meets the free surface. The effect of free surface on DSIF is investigated. The effect of waves traveling in thickness direction is explained. It is possible to estimate accurately the critical intersection angle of the crack front with the free surface at which square-root singularity is restored at the crack vertex under step loading. A new partitioning scheme is proposed for spatial integration of elastodynamic kernels.     

基于X-SBFEM的裂纹体非网格重剖分耦合模型研究

扩展有限元法利用了非网格重剖分技术,但需要基于裂尖解析解构造复杂的插值基函数,计算精度受网格疏密和插值基函数等因素影响。比例边界有限元法则在求解无限域和裂尖奇异性问题优势明显,两者衔接于有限元法理论内,可建立一种结合二者优势的断裂耦合数值模型。该文从虚功原理出发,利用位移协调与力平衡机制,提出了一种断裂计算的新方法X-SBFEM,达到了扩展有限元模拟裂纹主体、比例边界有限元模拟裂尖的目的。在数值算例中,通过边裂纹和混合型裂纹的应力强度因子计算,并与理论解对比,验证了该方法的准确性和有效性。     

An investigation of dynamic interaction between multiple cracks and inclusions by TDBEM

In this paper, the dynamic interaction between multiple inclusions and cracks is studied by the time-domain boundary element method (TDBEM). To deal with this problem, two kinds of time-domain boundary integral equations together with the sub-region technique are applied. The cracked solid is divided into homogeneous and isotropic sub-regions bounded by the interfaces between the inclusions and the matrix. The non-hypersingular traction boundary integral equations are applied on the crack-surfaces; while the traditional displacement boundary integral equations are used on the interfaces and the exterior boundaries. In the numerical solution procedure, square-root shape functions are adopted for the crack-opening-displacements to describe the proper asymptotic behavior in the vicinity of the crack-tips. Numerical results for dynamic stress intensity factors are presented for various cases. The effects of the inclusion position, material combinations and multiple micro-cracks on the dynamic stress intensity factors are discussed.     

Dynamic crack propagation analysis of orthotropic media by the extended finite element method

Dynamic crack propagation of composites is investigated in this paper based on the recent advances and development of orthotropic enrichment functions within the framework of partition of unity and the extended finite element method (XFEM). The method allows for analysis of the whole crack propagation pattern on an unaltered finite element mesh, defined independent of the existence of any predefined crack or its propagation path. A relatively simple, though efficient formulation is implemented, which consists of using a dynamic crack initiation toughness, a crack orientation along the maximum circumferential stress, and a simple equation to presume the crack speed. Dynamic stress intensity factors (DSIFs) are evaluated by means of the domain separation integral method. The governing elastodynamics equation is first transformed into a standard weak formulation and is then discretized into an XFEM system of time dependent equations, to be solved by the unconditionally stable Newmark time integration scheme. A number of benchmark and test problems are simulated and the results are compared with available reference results.     

Delamination analysis of composites by new orthotropic bimaterial extended finite element method

The extended finite element method (XFEM) is further improved for fracture analysis of composite laminates containing interlaminar delaminations. New set of bimaterial orthotropic enrichment functions are developed and utilized in XFEM analysis of linear‐elastic fracture mechanics of layered composites. Interlaminar crack‐tip enrichment functions are derived from analytical asymptotic displacement fields around a traction‐free interfacial crack. Also, heaviside and weak discontinuity enrichment functions are utilized in modeling discontinuous fields across interface cracks and bimaterial weak discontinuities, respectively. In this procedure, elements containing a crack‐tip or strong/weak discontinuities are not required to conform to those geometries. In addition, the same mesh can be used to analyze different interlaminar cracks or delamination propagation. The domain interaction integral approach is also adopted in order to numerically evaluate the mixed‐mode stress intensity factors. A number of benchmark tests are simulated to assess the performance of the proposed approach and the results are compared with available reference results. Copyright © 2011 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.     

A numerical method for a void–crack interaction under cyclic loads

This paper presents a numerical method to model a general system containing cracks and voids in an infinite elastic plate under remote cyclic loads. By extending Bueckner’s principle suited for a crack to a general system containing cracks and voids, the original problem is divided into a homogeneous problem (the one without cracks and voids) subjected to remote loads and a void-crack problem in an unloaded body with applied tractions on the surfaces of cracks and voids. Thus, the results in terms of the stress intensity factors can be calculated by considering the latter problem, which is analyzed easily by using the hybrid displacement discontinuity method (a boundary element method). Further, a fatigue growth technique of a mixed-mode crack is combined with the numerical approach to simulating a void–crack interaction problem under cyclic loads. Test examples are included to illustrate that the numerical method is very simple and effective for analyzing a void–crack interaction problem.     

Stress intensity factors evaluation for through-transverse crack in slab track system under vehicle dynamic load

The stress intensity factors (SIFs) for through-transverse crack in the China Railway Track System (CRTS II) slab track system under vehicle dynamic load are evaluated in this paper. A coupled dynamic model of a half-vehicle and the slab track is presented in which the half-vehicle is treated as a 18-degree-of-freedom multi-body system. The slab track is modeled as two continuous Bernoulli–Euler beams supported by a series of elastic rectangle plates on a viscoelastic foundation. The model is applied to calculate the vertical and lateral dynamic wheel–rail forces. A three-dimensional finite element model of the slab track system is then established in which the through-transverse crack at the bottom of concrete base is created by using extended finite element method (XFEM). The wheel–rail forces obtained by the vehicle-track dynamics calculation are utilized as the inputs to finite element model, and then the values of dynamic SIFs at the crack-tip are extracted from the XFEM solution by domain based interaction integral approach. The influences of subgrade modulus, crack length, crack angle, friction coefficient between cracked surfaces, and friction coefficient between faces of concrete base and subgrade on dynamic SIFs are investigated in detail. The analysis indicates that the subgrade modulus, crack length and crack angle have great effects on dynamic SIFs at the crack-tip, while both of the friction coefficients have negligible influences on variations of dynamic SIFs. Also the statistical characteristics of varying SIFs due to random wheel–rail forces are studied and results reveal that the distributions of dynamic SIFs follow an approximately Gaussian distribution with different mean values and standard deviations. The numerical results obtained are very useful in the maintenance of the slab track system.     

扩展有限元法在裂纹扩展问题中的应用

扩展有限元法(Extended finite element method,XFEM)是近几年发展起来的数值方法,属于传统有限元法的扩展,具有区别于传统有限元法的优点,在求解不连续断裂问题上具有更高的精度及效率。本文针对影响裂纹扩展的主要因素进行探讨,继而介绍扩展有限元的基本原理,并对其在裂纹扩展中的应用进行综述,同时对该方法的下一步研究进行了展望。     

Dynamic stress intensity factors around two parallel cracks in a functionally graded layer bonded to dissimilar half-planes subjected to anti-plane incident harmonic stress waves

The time-harmonic problem of determining the stress field around two parallel cracks in functionally graded materials (FGMs) is studied. The Fourier transform technique is used to reduce the boundary conditions to four simultaneous integral equations which are then solved by expanding the differences of crack surface displacements in a series. The unknown coefficients in the series are obtained by the Schmidt method. Numerical calculations are carried out for dynamic stress intensity factors (DSIF) in FGMs.     

The integral equation formulations of an infinite elastic medium containing inclusions, cracks and rigid lines

The integral equation formulations of an infinite homogeneous isotropic medium containing various inclusions, cracks and rigid lines are presented. The present integral equation formulations contain the displacements (no tractions) over the inclusion-matrix interfaces, the discontinuous displacements over crack surfaces and the axial and the shear forces along rigid-line inclusions. Besides, the sub-domain boundary element method is also used in the present research. Numerical results from the present method and the sub-domain boundary element method are compared and discussed.     

Representation of singular fields without asymptotic enrichment in the extended finite element method

In this paper, we replace the asymptotic enrichments around the crack tip in the extended finite element method (XFEM) with the semi‐analytical solution obtained by the scaled boundary finite element method (SBFEM). The proposed method does not require special numerical integration technique to compute the stiffness matrix, and it improves the capability of the XFEM to model cracks in homogeneous and/or heterogeneous materials without a priori knowledge of the asymptotic solutions. A Heaviside enrichment is used to represent the jump across the discontinuity surface. We call the method as the extended SBFEM. Numerical results presented for a few benchmark problems in the context of linear elastic fracture mechanics show that the proposed method yields accurate results with improved condition number. A simple code is annexed to compute the terms in the stiffness matrix, which can easily be integrated in any existing FEM/XFEM code. Copyright © 2013 John Wiley & Sons, Ltd.