陸野，周克棟，赫雷，李峻松，黃雪鷹

（1.南京理工大學(xué)機械工程學(xué)院，江蘇南京210094；2.中國兵器工業(yè)第208研究所，北京102202；3.63856部隊，吉林白城137001）

坡膛結(jié)構(gòu)參數(shù)對槍械內(nèi)彈道擠進時期的影響研究

陸野1，周克棟1，赫雷1，李峻松2，黃雪鷹3

（1.南京理工大學(xué)機械工程學(xué)院，江蘇南京210094；2.中國兵器工業(yè)第208研究所，北京102202；3.63856部隊，吉林白城137001）

槍械射擊過程中槍管坡膛處工作環(huán)境惡劣，為了研究槍管坡膛角度對擠進過程坡膛處受力的影響，建立了不同坡膛角下考慮槍管及彈頭結(jié)構(gòu)特性、本構(gòu)非線性等因素的三維有限元模型，分析了不同坡膛角對彈頭擠進過程的影響，獲得了不同坡膛角下擠進阻力隨擠進位移的變化數(shù)據(jù)；建立了擠進阻力的響應(yīng)面模型，基于上述數(shù)據(jù)采用Hermite多項式，求解獲得了擠進阻力以坡膛角和擠進位移為變量的計算公式；提出建立了考慮擠進阻力的彈頭擠進過程的動力學(xué)模型，編程計算了槍彈擠進過程中的擠進壓力，獲得了某大口徑機槍滿足彈頭初速條件的坡膛角度取值范圍為0.11°～1.13°，進而得到了緩減槍管坡膛受力、保證彈頭最高初速、滿足坡膛角設(shè)計范圍的坡膛角度最優(yōu)解為0.56°。

兵器科學(xué)與技術(shù)；槍管；彈頭擠進；坡膛；接觸；沖擊

0 引言

自動武器在連發(fā)射擊過程中，槍管經(jīng)歷著火藥燃氣高溫、高壓、彈頭高速、高加速度的高頻熱壓耦合沖擊、機械沖擊等作用，這種工況嚴重影響了槍管壽命。實踐表明，槍管坡膛處工作環(huán)境最惡劣，最容易先發(fā)生破壞［1］。目前，我國自動武器壽命偏低已成為制約我國新一代槍械性能提高的瓶頸問題。實際工況下的自動武器彈頭擠進以及內(nèi)彈道過程很難通過實驗來得到準(zhǔn)確的數(shù)據(jù)，而彈頭擠進過程又是影響槍管壽命的主要因素之一。目前，關(guān)于自動武器彈頭擠進過程的國內(nèi)外的研究報道還不是很多。文獻［2-4］對彈頭擠進過程進行了初步的分析與研究，但均未涉及彈頭擠進過程中膛線壓痕的形成過程及變形特點。文獻［5］建立了彈頭擠進槍管的有限元分析模型，通過數(shù)值模擬研究了鉛芯彈頭的擠進過程，展示了彈頭上膛線壓痕的形成過程和材料的流動情況，分析了擠進前后彈頭殼和鉛芯的變形特征，以及擠進結(jié)束后彈頭的殘余應(yīng)力。文獻［6］建立了彈帶材料基于連續(xù)介質(zhì)損傷力學(xué)的本構(gòu)模型，同時考慮了經(jīng)典內(nèi)彈道方程組和彈帶擠進過程的耦合效應(yīng)，對兩種坡膛結(jié)構(gòu)下彈帶擠進過程進行了數(shù)值模擬。上述研究均進行了三維數(shù)值仿真分析，但并沒有考慮彈頭實際內(nèi)部結(jié)構(gòu)對仿真結(jié)果的影響，同時也沒有對擠進阻力、擠進壓力等擠進參數(shù)隨槍管內(nèi)膛結(jié)構(gòu)參數(shù)變化的規(guī)律進行相應(yīng)的分析。本文針對某大口徑機槍的擠進內(nèi)彈道時期，構(gòu)建了槍彈擠進過程中的高沖擊、強摩擦和大變形接觸下的有限元模型，模擬了彈頭的擠進過程，分析了6種不同坡膛結(jié)構(gòu)對彈頭擠進過程的影響，解算了擠進阻力以坡膛角和擠進位移為變量的計算公式；通過建立考慮擠進阻力的彈頭擠進過程動力學(xué)模型，計算了變?nèi)轄顟B(tài)下槍彈擠進過程中的膛內(nèi)火藥燃氣平均壓力，分析了彈頭初速隨擠進壓力、擠進壓力隨坡膛角的變化規(guī)律，獲得了滿足彈頭初速條件的擠進壓力與坡膛角度取值范圍，進而得到了緩減槍管坡膛受力、保證彈頭最高初速、滿足坡膛角設(shè)計范圍的坡膛角度最優(yōu)解。研究結(jié)果對提高槍管壽命具有一定的指導(dǎo)意義。

