基于AUTODYN的壓縮空氣彈射內彈道研究

摘 要：以超近程防御武器系統壓縮空氣發射裝置為背景，利用流固耦合軟件AUYTODYN建立了不同空氣壓力、不同泄流面積、不同內彈道長度下的壓縮空氣彈射內彈道模型，通過對模型流固耦合數值仿真結果的分析，得到了彈藥速度與時間、彈藥速度與位移、彈藥速度與空氣壓力、彈藥加速度與時間的關系，并在此基礎上確定了發射裝置所用的壓縮氣體的體積與壓力以及發射管的內彈道長度。

關鍵詞：內彈道；壓縮空氣；彈射；流固耦合；AUTODY軟件

中圖分類號：TJ768 文獻標識碼：A 文章編號：1673-5048（2014）05-0046-04

0 引 言

超近程防御武器系統是一種用于打擊距離被保護目標20～200m范圍內的末敏彈、子母彈等來襲彈藥武器系統。其采用的壓縮空氣彈射方式，避免了采用自力發射方式時彈藥推進劑燃燒引起燃氣射流沖擊問題，既不會產生較大的紅外特征，又提高了彈藥的初速與有效載荷，同時還滿足超近程防御作戰紅外特征小、反應速度快的要求[1]。

壓縮空氣彈射系統最早應用于魚雷的發射，二戰后隨著彈道導彈和超視距空空導彈的發展，逐漸應用于戰略導彈和中距空空導彈的發射[2]。目前，人們對壓縮空氣彈射系統的內彈道研究已經取得了大量的研究成果，喬汝椿提出了采用壓縮空氣發射方式的魚雷發射裝置管內魚雷運動參數的估算方法[3]，許斌采用SIMULINK和ADAMS的聯合仿真，研究了機載導彈彈射式發射動力學行為[4]，廖振強建立了高壓氮氣驅動下彈射裝置的數學仿真模型，分析了導彈的分離參數和分離姿態[5]。

本文在前人研究的基礎上，以用于超近程防御系統的壓縮空氣發射裝置為應用背景，在炮口初速要求120m/s的條件下，利用AUTODYN軟件，建立壓縮空氣發射裝置內彈道模型，根據模型的數值仿真研究壓縮空氣壓力、發射裝置部分結構參數與彈丸彈射后初速的關系，為發射裝置進行結構設計提供依據。

1 發射裝置結構與內彈道模型

參考其他裝置，壓縮空氣發射裝置結構設計如圖1所示，該裝置采用電磁鐵閉鎖機構控制發射，其原理是：彈藥在閉鎖體的作用下緊貼壓縮氣瓶泄流口處的薄膜膜片，以此控制高壓氣體。在發射時，電磁鐵通電，銜鐵向左移動，手柄和閉鎖體失去銜鐵約束，向上運動，此時閉鎖體無法繼續約束彈藥，壓縮空氣沖破膜片，進入發射管的內彈道部分，推動彈藥運動。

考慮到平時維護以及重復使用的需要，將方案的使用過程設計為：發射管后部即為存儲壓縮空氣的空間，在平時不儲存壓縮空氣，進入戰備狀態時，先在泄流口處貼上薄膜膜片，然后裝填彈藥，使彈藥緊貼薄膜膜片，壓下手柄使閉鎖體緊貼彈藥，電磁鐵銜鐵壓住手柄，閉鎖力開始作用，彈藥受到約束。此時，開始對發射管后部壓縮空氣存儲區域充氣，充氣完畢后，系統進入待機狀態。

1.1 發射裝置內彈道模型相關參數的確定

在上述裝置基礎上，已知彈藥質量：4kg；炮口初速：120m/s；能量利用系數η：0.15；定容比熱容：717J/（kg·K）；發射閥節流后的空氣溫度：266K。則發射該彈藥所需的壓縮空氣質量按下式計算：

mB0=mTυ2Tc 100.1mm，則罐體長度L=78mm。

壓縮氣瓶尺寸數量級和大概尺寸范圍由此確定，該數據可作為AUTODYN中進行彈道估算的初始數據[6]。

1.2 分析內彈道性能影響的模型

