采用显微网格法研究韧性材料平面应力状态下裂尖附近的应变场

用显微网格数字图象处理方法，测量了韧性材料平面应力条件下Ｉ型单边裂纹尖端附近的应变场。实验结果表明：在裂尖两侧存在与裂纹方向夹角约为５０°的大变形带。最大应变位于离裂尖一定距离的大变形带上，随着载荷增加，最大应变的位置逐渐靠近裂纹尖端。应变峰值轨迹线在裂纹顶端呈放射状，载荷增加时，应变峰值的轨迹线位置并不改变。     

使用焦散线法与光弹法测定三维裂纹混合型应力强度因子

本文介绍了将焦散线法与应力冻结，“解冻”技术相结合，使焦散线法用于测量三维体内部裂纹前缘应力强度因子的实验方法，并针对复杂应力状态下三维裂纹前缘的不同应变奇异场，合理地综合运用焦散线法与三维光弹法，实际测量与分离了三维裂纹前缘混合型应力强度因子。     

应变损伤材料I型动态扩展裂纹尖端场的渐近分析

研究了应变损伤材料Ｉ型动态扩展的裂纹尖端场。假定材料服从Ｊ２流动理论，且损伤规律以幂律应变软化的规律给出。对于塑性区引进了应力函数φ，ψ０借助于动力学方程的分析，给出了渐近方程及数值解。结果表明，对于可压缩材料Ｉ型平面应变尖端场是完全由塑性区组成，没有弹性卸载区。在裂纹尖端附近，应力和应变分别具有如下的奇异性：σ￣（ｌｎＲ／ｒ）＾－ｎ／ｎ＋１，ε￣（ｌｎＲ／ｒ）＾１／ｎ＋１。     

考虑裂纹闭合效应的岩石损伤本构关系

岩石中的预存裂纹只有在一定的法向压应力即裂纹闭合应力的作用下才可能闭合，其闭合过程与其方位和外加应力场有关，并且，即使对于裂纹已经完全闭合的岩石，如果裂纹闭合应力不同，则岩石的应力应交关系也不相同。本文建立了考虑裂纹闭合效应的岩石细观损伤力学模型，分析了裂纹闭合应力对岩石损伤演化过程和应力应变关系的影响。数值结果表明裂纹闭合应力显著地改变岩石的应力应变关系，表现为随裂纹闭合应力的增加，岩石的轴向应变变化较小，侧向应变和体积应变则大为增加。     

用激光散斑法和云纹法测量含裂纹薄板的裂纹张开位移

平面应力的弹塑性断裂理论的研究一直受到广泛重视,正在深入发展,其中很多工作是通过研究裂纹张开位移和裂纹附近的应变场进行的。用光学方法对裂纹张开位移和裂纹周围的变形进行测量具有显著的优越性,不但测量精度高,而且能够把整个裂纹包括裂纹顶端的张开位移都同时测量出来,并得到裂纹周围的位移全场分布。近年来激光散斑法在实验应力分析中发展很快,已被用来测量裂纹张开位移和裂纹周围的应变场,其     

含微裂纹弹性体的应力应变关系

本义建立了考虑裂纹闭合和裂纹表面摩擦影响的含微裂纹弹性体的应力应变关系，给出了柔度张量增量的显式表达式。对于二维平面应力和平面应变状态，给出了等效工程弹性系数。数值计算结果表明，裂纹闭合和裂纹面摩擦对裂纹体的应力应变关系和等效工程弹性系数有重要影响。     

含微裂纹弹性体的应力应变关系

本义建立了考虑裂纹闭合和裂纹表面摩擦影响的含微裂纹弹性体的应力应变关系，给出了柔度张量增量的显式表达式。对于二维平面应力和平面应变状态，给出了等效工程弹性系数。数值计算结果表明，裂纹闭合和裂纹面摩擦对裂纹体的应力应变关系和等效工程弹性系数有重要影响。     

扩展裂纹尖端塑性场

本文通过对幂硬化材料中平面应变Ｉ型裂纹的扩展过程进行精细的弹塑性有限元计算，给出扩展裂纹尖端附近环形区域内弹塑性场的分布。首次提出适用于扩展裂纹尖端环形区域的三项解。其中首项为ＨＲＲ奇异解；第二项反映三轴应力的强弱；第三项与ＨＲＲ奇异性项相比还含有线性项。并指出：扩展裂纹尖端环形区域弹塑性应力应变场的分布和强弱可由Ｊ－Ｑ－ｋ２三参量刻划，此结论适用于不同试样几何。不同材料硬化指数以及由小范围屈服至     

有限宽板孔边弹塑性变形测试

