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.     

Axially symmetric thermal stress of an external circular crack under general thermal conditions

Summary The paper presents a solution for the linear thermoelastic problem of determining axisymmetric stress and displacement fields in an isotropic elastic solid of infinite extent weakened by an external circular crack under general mechanical loadings and general thermal conditions. The mechanical loadings and thermal conditions applied on the crack faces are axisymmetric, being non-symmetric about the crack plane. In similar lines of [7], equations of equilibrium of an elastic solid conducting heat have been solved using Hankel transforms and Abel operators of the first kind. Expressions for stress, displacement, temperature and heat flux functions are obtained in terms of Abel transforms of the first kind of the jumps of stress, displacement, temperature and heat flux at the crack plane. Two types of thermal conditions, that is, general surface temperatures and general heat flux on faces of the crack are considered. In both the cases, closed form solutions have been obtained for the unknown functions solving Abel type of integral equations. Explicit expressions for stresses, displacements, temperature fields, stress intensity factors have been obtained. Two special cases of thermal conditions in which: (i) crack faces are subjected to constant non-symmetric temperatures over a circular ring area, (ii) crack faces are subjected to constant non-symmetric heat flux over a circular ring area, have been considered. In some special cases, results have been compared with those from the literature.     

FRACTURE ANALYSIS OF AN ANNULAR CRACK IN A PIEZOELECTRIC LAYER

A flat annular crack in a piezoelectric layer subjected to electroelastic loadings is investigated under electrically impermeable boundary condition on the crack surface. Using Hankel transform technique, the mixed boundary value problem is reduced to a system of singular integral equations. With the aid of Gauss-Chebyshev integration technique, the integral equations are further reduced to a system of algebraic equations. The field intensity factor and energy release rate are determined. Numerical results reveal the effects of electric loadings and crack configuration on crack propagation and growth. The results seem useful for design of the piezoelectric structures and devices of high performance.     

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.     

Anti-plane analysis of multiple cracks originating from a circular hole in a magnetoelectroelastic solid

The anti-plane problem of multiple cracks originating from a circular hole in a magnetoelectroelastic solid is investigated under remotely uniform anti-plane mechanical loading and in-plane electromagnetic loadings. The boundary value problem is reduced to a Cauchy integral equation by a new mapping function and the complex variable method, which is further solved exactly. The analytic expressions of the complex potentials, the field intensity factors and the energy release rate are derived in closed-form under the assumption that the surfaces of the cracks and hole are both electrically and magnetically impermeable. The effects of crack configurations and combined loadings on the energy release rate are shown graphically. Several useful results which may have potential applications to the design and fracture analysis of magnetoelectroelastic structures are given.     

三维动态裂纹问题的超奇异积分方程法

基于弹性材料的动态基本方程，结合广义Betti-Rayleigh互易等式与时域下的边界积分方程，推导得到时域下的超奇异积分方程组。引入Laplace域下的动态基本解，将经过主部分析的积分核函数分解为静态和动态部分，其中动态积分核不具有奇异性。在裂纹前沿附近单元，采用与理论分析一致的平方根位移模型。结合Lubich时间卷积实现拉氏变换，采用配置点法计算超奇异积分，获得问题的数值解。并针对椭圆裂纹算例编写Fortran程序，得到冲击荷载作用下张开型裂纹的动态应力强度因子变化规律，数值结果稳定且收敛速度快。     

Stress intensity factors of the eccentric and edge cracked cylinders

Using the single crack solution and the regular solution of harmonic function,thetorsion problem of a cracked cylinder is reduced to solving a set of mixed-type integralequations which can be solved by combining the numerical method of singular integralequation with the boundary element method.Several numerical examples arecalculated and the stress intensity factors are obtained.     

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

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

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.     

Dynamic fracture analysis of an annular interfacial crack between dissimilar magnetoelectroelastic layers

Transient response of an annular interfacial crack between dissimilar magnetoelectroelastic layers under impacts is investigated. On the crack surface, magnetoelectrically impermeable boundary condition is adopted. Using Laplace and Hankel transform techniques, the mixed boundary value problem is reduced to a system of singular integral equations. The integral equations are further reduced to a system of algebraic equations with the aid of Jacobi polynomials. The dynamic field intensity factor and dynamic energy release rate are determined. Numerical results reveal the effects of electric or magnetic loadings and material parameters of composite on crack propagation and growth.     

Dynamic penny-shaped cracks in multilayer sandwich composites

A point force method is proposed for obtaining the dynamic elastic response of a multilayer sandwich composite in the presence of a penny-shaped crack under a harmonic loading. The sandwich composite is a multilayered solid whose lower half is the mirror image of the upper half with the center plane as the mirror. The crack is lying on the mirror plane of the composite. The solution of the mode I dynamic crack problem is formulated by integrating the Green’s function of a time-harmonic surface normal point force over the crack surface with an unknown point force distribution. The dual integral equations of the unknown point force distribution are established by considering the boundary conditions, which can be reduced to a Fredholm integral equation of the second kind. A complete solution of the crack problem under consideration can be obtained by solving this Fredholm integral equation. It will be shown that the results obtained by this approach are the same as some existing solutions.     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

