横观各向同性三维热弹性力学通解及其势理论法

通过引入两个位移函数，对用位移表达的运动平衡方程作了简化．利用算子理论，严格地导出了横观各向同性非耦合热弹性动力学问题的通解．对于静力学问题，通解的形式可进一步简化成用4个准调和函数来表示．具体考察了横观各向同性体内平面裂纹上下表面有对称分布温度作用的问题，推广了势理论方法，导出了一个积分方程和一个微分-积分方程．针对币状裂纹表面受均布温度作用情形，给出了具体的解。     

A penny-shaped interface crack between a functionally graded piezoelectric layer and a homogeneous piezoelectric layer

The problem of a penny-shaped interface crack between a functionally graded piezoelectric layer and a homogeneous piezoelectric layer is investigated. The surfaces of the composite structure are subjected to both mechanical and electrical loads. The crack surfaces are assumed to be electrically impermeable. Integral transform method is employed to reduce the problem to a Fredholm integral equation of the second kind. The stress intensity factor, electric displacement intensity factor and energy release rate are derived, some typical numerical results are plotted graphically. The effects of electrical loads, material nonhomogeneity and crack configuration on the fracture behaviors of the cracked composite structure are analyzed in detail.     

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

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

Three-dimensional continuum analysis for a unidirectional composite with a broken fiber

The three-dimensional problem of a periodic unidirectional composite with a penny-shaped crack traversing one of the fibers is analyzed by the continuum equations of elasticity. The solution of the crack problem is represented by a superposition of weighted unit normal displacement jump solutions, everyone of which forms a Green’s function. The Green’s functions for the unbounded periodic composite are obtained by the combined use of the representative cell method and the higher-order theory. The representative cell method, based on the triple discrete Fourier transform, allows the reduction of the problem of an infinite domain to a problem of a finite one in the transform space. This problem is solved by the higher-order theory according to which the transformed displacement vector is expressed by a second order expansion in terms of local coordinates, in conjunction with the equilibrium equations and the relevant boundary conditions. The actual elastic field is obtained by a numerical evaluation of the inverse transform. The accuracy of the suggested approach is verified by a comparison with the exact analytical solution for a penny-shaped crack embedded in a homogeneous medium. Results for a unidirectional composite with a broken fiber are given for various fiber volume fractions and fiber-to-matrix stiffness ratios. It is shown that for certain parameter combinations the use of the average stress in the fiber, as it is employed in the framework of the shear lag approach, for the prediction of composite’s strength, leads to an over estimation. To this end, the concept of “point stress concentration factor” is introduced to characterize the strength of the composite with a broken fiber. Several generalizations of the proposed approach are offered.     

Penny-shaped interfacial crack between dissimilar magnetoelectroelastic layers

A penny-shaped interfacial crack between dissimilar magnetoelectroelastic layers subjected to magnetoelectromechanical loads is investigated,where the magnetoelectrically impermeable crack surface condition is adopted. By using Hankel transform technique,the mixed boundary value problem is firstly reduced to a system of singular integral equations,which are further reduced to a system of algebraic equations. The field intensity factors and energy release rate are finally derived. Numerical results elucidate the eects of crack configuration,electric and/or magnetic loads,and material parameters of the magnetoelectroelastic layers on crack propagation and growth. This work should be useful for the design of magnetoelectroelastic composite structures.     

A penny-shaped crack in a transversely isotropic piezoelectric layer

In this paper, we develop a model to treat penny-shaped crack configuration in a piezoelectric layer of finite thickness. The piezoelectric layer is subjected to axially symmetric mechanical and electrical loads. Hankel transform technique is used to reduce the problem to the solution of a system of integral equations. A numerical solution for the crack tip fields is obtained for different crack radius and crack position.     

Fracture behavior of a bonded magneto-electro-elastic rectangular plate with an interface crack

