应用半权函数法计算界面裂纹的应力强度因子

应用半权函数法求解双材料界面裂纹的应力强度因子，得到以半权函数对参考位移与应力加权积分的形式表示的应力强度因子。针对特征值为复数λ的双材料界面裂纹裂尖应力和位移场，设置与之对应特征值为-λ的位移函数，即半权函数。半权函数的应力函数满足平衡方程，应力应变关系，界面的连续条件以及在裂纹面上面力为0；半权函数与裂纹体的几何尺寸无关，对边界条件没有要求。由功的互等定理得到应力强度因子KⅠ和KⅡ的积分形式表达式。本文计算了多种情况下界面裂纹应力强度因子的算例，与文献结果符合得很好。由于裂尖应力的振荡奇异性已经在积分中避免，只需考虑绕裂尖远场的任意路径上位移和应力，即使采用该路径上较粗糙的参考解也可以得到较精确的结果。     

考虑夹杂相互作用的复合陶瓷夹杂界面的断裂分析

复合材料中夹杂含量较高时,夹杂间的相互作用能显著改变材料细观应力应变场分布,基体和夹杂中的平均应力应变水平也会发生较大变化,导致复合材料强度等力学性能发生显著变化. 为修正单一夹杂模型运用在实际材料中的误差,基于相互作用直推估计法,建立一种考虑含夹杂相互作用的夹杂界面裂纹开裂模型. 首先根据相互作用直推估计法,得到残余应力和外载应力共同作用下夹杂中的平均应力,再计算无限大基体中相同的夹杂达到相同应力场时的等效加载应力,将此加载应力作为含界面裂纹夹杂的等效应力边界条件,在此边界条件下求得界面裂纹尖端的应力强度因子,进而得到界面裂纹开裂的极限加载条件,并分析了夹杂弹性性能、含量、热残余应力、夹杂尺寸等因素对界面裂纹开裂条件的影响. 结果表明,方法能够有效修正单夹杂模型运用在实际材料中的误差,较大的残余应力对界面裂纹开裂有重要的影响,夹杂刚度的影响并非单调且比较复杂;在残余应力较小时,降低柔性夹杂刚度或者增大刚性夹杂刚度都有利于提高材料强度;扩大夹杂尺寸将导致裂纹开裂极限应力显著降低,从而降低材料强度.      

双轴载荷作用下源于椭圆孔的分支裂纹的一种边界元分析

利用一种边界元方法来研究双轴载荷作用下无限大板中源于椭圆孔的分支裂纹．该边界元方法由Crouch与Starfied建立的常位移不连续单元和笔者提出的裂尖位移不连续单元构成．在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界,文中算例说明本数值方法对计算平面弹性裂纹的应力强度因子是非常有效的。该文对双轴载荷作用下无限大板中源于椭圆孔的分支裂纹的数值结果进一步证实本数值方法对计算复杂裂纹的应力强度因子的有效性,同时该数值结果可以揭示双轴载荷及裂纹体几何对应力强度因子的影响。     

圆柱形界面相裂纹在扭转载荷作用下的动态断裂

本文研究了位于界面相中的圆柱形界面裂纹的扭转冲击问题.采用Laplace、Fourier变换和位错密度函数将混合边值问题转化为求解Cauchy核奇异积分方程,利用Laplace数值反演技术计算了动态应力强度因子.讨论了材料特性和结构的几何尺寸对动态应力强度因子的影响.结果表明,随着界面相厚度的增加,无量纲化的动态应力强度因子减小.当裂纹靠近剪切弹性模量大的材料时,无量纲化的动态应力强度因子增大,反之减小.界面相两侧不同的材料组合对裂尖动态应力强度因子的影响是随着剪切弹性模量和质量密度的比值的增加而减小.界面相中裂纹长度对裂尖动态应力强度因子的影响比其他因素的影响大.     

双压电体界面上的电偶极子和裂纹5

求得了压电体双材料界面上的孤立二维电偶极子的解析解，结果表明某点电偶极子激发的应力－电位移场与该点到电偶极子的距离的平方成反比。研究了压电体双材料界面上的电偶极子对裂纹的作用，得到了问题的闭合解。在电偶极子的作用下，界面裂纹裂尖近区应力－电位移仍具有r^-1/2 iεα的振荡奇异性，文中求得裂尖应力强度因子，当电偶极子距裂尖距离ρ很近时，裂尖应力强度因子与ρ^-3/2-iεα成比例。     

