梯度非均匀有限层-基结构界面裂纹动态断裂的数值分析

研究了具有一个界面裂纹的有限尺寸梯度非均匀层-基结构在面内冲击载荷作用下的动态响应问题;提出了分析这类问题的一种数值方法,其核心是计算具有一定宽度的裂缝位移场,通过对到缝端距离及缝宽的双重插值求裂尖处动态应力强度因子。文中对涂层是弹性模量及密度连续变化的梯度非均匀材料,基体是各向同性均匀材料的层-基结构作了具体计算,数值结果清楚地显示出动态应力强度因子的时程变化规律,及有关参数(梯度参数,层厚,层长)对它的影响。     

瓷修复体界面断裂行为的模拟实验研究

本文利用云纹干涉法和云纹干涉－－有限元混合法，对瓷修复体的模拟双材料模型界面断裂问题进行了实验研究。用云纹干涉和数字错位云纹干涉法测量带边裂纹的双材料四点简支梁在剪切作用下界面表面的剪应变分布及界面两侧局部表面的位移场，实验表明，由于界面两两侧材料力学性质不同，表现出界面剪切断裂问题的非称性和裂尖附近复合型断裂的特点；用云纹干涉法和有限元法相结合的混合法对粘接界面角点应力奇异性进行研究，并对角点附近应力应变场作了分析，得到了应力奇异指数与边界楔角，载荷的关系，证明了用界面应力强度因子Kf来描述界面端部区域应力分布的公式，并得到了双材料界面端部区域的应力应变分布情况。本文的实验结果为进一步研究口腔金瓷修复体界面的优化设计提供了基础，同时也说明云纹干涉法对于双材料界面断裂行为的研究是有效的。     

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

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

有限尺寸下双材料界面附近裂纹位置对裂纹尖端的影响

随着复合材料的应用和发展,不同材料组成的界面结构越来越受到人们的重视。界面层两侧材料的性能相异会引起材料界面端奇异性,同时界面和界面附近存在裂纹会引起裂尖处的应力奇异性。因此双材料界面附近的力学分析是比较复杂的。本文建立双材料直角界面模型,在材料界面附近预设初始裂纹,计算了有限材料尺寸对界面应力场及其附近裂纹应力强度因子的影响。运用弹性力学中的 Goursat 公式求得直角界面端在有限尺寸下的应力场以及其应力强度系数。通过叠加原理和格林函数法进一步得到在直角界面端附近的裂纹尖端应力强度因子。计算结果表明,在适当范围内改变材料内裂纹与界面之间的距离,界面附近裂纹尖端的应力强度因子随着裂纹与界面距离的增加而减少,并且逐渐趋于稳定。分析结果可以为预测双材料结构复合材料界面失效位置提供参考。     

混合边界多裂隙体在简谐变温场作用下的断裂问题（张开型）

给出了无限平面作用有简谐变化的点热源时的位移场、应力场基本解、用间接法构造出混合边值多裂隙体在简谐变温场作用下的热断裂问题的边界积分方程,并离散求解.数值结果表明,该方法求解多裂隙体的简谐热断裂问题精度好,计算工作量少.文中计算了含边界裂缝的平板、含三条平行裂缝的平板在简谐变温场作用下缝端应力强度因子的变化过程,并与实验结果进行了比较,两者吻合良好.     

界面端附近裂纹的应力强度因子

结合材料的断裂形式可分为从界面端产生裂纹（沿界面或向母材内部层折）然后断裂与稍稍离开界面端处产生裂纹然后断裂这两种情况，在金属／陶瓷类结合材料中，后者出现的概率更大，本文利用结合材料界面端的奇异应力场和叠加原理，给出了界面端附近裂纹的应力强度因子近似计算公式，并用边界元数值计算验证了其有效性。     

双材料自由边界面端部问题的数值分析

对各向异性双材料自由边界面端部奇异性场问题进行了研究,利用有限元分析法所得到的各向异性双材料自由边界面端部的应力奇异性指数以及角分布函数,构造了一个自由边界面端部单元,据此建立了自由边界面端部奇异性场的杂交应力模型,并结合Hellinger-Reissner变分原理导出应力杂交元方程,建立了求解平面各向异性材料裂纹尖端问题的杂交元计算模型.与四节点单元相结合,提出一种求解自由边界面端部广义应力强度因子的杂交元法.考核例结果表明:本文方法的数值解精度高,可应用于各向异性材料双材料自由边界面端部问题.     

一种确定结合材料界面端应力强度因子的数值方法

基于弹性力平面问题的基本方程，给出了结合材料界面端的应力奇异性特征方程以及位移场和奇异应力场。提出了一种确定结合材料界面端应力强度因子的数值外插方法。对界面端区域进行了有限元网格单元划分。经过具体实例检验进一步确定了求解应力强度因子的最佳方向，该数值外插法的计算结果精度符合工程应用的要求，为工程材料强度的评价提供了有效的计算途径。     

