SINGULAR BEHAVIORS OF INTERFACIAL CRACKS IN 2D MAGNETOELECTROELASTIC BIMATERIALS

The singular characteristics of stress, electric displacement and magnetic induction fields near the tip of impermeable interracial cracks in two-dimensional magnetoelectroelastic bimaterials are studied using the generalized Stroh formalism. Two types of singularities are obtained： one is the oscillating singularity 1/2±iε, the other is the non-oscillating singularity 1/2±κ. It is found that the non-zero parameters ε and κ cannot coexist for one transversely isotropic MEE bimaterial, a similar result is obtained for transversely isotropic piezoelectric bimaterials.     

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.     

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.     

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.     

A unified definition for stress intensity factors of interface corners and cracks

Based upon linear fracture mechanics, it is well known that the singular order of stresses near the crack tip in homogeneous materials is a constant value −1/2, which is nothing to do with the material properties. For the interface cracks between two dissimilar materials, the near tip stresses are oscillatory due to the order of singularity being −1/2 ± iε and −1/2. The oscillation index ε is a constant related to the elastic properties of both materials. While for the general interface corners, their singular orders depend on the corner angle as well as the elastic properties of the materials. Owing to the difference of the singular orders of homogeneous cracks, interface cracks and interface corners, their associated stress intensity factors are usually defined separately and even not compatibly. Since homogenous cracks and interface cracks are just special cases of interface corners, in order to build a direct connection among them a unified definition for their stress intensity factors is proposed in this paper. Based upon the analytical solutions obtained previously for the multibonded anisotropic wedges, the near tip solutions for the general interface corners have been divided into five different categories depending on whether the singular order is distinct or repeated, real or complex. To provide a stable and efficient computing approach for the general mixed-mode stress intensity factors, the path-independent H-integral based on reciprocal theorem of Betti and Rayleigh is established in this paper. The complementary solutions needed for calculation of H-integral are also provided in this paper. To illustrate our results, several different kinds of examples are shown such as cracks in homogenous isotropic or anisotropic materials, central or edge notches in isotropic materials, interface cracks and interface corners between two dissimilar materials.     

Interface crack-tip generalized stress field and stress intensity factors in transversely isotropic piezoelectric bimaterials

Transversely isotropic piezoelectric (TIP) bimaterials with an impermeable interface crack have been classified [Int. J. Frac. 119 (2003) L41] into two classes corresponding to the vanishing of the two singularity parameters or κ. It is shown in the present paper that the related eigenvalue problems for either =0 or κ=0 are not degenerate. The crack-tip generalized stress fields are obtained subsequently. A new definition of crack-tip intensity factors is presented for interface cracks in practical TIP bimaterial of practical interest. Such defined intensity factors are real numbers, which dominate the maximum crack-tip stress singularity and do not generate any phase angle change under any dimension system transformation for physical quantities.     

Image singularities of Green’s functions for isotropic elastic bimaterials subjected to concentrated forces and dislocations

Green’s functions for isotropic materials in the two-dimensional problem for elastic bimaterials with perfectly bonded interface are reexamined in the present study. Although the Green’s function for an isotropic elastic bimaterial subjected to a line force or a line dislocation has been discussed by many authors, the physical meaning and the structure of the solution are not clear. In this investigation, the Green’s function for an elastic bimaterial is shown to consist of eight Green’s functions for a homogeneous infinite plane. One of the novel features is that Green’s functions for bimaterials can be expressed directly by knowing Green’s functions for the infinite plane. If the applied load is located in material 1, the solution for the half-plane of material 1 is constructed with the help of five Green’s functions corresponding to the infinite plane. However, the solution for the half-plane of material 2 only consists of three Green’s functions for the infinite plane. One of the five Green’s functions of material 1 and all the three Green’s functions of material 2 have their singularities located in the half-plane where the load is applied, and the other four image singularities of material 1 are located outside the half-plane at the same distance from the interface as that of the applied load. The nature and magnitude of the image singularities for both materials are presented explicitly from the principle of superposition, and classified according to different loads. It is known that for the problem of anisotropic bimaterials subjected to concentrated forces and dislocations, the image singularities are simply concentrated forces and dislocations with the stress singularity of order O(1/r). However, higher orders (O(1/r²) and O(1/r³)) of stress singularities are found to exist in this study for isotropic bimaterials. The highest order of the stress singularity is O(1/r³) for the image singularities of material 1, and is O(1/r²) for material 2. Using the present solution, Green’s functions associated with the problems of elastic half-plane with free and rigidly fixed boundaries, for homogeneous isotropic elastic solid, are obtained as special cases.     

