彈丸侵徹凍土實驗及仿真

王麒杰,王 健,劉 銳

(1.南京理工大學 能源與動力工程學院,江蘇 南京 210094;2.湖北航天飛行器研究所,湖北 武漢 430000)

凍土指溫度在0 ℃或0 ℃以下,并含有冰晶體的各種巖土和土壤,其由冰、未凍水、礦物材料和孔隙氣體構成[1]。我國受凍土影響區域約占國土面積的22%,主要分布在青藏高原和東北的高緯度地區,如大小興安嶺等區域[2]。在凍土區域,不論開展民事工程還是軍事工程都需要對凍土力學性質有一定的研究。開展民事工程前需要考慮對凍土進行挖掘、鉆探等作業需求;開展軍事工程前需要考慮彈丸對凍土的侵徹作用,分析凍土對防護工程的抗侵徹能力。因此,針對凍土結構的毀傷研究,對資源開采、建筑工程還是高寒環境作戰都具有重要意義。

在凍土侵徹方面,馬芹永[3]分析了影響凍土爆破性的因素,研究了不同溫度凍土的爆破塊度分布規律。劉志杰等[4]通過大量模型試驗,得到了凍土爆破指標與凍土溫度、含水率等主要參數的關系。陳翰等[5]采用理想彈塑性模型對凍土進行數值模擬,對比分析紫銅和鋁合金兩種材料所形成的射流對凍土材料的侵徹作用,研究發現不同炸高下桿式射流對凍土靶板侵徹效果影響較小。王軍[6]通過LS-DYNA對桿式射流侵徹凍土進行了數值仿真。

本文以7.62 mm制式彈丸為對象,開展了彈丸侵徹凍土的實驗研究和有限元仿真,分析凍土溫度和含水量的影響,對試驗結果與數值計算結果進行對比分析,得到凍土材料在不同條件下侵徹規律。

1 彈丸侵徹射擊試驗

1.1 試驗準備

由于凍土試樣的特殊性,需要嚴格保證試驗時的環境溫度應該低于0 ℃,防止溫度較高導致試樣融化以及產生溫度應力。因此,實驗在室內靶道開展,既可以調節試驗場地的溫度,又便于放置各類試驗器材。侵徹試驗所需的設備器材包括:

①56式自動步槍一支,7.62 mm制式標準步槍彈若干,彈丸初速720 m/s;

②高速攝像機(15 000 fps)一臺,紅外攝像機(50 Hz)一臺,計算機兩臺;

③低溫箱一臺、模具數套、靶架及試驗臺一套、石膏粉/染色石膏粉若干、標簽紙若干。

1.2 凍土靶體制作

如圖2所示,射擊試驗采用的凍土靶體由高低溫試驗箱制作,靶體模具選用高強度透明亞克力板材以方便拍攝侵徹過程,試樣模具尺寸為Φ15 cm×40 cm,模具質量為700 g,內容積為7 068.6 cm3。試樣參數如表1所示。

圖1 試驗所用56式半自動步槍及7.62 mm彈丸Fig.1 Type 56 semi-automatic rifle and 7.62 mm projectile used in the experiment

圖2 試樣制備Fig.2 Sample preparation

表1 試樣參數表Table 1 Parameters of the tested samples

1.3 實驗步驟

①開啟高速攝像機及紅外攝像機,將靶體安裝固定在靶架上,拍照記錄;

②射手進入射擊位置,其余試驗人員撤離至安全位置,做好隱蔽;

③射擊結束,采集試驗圖像與數據,拍照記錄試驗后冰體,保存并處理高速攝像機及紅外攝像機的拍攝圖像;

④清理試驗現場,準備下一次試驗;

⑤清理試驗場地,剩余彈藥入庫保存。

1.4 實驗結果分析

①侵徹0 ℃凍土。

彈丸侵徹0 ℃凍土的結果如圖3所示,侵徹深度為360 mm,彈道初始段直徑為20 mm,彈道平直未出現彎曲;在侵徹深度為130 mm處彈道直徑擴大,出現空腔,空腔直徑為49.5 mm,此時彈丸失穩出現翻轉。

圖3 彈丸侵徹常溫凍土靶板彈道Fig.3 Ballistics of projectile penetrating 0 ℃ frozen soil

②侵徹-10 ℃凍土。

彈丸侵徹-10 ℃凍土后的靶板情況如圖4 所示,著靶瞬間高速相機記錄的開坑情況如圖5所示。彈丸在凍土靶板表面開坑,凍土靶板表面土樣在開坑瞬間向外噴出,彈丸穿透凍土,土樣表面明顯崩裂,試樣前端被明顯破壞,后端未被明顯破壞。

