王麗群， 楊國來， 劉俊民， 葛建立

(1.南京理工大學(xué) 機械工程學(xué)院， 江蘇 南京 210094；2.內(nèi)蒙古北方重工業(yè)集團有限公司， 內(nèi)蒙古 包頭 014033)

?

面向火炮射擊密集度的隨機因素穩(wěn)健設(shè)計

王麗群1， 楊國來1， 劉俊民2， 葛建立1

(1.南京理工大學(xué) 機械工程學(xué)院， 江蘇 南京 210094；2.內(nèi)蒙古北方重工業(yè)集團有限公司， 內(nèi)蒙古 包頭 014033)

為了實現(xiàn)對射擊密集度的有效控制，解決考慮隨機因素的射擊密集度優(yōu)化設(shè)計問題，利用基于隨機模型和隨機優(yōu)化的穩(wěn)健設(shè)計理論，提出一種指標要求導(dǎo)向型的隨機因素參數(shù)區(qū)間計算方法?；诿商乜迥M，采用六自由度外彈道模型，利用某型大口徑榴彈炮數(shù)據(jù)求解射擊密集度并研究其統(tǒng)計特性。利用不合格概率、靈敏度指數(shù)兩個穩(wěn)健目標準則，構(gòu)建射擊密集度穩(wěn)健設(shè)計隨機模型。采用帶精英策略非支配排序遺傳算法并結(jié)合反向傳播神經(jīng)網(wǎng)絡(luò)代理模型技術(shù)對穩(wěn)健設(shè)計隨機模型進行優(yōu)化求解，確定符合指標要求的隨機因素參數(shù)區(qū)間。對比分析結(jié)果表明，該方法可同時保證射擊密集度指標的最優(yōu)性和穩(wěn)健性，對射擊密集度的預(yù)測具有一定可行性。

兵器科學(xué)與技術(shù)； 射擊密集度； 穩(wěn)健設(shè)計； 隨機因素； 多目標遺傳算法； 神經(jīng)網(wǎng)絡(luò)

0 引言

射擊密集度是火炮武器系統(tǒng)的重要戰(zhàn)術(shù)技術(shù)指標，也是射擊精度的重要組成部分。在火炮武器系統(tǒng)研制過程中，經(jīng)常發(fā)生射擊密集度達不到指標要求的情況，這主要是因為不能有效地控制影響射擊密集度的設(shè)計參量。而射擊密集度是一個涉及彈、炮、藥、氣象環(huán)境等眾多因素的系統(tǒng)問題，彈丸起始擾動、彈藥誤差、氣象條件等都將影響射擊密集度，并且這些因素是典型的隨機因素，均在一定的誤差范圍內(nèi)隨機變化。因此，找到最佳的隨機參數(shù)組合，確定各隨機因素的合理參數(shù)區(qū)間，可以為火炮關(guān)鍵結(jié)構(gòu)參數(shù)優(yōu)化設(shè)計以及彈藥的加工提供理論參考，進而在火炮的總體設(shè)計階段掌握影響射擊密集度的設(shè)計因素，實現(xiàn)對射擊密集度的有效控制。

目前對于火炮射擊密集度的研究大多集中在以炮口擾動為目標[1-2]的發(fā)射動力學(xué)分析，而研究火炮發(fā)射全過程中隨機因素較少。對于隨機因素的影響分析，曹寧等[3]研究了彈丸初速、彈丸質(zhì)量和彈丸偏心距等因素對車載炮射擊密集度的影響，建立了影響程度逐步回歸分析算法。吳宏等[4]采用均勻設(shè)計與逐步回歸結(jié)合的方法分析了彈炮間隙、質(zhì)量偏心、動不平衡角、發(fā)射藥量、擠進壓力、彈丸質(zhì)量和風(fēng)速對射擊密集度的影響。現(xiàn)有對射擊密集度的處理，均是采用誤差合成或多次計算取均值的方法，將射擊密集度這一隨機模型轉(zhuǎn)化為確定性模型求解，隨機變量或射擊密集度的隨機性沒有得到很好體現(xiàn)。同時，現(xiàn)有射擊密集度隨機因素的研究均不涉及參數(shù)優(yōu)化設(shè)計，無法根據(jù)指標要求直接計算給出合適的、富有理論依據(jù)的參數(shù)區(qū)間。針對以上問題，本文從“彈丸出炮口運動散布—射擊密集度”過程出發(fā)。采用基于隨機模型和隨機優(yōu)化的穩(wěn)健設(shè)計理論，先構(gòu)建射擊密集度隨機模型，再采用帶精英策略非支配排序遺傳算法(NSGA-Ⅱ)并結(jié)合反向傳播(BP)神經(jīng)網(wǎng)絡(luò)對隨機模型進行優(yōu)化求解，確定了符合指標要求的隨機因素參數(shù)區(qū)間，探討了考慮隨機因素的射擊密集度優(yōu)化設(shè)計問題。

