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

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

压电材料中椭圆片状裂纹问题精确解的一种方法

在线性压电陶瓷本构关系和裂纹边界绝缘的框架下,用超奇异积分方程的方法对椭圆类片状裂纹问题进行了重新研究.超奇异积分方程中的未知位移间断和电势间断近似地表示为基本密度函数与多项式之积,其中基本密度函数反映了椭圆片状裂纹前沿电弹性场的奇异性,而多项式在均布载荷作用下可用一个常数来表达.引入椭球坐标系后,得到了均布载荷作用下未知位移间断和电势间断的解析解.使用这些解析解和电弹性场强度的定义,得到了裂纹前沿Ⅰ型、Ⅱ型和Ⅲ型应力强度因子以及电位移强度因子的精确表达式.法向均布载荷作用下的结果与现有精确解完全一致,切向均布载荷作用下的结果则尚未见有其它报道.     

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

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

平行于异材界面的三维椭圆平片裂纹分析

讨论了拉伸载荷作用下平行于两相材料界面的椭圆平片裂纹问题．首先,使用有限部积分概念和两相材料界面完全接合时的点力基本解导出了一组以裂纹表面位移差为未知函数的超奇异积分方程组．该组方程表明,此时三种裂纹模型同时存在;其次,在数值求解该组方程的过程中,未知函数裂纹表面位移差被近似为位移差的基本密度函数与多项式之积．基本密度函数反映了裂纹前沿应力奇性性态;最后,以拉伸载荷为例,讨论了椭圆平片裂纹与界面的距离、裂纹形状比和不同材料组合对应力强度因子的影响,并以图表形式给出。     

圆形域多圆孔多裂纹反平面问题研究

本文运用复变函数及积分方程方法，求解了圆形域多圆孔多裂纹反平面问题，建立了两种类型的基本解。复叠加原理和所得的基本解并沿国圆孔和裂纹表面取待定的基本解密度函数，可得到一组以基本解密度函数为未知函数的Ｆｒｅｄｈｏｌｍ积分方程。通过该积分方程组的数值可以得到密度函数的离散值，进而得到了裂纹尖端的应力强度因子。     

反平面圆形夹杂和多圆孔多裂纹相互作用问题

动用复变函数及积分方法方法求解了反平面圆形夹杂和多圆孔多裂纹相互作用问题。为解决该问题,建立了两种类型的基本解。利用叠加原理和所得的基本解没圆孔和裂纹表面取待定的基本解密度函数,可得一组Ｆｒｅｄｈｏｌｍ积分方程,通过积分方程组的数值求解,可以得到密度函数的离散值,进而得到应力强度因子。     

界面上圆形衬砌结构对平面SH波散射

研究界面上的圆形衬砌结构对平面SH波散射与动应力集中问题.在一个含有半圆形衬砌缺口的弹性半空间水平面上,Green函数是受时间谐和的出平面线源载荷作用的位移基本解.采用沿界面“剖分”圆形衬砌结构的方法,并利用界面连续性条件建立起问题的定解积分方程组,进而得到圆形衬砌上的动应力集中解.最后给出了关于界面圆形衬砌结构上动应力集中系数的数值结果,并对界面圆形衬砌结构的动应力集中系数的影响进行了讨论.     

压电螺位错与含界面裂纹圆形涂层夹杂的干涉

研究位于基体或夹杂中任意点的压电螺型位错与含界面裂纹圆形涂层夹杂的电弹耦合干 涉问题. 运用复变函数方法,获得了基体,涂层和夹杂中复势函数的一般解答. 典型例 子给出了界面含有一条裂纹时,复势函数的精确级数形式解. 基于已获得的复势函数和广 义Peach-Koehler公式,计算了作用在位错上的像力. 讨论了裂纹几何条件,涂层厚度和材 料特性对位错平衡位置的影响规律. 结果表明,界面裂纹对涂层夹杂附近的位错运动有很大 的影响效应,含界面裂纹涂层夹杂对位错的捕获能力强于完整粘结情况;并发现界面裂纹长 度和涂层材料常数达到某一个临界值时可以改变像力的方向. 解答的特殊情形包含了以 往文献的几个结果.     

SH 波入射时半空间界面裂纹与圆形衬砌的相互作用

利用复变函数和Green函数法研究了双相介质半空间界面裂纹及界面附近圆形衬砌对SH 波的散射与动应力集中问题。首先,采用映像思想构造满足自由边界条件的散射波表达式,进而求解所需的Green函数;其次,采用裂纹切割技术构造裂纹,并根据连续性条件建立了求解该问题的无穷代数方程组;最后,给出了不同入射波数时界面裂纹与衬砌的相互作用。结果表明,裂纹的存在显著放大了衬砌界面的动应力集中。     

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

