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

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

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

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

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

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

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

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

变幅载荷下三维裂纹体的疲劳寿命估算方法探讨

1.引言变幅载荷下三维裂纹扩展的疲劳寿命估算是十分麻烦的.而三维裂纹体的应力强度因子计算的复杂性,使估算工作更加困难.为了适应三维裂纹体的工程分析,本文采用了权函数.G.C.Sih;J R.Rice.P.M.Besuner和,T.A.Cruse等把权函数用于工程断裂分析中.由于他们采用的权函数只能对内埋或半椭圆表面裂纹的长、短轴方向进行分析,不能求解三维裂纹前沿任意点处的应力强度因子,故无法对三维裂纹进行较合理的疲劳分析.本文采用文[9]形式的权函数求解应力强度因子,并引入当量裂纹长度的概念,导出裂纹扩展过程中△K_I～△c的近似关系,进而采用广义Willenborg模型估算工程构件中表面裂纹在变幅载荷下的剩余疲劳寿命. 2.基本原理和方法     

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

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

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

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

三维热权函数法和多虚拟裂纹扩展技术

热权函数法直接利用温度场与热权函数的乘积的积分来求应力强度因子(SIF)的过渡过程，它可以免除对每一时刻进行热弹性力学有限元或边界元应力分析，计算效率大大提高。本文给出了三维热权函数法的基本方程，并提出了求解三维热权函数法基本方程的多虚拟裂纹扩展法(MVCE法)。在MVCE法中，可以引入无穷多个虚拟裂纹扩展模式；虚拟裂纹扩展模式与应力强度因子的插值直接相联系，所得到的方程组的系数矩阵是一个三对角矩阵，具有良好的计算性能。它对于裂纹前缘SIF分布急剧变化的情况，有良好的数值模拟能力。实例计算表明，MVCE法具有权高的计算效率，并具有很高的计算精度。     

围压与径向荷载共同作用下巴西盘裂纹应力强度因子的解析解

利用权函数法推导了围压和径向荷载共同作用下,考虑裂纹面摩擦的预制裂纹巴西盘应力强度因子计算公式,从理论上分析了围压、径向荷载和裂纹面摩擦对巴西盘应力强度因子的影响。结果表明,围压对I型应力强度因子有很大影响,I型应力强度因子随围压的增大而减小。当裂纹面闭合后围压和摩擦系数对II型应力强度因子同样具有显著影响,考虑裂纹面有效剪应力的权函数法理论解与有限元数值解相吻合,表明理论分析的正确性。     

围压与径向荷载共同作用下巴西盘裂纹应力强度因子的解析解

利用权函数法推导了围压和径向荷载共同作用下,考虑裂纹面摩擦的预制裂纹巴西盘应力强度因子计算公式,从理论上分析了围压、径向荷载和裂纹面摩擦对巴西盘应力强度因子的影响。结果表明,围压对I型应力强度因子有很大影响,I型应力强度因子随围压的增大而减小。当裂纹面闭合后围压和摩擦系数对II型应力强度因子同样具有显著影响,考虑裂纹面有效剪应力的权函数法理论解与有限元数值解相吻合,表明理论分析的正确性。     

Transient dynamic responses of a cracked solid subjected to in-plane loadings

Transient dynamic responses of an elastic cracked solid subjected to in-plane surface loadings are investigated in this study. Two vertical cracks, a surface-breaking crack and a sub-surface crack, are considered. The frequency responses of the plane strain problem are calculated by the computational mechanics combining the finite element method with the boundary integral equation. The finite element method is used for the near-field enclosing the crack, while the boundary integral equation is applied for the far-field to satisfy the Sommerfeld radiation condition. The transient responses are then obtained using fast Fourier transform. Surface displacements, crack opening displacements, and dynamic stress intensity factors are presented to show the significant effects of the cracks. The interaction between the elastic waves and the cracks as well as the mode conversion phenomena can be observed and understood through the numerical simulations.     

Mixed-mode fracture analysis of orthotropic functionally graded materials under mechanical and thermal loads