Cylindrical interface crack in a hollow layered functionally graded cylinder under static torsion

A hollow functionally graded composite cylinder under static torsion, which consists of an inner and outer elastic circular tube with a cylindrical interface crack, is studied in this work. By utilizing Fourier integral transform method, the mixed boundary value problem is reduced to a Cauchy singular integral equation, from which the numerical results of the stress intensity factor are obtained by the Lobatto–Chebyshev quadrature technique. Numerical results demonstrate the coupled effects of geometrical, physical, and functionally graded parameters on the interfacial fracture behavior.     

Transient dynamic responses of a cracked solid subjected to in-plane loadings

Transient dynamic responses of an elastic cracked solid subjected to in-plane surface loadings are investigated in this study. Two vertical cracks, a surface-breaking crack and a sub-surface crack, are considered. The frequency responses of the plane strain problem are calculated by the computational mechanics combining the finite element method with the boundary integral equation. The finite element method is used for the near-field enclosing the crack, while the boundary integral equation is applied for the far-field to satisfy the Sommerfeld radiation condition. The transient responses are then obtained using fast Fourier transform. Surface displacements, crack opening displacements, and dynamic stress intensity factors are presented to show the significant effects of the cracks. The interaction between the elastic waves and the cracks as well as the mode conversion phenomena can be observed and understood through the numerical simulations.     

Dynamic analysis of interacting coplanar cracks in a half space with a clamped boundary condition using boundary integral equations

The three-dimensional dynamic problem of coplanar circular cracks in an elastic half-space with a clamped boundary condition is considered. The crack faces are subjected to harmonic loads. The problem is reduced to a system of two-dimensional boundary integral equations of the type of the Helmholtz potential for unknown discontinuities in the displacements of the opposite faces of the cracks. The stress intensity factors at the crack contours are obtained and discussed.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 1, pp. 153–159, January–February, 2005     

Analysis of anisotropic sector with a radial crack under anti-plane shear loading

In this paper, the anti-plane shear deformation of an anisotropic sector with a radial crack is investigated. The traction–traction boundary conditions are imposed on the radial edges and the traction-free condition is considered on the circular segment of the sector. A novel mathematical technique is employed for the solution of the problem. This technique consists of the use of some recently proposed finite complex transforms (Shahani, 1999), which have complex analogies to the standard finite Mellin transforms of the first and second kinds. However, it is essential to state the traction-free condition of the crack faces in the form of a singular integral equation which is done in this paper by describing an exact analytical method. The resultant dual integral equations are solved numerically to determine the stress intensity factors at the crack tips. In the special cases, the obtained results coincide with those cited in the literature.     

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

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

Circular welding problems with a crack

In this paper,the problem of an infinite plane with acircular hole welded by a nearly circular plate with a crackof different material is considered.The problem is tran-sformed into solving a certain boundary value problem ofanalytic functions and then reduced to solving a singularintegral equation along the crack.The formulas and somenumerical results of the factors of stress intensityfor the cases Mode Ⅰ and Mode Ⅱ are obtained at the endof this paper.     

STUDY ON DYNAMIC STRESS INTENSITY FACTORS OF DISK WITH A RADIAL EDGE CRACK SUBJECTED TO EXTERNAL IMPULSIVE PRESSURE

A dynamic weight function method is presented for dynamic stress intensity factors of circular disk with a radial edge crack under external impulsive pressure. The dynamic stresses in a circular disk are solved under abrupt step external pressure using the eigenfunction method. The solution consists of a quasi-static solution satisfying inhomogeneous boundary conditions and a dynamic solution satisfying homogeneous boundary conditions. By making use of Fourier-Bessel series expansion, the history and distribution of dynamic stresses in the circular disk are derived. Furthermore, the equation for stress intensity factors under uniform pressure is used as the reference case, the weight function equation for the circular disk containing an edge crack is worked out, and the dynamic stress intensity factor equation for the circular disk containing a radial edge crack can be given. The results indicate that the stress intensity factors under sudden step external pressure vary periodically with time, and the ratio of the maximum value of dynamic stress intensity factors to the corresponding static value is about 2.0.     

A crack of arbitrary shape inside an infinite transversely isotropic cylinder under arbitrary load

Summary  In the first part of the article an infinite circular cylinder is considered, made of transversely isotropic elastic material and weakened by a plane crack perpendicular to its axis O z. The crack is opened by an arbitrary normal stress. The second part is devoted to the same crack loaded by an arbitrary tangential stress. The complete solution in both cases is presented as a sum of the solution of a similar problem of a crack in an infinite space and an integral transform term, the parameters of which are determined from a set of linear algebraic equations derived from the boundary conditions. Governing integral equations with respect to the yet unknown crack displacement discontinuities are obtained. In the case of a circular crack, these equations can be inverted and solved by the method of consecutive interations. Received 30 November 2000; accepted for publication 3 May 2001