双压电体界面上的电偶极子和裂纹

求得了压电体双材料界面上的孤立二维电偶极子的解析解,结果表明某点电偶极子激发的应力一电位移场与该点到电偶极子的距离的平方成反比.研究了压电体双材料界面上的电偶极子对裂纹的作用,得到了问题的闭合解.在电偶极子的作用下,界面裂纹裂尖近区应力-电位移仍具有r-1/2+iεα的振荡奇异性,文中求得裂尖应力强度因子,当电偶极子距裂尖距离ρ很近时,裂尖应力强度因子与ρ-3/2-iεα成比例.     

相变滞后对裂纹屏蔽效应的数值分析

本文针对奥氏体-马氏体双相材料,研究裂纹尖端区弥散分布的奥氏体颗粒在应变诱发时发生的相变对裂纹的屏蔽效应。鉴于实验中已发现的不同相变滞后对裂纹屏蔽效应的不同影响,本研究通过在裂纹尖端区不同位置嵌入相变颗粒,考虑到裂纹尖端区应力应变场的奇异分布及其诱发的相变,将裂纹尖端区相变滞后问题转化为相变颗粒在裂纹尖端区的位置问题。计及奥氏体-马氏体相变的体积膨胀效应进行了平面应力裂纹问题数值模拟,得到了单个相变夹杂对裂纹屏蔽效应的影响规律。结果表明：裂纹尖端区相变夹杂的位置对裂纹的屏蔽效应在距裂尖2倍夹杂直径以内影响极大,且以裂尖86度方向为界。其影响规律与McMeeking 和 Evans理论预言的60度方向不同。     

I型裂纹尖端约束应力区模型及其解析解

考虑了I型裂纹尖端损伤区域内三种不同的约束应力分布形式,即右三角分布形式(情况A)、均匀分布形式(情况B)、左三角分布形式(情况C),并采用复变函数方法求得了应力强度因子与裂纹张开位移的解析解;在此基础上,通过数值计算得到了应力强度因子和裂纹张开位移随约束应力区长度、约束应力大小以及分布形式的变化规律。研究结果表明:随裂尖材料损伤程度的增加,裂尖损伤区内约束应力减小,应力强度因子和裂纹张开位移增大;约束应力的分布形式对应力强度因子和裂纹张开位移有显著影响;相对于其他区域,约束应力对裂纹尖端区域裂纹张开位移的影响较大。然而,对于裂尖损伤区域的形成与作用荷载、材料性质、构件几何尺寸之间的关系,还需要进行更为深入的研究。     

关于界面下裂纹对界面裂纹的屏蔽参数

讨论了界面下裂纹对界面裂干涉的合理屏蔽参数问题。对于一条在远场载荷作用下、受到附近界面下裂纹干涉的界面裂纹,采用伪力法计算不同长度单位时界面裂纹的G／G0、k1/k10和k11/k110以及K1／K10和K11／K110（G是能量释放率、k1 ik110是通常定义的应力强度因子、K10＋iK110是含标定长度的应力强度因子;G0、k10 ik110、K10＋iK110对应单一界面裂纹情况）。结果表明,G／G0、K1／K10和K11／K110是较为合理的描述界面下裂纹对界面裂纹屏蔽的参数。     

干涉问题中T应力与各向异性的作用

采用离散模型（包括半无限主裂纹和近尖微裂纹）研究了各向异性材料主微裂纹干涉问题中T应力对主裂尖参数的影响,并且与相同情况下各向同性材料的结果进行了比较,比较结果列于文中各图。研究结果表明,在各向异性材料和各向同性材料中T应力对主裂尖应力强度因子的影响趋势是相似的,但是由于T应力与材料各向异性性质的共同作用,使两种情况下T应力对主裂尖参数的影响结果存在着明显的偏差。     

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 OF A DISLOCATION DIPOLE WITH A MODE Ⅲ INTERFACE CRACK

