筒壁溝槽式火炮制退機的流場數值仿真

摘"要：為了研究某火炮筒壁溝槽式制退機在火炮后坐運動過程中內部流體的流動狀態，研究不同后坐速度下制退力的特性，應用Fluent軟件建立了制退機的等效流場模型，由工作腔壓力和流量，推導出液壓阻力系數K并分析流速、溝槽面積和液壓阻力系數K的關系。對比分析仿真和實驗數據表明：無論活塞桿勻速還是變速運動，同截面上兩不同流道的流液速度幾乎相等，且同截面不同流道的流量比近似等于各流道的面積比。

關鍵詞：火炮；制退機；流場；仿真

中圖分類號：TJ302""文獻標志碼：B""文章編號：1671-5276（2024）02-0153-04

Numerical Simulation of Internal Flow Field of Barrel Grooved Gun Brake

XU Mingke

（School of Mechanical Engineering，Nanjing University of Science and Technology，Nanjing 210094，China）

Abstract：In order to study the flow state of the internal fluid of the grooved recoil brake on the barrel wall of an artillery during the recoil motion of the artillery and the characteristics of the recoil force at different recoil speeds， the equivalent flow field model of the recoil brake was established using Fluent software. From the pressure and flow rate of the working chamber， the hydraulic pressure resistance coefficient K was derived， and the relationship between the flow rate， groove area and hydraulic resistance coefficient K was analyzed. The comparative analysis of simulation and experimental data shows that no matter the piston rod moves at a constant speed or a variable speed， the flow velocity of two different channels on the same section is almost equal， and the flow ratio of different channels on the same section is approximately equal to the area ratio of each channel.

Keywords：artillery；brake retreat machine；flow field；simulation

0"引言

在現代戰爭中，功率增大、機動性提高是傳統火炮的發展趨勢［1］，而體積小、質量輕、射程遠、打擊類型多是五人戰車發展的重要方向［2］。隨著戰爭多元化，裝備火炮的威力不斷提高，其火炮裝藥量及后坐力也隨之增大。軟后坐火炮可以大幅減小炮架受力，從而使火炮設計得更輕巧，提高火炮射擊穩定性以及綜合指標［3］。

制退機是軟后坐火炮的關鍵組成部分，理論上軟后坐技術可將平均后坐阻力減至常規火炮的1/4，實際中可降至1/3到1/2［4］。制退機是火炮反后坐裝置的重要組成部分，是一種液壓阻尼器。在火炮后坐的運動過程中，制退機通過內部液體受到擠壓對火炮產生阻力，同時消耗后坐動能，將后坐運動限制在規定的長度內并控制后坐和復進運動的規律［5］，對提高火炮的射擊穩定性有重要作用。制退機可分為節制桿式、活塞溝槽式、筒壁溝槽式［6］。溝槽式的制退機力與溝槽的面積有關，目前工程上一般用帶損耗項的伯努利方程進行求解計算［7］。

由于制退機內部結構的復雜性以及工作時間的短暫性，使得制退機的工作規律難以控制［8］。制退機內部流場的規律是決定制退機工作性能的重中之重。因此制退液在工作時的運動狀態一直以來都是制退機相關研究方面的重點［9］。流體仿真作為設計驗證手段能較好地與試驗數據吻合［10］。

1"某帶流口的筒壁溝槽式制退機原理

1.1"制退機的結構

制退機的結構形式有很多種，通常按流液孔的形成方式分為節制桿式、溝槽式、活門式、轉閥式和多孔襯筒式。本文研究對象采用帶流口的筒壁溝槽式制退機，如圖1所示。該制退機的內筒壁上刻有6條溝槽，深度隨后坐距離而變化，且在活塞后有一前擋環，環內部周向環列有4個面積相等的中心流口。

筒壁溝槽深度隨后坐位移的增大從5.5mm線性減小到0.5mm。經計算，其對應單道溝槽流液通道面積計算公式為：

R=6+［（－5.5/900）v］t+75（1）

S=40（R2－400－752－202）－20R2－400+

［2arcsinθ（20/R）×180/π］πR2360－72.157 86（2）

式中：v為活塞后坐的速度；R為溝槽內表面的半徑；S為單條溝槽的面積。溝槽面積是關于后坐位移一一對應關系的函數，當給定一定速度時，式（2）則變為溝槽面積與時間的對應函數。

當活塞隨制退桿向后后坐運動時，工作腔內的液體受到擠壓，通過流口和溝槽兩條流道向中間腔流動，同時中間腔內的液體也受到擠壓向非工作腔流動，因此該制退機相較于以往制退機內部流體運動更加復雜。

1.2"制退機的工作原理及數學模型

制退機流場邊界運動復雜，運動速度快且幅度大，局部流動通道狹窄且形狀不規則，這給制退機的研究工作帶來了巨大的挑戰［11］。筒壁溝槽式制退機在火炮后坐過程中需要起到阻尼作用時其內部有兩條流液路徑，分別是刻在制退筒內壁上的溝槽和與活塞桿相固連的前擋環上的流口，其內部流液原理簡圖如圖2所示。