Ⅰ-Ⅱ复合型裂缝断裂角的分层剪滞模型与分析

采用修正的剪滞理论建立了岩石、混凝土等准脆性材料的Ⅰ-Ⅱ复合型裂缝在单向拉伸荷载作用下的计算模型.得到了与实验相吻合且优于传统S判据的断裂角.通过对远场应力、斜裂缝区应力以及子层位秽的合理简化,得到了求解剪滞分析模型的边界条件,进而得到了含斜裂缝的各子层位移分布函数.引入最大应力集中因子,对Ⅰ-Ⅱ复合型裂缝前缘应力场进行简化;基于斜裂缝沿最大应力集中因子方向扩展,得到裂缝的断裂角.根据斜裂缝的应力分布,设置不同的子层分区,得到了更为细化的位移分布模式.通过对计算数据的分析,针对单向拉伸荷载作用下的Ⅰ-Ⅱ复合型裂缝,建立了按应力场分区设置子层的分层剪滞模型,得到更为精确的斜裂缝断裂角.     

双材料界面裂纹复应力强度因子的正则化边界元法

双材料界面裂纹渐近位移和应力场表现出剧烈的振荡特性, 许多用于表征经典平方根($r^{1/2})$和负平方根($r^{-1/2})$渐近物理场的传统数值方法失效, 给界面裂纹复应力强度因子($K_{1} +{i}K_{2} )$的精确求解增加了难度. 引入一种含有复振荡因子的新型"特殊裂尖单元", 可精确表征裂纹尖端渐近位移和应力场的振荡特性, 在避免裂尖区域高密度网格剖分的情况下, 可实现双材料界面裂纹复应力强度因子的精确求解. 此外, 结合边界元法中计算近奇异积分的正则化算法, 成功求解了大尺寸比(超薄)双材料界面裂纹的断裂力学参数. 数值算例表明, 所提算法稳定, 效率高, 在不增加计算量的前提下, 显著提高了裂尖近场力学参量和断裂力学参数的求解精度和数值稳定性.      

黏弹性界面断裂分析的增量“加料”有限元法

提出了一种适用于黏弹性界面裂纹问题的增量“加料” 有限元方法. 利用弹性界面裂纹尖端位移场的解答,通过对应原理和拉普拉斯逆变换近似方法,得到了黏弹性界面裂纹的尖端位移场. 用该位移场构造了黏弹性界面裂纹“加料” 单元和过渡单元位移模式,推导了增量“加料” 有限元方程,求解有限元方程可获得应力强度因子和应变能释放率等断裂参量. 建立了典型黏弹性界面裂纹平面问题“加料” 有限元模型,计算结果表明,对于弹性/黏弹性界面裂纹和黏弹性/黏弹性界面裂纹,该方法都能得到相当精确地断裂参量,并能很好地反映蠕变和松弛特性,可推广应用于黏弹性界面断裂问题的计算分析.      

光弹法实测正交异性双材料界面裂纹断裂参数

采用光弹贴片法实测正交异性双材料界面裂纹尖端区域的应力应变场, 进而求出界面裂纹的断裂力学参量. 将正交异性双材料板加工成拉伸试件,在聚碳酸酯贴片 的单侧表面镀金属铝膜,以提高贴片的反射效率. 沿贴片后的双材料界面预制裂缝,逐渐加 大载荷,得到一系列清晰的等差线条纹图. 利用正交异性双材料界面裂纹尖端应力分量表达 式计算出应力强度因子. 实验表明,光弹贴片法可有效地分析正交异性双材料界面裂纹问题.     

Interface crack problems for mode Ⅱ of double dissimilar orthotropic composite materials

The fracture problems near the interface crack tip for mode Ⅱ of double dissimilar orthotropic composite materials are studied. The mechanical models of interface crack for mode Ⅱ are given. By translating the governing equations into the generalized bi-harmonic equations,the stress functions containing two stress singularity exponents are derived with the help of a complex function method. Based on the boundary conditions,a system of non-homogeneous linear equations is found. Two real stress singularity exponents are determined be solving this system under appropriate conditions about himaterial engineering parameters. According to the uniqueness theorem of limit,both the formulae of stress intensity factors and theoretical solutions of stress field near the interface crack tip are derived. When the two orthotropic materials are the same,the stress singularity exponents,stress intensity factors and stresses for mode Ⅱ crack of the orthotropic single material are obtained.     

Interface crack problems for mode II of double dissimilar orthotropic composite materials

The fracture problems near the similar orthotropic composite materials are interface crack tip for mode Ⅱ of double disstudied. The mechanical models of interface crack for mode Ⅱ are given. By translating the governing equations into the generalized hi-harmonic equations, the stress functions containing two stress singularity exponents are derived with the help of a complex function method. Based on the boundary conditions, a system of non-homogeneous linear equations is found. Two real stress singularity exponents are determined be solving this system under appropriate conditions about bimaterial engineering parameters. According to the uniqueness theorem of limit, both the formulae of stress intensity factors and theoretical solutions of stress field near the interface crack tip are derived. When the two orthotropic materials are the same, the stress singularity exponents, stress intensity factors and stresses for mode II crack of the orthotropic single material are obtained.     