本文研究了圆形孔洞内衬砌与孔洞部分脱胶时对ＳＨ波的散射．将脱胶区看作表面不相接触的弧形界面裂纹，利用波函数展开法，并引入裂纹面的位错密度函数，将问题归结为一组奇异积分方程．通过数值计算获得了动应力强度因子（ＤＳＩＦ）和远场位移及散射截面（ＳＣＳ）．结果显示：由于脱胶，ＤＳＩＦ和ＳＣＳ在较低的频率上发生共振．     

Analysis of stress intensity factors of a ring-shaped interface crack

In this paper, numerical solutions of singular integral equations are discussed in the analysis of axi-symmetric interface cracks under torsion and tension. The problems of a ring-shaped interface crack are formulated in terms of a system of singular integral equations on the basis of the body force method. In the numerical analysis, unknown body force densities are approximated by the products of the fundamental density functions and power series, where the fundamental densities are chosen to express a two-dimensional interface crack exactly. The accuracy of the present analysis is verified by comparing the present results with the results obtained by other researchers for the limiting cases of the geometries. The calculation shows that the present method gives rapidly converging numerical results for those problems as well as for ordinary crack problems in homogeneous material. The stress intensity factors of a ring-shaped interface crack are shown in tables and charts with varying the material combinations and also geometrical conditions.     

Controlling parameter of the stress intensity factors for a planar interfacial crack in three-dimensional bimaterials

Numerical solutions of singular integral equations are discussed in the analysis of a planar rectangular interfacial crack in three-dimensional bimaterials subjected to tension. The problem is formulated as a system of singular integral equations on the basis of the body force method. In the numerical analysis, unknown body force densities are approximated by the products of the fundamental density functions and power series, where the fundamental density functions are chosen to express singular behavior along the crack front of the interface crack exactly. The calculation shows that the present method gives smooth variations of stress intensity factors along the crack front for various aspect ratios. The present method gives rapidly converging numerical results and highly satisfied boundary conditions throughout the crack boundary. The stress intensity factors are given with varying the material combination and aspect ratio of the crack. It is found that the stress intensity factors K_I and K_II are determined by the bimaterial constant ε alone, independent of elastic modulus ratio and Poisson’s ratio.     

Mode I fracture analysis of a piezoelectric strip based on real fundamental solutions

In existing papers, mode I crack problems of piezoelectric ceramics are generally solved in complex domain because of the complex fundamental solutions of in-plane piezoelectric governing equations. In fact, these problems can alternatively be analyzed in real number field by recasting the solutions in real form instead. The main purpose of the present work is to develop such real fundamental solutions by detailed eigenvalue and eigenvector analyses. As an example of application, the widely studied fracture problem of a piezoelectric strip with a center-situated crack under mode I loading condition is then revisited based on the real fundamental solutions. Mixed boundary value conditions of the crack are transformed into Cauchy singular integral equations, which are then solved numerically to get fracture parameters including the energy release rate and intensity factors. Convergence behaviors of the kernel functions are surveyed. Theoretical derivation and computation are validated by the exact solution in a special case. The effect of a combined geometrical parameter on the crack is also investigated.     

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.     

Stress intensity factors for an edge interface crack in a bonded semi-infinite plate for arbitrary material combination

Although a lot of interface crack problems were previously treated, few solutions are available under arbitrary crack lengths and material combinations. In this paper the stress intensity factors of an edge interface crack in a bonded strip are considered under tension with varying the crack length and material combinations systematically. Then, the limiting solutions are provided for an edge interface crack in a bonded semi-infinite plate under arbitrary material combinations. In order to calculate the stress intensity factors accurately, exact solutions in an infinite bonded plate are also considered to produce proportional singular stress fields in the analysis of FEM by superposing specific tensile and shear stresses at infinity. The details of this new numerical solution are described with clarifying the effect of the element size on the stress intensity factor. It is found that for the edge interface crack the normalized stress intensity factors are not always finite depending upon Dunders’ parameters. This behavior can be explained from the condition of the singular stress at the end of bonded strip. Convenient formulas are also given by fitting the computed results.     

Interaction between an interface crack and subinterface microcracks in metal/piezoelectric bimaterials

The J-integral analysis is presented for the interaction problem between a semi-infinite interface crack and subinterface matrix microcracks in dissimilar anisotropic materials. After deriving the fundamental solutions for an interface crack subjected to different loads and the fundamental solutions for an edge dislocation beneath the interface, the interaction problem is deduced to a system of singular integral equations with the aid of a superimposing technique. The integral equations are then solved numerically and a conservation law among three values of the J-integral is presented, which are induced from the interface crack tip, the microcracks and the remote field, respectively. The conservation law not only provides a necessary condition to confirm the numerical results derived, but also reveals that the microcrack shielding effect in such materials could be considered as a redistribution of the remote J-integral. It is this redistribution that does lead to the phenomenological shielding effect.     

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

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