Hypersingular integral equation method for a three-dimensional crack in anisotropic electro-magneto-elastic bimaterials

Using Green’s functions, the extended general displacement solutions of a three-dimensional crack problem in anisotropic electro-magneto-elastic (EME) bimaterials under extended loads are analyzed by the boundary element method. Then, the crack problem is reduced to solving a set of hypersingular integral equations (HIE) coupled with boundary integral equations. The singularity of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of HIE, and the exact analytical solutions of the extended singular stresses and extended stress intensity factors (SIFs) near the crack front in anisotropic EME bimaterials are given. Also, the numerical method of the HIE for a rectangular crack subjected to extended loads is put forward with the extended crack opening dislocation approximated by the product of basic density functions and polynomials. At last, numerical solutions of the extended SIFs of some examples are obtained.     

Mixed-mode stress intensity factors of 3D interface crack in fully coupled electromagnetothermoelastic multiphase composites

This contribution presents an extended hypersingular intergro-differential equation (E-HIDE) method for modeling the 3D interface crack problem in fully coupled electromagnetothermoelastic anisotropic multiphase composites under extended electro-magneto-thermo-elastic coupled loads through theoretical analysis and numerical simulations. First, based on the extended boundary element method, the 3D interface crack problem is reduced to solving a set of E-HIDEs coupled with extended boundary integral equations, in which the unknown functions are the extended displacement discontinuities. Then, the behavior of the extended singular stress indices around the interface crack front terminating at the interface is analyzed by the extended main-part analysis. The extended stress intensity factors near the crack front are defined. In addition, a numerical method for a 3D interface crack problem subjected to extended loads is proposed, in which the extended displacement discontinuities are approximated by the product of basic density functions and polynomials. Finally, the radiation distribution of extended stress intensity factors at the interface crack surface are calculated, and the results are presented toward demonstrating the applicability of the proposed method.     

Application of hypersingular integral equation method to three-dimensional crack in electromagnetothermoelastic multiphase composites

A three-dimensional crack problem in electromagnetothermoelastic multiphase composites (EMTE-MCs) under extended loads is investigated in this paper. Using Green’s functions, the extended general displacement solutions are obtained by the boundary element method. This crack problem is reduced to solving a set of hypersingular integral equations coupled with boundary integral equations, in which the unknown functions are the extended displacement discontinuities. Then, the behavior of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of hypersingular integral equations. Analytical solutions for the extended singular stresses, the extended stress intensity factors (SIFs) and the extended energy release rate near the crack front in EMTE-MCs are provided. Also, a numerical method of the hypersingular integral equations for a rectangular crack subjected to extended loads is put forward with the extended displacement discontinuities approximated by the product of basic density functions and polynomials. In addition, distributions of extended SIFs varying with the shape of the crack are presented. The results show that the present method accurately yields smooth variations of extended SIFs along the crack front.     

Analysis of multiple interfacial cracks in three-dimensional bimaterials using hypersingular integro-differential equation method

By using the concept of finite-part integral, a set of hypersingular integro-differential equations for multiple interracial cracks in a three-dimensional infinite bimaterial subjected to arbitrary loads is derived. In the numerical analysis, unknown displacement discontinuities are approximated with the products of the fundamental density functions and power series. The fundamental functions are chosen to express a two-dimensional interface crack rigorously. As illustrative examples, the stress intensity factors for two rectangular interface cracks are calculated for various spacing, crack shape and elastic constants. It is shown that the stress intensity factors decrease with the crack spacing.     

Three-dimensional dynamic Green’s functions in transversely isotropic bi-materials

By virtue of a complete representation using two displacement potentials, an analytical derivation of the elastodynamic Green’s functions for a linear elastic transversely isotropic bi-material full-space is presented. Three-dimensional point-load Green’s functions for stresses and displacements are given in complex-plane line-integral representations. The formulation includes a complete set of transformed stress–potential and displacement–potential relations, within the framework of Fourier expansions and Hankel integral transforms, that is useful in a variety of elastodynamic as well as elastostatic problems. For numerical computation of the integrals, a robust and effective methodology is laid out which gives the necessary account of the presence of singularities including branch points and pole on the path of integration. As illustrations, the present Green’s functions are analytically degenerated to the special cases such as half-space, surface and full-space Green’s functions. Some typical numerical examples are also given to show the general features of the bi-material Green’s functions.     

Anti-plane fracture analysis for the weak-discontinuous interface in a non-homogeneous piezoelectric bi-material structure

