錐角比對(duì)雙錐藥型罩射流成型影響的數(shù)值模擬

郝 鋼，李艷飛，劉裕濤

(1.陸軍裝備部駐沈陽(yáng)地區(qū)第二軍事代表室，遼寧 沈陽(yáng) 110043；2.遼沈工業(yè)集團(tuán)有限公司，遼寧 沈陽(yáng) 110045)

0 引言

聚能射流是一種基于聚能效應(yīng)形成的毀傷元,主要用于打擊裝甲目標(biāo),隨著裝甲防護(hù)能力的不斷提升,單錐藥型罩結(jié)構(gòu)已不能滿足目標(biāo)打擊要求,因此,雙錐藥型罩應(yīng)運(yùn)而生,國(guó)內(nèi)外許多學(xué)者對(duì)雙錐藥型罩的射流成型進(jìn)行了研究。

Walters等[1]建立了雙錐形藥型罩射流的侵徹計(jì)算模型,提出了所形成射流具有雙線性速度分布特性;楊曉紅等[2]通過(guò)試驗(yàn)獲得了雙錐形藥型罩在不同炸高下的靜破甲深度,并利用數(shù)值模擬和理論分析得到了雙錐藥型罩的最佳炸高;趙海平等[3]設(shè)計(jì)了一種錐-錐結(jié)合的藥型罩,通過(guò)改變上、下錐罩的長(zhǎng)度比例和下錐角,分析了組合藥型罩形成射流的規(guī)律;劉宏杰等[4]基于環(huán)形切割聚能裝藥戰(zhàn)斗部及灰關(guān)聯(lián)理論,通過(guò)數(shù)值模擬方法對(duì)環(huán)形雙錐罩聚能裝藥結(jié)構(gòu)進(jìn)行了優(yōu)化設(shè)計(jì);陳闖等[5]利用沖擊波Hugoniot關(guān)系修正伯努利方程,結(jié)合改進(jìn)的PER理論建立了考慮沖擊波、射流速度分布及射流狀態(tài)等影響因素的穿深計(jì)算模型;劉潤(rùn)滋等[6]研究了不同罩高的變壁厚和等壁厚雙錐藥型罩射流成型特性;王子明等[7]通過(guò)優(yōu)化雙錐藥型罩,提出了一種輔助型雙錐藥型罩,并通過(guò)數(shù)值模擬證明了優(yōu)化結(jié)構(gòu)在射流成型方面的優(yōu)勢(shì),達(dá)到了提升射流侵徹性能的目的;王珞冰等[8]對(duì)雙錐藥型罩結(jié)構(gòu)進(jìn)行了優(yōu)化設(shè)計(jì),得到了一種具有較優(yōu)結(jié)構(gòu)的大口徑雙錐藥型罩聚能裝藥結(jié)構(gòu)。由此可以看出,雙錐藥型罩的結(jié)構(gòu)參數(shù)對(duì)射流成型的影響較大,而錐角比是其中重要的影響因素,因此,本文針對(duì)錐角比對(duì)射流成型的影響展開(kāi)研究。

1 數(shù)值模擬模型建立

1.1 雙錐藥型罩聚能裝藥仿真模型

雙錐藥型罩聚能裝藥結(jié)構(gòu)如圖1所示。藥型罩直徑為40 mm,裝藥高度為60 mm,壁厚為1.2 mm,α為上錐角,β為下錐角。選用ANSYS/LS-DYNA進(jìn)行數(shù)值仿真,基于Euler進(jìn)行網(wǎng)格劃分,單元使用多物質(zhì)ALE算法,同時(shí)為了提高計(jì)算效率,僅建立1/4模型,并通過(guò)對(duì)稱邊界條件保證數(shù)值模擬的準(zhǔn)確性。建立的雙錐藥型罩聚能裝藥仿真模型如圖2所示。

圖1 雙錐藥型罩聚能裝藥結(jié)構(gòu)

圖2 建立的雙錐藥型罩聚能裝藥仿真模型

1.2 材料模型及參數(shù)

仿真時(shí)所采用的炸藥為B炸藥,使用High Explosive Burn模型描述,炸藥密度為1.72 g/cm3,爆速為8 425 m/s,爆壓為30.4 GPa;藥型罩材料選擇紫銅,采用Steinberg材料模型和Gruneisen狀態(tài)方程描述,材料密度為8.93 g/cm3,剪切模量為47.7 GPa,狀態(tài)方程中的常數(shù)N、B、F、A分別為0.45、2.83、0.01以及63.5;空氣采用Null材料模型和Linear Polynomial狀態(tài)方程描述,其密度為1.293 g/cm3,狀態(tài)方程中的常數(shù)C4和C5均為0.4。

