The scattering of general SH plane wave by interface crack between two dissimilar viscoelastic bodies

The scattering of general SH plane wave by an interface crack between two dissimilar viscoelastic bodies is studied and the dynamic stress intensity factor at the crack-tip is computed. The scattering problem can be decomposed into two problems: one is the reflection and refraction problem of general SH plane waves at perfect interface (with no crack); another is the scattering problem due to the existence of crack. For the first problem, the viscoelastic wave equation, displacement and stress continuity conditions across the interface are used to obtain the shear stress distribution at the interface. For the second problem, the integral transformation method is used to reduce the scattering problem into dual integral equations. Then, the dual integral equations are transformed into the Cauchy singular integral equation of first kind by introduction of the crack dislocation density function. Finally, the singular integral equation is solved by Kurtz's piecewise continuous function method. As a consequence, the crack opening displacement and dynamic stress intensity factor are obtained. At the end of the paper, a numerical example is given. The effects of incident angle, incident frequency and viscoelastic material parameters are analyzed. It is found that there is a frequency region for viscoelastic material within which the viscoelastic effects cannot be ignored. This work was supported by the National Natural Science Foundation of China (No.19772064) and by the project of CAS KJ 951-1-20     

平面波对有部分脱胶衬砌的圆形孔洞的散射

将波函数展开法与奇异积分方程技术相结合研究了平面波对有部分脱胶衬砌的圆形孔洞的散射。将脱离区看作弧形裂纹并忽略裂纹表面的相互作用。将衬砌和基体中的波场展开成Fou-rier-Bessel级数,利用混合边界条件得到一组对偶级数方程组并进一步转化成Hilbert奇异积分方程。数值求解给出了脱离区大小和衬砌厚度对动应力强度因子(DSIF)和散射截面(SCS)的影响。结果显示由于脱胶,动应力强度因子和散射截面呈现明显的低频共振特性。     

P波对界面部分脱胶的刚性圆柱夹杂物的散射

本文求解了弹性Ｐ波对界面部分脱胶的可动刚性圆柱夹杂物的散射问题。将脱胶区看作表面不相接触的弧形界面裂纹，借助波函数展开法并利用边界条件将问题转化为一组对偶级数方程。然后通过引入裂纹面的位错密度函数，将其化为一组具有Ｈｉｌｂｅｒｔ核的第二类奇异积分方程，并进一步化为Ｃａｕｃｈｙ型奇异积分方程组，数值求解方程组可获得动应力强度因子，夹杂物刚体振动位移和散射截面等重要参量。结果显示该类结构在较低的频率上发生共振，此低频共振特性与脱胶区大小，入射波方向、材料组合等多种参数有关。与已有方法相比，本文的方法更具一般性，适用于任意材料组合。     

SHEAR WAVE SCATTERING FROM A PARTIALLY DEBONDED PIEZOELECTRIC CYLINDRICAL INCLUSION

I. INTRODUCTION Owing to the intrinsic coupling characteristics between electric and elastic behaviors, piezoelectricmaterials have been used widely in technology such as transducers, actuators, sensors, etc. Studieson electroelastic problems of a piezo…     

Elastic wave scattering from a penny-shaped interface crack in coated materials

This paper is concerned with the elastic wave scattering induced by a penny-shaped interface crack in coated materials. Using the integral transform, the problem of wave scattering is reduced to a set of singular integral equations in matrix form. The singular integral equations are solved by the asymptotic analysis and contour integral technique, and the expressions for the stress and displacement as well as the dynamic stress intensity factors (SIFs) are obtained. Using numerical analysis, this approach is verified by the finite element method (FEM), and the numerical results agree well with the theoretical results. For various crack sizes and material combinations, the relations between the SIFs and the incident frequency are analyzed, and the amplitudes of the crack opening displacements (CODs) are plotted versus incident wavenumber. The investigation provides a theoretical basis for the dynamic failure analysis and nondestructive evaluation of coated materials.     

Interaction between a compliant disk-shaped inclusion and a crack upon incidence of an elastic wave

