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

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

基于权函数法的核主泵飞轮应力强度因子研究

采用权函数法确定了含裂纹飞轮在离心力和接触压力作用下的应力强度因子,计算了在同步转速工况下裂纹尖端应力强度因子的值,并与解析法和有限元法计算结果进行了对比。结果表明:权函数法与解析法的误差在3%以内,与有限元法的误差在1%以内,验证了权函数法计算应力强度因子的准确性高;在结构不变的情况下,权函数法可以快速求解不同载荷条件、不同长度裂纹的应力强度因子。通过控制变量法研究了不同参数对应力强度因子的影响,结果表明,飞轮裂纹尖端总应力强度因子随着裂纹尺寸、旋转转速、飞轮尺寸外径与内径比值的增大而增大。     

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

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

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

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

应用通用权函数法计算热机载荷作用下三维界面裂纹的应力强度因子

热机载荷共同作用下双材料、复合材料中的裂纹扩展往往发生在界面处,并且工程中实际遇到的裂纹大多数是三维裂纹,传统的求解热冲击及机械载荷共同作用下界面裂纹应力强度因子的数值方法如有限元、边界元法计算量大,计算效率低。由于通用权函数仅仅与裂纹体的几何形状有关,与载荷、时间无关,求解应力强度因子时避免了反复的应力分析,计算效率大大提高, 通用权函数法非常适合计算复杂冲击载荷下应力强度因子分布的过渡过程。根据Betti互易原理,本文推导出了三维界面裂纹问题通用权函数法的普遍表达式,并给出了热机载荷共同作用下三维界面I型、Ⅱ型和Ⅲ型裂纹问题通用权函数法的有限元格式. 通过实例计算比较,表明此方法得到的结果可以达到满意的工程应用精度。     

压缩荷载作用下多裂纹的相互干涉

结合裂面边界条件和位移单值条件得到裂面上有伪集中力作用时应力函数的复势解.同时利用"伪力法",在叠加原理的基础上,借助Chebyshev数值积分法求得"伪力"的大小,给出了压缩载荷作用下多裂纹相互干涉的解法,讨论了裂纹之间的干涉效应.算例计算了两个裂纹之间的干涉及共线裂纹的相互作用,并分析了裂纹间的距离、夹角等因素对应力强度因子的影响.     

含径向裂纹系的圆柱的弯曲与扭转

圆柱中任意径向裂纹系的Saint-Venant弯曲与扭转,还无一般的解析解法,本文采用裂纹面二侧应力差和位移差的混合边界条件提法,在求解了一组三节积分方程和一个Neumann问题后,精确地求得到了单裂纹基本解,利用此解给出了解决这类问题的一般方法。文中对二条非共线的等长边界裂纹的应力强度因子和柱的扭转刚度作了数值计算,扭转刚度的结果与L.A.Wigglesworth的完全一致。     

热机载荷作用下三维界面裂纹问题的权函数法分析

热机载荷共同作用下双材料和复合材料中的裂纹扩展往往发生在界面处,并且工程中实际遇到的裂纹大多数是三维裂纹.由于通用权函数仅仅与裂纹体的几何形状有关,与载荷、时间无关,因此在求解复杂冲击载荷下界面裂纹应力强度因子随时间的变化过程时,避免了反复的应力分析,计算效率得到提高.根据Betti互易原理,论文推导出三维界面裂纹问题通用权函数法的普遍表达式,并给出了热机载荷共同作用下三维界面Ⅰ型、Ⅱ型和Ⅲ型裂纹问题通用权函数法的有限元格式.通过与实例计算比较,表明此方法得到的结果可以达到满意的工程应用精度.     

动态扩展裂纹的若干反平面问题的研究

采用复变函数论,对反平面条件下的动态裂纹扩展问题进行研究。通过自相似函数的方法可以获得解析解的一般表达式。应用该法可以很容易地将所讨论的问题转化为Riemann—Hilbert问题,并可以相当简单地得到问题的闭合解。文中分别对裂纹面受均布载荷、坐标原点受集中增加载荷、坐标原点受瞬时冲击载荷以及裂纹面受运动集中载荷Px／t作用下的动态裂纹扩展问题进行求解,得到了裂纹扩展位移、裂纹尖端的应力和动态应力强度因子的解析解。应用该解并通过叠加原理,就可以求得任意复杂问题的解。     

