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

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

横观各向同性压电材料中裂纹问题的边界元法

采用Somigiliana公式给出了三维横观各向同性压电材料中的非渗漏裂纹问题的一般解和超奇异积分方程,其中未知函数为裂纹面上的位移间断和电势间断．在此基础上,使用有限部积分和边界元结合的方法,建立了超奇异积分方程的数值求解方法,并给出了一些典型数值算例的应力强度因子和电位移强度因子的数值结果,结果令人满意．     

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

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

三维体内含垂直于双材料界面混合型裂纹的超奇异积分解法

利用双材料位移基本解和Somigliana公式,将三维体内含垂直于双材料界面混合型裂纹问题归结为求解一组超奇异积分方程。使用主部分析法,通过对裂纹前沿应力奇性的分析,得到用裂纹面位移间断表示的应力强度因子的计算公式,进而利用超奇异积分方程未知解的理论分析结果和有限部积分理论,给出了超奇异积分方程的数值求解方法。最后,对典型算例的应力强度因子做了计算,并讨论了应力强度因子数值结果的收敛性及其随各参数变化的规律。     

压电材料中三维I型断裂力学分析

讨论了不可导通情况下三维横观各向同性压电材料中受拉伸和电载荷作用的平片裂纹I型断裂力学问题. 使用有限部积分概念，从三维线性压电理论出发，严格得到了一组以裂纹面位移间断和电势间断为未知变量的超奇异积分方程组；应用二维超奇异积分的主部分析法，从理论上分析得到了裂纹前沿应力和电势奇性指数以及应力和电位移奇性场，从而找到了以裂纹面位移间断和电势间断表示的应力和电位移强度因子、能量释放率表达式；为所得到的超奇异积分方程组建立了数值法，并用此计算了若干典型的平片裂纹问题，数值结果令人满意.     

双材料中平片裂纹问题的超奇异积分方程解法

利用三维断裂力学的超奇异积分方程方法,对双材料空间中重直于界面的平片裂纹Ⅰ型问题进行了研究。首先根据双材料空间的弹性力学基本解,使用边界积分方程方法,在有限部积分的意义下导出了以裂纹面位罗间断为未知函数的超奇异积分方程,并为其建立了数值法。在此基础上,讨论了用裂纹面位移问题计算应力强度因子的方法。最后用此计算了几个典型的Ⅰ型下片裂纹问题的应力强度因子,其数值结果令人满意。     

磁电热弹耦合材料三维多裂纹超奇异积分法

用超奇异积分方程法将多场耦合载荷作用下磁电热弹耦合材料内含任意形状和位置三维多裂纹问题转化为求解一以广义位移间断为未知函数的超奇异积分方程组问题,退化得到内含任意形状平行三维多裂纹问题的超奇异积分方程组;推导出平行三维多裂纹问题的裂纹前沿广义奇异应力场解析表达式、定义了广义(应力、应变能)强度因子和广义能量释放率;应用有限部积分概念及体积力法,为超奇异积分方程组建立了数值求解方法,编制了FORTRAN程序,以平行双裂纹为例,通过典型算例,研究了广义(应力、应变能)强度因子随裂纹位置、裂纹形状及材料参数变化规律,得到裂纹断裂评定准则. 最后,分析了裂纹间干扰、屏蔽作用及其在工程实际中的应用.      

三维界面裂纹的奇性应力场和应力强度因子分析

使用有限部积分概念和极限方法得到了三维平片界面裂纹的超奇异积分－微分方程组后,进一步利用二维超奇异积分主部分析方法,对裂纹前沿的应力场作了理论分析,并获得了其奇性应力场和裂纹面位移间断表示复位应力强度因子的精确表达式,为三维平片界面裂纹的超奇异积分－微分方程组的求解建立了数值方法,并分析了界面椭圆平片裂纹问题,和现有解比较,所得数值结果令人满意。     

三维无限体中两平行平片裂纹相互干扰的超奇异积分方程解法

本文利用三维断裂力学的超奇异积分方程求解理论，对三维无限体中两平行平片裂纹在任意载荷作用下的相互干扰问题作了研究。首先导出了以裂纹面移间断（位借）为未知函数的超奇异积分方程组，然后为其建立了有限积分边界元法；在此基础上，讨论用了裂纹面位移间断计算应力强度因子的方法，最后用此计算了两平行平片裂纹的相对位置对裂前沿应力强度因子的影响，其数值结果令人满意。     

