带裂纹厚壁圆筒应力强度因子的几种计算方法.docVIP

第26卷第4期南　京　理　工　大　学　学　报 Vol.26No.4带裂纹厚壁圆筒应力强度因子的 几种计算方法Ξ 陈爱军　徐　诚① (,ΞΞ 摘要　。该文在考虑裂纹尖端应力应变,,并且利用边界配置法的结果比较这2种方法的精度。同,关键词　应力强度因子,有限元分析,厚壁壳体 分类号　O346.1,TJ301 承受内压的厚壁筒在工程中广泛使用,如:工业管道、枪管炮管等等。厚壁筒在工作时,由于高压的循环作用,内壁常产生裂纹,继续工作,会由于裂纹的扩展而导致破坏,为了工作安全,一般采用断裂力学理论来分析,因此计算厚壁筒在内压p作用下的应力强度因子KIp具有很大的实际价值。确定应力强度因子的方法有很多[1～3],如:数值法、解析法、实验法等。解析法能解决的问题比较少,数值法则在工程中被广泛采用。数值法求解应力强度因子的方法有差分法、边界元、有限元法、边界配位法、权函数法等等。有限元法由于不受所研究的裂纹体几何形状和所受载荷复杂性的限制,在断裂力学中有着非常广泛的应用。本文通过有限单元法计算了裂纹体裂纹尖端附近的应力应变场,用位移法和应力法分别计算了静载下裂纹尖端的应力强度因子,并作了一些相应的分析和研究。 1　位移法求解应力强度因子的原理 如图1所示,厚壁筒的内壁考虑为具有深度为a的穿透型边裂纹。对于I型裂纹,取定Westergaard应力函数,经推导计算得裂纹尖端区域应力场的位移分量表达式为[4] U=4G (2k-1)cos(2k+1)sin-cos22θ-sin22V=4G(1) Ξ 　　Ξ ΞΞ收稿日期:2001-03-30　国防科技预研行业基金资助项目陈爱军　男　30岁　博士生 总第125期　　陈爱军　徐　诚　胡小秋　带裂纹厚壁圆筒应力强度因子的几种计算方法431式中,k=1+μ平面应力,本文以后的公式都是按平面应变公式而得到的。 3-4μ平面应变 图1　带裂纹厚壁筒模型图Fig.1　 Modelofthickwall2cracktip cylinder图2),得到裂纹面附近各点的位移公式=0V= 从而得　KI=?4KI(1-μ2)πE(2)(3)2r4(1-μ) 式中,V可以由有限元方法计算得到。 如图3,裂纹尖端处r=0。考虑到有限元计算位移 值的精度,计算裂纹尖端处的应力强度因子KI时,可从 裂纹面上取3个不同点rA、rB和rC,代入(3)式,可得 KIA=2rA4(1-μ)(4)KIB= KIC=2rB4(1-μ)2rC4(1-μ) 在r很小的情况下,利用拉格朗日插值方法,可求图3　求解裂纹尖端处应力强度因子示意图Fig.3　Solutionmodelofstress intensityfactorofcracktip +?K(rC-rA)(rC-rB)(5)得精度良好的裂纹尖端(r=0)处的应力强度因子?K?K　　KIr=0=(+)()()(rA-rBrA-rCrB- rArB-rC) 2　应力法求解应力强度因子原理 应力法求解应力强度因子与位移法很相似。首先求出裂纹尖端附近的应力场,再按应力与应力强度因子的关系求出应力强度因子。I型裂纹尖端附近的应力场公式为 σx= σy=πrcoscos(1-sin)?sin222(1+sin)?sin222πr裂纹平面的延长线面上,θ=0°,则可得(6) 432南　京　理　工　大　学　学　报　　　　　　　　　第26卷第4期 πr　　KI=σy?(7) 用有限元方法求得不同点的应力σy,按(7)式求出相应点的KI值,再利用前面位移法类似的插值方法可求得裂纹尖端的应力强度因子。 3　单元网格的划分和等参奇异元的应用、计算结果及其分析 厚壁筒处于平面应变状态,由于对称性,取厚壁筒的一半作为研究对象。力学可知[5～6],在裂纹尖端应变具有r-的奇异性,则不 能很好地反映这种奇异性。Barsoum,均4 存在的奇异性。 ,而在裂纹尖端附近采用退化三r2 ,。网格的划分见图4所示。图4　有限元模型图图5　加载曲线 Fig.4　SolutionmodelofstressintensityfactorsofcracktipFig.5　Theloadingcurve本文计算的厚壁圆筒的内外半径比R1/R2=0.444,材料为CrNiMoV钢,材料性能为E=200214MPa,σs=128kg/mm,μ=0.27,裂纹为单边静止裂纹,采用各向同性材料线2 弹性模型。本文研究计算静态应力强度因子所采用的加载曲线如图5。加载时间总长t=3s;计算时间步数为3,步长为1s。 利用有限元的计算结果,采用应力法和位移法计算了厚壁筒尺寸不变的条件下(R1=50mm, R2=113mm),裂纹深度a变化下的静态应力强度因子KI,见图6。以及a/d= 0.1时