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

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

内部压力作用下无限大板中三角孔-边裂纹应力强度因子研究

利用已有文献中的杂交位移不连续边界元法,重点研究了内部压力作用下无限大板中三角孔-边裂纹问题;通过改变孔的几何参数,分析了孔的几何参数对应力强度因子的影响。结果表明:孔对源于其裂纹的应力强度因子具有屏蔽影响和放大影响;当尺寸参数ad ≥adc(adc为某一定值)时,三角孔对源于其裂纹的应力强度因子具有屏蔽影响,并且三角孔尺寸越接近裂纹尺寸,这种屏蔽影响越强烈;当参数 ad≤adc时,三角孔对源于其裂纹的应力强度因子具有放大影响,并且在参数 ad=adm(adm为某一定值)处,这种放大影响达到最大。本文所得结果在工程上具有重要意义。     

改进型扩展比例边界有限元法

结合了扩展有限元法(extended finite elementmethods,XFEM)和比例边界有限元法(scaled boundary finite elementmethods,SBFEM)的主要优点,提出了一种改进型扩展比例边界有限元法(improvedextended scaled boundary finite elementmethods,$i$XSBFEM),为断裂问题模拟提供了一条新的途径.类似XFEM,采用两个正交的水平集函数表征材料内部裂纹面,并基于水平集函数判断单元切割类型;将被裂纹切割的单元作为SBFE的子域处理,采用SBFEM求解单元刚度矩阵,从而避免了XFEM中求解不连续单元刚度矩阵需要进一步进行单元子划分的缺陷;同时,借助XFEM的主要思想,将裂纹与单元边界交点的真实位移作为单元结点的附加自由度考虑,赋予了单元结点附加自由度明确的物理意义,可以直接根据位移求解结果得出裂纹与单元边界交点的位移;对于含有裂尖的单元,选取围绕裂尖单元一圈的若干层单元作为超级单元,并将此超级单元作为SBFE的一个子域求解刚度矩阵,超级单元内部的结点位移可通过SBFE的位移模式求解得到,应力强度因子可基于裂尖处的奇异位移(应力)直接获得,无需借助其他的数值方法.最后,通过若干数值算例验证了建议的$i$XSBFEM的有效性,相比于常规XFEM,$i$XSBFEM的基于位移范数的相对误差收敛性较好;采用$i$XSBFEM通过应力法和位移法直接计算得到的裂尖应力强度因子均与解析解吻合\较好.      

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

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

基于修正权函数的多裂纹无网格模拟

本文采用了一种基于不连续场修正权函数的无网格方法来处理二维平面多裂纹问题。相较于传统的无网格断裂不连续场和奇异场模拟方法,修正权函数法算法简便易实现。采用修正权函数处理多裂纹时,只需要对每一段裂纹周围节点的权函数进行修正,就能同时模拟多裂纹不连续位移场和多裂尖奇异场。本文采用基于不连续场修正权函数的无单元Galerkin方法（EFGM）,对Y型裂纹板、十字型裂纹板和孔边双裂纹板进行了分析。数值结果表明,在不引入扩展基函数情况下,通过修正权函数法能够得到精度较高的应力强度因子解,能较好地拟合多裂纹的裂尖奇异场。     

I型裂纹尖端约束应力区模型及其解析解

考虑了I型裂纹尖端损伤区域内三种不同的约束应力分布形式,即右三角分布形式(情况A)、均匀分布形式(情况B)、左三角分布形式(情况C),并采用复变函数方法求得了应力强度因子与裂纹张开位移的解析解;在此基础上,通过数值计算得到了应力强度因子和裂纹张开位移随约束应力区长度、约束应力大小以及分布形式的变化规律。研究结果表明:随裂尖材料损伤程度的增加,裂尖损伤区内约束应力减小,应力强度因子和裂纹张开位移增大;约束应力的分布形式对应力强度因子和裂纹张开位移有显著影响;相对于其他区域,约束应力对裂纹尖端区域裂纹张开位移的影响较大。然而,对于裂尖损伤区域的形成与作用荷载、材料性质、构件几何尺寸之间的关系,还需要进行更为深入的研究。     

任意多椭圆孔多裂纹无限大各向异性板应力强度因子求解的一种新方法

采用各向异性体平面弹性理论中的复势方法,应用保角变换技术,以F aber级数为工具,导出含任意多椭圆孔和裂纹群无限大各向异性板在远场载荷作用下其应力场和位移场的级数解,并在此基础上利用断裂力学方法确定裂纹尖端的应力强度因子,通过算例讨论了材料参数及裂纹、孔的尺寸等对应力强度因子的影响规律,得出了一些有益的结论。数值结果表明本文方法具有计算精度高、收敛速度快、方便快捷等优点,有利于全面系统地研究各参数对结构断裂性能的影响。     

