陳朝君，鄭建國，丁傳俊

（南京理工大學(xué)機(jī)械工程學(xué)院，江蘇南京 210094）

0 引言

制退機(jī)是火炮上控制受力和運(yùn)動的重要部件，通過擠壓制退液使其流過流液孔以達(dá)到耗能的效果，從而保證火炮射擊穩(wěn)定性。由于火炮射擊時，制退機(jī)后坐壓力達(dá)到20 MPa左右，后坐速度達(dá)到13 m/s左右，整個后坐時間持續(xù)0.2 s左右［1］，所以其內(nèi)部的流動規(guī)律比較復(fù)雜。

目前國內(nèi)的研究主要是以仿真為主，其中鄭建國［2］等建立二維軸對稱模型，并使用“仿運(yùn)動壁面法”，對制退機(jī)內(nèi)部流場進(jìn)行了仿真計算，并得出液壓阻力系數(shù)。張曉東［4-6］等通過建立真實制退機(jī)模型，并使用動網(wǎng)格技術(shù)，得到了較好的結(jié)果。范永［3］等建立了三維簡化模型，并采用動網(wǎng)格技術(shù)得出了制退機(jī)工作腔最大壓力和最大流速。但是以上研究和制退機(jī)實際的工作情況還有較大差距，主要的原因有:1）模型過于簡化，無法反應(yīng)制退機(jī)內(nèi)部真實結(jié)構(gòu);2）采用不適當(dāng)?shù)倪吔鐥l件，以上文獻(xiàn)將制退機(jī)非工作腔簡化為無端蓋管道，充滿制退液且出口壓力為零。3）制退機(jī)的非工作腔并非一直存在真空，而是一個壓力震蕩的過程。

本文建立了與實際模型主要特征相近的簡化物理模型，放棄以往采用零壓力出口的方法，利用fluent軟件網(wǎng)格動態(tài)層更新方法［7］，對制退機(jī)的內(nèi)部流場進(jìn)行了非定常仿真，并得到了制退機(jī)內(nèi)部工作腔、非工作腔的壓力、速度云圖以及非工作腔射流的速度曲線。仿真結(jié)果顯示網(wǎng)格運(yùn)動過程和制退機(jī)實際運(yùn)動過程一致，體現(xiàn)了制退機(jī)內(nèi)部真實流動情況，為下一步結(jié)構(gòu)的優(yōu)化設(shè)計和故障分析提供了理論依據(jù)。

1 制退機(jī)簡化模型以及計算模型的處理

火炮發(fā)射時，后座部分沿導(dǎo)軌向后運(yùn)動，制退桿活塞擠壓Ⅰ腔內(nèi)的制退液，導(dǎo)致Ⅰ腔工作壓力增大，使制退液分為兩股液流，其中一股流入非工作腔Ⅱ，另一股液流流入Ⅲ腔。同時由于制退桿的不斷被抽出，Ⅱ腔會出現(xiàn)真空。Ⅲ腔由于始終充滿液體不會產(chǎn)生真空。

圖1 制退機(jī)結(jié)構(gòu)簡圖

制退機(jī)內(nèi)部的流動屬于三維非定常流動，有運(yùn)動邊界且存在高速真空射流。為了簡化計算，對制退機(jī)的計算模型進(jìn)行了簡化:首先將三維對稱結(jié)構(gòu)簡化為二維軸對稱結(jié)構(gòu);其次簡化節(jié)制桿尺寸和調(diào)速筒內(nèi)部結(jié)構(gòu);同時將制退桿活塞簡化為環(huán)形通道，但傾角為實際尺寸。簡化后的模型如圖2所示，在基本保持制退機(jī)結(jié)構(gòu)尺寸的情況下，可以反映制退機(jī)內(nèi)部流場的真實情況。

圖2 簡化后的計算模型

1）制退機(jī)的運(yùn)動邊界采用導(dǎo)入profile文件的方法驅(qū)動邊界運(yùn)動，使得制退桿和制退活塞按照實測炮管后坐速度抽出。在抽出的過程中，運(yùn)動區(qū)域的網(wǎng)格采用動態(tài)層變方法自動更新。例如，當(dāng)工作腔內(nèi)的網(wǎng)格受到壓縮時，網(wǎng)格會壓縮變形繼而和鄰近層的網(wǎng)格合并;非工作腔內(nèi)的網(wǎng)格則受到拉伸時，網(wǎng)格會拉伸然后分裂生成新的網(wǎng)格。

