徐赫陽(yáng)，喬相信，郭克強(qiáng)，李文浩，彭琳茜

(沈陽(yáng)理工大學(xué) 裝備工程學(xué)院，沈陽(yáng) 110159)

?

某線膛火炮戰(zhàn)斗部發(fā)射強(qiáng)度數(shù)值模擬

徐赫陽(yáng)，喬相信，郭克強(qiáng)，李文浩，彭琳茜

(沈陽(yáng)理工大學(xué) 裝備工程學(xué)院，沈陽(yáng) 110159)

為研究某新型戰(zhàn)斗部的發(fā)射強(qiáng)度，采用LS-DYNA軟件建立戰(zhàn)斗部在膛內(nèi)運(yùn)動(dòng)模型，對(duì)該模型進(jìn)行了發(fā)射強(qiáng)度仿真和分析，得到了戰(zhàn)斗部在不同運(yùn)動(dòng)階段的應(yīng)力、應(yīng)變曲線。分析第一、第二臨界狀態(tài)時(shí)戰(zhàn)斗部的應(yīng)力、應(yīng)變值及相應(yīng)強(qiáng)度準(zhǔn)則，并進(jìn)行了實(shí)驗(yàn)驗(yàn)證，模擬結(jié)果與實(shí)彈實(shí)驗(yàn)結(jié)果基本吻合。結(jié)果表明，戰(zhàn)斗部實(shí)驗(yàn)和數(shù)值模擬結(jié)果均滿足其強(qiáng)度條件，模擬數(shù)據(jù)可為其結(jié)構(gòu)改進(jìn)和同類戰(zhàn)斗部強(qiáng)度分析提供參考。

彈藥工程；發(fā)射強(qiáng)度；數(shù)值模擬

彈丸在發(fā)射過(guò)程中沿炮膛運(yùn)動(dòng)，在多種復(fù)雜載荷的共同作用下，會(huì)使戰(zhàn)斗部殼體材料產(chǎn)生應(yīng)力和變形。如果戰(zhàn)斗部的變形超過(guò)其許用限度，會(huì)影響戰(zhàn)斗部在膛內(nèi)的正常運(yùn)動(dòng)，甚至造成彈體破裂、發(fā)生膛炸。為保證戰(zhàn)斗部發(fā)射的安全性，在戰(zhàn)斗部結(jié)構(gòu)設(shè)計(jì)時(shí)均需要對(duì)其進(jìn)行強(qiáng)度分析。

前蘇聯(lián)學(xué)者曾提出一種確定彈帶壓力的方法，但由于不能直接計(jì)算而不被廣泛使用。布林克方法將彈體簡(jiǎn)化為無(wú)限長(zhǎng)厚壁圓筒[1]，能有效校核結(jié)構(gòu)相對(duì)簡(jiǎn)單的圓柱部及定心部區(qū)域。國(guó)內(nèi)學(xué)者也做了大量的研究工作，并給出相關(guān)的工程近似公式。但由于戰(zhàn)斗部結(jié)構(gòu)的不規(guī)則性和載荷的復(fù)雜性，且允許出現(xiàn)一定程度的塑性變形，使得彈丸強(qiáng)度分析大多停留在實(shí)彈射擊實(shí)驗(yàn)的基礎(chǔ)上，常規(guī)方法無(wú)法很好的解答[2]。

本文借助有限元軟件對(duì)戰(zhàn)斗部膛內(nèi)運(yùn)動(dòng)過(guò)程進(jìn)行分析研究，校核其結(jié)構(gòu)強(qiáng)度，并進(jìn)行了實(shí)驗(yàn)驗(yàn)證，為該彈的改進(jìn)和研發(fā)提供參考。

1　載荷及失效分析[1]

1.1載荷分析

彈丸在發(fā)射過(guò)程中載荷主要有：火藥氣體壓力、彈帶壓力、導(dǎo)轉(zhuǎn)側(cè)力、慣性力、裝填物壓力、摩擦力、不均衡力等。其中火藥氣體壓力是最主要載荷，是推動(dòng)彈丸在膛內(nèi)運(yùn)動(dòng)的主要因素，對(duì)于線膛炮所配用的旋轉(zhuǎn)穩(wěn)定彈，由于有彈帶的密封作用，火藥氣體幾乎完全作用在彈帶后部的彈尾部，對(duì)彈尾部的結(jié)構(gòu)強(qiáng)度有很大影響。其計(jì)算公式：

