A piezoelectric-elastic body with inhomogeneities and a crack under plane electrical and anti-plane mechanical loads

Summary  In this paper, we study a two-dimensional electroelastic problem of an infinite piezoelectric body with two circular piezoelectric inhomogeneities, one of which contains a crack. We formulate the stress intensity factor (SIF) analytically and investigate it numerically. The problem is solved based on Bueckner's principle, and is reduced to a problem of a singular integral equation of the first kind with respect to the distribution function of screw dislocation. The effect of interaction between the two inhomogeneities and the crack on the electroelastic field as well as the control of the SIF by electrical loads is investigated. Received 18 April 2000; accepted for publication 24 October 2000     

Anti-plane stress intensity, energy release and energy density at crack tips in a functionally graded strip with linearly varying properties

The fracture problem of a crack in a functionally graded strip with its properties varying in a linear form along the strip thickness under an anti-plane load is considered. The embedded anti-plane crack is located in the middle of strip half way through the thickness. The third mode stress intensity factor is derived using two different methods. In the first method, by employing Fourier integral transforms, the governing equation is converted to a singular integral equation, which is subsequently solved numerically by the collocation method based on Chebyshev polynomials. Then, the problem is solved by means of finite element method in which quadrilateral 8-node singular elements around each crack tip are used. After inspecting the validity of the solution technique, effects of crack geometry and non-homogeneous material parameter on the stress intensity, energy release and energy density are studied and the results of analytical and FEM solutions are compared.     

Griffith crack moving in a piezoelectric strip

Summary  The dynamic problem of an impermeable crack of constant length 2a propagating along a piezoelectric ceramic strip is considered under the action of uniform anti-plane shear stress and uniform electric field. The integral transform technique is employed to reduce the mixed-boundary-value problem to a singular integral equation. For the case of a crack moving in the mid-plane, explicit analytic expressions for the electroelastic field and the field intensity factors are obtained, while for an eccentric crack moving along a piezoelectric strip, numerical results are determined via the Lobatto–Chebyshev collocation method for solving a resulting singular integral equation. The results reveal that the electric-displacement intensity factor is independent of the crack velocity, while other field intensity factors depend on the crack velocity when referred to the moving coordinate system. If the crack velocity vanishes, the present results reduce to those for a stationary crack in a piezoelectric strip. In contrast to the results for a stationary crack, applied stress gives rise to a singular electric field and applied electric field results in a singular stress for a moving crack in a piezoelectric strip. Received 14 August 2001; accepted for publication 24 September 2002 The author is indebted to the AAM Reviewers for their helpful suggestions for improving this paper. The work was supported by the National Natural Science Foundation of China under Grant 70272043.     

黏弹性材料界面裂纹应力场奇性分析

研究两半无限大黏弹性体间Griffith界面裂纹在简谐载荷作用下裂纹尖端动应力场的奇异特性.通过引入裂纹张开位移和裂纹位错密度函数,相应的混合边值问题归结为一组耦合的奇异积分方程.渐近分析表明裂尖动应力场的奇异特征完全包含在奇异积分方程的基本解中.通过对基本解的深入分析发现黏弹性材料界面裂纹裂尖动应力场具有与材料参数和外载荷频率相关的振荡奇异特性.以标准线性固体黏弹材料为例讨论了材料参数和载荷频率对奇性指数和振荡指数的影响.     

A long circular cylinder with a circumferential edge crack subjected to a uniform shearing stress

This paper contains an analysis of the stress distribution in a long circular cylinder of isotropic elastic material with a circumferential edge crack when it is deformed by the application of a uniform shearing stress. The crack with its center on the axis of the cylinder lies on the plane perpendicular to that axis, and the cylindrical surface is stress-free. By making a suitable representation of the stress function for the problem, the problem is reduced to the solution of a pair of singular integral equations. This pair of singular integral equations is solved numerically, and the stress intensity factor due to the effect of the crack size is tabulated.     

Interaction of general plane P wave and cylindrical inclusion partially debonded from its viscoelastic matrix

The interaction of a general plane P wave and an elastic cylindrical inclusion of infinite length partially debonded from its surrounding viscoelastic matrix of infinite extension is investigated. The debonded region is modeled as an arc-shaped interface crack between inclusion and matrix with non-contacting faces. With wave functions expansion and singular integral equation technique, the interaction problem is reduced to a set of simultaneous singular integral equations of crack dislocation density function. By analysis of the fundamental solution of the singular integral equation, it is found that dynamic stress field at the crack tip is oscillatory singular, which is related to the frequency of incident wave. The singular integral equations are solved numerically, and the crack open displacement and dynamic stress intensity factor are evaluated for various incident angles and frequencies. The project supported by the National Natural Science Foundation of China (19872002) and Climbing Foundation of Northern Jiaotong University     

Closed-form solution for two collinear mode-III cracks in an orthotropic elastic strip of finite width

The problem of an orthotropic strip containing two collinear cracks normal to the strip boundaries is considered. The Fourier series method is used to reduce the associated boundary value problem to triple series equations, then to a singular integral equation, which can be solved analytically. Under remote uniform antiplane shear loading, the stress field and the crack sliding displacement are determined analytically and stress intensity factors are also given in a closed form.     

An Interface Crack for a Functionally Graded Strip Sandwiched Between Two Homogeneous Layers of Finite Thickness