Mixed-mode fracture problems of orthotropic functionally graded materials (FGMs) are examined under mechanical and thermal loading conditions. In the case of mechanical loading, an embedded crack in an orthotropic FGM layer is considered. The crack is assumed to be loaded by arbitrary normal and shear tractions that are applied to its surfaces. An analytical solution based on the singular integral equations and a numerical approach based on the enriched finite elements are developed to evaluate the mixed-mode stress intensity factors and the energy release rate under the given mechanical loading conditions. The use of this dual approach methodology allowed the verifications of both methods leading to a highly accurate numerical predictive capability to assess the effects of material orthotropy and nonhomogeneity constants on the crack tip parameters. In the case of thermal loading, the response of periodic cracks in an orthotropic FGM layer subjected to transient thermal stresses is examined by means of the developed enriched finite element method. The results presented for the thermally loaded layer illustrate the influences of the material property gradation profiles and crack periodicity on the transient fracture mechanics parameters.     

Crack kinking in sandwich structures under three-point bending

Sandwich beams under three-point bending containing cracks in the core material very close to the upper skin interface are investigated. The cracks considered parallel or with an imperceptible inclination to the longitudinal beam axis and at different distances from the upper skin interface, are analyzed with static non-linear elastic two and three dimensional finite element analyses. From the proposed analyses stress intensity factors are calculated using linear elastic fracture mechanics. It is demonstrated that the crack propagation on the compression side of the core is mainly subjected in shear. The strain energy density criterion is used in order to determine the angle of kinking of the crack into the core.     

Multiple axial cracks in a coated hollow cylinder due to thermal shock

The transient thermal stress problem of an inner-surface-coated hollow cylinder with multiple pre-existing surface cracks contained in the coating is considered. The transient temperature, induced thermal stress, and the crack tip stress intensity factor (SIF) are calculated for the cylinder via finite element method (FEM), which is exposed to convective cooling from the inner surface. As an example, the material pair of a chromium coating and an underlying steel substrate 30CrNi2MoVA is particularly evaluated. Numerical results are obtained for the stress intensity factors as a function of normalized quantities such as time, crack length, convection severity, material constants and crack spacing.     

混凝土三维裂纹应力强度因子的简化算法

根据混凝土特性,本文提出了一种计算脆性材料三维裂缝应力强度因子的有限元方法,用以分析混凝土结构物中裂缝的稳定性及危害程度.其特点是不需要在裂缝尖端附近区域划分过细的网络,也不用引入特殊尖端单元,便于在工程中使用.本方法曾用于某拱坝坝肩裂缝稳定判别及应力重分布计算,效果良好.     

内部压力作用下矩形板中源于椭圆孔的分支裂纹应力强度因子的一种数值分析

应用一种边界元方法来研究内部压力作用下矩形板中源于椭圆孔的分支裂纹。该边界元方法由Crouch与Starfied建立的常位移不连续单元和笔者最近提出的裂尖位移不连续单元构成。在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界。本数值结果进一步证实这种数值方法对计算有限大板中复杂裂纹的应力强度因子的有效性,同时该数值结果可以揭示裂纹体几何对应力强度因子的影响。     

Three-dimensional mixed-mode stress intensity factors for cracks in functionally graded materials using enriched finite elements

Three-dimensional enriched finite elements are used to compute mixed-mode stress intensity factors (SIFs) for three-dimensional cracks in elastic functionally graded materials (FGMs) that are subject to general mixed-mode loading and constraint conditions. The method, which advantageously does not require special mesh configuration/modifications and post-processing of finite element results, is an enhancement of previous developments applied so far on isotropic homogeneous and isotropic interface cracks. The spatial variation of FGM material properties is taken into account at the level of element integration points. To validate the developed method, two- and three-dimensional mixed-mode fracture problems are selected from the literature for comparison. Two-dimensional cases include: inclined central crack in a large FGM medium under uniform tensile strain and stress loadings, a slanted crack in a finite-size FGM plate under exponentially varying tensile stress loading and an edge crack in a finite-size plate under shear traction load. The three-dimensional example models a deflected surface crack in a finite-size FGM plate under uniform tensile stress loading. Comparisons between current results and those from analytical and other numerical methods yield good agreement. Thus, it is concluded that the developed three-dimensional enriched finite elements are capable of accurately computing mixed-mode fracture parameters for cracks in FGMs.     

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.     

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.