質(zhì)偏子彈裝配方位對(duì)母彈徑向質(zhì)偏的影響*

蔡泓杰，聞　泉，王雨時(shí)，張志彪（南京理工大學(xué)機(jī)械工程學(xué)院，南京　210094）

質(zhì)偏子彈裝配方位對(duì)母彈徑向質(zhì)偏的影響*

蔡泓杰，聞泉，王雨時(shí)，張志彪
（南京理工大學(xué)機(jī)械工程學(xué)院，南京210094）

針對(duì)質(zhì)偏子彈尺寸變化和裝配方位對(duì)母彈徑向質(zhì)偏的影響問(wèn)題，先通過(guò)Solidworks軟件得到子彈尺寸變化造成的子彈極限質(zhì)偏，再通過(guò)蒙特卡羅法模擬質(zhì)偏子彈裝配方位對(duì)母彈徑向質(zhì)偏的影響，得到不同質(zhì)量、質(zhì)心分布規(guī)律子彈不同裝配方位下裝入母彈后的母彈質(zhì)心位置分布。為了減小母彈徑向質(zhì)偏，子彈固定方位裝配時(shí)，應(yīng)盡量減小子彈徑向質(zhì)偏量；而子彈隨機(jī)方位裝配時(shí)，子彈質(zhì)偏量不需要特殊控制。

子母彈；蒙特卡羅法；設(shè)計(jì)原則；質(zhì)偏；裝配方位

0　引言

母彈徑向質(zhì)偏會(huì)影響子母彈外彈道性能。母彈彈道變化將引起子彈拋射點(diǎn)諸元變化，進(jìn)而影響子彈分布和子母彈殺傷威力［1］。

母彈質(zhì)心計(jì)算時(shí)，一般將子彈質(zhì)心位置和質(zhì)量假設(shè)為恒定值。但實(shí)際上每一枚子彈質(zhì)心位置在一定范圍內(nèi)都是隨機(jī)變化的，其質(zhì)量也如此，這就會(huì)對(duì)母彈質(zhì)心位置造成影響，因此有必要研究母彈徑向質(zhì)偏。計(jì)算機(jī)模擬為研究該問(wèn)題提供了便捷途徑。文獻(xiàn)［2］以某導(dǎo)彈戰(zhàn)斗部部件為例分析了尺寸誤差對(duì)質(zhì)心偏移的影響，對(duì)比分析了采用極值法與蒙特卡羅法計(jì)算質(zhì)偏的結(jié)果，表明采用蒙特卡羅法更符合實(shí)際。文獻(xiàn)［3］分析了影響母彈質(zhì)心徑向位置的因素，包括子彈串徑向質(zhì)心分布、母彈結(jié)構(gòu)件尺寸變化、母彈及母彈結(jié)構(gòu)件同軸度、母彈內(nèi)各連接螺紋配合間隙等，并模擬了子母彈徑向質(zhì)心位置分布，結(jié)果表明質(zhì)心位置分布范圍很小但不能視為零。

現(xiàn)有母彈徑向質(zhì)心位置分析都未將子彈尺寸變化引起的徑向質(zhì)心變化考慮在內(nèi)。文中針對(duì)此問(wèn)題，對(duì)子彈尺寸變化造成的子彈徑向和軸向質(zhì)偏以及質(zhì)偏子彈裝配方位變化對(duì)母彈徑向質(zhì)偏的影響進(jìn)行了研究。

1　蒙特卡羅方法簡(jiǎn)介

蒙特卡羅方法又稱統(tǒng)計(jì)模擬法，是一種使用隨機(jī)數(shù)（或更常見(jiàn)的偽隨機(jī)數(shù)）來(lái)解決問(wèn)題的一種數(shù)值計(jì)算方法。其實(shí)質(zhì)是通過(guò)對(duì)建立的數(shù)學(xué)模型進(jìn)行大量隨機(jī)試驗(yàn)，利用概率論求得原始問(wèn)題的近似解。蒙特卡羅方法解決問(wèn)題的過(guò)程為：構(gòu)造或描述概率過(guò)程，實(shí)現(xiàn)從已知概率分布抽樣，建立各種估計(jì)量并得到結(jié)果。該方法已廣泛應(yīng)用于自然科學(xué)、工程技術(shù)和經(jīng)濟(jì)等各領(lǐng)域。當(dāng)試驗(yàn)費(fèi)用較大、試驗(yàn)危害較大或不可能實(shí)現(xiàn)時(shí)，該方法非常實(shí)用［4］。