圖4 侵徹后-10 ℃凍土靶板Fig.4 Permafrost target plate at -10 ℃ after penetration

圖5 開坑瞬間高速攝像Fig.5 High speed video recording at the moment of opening the pit

如圖5所示,著靶瞬間彈丸溫度高于靶板,圖像中顯示為凍土靶板前小塊光斑,隨后彈丸在靶板表面開坑,靶板前端面凍土被擊散飛出。

③侵徹-20 ℃凍土。

彈丸侵徹-20 ℃凍土的結果如圖6所示,入射處彈孔直徑為9 mm,侵徹深度為95 mm,彈道軌跡出現偏斜未出現空腔,到侵徹結束彈丸被甲與鋼芯分離,彈頭保持向下。

試驗結束后,彈體的破壞情況對比如圖7所示。彈體包括銅材料披甲以及鋼材料彈芯兩個部分,0 ℃下彈體1與-10 ℃下彈體2在侵徹后均保持彈形完整,-20 ℃下彈丸在侵徹結束后被甲與彈芯分離,被甲變形嚴重,彈芯和彈頭被磨平。

圖6 侵徹后-20 ℃凍土靶板Fig.6 Permafrost target plate at -20 ℃ after penetration

圖7 試驗后彈體Fig.7 Projectiles after test

2 數值計算

2.1 有限元計算模型

計算過程采用 Lagrange算法,采用六面體網格對模型進行劃分,并對彈體的彈頭部位進行一定程度的加密,彈丸質量為7.9 g,彈體有限元模型如圖8所示。凍土材料的侵徹毀傷數值計算模型如圖9所示。彈體整體與凍土材料之間采用面面侵蝕接觸算法 ERODING_SURFACE_TO_SURFACE。被甲與鋼芯之間采用自動面面接觸算法AUTOMATIC_ SURFACE_TO_SURFACE。冰體材料上添加邊界條件BOUNDARY_SPC_SET及BOUNDARY_NON_REFLECTING,計算所采用的單位制為cm-g-μs。

圖8 彈體有限元模型Fig.8 Finite element model of projectile

圖9 彈-靶有限元模型Fig.9 Finite element model of projectile-target

2.2 彈體與凍土材料屬性

數值模擬計算中使用的7.62 mm口徑步槍彈銅質被甲材料采用Johnson-Cook本構模型。對于大變形、高應變率下和高溫條件下的材料屬性,Johnson-Cook本構模型能夠進行有效的描述[7],其在大應變情況下的本構關系方程為:

(1)

Gruneisen狀態方程定義的壓縮材料壓力為:

(2)

式中:p為壓力,μ為體積變化率,C為νs-νp曲線的截距,S1,S2和S3為νs-νp曲線斜率系數,γ0為Grunrisen常數,a為對γ0常數的一階體積修正,E為內能。C,S1,S2,S3,γ0和a為常數,且均為狀態方程參數[8]。

彈體材料模型與狀態方程的參數參考文獻[9],凍土材料參考文獻[10],選用Platic-Kinematic本構模型對凍土材料進行數值模擬。

2.3 計算模型驗證

圖10和圖11分別為彈丸侵徹0 ℃和-10 ℃靶板試驗數值計算結果,圖12為彈丸侵徹-20 ℃下靶板結束后試驗與計算剩余彈丸對比圖。

圖10 彈丸侵徹0 ℃靶板試驗與計算彈道對比Fig.10 Comparison between test and computational trajectory of projectile penetration at 0 ℃

圖11 彈丸侵徹-10 ℃下靶板試驗數值計算云圖Fig.11 Simulation result of projectile penetration at -10 ℃

圖12 侵徹-20 ℃下靶板后試驗與計算剩余彈丸對比Fig.12 Comparison between the test and simulation of remaining projectiles after penetration at -20 ℃

分別對0 ℃、-10 ℃和-20 ℃三組土樣靶板進行了試驗與計算對比驗證,計算出的侵徹深度、彈丸剩余質量以及侵徹結束后彈丸姿態與試驗結果對比如表2所示,數值計算的侵徹深度與彈丸剩余質量能夠很好地符合試驗結果,侵徹結束后彈丸姿態均出現翻轉,與試驗結果相同。對比結果表明模型較為可靠。

表2 試驗與計算對比驗證Table 2 Comparison and verification of experiments and calculations

3 溫度對侵徹凍土的影響

