某迫擊炮裝彈下滑過程流場仿真分析

2015-12-23 05:25:16馬杰,周中佳,卜華濤等

兵器裝備工程學報 2015年1期

【裝備理論與裝備技術】

某迫擊炮裝彈下滑過程流場仿真分析

馬杰，周中佳，卜華濤，何永

(南京理工大學機械工程學院，南京210094)

摘要：迫擊炮采用炮口裝填，裝填過程中迫彈在重力作用下下滑，使基本藥管與膛底的擊針相撞而發火。迫彈下滑的速度與炮彈定心部和膛內壁的間隙大小相關，間隙過大，火藥燃燒的一致性變差，炮彈出炮口的初始擾動增大，因此彈丸散度將加大；間隙過小，一方面不能保證在小發射角時的擊發動能，另一方面無法保證發射射速。采用Fluent流場分析軟件研究了迫擊炮彈、炮間隙對裝彈中下滑過程的影響，對幾種不同間隙和射角狀況下彈丸下滑過程的流場和彈丸運動規律進行了仿真分析，并開展了相關的討論。

關鍵詞：迫擊炮；流場；迫彈；間隙；耦合

收稿日期：2014-09-20

作者簡介：馬杰(1989—)，男，碩士研究生，主要從事計算流體力學研究。

doi:10.11809/scbgxb2015.01.019

中圖分類號：TJ31

文章編號：1006-0707(2015)01-0068-04

本文引用格式：馬杰，周中佳，卜華濤，等.某迫擊炮裝彈下滑過程流場仿真分析[J].四川兵工學報，2015(1):68-71.

Citationformat:MAJie,ZHOUZhong-jia,BUHua-tao,etal.NumericalSimulationandAnalysisofMortar’sInsideFlowFieldDuringLoading[J].JournalofSichuanOrdnance,2015(1):68-71.

NumericalSimulationandAnalysisofMortar’s

InsideFlowFieldDuringLoading

MAJie,ZHOUZhong-jia,BUHua-tao,HEYong

(MechanicalEngineeringCollege,NanjingUniversityofScienceandTechnology,Nanjing210094,China)

Abstract:The mortar uses the muzzle loading. During the loading, the mortar ammunition will speed down under the effect of gravity, during which the basic medicine impacts hammer needle at innersole to produce fire. The speed at innersole is related to the space between the mortar ammunition with gun barrel. The excessive gap seriously hampers the consistency of gunpowder combustion which increases the initial disturbance in muzzle, so the projectile dispersion will increase. At the same time, the space is too small to guarantee the firing kinetic energy and the firing rate. In order to study the influence of gap and structure to the loading process, the flow field around the pill ,the loading process in different gap and angle of departure were simulated with fluent, and related analyses were discussed at the same time.

Keywords:mortar;flowfield;mortarammunition;space;coupling

迫擊炮以其體積小、質量輕、使用和攜帶方便等優點，在未來戰爭中仍將發揮重要的作用。如何保證迫擊炮裝填時的狀態穩定和安全、可靠性一直是研究人員迫切需要解決的問題，本文將從流場分析的角度對炮膛間隙對迫彈裝填的影響開展研究，并進行相關的討論。本文考慮炮膛與定心部間隙造成的氣體流失以及彈丸附近湍流對裝彈過程的影響。對迫擊炮膛內氣體建立二維雷諾平均方程，使用k-ε二方程標準湍流模型，采用Roe的FirstOrderUpwind,應用動網格技術模擬了整個彈丸的下滑過程。

1迫擊炮裝填過程分析

裝填時，迫擊炮彈在下滑的過程中受到重力、膛內氣體阻力和炮彈與身管之間的摩擦力作用。其下滑運動方程為

(1)

式中：m為彈丸質量；f為摩擦系數；po為大氣壓力；px為彈丸下滑中膛內氣體壓力；θ0迫擊炮的發射角；S為迫擊炮彈的橫截面積。近似地，在迫擊炮設計理論中提出了一種假定的下滑速度與時間關系，如圖1所示。

圖1　迫彈假定下滑速度-時間曲線

如圖1所示，假定迫彈開始為勻加速運動，當運動至0.3 s時，速度到最大值，此后迫彈在膛內以近似勻速運動，因而可得假定下滑速度

(2)

μ=αβ

(3)

式中：f為摩擦系數取0.2；Δ為炮彈與膛壁之間直徑上的間隙；a為消耗系數，表示縫隙出口處氣體密度為大氣密度的倍數；β為氣體流出的速度系數[1-2]。μ約在0.6～0.8范圍內，式(2)表明間隙的大小直接影響下滑速度，一般來說間隙增大，膛內氣體容易排出，加速度較大，當間隙足夠大是不出現近似等速運動，間隙減小時，氣體不容易排出，氣體受壓縮后，下滑運動出現較為明顯的振動現象，如圖2所示。

圖2　不同間隙大小情況下的迫彈下滑曲線

根據迫擊炮設計理論計算可得，120 mm迫擊炮的最佳間隙為0.6～0.8。

2迫彈裝填過程氣動數學模型

2.1控制方程

考慮到整個裝填過程，膛內間隙處的氣體馬赫數從0增長到1，同時彈帶前部的氣體經過加速，氣體馬赫數>1,所以氣體視為可壓縮氣體。根據以上假設建立迫彈流場二維軸對稱守恒形式的(N-S)方程

(4)

(5)

