映射函数法求方形孔口角点裂纹的应力强度因子

用映身函数法和Ｍｕｓｋｈｅｌｉｓｈｖｉｌｉ方程对无限大板内方形孔口角点裂纹在不同外载下进行应力分析，决定了裂纹尖端的应力强度因子随裂纹长ａ和孔构形尺寸Ｌ的变化规律，探讨了外载变化对Ｋ因子的影响，并使刚度导数法（有限元法）进行了计算，结果表明，在板宽Ｗ和Ｌ之比等于或大于５时，两种方法的误差在６％以内。     

界面裂纹问题中的通用权函数法

热载荷和机械载荷共同作用下复合材料中的裂纹扩展往往发生在界面处.传统求解热冲击及机械载荷共同作用下界面裂纹尖端的应力强度因子的数值方法(如有限元、边界元法等),计算工作量大、效率低.通用权函数与时间无关,运用通用权函数法可以免除对每个时刻的应力分析,计算效率可得到很大提高.本文将通用权函数法推广到求解热载荷和机械载荷共同作用下界面裂纹尖端的应力强度因子过渡过程的问题中,推导出求解平面双材料界面裂纹问题应力强度因子的通用权函数法计算格式.基于此格式,计算热载荷和机械载荷共同作用下界面裂纹尖端的应力强度因子.通过实例计算比较,表明此方法得到的结果可以达到与相互作用积分法相当的工程应用精度.最后,应用此方法研究了热障涂层受热冲击及表面力共同作用时裂纹长度以及涂层厚度对应力强度因子的影响.结果表明:在一定边界条件下,当热障涂层中存在边缘裂纹时,随着涂层厚度的增加,更容易导致裂纹的扩展和涂层的剥落.     

裂纹体分析的权函数理论与应用: 回顾和展望

断裂力学是工程材料和结构的疲劳与断裂分析、损伤容限设计和结构完整性评定的理论基础. 应力强度因子作为线弹性裂纹尖端奇异场的单一表征参量和裂纹扩展驱动力, 在裂纹体的断裂力学分析中发挥着关键作用. 权函数法为复杂受载裂纹体的应力强度因子求解计算提供了强有力的解析工具, 不但具有远高于各类数值解法的计算效率, 而且精度可靠, 使用方便. 本文结合笔者团队在权函数法方面的长期研究工作, 对该方法自20世纪70年代初提出至今半个世纪以来, 国际断裂界在二维和三维权函数理论与应用方面的主要研究进展作了回顾和评述, 并对其未来发展提出了展望. 主要内容涵盖: 当前国际断裂界广泛应用的3种二维裂纹解析权函数法简介和以格林函数为基准的验证评价; 三维裂纹问题的片条合成权函数法和点载荷权函数法; 权函数法在复杂受载裂纹体的应力强度因子和裂纹张开位移等关键力学参量计算、内聚力/桥连等裂纹模型分析、共线多裂纹权函数理论及其在剩余强度预测等方面的应用, 以及复杂裂纹几何的工程化权函数分析和权函数法的反向应用问题.      

半无限板边缘裂纹的权函数解法与评价

权函数法是求解裂纹体在任意受载条件下的应力强度因子和裂纹面位移等断裂力学参量的高效、高精度方法,与有限元等数值方法相比,在求解效率和可靠性方面均具有明显优势.针对半无限板边缘裂纹,系统分析了在国际断裂力学界较有代表性的Wu-Carlsson、Glinka-Shen和Fett-Munz三种解析形式的权函数法,进而以在远端均匀加载下的半无限板边缘裂纹面位移Wigglesworth解析解导得的权函数及其对应的格林函数解(即裂纹面受一对单位集中力作用下的应力强度因子)为基准,沿整个裂纹长度对3种权函数的精度逐点进行比较,并与文献中基于其他方法求得的权函数做了广泛对比,包括Bueckner,Hartranft-Sih以及Wigglesworth利用不同解析方法推导出的高精度的权函数.研究了3种参考载荷(均布/正反向线性分布应力、集中力)及其不同组合,以及裂纹嘴位移的几何条件对权函数精度的影响.结果表明,基于一种参考载荷下的裂纹面张开位移比基于两种参考载荷下的应力强度因子所得到的权函数具有更高的精度,而且后一种方法的精度明显受到所选参考载荷组合的影响;裂纹面位移在裂纹嘴处三阶导数等于零的条件对基于一个参考解的权函数精度的改进效果较小.最后给出了利用各种权函数方法计算得到的4种载荷条件下的应力强度因子,并对结果进行了比较.     

含电饱和区的压电双材料反平面界面Yoffe型运动裂纹问题研究

本文利用Fourier变换及Copson求解方法,得到了压电双材料中Yoffe型运动裂纹在裂纹尖端含有条状电饱和区条件下的相关解析解.结果表明电饱和尺寸只与电荷载有关,而与裂纹扩展速度无关;裂纹尖端的应力强度因子及电势跳变不仅与荷载及材料参数有关,而且还受到速度的影响,其中应力强度因子随着速度的增大而增大,而电势跳变随着速度的变化呈现递减趋势.     

