孫曉雄，田曉麗，沈劍，白晨，張軼夫

(中北大學(xué) 機電工程學(xué)院，山西 太原 030051)

火炮武器系統(tǒng)具有初速快、射速快、適應(yīng)性強以及可靠性強等優(yōu)點，在現(xiàn)代戰(zhàn)爭中大放異彩[1-3]。藥筒是火炮彈藥的重要組成部分之一，對提高火炮的發(fā)射速度，保護發(fā)射裝藥，減少火藥對炮膛的燒蝕，以及提高火炮的使用時間等方面起著關(guān)鍵作用。大口徑火炮所用的藥筒仍為鋼制結(jié)構(gòu)藥筒,在射擊中，由于藥筒的筒體強度低，退殼性能差，如若不能順利完成抽殼動作，造成卡殼，這就未達到藥筒的設(shè)計要求，也勢必會對火炮的連續(xù)射擊造成負面影響[4]。

針對藥筒存在的上述問題，王明廣等[5]對某火炮的抽殼過程進行數(shù)值仿真，分析了彈殼彈膛系統(tǒng)溫度差異對抽殼力的影響。譚波等[6]重點剖析了艦炮抽筒模板下移和抽筒左右方位的改變對抽殼產(chǎn)生的影響。楊吉林等[7]使用ANSYS軟件構(gòu)建了焊接鋼制結(jié)構(gòu)藥筒模型，研究溫度和壓力對藥筒力學(xué)特性造成的變化。李翔等[8]對比了貼膛時間與不同膛壓的關(guān)系，計算提取了艦炮不同區(qū)域的殘余正壓力，并分析了二者對抽殼的影響。倪路瑤等[9]研究了彈藥發(fā)射后，藥筒在膛壓作用下的變形情況，通過動力學(xué)仿真，擬合出了抽殼力曲線。

筆者基于中北大學(xué)開發(fā)的藥筒有限元分析系統(tǒng)CCAE(Cartridge Computer Aided Engineering)，分析藥筒劃分不同區(qū)域?qū)Τ闅ば阅芩鶐淼牟町?，研究藥筒在不同材料、不同溫度以及同一種材料不同材料參數(shù)對抽殼力的影響。

1 藥筒抽殼過程理論分析

1.1 有限元計算理論模型

焊接鋼質(zhì)藥筒是由直焊縫與環(huán)焊縫組合而成的錐體，其結(jié)構(gòu)和身管是帶有一定錐度的圓筒體，如圖1所示。

為了簡化問題，提出下述假設(shè)：藥筒所承受的載荷是平均載荷，僅是時間函數(shù)，在軸向和徑向都無變化；藥筒與身管是軸對稱模型；忽略藥筒底部中心點火孔的影響；忽略身管的后坐力在發(fā)射過程中的影響。基于以上假設(shè)，根據(jù)藥筒和身管結(jié)構(gòu)的對稱性，彈丸發(fā)射過程中外部施加載荷作用的對稱性，將圖1構(gòu)建的立體模型精簡為二維對稱的平面模型，如圖2所示。

在CCAE有限元計算中，根據(jù)變分原理，可推得彈塑性有限元求解方程[10]：

(1)

式中：[B]表示幾何矩陣；Δσij表示應(yīng)力張量；[N]是形函數(shù)矩陣；Sσ為身管或藥筒的表面積；Δqi表示身管和藥筒表面所承受的面積力密度；{pi}表示集中力；{W}表示體積力密度。

根據(jù)Von-Mises準則得

(2)

式中：Δσij為歐拉應(yīng)力張量；E為彈性模量；μ為波桑比；δik，δjl，δkl為克羅內(nèi)爾記號；Δεkl為格林應(yīng)變張量；dλ為應(yīng)變張量與應(yīng)力偏量間的比例因子；σ′ij為應(yīng)力偏量；α*是載荷性質(zhì)的判斷因子，取1。

依據(jù)Prandtl-Reuss彈塑性應(yīng)力-應(yīng)變關(guān)系得

(3)

式中：H為硬化系數(shù)；G為剪切彈性模量。從而可得到增量形式的彈塑性應(yīng)力-應(yīng)變關(guān)系方程：

(4)

1.2 藥筒抽殼過程分析

藥筒在發(fā)射時，根據(jù)火藥氣體壓力的變化和藥筒的工作情況，可將藥筒的作用過程分為4個時期。

第1時期：底火發(fā)火引燃發(fā)射藥，膛內(nèi)壓力開始增加，迫使藥筒變形，直到與身管相接觸。接著藥筒開始彈性變形，而后為塑性變形，藥筒材料性能將發(fā)生強化現(xiàn)象。

第2時期：即從藥筒與身管接觸到膛壓最大值為止，此時期藥筒與身管一起變形，對身管是彈性變形，對藥筒則是塑性變形，藥筒與身管達到更緊的貼合，藥筒壁受到徑向壓力，產(chǎn)生壓縮變形。