任意多孔多裂纹有限大板的应力强度因子分析

采用各向异性体平面弹性理论中的复势方法,以Faber级数为工具,应用保角映射技术和最小二乘边界配点法,导出内边界条件精确满足,外边界条件近似满足的含多椭圆孔及裂纹群有限大板在任意载荷作用下的应力场、位移场的级数解,建立了任意多椭圆孔及裂纹群有限大板应力强度因子的有效分析方法,讨论了各参数对裂尖应力强度因子及孔边应力集中的影响．数值结果表明,该方法具有计算精度高、收敛速度快、方便快捷等优点,有利于全面系统地研究各参数对结构断裂性能的影响．     

夹杂对两相材料界面裂纹的影响

根据界面上应力和位移的连续条件,得到了单向拉伸状态下,含有椭圆夹杂的无限大双材料组合板的复势解。进一步通过求解Hilbert问题,得到了含有夹杂和半无限界面裂纹的无限大板的应力场,并由此给出了裂尖的应力强度因子K。计算了夹杂的形状、夹杂的位置、夹杂的材料选取以及上、下半平面材料与夹杂材料的不同组合对裂尖应力强度的影响。计算结果表明夹杂到裂尖的距离和夹杂材料的性质对K影响较大,对于不同材料组合,该影响有较大差异。夹杂距裂尖较近时,会对K产生明显屏蔽作用,随着夹杂远离裂尖,对K的影响也逐渐减小。另外,软夹杂对K有屏蔽作用,硬夹杂对K有反屏蔽作用,而夹杂形状对K几乎没有影响。     

任意多椭圆孔多裂纹无限大板的应力强度因子与M积分

基于弹性力学中平面问题的复势方法,应用保角映射技术,以Faber级数为工具,导出含任意多椭圆孔及裂纹群无限大板在任意载荷作用下其应力场和位移场的级数解,并在此基础上计算了任意多裂纹板的应力强度因子和M积分,数值结果表明,该方法具有计算精度高、收敛速度快、方便快捷等解析法特有的优点。通过算例分析了不同裂纹倾角时M积分值随载荷方向的变化关系,并讨论了裂纹长度、裂纹间距及裂纹倾角等参数对应力强度因子的影响规律,获得了一些重要结论.     

An effective boundary element method for analysis of crack problems in a plane elastic plate

A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.     

A special crack tip displacement discontinuity element

Based on the analytical solution to the problem of a constant discontinuity in displacement over a finite line segment in the x, y plane of an infinite elastic solid and the note of the crack tip element by Crouch, in the present paper, the special crack tip displacement discontinuity element is developed. Further the analytical formulas for the stress intensity factors of crack problems in general plane elasticity are given. In the boundary element implementation the special crack tip displacement discontinuity element is placed locally at each crack tip on top of the non-singular constant displacement discontinuity elements that cover the entire crack surface. Numerical results show that the displacement discontinuity modeling technique of a crack presented in this paper is very effective.     

A numerical analysis of cracks emanating from an elliptical hole in a 2-D elasticity plate

This paper is concerned with the stress intensity factors (SIFs) of cracks emanating from an elliptical hole in an infinite or a finite plate under biaxial loads by using a boundary element method, which consists of the non-singular displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements due to the author. In the boundary element implementation the left or the right crack-tip element is placed locally at the corresponding left or right crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. A few numerical examples are included to show that the present approach is very efficient and accurate for the calculating the SIFs of crack problems in an infinite or a finite plate. The present numerical results of cracks emanating from an elliptical hole under biaxial loads can reveal the effect of the elliptical aspect ratio and the transverse load on the SIFs.     

FATIGUE GROWTH MODELING OF MIXED-MODE CRACK IN PLANE ELASTIC MEDIA

This paper presents an extension of a displacement discontinuity method with cracktip elements （a boundary element method） proposed by the author for fatigue crack growth analysis in plane elastic media under mixed-mode conditions. The boundary element method consists of the non-singular displacement discontinuity elements presented by Crouch and Starfield and the crack-tip displacement discontinuity elements due to the author. In the boundary element implementation the left or right crack-tip element is placed locally at the corresponding left or right crack tip on top of the non-singular displacement discontinuity elements that cover the entire crack surface and the other boundaries. Crack growth is simulated with an incremental crack extension analysis based on the maximum circumferential stress criterion. In the numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not required because of an intrinsic feature of the numerical approach. Crack growth is modeled by adding new boundary elements on the incremental crack extension to the previous crack boundaries. At the same time, the element characteristics of some related elements are adjusted according to the manner in which the boundary element method is implemented. As an example, the fatigue growth process of cracks emanating from a circular hole in a plane elastic plate is simulated using the numerical simulation approach.     