1 槍彈擠進阻力理論分析

槍彈擠進過程中，彈頭受到槍管坡膛、陰線和導(dǎo)轉(zhuǎn)側(cè)的作用力，分別為坡膛阻力Fc，陰線阻力Fs和導(dǎo)轉(zhuǎn)阻力Fd，這三部分的合力構(gòu)成了擠進阻力F.對擠進過程彈頭受力分析如圖1所示，其中η為纏角，β為坡膛角，α為導(dǎo)轉(zhuǎn)側(cè)角度。

令在接觸面上產(chǎn)生的接觸應(yīng)力分別為:σc、子c，σs、子s和σd、子d，其中法向應(yīng)力和切向應(yīng)力通過摩擦系數(shù)υ呈比例關(guān)系，即

圖1 擠進過程彈頭受力分析圖Fig.1 Mechanical analysis of projectile during engraving

由此可得沿彈軸方向的擠進阻力分別為

彈頭擠進膛線的過程時間極短，擠進結(jié)束后彈頭所獲得的速度非常小，彈頭擠進過程彈頭和坡膛的摩擦特性變化很小，可認為摩擦系數(shù)υ不變，由上述理論分析可知，在其余參數(shù)一定的情況下，擠進阻力僅與坡膛角和擠進位移有關(guān)。為了更加直觀地表述擠進阻力與坡膛角和擠進位移的數(shù)值關(guān)系，本文通過高精度有限元仿真獲得擠進阻力與坡膛角的對應(yīng)數(shù)據(jù)，采用隨機響應(yīng)面法建立擠進阻力與坡膛角及擠進位移的代理模型，計算擠進阻力隨坡膛角及擠進位移變化的函數(shù)表達式。

2 擠進過程有限元分析

2.1 有限元模型的建立

本文以某大口徑機槍槍管及彈頭為對象，槍管壁厚、膛線，彈頭內(nèi)部等保持實際結(jié)構(gòu)，如圖2所示，考慮彈頭的柔性效應(yīng)，采用有限元前處理軟件Hypermesh分別對槍管、彈頭主體和彈帶劃分網(wǎng)格，并對關(guān)鍵部位進行精細化處理，如圖3所示。

圖2 槍管與彈頭的三維實體模型Fig.2 3D solid model of barrel and projectile

圖3 槍管與彈頭的有限元模型Fig.3 Finite element models of barrel and projectile

2.2 加載與求解

由擊發(fā)底火到彈頭全部嵌入膛線的過程稱為擠進膛線時期，此過程由于彈頭移動量不大，可近似認為火藥定容燃燒，引用火藥氣體定容燃燒狀態(tài)方程為

式中:f為火藥力；Δ為裝填密度；ψ為火藥已燃燒的相對質(zhì)量；δm為火藥密度；α為余容；pB為點火藥壓力。

引用燃氣生成方程與燃燒速度方程為

式中:z為火藥已燃燒的相對厚度；n為燃速指數(shù)；u1為燃速系數(shù)；e1為火藥原厚度；χ、λ、μ均為藥形系數(shù)。

聯(lián)立方程（8）式、（9）式，代入某大口徑槍彈的火藥特性參數(shù)，通過MATLAB編程即可計算出擠進膛線時期的彈底壓力曲線施加于彈頭尾端面，即為力邊界條件。設(shè)置槍管楔形凹槽處為固定約束，該條件為位移邊界條件。

本有限元模型應(yīng)用Johoson-Cook材料本構(gòu)關(guān)系［7］，采用適用于求解復(fù)雜的非線性、大變形的顯式動力學(xué)問題的顯式非線性求解器LS-DYNA進行計算，針對擠進過程的特征，采用自動接觸模型算法，設(shè)置罰函數(shù)因子保持接觸界面的協(xié)調(diào)性；考慮到彈頭擠進過程中可能出現(xiàn)網(wǎng)格畸變現(xiàn)象，在計算過程中采用ALE自適應(yīng)網(wǎng)格技術(shù)，確保計算的精度及效率。