1 穩(wěn)健設(shè)計基本思想

穩(wěn)健設(shè)計的基本思想是在已知輸出空間下，找出最合適的輸入空間，使產(chǎn)品的質(zhì)量特性既限定在給定的輸出空間內(nèi)，它的均值又與目標值的差異為最小[5-6]。對于含有隨機因素的隨機問題，輸入空間指隨機設(shè)計變量的均值μx及其容差Δx-、Δx+. 若隨機設(shè)計變量服從正態(tài)分布，則輸入空間為其均值μx和標準差σ. 圖1展示了穩(wěn)健優(yōu)化設(shè)計解與一般優(yōu)化設(shè)計解的區(qū)別，其中g(shù)i(i=1，2，…，n)為隨機性約束函數(shù)，ai(i=1，2，…，n)為預(yù)先設(shè)定應(yīng)滿足的概率值。

圖1 穩(wěn)健設(shè)計最優(yōu)解Fig.1 Optimal solution for robust design of random model

2 射擊密集度穩(wěn)健設(shè)計隨機模型構(gòu)建

2.1 隨機因素的分類、選取與隨機性

工程實際中，隨機因素可以根據(jù)其不同特點分為可控因素和不可控因素兩類?？煽匾蛩匾卜Q隨機設(shè)計變量，是指在實際中可以控制均值和容差的改變減小輸出特性的波動和偏差的因素。數(shù)學(xué)表示為

xT=(x1,x2,…，xn)∈(Ω,ζ,P)?Rn.

(1)

不可控因素也稱隨機參數(shù)，是指對輸出特性有影響，而在實際操作中難以人為控制的因素，例如使用條件、操作人員、環(huán)境因素等。數(shù)學(xué)表示為

zT=(z1,z2,…，zk)∈(Ω,ζ,P)?Rk.

(2)

(1)式、(2)式中：x為n維隨機設(shè)計變量；z為k維隨機參數(shù)；(Ω,ζ,P)為概率空間，Ω為隨機事件，ζ為事件的全體，P為事件發(fā)生的概率，表示x和z都是屬于概率空間的向量。

對射擊密集度有影響的隨機因素描述為某種概率分布近似描述其隨機特性。在缺乏必要的統(tǒng)計試驗數(shù)據(jù)時，隨機設(shè)計變量假定為正態(tài)分布，并以隨機因素的均值μx和標準差σx作為設(shè)計變量，根據(jù)正態(tài)分布“3σ”原則抽取隨機數(shù)[7]。在本文的研究中，主要選用正態(tài)分布和隨機分布兩種分布類型，各隨機因素的具體概率分布類型見表1.

表1 隨機因素概率分布類型

2.2 射擊密集度的計算與隨機性分析

為了獲得射擊密集度的隨機分布類型，以某型大口徑榴彈炮數(shù)據(jù)編寫射擊密集度計算程序，并采用一組數(shù)據(jù)重復(fù)計算1 500次，對縱向密集度EX與橫向密集度EZ進行統(tǒng)計分析。

表2、圖2和圖3的統(tǒng)計結(jié)果是基于蒙特卡洛法模擬的計算方法，采用六自由度外彈道模型[8]，利用計算機隨機抽樣的方法求得的射擊密集度數(shù)據(jù)統(tǒng)計平均值。

表2 數(shù)據(jù)分布統(tǒng)計表

圖2 橫向密集度EZ統(tǒng)計分布直方圖Fig.2 Statistical distribution histogram of EZ

