王鈺婷, 黃正祥, 祖旭東, 馬彬, 肖強(qiáng)強(qiáng), 賈鑫
(南京理工大學(xué) 機(jī)械工程學(xué)院, 江蘇 南京 210094)
聚能裝藥具有強(qiáng)大的威力,且不限制其搭載平臺(tái)到達(dá)目標(biāo)時(shí)的存速,是空中武器平臺(tái)理想的戰(zhàn)斗部。非圓截面外形具有良好的氣動(dòng)力學(xué)特性,是空中武器平臺(tái)常見的彈體布局形式。因此,非圓截面空中武器平臺(tái)對(duì)聚能裝藥具有強(qiáng)烈的應(yīng)用需求,采用傳統(tǒng)旋轉(zhuǎn)對(duì)稱結(jié)構(gòu)聚能裝藥,雖然保證了破甲威力,但較低的炸藥裝填率使得破片威力和沖擊波威力受到限制,且增加了空間布局難度。目前,聚能裝藥的研究以旋轉(zhuǎn)對(duì)稱結(jié)構(gòu)為主,其射流成型及侵徹理論難以滿足非圓截面聚能裝藥,迫切需要開展非圓截面聚能裝藥射流成型與侵徹方面的研究。
截止目前,已有部分學(xué)者開展了非圓截面聚能裝藥研究。Li等、李硯東等研究發(fā)現(xiàn)半球形藥型罩加裝類圓矩形截面裝藥的聚能裝藥在多點(diǎn)起爆下,將形成一種水滴狀的高速成型彈丸。Cullis等研究的藥型罩為金字塔結(jié)構(gòu)的方形截面聚能裝藥和Wang等研究的藥型罩為旋轉(zhuǎn)對(duì)稱結(jié)構(gòu)的方形截面聚能裝藥在中心點(diǎn)起爆下形成的射流都具有如下特征:靠近杵體的射流截面形狀為星形,截面位置越靠近射流頭部,其形狀越趨近于圓形。Stewart等在多模戰(zhàn)斗部中嵌入了非旋轉(zhuǎn)對(duì)稱結(jié)構(gòu)的聚能裝藥,在偏心起爆下形成偏斜的射流。Wang等通過脈沖X光攝影試驗(yàn)研究了藥型罩為旋轉(zhuǎn)對(duì)稱結(jié)構(gòu)的方形截面聚能裝藥的射流成型特性,發(fā)現(xiàn)射流由凝聚部分和非凝聚部分組成,非凝聚部分的射流由一束沿著對(duì)稱軸運(yùn)動(dòng)的主體流體和四束沿著對(duì)稱面離軸運(yùn)動(dòng)的離軸流體組成。王鈺婷等研究了藥型罩為旋轉(zhuǎn)對(duì)稱結(jié)構(gòu)的橢圓形截面聚能裝藥的射流成型及對(duì)靶板的毀傷特性,發(fā)現(xiàn)射流在拉伸運(yùn)動(dòng)后期出現(xiàn)非凝聚現(xiàn)象,非凝聚的射流由關(guān)于截面長軸面對(duì)稱分布的兩束流體組成,射流侵徹能力大大降低。
目前,四邊形截面彈體布局已廣泛應(yīng)用于各國研制的新型彈箭上,如德國金牛座KEPD系列導(dǎo)彈,因此,相應(yīng)四邊形截面聚能戰(zhàn)斗部的研究十分必要。本文通過脈沖X光攝影試驗(yàn)研究等腰梯形截面聚能裝藥的射流成型特性,通過開展炸高為80 mm的侵徹深度(DOP)試驗(yàn)研究等腰梯形截面聚能裝藥對(duì)靶板的毀傷特性,基于數(shù)值模擬給出等腰梯形截面聚能裝藥形成準(zhǔn)直、凝聚射流的方法。研究結(jié)論對(duì)于等腰梯形截面聚能戰(zhàn)斗部的應(yīng)用具有參考價(jià)值。
等腰梯形截面聚能裝藥結(jié)構(gòu)如圖1所示,為旋轉(zhuǎn)成體藥型罩加裝等腰梯形截面裝藥。藥型罩采用錐角為60°的單錐結(jié)構(gòu),大端直徑54.2 mm,高度42 mm,厚度1 mm。裝藥截面銳角60°,裝藥截面內(nèi)切圓直徑56 mm,裝藥高度均為73.3 mm。、分別為裝藥剖面的裝藥半徑。
圖1 聚能裝藥Fig.1 Shaped charge
藥型罩材料為紫銅,炸藥類型為鑄裝B炸藥。為簡(jiǎn)化加工難度,鑄藥過程中采用3 mm的尼龍殼體以保證裝藥特征尺寸。尼龍殼體作為低阻抗材料,對(duì)爆轟波的傳播歷程沒有影響,可視作無殼聚能裝藥。
脈沖X光攝影試驗(yàn)布局如圖2所示,試驗(yàn)采用兩臺(tái)脈沖X光機(jī)組合拍攝,X射線管A和X射線管B呈45°匯交。聚能裝藥垂直放置于兩臺(tái)X射線管的匯交區(qū)域。試驗(yàn)過程中,設(shè)置X射線管B先工作,X射線管A后工作,獲得射流在不同時(shí)刻的X光照片。
圖2 脈沖X光攝影試驗(yàn)布局示意圖Fig.2 X-ray experimental configuration and setup
脈沖X光攝影試驗(yàn)的各裝置相對(duì)位置布設(shè)如圖3所示,其中,X射線管A拍攝等腰梯形截面聚能裝藥的對(duì)稱面視圖;X射線管B設(shè)置于截面下底一側(cè),拍攝45°側(cè)視圖。
圖3 各裝置相對(duì)位置布設(shè)圖Fig.3 Relative position of configurations in X-ray experiment
DOP試驗(yàn)現(xiàn)場(chǎng)布置如圖4所示,包括雷管、雷管座、擴(kuò)爆藥、聚能裝藥、炸高筒以及靶板。炸高為80 mm,炸高筒由硬紙板制成;靶板材料為45號(hào)鋼。試驗(yàn)前,標(biāo)記靶板與聚能裝藥的相對(duì)位置,其中靶板的標(biāo)記線對(duì)應(yīng)截面的下底中點(diǎn)。
圖4 DOP試驗(yàn)現(xiàn)場(chǎng)布置圖Fig.4 Experimental setup of DOP test
起爆后27.7 μs時(shí)刻射流的45°側(cè)視圖以及起爆后38.2 μs時(shí)刻射流的對(duì)稱面視圖如圖5所示。通過兩時(shí)刻射流的空間位置可以計(jì)算得到射流的頭部平均速度為6.41 km/s,尾部平均速度為3.24 km/s。從圖5中可以觀察到:射流明顯偏離軸線,隨著射流的軸向拉伸,偏離軸線的最大距離由10 mm增加至17 mm;射流后半部分呈非凝聚態(tài),非凝聚的射流由兩束偏離藥型罩軸線的流體構(gòu)成,一束與藥型罩軸線距離較遠(yuǎn)且尺寸較大,一束與藥型罩軸線距離較近且尺寸較小,隨著射流的軸向拉伸,兩束流體的間距增加。
圖5 脈沖X光圖像Fig.5 X-ray photos of jets
圓形截面聚能裝藥的脈沖X光圖像顯示形成的崩落環(huán)呈圓環(huán)狀,這歸因于崩落環(huán)上的藥型罩微元壓垮速度相同。由圖5可知,等腰梯形截面聚能裝藥形成的崩落環(huán)45°側(cè)視圖呈不對(duì)稱花瓣?duì)睿瑢?duì)稱面視圖呈線性狀,表明崩落環(huán)上的藥型罩微元具有不同的壓垮速度。
圖6所示為DOP試驗(yàn)后形成的靶板入孔照片,具體侵徹入孔數(shù)據(jù)如表1所示。在靶板表面觀察到一個(gè)初始侵徹孔和一個(gè)偏移侵徹孔;偏移侵徹孔分布在沿著聚能裝藥對(duì)稱面的上底一側(cè)。初始侵徹孔平均直徑為27.5 mm;第1發(fā)試驗(yàn)由于靶板缺失,難以測(cè)量偏移侵徹孔數(shù)據(jù),第2發(fā)試驗(yàn)的偏移侵徹孔直徑為16.5 mm,侵徹孔中心偏移距離為19.4 mm。
圖6 靶板入孔照片F(xiàn)ig.6 Photographs of penetrating hole on target
表1 DOP試驗(yàn)靶板入孔數(shù)據(jù)
圖7所示為DOP試驗(yàn)后形成的侵徹通道照片,具體侵徹深度數(shù)據(jù)如表2所示。初始侵徹通道平均深度為50 mm,偏移侵徹通道平均深度為47 mm。藥型罩結(jié)構(gòu)及炸藥類型相同的直徑為56 mm圓形截面聚能裝藥的侵徹深度為150 mm,與之相比,等腰梯形截面聚能裝藥的侵徹深度下降了66.7%。
圖7 侵徹通道照片F(xiàn)ig.7 Photographs of penetration channel
表2 侵徹深度數(shù)據(jù)
由圖7可知,炸高為80 mm時(shí),侵徹通道壁面較為光滑,表明射流在該炸高下并未明顯斷裂失穩(wěn)。由脈沖X光攝影試驗(yàn)可知,等腰梯形截面聚能裝藥的射流存在橫向速度。橫向速度較小的射流可以到達(dá)侵徹孔底部繼續(xù)增加侵徹深度,橫向速度較大的射流在到達(dá)靶板前,積累了較大的橫向位移,到達(dá)靶板時(shí)橫向位移大于初始侵徹孔直徑,形成偏移侵徹通道,這是造成侵徹深度下降的主要原因。
為更清晰地認(rèn)識(shí)等腰梯形截面聚能裝藥的射流成型,采用ANSYS-LSDYNA有限元軟件建立等腰梯形截面聚能裝藥的三維有限元模型,如圖8所示。由圖8可知:有限元模型由炸藥、藥型罩和空氣3部分組成;采用變尺寸網(wǎng)格,藥型罩母線方向網(wǎng)格尺寸為1.0~1.2 mm,炸藥、藥型罩和空氣采用歐拉共節(jié)點(diǎn)網(wǎng)格,起爆方式為上表面中心點(diǎn)起爆。
圖8 數(shù)值計(jì)算模型Fig.8 Numerical calculation model
炸藥類型為B炸藥,采用高能炸藥爆轟模型及Jones-Wilkins-Lee(JWL)狀態(tài)方程描述,主要參數(shù)如表3所示。藥型罩材料為紫銅,采用Johnson-Cook材料模型和Gruneisen狀態(tài)方程描述,主要參數(shù)如表4所示。
表3 B炸藥參數(shù)[12]
表4 紫銅材料參數(shù)[13]
數(shù)值模擬結(jié)果與脈沖X光圖像對(duì)比如圖9所示,具體參數(shù)對(duì)比如表5所示。由圖9可知數(shù)值模擬與試驗(yàn)的射流形態(tài)相仿:數(shù)值模擬得到射流頭部速度和尾部速度分別為6.23 km/s和3.01 km/s,與試驗(yàn)結(jié)果誤差分別為-2.4%和-7.5%;起爆后27.7 μs和38.2 μs時(shí)刻,射流長度的數(shù)值模擬值與試驗(yàn)值吻合程度較好,最大誤差為-4.2%;38.2 μs時(shí)刻,數(shù)值模擬得到射流偏離軸線的最大距離為11 mm,與試驗(yàn)結(jié)果誤差為35.3%,且數(shù)值模擬的位置比試驗(yàn)更接近射流尾部。分析原因,可能是等腰梯形截面聚能裝藥加工難度較大,使得試驗(yàn)中射流的橫向速度分布與表征理想情況的數(shù)值模擬結(jié)果存在一定的差異。但總體上數(shù)值模擬可較好地反映等腰梯形截面聚能裝藥的射流成型過程。
圖9 數(shù)值模擬結(jié)果(左)與脈沖X光照片(右)對(duì)比Fig.9 Comparison of simulated results (left) and X-ray images (right)
表5 數(shù)值模擬及X光試驗(yàn)結(jié)果對(duì)比
3.3.1 爆轟波及表面稀疏波的傳播
以距離裝藥上表面10 mm的裝藥截面為例,分析裝藥截面上的壓力變化過程,如圖10所示。由圖10可見:裝藥起爆后,爆轟波以球形波陣面?zhèn)鞑ァF鸨? μs時(shí)刻,爆轟波還未傳播至裝藥截面的邊緣,裝藥截面上爆轟產(chǎn)物壓力呈旋轉(zhuǎn)對(duì)稱分布;起爆后4 μs時(shí)刻,爆轟波傳播至剖面半徑較小的裝藥邊緣,表面稀疏波向爆轟產(chǎn)物傳入,使得爆轟產(chǎn)物的壓力降低;截面銳角平分線附近的裝藥半徑較大,爆轟波總是較晚傳播至截面銳角平分線附近的裝藥邊緣,截面銳角平分線附近爆轟產(chǎn)物的壓力總是減小的更慢;起爆后5 μs時(shí)刻,截面銳角平分線附近的爆轟產(chǎn)物依然具有較高的壓力。
圖10 距離裝藥上表面10 mm裝藥截面上的壓力變化Fig.10 Pressure change on the charge cross-section at 10 mm from the upper surface
3.3.2 藥型罩壓垮
起爆后8 μs和14 μs時(shí)刻藥型罩表面壓力云圖如圖11所示。由圖11可見:8 μs時(shí),表面稀疏波還未傳播至藥型罩表面,藥型罩表面壓力呈旋轉(zhuǎn)對(duì)稱分布;14 μs時(shí),表面稀疏波已陸續(xù)到達(dá)藥型罩表面,表面稀疏波傳入截面銳角平分線附近爆轟產(chǎn)物的時(shí)間總是較晚,則截面銳角平分線附近的藥型罩表面壓力也總是減小的更慢。
圖11 起爆后8 μs和14 μs時(shí)刻藥型罩表面壓力云圖 Fig.11 Pressure on liner surface at 8 μs and 14 μs after detonation
起爆后8 μs和14 μs時(shí)刻藥型罩微元速度云圖如圖12所示。由圖12可見:藥型罩壓垮初期,由于藥型罩表面壓力呈旋轉(zhuǎn)對(duì)稱分布,藥型罩微元速度也呈旋轉(zhuǎn)對(duì)稱分布;藥型罩壓垮后期,由于截面銳角平分線附近的藥型罩表面壓力總是減小得更慢,則截面銳角平分線附近的藥型罩微元速度總是增加得更快。
圖12 起爆后8 μs和14 μs藥型罩微元速度云圖Fig.12 Collapse velocity of liner elements at 8 μs and 14 μs after detonation
3.3.3 射流成型
起爆后不同時(shí)刻等腰梯形截面聚能裝藥射流的形態(tài)及速度分布如圖13所示。由圖13可見:炸藥起爆后20 μs時(shí),射流基本形成,觀察到射流對(duì)稱面上的尺寸略微大于其垂直面上的尺寸;射流速度呈梯度分布;射流存在沿著對(duì)稱面的橫向速度,且橫向速度具有速度梯度;隨著時(shí)間增長,射流沿軸向拉長的同時(shí),沿對(duì)稱面的尺寸也逐漸增加,最后發(fā)生橫向頸縮,出現(xiàn)非凝聚現(xiàn)象;非凝聚的射流由一束橫向速度較大的流體和一束橫向速度較小的流體構(gòu)成,橫向速度的差異使得兩束流體的間距隨時(shí)間增加而增加。
圖13 射流形態(tài)及速度分布圖Fig.13 Morphology and velocity distribution of jet
射流的成型特性與藥型罩壓垮特性息息相關(guān)。對(duì)于傳統(tǒng)旋轉(zhuǎn)對(duì)稱結(jié)構(gòu)聚能裝藥,同一截面的藥型罩微元以相同的壓垮速度到達(dá)軸線的同一點(diǎn)碰撞形成射流微元,因此射流沿著軸線運(yùn)動(dòng)。等腰梯形截面聚能裝藥射流的特性也可以從藥型罩壓垮特性中得到解釋。
等腰梯形截面聚能裝藥的藥型罩壓垮是一個(gè)三維問題,為便于闡述,將藥型罩壓垮過程投影到二維截面上開展描述,如圖14所示。圖14中,、為沿裝藥截面對(duì)稱軸的坐標(biāo)軸,、為垂直于裝藥截面對(duì)稱軸的坐標(biāo)軸,坐標(biāo)軸與數(shù)值計(jì)算模型相同。區(qū)域、′內(nèi),剖面上裝藥半徑、相同,同一剖面藥型罩微元經(jīng)歷相同的壓垮歷程;區(qū)域、′和區(qū)域、′內(nèi),剖面上裝藥半徑、不同,同一剖面的藥型罩微元的壓垮歷程存在差異。
圖14 藥型罩壓垮過程二維投影示意圖Fig.14 Schematic diagram of two-dimensional projection of liner collapse
靠近頂部的藥型罩微元碰撞形成射流的頭部,表面稀疏波傳播至這些藥型罩微元表面需要較長的時(shí)間。運(yùn)動(dòng)初始階段,表面稀疏波還未傳播至這些藥型罩微元表面,藥型罩微元表面壓力的衰減僅來源于藥型罩微元運(yùn)動(dòng)對(duì)其表面爆轟產(chǎn)物的稀疏擾動(dòng),藥型罩具有旋轉(zhuǎn)對(duì)稱結(jié)構(gòu)。因此,同一截面的藥型罩微元具有相同的壓垮速度,在空間上呈圓環(huán)分布。運(yùn)動(dòng)最后階段,表面稀疏波陸續(xù)傳播至這些藥型罩微元表面,傳入時(shí)間的差異使得壓垮速度出現(xiàn)輕微差異,在空間上的分布偏離圓環(huán)形狀,如圖14(a)所示,但由于存在差異的時(shí)間較短,對(duì)形成的射流微元特性影響較小,射流頭部的橫向速度較小,橫向速度梯度不明顯。
表面稀疏波很快陸續(xù)到達(dá)靠近藥型罩底部的藥型罩微元,運(yùn)動(dòng)初始階段,同一截面的藥型罩微元壓垮歷程就開始呈現(xiàn)差異,如圖14(b)所示。區(qū)域、′內(nèi)的藥型罩微元壓垮至軸線碰撞形成射流時(shí),區(qū)域、′以及區(qū)域、′內(nèi)的藥型罩微元與軸線還具有一定距離,此外,同一截面的藥型罩微元最終到達(dá)軸線的位置也存在差異,這些因素都將使得形成的射流微元具有不可忽視的橫向速度。當(dāng)藥型罩微元壓合速度差異較大時(shí),碰撞時(shí)間及碰撞位置的差異明顯,射流微元橫向速度具有速度梯度。
越靠近藥型罩底部,表面稀疏波傳播至藥型罩微元表面的時(shí)間越早,同一截面藥型罩微元壓垮速度呈現(xiàn)差異的時(shí)間也越早,射流微元的橫向速度越大,橫向速度梯度也越大。因此,觀察到射流尾部的非凝聚程度最為明顯。
3.3.4 內(nèi)切圓尺寸對(duì)藥型罩壓垮的影響
由3.3.1節(jié)~3.3.3節(jié)可知,表面稀疏波傳播至藥型罩表面的時(shí)間直接影響藥型罩的壓垮特性,進(jìn)而影響射流成型特性??梢酝ㄟ^增加裝藥截面內(nèi)切圓直徑與藥型罩口徑的相對(duì)大小來延遲表面稀疏波傳入的時(shí)間,以此弱化表面稀疏波傳入時(shí)間的差異對(duì)藥型罩壓垮及射流成型的影響??赏ㄟ^保持藥型罩結(jié)構(gòu)不變、增加裝藥截面內(nèi)切圓直徑,或保持裝藥截面內(nèi)切圓直徑、減小藥型罩口徑來實(shí)現(xiàn)上述目的。本研究將分析藥型罩結(jié)構(gòu)不變,裝藥截面內(nèi)切圓直徑的增加對(duì)藥型罩壓垮以及射流成型特性的影響。
圖15所示為等腰梯形截面聚能裝藥的藥型罩結(jié)構(gòu)不變、裝藥截面內(nèi)切圓直徑不同時(shí),藥型罩微元在起爆后12 μs和18 μs時(shí)刻的速度分布。由圖15可以看到,當(dāng)裝藥截面內(nèi)切圓直徑大于72.8 mm(56 mm的1.3倍)時(shí),12 μs時(shí)刻藥型罩微元速度呈旋轉(zhuǎn)對(duì)稱分布。當(dāng)裝藥截面內(nèi)切圓直徑增加至 84.0 mm(56 mm的1.5倍)時(shí),18 μs時(shí)刻藥型罩微元速度也呈旋轉(zhuǎn)對(duì)稱分布??梢哉J(rèn)為,裝藥截面內(nèi)切圓直徑為84 mm時(shí),藥型罩的壓垮具有旋轉(zhuǎn)對(duì)稱性。
圖15 內(nèi)切圓尺寸對(duì)藥型罩壓垮的影響Fig.15 Influence of the inscribed circle size on liner collapse
3.3.5 內(nèi)切圓尺寸對(duì)射流成型的影響
圖16展示了等腰梯形截面聚能裝藥的藥型罩結(jié)構(gòu)不變、裝藥截面內(nèi)切圓直徑不同時(shí),不同時(shí)刻射流的形態(tài)及橫向速度分布。由圖16可以看到:裝藥截面內(nèi)切圓直徑為56 mm時(shí),從起爆后30 μs到50 μs,射流最大橫向速度的變化區(qū)間為[570 m/s,430 m/s];裝藥截面內(nèi)切圓直徑增加至61.6 mm(56 mm的1.1倍)時(shí),射流的最大橫向速度變化區(qū)間為[325 m/s,266 m/s],射流橫向速度有所降低。此外,由圖16還可看出:射流的橫向速度梯度減小,射流的非凝聚現(xiàn)象得到明顯改善;裝藥截面內(nèi)切圓直徑增加至67.2 mm(56 mm的1.2倍)時(shí),射流的最大橫向速度變化區(qū)間降低至[170 m/s,159 m/s], 50 μs 內(nèi)射流一直呈凝聚態(tài);隨著裝藥截面內(nèi)切圓直徑的進(jìn)一步增加,射流橫向速度進(jìn)一步減小,裝藥截面內(nèi)切圓直徑增加至84.0 mm(56 mm的1.5倍)時(shí),射流的最大橫向速度變化區(qū)間為[50 m/s,25 m/s],可以認(rèn)為射流不具備明顯的橫向速度。
圖16 內(nèi)切圓尺寸對(duì)射流成型的影響Fig.16 Influence of the inscribed circle size on jet formation
本文以截面內(nèi)切圓直徑為56 mm,截面銳角為60°的等腰梯形截面聚能裝藥為研究對(duì)象,開展了脈沖X光攝影試驗(yàn),以及炸高為80 mm的DOP試驗(yàn)。利用ANSYS/LSDYNA有限元軟件對(duì)相關(guān)等腰梯形截面聚能裝藥的射流成型進(jìn)行了數(shù)值模擬。得出主要結(jié)論如下:
1)等腰梯形截面聚能裝藥形成的射流具有橫向速度,且橫向速度呈梯度分布。除射流頭部外,其余部分在拉伸運(yùn)動(dòng)后期出現(xiàn)非凝聚現(xiàn)象。非凝聚的射流由兩束偏離軸線的流體構(gòu)成,一束橫向速度較大且尺寸較大,一束橫向速度較小且尺寸較小。
2)炸高為80 mm時(shí),等腰梯形截面聚能裝藥的平均侵徹深度為50 mm,射流的橫向速度明顯降低了射流的侵徹能力。
3)保持藥型罩結(jié)構(gòu)不變、增加裝藥截面內(nèi)切圓直徑可以延遲表面稀疏波傳入的時(shí)間,從而弱化裝藥結(jié)構(gòu)非旋轉(zhuǎn)對(duì)稱性對(duì)射流成型的影響。對(duì)于本文的藥型罩結(jié)構(gòu),裝藥截面內(nèi)切圓直徑增加至67.2 mm(56 mm 的1.2倍時(shí))時(shí),射流的非凝聚現(xiàn)象消失。裝藥截面內(nèi)切圓直徑增加至84.0 mm(56 mm的1.5倍時(shí))時(shí),射流不具備明顯的橫向速度。