Stress intensity factors for asymmetric branched cracks in plane extension by using crack-tip displacement discontinuity elements

Stress intensity factors are important in the analysis of cracked materials. They are directly related to the fracture propagation and fatigue crack growth criteria. Based on the analytical solution (Crouch, S.L., 1976. Solution of plane elasticity problems by displacement discontinuity method, Int. J. Numer. Methods Eng. 10, pp. 301–343; Crouch, S.L., Starfield, A.M., 1983. Boundary Element Method in Solid Mechanics, with Application in Rock Mechanics and Geological Mechanics, London, Geore Allon and Unwin, Bonton, Sydney) to the problem of a constant discontinuity in displacement over a finite line segment in the x, y plane of an infinite elastic solid, recently, the crack-tip displacement discontinuity element which can be classified as the left and right crack-tip displacement discontinuity elements are developed by the author Yan, X., (in press. A special crack-tip displacement discontinuity element, Mechanics Research Communications) to model the crack-tip fields to more accurately compute the stress intensity factors of cracks in general plane elasticity. In the boundary element implementation the left or the right crack-tip displacement discontinuity element is placed locally at the corresponding left or right crack tip on top of the ordinary non-singular displacement discontinuity elements that cover the entire crack surface and the other boundaries. To prove further the efficiency of the suggested approach and provide more results of the stress intensity factors, in this study, analysis of an asymmetric branched crack bifurcated from a main crack in plane extension is carried out.     

A numerical analysis of perpendicular cracks under general in-plane loading with a hybrid displacement discontinuity method

In this paper, a numerical analysis of perpendicular cracks under general in-plane loading is performed by using a hybrid displacement discontinuity method which consists of the non-singular displacement discontinuity element presented by Crouch and Starfied and the crack tip displacement discontinuity elements by the author. In the boundary element implementation the left or the right crack tip displacement discontinuity element is placed locally at corresponding left or right crack tip on top of the ordinary non-singular displacement discontinuity elements that cover the entire crack surface and the other boundary. The present numerical results show that the numerical approach is simple, yet very accurate for calculating numerically stress intensity factors for perpendicular cracks under general in-plane loading.     

A boundary element modeling of fatigue crack growth in a plane elastic plate

This paper presents an extension of a boundary element method to fatigue growth analysis of mixed-mode cracked plane elastic bodies. The method consists of the non-singular displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity element due to the author. In the boundary element implementation the left or the right crack-tip element is placed locally at the corresponding left or right crack tip on top of non-singular displacement discontinuity elements that cover the entire crack surface and the other boundaries. Crack growth is simulated with an incremental crack extension analysis based on the modified maximum strain energy density criterion. In numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not required because of an intrinsic feature of the boundary element method. Crack growth is simulated by adding new boundary elements on the incremental crack extension to the previous crack boundaries. At the same time, the element characters of some related elements are adjusted according to the manner in which the boundary element method is implemented. Some numerical results of fatigue growth in a plane elastic plate with a center-inclined crack under uniaxial cyclic loading are given.     

The analysis of crack problems with non-local elasticity

IntroductionTheclassicalconhnuummechanicshasbeenusedtosolvemanyproblemsinmacrofracturemechanics,butencountersdifficulheswhentheeffectofITilcrocharacteristicdimensionshouldbetakenintoaccount.Thestressfieldverynearthecracktipisstillnotclear.Somephenomenaofshortcrackscannotbeexplained["']andsomemechanismoffracturehasnotbeensolvedyet.Thenon-localelashcitytheoryseemsattractivetotheseproblems.Thetheoryofnon-localelasticity,establishedanddevelopedbyEringenetal[3),connectstheclassicalcontinuummechan…     

压电材料反平面界面裂纹的杂交元分析

基于奇异性电弹场数值特征解开发了一种新型反平面界面裂纹尖端单元。将新型单元与四节点压电P-S单元组装,求解从绝缘到导通的任意电边界条件下,压电结构反平面界面裂纹尖端电弹场的数值解。考察了层厚、载荷类型和裂纹面间电边界条件等对反平面界面裂纹尖端断裂参数的影响。算例证明新型单元能使P-S单元数显著降低,计算结果更为精确。