扩展边界元法分析切口和裂纹结构应力场的准确性

在线弹性理论中,切口/裂纹结构尖端区域存在奇异应力场,数值方法不易求解.论文建立的扩展边界元法(XBEM)对围绕尖端区域位移函数采用自尖端径向距离r的渐近级数展开式表达,其级数项的幅值系数作为基本未知量,而外部区域采用常规边界元法离散方程.两者方程联立求解可获得切口和裂纹结构完整的位移和应力场.扩展边界元法具有半解析法特征,适用于一般的切口和裂纹结构应力场分析,其解可精细描述从尖端区域到整体结构区域的应力场.作者研制了扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明扩展边界元法求解切口和裂纹结构应力场的准确性和有效性.     

ANALYSIS OF 3-D NOTCHED/CRACKED STRUCTURES BY USING EXTENDED BOUNDARY ELEMENT METHOD 1)

在线弹性理论中,三维 V 形切口/裂纹结构尖端区域存在多重应力奇异性,常规数值方法不易求解. 本文提出和建立了三维扩展边界元法 (XBEM),用于分析三维线弹性 V 形切口/裂纹结构完整的位移和应力场. 先将三维线弹性 V 形切口/裂纹结构分为尖端小扇形柱和挖去小扇形柱后的外围结构. 尖端小扇形柱内的位移函数采用自尖端径向距离 $r$ 的渐近级数展开式表达,其中尖端区域的应力奇异指数、位移和应力特征角函数通过插值矩阵法获得. 而级数展开式各项的幅值系数作为基本未知量. 挖去扇形域后的外围结构采用常规边界元法分析. 两者方程联立求解可获得三维 V 形切口/裂纹结构完整的位移和应力场,包括切口/裂纹尖端区域精细的应力场. 扩展边界元法具有半解析法特征,适用于一般三维 V 形切口/裂纹结构完整位移场和应力场的分析,其解可精细描述从尖端区域到整体结构区域的完整应力场. 作者研制了三维扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明了扩展边界元法求解三维 V 形切口/裂纹结构完整应力场的准确性和有效性.     

三维切口/裂纹结构的扩展边界元法分析

在线弹性理论中,三维 V 形切口/裂纹结构尖端区域存在多重应力奇异性,常规数值方法不易求解. 本文提出和建立了三维扩展边界元法 (XBEM),用于分析三维线弹性 V 形切口/裂纹结构完整的位移和应力场. 先将三维线弹性 V 形切口/裂纹结构分为尖端小扇形柱和挖去小扇形柱后的外围结构. 尖端小扇形柱内的位移函数采用自尖端径向距离 $r$ 的渐近级数展开式表达,其中尖端区域的应力奇异指数、位移和应力特征角函数通过插值矩阵法获得. 而级数展开式各项的幅值系数作为基本未知量. 挖去扇形域后的外围结构采用常规边界元法分析. 两者方程联立求解可获得三维 V 形切口/裂纹结构完整的位移和应力场,包括切口/裂纹尖端区域精细的应力场. 扩展边界元法具有半解析法特征,适用于一般三维 V 形切口/裂纹结构完整位移场和应力场的分析,其解可精细描述从尖端区域到整体结构区域的完整应力场. 作者研制了三维扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明了扩展边界元法求解三维 V 形切口/裂纹结构完整应力场的准确性和有效性.      

Computations of modes I and II stress intensity factors of sharp notched plates under in-plane shear and bending loading by the fractal-like finite element method

The fractal-like finite element method (FFEM) is used to compute the stress intensity factors (SIFs) for different configurations of cracked/notched plates subject to in-plane shear and bending loading conditions. In the FFEM, the large number of unknown variables in the singular region around a notch tip is reduced to a small set of generalised co-ordinates by performing a fractal transformation using global interpolation functions. The use of exact analytical solutions of the displacement field around a notch tip as the global interpolation functions reduces the computational cost significantly and neither post-processing technique to extract SIFs nor special singular elements to model the singular region are required. The results of numerical examples of various configurations of cracked/notched plates are presented and validated via published data. Also, new results for cracked/notched plate problems are presented. These results demonstrate the accuracy and efficiency of the FFEM to compute the SIFs for notch problems under in-plane shear and bending loading conditions.     

