非均匀材料界面裂纹的Cell-Based光滑有限元法

为提高非均匀材料界面裂纹尖端断裂参数的求解精度,基于非均匀材料界面断裂力学、Cell-Based光滑有限元(Cell-SFEM)和非均匀材料的互交作用积分法,提出了求解非均匀材料界面裂纹尖端断裂参数的CellBased光滑有限元法,推导了基于Cell-Based光滑有限元法的非均匀材料的互交作用积分法,对非均匀材料间的界面裂纹尖端处正则应力强度因子进行了求解,并与参考解进行了比较,讨论了互交积分区域大小和光滑子元个数与正则应力强度因子的关系。数值算例结果表明:本方法具有很高的计算精度,对积分区域大小不敏感,可为设计、制造抗破坏非均匀材料提供依据。     

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

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

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

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

压电材料渗透型反平面界面裂纹的奇异因子

本文用复变函数解析延展原理,研究了集中载荷作用下的不同压电材料反平面应变 状态的电渗透型界面裂纹的耦合场:对单个裂纹,给出了封闭形式的复函数解和场强度因子。 结果表明,在裂尖处耦合场有(1/2)阶的奇异性。     

压电材料中反平面渗透裂纹的解

用复变函数的保角映射法,采用可渗透边界条件,研究了含裂纹的无限大压电材料在平面内电场和反平面荷载作用下的耦合场,得到了精确的解和场强度因子以及能量释放率。结果表明,电场强度在裂尖没有奇异性,应变、应力、电位移具有1/2阶的奇异性,能量释放率总是正的。     

加层压电条界面裂纹的稳态扩展

研究当压电条同时与两个不同材料的弹性条粘接在一起,在反平面机械载荷及面内电载荷联合作用下,长度不变的有限Griffith 界面裂纹沿加层压电条界面以常速稳态扩展时裂尖的动态断裂问题.应用Fourier积分变换将问题化为以第二类Fredholm积分方程表示的对偶积分方程,导出了相应的动应力强度因子表达式.给出了动应力强度因子与裂纹传播速度、裂纹长度、压电条及弹性条厚度、电荷载大小及方向的关系曲线.研究结果对结构设计及结构失效的预防具有理论和应用价值.     

含界面边裂纹压电材料反平面问题的应力强度因子

研究了含界面边裂纹的不同压电介质组成的复合材料在反平面荷载和平面内电场作用下的电弹场,得到了级数形式的基本解和应力强度因子,最后用边界配置法求解了应力强度因子.结果表明,在外加剪切荷载的作用下,应力强度因子与外加电场无关.     

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

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

混凝土裂纹扩展过程模拟的扩展有限元法研究

扩展有限元法(theextendedFiniteElementMethod,XFEM)为数值模拟结构裂纹扩展过程提供了一条有效途径。该文介绍了用扩展有限元法对混凝土结构裂纹扩展过程进行数值模拟的实现方法。采用虚拟裂缝模型模拟混凝土非线性断裂行为,针对二维四边形单元推导了详细的有限元列式。采用3种方案对非线性方程系统进行求解,分析了其求解思路并概括了其求解步骤。通过对带初始边缘裂纹的单向拉伸混凝土板的数值模拟,对3种求解方案的计算结果进行了分析和讨论。     

含裂纹的矩形截面压电材料反平面问题的应力场和电场

研究了含裂纹的矩形截面的压电材料在平面内电场和反平面荷载作用下的问题.得到了满足拉普拉斯方程、电渗透裂纹面边界条件的位移函数解和电势函数解,从而得到了电场和弹性场的基本解.最后,用边界配置法计算了应力强度因子和能量释放率.结果表明,这种半解析半数值的方法计算简便,而且具有足够的精确性和广泛的应用性.     

Linear smoothed extended finite element method

The extended finite element method was introduced in 1999 to treat problems involving discontinuities with no or minimal remeshing through appropriate enrichment functions. This enables elements to be split by a discontinuity, strong or weak, and hence requires the integration of discontinuous functions or functions with discontinuous derivatives over elementary volumes. A variety of approaches have been proposed to facilitate these special types of numerical integration, which have been shown to have a large impact on the accuracy and the convergence of the numerical solution. The smoothed extended finite element method (XFEM), for example, makes numerical integration elegant and simple by transforming volume integrals into surface integrals. However, it was reported in the literature that the strain smoothing is inaccurate when non‐polynomial functions are in the basis. In this paper, we investigate the benefits of a recently developed Linear smoothing procedure which provides better approximation to higher‐order polynomial fields in the basis. Some benchmark problems in the context of linear elastic fracture mechanics are solved and the results are compared with existing approaches. We observe that the stress intensity factors computed through the proposed linear smoothed XFEM is more accurate than that obtained through smoothed XFEM.     

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.     

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.     

Computational modeling of strong and weak discontinuities using the extended finite element method

In this article, the extended finite element method is employed to solve problems, including weak and strong discontinuities. To this end, a level set framework is used to represent the discontinuities location, and the Heaviside and Branch function are included in the standard finite element method. The case of two arbitrary curved cracks is solved numerically and stress intensity factor (SIF) values at the crack tips are calculated based on the evaluation of the crack tip opening displacement. Afterwards, J-integral methodology is adopted to evaluate the SIFs for isotropic and anisotropic bi-material interface crack problems. Numerical results are verified with those presented in the literature.     

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

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

A contact algorithm for cohesive cracks in the extended finite element method

Cracks with quasibrittle behavior are extremely common in engineering structures. The modeling of cohesive cracks involves strong nonlinearity in the contact, material, and complex transition between contact and cohesive forces. In this article, we propose a novel contact algorithm for cohesive cracks in the framework of the extended finite element method. A cohesive-contact constitutive model is introduced to characterize the complex mechanical behavior of the fracture process zone. To avoid the stress oscillations and ill-conditioned system matrix that often occur in the conventional contact approach, the proposed algorithm employs a special dual Lagrange multiplier to impose the contact constraint. This Lagrange multiplier is constructed by means of the area-weighted average and biorthogonality conditions at the element level. The system matrix can be condensed into a positive definite matrix with an unchanged size at a very low computational cost. In addition, we illustrate solving the cohesive crack contact problem using a novel iteration strategy. Several numerical experiments are performed to illustrate the efficiency and high-quality results of our method in contact analysis of cohesive cracks.     

Fracture mechanics analysis using the wavelet Galerkin method and extended finite element method

This paper presents fracture mechanics analysis using the wavelet Galerkin method and extended finite element method. The wavelet Galerkin method is a new methodology to solve partial differential equations where scaling/wavelet functions are used as basis functions. In solid/structural analyses, the analysis domain is divided into equally spaced structured cells and scaling functions are periodically placed throughout the domain. To improve accuracy, wavelet functions are superposed on the scaling functions within a region having a high stress concentration, such as near a hole or notch. Thus, the method can be considered a refinement technique in fixed‐grid approaches. However, because the basis functions are assumed to be continuous in applications of the wavelet Galerkin method, there are difficulties in treating displacement discontinuities across the crack surface. In the present research, we introduce enrichment functions in the wavelet Galerkin formulation to take into account the discontinuous displacements and high stress concentration around the crack tip by applying the concept of the extended finite element method. This paper presents the mathematical formulation and numerical implementation of the proposed technique. As numerical examples, stress intensity factor evaluations and crack propagation analyses for two‐dimensional cracks are presented. Copyright © 2012 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.