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

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

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.     

角度非均匀材料平面V形切口应力奇性分析

提出了一种确定角度非均匀材料平面V形切口尖端应力奇性指数的有效方法。首先,在弹性力学基本方程中引入V形切口尖端位移场的级数渐近展开,建立以位移为特征函数的变系数和非线性微分方程组。然后,采用微分求积法(DQM)求解微分方程组,可得到多阶应力奇性指数及其相对应的特征函数,该法具有公式简单、编程方便、计算量少和精度高等优点,可处理任意开口角度和任意材料组合的V形切口。典型算例验证了微分求积法的有效性和精确性。     

A new boundary element approach of modeling singular stress fields of plane V-notch problems

In this paper, a new boundary element (BE) approach is proposed to determine the singular stress field in plane V-notch structures. The method is based on an asymptotic expansion of the stresses in a small region around a notch tip and application of the conventional BE in the remaining region of the structure. The evaluation of stress singularities at a notch tip is transformed into an eigenvalue problem of ordinary differential equations that is solved by the interpolating matrix method in order to obtain singularity orders (degrees) and associated eigen-functions of the V-notch. The combination of the eigen-analysis for the small region and the conventional BE analysis for the remaining part of the structure results in both the singular stress field near the notch tip and the notch stress intensity factors (SIFs).Examples are given for V-notch plates made of isotropic materials. Comparisons and parametric studies on stresses and notch SIFs are carried out for various V-notch plates. The studies show that the new approach is accurate and effective in simulating singular stress fields in V-notch/crack structures.     

橡胶楔体与刚性缺口接触的大变形分析

利用一种新的橡胶材料应变能函数，对橡胶楔体与刚性缺口接触大变形问题进行了分析。得到了接触尖点附近变形的奇异性特征，给出了奇异性指数与材料常数、橡胶楔体角度、刚性缺口角度之间的关系式。同时编制了大变形有限元程序，计算得到了与理论解一致的结论。     

Path-independent integral for the sharp V-notch in longitudinal shear problem

By applying Noether’s theorem to the elastic energy density in longitudinal shear problem, it is shown that its symmetry-transformations of material space can be expressed by the real and imaginary parts of an analytic function. This kind of the symmetry-transformations leads to the existence of a conservation law in material space, which does not belong to trivial conservation laws and whose divergence-free expression gives a path-independent integral. It is found that by adjusting the analytic function, a finite value can be obtained from this path-independent integral calculated around the material point with any order singularity. For a sharp V-notch placed on the edge of homogenous materials and/or the interface of bi-materials, application shows that the finite value obtained from this path-independent integral is directly related to the notch stress intensity factor (NSIF) and does not depend on the location of integral endpoints chosen respectively along two traction-free surfaces of which form a notch opening angle. Usability is presented in an example to estimate the NSIF of a bi-material plate.     

Higher order fields for damaged nonlinear antiplane shear notch, crack and inclusion problems

Damaged nonlinear antiplane shear problems with a variety of singularities are studied analytically. A deformation plasticity theory coupled with damage is employed in analysis. The effect of microscopic damage is considered in terms of continuum damage mechanics approach. An exact solution for the general damaged nonlinear singular antiplane shear problem is derived in the stress plane by means of a hodograph transformation, then corresponding higher order asymptotic solutions are obtained by reversing the stress plane solution to the physical plane. As example, traction free sharp notch and crack, rigid sharp wedge and flat inclusion, and mixed boundary sharp notch problems are investigated, respectively. Consequently, higher order fields are obtained, in which analytical expressions of the dominant and second order singularity exponents and angular distribution functions of the near tip fields are derived. Effects of the damage and hardening exponents of materials and the geometric angle of notch/wedge on the near tip quantities are discussed in detail. It is found that damage leads to a weaker dominant singularity of stress, but to little stronger singularities of the dominant and second order terms of strain compared to that for undamaged material. It is also seen that damage has important effect on the angular distribution functions of the near tip stress and strain fields. As special cases, higher order analytical solutions of the crack and rigid flat inclusion tip fields are obtained, respectively, by reducing the notch/wedge tip solutions. Effects of damage and hardening exponents on the dominant and second order terms in the solutions of the crack and inclusion tip fields are discussed.     

单材料V型缺口尖端振荡性奇异应力场产生的条件

单材料V型缺口附近应力场存在奇异性,Williams在1952年针对不同边界条件下所产生的奇异性进行了讨论,结论表明,边界条件和材料的泊松比对奇异 均有影响,本文对Williams所提出的第三种边界条件（一边自由,一边固支）研究后发现,缺口尖端附近应力不仅存在幂次奇异,而且还会出现振荡性,振荡指数大小依赖于缺口角度和泊松比。     

Analytical representation of the non-square-root singular stress field at a finite angle sharp notch

The stress field near the tip of a finite angle sharp notch is singular. However, unlike a crack, the order of the singularity at the notch tip is less than one-half. Under tensile loading, such a singularity is characterized by a generalized stress intensity factor which is analogous to the mode I stress intensity factor used in fracture mechanics, but which has order less than one-half. By using a cohesive zone model for a notional crack emanating from the notch tip, we relate the critical value of the generalized stress intensity factor to the fracture toughness. The results show that this relation depends not only on the notch angle, but also on the maximum stress of the cohesive zone model. As expected the dependence on that maximum stress vanishes as the notch angle approaches zero. The results of this analysis compare very well with a numerical (finite element) analysis in the literature. For mixed-mode loading the limits of applicability of using a mode I failure criterion are explored.     