為了驗證高壓空氣壓力對內彈道性能的影響，在AUTODYN軟件中建立模型1以驗證不同壓力的高壓空氣對炮口初速等內彈道性能的影響[7]。為了驗證高壓氣瓶泄流面積對內彈道性能影響，在AUTODYN軟件中建立模型2以驗證不同泄流面積下的高壓空氣對炮口初速等內彈道性能的影響。

模型1的特點是考慮因素較多，除了要考慮在模型中填充空氣的壓縮，實際情況下高壓空氣由發射氣瓶通過口徑較小的通道流入發射管內彈道部分時的壓縮，還要考慮發射管與彈藥之間的摩擦。而模型2則主要考慮泄漏，相較簡單。

在模型簡化中，為了便于計算，將彈藥模型簡化為圓柱形，彈藥簡化后的狀況如圖2所示。=100.1mm。

不同的是在模型1的高壓空氣區域中在泄流直徑為20mm基礎上，分別填充進34MPa，35MPa，36 MPa的高壓空氣。而模型2在高壓空氣為35MPa基礎上，泄流直徑分別為20mm，24mm和28mm。

在彈管間隙問題上，由于壓縮空氣發射裝置和火炮、火箭炮等發射系統不同，在野戰條件下不存在由火箭發動機燃燒室高溫引起的彈體徑向彈性變形、炮管與彈藥之間存在溫度差、連續射擊過程中產生積碳等問題，因此能明顯影響壓縮空氣發射裝置彈管間隙的因素是彈藥幾何形狀誤差。彈藥幾何形狀誤差，主要受制于構成彈藥的各個部件的連接螺紋同心度以及彈藥各個構件的軸線彎曲[8]。并且在彈藥的各個構件中，戰斗部往往壁厚較大，剛度較大，不易產生幾何彎曲，而超近程防御武器系統采用壓縮空氣發射方式，彈藥上沒有火箭發動機或者其他薄壁結構，所以用于超近程防御武器系統的彈藥往往剛度較大，其彈管間隙的設定可明顯小于火箭炮等類武器的彈管間隙。通過比較相關武器的彈管間隙，將模型彈管間隙設定為0.1mm。

1.3 流固耦合模型的建立

由于涉及到流體與固體的相互作用，因此涉及到流固耦合，可以在AUTODYN軟件中實現這一過程。模型的建立過程：建立六個部件模型，分別是發射管模型、彈藥模型、發射管上邊界模型、發射管下邊界模型、泄流通道上邊界模型、泄流通道下邊界模型。建立的模型如圖3所示。

發射管模型由Euler單元建立，填充進空氣，在高壓空氣氣瓶區域填充進高壓空氣，在炮口設置出流邊界。彈藥模型由Lagrange單元建立，填充進與彈藥密度相同的材料。發射管上、下邊界模型由Lagrange單元建立，放置在發射管上邊界，在仿真過程中用于約束彈藥由于彈管間隙在+Y方向的運動。

泄流通道上邊界模型和泄流通道下邊界模型均由Lagrange單元建立，放置在模擬的泄流通道的兩邊，在仿真過程中用于約束空氣流動，模擬出泄流通道的效果。

在AUTODYN軟件分析中的可用耦合分為Lagrange-Lagrange耦合和Euler-Lagrange耦合兩種。在仿真模型中，Lagrange-Lagrange耦合主要用于彈藥和發射管邊界的接觸和摩擦分析，Euler-Lagrange耦合主要用于高壓空氣推動彈藥運動的流固耦合分析[9]。

在Lagrange-Lagrange耦合的設定中，設置彈藥模型與上下發射管模型的摩擦系數為f=0.15，以模擬彈藥在發射管中所受的來自發射管壁面的摩擦力。

Lagrange-Lagrange耦合要求兩個零件間的最小距離不得小于最小網格尺寸的1/10，在本模型中，彈管間隙為0.1mm，單邊距離為0.05mm，所以剖分網格時須將最小網格尺寸設定為小于0.5mm。

在Euler-Lagrange耦合類型的設定中，Euler-Lagrange耦合設定為全自動耦合[10]。

2 數值仿真結果及分析

2.1 模型1的仿真結果

針對不同壓力的高壓空氣對炮口初速等內彈道性能的影響，得到模型1的速度—時間以及速度-位移曲線如圖4所示。