2 數(shù)值模擬結(jié)果與分析

為了研究錐角比對(duì)雙錐藥型罩射流成型的影響,結(jié)合目前單錐藥型罩和雙錐藥型罩在射流成型研究時(shí)常用的錐角參數(shù),選擇上錐角分別為30°、40°和50°作為基礎(chǔ)角度,下錐角與上錐角的錐角比分別為1∶1.5、1∶1.8、1∶2、1∶2.2和1∶2.5,并以射流頭部最大速度、射流動(dòng)能以及射流長(zhǎng)度作為射流成型質(zhì)量的評(píng)價(jià)指標(biāo),以確定錐角比對(duì)射流成型質(zhì)量影響的趨勢(shì)。

2.1 射流頭部最大速度

通過(guò)數(shù)值模擬得到雙錐藥型罩射流頭部最大速度隨錐角比的變化情況和射流頭部最大速度隨上錐角的變化情況,分別如圖3、圖4所示。

圖3 射流頭部最大速度隨錐角比的變化情況

圖4 射流頭部最大速度隨上錐角的變化情況

由圖3可知:不同上錐角藥型罩射流頭部最大速度的變化趨勢(shì)均為隨著錐角比的增大而降低;上錐角為30°、錐角比為1∶1.5和1∶1.8時(shí),射流頭部最大速度分別為8 088 m/s和8 293 m/s,隨后隨著錐角比的增大射流頭部最大速度逐漸降低,在錐角比為1∶2.5時(shí)降為最小值7 905 m/s。

由圖4可知:同一雙錐藥型罩錐角比下,隨著上錐角的增大,射流頭部最大速度逐漸降低。

2.2 射流動(dòng)能變化

射流動(dòng)能隨雙錐藥型罩錐角比的變化情況和隨雙錐藥型罩上錐角的變化情況如圖5和圖6所示。

由圖5可知:當(dāng)雙錐藥型罩上錐角為30°時(shí),射流動(dòng)能隨著錐角比的增加而降低,在錐角比為1∶1.5時(shí),射流動(dòng)能最大值為7 kJ;上錐角為40°時(shí),射流動(dòng)能變化趨勢(shì)為先減小后增大,最小值出現(xiàn)在錐角比為1∶2時(shí),為6.51 kJ;上錐角為50°時(shí),射流動(dòng)能先隨著錐角比的增大而不斷增大,在錐角比為1∶2.2時(shí)達(dá)到最大值6.97 kJ,此時(shí)射流最大動(dòng)能接近上錐角為30°時(shí)的最大值,隨后射流動(dòng)能下降。由此可以看出,錐角比使藥型罩結(jié)構(gòu)發(fā)生了改變,從而影響了藥型罩中爆轟波的傳播,使不同藥型罩所形成的射流速度分布不同,進(jìn)而影響到射流動(dòng)能。

由圖6可知:當(dāng)錐角比為1∶1.5和1∶1.8時(shí),隨著上錐角的增大,射流動(dòng)能逐漸降低;當(dāng)錐角比大于1∶2后,上錐角為50°時(shí),射流具有的動(dòng)能最大。

2.3 射流長(zhǎng)度和有效射流長(zhǎng)度

雙錐藥型罩上錐角為30°、40°和50°時(shí)在不同錐角比下的射流長(zhǎng)度和有效射流長(zhǎng)度如表1所示。由表1可知:射流長(zhǎng)度和有效射流長(zhǎng)度均隨著錐角比的增加而降低,在一般情況下有效射流占有大比重的射流動(dòng)能,因此錐角比的增加會(huì)對(duì)射流的侵徹能力產(chǎn)生不利影響。

表1 不同錐角比下的射流長(zhǎng)度和有效射流長(zhǎng)度 mm

3 結(jié)論

通過(guò)對(duì)雙錐藥型罩聚能裝藥結(jié)構(gòu)進(jìn)行數(shù)值模擬,得到以下結(jié)論:

(1) 上錐角越小,射流頭部最大速度越大,錐角比越大,射流頭部最大速度越低。

(2) 上錐角較小時(shí)產(chǎn)生的射流動(dòng)能隨錐角比增加而降低,上錐角較大時(shí)產(chǎn)生的射流動(dòng)能隨錐角比增加而增加。

(3) 隨著錐角比的增加,射流長(zhǎng)度和有效射流長(zhǎng)度均降低,從而會(huì)對(duì)射流的侵徹能力產(chǎn)生不利影響。