插值矩阵法分析双材料平面V形切口奇异阶

对二维V形切口问题提出奇异阶分析的一个新方法.首先,以V形切口尖端附近位移场沿其径向渐近展开为基础,将其线弹性理论控制方程转换成切口尖端附近关于周向变量的常微分方程组特征值问题,然后将数值求解两点边值问题的插值矩阵法进一步拓展为求解一般常微分方程组特征值问题,插值矩阵法是在离散节点上采用微分方程中待求函数的最高阶导数作为基本未知量.由此,V形切口的应力奇性阶问题通过插值矩阵法获得,同时相应的切口附近位移场和应力场特征向量一并求出.     

垂直于界面V型缺口尖端的应力奇异性及其应用

基于双材料垂直于界面V型缺口理论,给出了单一材料和双材料裂纹问题、V型缺口问题应力强度因子的统一定义,得到了应力外推法计算双材料K_I的公式,数值算例验证了本文方法的有效性.以双材料单向拉伸和三点弯曲模型为对象,深入研究了双材料中弹性模量、泊松比、缺口深度、缺口张角对缺口尖端奇异应力场的影响,获得了一定范围内各种参数变化对缺口尖端奇异应力场的影响规律,为异体材料形成的V型缺口在应力断料中的应用提供了必要的参考依据.     

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.     

粘弹性双材料界面裂缝的缝端位移场及动态断裂分析

本文研究的是经常在实际工程中遇到的粘弹性双材料界面裂缝的动断裂问题.由于粘弹性自身的复杂性,使得粘弹性双材料界面裂缝缝端应力的奇异性较弹性呈现出更为复杂的形式,从而使动断裂问题的分析变得更为困难.根据此情况,本文采用复阻尼理论反映粘弹性体的运动规律,用复势理论和平面问题复变函数解答的科洛索夫公式推导了粘弹性双材料界面裂缝缝端位移场及动态应力强度因子的求解公式,利用特解边界元进行了粘弹性双域耦合动力响应计算,按求得的公式用位移外推法计算了单边裂纹板在动荷载作用下的动态应力强度因子.分析了粘性,弹模比和缝长对动态应力强度因子的影响,得出了一些有益的结论.     

Higher order asymptotic fields at the tip of a sharp V-notch in a power-hardening material

The higher order asymptotic fields at the tip of a sharp V-notch in a power-hardening material for plane strain problem of Mode I are derived. The order hierarchy in powers ofr for various hardening exponentsn and notch angles β is obtained. The angular distributions of stress for several cases are plotted. The self-similarity behavior between the higher order terms is noticed. It is found that the terms with higher order can be neglected for the V-notch angle β>45°. Project supported by the National Natural Science Foundation of China (Nos. 10132010 and 10072033).     

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 形切口/裂纹结构完整应力场的准确性和有效性.      

Stress singularities at crack corners

In this paper the stress and displacement fields near an embedded crack corner in a linear elastic medium are analytically computed. The conical-spherical coordinate system is introduced to solve this problem. It is observed that the strength of the stress singularity depends on the angle of the crack corner. The singularity becomes weaker, varying from r ^-1 to r ⁰, as the angle of the crack corner varies from 360° to 0°. Both symmetric and skew-symmetric loadings give the same variation of the behavior of the stress singularity. It is also found that the order of the singularity is independent of the Poisson's ratio, unlike the corner cracks at a free surface where Poisson's ratio affects the results.     

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

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

反平面V形切口塑性应力奇异性分析

论文提出了用插值矩阵法计算幂硬化塑性材料反平面V形切口和裂纹尖端区域的应力奇异性.首先在切口和裂纹尖端区域采用自尖端径向度量的渐近位移场假设,将其代入塑性全量理论的基本微分方程后,推导出包含应力奇异性特征指数和特征角函数的非线性常微分方程特征值问题.然后采用插值矩阵法迭代求解导出的控制方程,得到一般的塑性材料反平面V形切口和裂纹的前若干阶应力奇异阶和相应的特征角函数,该法的重要优点是以上求解的特征角函数和它们各阶导函数具有同阶精度,并且一次性地求出前若干阶特征对.同时,插值矩阵法计算量小,易于和其他方法联合使用,这些优点在后续求解尖端区域完全应力场非常优越.论文方法的计算结果与现有结果对照,发现吻合良好,表明了论文方法的有效性.     

压电切口张开角和深度对其尖端力电损伤场的影响

基于三维各向异性压电损伤本构理论,导出了广义平面应力问题的损伤本构方程,并据此分析了压电薄板板边V形切口尖端附近的力电损伤,研究了切口张开角和深度对切口尖端力电损伤的影响规律．结果发现：和张开角对切口尖端损伤的影响相比,深度的影响更为明显;在张开角对切口尖端力损伤的影响规律方面,压电材料与一般弹塑性材料存在明显差异,原因在于压电切口尖端力电载荷比会随着深度的改变发生很大变化;不同深度下张开角与切口尖端力、电损伤关系曲线随着张开角的增大由发散逐渐会聚,不同张开角下深度与切口尖端力、电损伤关系曲线随着切口加深由会聚逐渐发散,并且电损伤曲线表现得更为明显．