表2中的偏度與峰度分別指樣本數(shù)據(jù)標準化變量的3階、4階中心矩。偏度反映分布的對稱性，接近0則認為分布是對稱的；峰度亦是衡量偏離正態(tài) 分布的尺度，正態(tài)分布的峰度為3. 圖2、圖3中的紅線表示以該組射擊密集度數(shù)據(jù)的均值與標準差構(gòu)建的正態(tài)分布曲線。

圖3 縱向密集度lnEX統(tǒng)計分布直方圖Fig.3 Statistical distribution histogram of lnEX

由表2和圖2可以看出，橫向密集度服從正態(tài)分布；由表2和圖3可見，縱向密集度樣本數(shù)據(jù)的對數(shù)近似服從正態(tài)分布，即服從對數(shù)正態(tài)分布。在確定射擊密集度的分布類型后，可以通過多次重復(fù)計算求解射擊密集度的分布參數(shù)(均值、方差)，確定其概率密度函數(shù)并在穩(wěn)健設(shè)計中使用。

2.3 射擊密集度穩(wěn)健設(shè)計目標準則

(3)

式中：q為指標個數(shù)。P{·}∈(0,1)的數(shù)值越小，則不合格率越小，說明輸出特性越穩(wěn)定。

目標準則之二是靈敏度指數(shù)SI. 由于在工程實際中，希望每一項技術(shù)特性的實際值與目標值的差異盡量小，即希望輸出特性更加優(yōu)質(zhì)[7]，因此采用靈敏度指數(shù)SI衡量輸出特性的優(yōu)質(zhì)性，使

(4)

式中：μyj為yj的統(tǒng)計均值，當μyj=yoj時，SI=0.

概率P的計算方法主要有隨機模擬法、數(shù)值積分法和1次2階矩法等，本文采用數(shù)值積分法進行求解。依托某型大口徑榴彈炮數(shù)據(jù)，設(shè)定縱向密集 度目標值EX=1/yX=1/350，同時設(shè)定容差ΔyX為20；設(shè)定橫向密集度目標值EZ=yZ=0.6 mil，由于橫向密集度一般符合指標要求，因此放大容差設(shè)定，ΔyZ為0.3 mil. 在每次穩(wěn)健設(shè)計中，為了兼顧計算精度與求解速度，對一組隨機設(shè)計變量重復(fù)計算100次，計算獲得EX和EZ的均值μEX、μEZ與標準差σEX、σEZ，根據(jù)2.2節(jié)求得縱向密集度與橫向密集度的分布類型，代入(5)式中，獲得EX和EZ的概率密度函數(shù)。

(5)

由于正態(tài)分布與對數(shù)正態(tài)分布均不可積，因此通過(6)式高斯—洛巴托積分方法，求解EX和EZ概率密集函數(shù)的4次等效積分，計算不合格概率P.

(6)

式中：xk為積分節(jié)點坐標；Ak為積分系數(shù)；n為積分所用項數(shù)。

由于積分區(qū)間不是[-1,1]，因此按照(7)式先行轉(zhuǎn)換積分區(qū)間。

(7)

通過(4)式的定義計算靈敏度指數(shù)SI.

3 基于多目標遺傳算法的隨機模型優(yōu)化

采用模擬搜索算法對隨機模型進行優(yōu)化是隨機變量優(yōu)化設(shè)計的一種簡便有效的方法。本文采用NSGA-II多目標遺傳算法對射擊密集度穩(wěn)健設(shè)計隨機模型進行優(yōu)化，即通過對設(shè)計變量的隨機模擬搜索，尋找穩(wěn)健設(shè)計的最優(yōu)解，以確定最佳參數(shù)組合和各隨機因素參數(shù)區(qū)間。NSGA-II多目標遺傳方法基于Pareto最優(yōu)概念，是一種啟發(fā)式模擬搜索算法，其采用非支配排序過程、精英保留策略和無參數(shù)小生境操作算子[9-10]，可以獲得設(shè)計變量的最優(yōu)解集。但由第2節(jié)可知，每次穩(wěn)健設(shè)計構(gòu)建隨機模型的過程需要調(diào)用射擊密集度計算程序進行一定數(shù)量的隨機模擬，然后進行約束準則與目標準則的求解。為了保證計算精度，重復(fù)計算次數(shù)很大，產(chǎn)生巨大計算費用而不便用于后序優(yōu)化。為了提高計算效率，采用基于代理模型的優(yōu)化策略。本文的研究中選用BP神經(jīng)網(wǎng)絡(luò)，它具有良好的學(xué)習(xí)能力和非線性映射能力，一個3層的BP神經(jīng)網(wǎng)絡(luò)可以任意逼近復(fù)雜的非線性函數(shù)。