压电体内孔边裂纹的应力强度因子

本文研究含有Ⅲ型孔边裂纹压电弹性体的反平面问题.根据Muskhelishvili的数学弹性力学理论,并利用保角变换和Cauchy积分的方法,对含有圆孔孔边单裂纹和双裂纹的压电弹性体分别进行了分析.基于电不可穿透裂纹模型,得到了在反平面剪力和面内电载荷的共同作用下裂纹尖端应力强度因子的解析解.最后,通过数值算例,讨论了应力强度因子随裂纹长度变化的规律.结果表明:应力强度因子随着裂纹和孔的相对尺寸的增加而增加,并且单边裂纹的应力强度因子要比双边裂纹的应力强度因子大.     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

DETERMINATION OF CRACK SURFACE DISPLACEMENTS FOR A RADIAL CRACK EMANATING FROM A SEMI-CIRCULAR NOTCH USING WEIGHT FUNCTION METHOD

A radial crack emanating from a semi-circular notch is of significant engineering importance. Accurate determination of key fracture mechanics parameters is essential for damage tolerance design and fatigue crack growth life predictions. The purpose of this paper is to provide an efficient and accurate closed-form weight function approach to the calculation of crack surface displacements for a radial crack emanating from a semi-circular notch in a semi-infinite plate.Results are presented for two load conditions: remote applied stress and uniform stress segment applied to crack surfaces. Based on a correction of stress intensity factor ratio, highly accurate analytical equations of crack surface displacements under the two load conditions are developed by fitting the data obtained with the weight function method. It is demonstrated that the WuCarlsson closed-form weight functions are very efficient, accurate and easy-to-use for calculating crack surface displacements for arbitrary load conditions. The method will facilitate fatigue crack closure and other fracture mechanics analyses where accurate crack surface displacements are required.     

General solutions of plane problem for power function curved cracks

A new exact and universal conformal mapping is proposed. Using Muskhelishvili's complex potential method, the plane elasticity problem of power function curved cracks is investigated with an arbitrary power of a natural number, and the general solutions of the stress intensity factors (SIFs) for mode I and mode II at the crack tip are obtained under the remotely uniform tensile loads. The present results can be reduced to the well-known solutions when the power of the function takes different natural numbers. Numerical examples are conducted to reveal the effects of the coefficient, the power, and the projected length along the x-axis of the power function curved crack on the SIFs for mode I and mode II.     

压电效应下一维六方准晶双材料孔边裂纹的反平面问题研究

利用复变函数方法和保角变换技术研究了压电效应下一维六方准晶双材料中圆孔边单裂纹的反平面问题．考虑电不可渗透型边界条件,运用保角变换和Stroh公式得到了弹性体受远场剪切力和面内电载荷作用下裂纹尖端应力强度因子和能量释放率的解析解. 数值算例分析了几何参数、远场受力、电位移载荷对能量释放率的影响．结果表明：裂纹长度、耦合系数和远场剪切力的减小可以抑制裂纹的扩展．不考虑电场时,声子场应力对能量释放率的影响较小．本文的研究结果可作为研究一维六方压电准晶双材料孔边裂纹问题的理论基础,同时为压电准晶及其复合材料的设计、制备、优化和性能评估提供理论依据．     

Complex variable approach in studying modified polarization saturation model in two-dimensional semipermeable piezoelectric media

A modified polarization saturation model is proposed and addressed mathematically using a complex variable approach in two-dimensional(2 D) semipermeable piezoelectric media. In this model, an existing polarization saturation(PS) model in 2D piezoelectric media is modified by considering a linearly varying saturated normal electric displacement load in place of a constant normal electric displacement load, applied on a saturated electric zone. A centre cracked infinite 2D piezoelectric domain subject to an arbitrary poling direction and in-plane electromechanical loadings is considered for the analytical and numerical studies. Here, the problem is mathematically modeled as a non-homogeneous Riemann-Hilbert problem in terms of unknown complex potential functions representing electric displacement and stress components. Having solved the Hilbert problem, the solutions to the saturated zone length, the crack opening displacement(COD), the crack opening potential(COP), and the local stress intensity factors(SIFs) are obtained in explicit forms. A numerical study is also presented for the proposed modified model, showing the effects of the saturation condition on the applied electrical loading, the saturation zone length, and the COP. The results of fracture parameters obtained from the proposed model are compared with the existing PS model subject to electrical loading, crack face conditions, and polarization angles.     

