资源描述：
《带裂纹厚壁圆筒应力强度因子的几种计算方法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、第26卷第4期南京理工大学学报Vol.26No.42002年8月JournalofNanjingUniversityofScienceandTechnologyAug.2002带裂纹厚壁圆筒应力强度因子的X几种计算方法XX①①陈爱军徐诚胡小秋(南京理工大学理学院,①机械工程学院,南京210094)摘要确定应力强度因子是断裂力学的重要内容。该文在考虑裂纹尖端应力应变奇异性的前提下,通过有限元的位移法和应力法分别计算了承受高压厚壁筒裂纹尖端处的应力强度因子,并且利用边界配置法的结果比较这2种方法的精度。同时,还研究应力强度因子随裂纹深
2、度和厚壁筒尺寸的变化规律关键词应力强度因子,有限元分析,厚壁壳体分类号O346.1,TJ301承受内压的厚壁筒在工程中广泛使用,如:工业管道、枪管炮管等等。厚壁筒在工作时,由于高压的循环作用,内壁常产生裂纹,继续工作,会由于裂纹的扩展而导致破坏,为了工作安全,一般采用断裂力学理论来分析,因此计算厚壁筒在内压p作用下的应力强度因子KIp[1～3]具有很大的实际价值。确定应力强度因子的方法有很多,如:数值法、解析法、实验法等。解析法能解决的问题比较少,数值法则在工程中被广泛采用。数值法求解应力强度因子的方法有差分法、边界元、有限元法、
3、边界配位法、权函数法等等。有限元法由于不受所研究的裂纹体几何形状和所受载荷复杂性的限制,在断裂力学中有着非常广泛的应用。本文通过有限单元法计算了裂纹体裂纹尖端附近的应力应变场,用位移法和应力法分别计算了静载下裂纹尖端的应力强度因子,并作了一些相应的分析和研究。1位移法求解应力强度因子的原理如图1所示,厚壁筒的内壁考虑为具有深度为a的穿透型边裂纹。对于I型裂纹,取定[4]Westergaard应力函数,经推导计算得裂纹尖端区域应力场的位移分量表达式为KIrθ3U=(2k-1)cos-cosθ4G2π22(1)KIrθ3V=(2k+1
4、)sin-sinθ4G2π22X收稿日期:2001-03-30X国防科技预研行业基金资助项目XX陈爱军男30岁博士生总第125期陈爱军徐诚胡小秋带裂纹厚壁圆筒应力强度因子的几种计算方法4313-μ平面应力式中,k=1+μ,本文以后的公式都是按平面应变公式而得到的。3-4μ平面应变图1带裂纹厚壁筒模型图图2裂纹尖端坐标系Fig.1ModelofthickwallFig.2Coordinatesofcracktipcylinderwithcrack图2为裂纹尖端坐标系。将θ=π代入(1)式,得到裂纹面附近各点的位移公式r21U=0V=
5、·4KI(1-μ)(2)2πE2πE·V从而得KI=2(3)r4(1-μ)式中,V可以由有限元方法计算得到。如图3,裂纹尖端处r=0。考虑到有限元计算位移值的精度,计算裂纹尖端处的应力强度因子KI时,可从裂纹面上取3个不同点rA、rB和rC,代入(3)式,可得2πE·VAKIA=·2rA4(1-μ)2πE·VBKIB=·2(4)rB4(1-μ)图3求解裂纹尖端处应力2πE·VC强度因子示意图KIC=·2rC4(1-μ)Fig.3Solutionmodelofstress在r很小的情况下,利用拉格朗日插值方法,可求intensity
6、factorofcracktip得精度良好的裂纹尖端(r=0)处的应力强度因子rBrC·KIArArC·KIBrArB·KICKIr=0=(r)(r)+(r)(r)+(r)(r)(5)A-rBA-rCB-rAB-rCC-rAC-rB2应力法求解应力强度因子原理应力法求解应力强度因子与位移法很相似。首先求出裂纹尖端附近的应力场,再按应力与应力强度因子的关系求出应力强度因子。I型裂纹尖端附近的应力场公式为KIθθ3σx=cos(1-sin·sinθ)2πr222(6)KIθθ3σy=cos(1+sin·sinθ)2πr222裂纹平面的
7、延长线面上,θ=0°,则可得432南京理工大学学报第26卷第4期KI=σy·2πr(7)用有限元方法求得不同点的应力σy,按(7)式求出相应点的KI值,再利用前面位移法类似的插值方法可求得裂纹尖端的应力强度因子。3单元网格的划分和等参奇异元的应用、计算结果及其分析厚壁筒处于平面应变状态,由于对称性,取厚壁筒的一半作为研究对象。由线弹性断裂1[5～6]-力学可知,在裂纹尖端应变具有r2的奇异性,如果利用常规的单元来划分网格,则不1能很好地反映这种奇异性。Barsoum证明了退化为三角形的中点元在单元内任意一点,均41存在的奇异性。本
8、文划分网络单元采用八节点等差单元,而在裂纹尖端附近采用退化三r1角形的中点奇异等参元,进行有限元计算。网格的划分见图4所示。4图4有限元模型图图5加载曲线Fig.4Solutionmodelofstressintensityfactorsofcra