Determination of V-notch SIFs in multi-material anisotropic wedges by digital correlation experiments

In this paper, theoretical formulations based on the Stroh’s complex function approach were used to find the displacement field and H-integral of a sharp V-notch formed from several anisotropic materials. Displacements from the image-correlation experiments are then substituted into the least-squares formulation to find V-notch stress intensity factors (SIFs) in multi-material anisotropic wedges. Validations using the H-integral indicate that the experimental SIFs evaluated from the proposed method of acceptable accuracy. The major advantage is that the proposed method only requires displacements inside the specimen, and displacements near the notch tip, specimen boundaries, or notch surfaces are not necessary.     

METHOD TO CALCULATE STRESS INTENSITY FACTOR OF V-NOTCH IN BI-MATERIALS

Based on Zak＇s stress function, the eigen-equation of stress singularity ofbi-materials with a V-notch was obtained. A new definition of stress intensity factor for a perpendicular interfacial V-notch of bi-material was put forward. The effects of shear modulus and Poisson＇s ratio of the matrix material and attaching material on eigen-values were analyzed. A generalized expression for calculating/（i of the perpendicular V-notch of bi-materials was obtained by means of stress extrapolation. Effects of notch depth, notch angle and Poisson＇s ratio of materials on the singular stress field near the tip of the V-notch were analyzed systematically with numerical simulations. As an example, a finite plate with double edge notches under uniaxial uniform tension was calculated by the method presented and the influence of the notch angle and Poisson＇s ratio on the stress singularity near the tip of notch was obtained.     

两相材料V形切口应力强度因子边界元分析

建立了边界元法计算两相材料粘结V形切口奇异应力场的新途径。在V形切口尖端挖出一小扇形,将该扇形弧线边界的位移和面力表示为有限项奇性指数和特征角函数的线性组合,其组合系数即为广义应力强度因子,将该组合回代到在被挖去小扇形后的剩余结构内建立的边界积分方程,离散后可求解出组合系数,获得两相材料粘结V形切口尖端的应力强度因子。算例证明了本文方法的有效性。     

V形切口应力强度因子的一种边界元分析方法

将V形切口结构分成围绕切口尖端的小扇形和剩余结构两部分. 尖端处扇形域应力场表示成关于尖端距离$\rho$的渐近级数展开式,从线弹性理论方程推导出了一组分析平面V形切口奇异性的常微分方程特征值问题,通过求解特征方程,得到前若干个奇性指数和相应的特征向量. 再将切口尖端的位移和应力表示为有限个奇性阶和特征向量的组合. 然后用边界元法分析挖去小扇形后的剩余结构. 将位移和应力的线性组合与边界积分方程联立,求解获得切口根部区域的应力场、应力幅值系数和整体结构的位移和应力. 从而准确计算出平面V形切口的奇异应力场和应力强度因子.      

基于奇异强度因子的V型切口应力场和N-SIF评估方法

根据线弹性断裂力学理论,V形切口处的应力场具有奇异性,应力值趋于无穷大,峰值应力不能直接用于评定疲劳强度。通过引入了奇异强度因子“as”,单边缺口应力分布和缺口应力强度因子(N-SIF)的半解析公式被推导。考虑张开角和几何尺寸等因素,基于奇异强度因子拟合得到了切口应力评估的简易公式,可用于切口应力场和N-SIF值的快速评估。将简易公式评估结果与有限元结果以及传统文献结果进行对比分析,结果表明,本文简易公式可以准确地预报拉伸载荷下单边V型切口角平分线上的应力场和N-SIF值,实现了切口试样应力场的快速评估。     

Fracture investigation at V-notch tip using coherent gradient sensing (CGS)