In this paper the anti-plane problem for an interface crack between two dissimilar magneto-electro-elastic plates subjected to anti-plane mechanical and in-plane magneto-electrical loads is investigated. The interface crack is assumed to be either magneto-electrically impermeable or permeable, and the position of the interface crack is arbitrary. The finite Fourier transform method is employed to reduce the mixed boundary-value problem to triple trigonometric series equations. The dislocation density functions and proper replacement of the variables are introduced to reduce these series equations to a standard Cauchy singular integral equation of the first kind. The resulting integral equation together with the corresponding single-valued condition is approximated as a system of linear algebra equations which can be easily solved. Field intensity factors and energy release rates are determined numerically and discussed in detail. Numerical results show the effects of crack configuration and loading combination parameters on the fracture behaviors of crack tips according to energy release rate criterion. The study of this problem is expected to have applications to the investigation of dynamic fracture properties of magneto-electro-elastic materials with cracks.     

横观各向同性材料的三维断裂力学问题

从三维横观各向同性材料弹性力学理论出发, 使用Hadamard有限部积分概念, 导出了三维状态下单位位移间断(位错)集度的基 本解. 在此基础上, 进一步运用极限理论, 将任意载荷作用下, 三维无限大横观各向 同性材料弹性体中, 含有一个位于弹性对称面内的任意形状的片状裂纹问题, 归结为求 解一组超奇异积分方程的问题. 通过二维超奇异积分的主部分析方法, 精确地求得了裂纹前沿光滑点附近的应力奇异指数和奇异应力场, 从而找到了以裂纹表面位移间断表示的应力强度因子表达式及裂纹局部扩展所提供 的能量释放率. 作为以上理论的实际应用，最后给出了一个圆形片状裂纹问题 的精确解例和一个正方形片状裂纹问题的数值解例. 对受轴对称法向均布载荷作用下圆形片状裂纹问题, 讨论了超奇异积分方程的精确求解方法, 并获得了位移间断和应力强度因子的封闭解, 此结果与现有理论解完全一致.     

Mixed mode axisymmetric annular cracks in a finite layer: Off angle crack initiation

The solutions of axisymmetric Volterra type climb and glide edge dislocations are obtained in a layer by means of the Hankel transforms. Utilizing the same procedure, Green’s function solution is obtained for a layer under self-equilibration normal ring traction. The distributed dislocation technique is used to construct integral equations for a system of co-axial annular cracks where the layer is under axisymmetric normal loads. These equations are solved numerically to obtain dislocation density on the cracks surfaces. The results are employed to determine stress intensity factors for annular and penny-shaped cracks and the interaction between two co-axial penny-shaped cracks is studied. Moreover, the stress intensity factors of the interacting cracks are determined such that they can be further used in conjunction with strain energy density (SED) failure criterion to obtain the possible direction of crack initiation that may not be apparent under mixed mode conditions.     

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.     

On the penny-shaped crack in inhomogeneous elastic materials under normal extension

The solution of the problem of a penny-shaped crack in an inhomogeneous material with elastic coefficients which are varying continuously along the direction perpendicular to the crack is examined in this paper. We studied the problem for an inhomogeneous material which satisfies the conditions of either torsional deformation and normal extension. A series form solution to the problem is proposed and analytical expressions for the first two terms of the series are obtained by using a Hankel transform technique. In the solution a homogeneous body is chosen as the reference so that inhomogeneous quantities are treated as being perturbed from the zeros reference solutions. Closed form expressions for the relevant stress intensity factors and the crack energy are derived and specific cases of the problem are also considered.     

On the half-infinite crack problem in thermo-electro-elasticity

This paper concerns itself with the fundamental solutions of the thermo-electro-elastic field in an infinite medium, weakened by a half-infinite plane crack with two identical point thermal loads applied on the crack surfaces. The corresponding mixed boundary value problem is solved by virtue of the potential theory method conjugated with the general solutions. The boundary governing equations are solved by using the results available in literature. Exact and complete three-dimensional (3D) fundamental solutions are presented in terms of elementary functions. The singularity at the crack tip is analyzed explicitly. The obtained solutions will be of high significance to the related BEM analysis.     

Dynamic response of an elastic medium containing crack

A fundamental solution for an infinite elastic medium containing a penny-shaped crack subjected to dynamic torsional surface tractions is attempted. A double Laplace–Hankel integral transform with respect to time and space is applied both to motion equation and boundary conditions yielding dual integral equations. The solution of the derived dual integral equations is based on an analytic procedure using theorems of Bessel functions and ordinary differential equations. The dynamic displacements’ field is obtained by inversion of the corresponding Laplace–Hankel transformed variable. Results of a representative example for a crack subjected to pulse surface tractions are obtained and discussed.     

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.     