2　質(zhì)偏子彈質(zhì)心極限位置

某型子母彈子彈由彈體、藥柱、子彈引信等零部件組成，設(shè)計(jì)平均質(zhì)量61 g。子彈橫截面為正六方形［5］，其內(nèi)切圓直徑25.35 mm，子彈高59 mm。將坐標(biāo)系原點(diǎn)設(shè)定在子彈體底面中心，子彈外形及坐標(biāo)系OUVW選擇如圖1所示（子彈內(nèi)部結(jié)構(gòu)不對(duì)稱，U軸負(fù)向裝有電池，W軸負(fù)向裝有作動(dòng)器）。

圖1　子彈外形及坐標(biāo)系選擇

理論上講，任一零部件的任一尺寸變化，都會(huì)影響到子彈質(zhì)量和質(zhì)心位置，但質(zhì)量較小零部件（如0.003 g的探針和0.04 g的接電片）上所有尺寸以及質(zhì)量較大零部件次要尺寸（如引信體螺紋孔直徑和深度）變化對(duì)子彈整體質(zhì)量和質(zhì)心位置影響很小，可忽略。下面只研究子彈質(zhì)量較大零部件對(duì)質(zhì)心影響較大的尺寸在公差范圍內(nèi)變化對(duì)子彈質(zhì)心在其徑向兩個(gè)正交坐標(biāo)軸方向上變化的影響。

利用Solidworks軟件按圖定平均尺寸在計(jì)算機(jī)上制作三維造型，得到子彈平均尺寸下的質(zhì)心坐標(biāo)。為了得到質(zhì)心在軸向（V向）的最大值，將與之相關(guān)的高度尺寸（如彈體高度）調(diào)整到最大；零件質(zhì)心高度若大于子彈質(zhì)心高度，則將其上的孔和槽的尺寸調(diào)整到最小，使其質(zhì)量增大，整體質(zhì)心上移；若小于子彈質(zhì)心高度，則作相反處理。調(diào)整引信體、彈體、藥柱等零部件22個(gè)主要尺寸極值后，得質(zhì)心軸向坐標(biāo)最大值和最小值分別為29.42 mm和28.83 mm。同理，為了得到質(zhì)心在徑向某坐標(biāo)軸上的最大值，將在該坐標(biāo)軸正向的孔和槽的尺寸調(diào)整到最小，裝配在其上的零件的外形尺寸調(diào)整為最大，使得質(zhì)心位置向正向移動(dòng)。在該坐標(biāo)軸負(fù)向的孔、槽及零件作相反處理。

圖定平均尺寸情況下得到的質(zhì)心坐標(biāo)以及通過(guò)上述方法調(diào)整尺寸極值后得到的子彈質(zhì)心變化范圍如表1所示。同理，也可得到子彈質(zhì)量分布范圍為（59.72 g，62.35 g）。

表1　調(diào)整尺寸后所得子彈質(zhì)心坐標(biāo)極值

3　利用蒙特卡羅法模擬質(zhì)偏子彈固定方位裝配時(shí)的母彈徑向質(zhì)偏

文中涉及的母彈由15層子彈組成，每層子彈有30枚。母彈坐標(biāo)系及母彈中某一層子彈的排列方式見(jiàn)圖2所示。假設(shè)子彈裝入母彈時(shí)，由于子彈引信保險(xiǎn)方位需求或其它原因，子彈均按固定方位裝配，具體方位為子彈U軸正向與母彈X軸正向同向，子彈W軸正向與母彈Z軸正向同向。假設(shè)母彈結(jié)構(gòu)件或其組合均關(guān)于坐標(biāo)原點(diǎn)對(duì)稱，則其質(zhì)心位于坐標(biāo)原點(diǎn)。因此在計(jì)算母彈徑向質(zhì)心位置時(shí)，只需考慮質(zhì)偏子彈對(duì)母彈徑向質(zhì)心變化的影響。