同震断层三维裂纹扩展波动时域超奇异积分法

应用波动时域超奇异积分法将P波、S波和磁电热弹多场耦合作用下同震断层任意形状三维裂纹扩展问题转化为求解以广义位移间断率为未知函数的超奇异积分方程组问题;定义了广义应力强度因子,得到裂纹前沿广义奇异应力增量解析表达式;应用波动时域有限部积分概念及体积力法,为超奇异积分方程组建立了数值求解方法,编制了FORTRAN程序,以三维矩形裂纹扩展问题为例,通过典型算例,研究了广义应力强度因子随裂纹位置变化规律;分析了同震断层裂纹扩展中力、磁、电场辐射规律.      

TORSIONAL VIBRATIONS OF RIGID CIRCULAR PLATE ON TRANSVERSELY ISOTROPIC SATURATED SOIL

An analytical method was presented for the torsional vibrations of a rigid disk resting on transversely isotropic saturated soil. By Hankel transform, the dynamic governing differential equations for transversely isotropic saturated poroelastic medium were solved. Considering the mixed boundary-value conditions, the dual integral equations of torsional vibrations of a rigid circular plate resting on transversely isotropic saturated soil were established. By appropriate transform, the dual integral equations were converted into a Fredholm integral equation of the second kind. Subsequently, the dynamic compliance coefficient, the torsional angular amplitude of the foundation and the contact shear stress were expressed explicitly. Selected examples were presented to analyse the influence of saturated soil's anisotropy on the foundation's vibrations.     

横观各向同性饱和地基上中厚圆板的非轴对称振动

研究横观各向同性饱和土地基上中厚弹性圆板的非轴对称振动问题。基于横观各向同性饱和介质Biot波动方程的一般解,按混合边值问题建立了饱和地基与弹性中厚圆板非轴对称动力相互作用的对偶积分方程,并将对偶积分方程转化为易于计算的第二类Fredholm积分方程;采用数值方法求解该积分方程。数值算例结果表明,当h/a>0.05时,饱和半空间体上中厚度圆板在不同频率下的振动特性与相应频率下的刚性板的振动特性基本相同,当h/a<0.05时,板中心的位移将随h/a的减小而增大。     

Stationary heat conduction in a macro- and microperiodically layered solid

Summary This paper deals with the stationary heat conduction in a solid consisting of planar, isotropic or transversely isotropic layers. The following cases are considered: (1) arbitrarily layered solid, (2) macroperiodically layered solid, consisting of equal pairs of different basic layers with finite thicknesses, (3) microperiodically layered solid, consisting of equal layer groups of two or more different basic layers of infinitesimal thicknesses. Assuming a perfect contact between the layers, exact solutions for the plane, the axisymmetric and the general three-dimensional problem of macro- and microperiodically layered semispaces, respectively, are derived using integral transforms and the transfer-matrix method. It is proved that the microperiodically layered solid with isotropic or transversely isotropic basic layers is equivalent to a homogeneous transversely isotropic solid. Received 31 March 1999; accepted for publication 21 May 1999     

Dynamic interaction of a pile with a transversely isotropic elastic half-space under transverse excitations

This investigation is concerned with a mathematical analysis of an elastic circular cylindrical pile embedded in a transversely isotropic half-space under lateral dynamic excitations. A combination of time-harmonic horizontal shear force and moment are applied at the top end of the pile. The boundary value problem is formulated by decomposing the pile-medium system into a fictitious pile and an extended transversely isotropic half-space. A Fredholm integral equation of the second kind governs the interaction problem, whose solution is then computed numerically. Selected results for dynamic compliance bending moment, displacement and slope profiles are presented for different transversely isotropic half-spaces to portray the influence of degree of anisotropy of the medium on various aspects of the solution.     

Non-axisymmetrical vibration of elastic circular plate on layered transversely isotropic saturated ground