The concept of weak discontinuity is extended to functionally graded piezoelectric bi-material interface, and fracture analysis for the weak discontinuous interface is performed by the methods of Fourier integral transform and Cauchy singular integral equation. Numerical results of the total energy release rate (TERR) and the mechanical strain energy release rate (MSERR) are obtained to show the effects of non-homogeneity parameters, geometrical parameters and loads. Parametric studies yield three conclusions: (1) To reduce the weak-discontinuity of the interface is beneficial to resisting interfacial fracture. The effect of the weak-discontinuity of the interface on TERR and MSERR still depends on the strip width. The wider the strip, the more sensitive the TERR and MSERR will be to the weak-discontinuity of the interface. (2) To predict the effect of electric load on crack propagation, MSERR is more appropriate than TERR to be used as a fracture parameter. To predict the effect of mechanical load on crack propagation, both of them could be used as fracture parameters, and MSERR is more conservative. (3) Mechanical load and negative electric displacement load would promote crack propagation, but positive electric displacement load would retard it. For the structure applied by combined mechanical and positive electric displacement loads, crack propagation may be impeded by appropriately selecting the strip width and the ratio of non-homogeneity parameters.     

Evaluation of stress intensity factors for bi-material interface cracks using displacement jump methods

The aim of the present work is to investigate the numerical modeling of interfacial cracks that may appear at the interface between two isotropic elastic materials. The extended finite element method is employed to analyze brittle and bi-material interfacial fatigue crack growth by computing the mixed mode stress intensity factors (SIF). Three different approaches are introduced to compute the SIFs. In the first one, mixed mode SIF is deduced from the computation of the contour integral as per the classical J-integral method, whereas a displacement method is used to evaluate the SIF by using either one or two displacement jumps located along the crack path in the second and third approaches. The displacement jump method is rather classical for mono-materials, but has to our knowledge not been used up to now for a bi-material. Hence, use of displacement jump for characterizing bi-material cracks constitutes the main contribution of the present study. Several benchmark tests including parametric studies are performed to show the effectiveness of these computational methodologies for SIF considering static and fatigue problems of bi-material structures. It is found that results based on the displacement jump methods are in a very good agreement with those of exact solutions, such as for the J-integral method, but with a larger domain of applicability and a better numerical efficiency (less time consuming and less spurious boundary effect).     

双相介质半空间中圆孔对透射SH波的稳态分析

研究了平面SH波在半空间双相弹性介质中的传播。通过Green函数和积分方程方法,按照复变函数描述,对透射波被圆孔散射的情况进行稳态分析。将双相介质半空间沿界面剖分为1/4空间介质Ⅰ和含圆孔的1/4空间介质Ⅱ,分别构造了介质Ⅰ和介质Ⅱ中反平面点源荷载的Green函数,按双相介质中平面SH波的处理方法,给出介质Ⅰ和介质Ⅱ中的平面位移波,两种介质之间的相互作用力与对应Green函数的乘积沿界面的积分与平面位移波叠加得到介质Ⅰ和介质Ⅱ中的全部位移场。按照界面的位移连续条件,定解积分方程组,得到问题的稳态解,并给出圆孔位置和介质参数对散射的影响。     

Interaction of a screw dislocation with a semi-infinite interfacial crack in a magneto-electro-elastic bi-material

The interaction between a screw dislocation and a semi-infinite interfacial crack in a transversely isotropic magneto-electro-elastic bi-material is investigated. The dislocation line is perpendicular to the isotropic basal plane of the bi-material. The elastic and electromagnetic fields induced by the dislocation are obtained through the use of the complex variable method together with the superposition scheme. The stress, electric displacement and magnetic intensity factors as well as the image exerted on the dislocation are given explicitly. We find that the intensity factors are expressed in terms of the so-called effective materials and the radial component of the image force is only dependent on the elastic modulus of the material with the dislocation. As an illustrative example, the bi-material that consists of piezoelectric and piezomagnetic phases is analyzed.     

Propagation of a mode-III interfacial crack in a piezoelectric–piezomagnetic bi-material