2）通過對模型內(nèi)部區(qū)域進(jìn)行分區(qū)來劃分網(wǎng)格，靜區(qū)域內(nèi)采用規(guī)則的四邊形網(wǎng)格，并且控制網(wǎng)格的疏密程度，使貼近壁面的網(wǎng)格更加精細(xì)，以此來代替流動的邊界層。對于不規(guī)則形狀的動區(qū)域則采用三角形非結(jié)構(gòu)網(wǎng)格，并嚴(yán)格控制網(wǎng)格品質(zhì)。動區(qū)域和靜區(qū)域之間采用interface進(jìn)行數(shù)據(jù)交換。

3）除了對稱軸為對稱邊界條件外，其他內(nèi)部邊界均作為壁面邊界條件處理。壁面有動壁面和靜壁面之分。

最后整個計算區(qū)域共生成89 184個網(wǎng)格，網(wǎng)格劃分足夠細(xì)密，可以使求解結(jié)果不存在網(wǎng)格依賴性。

在此將制退機(jī)的計算模型采用二維軸對稱模型的一半，制退機(jī)劃分的網(wǎng)格模型如圖3所示，其中制退機(jī)流場計算域網(wǎng)格如圖3（a），其中網(wǎng)格的局部放大圖如圖3（b）、（c）。為了數(shù)值模擬的有效性必須保證網(wǎng)格的品質(zhì)，因此采用以下原則:

圖3 制退機(jī)網(wǎng)格模型

1）劃分計算模型的網(wǎng)格采用分區(qū)劃分網(wǎng)格;

2）根據(jù)制退機(jī)結(jié)構(gòu)的復(fù)雜情況，對網(wǎng)格進(jìn)行加密處理;

3）布置網(wǎng)格時，結(jié)構(gòu)規(guī)則的地方采用四邊形網(wǎng)格，結(jié)構(gòu)不規(guī)則的地方采用三角形網(wǎng)格;

4）計算過程中各子區(qū)域數(shù)據(jù)通過滑移網(wǎng)格界面進(jìn)行數(shù)據(jù)交換。

2 計算采用的數(shù)值方法