假設制退液為不可壓流體，對于后坐過程，有

Avp=Sgvg+Slvl（3）

式中：A為活塞受壓面積；vp為活塞絕對速度；Sg為中心流口面積；vg為制退液流過流口的平均絕對速度；Sl為溝槽總面積；vl為制退液流過溝槽的平均絕對速度。運用伯努利方程（流體的能量守恒定律）建立制退機工作腔的壓力計算公式。由圖2可見，筒壁溝槽通道分為兩段，分別是由工作腔內的液體分流流過的溝槽通道和最終匯流到非工作腔的溝槽通道，沿后坐流液路徑，根據伯努利公式有：

z1g+P1/ρ+v20/2=z2g+P2/ρ+v21/2+Hr12（4）

z3g+P1/ρ+v20/2=z4g+P2/ρ+v22/2+Hr34（5）

z5g+P2/ρ+v21/2+v22/2=z6g+P3/ρ+v23/2+Hr56（6）

式中：ρ為制退液密度；下標1代表工作腔；下標2代表前擋環與活塞中間的腔。式（4）和式（5）是液體流過工作腔的分流公式，式（6）是液體流入非工作腔的匯流公式。

由于工作腔和非工作腔重力勢能項變化相對較小，近似相等。假設工作腔液體的流速v0為0，非工作腔抽真空，其壓強P3為0，并假設損耗與動能成正比，即

Hr12=ξ1v21/2，Hr34=ξ2v22/2，Hr56=ξ3v23/2

經化簡得到：

P1/ρ=P2/ρ+K1v21/2（7）

P1/ρ=P2/ρ+K2v22/2（8）

P2/ρ+v21/2+v22/2=P3/ρ+K3v23/2（9）

式中K1、K2、K3分別為從工作腔分流到溝槽、分流到流口和匯流到非工作腔流道的液壓阻力系數。有研究者處理燃油時將供給管路中的K視為面積的函數［12］，實際上其與流動速度和流道的結構特征有關，且其值可在一定范圍內發生變化。

故工作腔壓力：

P1=ρ（K1v21/2+K3v23/2－v21/2－v22/2）（10）

或

P1=ρ（K2v22/2+K3v23/2－v21/2－v22/2）（11）

取式（10）計算制退機力：

F=A1ρ（K1v21/2+K3v23/2－v21/2－v22/2）（12）

以上所有公式中，v1、v2和v3由體積流量q1、q2和q3與對應截面的流液通道面積的比值計算，P1、P2和P3由Fluent監測得到。因此，知道了制退液的流速即可計算制退機力。有研究人員將帶旁路的閥控制退機的主流支流流速比值設為液壓阻力系數的1/2次方的反比，但是針對于氣體的流量分配研究表明，流量分配與過液面幾何因素和流速有關［13］。

2"制退機內部流體仿真

2.1"有限元建模及仿真設置

火炮后坐時，制退桿帶動活塞沿搖架導軌向后運動，液體受到擠壓后從工作腔流入非工作腔，在擠壓過程中非工作腔的空間不斷增大，而工作腔不斷被壓縮減?。?4］。利用有限元軟件ANSYS ICEM分別繪制制退機內腔和六道溝槽網格模型，最終合并成制退機流體腔網格模型。網格類型采用六面體網格，大多數網格尺寸約為3mm。由于整體尺寸較大和網格局部加密，制退機流體腔網格模型共包含600多萬個網格單元和600多萬個網格節點。

將建立的網格模型文件輸出格式設為msh文件類型，在Fluent中將網格模型導入，制退液密度取1 100kg/m3，制退液黏性系數取0.003 3kg/（m·s）。采用瞬態仿真的方法，令制退機以不同速度［2，4，6，8］m/s勻速后坐，將非工作腔端面設為壓強出口，出口高斯壓強設為0。監測前擋環工作面壓強P1和活塞工作面壓強P2，動網格設置如圖3所示。求解器類型為壓力基，湍流模型為標準k-ε模型，離散算法為SIMPLE，殘差最大值取0.001，制退液為不可壓液體，提交計算。

2.2"仿真結果分析

在傳統的計算方法中，K值的選取是根據經驗，參考同類型制退機，在一定范圍內選取，這導致計算結果和實際測量結果存在較大的誤差。采用數值方法對制退機內流場進行模擬計算，避免了選取制退機液壓阻力系數K，使液壓阻力更接近真實值［15］。影響筒壁溝槽式制退機主流流場變化的主要參數為制退機的結構特征、尺寸和后坐運動速度，其中結構特征和尺寸是固定的，但是在后坐運動過程中，后坐速度是不斷變化的［16］。活塞工作面壓強與溝槽總面積的對應關系如圖4所示（本刊黑白印刷，相關疑問請咨詢作者）。將仿真得到的相關數據代入式（7）、式（8）和式（9）中計算得到K1、K2和K3與溝槽面積的對應關系，如圖5所示。