Interaction between a screw dislocation dipole and a mode Ⅲ interface crack is investigated. By using the complex variable method, the closed form solutions for complex potentials are obtained when a screw dislocation dipole lies inside a medium. The stress fields and the stress intensity factors at the tip of the interface crack produced by the screw dislocation dipole are given. The influence of the orientation, the dipole arm and the location of the screw dislocation dipole as well as the material mismatch on the stress intensity factors is discussed. zThe image force and the image torque acting on the screw dislocation dipole center are also calculated. The mechanical equilibrium position of the screw dislocation dipole is examined for various material property combinations and crack geometries. The results indicate that the shielding or anti-shielding effect on the stress intensity factor increases abruptly when the dislocation dipole approaches the tip of the crack. Additionally, the disturbation of the interface crack on the motion of the dislocation dipole is also significant.     

INTERACTION OF A DISLOCATION DIPOLE WITH A MODE III INTERFACE CRACK

Interaction between a screw dislocation dipole and a mode III interface crack is investigated. By using the complex variable method, the closed form solutions for complex potentials are obtained when a screw dislocation dipole lies inside a medium. The stress fields and the stress intensity factors at the tip of the interface crack produced by the screw dislocation dipole are given. The influence of the orientation, the dipole arm and the location of the screw dislocation dipole as well as the material mismatch on the stress intensity factors is discussed. The image force and the image torque acting on the screw dislocation dipole center are also calculated. The mechanical equilibrium position of the screw dislocation dipole is examined for various material property combinations and crack geometries. The results indicate that the shielding or anti-shielding effect on the stress intensity factor increases abruptly when the dislocation dipole approaches the tip of the crack. Additionally, the disturbation of the interface crack on the motion of the dislocation dipole is also significant.     

位错偶极子和Ⅲ型界面裂纹的干涉效应

研究了多晶体材料中螺型位错偶极子和界面裂纹的弹性干涉作用.利用复变函数方法,得到了该问题复势函数的封闭形式解答.求出了由位错偶极子诱导的应力场和裂纹尖端应力强度应子,分析了偶极子的方向,偶臂和位置以及材料失配对应力强度因子的影响.推导了作用在螺型位错偶极子中心的像力和力偶矩,并讨论了界面裂纹几何条件和不同材料特征组合对位错偶极子平衡位置的影响规律.结果表明,裂纹尖端的螺型位错偶极子对应力强度因子会产生强烈的屏蔽或反屏蔽效应.同时,界面裂纹对螺型位错偶极子在材料中运动有很强的扰动作用.     

Some problems in the antiplane shear deformation of bi-material wedges

The antiplane shear deformation of a bi-material wedge with finite radius is studied in this paper. Depending upon the boundary condition prescribed on the circular segment of the wedge, traction or displacement, two problems are analyzed. In each problem two different cases of boundary conditions on the radial edges of the composite wedge are considered. The radial boundary data are: traction–displacement and traction–traction. The solution of governing differential equations is accomplished by means of finite Mellin transforms. The closed form solutions are obtained for displacement and stress fields in the entire domain. The geometric singularities of stress fields are observed to be dependent on material property, in general. However, in the special case of equal apex angles in the traction–traction problem, this dependency ceases to exist and the geometric singularity shows dependency only upon the apex angle. A result which is in agreement with that cited in the literature for bi-material wedges with infinite radii. In part II of the paper, Antiplane shear deformation of bi-material circular media containing an interfacial edge crack is considered. As a special case of bi-material wedges studied in part I of the paper, explicit expressions are derived for the stress intensity factor at the tip of an edge crack lying at the interface of the bi-material media. It is seen that in general, the stress intensity factor is a function of material property. However, in special cases of traction–traction problem, i.e., similar materials and also equal apex angles, the stress intensity factor becomes independent of material property and the result coincides with the results in the literature.     

Analysis of stress intensity factor in orthotropic bi-material mixed interface crack

Adopting the complex function approach, the paper studies the stress intensity factor in orthotropic bi-material interface cracks under mixed loads. With consideration of the boundary conditions, a new stress function is introduced to transform the problem of bi-material interface crack into a boundary value problem of partial differential equations. Two sets of non-homogeneous linear equations with 16 unknowns are constructed. By solving the equations, the expressions for the real bi-material elastic constant εt and the real stress singularity exponents λt are obtained with the bi-material engineering parameters satisfying certain conditions. By the uniqueness theorem of limit,undetermined coefficients are determined, and thus the bi-material stress intensity factor in mixed cracks is obtained. The bi-material stress intensity factor characterizes features of mixed cracks. When orthotropic bi-materials are of the same material, the degenerate solution to the stress intensity factor in mixed bi-material interface cracks is in complete agreement with the present classic conclusion. The relationship between the bi-material stress intensity factor and the ratio of bi-material shear modulus and the relationship between the bi-material stress intensity factor and the ratio of bi-material Young's modulus are given in the numerical analysis.     