The propagation of harmonic elastic wave in an infinite three-dimensional matrix containing an interacting low-rigidity disk-shaped inclusion and a crack. The problem is reduced to a system of boundary integral equations for functions that characterize jumps of displacements on the inclusion and crack. The unknown functions are determined by numerical solution of the system of boundary integral equations. For the symmetric problem, graphs are given of the dynamic stress intensity factors in the vicinity of the circular inclusion and the crack on the wavenumber for different distances between them and different compliance parameters of the inclusion.     

Scattering of SHPlease define abbreviation SH waves by arc-shaped interface cracks between a cylindrical magneto-electro-elastic inclusion and matrix: near fields

Summary In this paper, the scattering of SH waves by a magneto-electro-elastic cylindrical inclusion partially debonded from its surrounding magneto-electro-elastic material is investigated by using the wavefunction expansion method and a singular integral equation technique. The debonding regions are modeled as multiple arc-shaped interface cracks with non-contacting faces. The magneto-electric impermeable boundary conditions are adopted. By expressing the scattered fields as wavefunction expansions with unknown coefficients, the mixed boundary-value problem is firstly reduced to a set of simultaneous dual-series equations. Then, dislocation density functions are introduced as unknowns to transform these dual-series equations to Cauchy singular integral equations of the first type,which can be numerically solved easily. The solution is valid for arbitrary number and size of the arc-shaped interface cracks. Finally, numerical results of the dynamic stress intensity factors are presented for the cases of one debond. The effects of incident direction, crack configuration and various material parameters on the dynamic stress intensity factors are discussed. The solution of this problem is expected to have applications in the investigation of dynamic fracture properties of magneto-electro-elastic materials with cracks.The work was supported by the National Natural Science Fund of China (Project No. 19772029) and the Research Fund for Doctors of Hebei Province, China (Project No. B2001213).     

Scattering of love waves by an interface crack between a piezoelectric layer and an elastic substrate

The scattering of Love waves by an interface crack between a piezoelectric layer and an elastic substrate is investigated by using the integral transform and singular integral equation techniques. The dynamic stress intensity factors of the left and the right crack tips are determined. It is found from numerical calculation that the dynamic response of the system depends significantly on the crack configuration, the material combination and the propagating direction of the incident wave. It is expected that specifying an appropriate material combination may retard the growth of the crack for a certain crack configuration. Project supported by the National Natural Science Foundation of China (No. 19891180), the Fundamental Research Foundation of Tsinghua University (JZ 2000.007) and the Fund of the Education Ministry of China.     

Scattering of harmonic waves by a disk-shaped rigid inclusion in a three-dimensional elastic matrix

We consider a three-dimensional problem on the interaction of harmonic waves with a thin rigid movable inclusion in an infinite elastic body. The problem is reduced to solving a system of two-dimensional boundary integral equations of Helmholtz potential type for the stress jump functions on the opposite surfaces of the inclusion. We propose a boundary element method for solving the integral equations on the basis of the regularization of their weakly singular kernels. Using the asymptotic relations between the amplitude-frequency characteristics of the wave farzone field and the obtained boundary stress jump functions, we determine the amplitudes of the shear plane wave scattering by a circular disk-shaped inclusion for various directions of the wave incident on the inclusion and for a broad range of wave numbers.     

薄膜涂层材料中币形界面裂纹的弹性波散射

研究了薄膜涂层材料中币形界面裂纹的弹性波散射问题,建立了含有币形界面裂纹的覆层半空间模型,采用Hankel积分变换,将裂纹对弹性波散射的问题转化为求解矩阵形式的奇异积分方程。结合渐近分析和围道积分技术得到积分方程的解,进一步推导了散射波的应力场和位移场,以及动应力强度因子的理论计算公式。在数值算例中,分析了不同材料组合和裂纹尺寸情况下动应力强度因子与入射波频率的关系,并给出了裂纹张开位移的结果。为薄膜涂层材料的动态破坏分析提供了一定的理论基础。     

Interaction between a circular inclusion and a symmetrically branched crack