The transient response of a semi-infinite mode-III interfacial crack propagating between piezoelectric (PE) and piezomagnetic (PM) half spaces is investigated in this paper. The integral transform method together with the Wiener–Hopf and Cagniard–de Hoop techniques is used to solve the mixed boundary value problem under consideration. The existence of generalized Maerfeld–Tournois interfacial wave is discussed and the solutions of the coupled fields are derived for four different cases of bulk shear wave velocity. The dynamic intensity factors of stress, electric displacement and magnetic induction as well as energy release rate (ERR) are obtained in explicit forms. The numerical results of the universal functions and dimensionless ERR for several different material combinations are presented and discussed in details. It is found that the Bleustein–Gulyaev (generalized Maerfeld–Tournois) waves dominate the dynamic characteristics of the interfacial crack propagation in PE–PM bi-material.     

Theoretical model of crack branching in magnetoelectric thermoelastic materials

Thermomagnetoelectroelastic crack branching of magnetoelectro thermoelastic materials is theoretically investigated based on Stroh formalism and continuous distribution of dislocation approach. The crack face boundary condition is assumed to be fully thermally, electrically and magnetically impermeable. Explicit Green’s functions for the interaction of a crack and a thermomagnetoelectroelastic dislocation (i.e., a thermal dislocation, a mechanical dislocation, an electric dipole and a magnetic dipole located at a same point) are presented. The problem is reduced to two sets of coupled singular integral equations with the thermal dislocation and magnetoelectroelastic dislocation densities along the branched crack line as the unknown variables. As a result, the formulations for the stress, electric displacement and magnetic induction intensity factors and energy release rate at the branched crack tip are expressed in terms of the dislocation density functions and the branch angle. Numerical results are presented to study the effect of applied thermal flux, electric field and magnetic field on the crack propagation path by using the maximum energy release rate criterion.     

Analysis method of planar cracks of arbitrary shape in the isotropic plane of a three-dimensional transversely isotropic magnetoelectroelastic medium

The hyper-singular boundary integral equation method of crack analysis in three-dimensional transversely isotropic magnetoelectroelastic media is proposed. Based on the fundamental solutions or Green’s functions of three-dimensional transversely isotropic magnetoelectroelastic media and the corresponding Somigliana identity, the boundary integral equations for a planar crack of arbitrary shape in the plane of isotropy are obtained in terms of the extended displacement discontinuities across crack faces. The extended displacement discontinuities include the displacement discontinuities, the electric potential discontinuity and the magnetic potential discontinuity, and correspondingly the extended tractions on crack face represent the conventional tractions, the electric displacement and the magnetic induction boundary values. The near crack tip fields and the intensity factors in terms of the extended displacement discontinuities are derived by boundary integral equation approach. A solution method is proposed by use of the analogy between the boundary integral equations of the magnetoelectroelastic media and the purely elastic materials. The influence of different electric and magnetic boundary conditions, i.e., electrically and magnetically impermeable and permeable conditions, electrically impermeable and magnetically permeable condition, and electrically permeable and magnetically impermeable condition, on the solutions is studied. The crack opening model is proposed to consider the real crack opening and the electric and magnetic fields in the crack cavity under combined mechanical-electric-magnetic loadings. An iteration approach is presented for the solution of the non-linear model. The exact solution is obtained for the case of uniformly applied loadings on the crack faces. Numerical results for a square crack under different electric and magnetic boundary conditions are displayed to demonstrate the proposed method.     

垂直磁压电材料界面三维裂纹的超奇异积分法

应用有限部积分概念和广义位移基本解,垂直于磁压电双材料界面三维复合型裂纹问题被转 化为求解一组以裂纹表面广义位移间断为未知函数的超奇异积分方程问题. 进而,通过主部 分析法精确地求得裂纹尖端光滑点附近的奇性应力场解析表达式. 然后,通过将裂纹表面 位移间断未知函数表达为位移间断基本密度函数与多项式之积,使用有限部积分法对超奇异 积分方程组建立了数值方法. 最后,通过典型算例计算,讨论了广义应力强度因子的变化规 律.     

模拟裂纹扩展的单位分解扩展无网格法

单位分解扩展无网格法(PUEM)是一种求解不连续问题的新型无网格方法.其基于单位分解思想,通过在传统无网格法的近似函数中加入扩展项来反映由裂纹所产生的不连续位移场.详细描述了水平集方法,PUEM不连续近似函数的构造及控制方程的离散.针对裂纹扩展问题,提出了一种十分简单的水平集更新算法;讨论了不同的节点数、高斯积分阶次以及围线积分区域对应力强度因子计算结果的影响,并给出了合理的参数;模拟了边裂纹和中心裂纹的扩展问题,并与XFEM的数值结果进行了比较.数值算例表明,本文方法具有较高的计算精度,是模拟裂纹扩展非常有效的方法,具有广阔的应用前景.     