An antisymmetric problem of a penny-shaped crack in a piezoelectric medium

Summary  An exact, three-dimensional analysis is developed for a penny-shaped crack in an infinite transversely isotropic piezoelectric medium. The crack is assumed to be parallel to the plane of isotropy, with its faces subjected to a couple of concentrated normal forces and a couple of point electric charges that are antisymmetric with respect to the crack plane. The fundamental solution of a concentrated force and a point charge acting on the surface of a piezoelectric half-space is employed to derive the integral equations for the general boundary value problem. For the above antisymmetric crack problem, complete expressions for the elastoelectric field are obtained. A numerical calculation is finally performed to show the piezoelectric effect in piezoelectric materials. It is noted here that the present analysis is an extension of Fabrikant's theory for elasticity. Received 30 August 1999; accepted for publication 1 March 2000     

复合型运动荷载作用下横观各向同性体三维裂纹的应力强度因子

讨论横观各向同性体中含一半平面裂纹,在裂纹面上作用有运动点荷载的三维复合型应力强度因子历史,通过积分变换技术,最终将问题归结为求解Wiener-Hopf型积分方程组,该文给出了求解这一类积分方程组的一般性方法,在此基础上,基于Abel定理和Cagniard-deHoop方法,求得Ⅱ、Ⅲ型复合应力强度因子的解析解,最后通过数值结果揭示了横观各向同性材料三维方尖端场的动态特性。     

CRACK TIP PLASTICITY OF A THERMALLY LOADED PENNY-SHAPED CRACK IN AN INFINITE SPACE OF 1D QC

The present work is concerned with a penny-shaped Dugdale crack embedded in an infinite space of one-dimensional(1D) hexagonal quasicrystals and subjected to two identical axisymmetric temperature loadings on the upper and lower crack surfaces. Applying Dugdale hypothesis to thermo-elastic results, the extent of the plastic zone at the crack tip is determined.The normal stress outside the plastic zone and crack surface displacement are derived in terms of special functions. For a uniform loading case, the corresponding results are presented by simplifying the preceding results. Numerical calculations are carried out to show the influence of some parameters.     

Thermal stress analysis of a three-dimensional anticrack in a transversely isotropic solid

A three-dimensional analysis is performed for an infinite transversely isotropic elastic body containing an insulated rigid sheet-like inclusion (an anticrack) in the isotropy plane under a remote perpendicularly uniform heat flow. A general solution scheme is presented for the resulting boundary-value problems. Accurate results are obtained by constructing suitable potential solutions and reducing the thermal problem to a mechanical analog for the corresponding isotropic problem. The governing boundary integral equation for a planar anticrack of arbitrary shape is obtained in terms of a normal stress discontinuity. As an illustration, a complete solution for a rigid circular inclusion is obtained in terms of elementary functions and analyzed. This solution is compared with that corresponding to a penny-shaped crack problem.     

含币形裂纹的压电体在点力和点电荷作用下的裂纹张开位移

本文利用Chen和Shioya给出的在横观各向同性压电无限体内币形裂纹上下表面作用对称法向点力和点电荷情形下的解,结合压电材料之功的互等定,用初等函数的形式给出了在压电无限体中任意一点作用任意点力和点电荷情形下币形裂纹的张开位移,并对PZT－4压电陶瓷和非压电材料作了计算分析。     

Propagation of a penny-shaped fluid-driven fracture in an impermeable rock: asymptotic solutions

This paper presents an analysis of the propagation of a penny-shaped hydraulic fracture in an impermeable elastic rock. The fracture is driven by an incompressible Newtonian fluid injected from a source at the center of the fracture. The fluid flow is modeled according to lubrication theory, while the elastic response is governed by a singular integral equation relating the crack opening and the fluid pressure. It is shown that the scaled equations contain only one parameter, a dimensionless toughness, which controls the regimes of fracture propagation. Asymptotic solutions for zero and large dimensionless toughness are constructed. Finally, the regimes of fracture propagation are analyzed by matching the asymptotic solutions with results of a numerical algorithm applicable to arbitrary toughness.