The interaction problem between a circular inclusion and a symmetrically branched crack embedded in an infinite elastic medium is solved. The branched crack is modeled as three straight cracks which intersect at a common point and each crack is treated as a continuous contribution of edge dislocations. Green's functions are used to reduce the problem into a system of singular equations consisting of the distributions of Burger's dislocation vectors as unknown functions through the superposition technique. The resulting integral equations are solved numerically by the method of Gauss-Chebychev quadrature. The proposed integral equation approach is first verified for two limiting cases against the literature. More effort is paid on the effect of inclusion on both the Mode I and Mode lI stress intensity factors at the branch tips. The effect of inclusion on the branching path is also investigated.     

Interaction between rigid-disc inclusion and penny-shaped crack under elastic time-harmonic wave incidence

Influence of a rigid-disc massive inclusion on a neighboring penny-shaped crack induced by the time-harmonic wave propagation in an infinite elastic matrix is investigated by the numerical solution of associated 3D elastodynamic problem. No restrictions on the mutual orientation of interacting objects and direction of wave incidence are assumed. The inclusion is perfectly bonded with a matrix and supposes the translations and rotations, the crack faces are load-free. Frequency-domain problem is reduced to a system of boundary integral equations (BIEs) relative to the interfacial stress jumps (ISJs) on the inclusion and the crack opening displacements (CODs). The subtraction technique in conjunction with mapping technique, under taking into account the structure of solution at the fronts of inclusion and crack, is applied for regularization of BIEs obtained. A discrete analogue of equations is constructed by using the collocation scheme. Numerical calculations are carried out for the grazing incidence of a plane P-wave on the crack, where the interacting inclusion is coplanar and perpendicular to the crack, and has the same radius. The shielding and amplification effects of inclusion are assessed by the analysis of mode-I stress intensity factor (SIF) in the crack vicinity depending on the wave number, incident wave direction, position of the crack front point, inclusion mass, crack-inclusion orientation and distance.     

Interaction between crack and elastic inclusion

ＩＮＴＥＲＡＣＴＩＯＮＢＥＴＷＥＥＮＣＲＡＣＫＡＮＤＥＬＡＳＴＩＣＩＮＣＬＵＳＩＯＮＺｈａｎｇＭｉｎｇ－ｈｕａｎ（张明焕）,ＴａｎｇＲｅｎ－ｊｉ（汤任基）（ＳｈａｎｇｈａｉＪｉａｏｔｏｎｇＵｎｉｖｅｒｓｉｔｙ,Ｓｈａｎｇｈａｉ,２０００３０,Ｐ．Ｒ．...     

Stress singularities for three-dimensional corners using the boundary integral equation method

The boundary integral equation method is developed to study three-dimensional asymptotic singular stress fields at vertices of a pyramidal notch or inclusion in an isotropic elastic space. Two-dimensional boundary integral equations are used for the infinite body with pyramidal notches and inclusions when either stresses or displacements are specified on its surface. Applying the Mellin integral transformation reduces the problem to one-dimensional singular integral equations over a closed, piece-wise smooth line. Using quadrature formulas for regular and singular integrals with Hilbert and logarithmic kernels, these integral equations are reduced to a homogeneous system of linear algebraic equations. Setting its determinant to zero provides a characteristic equation for the determination of the stress singularity power. Numerical results are obtained and compared against known eigenvalues from the literature for an infinite region with a conical notch or inclusion, for a Fichera vertex, and for a half-space with a wedge-shaped notch or inclusion.     

Dynamic interaction between a matrix crack and a circular inhomogeneity with a distinct interphase

This article provides a theoretical treatment of the dynamic interaction between a matrix crack and an arbitrarily located circular inhomogeneity with a distinct interphase under antiplane loading. The matrix⧹inhomogeneity interphase is characterized by a linear spring model. The theoretical formulations governing the steady state problem are based upon the use of integral transform techniques, Bessel function expansions and a Pseudo-incident wave technique. The closed form expression for the resulting stress intensity factor at the matrix crack is obtained by solving the appropriate singular integral equations using Chebyshev polynomials. Typical examples are provided to show the effect of the location of the inhomogeneity, the material combination and the interface property upon the dynamic stress intensity factor of the matrix crack.     

ANTIPLANE CIRCULAR INCLUSION WITH A CURVED CRACK CROSSING THE BOUNDARY

