多层材料结构的界面裂纹尖端复应力强度因子

本文建立了一种多层材料复合结构的界面裂纹问题分析模型。当两种材料之间插入第三种薄层弹性材料，裂纹位于第三种材料与第一或第二种弹性材料的界面上，且插页材料３＾＃的厚度相对于裂纹尺寸或平面内其他尺寸很小时，可以得到该问题裂纹尖端的复应力强度因子通式。本文用有限元法对结果进行了数值验证，并进行了有关问题的讨论。     

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

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

界面裂纹问题中的弹性T项和应力强度因子

研究两相材料有限板含单边界面裂纹的断裂力学特性，对不同的材料组合用广义变分法分析了不同尺寸试件和裂纹长度下的应力强度因子和弹性T项，讨论了材料特性对应力强度因子和弹性T项的作用．分析了试件尺寸和裂纹长度对应力强度因子和弹性T项的影响．     

薄膜涂层材料中币形界面裂纹的弹性波散射

研究了薄膜涂层材料中币形界面裂纹的弹性波散射问题,建立了含有币形界面裂纹的覆层半空间模型,采用Hankel积分变换,将裂纹对弹性波散射的问题转化为求解矩阵形式的奇异积分方程。结合渐近分析和围道积分技术得到积分方程的解,进一步推导了散射波的应力场和位移场,以及动应力强度因子的理论计算公式。在数值算例中,分析了不同材料组合和裂纹尺寸情况下动应力强度因子与入射波频率的关系,并给出了裂纹张开位移的结果。为薄膜涂层材料的动态破坏分析提供了一定的理论基础。     

热释电材料问题的通解与界面裂纹

该文讨论了热释电材料中的热弹性问题的一般解，进而求解了共线界面裂纹问题．利用Ｓｔｒｏｈ方法，把热释电材料的热弹性界面裂纹问题化为一向量形式的Ｈｉｌｂｅｒｔ问题，求出这一Ｈｉｌｂｅｒｔ问题的通解，进而求得了热释电材料热弹性界面裂纹的闭合解，得到了温度、热流、位移、电势、应力和电位移的全场解，得到了裂纹张开位移及电势差的精确表达式．在此基础上，还求得了均匀热释电体中单个热弹性裂纹裂尖场，单个界面裂纹裂尖场以及点热源与界面裂纹的作用．此外，该文还对界面裂纹顶点附近的端部场作了渐近分析．     

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

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

纤维端部的界面裂纹分析

基于弹性力学空间轴对称问题的通解,研究了短纤维增强复合材料中纤维端部的轴对称币形和柱形界面裂纹尖端的应力奇异性,得到了裂纹尖端附近的奇异应力场．研究结果表明,这两种轴对称界面裂纹尖端的应力奇异性相同,并且与平面应变状态下相应模型的应力奇异性一致,材料性能对裂纹尖端附近奇异应力场的影响可用三个组合参数描述     

不同材料界面共圆弧裂纹的反平面问题

本文首次将不同弹性材料界面共圆弧裂纹版平面问题,化为解析函数边值问题,获得了一般解答,由此求出了几种典型情况的精确解,算出了应力强度因子。当两种材料相同时,本文结果与文[S]完全吻合。     

界面裂纹尖端有限变形弹性场的渐近分析

采用完全非线性弹性理论，研究了一类新的可压缩超弹性材料形成的界面裂纹问题，给出了平面应变条件下裂尖场的渐近解．揭示了界面裂纹尖端场的变形特征．     

有限宽的粘接SANDWICH型正交各向异性板条具有界面裂纹或断裂的中间板条

本文研究了有限宽、粘接的对称SANDW(?)CH型正交各向异性板条的静裂纹问题.在中间板条有内部裂纹和完全断裂的两种情形,解法和应力奇异性分析的过程都和板条为各向同性时相似;但在界面裂纹时,却归结为解一组与各向同性粘接板条不同的二类柯西型奇异积分方程.此时,各向同性粘接板条界面裂纹的应力强度因子的定义已不再适用.本文提出一种广义的应力强度因子定义,并给出上述三种裂纹问题的算例,计算裂纹长度、板条宽度或弹性常数对应力强度因子的影响.     

One model of fatigue crack growth

A model for crack growth is proposed based on studies of the variation in the curvature radius at the crack tip during cyclic loading. Relations are obtained between mechanical material characteristics, crack geometry, and the rate of crack growth in a structure under cyclic loading. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 4, pp. 167–175, July–August, 2009.     

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.     

Comparative study of an interface crack for different wedge-interface models

Summary  An interface crack problem is investigated under various assumptions on an interface between two elastic materials. The interface is modeled by an additional third structure (thin elastic wedge of differing elastic properties) matching the bonded materials, or by introducing special boundary conditions on the crack line ahead. The main emphasis of the paper is placed on a comparison of the asymptotic expansion of the elastic solutions near the crack tip obtained for the different models. In particular, the behaviour of the stress singularity exponent and the generalized SIF are discussed. Numerical examples are presented. Received 16 August 2000; accepted for publication 26 May 2001     

Field structure at mode III dynamically propagating crack tip in elastic-viscoplastic materials