The thermoelastic equilibrium of a transversally isotropic medium with an elliptic crack under symmetric loading

The stress intensity factors (SIFs) are evaluated for flat elliptical cracks located in a transversally isotropic material (cracks are assumed perpendicular to the transtropy axis) under an arbitrary load and symmetric temperature. The SIFs for an elliptical crack in a transversally isotropic medium are determined using the formulas (derived by the author in his previous studies) of transition from an isotropic to transversally isotropic material and the relative problem for an isotropic medium. It is proved that these formulas can be employed for an arbitrary homogeneous transversally isotropic material (no matter whether the roots of some characteristic equation of the material are real or complex) with an arbitrary flat crack or a system of coplanar flat cracks, including elliptical ones, under an arbitrary load and symmetric temperature. A transversally isotropic material with two coplanar elliptical cracks is considered as an illustrative example. The dependences of the SIFs on the parameters of cracks and their arrangement at a decreasing temperature are presented. S. P. Timoshenko Institute of Mechanics, National Academy of Sciences of Ukraine Kiev. Translated from Prikladnaya Mekhanika, Vol. 36, No. 4, pp. 96–105, April, 2000.     

Controlling Stress Intensity Factors During a Fatigue Crack Propagation Using Digital Image Correlation and a Load Shedding Procedure

Digital Image Correlation is used for controlling load shedding fatigue crack propagation. A specific algorithm is used to perform Stress Intensity Factors (SIFs) and crack length estimation in real time. Crack length measurements are validated by comparison with potential drop technique. SIFs results are compared with more common techniques using standard analytical formula considering confined plasticity at the crack tip. The proposed non-contact method is shown to be a powerful tool to control crack propagation.     

Elastic wave scattering from a penny-shaped interface crack in coated materials

This paper is concerned with the elastic wave scattering induced by a penny-shaped interface crack in coated materials. Using the integral transform, the problem of wave scattering is reduced to a set of singular integral equations in matrix form. The singular integral equations are solved by the asymptotic analysis and contour integral technique, and the expressions for the stress and displacement as well as the dynamic stress intensity factors (SIFs) are obtained. Using numerical analysis, this approach is verified by the finite element method (FEM), and the numerical results agree well with the theoretical results. For various crack sizes and material combinations, the relations between the SIFs and the incident frequency are analyzed, and the amplitudes of the crack opening displacements (CODs) are plotted versus incident wavenumber. The investigation provides a theoretical basis for the dynamic failure analysis and nondestructive evaluation of coated materials.     

Measurement of energy release rate and energy flux of rapidly bifurcating crack in Homalite 100 and Araldite B by high-speed holographic microscopy

High-speed holographic microscopy is applied to take three successive photographs of fast propagating cracks in Homalite 100 or in Araldite B at the moment of bifurcation. Crack speed at bifurcation is about 540 m/s on Homalite 100, and about 450 m/s on Araldite B. From the photographs, crack speeds immediately before and after bifurcation are obtained, and it is found that discontinuous change of crack speed does not exist at the moment of bifurcation in the case of Homalite 100, but exists in the case of Araldite B. From the photographs, crack opening displacement (COD) is also measured along the cracks as a function of distance r from the crack tips. The measurement results show that the CODs are proportional to √r before bifurcation. After bifurcation, the CODs of mother cracks are proportional to √r, though the CODs of branch cracks are not always proportional to √r. The energy release rate is obtained from the measured CODs, and it is found that energy release rate is continuous at bifurcation point in both cases of Homalite 100 and Araldite B. Energy flux that shows the energy flow toward a crack tip is also obtained.