3.1 優(yōu)化設(shè)計變量與優(yōu)化數(shù)學(xué)模型

表3 隨機設(shè)計變量及取值范圍

基于射擊密集度穩(wěn)健設(shè)計隨機模型，建立以彈丸起始擾動參數(shù)和彈丸結(jié)構(gòu)參數(shù)誤差為隨機設(shè)計變量，以射擊密集度指標要求為約束條件，穩(wěn)健目標準則為優(yōu)化目標的數(shù)學(xué)模型為

s.t.

(8)

3.2 BP神經(jīng)網(wǎng)絡(luò)代理模型構(gòu)建與可信度評估

參見圖4，本文構(gòu)造BP神經(jīng)網(wǎng)絡(luò)代理模型經(jīng)過3個步驟：首先依托某型大口徑榴彈炮數(shù)據(jù)，采用最優(yōu)拉丁超立方試驗設(shè)計安排1 000組試驗樣本；然后代入射擊密集度穩(wěn)健設(shè)計隨機模型中，得到不合格率P、靈敏度指數(shù)SI、縱向密集度均值EX、橫向密集度均值EZ作為輸出，獲得輸入到輸出的訓(xùn)練樣本；利用Matlab神經(jīng)網(wǎng)絡(luò)工具箱訓(xùn)練一個3層的BP神經(jīng)網(wǎng)絡(luò)。

圖4 BP神經(jīng)網(wǎng)絡(luò)代理模型構(gòu)建過程Fig.4 Model building process of BP neural network

為了提高訓(xùn)練精度，在訓(xùn)練前進行樣本數(shù)據(jù)的歸一化處理，消除各設(shè)計變量量綱與量級之間的差別。輸入層節(jié)點數(shù)設(shè)為15，隱含層節(jié)點數(shù)經(jīng)反復(fù)試算設(shè)定為30，其傳遞函數(shù)采用tansig函數(shù)；輸出層節(jié)點數(shù)為4，其傳遞函數(shù)采用purelin函數(shù)。采用精度較高的trainlm訓(xùn)練函數(shù)和learngdm學(xué)習(xí)函數(shù)對神經(jīng)網(wǎng)絡(luò)進行訓(xùn)練。

采用最優(yōu)拉丁方試驗設(shè)計生成100組驗證樣本，帶入原射擊密集度穩(wěn)健設(shè)計隨機模型中獲得相應(yīng)輸出，計算確定性系數(shù)R2檢驗神經(jīng)網(wǎng)絡(luò)模型可信度。

(9)

表4 神經(jīng)網(wǎng)絡(luò)可信度檢驗

3.3 優(yōu)化結(jié)果與分析

優(yōu)化計算得到隨機設(shè)計變量的Pareto最優(yōu)解集，優(yōu)化結(jié)果的Pareto前沿如圖5所示。

圖5 優(yōu)化結(jié)果Pareto前沿Fig.5 Pareto front

由圖5可見，Pareto前沿中的設(shè)計點，穩(wěn)健目標準則均得到不同程度的優(yōu)化。Pareto前沿左側(cè)的設(shè)計P較小，輸出特性較為穩(wěn)定；右側(cè)的設(shè)計SI較小，輸出特性更加優(yōu)質(zhì)；中部的設(shè)計較為折中。在Pareto前沿的左、中、右側(cè)分別選擇一組設(shè)計方案，即圖5中A、B、C設(shè)計點，以使選擇的方案具有代表性。帶入原射擊密集度穩(wěn)健設(shè)計模型中進行詳細的分析計算，并與采用遺傳算法優(yōu)化(不含穩(wěn)健設(shè)計)得到的普通優(yōu)化結(jié)果進行對比，計算結(jié)果見表5.

表5 不同設(shè)計的對比