圖2　子彈排列方式及母彈坐標(biāo)系

母彈質(zhì)心計(jì)算公式如下：

式中：Xm和Zm為母彈質(zhì)心徑向坐標(biāo)值；mi為子彈質(zhì)量；xi，和zi均為子彈質(zhì)心坐標(biāo)；M為母彈減去子彈的結(jié)構(gòu)件質(zhì)量之和，M=10.5 kg。參考某高炮彈頭觸發(fā)引信機(jī)械零件質(zhì)量分布特性［6］研究結(jié)果，約有半數(shù)種類的引信零件質(zhì)量特性服從正態(tài)分布或威布爾分布，整體引信質(zhì)量服從威布爾分布。文中分別假設(shè)子彈質(zhì)量服從均值為平均質(zhì)量、標(biāo)準(zhǔn)差分別為六分之一質(zhì)量變化量和十分之一質(zhì)量變化量的正態(tài)分布和尺度參數(shù)為61.85、形狀參數(shù)為435的威布爾分布（使得滿足該參數(shù)分布的子彈質(zhì)量極值與上文中得到的子彈質(zhì)量分布范圍相符），假設(shè)子彈質(zhì)心的偏移量滿足均值為坐標(biāo)范圍平均值，標(biāo)準(zhǔn)差分別為六分之一質(zhì)心變化量、十分之一質(zhì)心變化量的正態(tài)分布。

計(jì)算機(jī)模擬母彈的徑向質(zhì)心坐標(biāo)時(shí)，先按假設(shè)的分布規(guī)律隨機(jī)生成450個(gè)子彈的質(zhì)心坐標(biāo)，將其轉(zhuǎn)換到母彈坐標(biāo)系，通過(guò)上述質(zhì)心計(jì)算公式即可求得一次模擬結(jié)果。重復(fù)上述過(guò)程106次。

假設(shè)子彈質(zhì)量服從均值為平均質(zhì)量、標(biāo)準(zhǔn)差為六分之一質(zhì)量變化量的正態(tài)分布，子彈質(zhì)心的偏移量服從標(biāo)準(zhǔn)差為六分之一質(zhì)心變化量的正態(tài)分布，則子彈徑向偏心距（徑向坐標(biāo)平方和的平方根）均值約為0.235。利用蒙特卡羅方法分別模擬104次、105次和106次，所得結(jié)果如表2所示。母彈質(zhì)心在徑向坐標(biāo)軸上的分布情況如圖3所示。對(duì)結(jié)果進(jìn)行χ2擬合優(yōu)度檢驗(yàn)和Lilliefors檢驗(yàn)，表明徑向質(zhì)心坐標(biāo)及偏心距均服從正態(tài)分布。

表2　應(yīng)用蒙特卡羅方法模擬不同次數(shù)得到的母彈質(zhì)心坐標(biāo)結(jié)果對(duì)比

圖3　模擬106次后母彈質(zhì)心分布

由表2得，模擬106次后母彈質(zhì)心坐標(biāo)數(shù)值已經(jīng)收斂，故下文均按106次模擬。

假設(shè)子彈質(zhì)量和質(zhì)心的徑向偏移量服從上文所列其它分布情況，模擬結(jié)果如表3所示。由于文中涉及的子彈偏心量很小，不能確定母彈偏心量小是否由此引起，所以將子彈在兩徑向坐標(biāo)軸上的偏心量分別放大不同倍數(shù)，放大倍數(shù)及結(jié)果（以質(zhì)量服從正態(tài)分布，標(biāo)準(zhǔn)差為六分之一質(zhì)量變化量；質(zhì)心服從正態(tài)分布，標(biāo)準(zhǔn)差為六分之一質(zhì)心變化量為例）如表4所示。