测量应变集中区域的弹塑性变形，对于研究材料损伤或微裂纹的产生，以预防宏观裂纹的产生及扩展具有重要的意义。为给航空发动机轮盘的损伤容限设计研究提供相关材料的力学实验参数，做了两个有限宽中心带孔试件的拉伸试验。试验中，采用单调逐级加、卸栽的循环加栽方式，并应用自动网格法和数字图像相关技术测量了孔用的全场位移分布。再应用最小二乘法将离散的位移分量拟合成二元多项式函数，求解出拉伸方向的应变分布，并给出孔边应力σ与孔边应变εy的关系曲线。     

裂纹纹尖弹塑性应力应变场的实验研究和数值分析

本文分别采用激光和白光DSCM(数字散斑相关测量)方法对聚碳酸脂材料的裂纹尖端弹塑性应力、应变场进行了实验研究，两种方法都取得了较好的结果,有限元计算采用了各向同性多线性随动强化的本构模型，较好地模拟了材料应力、应变曲线，并对裂纹尖端进行了有限元网格的细化,通过对实验结果和有限元计算结果的比较，可以发现两者具有相近的趋势。     

裂纹在冲击载荷作用下起裂的临界载荷面

采用有限元方法研究裂纹在I型短脉冲载荷作用下应力强度因子随时间的变化 ,用应力强度因子的初始上升时间Tr 对时间坐标无量纲化 ,对应力强度因子初始上升段进行曲线拟合 ,得到了上升段的曲线表达式。运用简单弹性梁理论和Lagrangian运动方程 ,获得载荷与时间对裂纹作用的关系式 ,结合有限元的结果 ,得到了上升时间Tr 的计算表达式 ,并进一步推出了裂纹在冲击载荷作用下起裂的临界载荷面。     

刚度微分法计算压电材料平面断裂问题

把计算应变能释放率的刚度微分法推广到压电材料平面断裂问题．在此基础上，利用压电材料平面断裂问题的有限元数值解作为真实场，用Ｓｏｓａ的平面问题裂端渐近解作为辅助场，由推广的交互Ｍ积分法求得了应力强度因子ＫＩ，ＫＩＩ和电位移强度因子ＫＩＶ．算例表明，计算结果与理论解符合得很好     

含孔von Karman板中非线性波散射与边值问题

基于ｖｏｎ Ｋａｒｍａｎ板大挠度弯曲理论,利用小参数摄动法,分析研究了含孔ｖｏｎＫａｒｍａｎ板的非线性波散射与动应力集中问题,其中一类可看成是薄板弯曲波动问题的控制方程。当有单频波入射时,由于弯曲应力与膜应力状态的非线性耦合,孔洞会产生高次谐波散射现象。建立了求解本问题的边界积分方程法,利用积分方程法交替求求这两类问题,最终可获得问题的近似分析解。     

A method for measuring mode I crack tip constraint under static and dynamic loading conditions

A novel experimental technique for measuring crack tipT-stress, and hence in-plane crack tip constraint, in elastic materials has been developed. The method exploits optimal positioning of stacked strain gage rosette near a mode I crack tip such that the influence of dominant singular strains is negated in order to determineT-stress accurately. The method is demonstrated for quasi-static and low-velocity impact loading conditions and two values of crack length to plate width ratios (a/W). By coupling this new method with the Dally-Sanford single strain gage method for measuring the mode I stress intensity factorK _I , the crack tip biaxiality parameter is also measured experimentally. Complementary small strain, static and dynamic finite element simulations are carried out under plane stress conditions. Time histories ofK _I andT-stress are computed by regression analysis of the displacement and stress fields, respectively. The experimental results are in good agreement with those obtained from numerical simulations. Preliminary data for critical values ofK _I and β for dynamic experiments involving epoxy specimens are reported. Dynamic crack initiation toughness shows an increasing trend as β becomes more negative at higher impact velocities.     

SIF-based fracture criterion for interface cracks

The complex stress intensity factor K governing the stress field of an interface crack tip may be split into two parts, i.e.,■ and s~(-iε), so that K = ■ s~(-iε), s is a characteristic length and ε is the oscillatory index. ■ has the same dimension as the classical stress intensity factor and characterizes the interface crack tip field. That means a criterion for interface cracks may be formulated directly with■, as Irwin(ASME J. Appl. Mech. 24:361–364, 1957) did in 1957 for the classical fracture mechanics. Then, for an interface crack,it is demonstrated that the quasi Mode I and Mode II tip fields can be defined and distinguished from the coupled mode tip fields. Built upon SIF-based fracture criteria for quasi Mode I and Mode II, the stress intensity factor(SIF)-based fracture criterion for mixed mode interface cracks is proposed and validated against existing experimental results.     