The weakly singular integral equation used to solve the problem of the curved crack crossing the boundary of the antiplane circular inclusion is presented. Using the principal part analysis method of the Cauchy type integral equation, the singular stress index at the intersection and the singular stress of angular regions near the intersection are obtained. By using the singular stress obtained, the stress intensity factor at the intersection is, defined. After the numerical solution of the integral equation, the stress intensity factors at the end points of the crack and intersection are obtainable. The research is supported by National Natural Science Foundation of China (No. 59879012) and is the project of Chinese Foundation of State Education Commission (No. 98024832).     

A new integral equation formulation of two-dimensional inclusion–crack problems

A new integral equation formulation of two-dimensional infinite isotropic medium (matrix) with various inclusions and cracks is presented in this paper. The proposed integral formulation only contains the unknown displacements on the inclusion–matrix interfaces and the discontinuous displacements over the cracks. In order to solve the inclusion–crack problems, the displacement integral equation is used when the source points are acting on the inclusion–matrix interfaces, whilst the stress integral equation is adopted when the source points are being on the crack surfaces. Thus, the resulting system of equations can be formulated so that the displacements on the inclusion–matrix interfaces and the discontinuous displacements over the cracks can be obtained. Based on one point formulation, the stress intensity factors at the crack tips can be achieved. Numerical results from the present method are in excellent agreement with those from the conventional boundary element method.     

Analysis of the Singularity for the Vertical Contact Problem of Line Crack-Inclusion

ＩｎｔｒｏｄｕｃｔｉｏｎＵｐｔｏｎｏｗ ,ｔｈｅｔｅｃｈｎｉｃａｌｌｉｔｅｒａｔｕｒｅｏｎｓｅｐａｒａｔｅｃｒａｃｋｓ,ｖｏｉｄｓ,ｉｎｃｌｕｓｉｏｎｓａｎｄｔｈｅｉｎｔｅｒａｃｔｉｏｎｓｂｅｔｗｅｅｎｃｒａｃｋｓａｎｄｉｎｃｌｕｓｉｏｎｓｈａｖｅｂｅｅｎｑｕｉｔｅｅｘｔｅｎｓｉｖｅ.Ｈｏｗｅｖｅｒ,ｔｈｅｃｏｎｔａｃｔｐｒｏｂｌｅｍｓｏｆｃｒａｃｋ＿ｉｎｃｌｕｓｉｏｎｄｏｎｏｔｓｅｅｍｔｏｂｅａｓｗｉｄｅｌｙｓｔｕｄｉｅｄ .Ｔｈｉｓｐａｐｅｒｃａｎｂｅｒｅｇａｒｄｅｄａｓｔｈｅｆｕｒｔｈｅｒｒｅｓｅ…     

Stress-intensity factors for 3-D dynamic loading of a cracked half-space

A half-space containing a surface-breaking crack of uniform depth is subjected to three-dimensional dynamic loading. The elastodynamic stress-analysis problem has been decomposed into two problems, which are symmetric and antisymmetric, respectively, relative to the plane of the crack. The formulation of each problem has been reduced to a system of singular integral equations of the first kind. The symmetric problem is governed by a single integral equation for the opening-mode dislocation density. A pair of coupled integral equations for the two sliding-mode dislocation densities govern the antisymmetric problem. The systems of integral equations are solved numerically. The stress-intensity factors are obtained directly from the dislocation densities. The formulation is valid for arbitrary 3-D loading of the half-space. As an example, an applied stress field corresponding to an incident Rayleigh surface wave has been considered. The dependence of the stress-intensity factors on the frequency, and on the angle of incidence, is displayed in a set of figures.     

三维裂纹分析的主值型面力边界积分方程法

给出了一组只包含Ｃａｕｃｈｙ主值积分、不含有强奇异积分的三维静动力边界积分方程及其应用于裂纹问题的具体列式，并给出了几何轴对称问题的相应半解析边界元求解方法，将三维问题降阶为一维数值问题．文中分析了无限、半无限介质中圆裂纹、平行圆裂纹系、球面裂纹等在静载及应力波作用下的静力或瞬态动力响应问题，求得了相应的应力强度因子．