Local deformation field and fracture characterization of mode I V-notch tip are studied using coherent gradient sensing (CGS). First, the governing equations that relate to the CGS measurements and the elastic solution at mode I V-notch tip are derived in terms of the stress intensity factor, material constant, notch angle and fringe order. Then, a series of CGS fringe patterns of mode I V-notch are simulated, and the effects of the notch angle on the shape and size of CGS fringe pattern are analyzed. Finally, the local deformation field and fracture characterization of mode I V-notch tip with different V-notch angles are experimentally investigated using three-point-bending specimen via CGS method. The CGS interference fringe patterns obtained from experiments and simulations show a good agreement. The stress intensity factor obtained from CGS measurements shows a good agreement with finite element results under K-dominant assumption.     

Equivalence of the notch stress intensity factor,tip opening displacement and energy release rate for a sharp V-notch

For an infinite elastic plane with a sharp V-notch under the action of symmetrically loading at infinity, the length of crack initiation ahead of the V-notch’s tip is estimated according to a modified Griffith approach. Applying a new conservation integral to the perfectly plastic strip (Dugdale model) ahead of the V-notch’s tip, the relationship between notch stress intensity factor (NSIF) and notch tip opening displacement (NTOD) is presented. Also, the relationship between NSIF and perfectly plastic strip size (PPSS) is found. Since there are three fracture parameters (NSIF, NTOD, and PPSS) with changeable notch opening angle in two basic relationships, it is necessary to select one critical parameter with changeable notch opening angle or two independent critical parameters, respectively. With the help of a characteristic length, it is found by this new conservation integral that the NSIF, NTOD and energy release rate are equivalent in the case of small-scale yielding. Especially, the characteristic length possesses clear physical meaning and, for example, depends on both the critical NSIF and SIF or both the NTOD and CTOD, respectively, in which SIF and CTOD are from the tip of a crack degenerated from the sharp V-notch. The dependence of NSIF on NTOD and PPSS is presented according to the equivalence, and the critical NSIF depending on the critical NTOD with a notch opening angle is also predicted.     

COMPUTATION OF STRESS INTENSITY FACTORS BY THE SUB-REGION MIXED FINITE ELEMENT METHOD OF LINES

Based on the sub-region generalized variational principle,a sub-region mixed ver- sion of the newly-developed semi-analytical‘finite element method of lines’(FEMOL)is pro- posed in this paper for accurate and efficient computation of stress intensity factors(SIFs)of two-dimensional notches/cracks.The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used,with the sought SIFs being among the unknown coefficients.The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements.A mixed system of ordinary differential equations(ODEs) and al- gebraic equations is derived via the sub-region generalized variational principle.A singularity removal technique that eliminates the stress parameters from the mixed equation system even- tually yields a standard FEMOL ODE system,the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver.A number of numerical examples,including bi-material notches/cracks in anti-plane and plane elasticity,are given to show the generally excellent performance of the proposed method.     

Reliability assessment of a bi-material notch: Strain energy density factor approach

Joints of different materials have many applications in structural engineering and microelectronics. In the present contribution the joint is modelled as a bi-material notch. The singular stress field near the notch tip is investigated. Depending on the notch geometry and materials, the stress field can have one or two singularities. It is shown that to study the problem of a crack onset at the notch, both terms have to be taken into account. Criteria for the direction and for crack nucleation are formulated. The approach utilizes the knowledge of the strain energy density factor distribution in a bi-material notch vicinity.     

辛解析奇异元在动荷载作用下V型切口问题中的应用

利用辛解析奇异单元,结合时域精细算法,研究了动荷载作用下的含平面V型切口问题。时域上,采用时域精细算法,并结合自适应算法控制展开项数,保证了计算精度。空间域上,切口尖端附近采用辛解析奇异单元,其余区域采用常规有限单元,避免了局部网格加密。本文使用的辛解析奇异单元不需要过渡单元和局部网格加密,且能够通过奇异单元内部的参数关系直接给出切口尖端的应力强度因子,不需要复杂的后处理过程。数值结果表明,本文方法具有良好的精度和稳定性,可以准确地计算动态应力强度因子。     

Crack nucleation from a single notch caused by stress-dependent surface reactions