数字梯度敏感方法及其在航空透明件断裂力学中的应用

本文将数字梯度敏感方法用于航空透明件断裂力学问题研究。首先,基于透明材料的弹性-光学效应,建立了透明件应力状态与光线穿过透明件后的偏转角之间的关系。在平面应力假设下,利用最小二乘拟合建立了I型裂纹尖端应力强度因子与光学偏转角的关系。其次,通过数字梯度敏感方法搭建非接触光学测试平台,开展了带单边裂纹的航空有机玻璃试件三点弯曲实验,应用数字梯度敏感方法提取了I型裂纹尖端应力强度因子。最后,通过选择不同计算子区域和步长大小,分析了数字梯度敏感方法中的子区域和步长选择对计算结果的影响。研究结果表明,数字梯度敏感方法实验所得应力强度因子与经验公式计算所得结果偏差小于10%,通过增加最小二乘拟合项数以及合理的子区域和步长选择可以减小数字梯度敏感方法计算应力强度因子误差。     

求解混合型裂纹应力强度因子的围线积分法

本文用复变函数理论推导出裂纹的辅助场，并用Betti功互等定理给出求解混合型裂纹应力强度因子的远场围绕积分法．此方法与积分路径的选择无关，用有限元法计算出远离裂纹尖端的位移场和应力场，就可通过计算绕裂端的围线积分，精确地给出混合型裂纹的应力强度因子KⅠ和KⅡ的数值解．     

有限厚度板穿透裂纹前缘附近三维弹性应力场分析

通过三维有限元计算来研究有限宽度、有限厚度含有穿透裂纹板的裂纹前缘应力场，从中找出应力强度因子与板的厚度、裂纹长度之间的关系，同时还分析了裂尖的三维约束程度和三维约束区的大小。分析结果表明：应力强度因子沿厚度的分布是不均匀的，应力强度因子的最大值及其位置与厚度有关；有限厚度板中面应力强度因子(KI)m-p及最大应力强度因子(KI)max均大于平面应力或平面应变的应力强度因子。对有限厚度裂纹问题，按平面应力或平面应变来考虑是不安全的；板中面的应力强度因子(KI)m-p及最大应力强度因子(KI)max是厚度B／a的函数；板的中面离面约束系数Tx最大，自由面(z=B)Tx=0。沿厚度方向裂尖附近的离面约束系数Tx也是z／B和B／a的函数，随着厚度的增加离面约束系数Tx增大，离中面越近离面约束系数Tx越大。Tx随着x的增大急剧减小，三维约束影响区域大小大约为板厚的一半，且裂纹长度a／W对应力强度因子沿厚度变化规律及Tx影响区域大小影响较小。     

Mode I and mode II crack tip asymptotic fields with strain gradient effects

The strain gradient effect becomes significant when the size of fracture process zone around a crack tip is comparable to the intrinsic material lengthl, typically of the order of microns. Using the new strain gradient deformation theory given by Chen and Wang, the asymptotic fields near a crack tip in an elastic-plastic material with strain gradient effects are investigated. It is established that the dominant strain field is irrotational. For mode I plane stress crack tip asymptotic field, the stress asymptotic field and the couple stress asymptotic field can not exist simultaneously. In the stress dominated asymptotic field, the angular distributions of stresses are consistent with the classical plane stress HRR field; In the couple stress dominated asymptotic field, the angular distributions of couple stresses are consistent with that obtained by Huang et al. For mode II plane stress and plane strain crack tip asymptotic fields, only the stress-dominated asymptotic fields exist. The couple stress asymptotic field is less singular than the stress asymptotic fields. The stress asymptotic fields are the same as mode II plane stress and plane strain HRR fields, respectively. The increase in stresses is not observed in strain gradient plasticity for mode I and mode II, because the present theory is based only on the rotational gradient of deformation and the crack tip asymptotic fields are irrotational and dominated by the stretching gradient. The project supported by the National Natural Science Foundation of China (19704100), National Natural Science Foundation of Chinese Academy of Sciences (KJ951-1-20), CAS K.C. Wong Post-doctoral Research Award Fund and Post-doctoral Science Fund of China     

Three-dimensional elliptic crack under impact loading

The dynamic stress intensity factor of a three-dimensional elliptic crack under impact loading is determined with the finite element method. The computation results can take into account the influence of time and the ratio of the wave speeds on the stress intensity factor. The present method is suitable not only for three-dimensional dynamic crack, but also for three-dimensional dynamic contact. Project supported by the National Natural Science Foundation of China (No. K19672007).