Finite-part integral and boundary element method to solve three-dimensional crack problems in piezoelectric materials

Using the hypersingular integral equation method based on body force method, a planar crack in a three-dimensional transversely isotropic piezoelectric solid under mechanical and electrical loads is analyzed. This crack problem is reduced to solve a set of hypersingular integral equations. Compare with the crack problems in elastic isotropic materials, it is shown that for the impermeable crack, the intensity factors for piezoelectric materials can be obtained from those for elastic isotropic materials. Based on the exact analytical solution of the singular stresses and electrical displacements near the crack front, the numerical method of the hypersingular integral equation is proposed by the finite-part integral method and boundary element method, which the square root models of the displacement and electric potential discontinuities in elements near the crack front are applied. Finally, the numerical solutions of the stress and electric field intensity factors of some examples are given.     

INTERFACIAL CRACK ANALYSIS IN THREE-DIMENSIONAL TRANSVERSELY ISOTROPIC BI-MATERIALS BY BOUNDARY INTEGRAL EQUATION METHOD

The integral-differential equations for three-dimensional planar interfacial cracks of arbitrary shape in transversely isotropic bimaterials were derived by virtue of the Somigliana identity and the fundamental solutions, in which the displacement discontinuities across the crack faces are the unknowns to be determined. The interface is parallel to both the planes of isotropy. The singular behaviors of displacement and stress near the crack border were analyzed and the stress singularity indexes were obtained by integral equation method. The stress intensity factors were expressed in terms of the displacement discontinuities. In the non-oscillatory case, the hyper-singular boundary integral-differential equations were reduced to hyper-singular boundary integral equations similar to those of homogeneously isotropic materials.     

Interaction of a conductive crack and of an electrode at a piezoelectric bimaterial interface

The interaction of a conductive crack and an electrode at a piezoelectric bi-material interface is studied. The bimaterial is subjected to an in-plane electrical field parallel to the interface and an anti-plane mechanical loading. The problem is formulated and reduced, via the application of sectionally analytic vector functions, to a combined Dirichlet–Riemann boundary value problem. Simple analytical expressions for the stress, the electric field, and their intensity factors as well as for the crack faces' displacement jump are derived. Our numerical results illustrate the proposed approach and permit to draw some conclusions on the crack–electrode interaction.     

横观各向同性压电材料中裂纹问题的边界元法

采用Somigiliana公式给出了三维横观各向同性压电材料中的非渗漏裂纹问题的一般解和超奇异积分方程,其中未知函数为裂纹面上的位移间断和电势间断．在此基础上,使用有限部积分和边界元结合的方法,建立了超奇异积分方程的数值求解方法,并给出了一些典型数值算例的应力强度因子和电位移强度因子的数值结果,结果令人满意．     

Two dimensional numerical simulation of crack kinking from an interface under dynamic loading by time domain boundary element method

The hybrid time-domain boundary element method, together with the multi-region technique, is applied to simulate the dynamic process of propagation and/or kinking of an interface crack in a two-dimensional bi-material. The whole bi-material is divided into two regions along the interface. The traditional displacement boundary integral equations are employed with respect to each region. However, when the crack kinks into the matrix material, the non-hypersingular traction boundary integral equations 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 controlled by the fracture criterion. Kinking of the interface crack is controlled by a criterion developed from the quasi-static one. Once the crack kinks into the matrix, its propagation is modeled by adding new elements of constant length to the moving crack-tip controlled by a criterion extended from the quasi-static maximum circumferential stress. The numerical results of the crack growth trajectory for different material combinations are computed and compared with the corresponding experimental results. Good agreement between numerical and experimental results implies that the present boundary element numerical method can provide an excellent simulation for the dynamic propagation and deflection of an interface crack.     

含圆形嵌体弹性平面中径向裂纹问题的超奇异积分方程方法

根据含圆形嵌体平面问题在极坐标下的弹性力学基本解,使用Betti互换定理,在有限部积分意义下将问题归结为两个以裂纹岸位移间断为基本未知量、对于Ⅰ型和Ⅱ型问题相互独立的超奇异积分方程,对含圆形嵌体弹性平面中的径向裂纹问题进行了研究.根据有限部积分原理,建立了问题的数值算法.计算结果表明,嵌体半径、裂纹位置及材料剪切弹性模量等都对裂纹应力强度因子具有较为明显的影响.