為確保彈頭擠進過程的響應(yīng)面代理模型的準(zhǔn)確性，本文基于現(xiàn)有槍管坡膛角的設(shè)計范圍，建立了6種不同坡膛角（如表1所示）的有限元模型，分別進行了仿真分析。針對彈頭被甲的變形與膛線刻痕的形成情況和擠進過程擠進阻力的變化情況，使用后處理軟件Ls-Prepost，對仿真計算的數(shù)據(jù)進行整理分析。

表1 6種不同方案坡膛角度列表Tab.1 Forcing cone angles of six different schemes

3 槍彈擠進過程代理模型的建立

本文采用隨機響應(yīng)面法構(gòu)建彈頭擠進過程的代理模型［8］。目前常用的數(shù)據(jù)擬合方法有Gauss法、Sina法、Lorentz法、Voigt法、Hermite法等，其中前面4種方法適用于單變量數(shù)據(jù)的擬合，而Hermite多項式法可用于多個變量數(shù)據(jù)的擬合，擬合精度高。本文基于有限元仿真獲得了不同坡膛角及擠進位移兩個變量的擠進阻力的變化數(shù)據(jù)，故采用Hermite多項式編寫MATLAB程序求解擠進阻力以坡膛角和擠進位移為變量的計算公式，分析擠進阻力的變化規(guī)律。輸出響應(yīng)量H采用如下Hermite多項式展開來表示:

式中:a0、ai1、ai1i2等為待定系數(shù)；xi1，xi2，…，xin為獨立的變量；q為變量的個數(shù)；φi為i階Hermite多項式。本文中輸出響應(yīng)量為擠進阻力F，獨立變量為坡膛角β和擠進位移x.

確定（10）式的主要任務(wù)是求解待定系數(shù)ai1，目前求解待定系數(shù)的方法很多，例如配點法、高效配點法和回歸法等。Isukkapalli等認為回歸法穩(wěn)健性較好［9］，本文采用逐步回歸的方式得到彈頭擠進過程的響應(yīng)面代理模型。

Hermite多項式的遞推關(guān)系為

4 擠進過程彈頭動力學(xué)方程的建立

4.1 擠進壓力的影響分析

擠進壓力是指彈頭全部嵌入膛線受到最大阻力時所對應(yīng)的膛內(nèi)火藥燃氣壓力。由于坡膛角的變化會引起擠進過程阻力、時間等參數(shù)發(fā)生變化，而火藥的燃燒對時間極其敏感，從而影響膛內(nèi)火藥燃氣壓力，使得整個內(nèi)彈道性能發(fā)生變化，因此為保證彈頭的初速性能，需要保證擠進壓力和坡膛角在一定范圍內(nèi)變化。為有效分析槍彈擠進過程的影響，并提出保證彈頭的初速性能，降低擠進阻力峰值的可行方案，本節(jié)在隨機響應(yīng)面法建立的槍彈擠進過程的代理模型基礎(chǔ)上，結(jié)合火藥變?nèi)萑紵^程和彈頭起始動力學(xué)過程來分析槍彈擠進動力學(xué)過程，建立了考慮擠進阻力的彈頭擠進過程動力學(xué)模型。

4.2 擠進過程動力學(xué)方程的建立

由于經(jīng)典內(nèi)彈道方程忽略了彈頭擠進膛線的過程，假定彈頭全部擠進膛線達到擠進壓力時，彈頭才開始運動，顯然誤差較大。本文通過將建立的擠進阻力代理模型與考慮了火藥氣體變?nèi)萑紵蛷楊^運動過程的經(jīng)典內(nèi)彈道方程組耦合，建立了考慮擠進阻力的彈頭擠進過程動力學(xué)模型，以此來仿真計算膛內(nèi)火藥燃氣壓力。新建立的槍彈擠進過程動力學(xué)方程如下:

式中:S為身管的等效截面積；u1為正比燃燒定律中的燃速系數(shù)；m為彈頭質(zhì)量；p為膛內(nèi)火藥燃氣壓力；F由槍彈擠進過程的代理模型求解確定；l為彈頭行程；lψ為藥室自由容積縮徑長；f為火藥力；ω為裝藥量；θ=k-1，k為比熱比；φ為忽略了擠進阻力所做功的次要功系數(shù)，本文中將擠進阻力所做功進行了詳細求解，而不是以往采用系數(shù)K3進行粗略的估算；K2、K4、K5分別為彈頭旋轉(zhuǎn)功、火藥氣體的運動功和武器后坐部分運動功占主要功的百分比。

