■ 胡智航 成遠礦業(yè)開發(fā)股份有限公司
■ 劉萬通 遼寧科技大學
■ 孟強 白爽 何茂林 成遠礦業(yè)開發(fā)股份有限公司
Luccioni B[1]和崔溦等[2]通過模擬對爆炸成坑進行了數值分析,土體密度對爆坑直徑影響較大,剪切模量影響較小。穆朝民等[3]研究炸藥在不同深度下爆炸成坑和應力波傳播規(guī)律,得出半封閉狀態(tài)下爆坑半徑經驗公式。賈永勝等[4]和范一鍇等[5]研究了藥量、深度、距離對爆坑成型規(guī)律。徐干成等[6]和岳松林等[7]通過模擬得出爆炸成坑相似規(guī)律。張宇等[8]利用Autodyn模擬爆炸成壕構建單兵掩體和實物實驗做對比,證明數值模擬的可靠性。景彤等[9]利用ANSYS/LS-DYNA 有限元軟件模擬低速侵徹土壤,在土質參數準確的情況下數值模擬具有一定可靠性
本文主要研究無人機投彈重力侵徹土壤,土中爆炸拋撒土顆粒滅火與爆炸成坑形成隔離帶阻隔滅火研究。
本文數值模擬采用SPH 算法,一種特殊的拉格朗日數值計算方法,離散化過程不需要網格劃分,進而不會出現網格纏結和網格單元畸變影響計算結果。SPH 算法最早應用于天體物理研究,后廣泛應用于大型變問題、沖擊破壞等。
SPH 算法通過質點的積分方程進行求解,同樣需要守恒方程與材料本構方程。該算法將特定空間離散成一系列的粒子,粒子之間沒有固定連接,每個粒子有質量、內能、速度、坐標等。采用SPH 方法離散的質量、動量和能量守恒控制方程如式下[10]:
式中:P 為壓力;V 為相對比容;E 為單位體積的內能。彈體材料采用STEEL 4340,用 Shock 狀態(tài)方程和 Johnson Cook 強度模型。參數見表1、表2。
表1 4340 鋼參數
表2 TNT 的JWL 狀態(tài)方程參數
投彈某靶場的土質為松軟的黃土,數值計算采用Autodyn 中的Sand[11]模型,選用線性壓實狀態(tài)方程,強度模型采用顆粒強度模型,該模型考慮了不同土介質的強度,同時可以描繪密度強化效應、剪切模量與密度的關系,該模型可以有效地描述爆坑壁密度變化[12]。
彈體侵徹土壤經驗公式別列贊公式:
由上述公式可知彈體侵徹土壤主要受到土壤介質的阻力作用,阻力的大小于彈體口徑、速度、結構、侵徹介質性質有關。
建立侵徹初速度為20m/s、動能相同的三種重力侵入模型,分析彈頭形狀侵入土體作用過程,以及彈頭形狀對重力侵入土體深度的影響規(guī)律.
三種不同形狀彈體依靠重力侵入土體作用過程圖,彈體以20m/s 的速度垂直侵徹土壤介質,在彈體侵入過程中,靠近彈頭的土體介質受到很強的沖擊,產生很大的徑向應力和切向應力,其值遠遠大于土壤的動態(tài)抗拉強度和抗壓強度,彈體周圍土體受到強烈壓縮;隨著彈體繼續(xù)侵入,周圍土壤持續(xù)向外部擠壓作用,土體被不斷壓密,侵徹耗能增加,故而在后期彈體動能急劇減少,侵徹速率降低。彈頭形狀的侵徹深度數據見表3。
表3 彈頭形狀的侵徹深度數據表
半球形彈頭的侵徹深度相比于卵形彈頭和圓錐形彈頭侵徹深度較淺,主要是由于半球形彈頭在最初的接觸面積較大,侵徹過程受力面積較大,土壤壓實面積較大,單位時間內耗能過快,導致侵徹深度較小。圓錐彈頭與卵形彈頭在侵徹過程中的受力較小,在實際實驗中可將彈頭設計為圓錐形彈頭或者卵形彈頭,重力侵入效果會更好。
投彈高度影響彈體侵徹土壤的初始速度,該速度直接影響侵徹深度,隨著侵徹深度增加,爆坑體積也在增加,但是達到有利深度后爆坑體積會減少。民用無人機飛行高度會有一定限制,研究投彈高度對侵徹影響,可以為后期實物實驗提供指導。以卵形彈頭作為侵徹模型,觀察彈體重力侵徹土壤,分析投彈高度對重力侵入深度的影響規(guī)律,為了簡化計算過程,忽略空氣阻力,彈體按自由落體運動計算接觸初始速度,具體工況安排如表4。
表4 投彈高度下侵徹初速度數據表
圖1 速度位移變化時程曲線
深度依次間隔0.5m,分析一定藥量下侵入深度對爆炸鼓包、拋擲、成坑及覆土范圍的影響。
根據實驗現象可以觀測到隨著侵徹深度的增加,成坑深度呈上升趨勢,但上升速率逐漸減小,下部土體受到爆炸作用逐漸密實,導致成坑深度的增加速率逐漸減小。成坑開口和覆土范圍隨著深度的增加呈現上升后下降的趨勢,侵徹深度較小時,炸藥能量部分作用于空氣,范圍較小;侵徹深度逐漸增加,爆炸拋擲土體增加,范圍逐漸增大;侵徹深度的進一步增大,炸藥爆炸導致土體鼓包運動,拋擲較少土體,導致兩者范圍下降。
綜上所述,彈頭形狀與投彈高度會影響直接影響侵徹深度,裝藥量一定條件下侵徹深度會影響爆坑直徑與爆炸覆土范圍。覆土范圍與覆土量會影響后期滅火效果,爆坑同時可以充當隔離溝,阻隔火體蔓延。
為研究爆炸覆土滅火效果,在某場地進行實物實驗,土壤以碎石土壤為主,人為制造火場,創(chuàng)造侵徹深度,分別設計兩種不同工況。彈頭形狀的侵徹深度數據見表5。
表5 彈頭形狀的侵徹深度數據表
引爆后巖土被拋撒而出,土顆粒經過拋撒最終周向覆蓋火圈,覆土面積大于火圈面積,火源半徑1m、2m 的火場均被撲滅,且無復燃現象。爆炸覆土可以有效撲覆土范圍內的火源,靜爆實驗證明爆炸覆土滅火的可行性。
實物試驗微觀現象無法觀測,在相同工況下進行模擬。(見圖2)模擬仿真同樣使用Autodyn 建立10m9m 土域,填充粒子數為800000,炸藥尺寸、位置與實物試驗一致。
圖2 觀測點密度、速度變化曲線
爆轟產物與爆生氣體迅速對爆坑底部土壤壓縮,炸藥附近的土壤密度與速度迅速增加到最值,爆坑底部土壤不斷被壓實,后期坑底土壤密度變化很?。浑S著爆炸載荷施加于爆坑側壁,側壁土壤同樣被迅速壓縮,速度、密度不斷增大,后期兩側土壤部分被拋撒而出,密度開始減少,被拋撒而出的土顆粒速度短時間保持一定速度周向飛散,由于周向的土壤阻擋,未被拋撒出的土壤速度不斷減少,最終兩壁側土壤速度趨近于0,表明土體不斷被壓縮,直至被壓實。
圖3 壓力云圖與密度云圖
根據實驗現象可以觀測到爆坑深度與爆坑直徑不斷增加,爆坑深度在14ms 時增長趨于平緩,最終穩(wěn)定于1m;爆坑直徑在17ms 時增長趨于平緩,最終穩(wěn)定于2.25m。開始時大部分爆炸載荷用于壓縮坑底,爆坑深度增加明顯,隨著時間的增加,爆炸載荷開始用于拋撒與壓縮壁側土壤,爆坑直徑增加明顯,爆坑深度趨于穩(wěn)定。
實物試驗(來源于山東特種工業(yè)集團)在某靶場進行,土質為松軟沙土,裝藥為4.3kg 的TNT,裝藥尺寸為Φ0.09m×0.41m,彈體侵徹土壤0.5m 后引爆。最終形成半徑為1.2m,深度為0.9m,覆土范圍約3.2m。(見圖4)
圖4 實物投彈爆炸覆土
如表6 所示,數值計算得到的爆坑直徑與爆坑深度與實物試驗比較一致,最大誤差為10%,表明數值建模方法與計算方法合理可靠。
表6 計算值與實際值對比
建模與上述一致,裝藥距離為1m,裝藥情況與觀測點設置如圖5,觀測點間隔為10cm。30ms 后形成爆坑直徑為3.3m,爆坑深度為0.82m,組合投彈可以有效增大爆坑直徑。應用于森林火災可以在迅速形成防火溝(阻火隔離帶)。
圖5 組合投彈模擬
1. 彈體形狀、投彈高度會影響侵徹深度,侵徹深度會影響爆坑尺寸與爆炸覆土范圍。
2.SPH 算法可以形象描述爆炸拋撒土體顆粒飛散,爆坑內土體密度。數值模擬仿真與實物試驗誤差較小,比較合理。
3. 組合投彈不僅可以擴大爆坑尺寸,同時可以迅速形成阻火隔離帶。