F=[ρμ ρμ2+p-τxxρυμ-τxy

(6)

G=[ρv ρυμ-τyxρμ2+p-τyy

(7)

(8)

氣體狀態方程

Pv=nRT

(9)

量熱狀態方程

e=cvT

(10)

式中：ρ為大氣密度；μ為膛內氣體的徑向流速；υ為膛內氣體的軸向流速；e為氣體單位比內能；p為壓強；f為作用在單位質量流體微團的體積力[3]。

Gk+Gb-ρε-YM+Sk

(11)

(12)

式中:Gk表示由層流速度梯度產生的湍流動能；Gb是由浮力產生的湍流動能；YM由于在可壓縮湍流中，過度的擴散產生的波動C1,C2,C3是常量；σk和σε是k方程和ε方程的湍流prandtl數；Sk,Se是用戶定義的[4]。

2.2計算區域網格劃分方法

計算區域網格的劃分，考慮到氣體和迫彈的流固耦合作用，采用結構化動態層模型：

Hmin>(1+as)*h0

(13)

Hmin

(14)

式中：ac為合并因子，取0.04；as為分割因子，取0.4；h0為理想單元高度，取0.005 m；Hmin為單元的最小高度。在網格運動條件下，流場內控制體發生變化，采用有限容積法來表示守恒方程[5]。

以下為計算區域分塊網格的劃分方法。炮口和膛底為靜態區域，彈丸和流體區域為剛體運動區域，剛體運動區域的運動通過編程定義，運用牛頓第二定律計算膛內氣體對彈丸的阻力。為了能夠捕捉到間隙的流場，對定心部與身管之間的間隙進行加密處理，如圖3所示。

3迫彈裝彈下滑過程流場模擬結果與分析

本文分別模擬和分析了發射角為45°,間隙為0.6mm、0.8mm、0.4mm的3種裝填下滑過程工況，通過fluent流體分析軟件計算并獲取了彈丸下滑過程的流場等運動狀態，圖4為流場速度等值線圖。

圖3　迫彈下滑模型示意圖

圖4　彈丸在5.939 9 e-1( s)時膛內的速度等值線

可以看出整個流場除了間隙處速度比較大，其他區域速度都接近0。為了進一步分析流場分布規律，對區域一、二進行放大處理，如圖5、圖6所示。

圖5　區域一的速度等值線

圖6　區域二的速度等值線

可以看出，膛內氣體在密封圈之前間隙逐漸減小的區域，壓力逐漸減小，速度逐漸增大，在進入密封圈區域之后，截面積增大，此時氣流馬赫數仍小于1，速度有略微的減小。之后氣體通過等效縮擴管的密封圈之后的區域，氣體加速為超音速氣體。

為了能夠清晰了解彈丸在下滑過程中，主要流場區域的流動特性，從fluent流體軟件中提取了0.539s的速度流線圖，如圖7所示。

圖7　0.539 s時彈丸附近的速度流線

如圖7所示，彈丸在整個下滑的過程中，除間隙附近外，整個流體區域的氣體流速較小，根據雷諾數計算公式可得此流體區域內雷諾數遠小于2 320，因此流動可視為層流運動。

從速度等值線圖得知，間隙處在t=0.4之后，彈丸間隙末端的氣體加速為超音速氣體，雷諾數在3e4左右遠大于湍流臨界值，通過圖8發現，在密封圈處的湍流動能急速增加。在出密封圈的附近湍流達到最大值約為3.25e3(m2/s2)。

圖8　0.539 s時彈丸密封圈附近的湍流動能等值線

膛內壓力是影響彈丸運動軌跡的主要因素之一，通過fluent流體計算軟件，得到了膛內靠近膛底某一點壓力-時間曲線，如圖9所示。

圖9　Δ=0.6靠近膛底部位的壓力-時間的變化曲線

可以看出，在整個裝填過程中，堂內壓力變化同樣經歷了緩慢-快速-緩慢的一個升壓過程。利用fluent軟件，得到了彈丸在不同間隙的情況下，發射角為45°的速度-時間曲線和彈丸受力-時間曲線圖，如圖10、圖11所示。

圖10　不同間隙下的v-t曲線

圖11　不同間隙下迫彈周圍整體的p-t曲線

從v-t曲線可以看出迫彈在經過近似勻加速后減速運動，通過流體分析的結果與理論分析的結果大體吻合。從p-t曲線可以看出，迫彈周圍整體壓力在經過一段時間緩慢增加后到0.4s左右進入了一個快速增加的階段。此時彈丸整體壓力與摩擦的合力仍小于重力，所以彈丸會進入一個加速度減小的加速階段。從速度曲線看出在t=0.5s左右速度達到最大值，此時膛內壓力達到臨界值，彈丸周圍整體壓力與摩擦力的合力與重力平衡，之后，重力小于阻力，彈丸進入加速度逐漸增大的減速階段。通過對模型和計算結果的分析，得出彈丸間隙附近的流體區域可以簡化為一個給定面積比的縮擴管，管后氣體壓力基本保持不變，進口氣體總壓從p0開始逐漸提高，流體區域從過度膨脹工作狀態的第三類開始，然后隨著膛內壓力的升高，依次出現第二類、第一類和完全膨脹工作狀態、不完全膨脹的工作狀態的流動情形。仿真結果表明間隙0.6mm，t=0.4s左右流體區域進入不完全膨脹的工作狀態，最小界面處為馬赫數為1的臨界狀態，根據縮擴管的氣體流量公式