5 算例

本文以某大口徑機槍的擠進過程為例，進行了建模、仿真與分析。

5.1 彈頭被甲膛線刻痕的形成過程分析

為更加直觀地反應(yīng)彈頭被甲膛線刻痕的形成過程，取膛線刻痕形成的階段特征圖，按時間順序展示，如圖4所示，隨著彈頭擠進槍管，彈頭被甲膛線刻痕由彈頭向彈尾延展，逐漸加深成形。不同坡膛結(jié)構(gòu)下刻痕形成過程、形狀基本一致。

圖4 膛線刻痕的形成過程Fig.4 Formation process of rifling indentation

5.2 擠進過程擠進阻力分析

通過對6種不同坡膛角時擠進過程的仿真結(jié)果中每個單元所受的阻力進行合成得到擠進阻力隨時間及空間變化的規(guī)律，如圖5、圖6所示。由擠進阻力隨時間變化圖圖5可知，隨著坡膛角度變大，擠進阻力達到的峰值越大；對所有坡膛角，擠進阻力在達到峰值后逐漸下降并趨于穩(wěn)定，這是由于擠進阻力包括克服彈頭材料變形的力和軸向摩擦力，擠進結(jié)束后克服彈頭材料變形的力消失，擠進阻力趨于的穩(wěn)定值為軸向摩擦力；隨著坡膛角度變大，擠進阻力達到峰值的時刻越靠前，這是由于隨著坡膛角度變大，擠進阻力中克服彈頭材料變形的力增大的速率變快，因此擠進阻力達到峰值的時間隨坡膛角增大而減少。由擠進阻力隨擠進位移變化圖圖6可知，擠進阻力的峰值隨著坡膛角度的增加而增加，彈頭與槍管膛線接觸后擠進阻力迅速增加，由于彈頭中間部位設(shè)有階梯凹槽，因此擠進阻力存在短暫下降階段。

圖5 不同坡膛角時擠進阻力隨時間變化圖Fig.5 Engraving resistance vs.time at different angles of forcing cone

圖6 不同坡膛角時擠進阻力隨擠進位移變化圖Fig.6 Engraving resistance vs.engraving displacement at different angles of forcing cone

基于隨機響應(yīng)面法構(gòu)建的彈頭擠進過程的代理模型與6組不同坡膛角對應(yīng)的擠進阻力隨位移的變化數(shù)據(jù)，采用MATLAB編寫計算程序，獲得擠進阻力以坡膛角和擠進位移為變量的計算公式如下:

矩陣A為常數(shù)矩陣，采用逐步回歸法通過MATLAB程序求解Hermite多項式中a0、ai1、ai1i2等待定系數(shù)獲得。

矩陣X中:X0=1；X1=beta（1）；X2=m（1）；X3=beta（2）；X4=m（2）；X5=beta（1）×m（1）；X6= beta（3）；X7=m（3）；X8=beta（2）×m（1）；X9= beta（1）×m（2）；X10=beta（3）×m（1）；X11=beta（2）× m（2）；X12=beta（1）×m（3）；X13=m（4）；X14= beta（3）×m（2）；X15=beta（2）×m（3）；X16=beta（1）× m（4）；X17=m（5）；X18=beta（3）×m（3）；X19= beta（2）×m（4）；X20=beta（1）×m（5）；X21=m（6）；X22=beta（3）×m（4）；X23=beta（2）×m（5）；X24= beta（1）×m（6）；X25=m（7）；X26=beta（3）×m（5）；X27=beta（2）×m（6）；X28=beta（1）×m（7）；X29= beta（3）×m（6）；X30=beta（2）×m（7）；X31= beta（3）×m（7）.其中:m（0）=1；m（1）=2x；m（2）=4x2-2；m（3）=8x3-12x；m（4）=16x4-48x2+12；m（5）=32x5-160x3+120x；m（6）= 2x（32x5-160x3+120x）-10×（16x4-48x2+12）；m（7）=2x［2x（32x5-160x3+120x）-10×（16x4-48x2+12）］-12×（32x5-160x3+120x）；beta（1）=2β；beta（2）=4β2-2；beta（3）=8β3-12β.

5.3 計算與分析

