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.