The surface of a solid under stress is unstable if there is mass exchange and transportation along the surface. A notch on the surface can be a preferred site for crack nucleation. This paper studies the evolution of a surface notch under stress dependant reaction. The surface is represented by a family of curves with many degrees of freedom, and the elastic field is solved using complex conformal mapping method. In the numerical simulations, the notch either deepens immediately and forms singular tip; or becomes perfectly flat after long time; or becomes blunter first, evolving slowly for a long time, then nucleates new surface instability and finally forms sharp tip; or becomes near flat first, then forms a bump on the surface. When the activation strain is small or the reaction is near equilibrium, the notch stability is mainly determined by the competition of strain energy and surface energy, and the sign of stress has little effect. When the surface reaction is away from equilibrium and activation strain is not small, stresses with different signs give totally different surface stability behaviors, depending on whether the solid is losing mass to or gaining mass from the environment.     

Stress concentration factors of periodic notches determined from the strain energy density

Stress concentration factors (SCFs) of a number of flat plates and round bars with periodic U- and V-notches are evaluated. Tension, bending and torsion loadings are considered in the investigation. The main objective of the investigation is to take advantage of the local strain energy density (SED) averaged on a control volume surrounding the tip of the middle notch and to estimate the SCF of each component by using a relatively coarse mesh. The unique advantage of SED method is the most prominent application of such a technique in the current study. Systematic FE simulations by considering a wide range of notch acuity and relative frequency of periodic U- and V-notch components are performed. More than two hundred and fifty models have been examined. The results of this study are compared with those provided by other researchers in the past and recent literature. Two new expressions of the notch depth reduction factor for the case of normal stresses (tension and bending) and torsion are also proposed to match the results from SED approach.     

Strain energy approach to compute stress intensity factors for isotropic homogeneous and bi-material V-notches

A strain energy approach (SEA) is developed to compute the general stress intensity factors (SIFs) for isotropic homogeneous and bi-material plates containing cracks and notches subject to mode I, II and III loading conditions. The approach is based on the strain energy of a control volume around the notch tip, which may be computed by using commercial finite element packages. The formulae are simple and easy to implement. Various numerical examples are presented and compared to corresponding published results or results that are computed using different numerical methods to demonstrate the accuracy of the SEA. Many of those results are new, especially for the cases of bi-material notches where the problem is quite complicated.     

孔洞对平面裂纹应力强度因子影响分析的广义参数Williams单元

本文采用圆形奇异区广义参数Williams单元（W单元）建立了中心裂纹与圆孔共存的平面应力模型,奇异区外围利用ABAQUS有限元软件自动网格离散技术与FORTRAN95编程前处理相结合,克服了自主编程中网格离散的局限性．算例分析了圆孔位置和几何参数对I-II混合型裂纹尖端应力强度因子(SIFs)的影响,并与扩展有限元法(XFEM)计算结果进行比较．结果表明：靠近圆孔一侧的裂尖SIFs大于远离圆孔一侧的裂尖SIFs;控制圆孔左边缘到裂纹中心的距离,则两侧裂尖SIFs随圆孔半径的增大而增大;圆孔中心与裂纹中心水平距离越远,圆孔对裂纹扩展的影响越小．同时,基于圆形奇异区的W单元直接计算得到的裂尖SIFs与扩展有限元法得到的解吻合较好,证明了W单元对奇异区离散形状不敏感,且具有高效率和高精度．     

A unified approach to problems of stress concentration near V-shaped notches with sharp and rounded tip

The paper proposes a unified approach to problems of stress concentration near notches with sharp and rounded tip based on the method of singular integral equations. A solution for an elastic region having a V-shaped notch with rounded tip of large curvature is first found. Then, the stress intensity factor at the tip of a sharp-tipped notch is calculated by passing to the limit. Numerical results are obtained for a slit and a square hole in an elastic plane and an edge notch in a half-plane __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 2, pp. 70–87, February 2007. For the centenary of the birth of G. N. Savin.     

Calculation of Notch H-Integrals Using Image Correlation Experiments

This study evaluated notch H-integrals as well as stress intensity factors (SIFs) using image-correlation experiments for anisotropic materials. First, complex displacement and stress functions are deduced into an H-integral equation. Displacements and stresses from image-correlation experiments are then substituted into the H-integral equation to evaluate the notch SIFs. Experimental results compared with finite element analyses show that the SIFs evaluated using the current method are acceptably accurate.