根據(jù)建立的擠進阻力代理模型與考慮了擠進阻力的擠進過程動力學(xué)方程，計算了6種不同坡膛角下膛內(nèi)火藥燃氣壓力隨擠進位移變化關(guān)系，如圖7所示，結(jié)果顯示，膛內(nèi)火藥燃氣壓力隨著擠進位移增加而單調(diào)增加，擠進位移為25 mm時彈頭完全嵌入膛線，此時的膛內(nèi)火藥燃氣壓力即擠進壓力，隨坡膛角的增加而變小。這是由于隨著坡膛角度的增加，相同的彈頭刻痕長度對應(yīng)的擠進時間減小，火藥燃燒時間減小導(dǎo)致壓力變小。

圖7 膛內(nèi)火藥燃氣壓力隨擠進位移的變化曲線Fig.7 Powder gas pressure vs.engraving displacement

基于擠進過程動力學(xué)方程，計算現(xiàn)有槍管長度下，彈頭出膛口時的速度隨坡膛角的變化曲線，如圖8所示。彈頭出膛口時的速度隨著坡膛角的變化先增加后降低，這是由于坡膛角過小，會引起擠進時間過長，火藥燃燒過于充分，坡膛角過大，會引起擠進時間過短，導(dǎo)致火藥燃燒不充分，兩種情況均會引起彈頭在內(nèi)彈道時期所受的彈后平均壓力降低，導(dǎo)致彈頭出膛口時的速度下降。由于某大口徑槍彈技術(shù)指標(biāo)要求初速大于810.0 m/s，在保證彈頭初速的前提下，坡膛角的取值范圍為0.11°～1.13°.

由于在陰線深度一定的條件下，坡膛角度越小，槍管坡膛五錐長度就越長，考慮到實際加工的可能性，坡膛五錐長度應(yīng)在一定的合理范圍內(nèi)，即滿足表1中所提出的槍管坡膛角的設(shè)計范圍。因此坡膛角最優(yōu)解的選擇應(yīng)滿足3個方面的條件，即相對于原坡膛角度擠進壓力減小、彈頭初速最高、在槍管坡膛角的設(shè)計范圍內(nèi)。

本文所研究槍管的實際坡膛角度為0.85°，由圖8可知，最高初速對應(yīng)的坡膛角為0.53°，根據(jù)表1中有槍管坡膛角的設(shè)計范圍，選取槍管的最佳坡膛角為0.56°；將實際坡膛角與本文所得最佳坡膛角0.56°分別采用有限元仿真方法分析和代入代理模型計算，獲得仿真與代理模型分別計算所得擠進阻力隨位移變化對比圖如圖9所示，顯然兩條仿真曲線均與代理模型計算得到的擠進阻力曲線具有較高的一致性，驗證了代理模型的有效性，同時坡膛角為0.56°時，擠進阻力較小，坡膛受力較為緩和。

圖8 彈頭出膛口時的速度隨坡膛角的變化曲線Fig.8 Muzzle velocity vs.forcing cone angle

圖9 擠進阻力隨擠進位移變化對比圖Fig.9 Engraving resistance vs.engraving displacement

6 結(jié)論

本文針對某大口徑槍械內(nèi)彈道擠進時期，進行了理論分析與數(shù)值仿真計算，得出以下結(jié)論:

1）對槍彈擠進過程的理論分析表明在槍管陽線寬度、纏角和導(dǎo)轉(zhuǎn)側(cè)角度等參數(shù)以及槍彈材料一定的情況下，擠進阻力僅與坡膛角和擠進位移有關(guān)。

2）通過建立不同坡膛角的有限元模型并進行仿真分析表明隨著坡膛角度變大，擠進阻力達到的峰值越大，擠進阻力達到峰值的時刻越靠前，在達到峰值后擠進阻力有所下降并趨于穩(wěn)定。

3）基于隨機響應(yīng)面法建立了擠進阻力與坡膛角及擠進位移的代理模型，數(shù)值解算得到了擠進阻力與坡膛角及擠進位移的函數(shù)表達式。

4）以經(jīng)典內(nèi)彈道方程組為基礎(chǔ)，提出了考慮擠進阻力的槍彈擠進過程動力學(xué)方程。

5）彈頭出膛口時的速度隨著坡膛角的變化先增加后降低，對本文研究的某大口徑機槍滿足彈頭初速條件的坡膛角度取值范圍為0.11°～1.13°；在槍管坡膛角6種設(shè)計方案中滿足彈頭初速最高和擠進壓力較小的最佳坡膛角為0.56°.

（

）

［1］ 歐學(xué)炳，殷仁龍，王學(xué)顏.自動武器結(jié)構(gòu)設(shè)計［M］.北京:北京理工大學(xué)出版社，1995. OU Xue-bing，YIN Ren-long，WANG Xue-yan.Structure design of automatic weapons［M］.Beijing:Beijing Institute of Technology Press，1995.（in Chinese）