An elastic-viscoplastic mechanics model is used to investigate asymptotically the mode Ⅲ dynamically propagating crack tip field in elastic-viscoplastic materials. The stress and strain fields at the crack tip possess the same power-law singularity under a linear-hardening condition. The singularity exponent is uniquely determined by the viscosity coefficient of the material. Numerical results indicate that the motion parameter of the crack propagating speed has little effect on the zone structure at the crack tip. The hardening coefficient dominates the structure of the crack-tip field. However, the secondary plastic zone has little influence on the field. The viscosity of the material dominates the strength of stress and strain fields at the crack tip while it does have certain influence on the crack-tip field structure. The dynamic crack-tip field degenerates into the relevant quasi-static solution when the crack moving speed is zero. The corresponding perfectly-plastic solution is recovered from the linear-hardening solution when the hardening coefficient becomes zero.     

Simulation of size effects by a phase field model for fracture

In phase field fracture models the value of the order parameter distinguishes between broken and undamaged material. At crack faces the order parameter interpolates smoothly between these two states of the material, which can be regarded as phases. The crack evolution follows implicitly from the time integration of an evolution equation of the order parameter, which is coupled to the mechanical field equations. Among other phenomena phase field fracture models are able to reproduce crack nucleation in initially sound materials. For a 1D setting it has been shown that crack nucleation is triggered by the loss of stability of the unfractured, spatially homogeneous solution, and that the stability point depends on the size of the considered structure. This work numerically investigates to which extend size effects are reproduced by the 2D phase field model. Exemplarily, a finite element study of the hole size effect is performed and the simulation results are compared to experimental data.     

Field structure at mode Ⅲ dynamically propagating crack tip in elastic-viscoplastic materials

An elastic-viscoplastic mechanics model is used to investigate asymptotically the mode Ⅲ dynamically propagating crack tip field in elastic-viscoplastic materials. The stress and strain fields at the crack tip possess the same power-law singularity under a linear-hardening condition. The singularity exponent is uniquely determined by the viscosity coefficient of the material. Numerical results indicate that the motion parameter of the crack propagating speed has little effect on the zone structure at the crack tip. The hardening coefficient dominates the structure of the crack-tip field. However, the secondary plastic zone has little influence on the field. The viscosity of the material dominates the strength of stress and strain fields at the crack tip while it does have certain influence on the crack-tip field structure. The dynamic crack-tip field degenerates into the relevant quasi-static solution when the crack moving speed is zero. The corresponding perfectly-plastic solution is recovered from the linear-hardening solution when the hardening coefficient becomes zero.     

一种随机裂纹结构的裂纹扩展路径分析方法

裂纹结构中存在大量不确定性因素,如裂纹长度、材料性质、外部载荷等,裂纹扩展路径的不确定性分析对研究随机裂纹结构损伤和断裂的力学特性并预测其性能及可靠性具有重要意义。本文提出了一种适应于混合载荷模式下随机裂纹结构的裂纹扩展路径分析方法。该方法考虑了裂纹长度、材料性质和外部载荷等的随机性,并通过蒙特卡洛方法对随机参数空间进行采样。采用比例边界有限元方法计算结构应力强度因子,进而模拟单次裂纹扩展路径。在此基础上,通过概率分析方法获得随机裂纹结构中裂纹扩展路径的统计特性。最后给出了两个数值算例验证了本文方法的有效性。     

Anti-plane stress intensity, energy release and energy density at crack tips in a functionally graded strip with linearly varying properties

The fracture problem of a crack in a functionally graded strip with its properties varying in a linear form along the strip thickness under an anti-plane load is considered. The embedded anti-plane crack is located in the middle of strip half way through the thickness. The third mode stress intensity factor is derived using two different methods. In the first method, by employing Fourier integral transforms, the governing equation is converted to a singular integral equation, which is subsequently solved numerically by the collocation method based on Chebyshev polynomials. Then, the problem is solved by means of finite element method in which quadrilateral 8-node singular elements around each crack tip are used. After inspecting the validity of the solution technique, effects of crack geometry and non-homogeneous material parameter on the stress intensity, energy release and energy density are studied and the results of analytical and FEM solutions are compared.     

Mode III fracture problem of two arbitrarily oriented cracks located within two bonded functionally graded material strips

This paper shows the anti-plane crack problem of two bonded functionally graded material (FGM) strips. Each strip contains an arbitrarily oriented crack. The material properties of the strips are assumed in exponential forms varied in the direction normal to the interface. After employing the Fourier transforms, the unknowns are solved from the interface conditions, boundary conditions and the condition on the crack surfaces. The problem can then be reduced to a system of singular integral equations, which are solved numerically by applying the Gauss-Chebyshev integration formula to obtain the stress intensity factors at the crack tips. In the discussions, several degenerated problems are considered to demonstrate the influence of the non-homogeneous parameters, crack orientations, edge effects and the crack interactions on the normalized intensity factors. In general, the factors are larger when crack tips are located in stronger material. Also, the factors increase as the crack is oriented in the direction normal to the interface. The conclusions made in this research can be used to evaluate the safety of two bonded strips once the cracks exist inside the structure.