仿真是在的正常氣溫下進(jìn)行，采用標(biāo)準(zhǔn)裝藥，0°射角。制退液設(shè)定為不可壓縮粘性流體，密度為1 160 kg·m－3動力粘度為0.01 Pa·s。在軟件中采用非耦合算法求解控制方程。采用標(biāo)準(zhǔn)湍流模型，壓力項采用PRESTO方法進(jìn)行離散，動量項采用一階迎風(fēng)格式離散，其他項均采用二階迎風(fēng)格式進(jìn)行離散。時間步長選擇為1e－06 s，整個迭代次數(shù)160 000，計算時間為0.16 s。`

3 仿真結(jié)果分析

仿真后得到各腔的壓力、速度、湍動能等分布云圖以及流線圖。圖4給出了制退機(jī)后坐速度和非工作腔射流速度曲線;圖5給出了制退機(jī)工作腔平均壓力曲線;圖6給出了制退機(jī)非工作腔的平均壓力。

圖4、圖5中可以看出當(dāng)制退機(jī)后坐速度達(dá)到最大時，制退機(jī)非工作腔內(nèi)的射流速度也達(dá)到最大值，即在0.013 6 s時非工作腔射流速度達(dá)到最大值105.23 m/s。而工作腔的最大壓力并不是出現(xiàn)在此時，在0.012 6 s時，工作腔平均壓力的最大值為8.24 MPa，在后坐速度最大時工作腔平均壓力為7.76 MPa，兩者相差并不大。

圖4和圖5中制退機(jī)后坐速度和非工作腔射流速度曲線呈現(xiàn)線性變化趨勢，制退機(jī)工作腔內(nèi)的平均壓力的變化趨勢則是近似二次變化的，這是符合設(shè)計理論的。

圖6非工作腔的平均壓力出現(xiàn)了短暫的負(fù)壓值，這是由于工作腔壓力急劇升高，制退液在高壓作用下通過流液孔進(jìn)入Ⅱ腔和Ⅲ腔，此時Ⅱ腔內(nèi)的壓力開始增大，隨著制退桿的不斷抽出，Ⅱ腔持續(xù)形成短暫的真空，在圖中即表現(xiàn)為壓力出現(xiàn)了短暫的負(fù)壓值，隨后Ⅱ腔內(nèi)的壓力基本保持低壓狀態(tài)。這是符合制退機(jī)工作情況的。

圖6 非工作腔平均壓力

圖7、圖8中可以看出制退機(jī)射流速度在流液孔附近迅速達(dá)到最大值，因此這股液流為主流。當(dāng)液流進(jìn)入Ⅱ腔后，由于其速度很高，當(dāng)它和Ⅱ腔內(nèi)的液體迅速混合后，在活塞附近立刻形成湍動能很高的漩渦，通過這一過程液體壓力、速度下降，制退液溫度升高，火炮動能由此轉(zhuǎn)化為液體的內(nèi)能，最終達(dá)到耗能制動的目的。另一股液流進(jìn)入Ⅲ腔，即為支流。由于Ⅲ腔在設(shè)計時，確保了其腔內(nèi)始終充滿液體，所以在后坐過程中，雖然制退桿不斷抽出，Ⅲ腔內(nèi)部空間增大，但是其內(nèi)部不會出現(xiàn)真空。

圖9中當(dāng)液體流過流液孔進(jìn)入Ⅲ腔，在經(jīng)過調(diào)速筒流液孔時，壓力、速度均下降。隨著時間的推移，制退機(jī)內(nèi)部的壓力、速度均不斷下降，漩渦的劇烈程度也有所下降。同時在Ⅱ腔，其內(nèi)部只會出現(xiàn)一個漩渦，這是因為非工作腔的流動類似于臺階流動，下游的漩渦會被較高速度的射流沖散，而右端如果采用零壓出口，在流動的后期，由于非工作腔內(nèi)部空間增大，腔內(nèi)壓力分布不均，則會出現(xiàn)兩個漩渦，這與理論不符。

圖9 不同時刻流線分布圖

4 結(jié)論

通過對制退機(jī)二維軸對稱模型的非定常仿真，較為真實地得出了通常實驗難以測定了流動數(shù)據(jù)，如壓力、速度以及真空射流速度等數(shù)據(jù)。而且通過分析認(rèn)為非工作腔只會出現(xiàn)一個漩渦，認(rèn)為在仿真中設(shè)置右端為零出口的邊界條件并不合適。因此，本文提出的仿真方法可以較為合理地仿真制退機(jī)內(nèi)部流動情況，可以為制退機(jī)的優(yōu)化設(shè)計和故障分析提供一定的理論基礎(chǔ)。

［1］高樹滋，陳運(yùn)生，張月林，等.火炮反后座裝置設(shè)計［M］.北京:兵器工業(yè)出版社，1995.

［2］鄭建國.火炮制退機(jī)流場的數(shù)值模擬［J］.力學(xué)與實踐，2001，23（2）:30-32.

［3］范永，劉樹華，曹廣群.基于動網(wǎng)格的某駐退機(jī)三維流場數(shù)值模擬與分析［J］.火炮發(fā)射與控制學(xué)報，2010.

［4］張曉東，張培林，傅建平，等.基于動網(wǎng)格的火炮制退機(jī)內(nèi)部流場數(shù)值模擬［J］.南京理工大學(xué)學(xué)報，2010，34（4）.

［5］張曉東，張培林，傅建平，等.基于二維模型的火炮沖擊運(yùn)動計算［J］.振動與沖擊，2011，30（2）.

［6］張曉東，張培林，傅建平，等.基于CFD與協(xié)同仿真的火炮后坐分析計算［J］.彈道學(xué)報，2010，22（3）.

［7］溫正，石良臣，任毅如.FLUENT流體計算應(yīng)用教程［M］.北京:清華大學(xué)出版社，2009.