厚壁筒內(nèi)裂紋深度與動態(tài)應(yīng)力強度因子的關(guān)系

1、第20卷第6期南京理工大學(xué)學(xué)報Vol.20No.6厚壁筒內(nèi)裂紋深度與動態(tài)應(yīng)力強度因子的關(guān)系楊宏偉馬大為王更新王志群(南京理工大學(xué)機械學(xué)院,摘要根據(jù)有限元法的基本理論,5種裂紋深度(610mm、1410、263210mm)的靜止邊裂紋的厚壁圓筒,設(shè)計了合理的計算動態(tài)應(yīng)力強度因5種裂紋深度下厚壁筒在2種載荷(突加和正弦)作用下的動態(tài)應(yīng)力強度因子,從而得到了一些動態(tài)應(yīng)力強度因子隨裂紋深度和時間變化的規(guī)律。該規(guī)律對厚壁筒的設(shè)計及工程應(yīng)用有一定的參考價值。關(guān)鍵詞裂紋,應(yīng)力強度因子,有限元法,厚壁殼體分類號O346.1工程中廣泛使用厚壁筒構(gòu)件,如石油化工設(shè)備中的高壓容器和管道、槍炮身管等。厚壁筒在工作時,

2、由于內(nèi)膛承受動態(tài)內(nèi)壓,故常在內(nèi)壁出現(xiàn)裂紋,繼續(xù)使用這些厚壁筒時,裂紋就要不斷擴展,從而造成破壞事故,故在進行斷裂力學(xué)的安全分析和疲勞壽命設(shè)計時,需要計算表面裂紋的應(yīng)力強度因子。目前,對厚壁筒在靜態(tài)壓力下和自增強殘余應(yīng)力下的應(yīng)力強度因子的計算研究已有不少成就123。然而,實際應(yīng)用中的厚壁筒,一般承受的是動態(tài)內(nèi)壓,故需研究在動態(tài)內(nèi)壓下厚壁筒的應(yīng)力強度因子。本文根據(jù)有限元理論,計算了厚壁筒在2種動態(tài)內(nèi)壓載荷作用下,不同裂紋深度的動態(tài)應(yīng)力強度因子,并分析了它們之間的變化關(guān)系。1計算方案111基本原理動態(tài)有限元程序是以方程+Cu(1)Mu+Ku=R(t)為基礎(chǔ)的。式中,M、C和K分別表示質(zhì)量、阻尼和剛度

3、矩陣,R(t)是外載荷列陣,u、分別是位移、u和u速度向量和加速度向量。求解每一時刻的動態(tài)方程(1),可求出各時刻的位移u,本文采用位移法求解裂紋尖端的應(yīng)力強度因子。本文于1996年1月27日收到國防科技預(yù)研行業(yè)基金項目楊宏偉男33歲博士生© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. 530南京理工大學(xué)學(xué)報第20卷第6期112裂紋尖端位移場與應(yīng)力強度因子的關(guān)系動態(tài)載荷下靜止裂紋尖端的應(yīng)力場和位移場表達式與靜載下的類似4,可用靜載下的位移與應(yīng)力強度因子的關(guān)系說明動

4、載下的位移與動態(tài)應(yīng)力強度因子的關(guān)系。經(jīng)推導(dǎo),有K1D(t)=)r4(1-(2)式中,V(t)為相應(yīng)時刻裂紋面上各點處垂直于裂紋面的位移,可用有限元法求得,外推到裂紋尖端r=0處,得到對應(yīng)時刻的裂紋尖端的動態(tài)應(yīng)力強度因子;r為所求位移點距裂紋尖端(2)式有r項,在計算時需處理好1r的奇異性問題。的距離;E為楊氏模量;為泊松比。113網(wǎng)格劃分與奇異性處理本文的研究對象為各向同性材料,考慮到對稱性,并假定其處于應(yīng)變狀態(tài)。網(wǎng)格采用八節(jié)點等參元,r形四分之一中點等參元,考慮了合理的優(yōu)化尺寸,1。圖1計算網(wǎng)格Fig.1Calculationgrids圖2載荷Fig.2Pressure2具體應(yīng)用研究對象為內(nèi)

5、半徑R1=6010mm,外半徑R2=12510mm的厚壁筒具有單邊靜止半橢圓表面裂紋,a =0127,E=200GPa,分別研究具有不同裂紋尺寸a的情況。c=015,211對網(wǎng)格合理性的檢驗分別計算了具有5種不同裂紋深度的厚壁筒承受同種內(nèi)壓p的靜態(tài)應(yīng)力強度因子K1S,并同時利用文獻1計算出相應(yīng)的靜態(tài)值K1S1,對兩者作了比較,其結(jié)果見表1,結(jié)果符合較好,這說明網(wǎng)格劃分是合理的。表1靜態(tài)應(yīng)力強度因子KIS的算法比較Table1ComparisonofcalculationalmethodsofK1S6101410201026106261383818621814514521S11 24904904

6、90490180110278149334181397164183155280110328148399174相對誤差(%)119016119015© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. 總第90期楊宏偉馬大為王更新王志群厚壁筒內(nèi)裂紋深度與動態(tài)應(yīng)力強度因子的關(guān)系531212動態(tài)應(yīng)力強度因子計算首先對程序能否用于動態(tài)計算進行了考核,其結(jié)果與文獻5作的比較表明程序設(shè)計合理,然后對以上5種裂紋深度模型分別施加突加載荷和正弦波內(nèi)壓2種動態(tài)載荷,見圖2(a)、(b)。最

7、大動載值為p0=490MPa,正弦波內(nèi)壓為p(t)=p0sin(t)。2種載荷的作用時間相6同,均為0t610ms。通過計算,得到分別在這2種載荷作用下,具有5種裂紋尺寸的厚壁筒的裂紋尖端的動態(tài)應(yīng)力強度因子,其結(jié)果見圖3和圖4。圖3沖擊載荷下的動態(tài)應(yīng)力強度因子Fig.3Dynamicstressintensityfactorsunderimpactpressure圖4正弦波載荷下的動態(tài)應(yīng)力強度因子Fig.4Dynamicstressintensityfactorsundersinepressure為了研究動態(tài)應(yīng)力強度因子的變化規(guī)律,表2分別列出了在2種載荷作用下5種不同裂紋尺寸的一些計算結(jié)果,

8、其中,B=,式中,K1Dmax表示最大動態(tài)應(yīng)力強度因子,t表示B達到K1S最大值的時刻,T表示整個載荷作用的時間。表2應(yīng)力強度因子比較Table2Comparisonofstressintensityfactorsa mmB突加載荷t ms510514515519610T msB正弦載荷t ms512515515518610T ms6101410201026103210117561160111573115301148161061061061061011222111371107911029019646106106106106103結(jié)果分析對于同種動態(tài)加載載荷,隨著裂紋深度a的增加。經(jīng)分析計算結(jié)果

9、,得如下規(guī)律:(1)對應(yīng)于某一時刻的裂紋尖端的K1D(t)也逐漸變大,且它的增長速率也在變大,在加載初期,K1D(t)與a和時間呈近似的線性關(guān)系,而與加載方式關(guān)系不大;(2)裂紋尖端的K1Dmax相對于各自的靜態(tài)值的比值在逐漸地減小,K1Dmax出現(xiàn)的時間在逐漸地增長,這主要是由于載荷作用在裂紋面和內(nèi)壁上,裂紋越深,應(yīng)力波傳播到裂紋尖端所用的時間越長之故;© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. 532南京理工大學(xué)學(xué)報第20卷第6期(3)對于正弦載荷,當(dāng)裂紋

10、深度達到一定值時,其動態(tài)應(yīng)力強度因子的最大值開始小于靜態(tài)值,這說明動態(tài)應(yīng)力強度因子的最大值與載荷的加載方式關(guān)系較緊密;(4)對2種不同的作用載荷,當(dāng)K1Dmax出現(xiàn)后,K1D(t)都會有不同程度的下降,這說明了K1D(t)隨時間變化的趨勢主要由其結(jié)構(gòu)本身的固有特性所決定;(5)對于同種裂紋尺寸的厚壁筒,加載方式不同對K1D(t)的變化趨勢影響不大。參考文獻1中國航空研究院主編.應(yīng)力強度因子手冊.北京:科學(xué)出版社,19812王志群,朱務(wù)學(xué),查子初.,1989,11(3):62663王志群,何金泉.,1989,13(4):1021074范天佑.,19905WL.echanicsoffracture

11、:Vol4.Leyden:NoordhoffInternationalRelationbetweenCrackDepthsandDynamicStressIntensityFactorsofThick-walledCylinderYangHongweiMaDaweiWangGengxinWangZhiqun(SchoolofMechanics,SchoolofSciences,NUST,Nanjing210094)ABSTRACTAccordingtothebasictheoryoffiniteelementmethods,usingdisplacementmethod,first,thest

12、ationarystressintensityfactorsunderthesameinternalpressurearecomputedforthick2walledcylinderswithfivekindscrackdepths(610mm,1410mm,2010mm,2610mmand3210mm),then,therationalplanofcomputingdynamicstressin2tensityfactorsisdesigned.Thedynamicstressintensityfactorsaboutthefivekindscrackdepthsarecomputedun

13、dertwokindsinternalpressures(impactandsine).Theruleofthestressintensityfactorsofthick2walledcylinderwiththechangeofcrackdepthsandtimeareobtained.Theruleisofsomereferencevaluetothedesignandengineeringofthethick2walledcylinder.KEYWORDScracks,stressintensityfactors,finiteelementmethods,thickwallshells 簡訊我校方大綱教授等4位同志榮獲“紅衫樹園丁獎”為了慶祝第12個教師節(jié),在徐州卷煙廠“紅杉樹教育獎勵委員會”的資助下,江蘇省教委在全省普通高校開