考虑表面效应时正三角形孔边裂纹反平面剪切问题的断裂力学分析

基于表面弹性理论和保角映射,研究了远场反平面剪切载荷作用下考虑表面效应时正三角形孔边裂纹问题的断裂性能。给出了孔边应力场解答,获得了裂纹尖端应力强度因子解析解答。数值算例讨论了应力强度因子随三角形孔尺寸、裂纹长度和表面性能的变化规律。结果表明：当三角形孔尺寸在在纳米量级时,无量纲应力强度因子受孔隙尺寸影响显著;随着三角形孔尺寸的增大,本文结果趋近于经典断裂理论解答;无量纲应力强度因子随裂纹长度的增加,数值先增大而后减小;裂纹相对长度较小时,表面效应影响较弱;应力强度因子的尺寸效应受表面性能影响显著。     

含孔边裂纹球壳的断裂分析

在依据Reissner理论得出的球壳裂纹尖端应力应变场展开式基础上,采用局部—整体分析法和权函数方法分别计算承受内压的含孔边裂纹球壳的应力强度因子.在有限元的模式中考虑剪切变形的影响,并对奇异元模式的应力应变场展开式的项数选择、奇异元最佳尺寸的选取进行了分析.本文计算和分析了在不同几何尺寸,不同开孔大小以及不同剪切参量条件下承受内压的含孔边裂纹球壳的应力强度因子及其变化规律.     

热冲击应力强度因子过渡过程分析的热权函数法

热权函数法利用温度与热权函数的乘积的积分来直接计算热冲击过程中裂纹尖端的应力强度因子过渡过程。热权函数与时间τ无关。由于免除了对每一时刻τ所需作的有限元或边界元分析,计算过程大大简化,计算效率得到极大的提高。本文将热权函数与有限元法直接耦合,给出了基于刚度阵导数法和Gauss-Chebyshev积分的热权函数计算格式。实例计算表明,本文给出的热权函数计算格式具有满意的计算精度。     

轴对称热冲击SIF过渡过程分析的热权函数法

热权函数法利用温度与热权函数的乘积的积分来直接计算热冲击过程中裂纹尖端的应力强度因子过渡过程。热权函数与时间τ无关,由于免除了对每一时刻τ所需作的有限元或边界元分析,计算过程大大简化,计算效率得到极大提高,本文将热权函数与有限元法直接耦合,给出了基于刚度阵导数法的轴对称问题的热权函数计算格式,实例计算表明,本文给出的热权函数计算格式具有满意的计算精度。     

裂纹面受荷载作用的应力强度因子的计算

基于比例边界有限元法计算了裂纹面有荷载作用情况下裂纹尖端的应力强度因子,给出了有限介质裂纹面作用荷载的比例边界有限元方程的基本求解过程.对于随径向坐标任意变化的一类面荷载的积分能够显式计算,不需要引入额外的近似;并将计算结果与解析解和数值结果进行对比,结果表明比例边界有限元法在计算裂纹面作用荷载时的应力强度因子是有效且精确的.此外,该方法可方便地处理各向异性材料裂纹问题,本文给出了正交各向异性矩形盘裂纹面受均布荷载情况的应力强度因子.     

WEIGHT FUNCTION FOR STRESS INTENSITY FACTORS IN ROTATING THICK-WALLED CYLINDER

The equation of stress intensity factors(SIF) of internally pressurized thick- walled cylinder was used as the reference case.SIF equation of rotating thick-walled cylinder containing a radial crack along the internal bore was presented in weight function method.The weight fumction formulas were worked out and can be used for all kinds of depth of cracks,rotating speed,material,size of thick-walled cylinder to calculate the stress intensity factors.The results indicated the validity and effectiveness of these formulas.Meanwhile,the rules of the stress intensity factors in rotating thick-walled cylinder with the change of crack depths and the ratio of outer radius to inner radius were studied.The studies are valuable to engineering application.     

Perturbation of a dynamic planar crack moving in a model viscoelastic solid

Weight functions, which give stress intensity factors in terms of applied loading, are constructed, for three-dimensional time-dependent loading of a semi-infinite crack, propagating at uniform speed. Both a model problem, governed by a scalar wave equation, and the full vectorial problem for Mode I loading, are considered. The medium through which the crack propagates is viscoelastic; the approach is general but explicit formulae are given when the medium is a Maxwell fluid. The weight functions are exploited to develop formulae for the first-order perturbations of stress intensity factors when the crack edge is no longer straight but becomes slightly wavy. Implications for stability, and for “crack front waves” in the case of the Mode I problem, are discussed.     

界面裂纹问题中的权函数方法

本文将Ｐａｒｉｓ等确定均匀材料中裂纹尖端应力强度因子的权函数方法推广应用到界面裂纹问题，给出了界面裂纹尖端附近或无限大体半无限界面裂纹问题的权函数的显式表达式。利用此权函数表达式可以很简便地求解界面裂纹尖端附近一些外来作用引起的应力强度因子，比如任意分布力、相变应变、位错和热等。作为一个算例，本文计算了界面一侧一个刃型位错引起的应力强度因子。     