The non-axisymmetrical vibration of elastic circular plate resting on a layered transversely isotropic saturated ground was studied.First,the 3-d dynamic equations in cylindrical coordinate for transversely isotropic saturated soils were transformed into a group of governing differential equations with 1-order by the technique of Fourier ex- panding with respect to azimuth,and the state equation is established by Hankel integral transform method,furthermore the transfer matrixes within layered media are derived based on the solutions of the state equation.Secondly,by the transfer matrixes,the general solutions of dynamic response for layered transversely isotropic saturated ground excited by an arbitrary harmonic force were established under the boundary conditions, drainage conditions on the surface of.ground as well as the contact conditions.Thirdly, the problem was led to a pair of dual integral equations describing the mixed boundary- value problem which can be reduced to the Fredholm integral equations of the second kind solved by numerical procedure easily.At the end of this paper,a numerical result concerning vertical and radical displacements both the surface of saturated ground and plate is evaluated.     

Analysis of a crack in a half-plane piezoelectric solid with traction-induction free boundary

In this paper, the problem of a crack embedded in a half-plane piezoelectric solid with traction-induction free boundary is analyzed. A system of singular integral equations is formulated for the materials with general anisotropic piezoelectric properties and for the crack with arbitrary orientation. The kernel functions developed are in complex form for general anisotropic piezoelectric materials and are then specialized to the case of transversely isotropic piezoelectric materials which are in real form. The obtained coupled mechanical and electric real kernel functions may be reduced to those kernel functions for purely elastic problems when the electric effects disappear. The system of singular integral equations is solved numerically and the coupling effects of the mechanical and electric phenomena are presented by the generalized stress intensity factors for transversely isotropic piezoelectric materials.     

Application of generalized images method to contact problems for a transversely isotropic elastic layer on a smooth half-space

The idea, first used by the author for the case of crack problems, is applied here to solve a contact problem for a transversely isotropic elastic layer resting on a smooth elastic half-space, made of a different transversely isotropic material. A rigid punch of arbitrary shape is pressed against the layer’s free surface. The governing integral equation is derived; it is mathematically equivalent to that of an electrostatic problem of an infinite row of coaxial charged disks in the shape of the domain of contact. The case of circular domain of contact is considered in detail. As a comparison, the method of integral transforms is also used to solve the problem. The main difference of our integral transform approach with the existing ones is in separating of our half-space solution from the integral transform terms. It is shown that both methods lead to the same results, thus giving a new interpretation to the integral transform as a sum of an infinite series of generalized images.     

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.     

Contact problems for several transversely isotropic elastic layers on a smooth elastic half-space

The idea of generalized images, first used by the author for the case of crack problems, is applied here to solve a contact problem for n transversely isotropic elastic layers, with smooth interfaces, resting on a smooth elastic half-space, made of a different transversely isotropic material. A rigid punch of arbitrary shape is pressed against the top layer’s free surface. The governing integral equation is derived for the case of two layers; it is mathematically equivalent to that of an electrostatic problem of an infinite row of coaxial charged disks in the shape of the domain of contact. This result is then generalized for an arbitrary number of layers. As a comparison, the method of integral transforms is also used to solve the problem. The main difference of our integral transform approach with the existing ones is in separating of our half-space solution from the integral transform terms. It is shown that both methods lead to the same results, thus giving a new interpretation to the integral transform as a sum of an infinite series of generalized images.     

Exact Electromagnetothermoelastic Solution for a Transversely Isotropic Piezoelectric Hollow Sphere Subjected to Arbitrary Thermal Shock

This paper presents analytical study for electromagnetothermoelastic transient behavior of a transversely isotropic hollow sphere, placed in a uniform magnetic field, subjected to arbitrary thermal shock. Exact solutions for the transient responses of stresses, perturbation of magnetic field vector, electric displacement and electric potential in the transversely isotropic piezoelectric hollow sphere are obtained by means of the Hankel transform, the Laplace transform and their inverse transforms. An interpolation method is used to solve the Volterra integral equation of the second kind caused by interactions among electric, magnetic, thermal and elastic fields. From the sample numerical calculations, it is seen that the present method is suitable for the transversely isotropic hollow sphere, placed in a uniform magnetic field, subjected to arbitrary thermal shock. Finally, the result can be used as a reference to solve other transient coupling problems of electromagnetothermoelasticity.