K1=4.259 49－0.005 58x+4.176 61×10－6x2－9.771 18×10－10x3（13）

K2=3.438 13－0.003 62x+3.132 88×10－6x2－9.989 17×10－10x3（14）

K3=5.848 10－0.010 29x+1.049 79×10－5x2－9.771 18×10－9x3（15）

不同后坐速度下，中心流口通道與溝槽通道流量比值與流液通道面積的比值關系如圖6所示，兩比值接近于1∶1，近似為一條線性直線。使用Matlab對火炮后坐部分擊發全過程的運動進行仿真，仿真得到活塞桿的后坐速度，將后坐速度用profile文件的形式導入Fluent里驅動活塞運動進行流體仿真，后坐速度如圖7所示。將仿真得到的數據經過計算，代入到式（12）計算，用Matlab仿真計算得到制退機力并與仿真監測的活塞工作面靜壓面積進行對比，結果如圖8所示。

根據以上仿真結果對比發現，制退液從工作腔分流流入溝槽流道的流量q1和中心流道的流量q2的流量比約等于溝槽流道的流液通道面積與中心流道的流液通道面積的比值，且流入兩流道時的流速近似相等。這說明仿真過程監測的制退機力與本文模型計算得到的力不論在數值還是變化趨勢上都比較接近。

3"結語

本文利用ANSYS ICEM對某制退機內部的流體模型進行網格劃分，利用Fluent進行流體仿真計算，推導出了該模型的液壓阻力系數K并將模型計算得到的制退機力與仿真監測的進行對比分析。對于筒壁溝槽式制退機，可通過勻速仿真的方法得到液壓阻力系數與面積的關系且適用于變速仿真。對制退機內流場的變化機理作進一步的仿真研究，同時得到了制退機液壓阻力曲線，避免了人工憑經驗選取液壓阻力系數K，提高了計算精度。通過對比仿真和計算結果，液壓阻力近似相等，對今后制退機的設計及優化起到一定的參考和借鑒作用。

參考文獻：

［1］ QIU M，SI P，SONG J，et al. Recoil reduction method of Gun with side to rear jet controlled by piston motion［J］. Symmetry，2021，13（3）：396.

［2］ LIU J B，HUANG H S，ZHU W F，et al. Research on high efficiency muzzle brake technology of small caliber automatic Gun［J］. Journal of Physics-Conference Series，2020，1507（3）：032003.

［3］ 秦凱. 某軟后坐火炮動力學研究［D］. 南京：南京理工大學，2020.

［4］ 談樂斌，侯保林，陳衛民. 降低火炮后坐力技術概述［J］. 火炮發射與控制學報，2006，27（4）：69-72.

［5］ 張相炎. 武器發射系統設計理論［M］.北京：國防工業出版社，2015.

［6］ 高樹滋， 陳運生， 張月林，等. 火炮反后坐裝置設計［M］. 北京：兵器工業出版社，1995.

［7］ 郝保臣，孫衛國，潘孝斌. 某火炮液壓緩沖器緩沖特性研究［J］. 兵器裝備工程學報，2018，39（9）：54-58.

［8］ SEFIANE K，BENIELLI D，STEINCHEN A. A new mechanism for pool boiling crisis，recoil instability and contact angle influence［J］. Colloids and Surfaces A：Physicochemical and Engineering Aspects，1998，142（2/3）：361-373.

［9］ 林志宸. 新型閥控式制退機性能研究［D］. 南京：南京理工大學，2020.

［10］ 張曉東，張培林，傅建平，等. 基于動網格的火炮制退機內部流場數值模擬［J］. 南京理工大學學報（自然科學版），2010，34（4）：533-536，546.

［11］ 孫也尊，秦俊奇，狄長春，等. 基于計算流體力學的火炮制退機故障仿真研究［J］. 計算機仿真，2014，31（2）：31-34，96.

［12］ CHO J H. Engineering design and analysis of flow distribution in multiple flow circuits with partial contractions［J］. Journal of Mechanical Science and Technology，2021，35（9）：3979-3987.

［13］ 楊春信，席兵. 穩定流動的流量分配計算方法［J］. 流體機械，2000，28（2）：38-39.

［14］ 肖本勇. 一種變阻力火炮制退技術研究［D］. 南京：南京理工大學，2017.

［15］ 范永，劉樹華，曹廣群. 基于動網格的某駐退機三維流場數值模擬與分析［J］. 火炮發射與控制學報，2010，31（4）：63-65.

［16］ 潘孝斌，宋彥明，談樂斌. 筒壁溝槽式制退機主流液壓阻力系數分析［J］. 振動與沖擊，2016，35（17）：146-150.

收稿日期：20220906