在泄流直徑20mm的條件下，34MPa，35MPa壓力的壓縮空氣加速過程不平穩，炮口初速無法滿足120m/s炮口初速的要求，36MPa壓力的壓縮空氣加速過程平穩，可使彈藥炮口初速達到116m/s。

2.2 模型2的仿真結果

針對不同高壓氣瓶泄流面積對內彈道性能影響，得到模型2的速度—時間以及速度—位移曲線如圖5所示。

在35MPa下，泄流直徑為20mm時，加速過程不穩定，炮口初速無法滿足120m/s的要求，當泄流直徑由20mm增加到24mm時，加速過程變得明顯平穩，炮口初速由不足100m/s增加到119 m/s，當泄流直徑由24mm增加到28mm時，炮口初速由119m/s增加到121m/s。

2.3 仿真結果中提取的速度—位移關系

在分析上述結果基礎上為了得到合適的炮管長度，選取合適模型的位移—速度曲線如圖6所示。

在泄流直徑28mm、壓縮空氣壓力35MPa的情況下，彈藥在800mm位移處即能達到120m/s的速度。

3 試驗對比

根據仿真數據，建立試驗模型，由于35MPa氣體較難獲得，本次試驗采取20MPa試驗氣體作為對比，僅作軟件仿真與試驗對比。

通過試驗測試，得到時間—加速度數據圖，然后在軟件內做出相同模型仿真，圖7為試驗數據與軟件模擬在相同尺寸數據上的結果對比。

通過對比發現軟件仿真數據偏大，與試驗有一定誤差，但在整體趨勢上保持一致，并且誤差在可接受范圍之內，因此軟件有一定可信度，并且在初步設計階段可以接受。

4 結 論

本文通過在AUTODYN軟件中建立壓縮空氣發射裝置的內彈道模型，獲得了壓縮空氣壓力、發射管泄流直徑、彈藥位移、炮口初速之間的關系，得到了如下結論：

（1）在泄流直徑20mm的條件下，34MPa，35 MPa壓力的壓縮空氣加速過程不平穩，增加壓縮空氣的泄流直徑可以解決彈藥加速過程不平穩的問題；

（2）壓縮空氣壓力35MPa，當泄流直徑由20 mm增加到24mm時，加速過程變得明顯平穩，炮口初速由不足100m/s增加到119m/s，在泄流直徑20mm變化到24mm時，炮口初速急劇提高。當泄流直徑由24mm增加到28mm時，炮口初速由119m/s增加到121m/s，提升并不高；

（3）在泄流直徑28mm、壓縮空氣壓力35MPa的情況下，彈藥在400mm位移處即能達到100m/s的速度，發射裝置的壓縮空氣泄流直徑可設定為28mm，壓縮空氣壓力可選取為35MPa，發射管長度可設定為800mm。

參考文獻：

[1]李廣裕.戰略導彈彈射技術的發展[J].國外導彈與航天運載器，1990（7）：38-49.

[2]任淼，王秀萍.國外空空導彈發展動態研究[J].航空兵器，2013（5）：12-17.

[3]喬汝春.輕型艦載魚雷發射裝置管內魚雷運動參數的估算法[J].魚雷技術，1999（3）：31-34.

[4]許斌，楊積東，劉廣，等.機載導彈彈射式發射建模與仿真[J].系統仿真學報，2001，23（7）：51-54.

[5]廖振強，王濤，何大平，等.拋放彈彈射機構優化設計[J].南京理工大學學報，2003，27（5）：573-577.

[6]石少卿，汪敏，孫波，等.AUTODYN工程動力分析及應用實例[M].北京：中國建筑工業出版社，2012.

[7]ANSYInc.，DocumentationforRelease12.1[Z].ANSY，Inc，2009.

[8]賈孟義.火箭彈與發射管最小彈管間隙的計算[J].兵工學報.彈箭分冊，1981（2）：81-92.

[9]劉玉磊.燃氣舵流固耦合傳熱數值分析[J].航空兵器，2013（3）：41-43.

[10]ChoiM，LeeJungRyul，KongCheolWon.Development ofaNumericalModelforanExpandingTubewithLinear ExplosiveUsingAUTODYN[J/OL].ShockandVibration，vol.2014，ArticleID436156，10pages，2014. doi：10.1155/2014/436156.http：∥dx.doi.org/10. 1155/2014/436156.