On the crack-tip stress singularity of interfacial cracks in transversely isotropic piezoelectric bimaterials

In this paper, characteristics of the interface crack-tip stress and electric displacement fields in transversely isotropic piezoelectric bimaterials are studied. The authors have proven, within the framework of the generalized Stroh formalism for piezoelectric bimaterials, that there is no coexistence of the parameters (oscillating) and κ (non-oscillating) in the interface crack-tip generalized stress field for all transversely isotropic piezoelectric bimaterials. This leads to the classification of piezoelectric bimaterials into one group that exhibits the oscillating property in the interface crack-tip generalized stress field and the other that does not. Fifteen (15) pair-combinations of six (6) piezoelectric materials PZT-4, PZT-5H, PZT-6B, PZT-7A, P-7, and BaTiO₃, which are commonly used in practice, are numerically analyzed in this study, and the results backup the above theoretical conclusions. Moreover, the associated eigenvectors for such material systems (with either =0 or κ=0) are also obtained numerically, and the result show that there still exist four linear independent associate eigenvectors for each bimaterial.     

Controlling parameter of the stress intensity factors for a planar interfacial crack in three-dimensional bimaterials

Numerical solutions of singular integral equations are discussed in the analysis of a planar rectangular interfacial crack in three-dimensional bimaterials subjected to tension. The problem is formulated as a system of singular integral equations on the basis of the body force method. In the numerical analysis, unknown body force densities are approximated by the products of the fundamental density functions and power series, where the fundamental density functions are chosen to express singular behavior along the crack front of the interface crack exactly. The calculation shows that the present method gives smooth variations of stress intensity factors along the crack front for various aspect ratios. The present method gives rapidly converging numerical results and highly satisfied boundary conditions throughout the crack boundary. The stress intensity factors are given with varying the material combination and aspect ratio of the crack. It is found that the stress intensity factors K_I and K_II are determined by the bimaterial constant ε alone, independent of elastic modulus ratio and Poisson’s ratio.     

双材料界面裂纹的改进广义有限差分法

采用数值方法进行断裂力学分析时,裂纹尖端奇异区域处理的好坏直接关系到最终断裂力学参数的求解精度。与传统均匀介质不同,复合材料界面裂纹渐近位移和应力场表现出剧烈的振荡特性,许多用于表征经典的平方根和负平方根物理场渐近性的传统方法也因此失效。论文提出了一种改进的广义有限差分法,该方法基于多元函数泰勒级数展开和移动最小二乘法的思想,将节点变量的各阶导数由相邻点集函数的加权线性累加来近似,具有无网格、无数值积分、数据准备简单、稀疏矩阵快速求解等优点。为提高该方法求解断裂力学问题的计算精度和数值稳定性,论文引入了裂尖奇异区域局部点簇的自动创建技术和一种基于局部点簇几何尺寸的矩阵正则化算法。数值算例表明,所提算法稳定,效率高,在不增加计算量的前提下,显著提高了裂尖近场力学参量和断裂力学参数的求解精度和数值稳定性。     

双材料界面裂纹复应力强度因子的正则化边界元法

双材料界面裂纹渐近位移和应力场表现出剧烈的振荡特性,许多用于表征经典平方根(r1/2)和负平方根(r?1/2)渐近物理场的传统数值方法失效,给界面裂纹复应力强度因子(K1+iK2)的精确求解增加了难度.引入一种含有复振荡因子的新型"特殊裂尖单元",可精确表征裂纹尖端渐近位移和应力场的振荡特性,在避免裂尖区域高密度网格剖...     

Two-dimensional Green's functions in anisotropic multiferroic bimaterials with a viscous interface

We derive, by virtue of the unified Stroh formalism, the extremely concise and elegant solutions for two-dimensional and (quasi-static) time-dependent Green's functions in anisotropic magnetoelectroelastic multiferroic bimaterials with a viscous interface subjected to an extended line force and an extended line dislocation located in the upper half-plane. It is found for the first time that, in the multiferroic bimaterial Green's functions, there are 25 static image singularities and 50 moving image singularities in the form of the extended line force and extended line dislocation in the upper or lower half-plane. It is further observed that, as time evolves, the moving image singularities, which originate from the locations of the static image singularities, will move further away from the viscous interface with explicit time-dependent locations. Moreover, explicit expression of the time-dependent image force on the extended line dislocation due to its interaction with the viscous interface is derived, which is also valid for mathematically degenerate materials. Several special cases are discussed in detail for the image force expression to illustrate the influence of the viscous interface on the mobility of the extended line dislocation, and various interesting features are observed. These Green's functions can not only be directly applied to the study of dislocation mobility in the novel multiferroic bimaterials, they can also be utilized as kernel functions in a boundary integral formulation to investigate more complicated boundary value problems where multiferroic materials/composites are involved.     