(1)

式中：P為炮膛壓力；φ為次要功系數(shù)；mω為發(fā)射藥質(zhì)量；m為彈丸質(zhì)量。

1.2臨界狀態(tài)

彈丸的發(fā)射過(guò)程將面臨三個(gè)危險(xiǎn)臨界狀態(tài)。

第一臨界狀態(tài)，即彈帶完全嵌入身管膛線時(shí)。彈帶區(qū)在徑向壓力的作用下徑向壓縮變形，除彈帶區(qū)域，其它區(qū)域的應(yīng)力、應(yīng)變值均較小。彈帶區(qū)變形如圖1所示。

圖1　第一臨界狀態(tài)彈體變形示意圖

第二臨界狀態(tài)，即膛內(nèi)壓力達(dá)到最大時(shí)。此時(shí)彈底壓力值、軸向過(guò)載達(dá)到峰值；彈底由于火藥氣體的作用產(chǎn)生凹陷變形，彈尾部和彈帶區(qū)產(chǎn)生徑向壓縮變形。在軸向慣性力作用下，彈丸頭部和圓柱部將產(chǎn)生軸向墩粗、徑向膨脹變形。彈體變形如圖2所示。

圖2　第二臨界狀態(tài)彈體變形示意圖

彈丸第三臨界狀態(tài)，即彈丸完全脫離炮口時(shí)。由于大部分載荷突然卸載，使彈體材料因彈性恢復(fù)而發(fā)生震顫。因此對(duì)于抗拉強(qiáng)度大大低于抗壓強(qiáng)度的脆性材料，必須考慮該狀態(tài)對(duì)彈體強(qiáng)度的影響。文中的殼體材料為35CrMnSiA鋼，有較好的韌性和抗拉性，故數(shù)值模擬時(shí)忽略該狀態(tài)對(duì)彈丸的影響。

1.3失效準(zhǔn)則

1)應(yīng)力準(zhǔn)則

強(qiáng)度條件為：σ<[σ0]

式中：σ為材料所受應(yīng)力；[σ0]為材料的許用應(yīng)力。

2)應(yīng)變準(zhǔn)則

強(qiáng)度條件為：2W*<[2W*]

式中：W*為材料的殘余變形量；[2W*]為允許的材料殘余變形量。

根據(jù)文獻(xiàn)[1]，彈體材料的殘余變形許用值如表1所示。

表1　彈體殘余變形許用量

2　數(shù)值模擬

2.1物理模型

戰(zhàn)斗部的結(jié)構(gòu)如圖3所示，由引信、頭螺、含能破片、彈丸殼體、推板、拋射藥、炸藥等組成?？紤]戰(zhàn)斗部零部件的高度軸對(duì)稱特性，實(shí)際計(jì)算時(shí)采用1/4模型，cm-g-μs單位制。

圖3　戰(zhàn)斗部結(jié)構(gòu)和火炮身管模型

2.2數(shù)值模型[3]

戰(zhàn)斗部各零部件及火炮身管均采用Lagrange算法。在戰(zhàn)斗部擠進(jìn)過(guò)程中，身管材料變形量很小，可以忽略不計(jì)，將其簡(jiǎn)化為剛性材料。參數(shù)分別為ρ=7.83g/cm3、E=210GPa、μ=0.3。

拋射藥和炸藥采用Elastic材料模型描述。彈丸引信、頭螺、含能破片殼體、含能材料、彈丸殼體和推板均采用Plastic-Kinematic材料模型，應(yīng)力應(yīng)變關(guān)系為

σy=[1+(ε/C)-q](σ0+βEεeff)

(2)

式中：ε為應(yīng)變率；εeff為有效塑性應(yīng)變；σ0、E、C為材料參數(shù)，表征材料敏感率特征。具體參數(shù)值見(jiàn)表2，未列出參數(shù)采用軟件自帶缺省值或忽略。

表2　各材料性能參數(shù)[4-5]

針對(duì)戰(zhàn)斗部的前后彈帶，考慮到其在隨彈丸一同擠進(jìn)身管膛線過(guò)程中會(huì)出現(xiàn)高強(qiáng)度、高應(yīng)變、高應(yīng)變率的特征，采用Johnson-Cook材料模型和Gruneisen狀態(tài)方程描述。Johnson-Cook模型的基本公式為

(3)

表3　彈帶材料性能參數(shù)[4-5]

2.3數(shù)值模擬結(jié)果及分析