凍土材料具有明顯的溫度效應,隨著溫度的降低凍土材料中冰顆粒的強度以及冰與凍土之間膠結強度增大,凍土強度隨之增強。為研究溫度對于凍土材料抗侵徹性能影響,選取0 ℃、-5 ℃、-10 ℃、-15 ℃、-20 ℃,5個代表溫度進行有限元仿真計算。

彈丸侵徹姿態均為垂直侵徹,侵徹速度為700 m/s,含水率與侵徹試驗相同,為19%。由圖13可知,彈丸侵徹0 ℃凍土材料時,彈體姿態出現偏轉,彈丸失穩;在侵徹其他溫度靶板時彈丸均保持垂直侵徹狀態。0 ℃下,彈丸未發生大角度偏轉,彈孔直徑略大于彈徑;-5 ℃下彈丸侵徹凍土靶板過程中隨著彈丸偏轉,彈孔直徑逐漸增大,彈道軌跡未出現大角度偏轉;-10 ℃、-15 ℃、-20 ℃條件下彈丸未出現偏轉。

圖13 不同溫度下彈丸侵徹彈道軌跡Fig.13 Trajectory of projectile penetration at different temperatures

圖14為溫度對彈丸侵徹深度時程和加速度時程。由圖可知,隨著凍土溫度的降低侵深逐漸減小,彈丸侵徹凍土最大深度由0 ℃時287 mm降至-20 ℃時的103 mm。侵徹初始階段彈丸加速度較大,隨著速度的降低,加速度減小。對侵徹深度H與凍土溫度TP的關系進行擬合,結果可得H=9.1TP+285.6。計算結果與擬合公式對比如圖15所示,侵徹深度和凍土溫度基本呈線性關系,擬合公式的擬合效果較好。

圖14 不同溫度下侵徹深度和加速度時程曲線Fig.14 Time history curves of penetration depth and acceleration at different temperatures

圖15 凍土溫度對侵深的影響Fig.15 Influence of permafrost temperature on depth invasion

4 含水率對侵徹凍土的影響

含水率η是凍土材料力學性能的重要影響因素,假設彈丸以速度700 m/s,對-10 ℃下不同含水率的凍土靶板進行垂直侵徹,計算結果如表3所示。

表3 不同含水率下彈丸侵徹凍土結果Table 3 Penetration results under different moisture cuts

凍土含水率為15%時,彈丸侵徹深度最大335.1 mm,隨著含水率增大,侵徹深度減小,到22%含水率時侵徹深度最小為150.9 mm,之后彈丸侵徹深度隨著含水率的增大而增大。彈丸所受加速度的變化規律與侵徹深度一致,當侵徹含水率為22%的靶板時,彈丸受到的加速度為2.65×106m/s2。彈丸在4種含水率下均為正姿態侵徹,彈丸偏轉角較小。

結果表明彈丸在侵徹不同含水率凍土靶板時,侵徹深度隨著凍土含水率的上升呈先降低后升高的變化趨勢。在本文工況下,含水率為15%時,侵徹深度最深。在達到最佳含水率之前,彈丸侵徹深度隨著含水率增大而減小,凍土靶板超過最佳含水率后,彈丸侵徹深度隨著含水率增大而增大,侵徹深度與含水率擬合關系為:H=2.0η2-95.8η+1 309.8,式中η為含水率值。計算結果和擬合值的對比如圖16所示。

圖16 含水率對侵徹深度的影響Fig.16 Effect of moisture content on penetration depth

5 結論

①對0 ℃、-10 ℃、-20 ℃下土壤進行實彈射擊試驗,試驗結果表明,彈丸侵徹深度隨著溫度降低而減小。彈丸在侵徹常溫靶板過程中出現翻轉,彈道可分為三個階段,第一階段彈丸保持旋轉穩定,彈道直徑略大于彈徑,彈道未出現偏轉;第二階段彈丸失穩翻轉產生空腔,彈道軌跡出現偏轉;第三階段彈丸翻轉結束,以任意姿態向下侵徹至動能衰減為0。

②彈丸侵徹不同溫度下凍土靶板,侵徹深度與靶板溫度成線性關系:H=9.1TP+285.6,0 ℃以上彈道軌跡出現空腔且彈丸出現翻轉,0 ℃以下彈道軌跡未出現空腔且彈丸偏轉角較小,隨溫度降低彈道軌跡偏轉減小。

③含水率對凍土靶板防護能力有顯著影響,含水率對彈丸侵徹深度影響可擬合為H=2.0η2-95.8η+1 309.8。