Magnetoelectroelastic media containing a row of moving shear cracks

The elastic, electric and magnetic fields created within transversely isotropic magnetoelectroelastic media around an infinite row of uniformly-moving, collinear, antiplane shear cracks are studied, using an extension of the powerful method of dislocation layers. This analysis additionally provides the solutions for a finite magnetoelectroelastic plate containing a single crack and a plate with an edge crack.     

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

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

Anti-plane shear crack normal to and terminating at the interface of two bonded piezoelectric ceramics

The electroelastic analysis of two bonded dissimilar piezoelectric ceramics with a crack perpendicular to and terminating at the interface is made. By using Fourier integral transform, the associated boundary value problem is reduced to a singular integral equation with generalized Cauchy kernel, the solution of which is given in closed form. Results are presented for a permeable crack under anti-plane shear loading and in-plane electric loading. Obtained results indicate that the electroelastic field near the crack tip in the homogeneous piezoelectric ceramic is dominated by a traditional inverse square-root singularity, while the electroelastic field near the crack tip at the interface exhibits the singularity of power law r^−α, r being distance from the interface crack tip and α depending on the material constants of a bi-piezoceramic. In particular, electric field has no singularity at the crack tip in a homogeneous solid, whereas it is singular around the interface crack tip. Numerical results are given graphically to show the effects of the material properties on the singularity order and field intensity factors.     

Line-integral representations for extended displacements,stresses, and interaction energy of arbitrary dislocation loops in transversely isotropic magneto-electro-elastic bimaterials

In addition to the hexagonal crystals of class 6 mm, many piezoelectric materials （e.g., BaTiO3）, piezomagnetic materials （e.g., CoFe2O4）, and multiferroic com-posite materials （e.g., BaTiO3-CoFe2O4 composites） also exhibit symmetry of transverse isotropy after poling, with the isotropic plane perpendicular to the poling direction. In this paper, simple and elegant line-integral expressions are derived for extended displace-ments, extended stresses, self-energy, and interaction energy of arbitrarily shaped, three-dimensional （3D） dislocation loops with a constant extended Burgers vector in trans-versely isotropic magneto-electro-elastic （MEE） bimaterials （i.e., joined half-spaces）. The derived solutions can also be simply reduced to those expressions for piezoelectric, piezo-magnetic, or purely elastic materials. Several numerical examples are given to show both the multi-field coupling effect and the interface/surface effect in transversely isotropic MEE materials.     

Line-integral representations for the elastic displacements,stresses and interaction energy of arbitrary dislocation loops in transversely isotropic bimaterials

The elastic displacements, stresses and interaction energy of arbitrarily shaped dislocation loops with general Burgers vectors in transversely isotropic bimaterials (i.e. joined half-spaces) are expressed in terms of simple line integrals for the first time. These expressions are very similar to their isotropic full-space counterparts in the literature and can be easily incorporated into three-dimensional (3D) dislocation dynamics (DD) simulations for hexagonal crystals with interfaces/surfaces. All possible degenerate cases, e.g. isotropic bimaterials and isotropic half-space, are considered in detail. The singularities intrinsic to the classical continuum theory of dislocations are removed by spreading the Burgers vector anisotropically around every point on the dislocation line according to three particular spreading functions. This non-singular treatment guarantees the equivalence among different versions of the energy formulae and their consistency with the stress formula presented in this paper. Several numerical examples are provided as verification of the derived dislocation solutions, which further show significant influence of material anisotropy and bimaterial interface on the elastic fields and interaction energy of dislocation loops.     

Finite element analysis of a bimaterial interface crack

In this investigation, the enriched element method developed by Benzley was extended to treat the stress analysis problem involving a bimaterial interface crack. Unlike crack problems in isotropic elasticity, where the stress singularity at the crack tip is of the inverse square root type, the interface crack contains an additional oscillatory singularity. Although the effect of this oscillatory characteristic is confined to a region very close to the crak tip, it nevertheless requires proper treatment in order to obtain accurate predictions on the stress intensity factors. Using appropriate crack tip stress and displacement expressions, the enriched element method can model the stress singularity for an interface crack exactly. The finite element implementation of this method has been made on the code APES. Stress intensity factor results predicted by the modified APES program compare favorably with those available in the literature. This indicates tha the enriched element technique provides an accurate and efficient numerical tool for the analysis of bimaterial interface crack problems.