第3時期：從最大膛壓到膛壓降到大氣壓為止，膛壓完全下降，使身管最終變形位置恢復(fù)到原來的位置。此時由于藥筒已經(jīng)產(chǎn)生了塑性變形，它不可能恢復(fù)到原來位置。繼續(xù)受熱從膛壓下降至大氣壓到抽殼完畢，此時期結(jié)束會出現(xiàn)兩種情況：一種是藥筒恢復(fù)后的尺寸小于身管的尺寸，因而產(chǎn)生間隙，對抽殼有利；另一種是藥筒恢復(fù)的尺寸仍大于身管尺寸，因而產(chǎn)生過盈，對抽殼不利。

第4時期：藥筒和身管形成最終間隙，藥筒順利抽殼。

由以上分析不難看出，在藥筒的作用中，只要能形成有利的最終間隙，會給抽殼帶來方便，抽殼力的大小又會影響最終間隙。

1.3 影響抽殼性能的因素分析

藥筒的抽殼性能主要與藥筒機械性能、初始間隙、火藥氣體壓力、藥筒與身管間摩擦力的大小等因素相關(guān)[11]。

1)藥筒機械性能的影響。根據(jù)金屬材料的機械性能可知，藥筒材料的強度極限越高，發(fā)射后還原的彈性變形就越大，即藥筒與身管形成的利于藥筒順利抽出的最終間隙越大，因而容易抽殼。

2)初始間隙的影響。初始間隙的增大將會使藥筒的最終間隙也增大，但是初始間隙過大會造成閉氣性能不良和增加藥筒縱向破裂的可能。

3)火藥氣體壓力的影響。火藥氣體壓力的大小，是決定藥筒抽殼性能和最終間隙的主要因素。當藥筒強度一定時，膛壓增加，藥筒的殘余變形變大，因而形成的最終間隙就小，甚至產(chǎn)生“過盈”而卡殼。

4)藥筒和身管摩擦力大小的影響。表面粗糙或者質(zhì)量不佳，如銹蝕、劃痕等造成摩擦阻力增加，從而影響抽殼。

金屬藥筒在射擊時的變形，是由大小不一的幾部分組成的，為適應(yīng)金屬藥筒抽殼性能的需要，有必要根據(jù)藥筒在膛內(nèi)變形的不同特點，將藥筒分成若干個區(qū)，以便進行抽殼分析。

采用分區(qū)的布局是為了方便顯示與比較不同區(qū)域的變形與位移，同時能夠更加直觀的比較藥筒不同區(qū)域與身管之間接觸作用的全過程，將藥筒區(qū)域劃分為4個區(qū)，分別為：藥筒口部、藥筒斜肩部、藥筒體部和藥筒底部。如圖3所示。

2 有限元計算前處理

2.1 有限元模型

藥筒壁厚很薄，藥筒口部、藥筒斜肩部和藥筒體部最厚的壁厚也只有2.3 mm，但整個藥筒長度可達700 mm，為了控制網(wǎng)格質(zhì)量，考慮藥筒和身管的結(jié)構(gòu)和尺寸存在很大的差異，在CCAE中將二者劃分成獨立的部分。考慮到身管在火炮發(fā)射過程變形較小，所以身管網(wǎng)格尺寸相對粗糙。網(wǎng)格劃分如圖4所示。

2.2 不同鋼材材料模型

在構(gòu)建藥筒和身管的具體材料參數(shù)時依據(jù)如下兩點：

1)藥筒包含4個材料區(qū)，每個區(qū)域內(nèi)材料具有一致性，藥筒的材料模型選用MAT_PLASTIC_KINEMATIC_TITLE，該模型適用于各向同性和隨動強化塑性模型，并考慮了速率效應(yīng)，是一個非常經(jīng)濟有效的模型。

2)整個身管選用材料模型為線彈性的MAT_ELASTIC_TITLE，其是一種各向同性的彈性材料，可用于梁和LS-DYNA中的實體單元，也就是材料的變形在材料外部的載荷去除后可全部恢復(fù)，應(yīng)力應(yīng)變呈線性關(guān)系。

2.3 不同物理參數(shù)材料模型

計算中對鋼材的材料參數(shù)的切線模量進行修改，通過統(tǒng)計藥筒各個分區(qū)與身管接觸面的殘余正壓力來計算最終的抽殼力，對比分析不同材料參數(shù)對抽殼力的影響情況。材料模型參數(shù)如表1所示，修改后的不同鋼材的材料參數(shù)如表2所示，其余部分參數(shù)與表1一致。

表1 20#/30#/45#鋼材料參數(shù)

表2 20#/30#/45#鋼修改材料參數(shù)

2.4 不同溫度下膛壓曲線