表3　子彈固定方位裝配時(shí)母彈質(zhì)心坐標(biāo)模擬結(jié)果

由表3和表4可得，母彈徑向質(zhì)偏均值與子彈徑向質(zhì)偏均值成正比，比值約為母彈總質(zhì)量與子彈質(zhì)量和之比；質(zhì)量分布特性和質(zhì)心分布規(guī)律影響偏心距的取值范圍，但無(wú)論何種情況，母彈質(zhì)心兩徑向坐標(biāo)軸上的坐標(biāo)值以及其偏心距均服從正態(tài)分布。因此，在用固定方位裝配子彈的情況下，應(yīng)盡量減小子彈質(zhì)心的徑向偏移量。

表4　子彈偏心量放大后母彈質(zhì)心坐標(biāo)模擬結(jié)果

4　利用蒙特卡羅法模擬質(zhì)偏子彈隨機(jī)方位裝配時(shí)的母彈徑向質(zhì)偏

對(duì)于文中研究的子母彈來(lái)說(shuō)，子彈橫截面為正六方形，繞子彈坐標(biāo)系的V軸旋轉(zhuǎn)60°、120°、180°、240° 和300°后均能裝入母彈。模擬實(shí)際情況，假設(shè)子彈裝入前繞V軸旋轉(zhuǎn)的角度θ服從上述6種情況的均勻分布。旋轉(zhuǎn)后子彈質(zhì)心坐標(biāo)U'和W'滿足公式U'= Ucosθ+Wsinθ，W'=-Usinθ+Wcosθ。由概率論與數(shù)理統(tǒng)計(jì)知識(shí)可知，假設(shè)變量X服從正態(tài)分布，則X的任意線性函數(shù)也服從正態(tài)分布［7］，當(dāng)θ為常數(shù)時(shí)，U' 和W'服從正態(tài)分布。旋轉(zhuǎn)上述角度后子彈質(zhì)心的坐標(biāo)滿足表5。

表5　旋轉(zhuǎn)后子彈質(zhì)心坐標(biāo)分布規(guī)律

計(jì)算機(jī)模擬母彈的徑向質(zhì)心坐標(biāo)時(shí)，先通過(guò)均勻分布隨機(jī)數(shù)模擬子彈旋轉(zhuǎn)角度，并按照表5中的分布規(guī)律隨機(jī)生成450個(gè)子彈的質(zhì)心坐標(biāo)，將其轉(zhuǎn)換到母彈坐標(biāo)系中，通過(guò)上述質(zhì)心計(jì)算公式即可求得一次模擬結(jié)果。重復(fù)上述過(guò)程106次。

假設(shè)質(zhì)量服從正態(tài)分布N（61，0.438 32），徑向質(zhì)心服從標(biāo)準(zhǔn)差為六分之一質(zhì)心變化量的正態(tài)分布，模擬結(jié)果見(jiàn)圖4所示。假設(shè)質(zhì)量服從威布爾分布W （61.85，435），徑向質(zhì)心服從標(biāo)準(zhǔn)差為六分之一質(zhì)心變化量的正態(tài)分布，模擬結(jié)果見(jiàn)圖5所示。上述情況下的母彈質(zhì)心坐標(biāo)求解結(jié)果如表6所示。

由于文中涉及的子彈偏心量很小，不能確定母彈偏心量小是否由此引起，所以將子彈在兩徑向坐標(biāo)軸上的偏心量分別放大不同倍數(shù)，放大倍數(shù)及結(jié)果（以質(zhì)量服從正態(tài)分布，標(biāo)準(zhǔn)差為六分之一質(zhì)量變化量為例）如表7所示。

圖4　質(zhì)量正態(tài)分布時(shí)母彈質(zhì)心徑向偏心距分布

圖5　質(zhì)量威布爾分布時(shí)母彈質(zhì)心徑向偏心距分布

表6　子彈隨機(jī)方位裝配時(shí)母彈質(zhì)心坐標(biāo)模擬結(jié)果

表7　子彈偏心量放大后母彈質(zhì)心坐標(biāo)模擬結(jié)果

