摘要：為了研究多層混凝土介質內爆炸相似律,基於量綱分析法建立了多層混凝土介質內爆炸相似模型。利用有限元分析軟體LS-DYNA分析了裝藥深度和縮尺比例對原型和縮比模型毀傷形態及毀傷效應參量相似性的影響,並以梯恩梯和含鋁炸藥進行原型和縮尺比例為1/2的縮比模型內爆炸燬傷相似性試驗。結果表明:不同裝藥深度下的原型與縮比模型的毀傷形態相同,毀傷效應參量符合相似規律;靶板邊界約束相同時,含鋁炸藥在原型和縮尺比例為1/2的縮比模型內爆炸燬傷具有相似性,縮比毀傷效應參量歸一化後與原型的偏差在8%以內;隨著縮尺比例的減小,毀傷效應參量與原型的偏離程度增大,可選用適當縮尺比例的縮比模型試驗結果預測原型的內爆炸燬傷效應。
關鍵詞： 內爆炸 內爆炸燬傷效應 多層混凝土介質 爆炸力學 相似定律 縮比試驗
反堅固目標彈藥對機場跑道、港口及高速公路等目標的內爆炸燬傷效應及威力評估,是兵器領域的研究熱點。上述目標通常是由面層高強度混凝土、貧混凝土基層及壓實土層組成的多層混凝土介質。在目前彈藥設計及威力評估過程中,試驗研究仍是不可或缺的重要方法。以相似理論為基礎的縮比模型試驗能節約成本,縮短試驗週期,成為研究爆炸燬傷效應的一種重要手段。因此,分析多層混凝土介質內爆炸相似律,建立原型與縮比模型的相似準則,使模型試驗能夠預測原型試驗的毀傷效果,具有重要現實意義。
關於混凝土或混凝土-土介質在爆炸載荷作用下的成坑相似律,部分學者透過相似理論建立了相似準則。Westine[1]分析了混凝土-土複合介質內爆炸燬傷效應的關鍵影響因素,基於量綱分析法建立了相似模型。Chabai[2]在深埋裝藥的爆炸成坑效應中考慮了重力影響,得到4次方根比例定律,即r∝Q1/4(Q為炸藥釋放的總能量,r為相似比)。楊振聲[3]和談慶明[4]分析了工程爆破模型律,研究表明當介質拋擲高度y1與最小抵抗線W為同一量級時,即最小抵抗線超過25m時需考慮重力影響。Wang等[5]分析了雙層混凝土介質中的爆炸相似律,得到考慮重力時模型律的兩種方法:慣性加速度法(採用離心機實現)和等效材料法(材料強度進行縮比)。嶽松林等[6]推導了岩石內爆炸成坑相似律,闡述了地下深埋爆炸效應的模型試驗方法,隨後徐小輝等[7]給出了模型試驗等效材料的相似條件。楊亞東等[8]和Yao等[9]分析了縮比係數對密閉空間內爆炸相似律的影響。Xie等[10]基於相似理論建立了土壤爆破彈坑引數的無量綱經驗公式。荊松吉等[11]採用有限元分析軟體LS-DYNA分析了多層混凝土介質的爆炸相似律,分析表明多層混凝土介質中存在相似律,但縮尺比例過大時,會透過面效應影響原型與縮比模型的相似性。
上述文獻主要是基於相似理論分析混凝土-土複合介質結構內爆炸相似律,對多層混凝土介質內爆炸相似律的研究較少。基於量綱分析法和數值模擬的混凝土內爆炸相似律分析都是將混凝土看成均質材料,需要透過考慮骨料的內爆炸試驗對相似律進行驗證。此外,由於含鋁炸藥具有高爆熱、高爆容和釋能時間長的優點,廣泛應用於航空武器彈藥,而含鋁炸藥中鋁粉與爆轟產物的二次反應表現出非理想性,在多層混凝土介質內爆炸燬傷是否滿足相似律,需要相關研究驗證。
本文利用量綱理論匯出多層混凝土介質內爆炸燬傷相似律,建立內爆炸相似理論模型。在理論分析基礎上,利用LS-DYNA軟體進行了不同裝藥深度和縮尺比例下原型和縮比模型的內爆炸數值模擬,並採用梯恩梯(TNT)和含鋁炸藥進行原型和縮尺比例為1/2的縮比模型內爆炸相似性驗證試驗。透過模擬和試驗分析了裝藥深度及鋁粉的二次反應釋能對毀傷效應參量、介質拋擲速度和位移相似性的影響,揭示了縮尺比例對內爆炸燬傷效應參量的影響,論證了模型試驗的可行性。所得結論可為多層混凝土介質內爆炸縮比模型試驗設計提供參考依據。
1、多層混凝土介質內爆炸相似模型
1.1毀傷效應參量及影響因素
裝藥在多層混凝土介質內爆炸燬傷效應參量主要有:實際彈坑深度H(實際彈坑底部到地表的距離);彈坑半徑Rc(混凝土道面上彈坑的平均半徑);有效毀傷半徑Reff(混凝土道面上環向裂紋的平均半徑);最大空腔半徑Rmax(在成坑模式中,Rmax一般等於Rc);鼓包隆起過程中的鼓包位移y和鼓包速度v.此外,還引入1個參量MODE來描述毀傷模式(成坑模式、隆起模式或隱坑模式)。
上述毀傷參量的主要影響因素來自裝藥、混凝土、土壤的材料性質及幾何尺寸等,包括:
1)裝藥:裝藥量w,裝藥密度ρe,裝藥爆熱Qe和爆轟產物膨脹指數γ;
2)多層混凝土介質靶:面層混凝土參數包括面層厚度hf、面層長度、寬度lf、反映慣性效應的材料密度ρf、反應強度效應的材料極限抗壓強度σf;碎石層引數包括碎石層厚度hm、介質材料密度ρm、極限抗壓強度σm;壓實土壤層包括介質密度ρs、聲速cs.聲速cs和密度ρs反映材料強度效應的影響。忽略多層混凝土介質材料的應變率效應,認為其他強度引數與抗壓強度呈正比;
3)交會引數:裝藥深度h和裝藥傾斜角度α.
1.2內爆炸燬傷量綱分析
在進行量綱分析前作如下假設[12]:1)不考慮介質應變率的影響,忽略介質的黏性和熱傳導;2)混凝土為各向同性的脆性材料;3)忽略爆轟產物初始壓力的洩露;4)由於炸藥裝藥深度較小,忽略重力影響;5)炸藥瞬時爆轟。表1和表2分別列出了7個毀傷效應參量和15個影響因素及其量綱(以M、L和T為基本量綱)。
表1毀傷效應參量的量綱
表2影響因素及其量綱
選取w、Qe和h為基本量,可將7個毀傷效應參量寫成無量綱形式:
根據!定理,以上公式可寫成如下無量綱關係式:
根據量綱分析相似定理,當原型和縮比模型(1)式中的!項對應相等,即相似條件為(2)式時,則πip=πim,下角標p表示原型,m表示模型。
當縮比試驗和原型試驗中採用的炸藥型別、密度及靶板材料、強度和密度相同,以及縮尺比例為λ時,求解(2)式可得相似模型如下:
由πip=πim可得原型和縮比模型的毀傷效應參量關係,如表3所示。
表3原型和縮比模型毀傷效應參量關係
根據表3中原型和縮比模型體系內的鼓包位移及速度關係,可知時間滿足tm/tp=λ.下面透過數值模擬檢驗相似模型的正確性,並分析裝藥深度和縮尺比例對相似性的影響。
2、內爆相似性數值模擬
2.1計算模型及材料引數
多層混凝土介質由40cm厚混凝土、40cm厚碎石層和120cm厚壓實土層組成,如圖1所示。其中混凝土面層由邊長為5m的正方形中心靶及寬度為1.5m的邊界靶組成,中心靶和邊界靶存在1cm的間隙。裝藥為TNT炸藥,位於靶板中心位置,並用壓實土填充炮孔。考慮到結構的對稱性,建立1/4模型,數值模擬採用多物質Euler材料與Lagrange結構相耦合的演算法,其中面層混凝土採用單點積分Lagrange六面體網格,其餘均採用Euler六面體網格,原型網格尺寸為3cm,縮比模型網格尺寸根據縮尺比例縮比,即縮尺比例為1/2的縮比模型網格尺寸為1.5cm.
TNT炸藥採用High＿Explosive＿Burn材料模型及JWL狀態方程描述。
圖11/4計算模型
式中:p為爆轟產物壓力;A、B、R1、R2和ω為表徵炸藥爆轟特徵的常數;E0為初始比內能;為相對比容。TNT炸藥材料引數如表4所示。
表4TNT炸藥材料引數
混凝土採用TCK模型,其主要特點是把動態脆性破裂過程看作拉伸損傷的累積,材料的壓縮響應看作理想彈塑性。透過使用者自定義模型材料介面將TCK模型嵌入LS-DYNA中,混凝土材料模型引數如表5所示。表5中:Ef為楊氏模量,Yf為屈服強度,ν為泊松比,kic為材料斷裂韌性,k為材料常數,Fs為拉伸失效應變,Fd為壓縮失效應變。
表5混凝土材料引數
碎石層採用流體彈塑性模型和Gruneisen狀態方程,其中ρm=2.2g/cm3,Ym=8MPa.壓實土層採用帶失效的土壤與泡沫材料模型,具體引數見文獻[13]。
2.2裝藥深度對相似性影響分析
本節建立原型(縮尺比例為1)和縮尺比例為1/2的縮比模型內爆炸燬傷計算模型,分析不同裝藥深度下原型與縮比模型毀傷效應參量(彈坑口部半徑、有效毀傷半徑、最大爆腔半徑和彈坑深度)的相似性。其中原型裝藥量為3kgTNT,裝藥深度分別為40cm、60cm和90cm,靶板尺寸如2.1節所示。縮比模型的長度尺寸減小一半,即面層和碎石層為20cm,壓實土層為60cm,裝藥量為0.375kg,裝藥深度分別為20cm、30cm和45cm,所有材料模型及引數保持不變。
2.2.1毀傷形態對比
3種裝藥深度下的原型和縮比模型毀傷效果如圖2所示,其中原型為40ms時的毀傷效果,縮比模型為20ms時的毀傷效果,此時毀傷形態基本形成。由圖2對比可知:不同裝藥深度下原型和縮比模型的毀傷特性較為一致,原型和縮比模型的面層混凝土均出現了破碎和裂紋,碎石層和壓實土層形成的爆腔形態相同,由於土壤強度較低,壓實土層的空腔較碎石層大。在裝藥深度H分別為40、60cm時,原型和縮比模型的面層出現徑向裂紋,裂紋的走向和條數較為一致。由此可以看出,不同裝藥深度下原型和縮比模型的毀傷形態具有相似性。
2.2.2毀傷過程對比
以原型裝藥深度60cm為例,對比原型和縮比模型毀傷過程中介質運動的相似性。原型選取混凝土面層表面距圓心分別為50cm、100cm和150cm的3個節點,縮比模型距圓心為25cm、50cm和75cm3個節點,分別對比不同位置處的豎直方向節點速度和位移。圖3和圖4分別為節點位移和節點速度時程曲線,由於位移和速度均與時間有關,為了方便比較相似性,將縮比模型的位移和速度曲線化為原型時間體系下的時程曲線。從圖3中可以看出,原型和縮比模型相似位置處的節點位移較為一致,不同位置處的節點位移滿足相似模型中的yp/ym=2的相似關係。從圖4中可以看出,不同位置處原型和縮比模型的節點速度時程曲線趨勢相同,只是距離炸藥較近的節點在末端速度差別較大,原因是距離炸藥較近時,由於衝擊波作用,會產生層裂和拋擲現象,造成節點速度有所差別,但兩個體系下的速度近似滿足vp/vm=1的相似關係。因此原型和縮比模型的毀傷過程滿足相似性。
圖2原型和縮比模型的毀傷效果
圖3節點位移時程曲線
圖4節點速度時程曲線
2.2.3毀傷效應參量的對比
多層混凝土介質的毀傷包括炸坑、爆腔、裂紋及隆起,由於毀傷效果的複雜性,難以精確判斷毀傷半徑。本文的有效毀傷半徑是指面層混凝土嚴重毀壞區域(包括炸坑和環向裂紋內部的隆起,沒有考慮徑向裂紋的影響)半徑,不同裝藥深度下原型和縮比模型的彈坑口部半徑Rc、有效毀傷半徑Reff、最大爆腔半徑Rmax和彈坑深度H(兩次測量平均值)如表6所示,原型裝藥深度為90cm時,沒有形成漏斗坑,由於計算模型的壓實土層厚度不夠,無法測量彈坑深度。從表6中可以看出,原型和縮比模型的彈坑口部半徑、有效毀傷半徑和彈坑深度的相似比差別不大,均在1∶0.48到1∶0.52之間,而不同炸點深度下的爆腔相似比差別較大,主要原因是網格數量較少引起的介質流動介面的誤差。由此可見原型和縮比模型的毀傷效應參量滿足2∶1的幾何相似性。
表6不同裝藥深度下的毀傷效應參量
2.3縮尺比例對相似性影響分析
分別進行縮尺比例為1、1/2、1/4和1/6的縮比模型內爆炸燬傷模擬,分析不同縮尺比例下的內爆炸燬傷效應參量是否滿足相似規律。裝藥量分別為3kg、375g、46.88g和13.89g(原型的1/λ3),裝藥深度分別為60cm、30cm、15cm和10cm,靶板尺寸按縮尺比例縮小,材料模型和引數不變。其中縮尺比例λ=1和λ=1/2的毀傷效果如圖2(b)和圖2(e)所示,縮尺比例λ=1/4和λ=1/6的毀傷效果如圖5所示,不同縮尺比例下的毀傷效應參量如表7所示。
圖6為不同縮尺比例下的毀傷效應參量歸一化(除以縮尺比例)後與原型毀傷參量相似性對比。
從圖5中可以看出,原型和縮尺比例為1/2的模型裂紋擴充套件及毀傷區形態相似,進一步增大縮尺比例時,毀傷區形態差異逐漸增大,表面無徑向擴充套件裂紋。從表7和圖6可以看出,原型與不同縮尺比例的縮比模型歸一化後彈坑口部半徑比為1∶0.98∶0.96∶0.72,有效毀傷半徑比為1∶0.98∶0.96∶1.08,最大爆腔半徑比為1∶1.06∶1.12∶1.2,彈坑深度比為1∶1.02∶1.01∶1.05.隨著縮尺比例的增加,毀傷效應參量與原型的偏差逐漸增大,相似性呈減小趨勢,λ=1/6時的縮比模型歸一化後與原型的毀傷參量相對誤差最大為28%.這是因為:1)縮尺比例逐漸變小時,毀傷效應參量的取值偏差逐漸增大;2)多層混凝土介質內爆炸燬傷存在層裂、崩落和層間相互作用等面效應[6],從能量角度來看,當縮尺比例逐漸減小時,用於面效應的能量佔總能量的比例逐漸增加,對相似性的破壞越嚴重;3)原型和縮比模型的靶板材料和強度相同,導致斷裂韌性相同,不符合相似律,因此縮比模型對靶板的斷裂抗力比原型大。由此可見,多層混凝土介質內爆炸縮比模型試驗的縮尺比例不宜過小,一般不能小於1/4.
圖5毀傷效果圖
表7不同縮尺比例下的毀傷參量
圖6縮尺比例對毀傷效應參量的影響
3、內爆相似性試驗驗證
由於數值模擬中將混凝土視為均質材料,沒有考慮混凝土中骨料對相似律的影響,模擬也無法驗證含鋁炸藥在多層混凝土介質內爆炸相似性,因此構建原型和縮尺比例為1/2的縮比模型,分別採用TNT和海薩爾(RDX/Al/binder)炸藥進行內爆炸相似性驗證試驗。
3.1試驗概況
原型靶為邊長5m的正方形,由40cm厚混凝土層、40cm厚碎石層和半無限厚壓實土壤層組成,縮比模型尺寸根據原型按縮尺比例λ=1/2進行縮比,如圖7所示,各塊靶標間由1cm厚木板隔開,靶標中間預留裝藥孔。裝藥分別為TNT和海薩爾,原型裝藥量為3kg,其中TNT裝藥深度為53cm,海薩爾裝藥深度為51cm.縮比模型裝藥量為375g,裝藥深度按縮尺比例縮比,炸藥裝填及起爆方式如圖8所示。
圖7試驗靶標
圖8裝藥裝填方式
3.2鼓包隆起過程相似性分析
以海薩爾炸藥爆炸燬傷原型和縮比模型為例,兩個體系相似時間處的鼓包隆起過程如圖9所示。從圖9中可以看出,由於鋁粉燃燒,在縮比模型中可以看到火光從炮孔中噴出,原型和縮比模型的鼓包隆起過程相似,均先產生徑向裂紋,後產生環向裂紋,在爆生氣體作用下,徑向裂紋和環向裂紋逐漸擴充套件變寬,最後形成梯形混凝土碎塊。
圖9鼓包隆起過程(海薩爾炸藥)
透過分析高速攝影機獲得的各種工況鼓包運動過程,得到TNT和海薩爾炸藥爆炸燬傷縮比模型和原型在不同時刻處的鼓包隆起高度和速度,分別如圖10和圖11所示,其中原型λ=1,縮比模型λ=1/2.
圖10鼓包隆起高度-時間曲線
圖11鼓包隆起平均速度-時間曲線
從圖10中可以看出,TNT和海薩爾炸藥在相似時間處的鼓包隆起高度近似滿足Hp/Hm=2,其中TNT的縮比模型隆起高度除以縮尺比例後與原型的最大誤差為8%,海薩爾炸藥的最大誤差為6%.從圖11中可以看出,原型和模型體系中的速度近似滿足vp/vm=1,其中TNT體系速度偏差較大的原因是因為縮比模型中有兩邊界約束較弱,應力波達到邊界後沒有反射,造成能量外洩,導致縮比模型的鼓包平均速度與原型相比較小。由此可以看出,在邊界約束相同時,原型與縮比模型的隆起過程符合相似規律。
3.3毀傷效應參量相似性分析
TNT和海薩爾炸藥的爆炸燬傷效果如圖12所示。從圖12中可以看出,由於TNT爆炸燬傷試驗的縮比模型靶標有兩邊界約束較弱,造成部分割槽域隆起高度較大,但縮比模型和原型的裂紋條數均為8條,擴充套件方向也大致相同,只是縮比模型由於邊界效應導致裂紋寬度較大。海薩爾炸藥的原型和縮比模型毀傷效果大致相同,均為成坑模式,徑向擴充套件裂紋均為9條,裂紋擴充套件方向也相同。由此可以看出,邊界約束相同時,原型和縮比模型的毀傷模式及裂紋擴充套件形態吻合較好。
圖12爆炸燬傷效果圖
TNT和海薩爾炸藥毀傷原型和縮比模型後的毀傷效應參量如表8所示。
表8毀傷效應參量的對比
注:Dc為彈坑口部直徑,Deff為環向裂紋直徑。
透過對比可知,TNT炸藥毀傷的原型與縮比模型彈坑口部直徑比約為1∶0.35,環形裂紋直徑比為1∶0.47,最大爆腔半徑比為1∶0.52,實際彈坑深度比為1∶0.47.除彈坑口部直徑外,縮比模型的毀傷參量(環向裂紋直徑、最大爆腔半徑和實際彈坑深度)除以縮尺比例後與原型的誤差均在6%以內。彈坑口部直徑偏差較大是因為靶板邊界效應的影響,縮比模型在衝擊波和爆轟氣體的作用下,部分靶塊隆起高度較大,導致面層混凝土崩落較少,彈坑口部直徑較小,使彈坑口部的相似比偏差較大。
海薩爾炸藥毀傷的原型與縮比模型毀傷模式相同,且原型和縮比模型邊界約束相同,其彈坑口部直徑比為1∶0.52,均無環向裂紋,最大爆腔半徑比為1∶0.70,實際彈坑深度比為1∶0.54.除最大爆腔半徑外,縮比模型的毀傷參量(彈坑口部直徑和實際彈坑深度)除以縮尺比例後與原型的誤差均在8%以內,最大爆腔半徑相似比的誤差較大的原因可能是因為混凝土層和碎石層的強度相差較大,在爆生氣體的作用下,兩層介質間出現了分層,使試驗測量的爆腔半徑包含層間間隙的長度,導致縮比模型的最大爆腔半徑偏大。另外,由於靶板養護的差別,試驗中材料的力學效能也不可能完全相同,會存在一定離散,造成試驗結果存在一定誤差。因此在邊界約束相同條件下,兩種炸藥的原型和縮比模型毀傷效應參量滿足幾何相似性。
4、結論
本文透過建立多層混凝土介質內爆炸燬傷效應相似模型,並對內爆炸燬傷效應參量和毀傷形態的相似性進行數值模擬和試驗研究。得出結論如下:
1)適當的縮尺比例下,不同裝藥深度下原型和縮比模型的毀傷形態相同,毀傷效應參量(彈坑口部半徑、有效毀傷半徑、最大爆腔半徑和彈坑深度)滿足幾何相似性,原型和縮比模型的鼓包隆起高度和速度服從相似規律。
2)在邊界約束條件相同時,含鋁炸藥在原型和縮比模型(λ=1/2)內爆炸燬傷破壞現象相同,毀傷效應參量符合相似性,除最大爆腔半徑外,毀傷特徵參量歸一化後與原型的誤差均在8%以內。
3)隨著縮尺比例的減小,縮比模型的毀傷形態與原型差異增大,毀傷效應參量與原型的偏離程度也逐漸增大,相似程度呈下降趨勢。工程上可採用縮尺比例不小於1/4的縮比模型來預測原型的內爆炸燬傷效應。