In this paper, the behavior of an interface crack for a functionally graded strip sandwiched between two homogeneous layers of finite thickness subjected to an uniform tension is resolved using a somewhat different approach, named the Schmidt method. The Fourier transform technique is applied and a mixed boundary value problem is reduced to two pairs of dual integral equations in which the unknown variables are the jumps of the displacements across the crack surface. To solve the dual integral equations, the jumps of the displacements across the crack surfaces are expanded in a series of Jacobi polynomials. This process is quite different from those adopted in previous works. Numerical examples are provided to show the effects of the crack length, the thickness of the material layer and the materials constants upon the stress intensity factor of the cracks. It can be obtained that the results of the present paper are the same as ones of the same problem that was solved by the singular integral equation method. As a special case, when the material properties are not continuous through the crack line, an approximate solution of the interface crack problem is also given under the assumption that the effect of the crack surface interference very near the crack tips is negligible. Contrary to the previous solution of the interface crack, it is found that the stress singularities of the present interface crack solution are the same as ones of the ordinary crack in homogenous materials.     

Variation of the stress intensity factor along the front of a 3-D rectangular crack subjected to mixed-mode load

Summary  The singular integral equation method is applied to the calculation of the stress intensity factor at the front of a rectangular crack subjected to mixed-mode load. The stress field induced by a body force doublet is used as a fundamental solution. The problem is formulated as a system of integral equations with r ⁻³-singularities. In solving the integral equations, unknown functions of body-force densities are approximated by the product of polynomial and fundamental densities. The fundamental densities are chosen to express two-dimensional cracks in an infinite body for the limiting cases of the aspect ratio of the rectangle. The present method yields rapidly converging numerical results and satisfies boundary conditions all over the crack boundary. A smooth distribution of the stress intensity factor along the crack front is presented for various crack shapes and different Poisson's ratio. Received 5 March 2002; accepted for publication 2 July 2002     

曲线裂纹和反平面圆形夹杂相交问题

建立了和反平面圆夹杂界面相交的曲线裂纹的弱奇异积分方程,利用Cauchy型奇异积分方程主部分析方法研究了穿过反平面圆夹杂界面的曲线裂纹在交点处的奇性应力指数以及交点处角形域内的奇性应力,并根据奇性应力定义了交点处的应力强度因子。通过对弱奇异积分方程的数值求解,可得裂纹端点和交点处的应力强度因子。     

Diffraction of torsional elastic waves by a peripheral edge crack around a spherical cavity

This paper deals with the problem of finding the stress distribution in the neighborhood of a peripheral edge crack in a spherical cavity. The crack is excited by a torsional standing wave.The problem is solved by using integral transforms and is reduced to the solution of a singular integral equation. The solution of this equation is obtained numerically by the method due to Erdogan, Gupta, and Cook, and the stress intensity factors are displayed graphically.     

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).     

Strip yield zones around ring-shaped crack in transversely isotropic thick layer

A ring-shaped crack under uniform load in an infinitely long elastic–perfectly plastic thick layer is considered. The problem is formulated for a transversely isotropic material by using integral transform technique. Due to the geometry of the configuration, Hankel integral transform technique was chosen and the problem was reduced to a singular integral equation which is solved numerically by using Gaussian Quadrature Formulae and the values were evaluated at discrete points. The plastic zone widths were obtained by using the plastic strip model after stress intensity factors were obtained. Numerical results are plotted for various ring-shaped crack sizes and transversely isotropic materials. It was found that the width of the plastic zone at the inner edge of the crack was greater than the outer one.     

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.     

Anti-plane fracture of a functionally graded material strip

This paper considers the anti-plane (or mode III) crack problem in a functionally graded material strip. The shear modulus of the strip is considered for a class of functional forms for which the equilibrium equation has an analytical solution. The problem is solved by means of singular integral equation technique. Both a single crack and a series of collinear cracks are studied. The results are tabulated and plotted to show the effect of the material nonhomogeneity and crack location on the stress intensity factors.     

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     

Interaction between curved crack and elastic inclusion in an infinite plate

Summary In this paper, the curved-crack problem for an infinite plate containing an elastic inclusion is considered. A fundamental solution is proposed, which corresponds to the stress field caused by a point dislocation in an infinite plate containing an elastic inclusion. By placing the distributed dislocation along the prospective site of the crack, and by using the resultant force function as the right-hand term in the equation, a weaker singular integral equation is obtainable. The equation is solved numerically, and the stress intensity factors at the crack tips are evaluated. Interaction between the curved crack and the elastic inclusion is analyzed. Received 8 October 1996; accepted for publication 27 March 1997     

半平面多边缘裂纹反平面问题的奇异积分方程

利用复变函数和奇异积分方程方法,求解弹性范围内半平面多边缘裂纹的反平面问题．提出了满足半平面边界自由的由分布位错密度表示的单边缘裂纹的基本解,此基本解由主要部分和辅助部分组成．将半平面多边缘裂纹问题看作是许多单边缘裂纹问题的叠加,建立了一组Cauchy型奇异积分方程．然后,利用半开型积分法则求解该奇异积分方程,得到了裂纹端处的应力强度因子．最后,给出了几个数值算例．     

Interaction between crack and elastic inclusion

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

Eccentric crack in a piezoelectric strip bonded to half planes

We consider the problem of determining the singular stresses and electric fields in a piezoelectric ceramic strip containing an eccentric Griffith crack off the centre line bonded to two elastic half planes under anti-plane shear loading using the continuous crack-face condition. Fourier transforms are used to reduce the problem to the solution of two pairs of dual integral equations, which are then expressed to a Fredholm integral equation of the second kind. Numerical values on the stress intensity factor and energy release rate are obtained.