劉文祥,譚書舜,劉冠蘭
(西北核技術(shù)研究所,陜西 西安710024)
美國Lawrence Livermore國家實驗建立了JASPER裝置[1],該裝置的一項關(guān)鍵技術(shù)是超快速密封技術(shù),利用包裹在管道外面的炸藥爆轟產(chǎn)生向內(nèi)聚合的爆炸載荷壓縮圓管,進而實現(xiàn)管道在100μs內(nèi)封閉。內(nèi)聚滑移爆轟壓縮圓管還在金屬管道焊接等領(lǐng)域得到應用[2-5]。工程分析在相關(guān)研究中起著非常重要的作用,可以提供快捷的參考,也可以加深對物理問題的認識。針對內(nèi)聚爆轟壓縮管道的研究較少,邵炳璜等[6]做了這方面的工作,但管道壓縮量較大時,分析解與數(shù)值模擬結(jié)果、實驗結(jié)果相差較大。
本文中,對內(nèi)聚滑移爆轟加載下圓管壓縮運動過程進行工程分析,分析結(jié)果由數(shù)值模擬結(jié)果驗證。
圖1為內(nèi)聚滑移爆轟壓縮圓管的示意圖,管道外面均勻包裹炸藥,炸藥由端部截面起爆,炸藥進行滑移爆轟,壓縮管道進行封閉。設r00為炸藥初始外半徑,r10為管壁的初始外半徑,r20為管壁的初始內(nèi)半徑,r1、r2為管壁的流動外半徑和內(nèi)半徑。ρ1為管道材料密度,ρ0為炸藥裝藥密度。
內(nèi)聚滑移爆轟壓縮圓管工程模型的假定條件:(1)忽略應力波在管道壁內(nèi)的傳播和反射過程,即管壁內(nèi)應力波傳播時間為零;(2)管道外壁爆炸載荷隨時間呈指數(shù)衰減,衰減特征時間與管壁最大徑向速度呈指數(shù)減少的關(guān)系;(3)僅考慮管道材料徑向運動,忽略軸向運動;(4)管道材料不可壓縮,為理想剛塑性材料,在管道壓縮瞬間周向應力達到屈服值;(5)管壁內(nèi)速度呈線性分布;(6)管道不發(fā)生層裂。
炸藥進行穩(wěn)定爆轟后,炸藥CJ爆轟壓力為:
問題的處理遵循從簡單到復雜的方法。首先假設管道壁為剛性壁面,以便得到未受管壁運動干擾情況下的管壁外載荷情況。根據(jù)爆炸載荷衰減特征,認為管道外壁面爆炸載荷隨時間呈指數(shù)衰減,則管壁外壓力和沖量可以表示為:
式中:T0為衰減特征時間。通過數(shù)值模擬研究發(fā)現(xiàn),衰減特征時間與裝藥厚度的關(guān)系為線性關(guān)系:
式中:a為常數(shù),只與炸藥性質(zhì)有關(guān)。
爆炸載荷驅(qū)動下管道向內(nèi)壓縮,這使衰減特征時間T0更短。由于管道壁運動速度在極短時間(相對于管道密閉時間來說)內(nèi)到達峰值,可以認為衰減特征時間與管道外壁速度峰值存在關(guān)聯(lián),修正后的衰減特征時間為:
式中:b、c和d為常數(shù)。利用管道外壁最大速度對爆炸載荷衰減特征時間進行修正,修正公式系數(shù)與管道材料有關(guān)。
設管道壁內(nèi)速度呈線性分布,存在關(guān)系:
式中:v1(t)為管壁壓縮過程中外表面的速度,v2(t)為管壁內(nèi)表面的速度。由管道材料的不可壓縮性得到:
進而得到:
由式(6)、(8),得到管道的動量表示式:
爆炸載荷加載下管道獲得動量,這些動量在后續(xù)的管道壓縮過程中不斷被管道變形的內(nèi)部應力(主要是周向應力)消耗掉。
假設管道材料為理想剛塑性材料,管道周向應力瞬間達到材料屈服強度Y0。為了簡化問題,可以把內(nèi)部應力等效為虛擬的、加載在管道外壁的載荷,方向與爆炸載荷相反:
內(nèi)部應力消耗的管道動量為:
依據(jù)動量守恒關(guān)系,爆炸載荷施加在管道外壁的沖量等于管道的動量與管道內(nèi)部應力消耗的動量之和。求解式(1)~(12),可得出管道的運動過程:
對爆炸壓縮管道的過程進行數(shù)值模擬,采用二維軸對稱模型(見圖2)。模型中炸藥為PETN,爆轟產(chǎn)物采用JWL狀態(tài)方程,初始密度ρ0=1 770kg/m3,A=617GPa,B=1.692 6GPa,R1=4.4,R2=1.2,ω=0.3,比內(nèi)能e0=5.7MJ/kg。管道材料為紫銅,初始密度ρ0=8 900kg/m3,屈服強度為120MPa。管道采用拉氏網(wǎng)格,網(wǎng)格尺寸為1.2mm×1.2mm,炸藥和空氣采用歐拉網(wǎng)格,網(wǎng)格尺寸為1mm×1mm,爆炸產(chǎn)物、空氣與管道之間相互作用采用流固耦合模擬。
圖2 數(shù)值模擬模型Fig.2 Numerical model
圖3為數(shù)值模擬和本文工程分析得到的管壁運動過程,vi為管道中部速度。數(shù)值模擬結(jié)果顯示,管壁速度曲線存在高頻振蕩,這是由上升時間極短的爆炸波加載在管壁形成應力波、應力波在壁內(nèi)來回反射引起的。數(shù)值模擬結(jié)果還顯示,管道內(nèi)壁速度為零之后外壁速度才為零,這是應力波到達內(nèi)外壁的時間有先后造成的。工程模型忽略了應力波在管道壁內(nèi)的傳播和反射過程,因此管壁運動曲線不存在這種振蕩,管內(nèi)外壁速度同時為零。數(shù)值模擬結(jié)果中管壁速度曲線顯示速度出現(xiàn)小于零的情況,這是數(shù)值模擬中管道內(nèi)壁碰撞后反彈造成的。從圖3中可以看出,數(shù)值模擬結(jié)果與工程分析結(jié)果基本吻合。
圖3 管壁運動過程的數(shù)值模擬結(jié)果和工程分析解Fig.3 Numerical results and analytical results of tube movement
通過工程模型可以快速得到不同參數(shù)管道封閉所需裝藥量以及管道封閉所用時間等,反過來也可以得到裝藥量與管道運動速度的關(guān)系。這里以內(nèi)壁絕對徑向撞擊速度500m/s(單邊內(nèi)壁徑向速度為250m/s)為標準,即認為該速度下管道有效封閉,管道封閉所需炸藥量和管道封閉所用時間見圖4~5。
圖4 管道封閉所需炸藥量Fig.4 Relationship between explosives needed to close tube and tube thickness
圖5 管道封閉所需時間Fig.5 Relationship between tube closing time and tube thickness
擬合數(shù)據(jù),得到管道封閉所需炸藥量和管道封閉所用時間的表達式為:
式中:h1為炸藥厚度,h2為管道壁厚,r20為管道內(nèi)半徑,T為管道封閉時間,D為炸藥爆速。
對內(nèi)聚滑移爆轟壓縮管道的管壁運動過程進行工程分析,并利用數(shù)值模擬進行了驗證,通過分析可以得出以下結(jié)論:
(1)工程分析結(jié)果與數(shù)值模擬結(jié)果相吻合,驗證工程模型是合理的,同時也說明模型的假設條件是物理過程中存在的真實條件或者近似條件。
(2)工程模型得到了不同管道參數(shù)下管道封閉所需炸藥量和封閉時間的表達式,這些結(jié)論可作為相關(guān)研究的參考。
[1]Braddy R W,Homes N,Konrad C H,et al.Joint actinide shock physics experimental research(JASPER)facility overview[C]∥52nd Meeting of the Aeroballistic Range Association.Québec City,Québec,Canada,2001.
[2]隋國發(fā),李金山,馬貝,等.炸藥量對雙層圓管爆炸焊接影響的數(shù)值模擬[J].材料科學與工藝,2010,18(6):741-745.Sui Guo-fa,Li Jin-shan,Ma Bei,et al.Simulation for the influences of explosive amount on explosive welding of two-layer tubes[J].Materials Science and Technology,2010,18(6):741-745.
[3]陳軍.滑移爆轟驅(qū)動金屬圓管及對碰區(qū)的研究[D].四川綿陽:中國工程物理研究院,2001.
[4]李雪梅,金孝剛,李大紅,等.液體炸藥滑移內(nèi)爆加載下鋼管的變形與層裂破壞研究[J].爆炸與沖擊,2003,23(6):523-528.Li Xue-mei,Jin Xiao-gang,Li Da-h(huán)ong,et al.Deformation and spallation of cylindrical steel tube loaded by sliding implosion of liquid explosive[J].Explosion and Shock Waves,2003,23(6):523-528.
[5]李雪梅,汪小松,王鵬來,等.內(nèi)爆和柱殼條件下無氧銅的層裂特性[J].爆炸與沖擊,2009,29(2):162-166.Li Xue-mei,Wang Xiao-song,Wang Peng-lai,et al.Spall of cylindrical copper by converging sliding detonation[J].Explosion and Shock Waves,2009,29(2):162-166..
[6]薛丙璜,張凱.爆炸焊接原理及其應用[M].遼寧大連:大連工學院出版社,1987.