圖4和圖5分別給出了上彈帶完全擠進(jìn)膛線和膛內(nèi)壓力最大時(shí)刻的等效應(yīng)力分布云圖，由于篇幅所限，其它時(shí)刻的等效應(yīng)力云圖不再列出。

圖4　應(yīng)力分布(彈帶完全擠入)

圖5　應(yīng)力分布(膛壓最大)

從計(jì)算結(jié)果來(lái)看，隨著彈帶不斷擠進(jìn)膛線，膛線對(duì)彈帶的壓力不斷增大，彈帶凹槽區(qū)域的應(yīng)力隨之增大，應(yīng)力波由彈帶凹槽處向兩側(cè)傳播，當(dāng)上彈帶完全擠進(jìn)膛線時(shí)，最大應(yīng)力可達(dá)到1382MPa，出現(xiàn)在彈帶凹槽處頂端。與此同時(shí)，彈帶與彈體的塑性應(yīng)變也不斷增大，彈帶最大塑性應(yīng)變?yōu)?.284mm，出現(xiàn)在彈帶頂部。

1931年9月，但采爾如約來(lái)到南京的中研院，帶來(lái)民族學(xué)標(biāo)本134種，另外還有圖標(biāo)120幅、各族照片28幅。民族博物館建設(shè)，因經(jīng)費(fèi)等各種原因耽擱下來(lái)，最終未能實(shí)現(xiàn)。為充分利用但采爾帶來(lái)的這些物品，蔡元培就專門(mén)指示在中研院辟出地方，設(shè)立了一個(gè)民族學(xué)標(biāo)本陳列室。遠(yuǎn)道而來(lái)的但采爾便留在中研院民族學(xué)組從事民族學(xué)研究。

隨著彈體的不斷運(yùn)動(dòng)，上、下彈帶全部進(jìn)入膛線，由擠壓彈帶而引起的徑向應(yīng)力開(kāi)始減弱，彈丸加速度逐漸增大，由火藥推力引起的軸向慣性力開(kāi)始占據(jù)主導(dǎo)，在2700～3000μs時(shí)間內(nèi)達(dá)到最大。彈丸頂部的含能破片及推板與彈體的接觸應(yīng)力不斷增加，含能破片區(qū)域有明顯的傳向彈底的應(yīng)力波，最大應(yīng)力為481MPa，出現(xiàn)在梅花形推板和彈體接觸區(qū)域。而彈體由于出現(xiàn)應(yīng)力集中點(diǎn)，最大應(yīng)力達(dá)到995MPa，除應(yīng)力集中點(diǎn)外其它主要受力部位應(yīng)力在400～600MPa之間，最大應(yīng)力出現(xiàn)在彈底靠近彈腔區(qū)域。

由于該戰(zhàn)斗部結(jié)構(gòu)較為復(fù)雜，為更精確分析殼體強(qiáng)度，了解各區(qū)域的應(yīng)力、應(yīng)變分布情況，在彈丸殼體壁厚突變和應(yīng)力集中處選取7個(gè)觀察點(diǎn)，如圖6所示，并提取發(fā)射過(guò)程中相應(yīng)單元的應(yīng)力、應(yīng)變時(shí)間歷程曲線，如圖7和圖8所示。

圖6　戰(zhàn)斗部殼體觀察點(diǎn)單元

圖7　應(yīng)力曲線

圖8　塑性應(yīng)變曲線

由圖7和圖8可以看出，900～1200μs彈丸上、下彈帶相繼擠進(jìn)膛線，面臨第一臨界狀態(tài)。彈帶凹槽中間處單元E應(yīng)力出現(xiàn)輕微波動(dòng)，其最大應(yīng)力遠(yuǎn)小于材料屈服強(qiáng)度(1275MPa)，未發(fā)生塑性變形。單元C、D應(yīng)力值出現(xiàn)劇烈波動(dòng)，其峰值高于殼體材料屈服強(qiáng)度，分別出現(xiàn)了0.0178mm和0.0032mm的輕微變形。在彈丸實(shí)際加工過(guò)程中，各棱角處會(huì)采用工藝倒角以避免此類因應(yīng)力集中現(xiàn)象而導(dǎo)致的變形，因此可以忽略其變形，即彈丸未發(fā)生變形。