由表5可知，基于隨機模型與隨機優(yōu)化的穩(wěn)健設(shè)計得到的3組優(yōu)化結(jié)果雖不相同，但都可以滿足預(yù)先設(shè)定的射擊密集度EX=1/350、EZ=0.6 mil的指標要求?；谶z傳算法未經(jīng)穩(wěn)健設(shè)計的普通優(yōu)化結(jié)果，在均值上可以滿足EX=1/350的目標，但偏離目標程度較大，這意味著無法實現(xiàn)精確的設(shè)計；過高的指標意味著過高的技術(shù)成本，這也不是我們所希望的；同時穩(wěn)健性較差，計算結(jié)果必將不穩(wěn)定。需要特別指出的是，表5中目標準則SI變化規(guī)律滿足圖5，目標準則P與圖5有一定差異，這一方面是代理模型精度造成的，另一方面由于3組設(shè)計均滿足約束條件，Pareto前沿中穩(wěn)健目標準則P結(jié)果相差并不明顯，約為0.005左右，這種差異是可接受的。

兼顧不合格概率P和靈敏度指數(shù)SI兩個穩(wěn)健目標準則，選擇C設(shè)計點作為最終的設(shè)計，從而確定縱向密集度EX=1/350、橫向密集度EZ=0.6 mil時各隨機因素的參數(shù)區(qū)間，數(shù)據(jù)圓整后如表6所示。

表6 隨機因素參數(shù)區(qū)間優(yōu)化結(jié)果

相似的，可以根據(jù)工程射擊密集度指標需要，設(shè)定不同的目標，例如縱向密集度EX為1/300、1/400，采用相同的方法就可以獲得不同指標下的隨機因素參數(shù)區(qū)間。

4 結(jié)論

本文從“彈丸出炮口運動散布- 射擊密集度”過程出發(fā)，采用基于隨機模型和隨機優(yōu)化的穩(wěn)健設(shè)計理論，探討了考慮隨機因素的射擊密集度優(yōu)化設(shè)計問題，以某型大口徑榴彈炮為例初步驗證了該方法的有效性。

目前僅研究了火炮發(fā)射過程的一個階段，若將該方法進一步擴展到火炮關(guān)鍵結(jié)構(gòu)參數(shù)誤差，從理論上揭示其與射擊密集度的傳遞機理及影響關(guān)系，將在總體設(shè)計階段實現(xiàn)對射擊密集度的有效控制。另一方面，火炮發(fā)射過程中的隨機因素目前尚缺乏相關(guān)統(tǒng)計試驗數(shù)據(jù)，未來若從試驗角度揭示隨機因素的分布類型，引入更加完整的氣象模型與數(shù)據(jù)，計算結(jié)果將更加準確。

References)

[1] 梁傳建, 楊國來, 王曉鋒. 基于神經(jīng)網(wǎng)絡(luò)和遺傳算法的火炮結(jié)構(gòu)動力學(xué)優(yōu)化[J]. 兵工學(xué)報, 2015, 36(5): 789-795. LIANG Chuan-jian, YANG Guo-lai, WANG Xiao-feng. Structural dynamics optimization of gun based on neural networks and genetic algorithms[J]. Acta Armamentarii, 2015, 36(5): 789-795. (in Chinese)

[2] 謝潤, 楊國來. 自行高炮行進間射擊炮口響應(yīng)特性研究[J]. 兵工學(xué)報, 2014, 35(8): 1158-1163. XIE Run, YANG Guo-lai. Research on response characteristics of muzzle for self-propelled antiaircraft gun firing on the move[J]. Acta Armamentarii, 2014, 35(8): 1158-1163. (in Chinese)

[3] 曹寧, 王曉鋒, 徐亞棟, 等. 基于逐步回歸法的車載炮射擊精度影響因素分析[J]. 南京理工大學(xué)學(xué)報, 2013, 37(4): 557-560. CAO Ning, WANG Xiao-feng, XU Ya-Dong, et al. Evaluation of influencing factors on firing accuracy of wheeled artillery based on stepwise regression analysis[J]. Journal of Nanjing University of Science and Technology, 2013, 37(4): 557-560. (in Chinese)