綜合表6和表7可以看出，母彈質(zhì)心徑向偏心距均值及最大值均遠(yuǎn)小于子彈徑向偏心距，其最小值趨向于零，且偏心距近似服從威布爾分布。因此，在用隨機(jī)方位（僅考慮子彈繞自身對(duì)稱軸的旋轉(zhuǎn)方位）裝配子彈的情況下，只需保證子彈質(zhì)偏不過(guò)大即可保證母彈質(zhì)偏較小。

5　結(jié)論

文中將蒙特卡羅方法用于分析質(zhì)偏子彈裝配方位對(duì)母彈徑向質(zhì)偏的影響。仿真結(jié)果表明，子彈固定方位裝配時(shí)，母彈質(zhì)心兩徑向坐標(biāo)軸上的坐標(biāo)均值與子彈相應(yīng)坐標(biāo)軸上的坐標(biāo)均值成正比。因此，必須在子彈設(shè)計(jì)時(shí)盡量減小子彈的徑向質(zhì)偏量。子彈隨機(jī)方位（僅考慮子彈繞自身對(duì)稱軸的旋轉(zhuǎn)方位）裝配時(shí)，母彈質(zhì)心徑向偏心距均值及最大值均遠(yuǎn)小于子彈徑向偏心距，只需保證子彈質(zhì)偏不過(guò)大即可。

［1］ 楊啟仁.子母彈飛行動(dòng)力學(xué)［M］.北京：國(guó)防工業(yè)出版社，1999：141-142.

［2］葛任偉，歐陽(yáng)勇，李佳.基于Monte Carlo方法的戰(zhàn)斗部部件質(zhì)心偏移量仿真分析［J］.中國(guó)機(jī)械工程，2012，23（7）：782-785.

［3］聞泉，王雨時(shí).子母彈徑向質(zhì)心位置的蒙特卡羅模擬［J］.彈箭與制導(dǎo)學(xué)報(bào)，2006，26（2）：403-405.

［4］朱平.擬蒙特卡羅方法的若干研究與應(yīng)用［D］.浙江：浙江大學(xué)，2010.

［5］馮彥哲，王雨時(shí).使母彈威力最大化的子母彈子彈橫截面形狀與布局分析［J］.彈箭與制導(dǎo)學(xué)報(bào)，2008，28 （3）：135-142.

［6］王曉鵬，王雨時(shí)，聞泉，等.某高炮彈頭觸發(fā)引信機(jī)械零件質(zhì)量分布特性［J］.四川兵工學(xué)報(bào)，2015，36（5）：63-68.

［7］盛驟，謝式千，潘承毅.概率論與數(shù)理統(tǒng)計(jì)［M］.北京：高等教育出版社，2001：65-66.

Effects of Assembling Orientation of Submunition with Centroid Excursion on Radial Centroid Excursion of Cargo Projectile

CAI Hongjie，WEN Quan，WANG Yushi，ZHANG Zhibiao
（School of Mechanical Engineering，Nanjing University of Science and Technology，Nanjing 210094，China）

Dimensional variation and assembling orientation of submunition with centroid excursion affect radial centroid excursion of cargo projectile.To analyze it，the effects of dimensional variation of submunition on ultimate centroid excursion of submunition was analyzed by Solidworks firstly.Then，the effects of assembling orientation of submunition with centroid excursion on radial centroid excursion of cargo projectile was analyzed using Monte-Carlo method.Submunitions obeyed different mass distribution and centroid distribution were assembled in the cargo projectile in different orientations and the radial centroid distribution of cargo projectile was calculated.To effectively decrease the radial centroid excursion of cargo projectile，the centroid excursion of submunition must be reduced when submunition was assembled in certain orientation，while it needn’t to be specially controlled when assembled in random orientation.

cargo projectile；Monte-Carlo method；design principle；centroid excursion；assembling orientation

10.15892/j.cnki.djzdxb.2016.01.022

TJ413.3

A

2015-03-31

蔡泓杰（1991-），男，江蘇無(wú)錫人，碩士研究生，研究方向：引信及彈藥技術(shù)。