Weight functions for a surface semi-elliptical crack in a finite thickness plate

Weight functions for the surface and the deepest point of a semi-elliptical crack in a finite thickness plate were derived from a general weight function and two references stress intensity factors. The weight functions were validated against finite element data by comparison of stress intensity factors calculated for several linear and non-linear stress fields. Differences were less than 2% for the surface point and 5% for the deepest point. The final weight functions are given in closed forms suitable for computer numerical integration. The weight functions appear to be particularly suitable for fatigue crack growth prediction of semi-elliptical cracks and fracture analysis of such cracks in complex stress fields.     

含圆形嵌体弹性平面中径向裂纹问题的超奇异积分方程方法

根据含圆形嵌体平面问题在极坐标下的弹性力学基本解,使用Betti互换定理,在有限部积分意义下将问题归结为两个以裂纹岸位移间断为基本未知量、对于Ⅰ型和Ⅱ型问题相互独立的超奇异积分方程,对含圆形嵌体弹性平面中的径向裂纹问题进行了研究.根据有限部积分原理,建立了问题的数值算法.计算结果表明,嵌体半径、裂纹位置及材料剪切弹性模量等都对裂纹应力强度因子具有较为明显的影响.     

Analysis of dynamic stress intensity factors of three-point bend specimen containing crack

A new formula is obtained to calculate dynamic stress intensity factors of the three-point bending specimen containing a single edge crack in this study. Firstly, the weight function for three-point bending specimen containing a single edge crack is derived from a general weight function form and two reference stress intensity factors, the coefficients of the weight function are given. Secondly, the history and distribution of dynamic stresses in uncracked three-point bending specimen are derived based on the vibration theory. Finally, the dynamic stress intensity factors equations for three-pointing specimen with a single edge crack subjected to impact loadings are obtained by the weight function method. The obtained formula is verified by the comparison with the numerical results of the finite element method (FEM). Good agreements have been achieved. The law of dynamic stress intensity factors of the three-point bending specimen under impact loadings varing with crack depths and loading rates is studied.     

Evaluation of stress intensity factors by multiple reference state weight function approach

The advantages of expressing stress intensity factors in terms of weight functions are widely appreciated. However, the main obstacle in the determination of weight functions has been the definition of the crack opening displacement (COD) field. There are several approaches currently used, the most common is assuming an expression to define COD in terms of the crack dimensions and stress state. A recent development in weight function application simplifies the tranditional stress intensity factor calculation. This development uses more than one reference stress intensity factor and associated stress field to eliminate the need to assume a COD field. This paper describes current application of COD for weight functions and explains the full advantage of adopting a multiple reference state (MRS) weight function approach.     

Analysis of dynamic stress intensity factors of three-point bend specimen containing crack

A new formula is obtained to calculate dynamic stress intensity factors of the three-point bending specimen containing a single edge crack in this study. Firstly, the weight function for three-point bending specimen containing a single edge crack is derived from a general weight function form and two reference stress intensity factors, the coefficients of the weight function are given. Secondly, the history and distribution of dynamic stresses in uncracked three-point bending specimen are derived based on the vibration theory. Finally, the dynamic stress intensity factors equations for three-pointing specimen with a single edge crack subjected to impact loadings are obtained by the weight function method. The obtained formula is verified by the comparison with the numerical results of the finite element method (FEM). Good agreements have been achieved. The law of dynamic stress intensity factors of the three-point bending specimen under impact loadings varing with crack depths and loading rates is studied.     

基于不连续场修正权的无网格断裂分析

本文采用了一种基于不连续场修正权函数的无网格方法来处理二维平面问题中的有限长裂纹。相较于目前常用的无网格裂纹不连续性处理方案,采用修正权函数处理裂纹附近不连续场时只需要对原权函数进行修正,算法简便易实现。本文采用基于不连续场修正权函数的无单元Galerkin方法（EFGM）,对在边界上施加I-II混合型裂纹位移场的斜裂纹板进行了数值分析。并与可视性准则、衍射法和透射法等不连续准则对比了裂尖位移场、应力场和应力强度因子解的数值精度。另外,本文还对这四种不连续准则形函数的计算效率进行了分析和比较。     

Pre-kinking of a moving crack in a magnetoelectroelastic material under in-plane loading

A constant moving crack in a magnetoelectroelastic material under in-plane mechanical, electric and magnetic loading is studied for impermeable crack surface boundary conditions. Fourier transform is employed to reduce the mixed boundary value problem of the crack to dual integral equations, which are solved exactly. Steady-state asymptotic fields near the crack tip are obtained in closed form and the corresponding field intensity factors are expressed explicitly. The crack speed influences the singular field distribution around the crack tip and the effects of electric and magnetic loading on the crack tip fields are discussed. The crack kinking phenomena is investigated using the maximum hoop stress intensity factor criterion. The magnitude of the maximum hoop stress intensity factor tends to increase as the crack speed increases.