Dynamic interaction between a sub-interface crack and the interface in a bi-material: time-domain BEM analysis

Summary Transient response of a sub-interface crack in a bi-material is studied with emphasis on the dynamic interaction between the crack and the interface, by combining the traditional time-domain displacement boundary element method (BEM) and the non-hypersingular traction BEM. Computations are performed for an unbounded bi-material with a crack subjected to impact tensile loading on its faces or incident impact waves and a bounded rectangular bi-material plate under remote impact tensile loading. Numerical results of the dynamic stress intensity factors (DSIFs) and dynamic interface tractions are presented for various material combinations and crack locations. It is shown that pronounced increases in DSIFs and the interface tractions may be caused in some cases because of the dynamic interaction between the crack and the interface.This work was initialized during the second author's stay at Institute of Mechanics, TU Darmstadt, Germany under the support of the Alexander von Humboldt Foundation. Discussion on the BEM formulation with Dr. Seelig is gratefully acknowledged. The first two authors are also grateful for the partial support by the China National Natural Science Foundation under Grant No. 10025211 and the NJTU Scientific Paper Fund (PD195).     

Shielding effect and emission condition of a screw dislocation near a blunt crack in elliptical inhomogeneity

The shielding effect and emission condition of a screw dislocation near a blunt crack in elastic elliptical inhomogeneity is dealt with. Utilizing the Muskhelishvili complex variable method, the explicit series form solutions of the complex potentials in the matrix and the inclusion regions are derived. The stress intensity factor and critical stress intensity factor for dislocation emission are also calculated. The influences of the orientation of the dislocation and morphology of the blunt crack as well as the material elastic dissimilarity upon the shielding effect and emission criterion are discussed in detail. As a result, numerical analysis and discussion show that the positive screw dislocation can reduce the stress intensity factor of the crack tip (shielding effect) only when it is located in the certain region. The shielding effect increases with the increase of the shear modulus of the matrix and the curvature radius of the blunt crack tip, but decreases with the increase of dislocation azimuth angle. The critical loads at infinity for dislocation emission increases with the increment of the emission angle and the curvature radius of the blunt crack tip, and the most probable angle for screw dislocation emission is zero. The present solutions contain previous results as the special cases.     

Antiplane interaction of line crack with an arbitrarily located elliptical inclusion

The antiplane problem of the interaction between a main crack and an arbitrarily located elastic elliptical inclusion near its tip is addressed in the current study. The analysis is based on the use of the complex potentials for the antiplane problem, Laurent series expansion method and an appropriate superposition scheme. The stress intensity factor at the main crack is obtained in a general series form. Explicit asymptotic solutions are also derived by using a perturbation technique and retaining the leading order terms in series expansion. The present solutions are shown to coincide with the Taylor expansion of exact solutions for special cases available in the literature. Discussed are changes in the crack tip stress intensity which can be enhanced or suppressed depending on the location of the elliptical inclusion. The explicit solutions provided herein are well suited for the further quantitative analysis of toughening mechanisms in ceramic composite materials.     

A crack perpendicular to the bimaterial interface in finite solid

The dislocation simulation method is used in this paper to derive the basic equations for a crack perpendicular to the bimaterial interface in a finite solid. The complete solutions to the problem, including the T stress and the stress intensity factors are obtained. The stress field characteristics are investigated in detail. It is found that when the crack is within a weaker material, the stress intensity factor is smaller than that in a homogeneous material and it decreases when the distance between the crack tip and interface decreases. When the crack is within a stiffer material, the stress intensity factor is larger than that in a homogeneous material and it increases when the distance between the crack tip and interface decreases. In both cases, the stress intensity factor will increase when the ratio of the size of a sample to the crack length decreases. A comparison of stress intensity factors between a finite problem and an infinite problem has been given also. The stress distribution ahead of the crack tip, which is near the interface, is shown in details and the T stress effect is considered.