[4] 吳宏, 趙有守. 隨機因素對彈丸射擊精度影響的分析方法[J]. 彈箭與制導(dǎo)學(xué)報, 2001, 21(3): 56-59. WU Hong, ZHAO You-shou. A method to study fire precision of projectile effected by random factors[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2001, 21 (3): 56-59. (in Chinese)

[5] 申立娟, 楊軍, 趙宇. 基于方差傳遞模型的飛機蒙皮拉形工藝穩(wěn)健設(shè)計[J]. 機械工程學(xué)報, 2011, 47(1): 145-151. SHEN Li-juan, YANG Jun, ZHAO Yu. Robust design of aircraft skin stretch forming process based on transmitted variance model[J]. Journal of Mechanical Engineering, 2011, 47(1): 145-151. (in Chinese)

[6] Huang B Q, Du X P. Analytical robustness assessment for robust design[J]. Structural and Multidisciplinary Optimization, 2007, 34(2): 123-137.

[7] 陳立周. 穩(wěn)健設(shè)計[M]. 北京: 機械工業(yè)出版社, 2000：255-263. CHEN Li-zhou. Robust design[M]. Beijing: China Machine Press, 2000: 255-263. (in Chinese)

[8] 郭錫福. 遠程火炮武器系統(tǒng)射擊精度分析[M]. 北京: 國防工業(yè)出版社, 2004: 38-60. GUO Xi-fu. Firing accuracy analysis for long range gun weapon systems[M]. Beijing: National Defense Industry Press, 2004: 38-60. (in Chinese)

[9] 蔡云驤, 周志勇, 徐小倩, 等. 一種基于精英策略非支配排序遺傳算法-Ⅱ的多層吸波涂層結(jié)構(gòu)設(shè)計方法[J].兵工學(xué)報, 2015, 36(8): 1574-1579. CAI Yun-xiang, ZHOU Zhi-yong, XU Xiao-qian, et al. A design method for multilayered microwave absorbing coating structure based on NSGA-II[J]. Acta Armamentarii, 2015, 36(8): 1574-1579. (in Chinese)

[10] 錢輝仲, 顧克秋, 彭迪, 等. 基于NSGA-II算法的超輕型火炮搖架多目標優(yōu)化設(shè)計[J]. 機械設(shè)計, 2012, 29(6): 36-39. QIAN Hui-zhong, GU Ke-qiu, PENG Di, et al. Multi- objective optimization design of ultralightweight artillery cradle based on NSGA-II algorithm[J]. Journal of Machine Design, 2012, 29(6): 36-39. (in Chinese)

Robust Design of Random Factors on Gun Firing Dispersion

WANG Li-qun1, YANG Guo-lai1, LIU Jun-min2, GE Jian-li1

(1.School of Mechanical Engineering，Nanjing University of Science and Technology，Nanjing 210094，Jiangsu，China； 2.Inner Mongolia North Heavy Industries Group Co., Ltd.，Baotou 014033，Inner Mongolia，China)

To study the effective ways to improve the firing dispersion and solve its optimal design problem considering random factors, a calculation method of parameters interval oriented by the predefined indicators is proposed by using robust design theory based on random model and random optimization. A large caliber howitzer data is used to calculate firing dispersion based on Monte Carlo simulation and six degrees of freedom exterior ballistics model, and the statistical properties of firing dispersion are analyzed. A random model of robust design is built by calculating the objective criteria of robust design including unqualified probability and sensitivity index. Based on multi-objective genetic algorithm (NSGA-II) and BP neural network, the random model is optimized to determine the parameters interval of random factors which meet the requirements. The analysis results show that the proposed method can ensure the optimality and robustness of firing dispersion, and is feasible for the prediction of firing dispersion.

ordnance science and technology; firing dispersion; robust design; random factor; multi-objective genetic algorithm; neural network

2016-05-13

國家“973”計劃項目(6132490303)；國家自然科學(xué)基金項目(11572158)；中央高?；究蒲袠I(yè)務(wù)費專項資金項目(30915118825)

王麗群(1992—)，男，博士研究生。E-mail：lqwangnjust101@163.com； 楊國來(1968—)，男，教授，博士生導(dǎo)師。E-mail：yyanggl@mail.njust.edu.cn

TJ301

A

1000-1093(2016)11-1983-06

10.3969/j.issn.1000-1093.2016.11.003