2850μs時(shí)，膛壓值達(dá)到最大，面臨第二臨界狀態(tài)。殼體壁厚突變處單元A、B和彈底拐角及彈底中心處單元F、G的應(yīng)力值均隨彈丸軸向慣性力的增大而增大，并達(dá)到各自峰值，但以上單元的應(yīng)力峰值均低于其屈服強(qiáng)度，處于線彈性階段。

3　實(shí)驗(yàn)驗(yàn)證

為驗(yàn)證數(shù)值模擬結(jié)果正確性，進(jìn)行了彈丸強(qiáng)度實(shí)驗(yàn)。彈丸配阻力引信，其結(jié)構(gòu)如圖9所示。

圖9　實(shí)驗(yàn)彈丸

以339.7MPa膛壓，采用某型彈道炮進(jìn)行強(qiáng)裝藥試射，實(shí)驗(yàn)射擊數(shù)量共2發(fā)。選定彈體上5個(gè)特殊位置，測(cè)量彈丸各位置發(fā)射前后的尺寸，結(jié)果如表4所示。

表4　射擊前后尺寸　mm

由表4可以看出，兩發(fā)彈丸在射擊前后各部位變形量均為0，即彈丸在射擊過(guò)程中殼體結(jié)構(gòu)未發(fā)生塑性變形，與模擬仿真結(jié)果相吻合。

4　結(jié)論

(1)彈丸在膛內(nèi)運(yùn)動(dòng)過(guò)程中，第一臨界狀態(tài)時(shí)，由于并未對(duì)幾何模型進(jìn)行倒角，致使彈帶凹槽頂部出現(xiàn)應(yīng)力集中現(xiàn)象，局部應(yīng)力達(dá)到1382MPa，超過(guò)材料屈服應(yīng)力1275MPa，排除該因素影響，彈丸各部位均未發(fā)生變形。第二臨界狀態(tài)時(shí)彈體最大應(yīng)力出現(xiàn)在彈底及靠近彈腔附近區(qū)域，其最大應(yīng)力值為573MPa，小于屈服應(yīng)力，處于彈性變形階段。

(2)實(shí)驗(yàn)結(jié)果表明彈體各部位均未發(fā)生塑性變形，模擬所得相關(guān)數(shù)據(jù)符合實(shí)際情況，證明該仿真模型較為合理，可以作為同類彈丸設(shè)計(jì)參考。

[1]魏惠之，朱鶴松，汪東暉，等.彈丸設(shè)計(jì)理論[M].北京：國(guó)防工業(yè)出版社，1985：43-47.

[2]卜仕彪.復(fù)合增程榴彈第一臨界狀態(tài)動(dòng)力學(xué)模型計(jì)算分析[D].南京：南京理工大學(xué)，2006：4-11.

[3]李裕春，時(shí)黨勇.ANSYS10.0/LS-DYNA基礎(chǔ)理論與工程實(shí)踐[M].北京：中國(guó)水利水電出版社，2006.

[4]米糧川，高樹(shù)滋，賴長(zhǎng)纓.內(nèi)彈道擠進(jìn)過(guò)程仿真研究[J].彈箭與制導(dǎo)學(xué)報(bào)，2012,32(1)：133-135.

[5]趙曉寧，何勇，張先鋒，等.A3鋼抗高速桿彈侵徹的數(shù)值模擬與實(shí)驗(yàn)研究[J].北京理工大學(xué)學(xué)報(bào)，2011，35(2)：165-166.

(責(zé)任編輯：趙麗琴)

The Numerical Simulation of Launching Intensity for Rifled Gun Warhead

XU Heyang，QIAO Xiangxin，GUO Keqiang，LI Wenhao，PENG Linxi

(Shenyang Ligong University，Shenyang 110159，China)

To study the launching intensity of a new warhead,LS-DYNA software is used to develop the motion model of warhead in the bore,using the software to make the launching intensity simulation and analysis for the model to get the stress and strain curve of warhead in the different stages of exercise,and to analyze the stress and strain data when the warhead is in the first and the second critical state and the corresponding value of the strength criterion,and to have a test.Simulation results with actual firing test results are basically consistent.The results show that the results of numerical simulation and experiment can meet the requirements of practical application,and the simulation data can provide reference for the structural improvement and the strength analysis of the same kind of warhead.

ammunition engineering；launching intensity；numerical simulation

2015-12-29

徐赫陽(yáng)(1992—)，男，碩士研究生；通訊作者：?jiǎn)滔嘈?1959—)，男，教授，研究方向：彈藥工程與爆炸技術(shù)。

TJ430

A