火炮發(fā)射時的火藥溫度能夠瞬間達3 000～4 000 K，藥筒內(nèi)壁會在1～10 ms的時間內(nèi)受到高溫高壓氣體的沖擊作用，使彈殼內(nèi)壁溫度瞬間升高，且短時間內(nèi)來不及傳遞，造成彈殼內(nèi)壁溫度高，膨脹大；而彈殼外壁溫度低，膨脹小。因此，有必要研究不同溫度下藥筒的抽殼性能，不同溫度下的膛壓曲線如圖5所示。

實驗選取了高溫(50 ℃)和常溫(15 ℃)2種工作溫度，將火藥氣體壓力以壓力載荷的形式加在藥筒內(nèi)壁，常溫條件下在3.31 ms時最大膛壓約為267 MPa，高溫條件下在3.52 ms時最大膛壓約為331 MPa。

3 計算結(jié)果

3.1 抽殼力計算公式

利用CCAE藥筒有限元軟件對藥筒抽殼力進行統(tǒng)計計算，其計算公式為

R=∑fiFi,

(5)

式中：R表示藥筒抽殼力；i表示接觸面編號；fi表示接觸面i處摩擦系數(shù)；Fi表示接觸面i處殘余正壓力。

3.2 藥筒不同分區(qū)抽殼力實驗結(jié)果

藥筒抽殼過程應(yīng)力云圖如圖6～8所示，因篇幅有限，只顯示20#鋼在常溫條件下，藥筒貼膛前、藥筒最大膛壓時和藥筒卸載后各部位應(yīng)力分布云圖。從數(shù)值模擬得到的結(jié)果可以看出，不論是藥筒貼膛前、藥筒最大膛壓時還是藥筒卸載后，藥筒口部對應(yīng)數(shù)值顯示的應(yīng)力云圖均為最大的，藥筒斜肩部次之，其中藥筒底部的應(yīng)力最小。

3.3 不同鋼材對抽殼力的影響

依照式(5)的抽殼力計算公式，統(tǒng)計得到3種鋼材在常溫條件下卸載后不同部位的抽殼力如圖9所示。從圖9可得，卸載后抽殼力主要分布在藥筒口部和藥筒斜肩部，30#鋼的抽殼力比20#鋼、45#鋼的抽殼力稍微略低，相對而言抽殼較為順利，但是藥筒總體的抽殼力仍然偏高。

3.4 不同物理參數(shù)材料模型對抽殼力的影響

依據(jù)實驗要求，藥筒總體抽殼力需在18 kN以下，因此對鋼材的切線模量進行修改。在常溫條件下，統(tǒng)計不同鋼材卸載后，未修改材料參數(shù)和修改材料參數(shù)的抽殼力對比如圖10所示。

從圖10可以很直觀地看出，修改材料部分參數(shù)后，3種鋼材的抽殼力下降很明顯，其中3種鋼材藥筒口部的抽殼力從大約20 kN下降到約10 kN,藥筒斜肩部抽殼力從大約10 kN下降到約5 kN，藥筒體部從3 kN降到約2 kN，藥筒底部變化不明顯。整個藥筒的抽殼力下降約50%左右，對于某些抽殼困難的彈殼來說，抽殼力下降意味著降低了卡殼的出現(xiàn)率。

依據(jù)仿真計算的藥筒各部位抽殼力數(shù)值可以得出，藥筒口部、藥筒斜肩部和藥筒體部采用30#鋼，藥筒底部采用20#鋼可以達到抽殼力最低的效果。

3.5 不同溫度對抽殼力的影響

不同溫度下，20#鋼以及修改20#鋼材料參數(shù)的抽殼力對比如圖11所示。從圖11可以看出無論是否修改材料參數(shù)，藥筒的各部位抽殼力都隨溫度的升高而升高。但是在高溫條件下，藥筒口部抽殼力從27.5 kN下降到19.25 kN，藥筒斜肩部抽殼力從14.47 kN下降到10.13 kN,藥筒體部和藥筒底部變化不明顯，修改部分材料參數(shù)仍然可以達到降低抽殼力的效果。

4 結(jié)論

通過對藥筒抽殼過程進行理論分析與有限元計算仿真發(fā)現(xiàn)，藥筒口部和藥筒斜肩部對抽殼力的影響最大。因此，在藥筒的設(shè)計中，可以考慮不同部位采用不同鋼材，從而達到順利抽殼的目的，結(jié)論如下：

1)換用不同鋼材對藥筒的抽殼力有一定的影響，

其中30#鋼相比20#鋼、45#鋼而言比較容易抽殼。

2)統(tǒng)計了3種不同鋼材對應(yīng)的修改材料的抽殼力，藥筒口部和斜肩部的抽殼力下降十分明顯。

3)藥筒各個部位的抽殼力都隨溫度升高而升高，溫度對抽殼力的影響不可忽視。