［2］ 丘爾巴諾夫E B.擠進時期內(nèi)彈道學(xué)與擠進壓力計算［M］.楊敬榮，譯.北京:國防工業(yè)出版社，1997. Qierbarov E B.Interior ballistics and engraving force calculation during engraving of projectile［M］.YANG Jing-rong，translated. Beijing:National Defense Industry Press，1997.（in Chinese）

［3］ 段吉員，王彥平.彈丸擠進壓力對平衡炮內(nèi)彈道性能的影響研究［J］.高壓物理學(xué)報，2006，20（2）:189-193. DUAN Ji-Yuan，WANG Yan-ping，et al.The effect of engraving force of projectile on interior ballistic performance of DAVIS Gun［J］.Chinese Journal of High Pressure Physics，2006，20（2）: 189-193.（in Chinese）

［4］ Joseph T S，Dipak K，Michael M.Small caliber modeling from design to manufacture to launch［C］∥23th International Symposium on Ballistics.Tarragona，Spain:International Ballistics Committee，2007:557-564.

［5］ 樊黎霞，何湘玥.彈丸擠進過程的有限元模擬與分析［J］.兵工學(xué)報，2011，32（8）:963-969. FAN Li-xia，HE Xiang-yue.Finite element simulation and process analysis of projectile entering into barrel［J］.Acta Armamentarii，2011，32（8）:963-969.（in Chinese）

［6］ 孫河洋，馬吉勝，李偉，等.坡膛結(jié)構(gòu)變化對火炮內(nèi)彈道性能影響的研究［J］.兵工學(xué)報，2012，33（6）:669-675. SUN He-yang，MA Ji-sheng，LI Wei，et al.Study on influence of bore structure on gun's interior ballistic performances［J］.Acta Armamentarii，2012，33（6）:669-675.（in Chinese）

［7］ Johnson G，Cook W.Fracture characteristics of three metals subjected to various strains，strain rates，temperatures and pressures［J］.Engineering Fracture Mechanics，1985，21:31-48.

［8］ Basaga H B，Kaymaz I，Bayraktar A.An improved response surface method for reliability analysis of structures［J］.Structural Engineering and Mechanics，2012，42（4）:175-189.

［9］ Isukapalli S S，Ray A.Stochastic response surface methods（SRSMs）for uncertainty propagation:application to environmental and biological systems［J］.Risk Analysis，1998，18（3）:351-363.

Influence of Structure Parameters of Forcing Cone on Small Arms Interior Ballistics During Engraving

LU Ye1，ZHOU Ke-dong1，HE Lei1，LI Jun-song2，HUANG Xue-ying3
（1.School of Mechanical Engineering，Nanjing University of Science and Technology，Nanjing 210094，Jiangsu，China；2.No.208 Research Institute of China Ordnance Industries，Beijing 102202，China；3.Unit 63856 of PLA，Baicheng 137001，Jilin，China）

The working condition of barrel forcing cone is worse during the firing process of small arms. To study the influence of the forcing cone angles on the stress of barrel forcing cone，the three-dimensional finite element models of barrel and projectile at different forcing cone angles are established，in which the structures of barrel and projectile，and nonlinear constitutive relations are considered.The effects of different forcing cone angles on the projectile engraving process are analyzed.The data of engraving resistance is obtained under the condition of different forcing cone angles.A response surface model of engraving resistance is established.The formula and the variation rules of projectile engraving resistance along with the forcing cone angle and engraving displacement are calculated based on the above data by using Hermite polynomial.A dynamic model of projectile engraving process is established in consideration of the engraving resistance.The engraving pressure during the projectile engraving process is solved by programming.The forcing cone angles among 0.11°and 1.13°，which satisfy the requirement of muzzlevelocity of the projectile，are obtained，and then the best forcing cone angle of 0.56°，which could reduce the barrel stress to ensure the highest velocity and meet the design range of forcing cone angles in case of satisfying the interior ballistic behaviors，is got.

ordnance science and technology；barrel；projectile engraving；forcing cone；contact；impact

TJ303.9

A

1000-1093（2015）07-1363-07

10.3969/j.issn.1000-1093.2015.07.028

2014-09-18

陸野（1988—），男，博士研究生。E-mail:luye_njust@163.com；周克棟（1964—），男，教授，博士生導(dǎo)師。E-mail:zkd81151@njust.edu.cn