Effects of fractal crack

Experimental results indicate that propagation paths of cracks in geomaterials are often irregular, producing rough fracture surfaces which are fractal. In this paper, crack tip motion along a fractal crack trace is discussed. A fractal kinking model of the crack extension path is established to describe irregular crack growth. The length, velocity and kinking effects of the fractal crack are analysed. A formula is derived to describe the effects of fractal crack propagation on the dynamic stress intensity factor and on crack velocity. Finally, expressions of stress and displacement fields near the fractal crack tip are given.     

Dynamic stress intensity factor KⅢ and dynamic crack propagation characteristics of anisotropic materials

Based on the mechanics of anisotropic materials,the dynamic propagation problem of a mode Ⅲ crack in an infinite anisotropic body is investigated.Stress,strain and displacement around the crack tip are expressed as an analytical complex function,which can be represented in power series.Constant coefficients of series are determined by boundary conditions.Expressions of dynamic stress intensity factors for a mode Ⅲ crack are obtained.Components of dynamic stress,dynamic strain and dynamic displacement around the crack tip are derived.Crack propagation characteristics are represented by the mechanical properties of the anisotropic materials,i.e.,crack propagation velocity M and the parameter α.The faster the crack velocity is,the greater the maximums of stress components and dynamic displacement components around the crack tip are.In particular,the parameter α affects stress and dynamic displacement around the crack tip.     

Dynamic stress intensity factor K III and dynamic crack propagation characteristics of anisotropic materials

Based on the mechanics of anisotropic materials, the dynamic propagation problem of a mode Ⅲ crack in an infinite anisotropic body is investigated. Stress, strain and displacement around the crack tip are expressed as an analytical complex function, which can be represented in power series. Constant coefficients of series are determined by boundary conditions. Expressions of dynamic stress intensity factors for a mode Ⅲ crack are obtained. Components of dynamic stress, dynamic strain and dynamic displacement around the crack tip are derived. Crack propagation characteristics are represented by the mechanical properties of the anisotropic materials, i.e., crack propagation velocity M and the parameter ~. The faster the crack velocity is, the greater the maximums of stress components and dynamic displacement components around the crack tip are. In particular, the parameter α affects stress and dynamic displacement around the crack tip.     

Singularities of an interface crack in electrostrictive materials

In the present work, the singularities of an interface crack between two dissimilar electrostrictive materials under electric loads are investigated. Within the framework of two-dimensional deformation, the problem is solved using the complex variable method. Three crack models, that is, permeable, impermeable and conducting crack models are considered individually. Complex potentials and intensity factors of total stresses are derived by considering both the Maxwell stresses in the surrounding space at infinity and inside the crack. It is found that, for the above three crack models, the singularities of total stress are the same as those in traditional bi-materials with an interface crack; however, the intensities of the total stress depend on the actual crack model used.     

The scattering of general SH plane wave by interface crack between two dissimilar viscoelastic bodies

The scattering of general SH plane wave by an interface crack between two dissimilar viscoelastic bodies is studied and the dynamic stress intensity factor at the crack-tip is computed. The scattering problem can be decomposed into two problems: one is the reflection and refraction problem of general SH plane waves at perfect interface (with no crack); another is the scattering problem due to the existence of crack. For the first problem, the viscoelastic wave equation, displacement and stress continuity conditions across the interface are used to obtain the shear stress distribution at the interface. For the second problem, the integral transformation method is used to reduce the scattering problem into dual integral equations. Then, the dual integral equations are transformed into the Cauchy singular integral equation of first kind by introduction of the crack dislocation density function. Finally, the singular integral equation is solved by Kurtz's piecewise continuous function method. As a consequence, the crack opening displacement and dynamic stress intensity factor are obtained. At the end of the paper, a numerical example is given. The effects of incident angle, incident frequency and viscoelastic material parameters are analyzed. It is found that there is a frequency region for viscoelastic material within which the viscoelastic effects cannot be ignored. This work was supported by the National Natural Science Foundation of China (No.19772064) and by the project of CAS KJ 951-1-20     

黏弹性体界面裂纹的冲击响应

研究两半无限大黏弹性体界面Griffith裂纹在反平面剪切突出载荷下,裂纹尖端动应力强度因子的时间响应,首先,运用积分变换方法将黏弹性混合黑社会问题化成变换域上的对偶积分方程,通过引入裂纹位错密度函数进一步化成Cauchy型奇异积分方程,运用分片连续函数法数值求解奇异积分方程,得到变换域内的动应力强度因子,再用Laplace积分变换数值反演方法,将变换域的解反演到时间域内,最终求得动应力强度因子的时间响应,并对黏弹性参数的影响进行分析。