迫擊炮炮箍位置優(yōu)化設(shè)計

林智文,葛建立,張鴻浩,王雪嫣,楊國來

(1.南京理工大學(xué) 機(jī)械工程學(xué)院,南京 210094; 2.北京特種機(jī)電研究所,北京 100012)

【裝備理論與裝備技術(shù)】

迫擊炮炮箍位置優(yōu)化設(shè)計

林智文1,葛建立1,張鴻浩2,王雪嫣1,楊國來1

(1.南京理工大學(xué) 機(jī)械工程學(xué)院,南京 210094; 2.北京特種機(jī)電研究所,北京 100012)

針對迫擊炮炮口擾動問題，通過建立某120 mm迫擊炮結(jié)構(gòu)動力學(xué)模型，采用有限元法、層次分析法、Kriging模型和遺傳算法相結(jié)合的方法，進(jìn)行了迫擊炮結(jié)構(gòu)動力學(xué)優(yōu)化。以試驗數(shù)據(jù)為訓(xùn)練樣本，建立了基于迫擊炮炮箍位置與炮口擾動之間的非線性映射關(guān)系。通過構(gòu)建優(yōu)化目標(biāo)函數(shù)，利用遺傳算法對目標(biāo)函數(shù)進(jìn)行求解，實現(xiàn)了迫擊炮炮口擾動量最小值的優(yōu)化。研究表明，在多工況下優(yōu)化后的迫擊炮炮箍位置，有效減少了炮口擾動，對于研究迫擊炮炮口擾動問題具有一定的參考意義。

迫擊炮；層次分析法；Kriging模型；結(jié)構(gòu)動力學(xué)優(yōu)化

迫擊炮在山地戰(zhàn)、城市戰(zhàn)中具有不可替代的作用。目前對迫擊炮的研究非?；钴S。迫擊炮的射擊精度是其研制過程中需要考慮的重要問題，而炮口擾動對射擊精度的影響很大。研究表明，座鈑、炮箍、緩沖機(jī)和炮架等零部件的結(jié)構(gòu)設(shè)計和布局，以及迫擊炮炮彈的質(zhì)心、彈炮間隙等都是影響迫擊炮炮口擾動的關(guān)鍵因素。其中，炮箍位置的選定更是迫擊炮總體布置的重要問題[1]。

隨著結(jié)構(gòu)優(yōu)化方法的日趨成熟，針對迫擊炮采用某種優(yōu)化方法，可以降低迫擊炮炮口擾動。鄔欣明[2]應(yīng)用虛擬樣機(jī)和聯(lián)合仿真技術(shù)研究了某車載自動迫擊炮的射擊穩(wěn)定性和炮口擾動問題。陳杰等[3]結(jié)合正交試驗的方法進(jìn)行了火炮上架多工況結(jié)構(gòu)優(yōu)化，使得整體結(jié)構(gòu)更加合理，各工況下的上架剛強(qiáng)度均有所改善。蘭鳳崇等[4]利用層次分析法確定權(quán)重比，針對方程式賽車的車身進(jìn)行了多工況結(jié)構(gòu)優(yōu)化，提高了各工況下的結(jié)構(gòu)性能。李小剛等[5]基于雙層更新Kriging模型和鄰域培植遺傳算法，對某型號大型汽輪發(fā)電機(jī)定子端部繞組錐環(huán)固定結(jié)構(gòu)進(jìn)行了動態(tài)特性穩(wěn)健優(yōu)化設(shè)計。當(dāng)前針對多工況結(jié)構(gòu)優(yōu)化問題的研究方法有很多，但是在迫擊炮領(lǐng)域均沒有對結(jié)構(gòu)多工況和各工況權(quán)重系數(shù)進(jìn)行深入研究。

為了優(yōu)化迫擊炮多工況條件下的炮口擾動，本研究采用層次分析法確定各工況權(quán)重系數(shù)?；贙riging模型和遺傳算法相結(jié)合的方法，對迫擊炮進(jìn)行結(jié)構(gòu)動力學(xué)優(yōu)化，利用有限元分析軟件ABAQUS對優(yōu)化前后的迫擊炮總體結(jié)構(gòu)進(jìn)行動力學(xué)數(shù)值計算，通過對比分析說明該研究方法的可行性。

1 迫擊炮結(jié)構(gòu)動力學(xué)建模

1.1 有限元模型

本文研究的對象是某120 mm迫擊炮，為了更好地模擬實際情況，本文采用帶彈的全炮模型進(jìn)行計算，土壤選擇試驗常用的中硬土。

針對迫擊炮的有限元建模，對主要零部件結(jié)構(gòu)盡量保持不變，為了提高建模和計算效率，對凸臺、圓角和圓孔等結(jié)構(gòu)進(jìn)行了簡化處理。在網(wǎng)格劃分中，迫擊炮炮身、炮架、炮彈結(jié)構(gòu)和土壤網(wǎng)格模型均采用了實體單元，而座鈑結(jié)構(gòu)、局部尾翼片等結(jié)構(gòu)則采用了殼單元。迫擊炮各零部件之間的復(fù)雜連接關(guān)系也需要妥善處理，通過在炮尾尾球中心建立一個節(jié)點與炮尾底面剛性連接，座鈑駐臼中心建立一個節(jié)點與駐臼四周剛性連接，兩節(jié)點間定義一個球鉸模擬炮尾與座鈑的鉸連接關(guān)系。通過建立面接觸，模擬炮架、座鈑與土壤之間的接觸關(guān)系。

所建立的迫擊炮全炮有限元模型如圖1所示。建立坐標(biāo)系方向如下：與炮身身管軸線重合的坐標(biāo)軸定為x軸，炮尾指向炮口方向為正方向，與身管軸線垂直且指向上方的方向為y軸正方向，z軸由右手定則確定。

圖1 迫擊炮全炮有限元模型

1.2 載荷與邊界條件

迫擊炮彈在膛內(nèi)運動的時間較短，運動和受力情況復(fù)雜，整個發(fā)射期間全炮主要受到火藥氣體壓力和重力的作用。在本文的研究中，火藥氣體壓力對迫擊炮的作用是在膛底隨時間變化的等效壓力，如圖2所示。而作用在炮彈上的火藥氣體壓力則是彈底隨時間變化的等效壓力，如圖3所示。重力載荷作為常力直接加載在模型中。在土壤上施加全約束邊界條件。

圖2 膛底等效壓力

圖3 彈底等效壓力

2 多工況炮口擾動優(yōu)化設(shè)計

2.1 優(yōu)化設(shè)計變量

為了研究迫擊炮炮口擾動問題，根據(jù)該型迫擊炮的試驗及實際運用情況，選取上下極限射角45°和80°，以及常用射角60°這3種不同射角工況進(jìn)行研究。在本文研究中，選擇迫擊炮炮箍位置作為優(yōu)化設(shè)計變量，其值為炮箍距炮口的直線距離，示意圖如圖4所示。根據(jù)相關(guān)設(shè)計手冊以及考慮實際操作使用情況，確定設(shè)計變量的初值和取值范圍如表1所示。

表1 設(shè)計變量初值和取值范圍 mm

圖4 炮箍位置示意圖

2.2 各工況權(quán)重系數(shù)確定

在研究迫擊炮多工況結(jié)構(gòu)優(yōu)化問題時，每一種工況都需要對應(yīng)一個權(quán)重系數(shù)。傳統(tǒng)上多工況權(quán)重系數(shù)的確定，往往是設(shè)計人員依據(jù)認(rèn)定的工況重要程度給出一個經(jīng)驗值，但采用這種方法給出的權(quán)重系數(shù)一般只是設(shè)計人員的定性判斷，無法精確。遇到工況數(shù)量較多的情況，設(shè)計人員無法依據(jù)各工況的重要度給出相適應(yīng)的權(quán)重系數(shù)。為了解決上述問題，本文利用層次分析法，研究迫擊炮結(jié)構(gòu)優(yōu)化各工況下的權(quán)重系數(shù)。

首先，依據(jù)設(shè)計人員所要研究的問題建立決策層次，然后兩兩比較不同工況間的重要程度，構(gòu)建一個決策矩陣A= (aij)n×n，n代表決策的數(shù)目(本文中表示工況數(shù)目)，aij是第i層決策數(shù)對比第j層決策數(shù)的重要程度。其中，決策矩陣如式(1)所示，aij的重要性參考定義如表2所示。

(1)

表2 決策矩陣重要性參考定義表

在構(gòu)造判斷矩陣A之后，求出相應(yīng)的特征向量并歸一化，即可得最終優(yōu)化目標(biāo)函數(shù)各工況的權(quán)重系數(shù)。在實際應(yīng)用中，由于客觀事物的復(fù)雜性，進(jìn)行兩兩對比判斷時常帶有主觀性和片面性。因此，在構(gòu)造判斷矩陣A之后，還需要進(jìn)行一致性檢驗，若不滿足一致性要求，還需要重新建立判斷矩陣。根據(jù)SAATY等[6-7]的一致性判斷準(zhǔn)則，每個判斷矩陣都用一個一致性比值(Consistency ratio,CR)來判斷矩陣一致性的程度。

如果一致性比值滿足式(2)，則一致性可以接受

(2)

表3 隨機(jī)一致性指標(biāo)R參考值

本文研究的迫擊炮炮口擾動，分別為45°、60°、80°3種不同射角工況，n=3。參考表2，并參考相關(guān)手冊，對各工況的重要程度進(jìn)行兩兩對比，得到判斷矩陣A

(3)

求得該矩陣的最大特征值為3.04，由n=3可得R=0.58，根據(jù)式(2)算得該矩陣的一致性比值為0.034<0.1，判定該矩陣的一致性可以接受，不需要重新建立判斷矩陣。最大特征值對應(yīng)的矢量為(0.260,0.634,0.106)T，因此45°、60°、80° 3種不同射角工況的權(quán)重系數(shù)分別為0.260、0.634、0.106。

2.3 Kriging模型

在迫擊炮炮口擾動問題的研究中，炮箍位置與炮口擾動之間存在著一種復(fù)雜的非線性映射關(guān)系，鑒于無法通過構(gòu)建普通函數(shù)關(guān)系式來表達(dá)，選用Kriging模型[8-9]反映這種關(guān)系。通過自身這種非常強(qiáng)大的非線性映射能力，Kriging模型可以在不參考任何經(jīng)驗公式的情況下，通過自主學(xué)習(xí)或訓(xùn)練，自動總結(jié)出數(shù)據(jù)之間的函數(shù)關(guān)系，是一種常用且有效的建模手段。Kriging模型假設(shè)響應(yīng)值與樣本點之間的關(guān)系式為

y(x) = f(x)+z(x)

(4)

其中：y(x)是待擬合的響應(yīng)函數(shù)；f(x)是多項式函數(shù)的響應(yīng)模型，相當(dāng)于對設(shè)計空間的一個全局模擬。一般情況下，將f(x) 簡化為一常數(shù)β。z(x)是一個隨機(jī)函數(shù)，通過樣本點插值獲得的局部偏差。

建立Kriging模型需要一系列訓(xùn)練樣本。適當(dāng)?shù)挠?xùn)練樣本數(shù)量以及合理的分布能夠使Kriging模型更加確切地表達(dá)出結(jié)構(gòu)的映射關(guān)系。因此，在45°、60°、80° 3種不同射角工況下，分別在設(shè)計變量的取值區(qū)間內(nèi)均勻地取34個值，構(gòu)成總數(shù)為102的輸入樣本。根據(jù)樣本的選擇對第1節(jié)所建立的迫擊炮動力學(xué)模型做出相應(yīng)的修改并進(jìn)行動力學(xué)分析，從而獲得3種不同射角下的迫擊炮動態(tài)響應(yīng)。本文優(yōu)化設(shè)計的主要目標(biāo)是降低迫擊炮炮口擾動，而彈丸出炮口時刻的炮口角位移和角速度可以作為衡量炮口擾動的主要指標(biāo)。通過提取出彈丸出炮口時刻的炮口回轉(zhuǎn)角位移θy、高低角位移θz、回轉(zhuǎn)角速度ωy和高低角速度ωz作為輸出樣本，可以獲得Kriging模型試驗樣本庫所需的輸入和輸出樣本。

利用Lophaven等[10]開發(fā)的DACE Matlab Kriging Toolbox工具箱，對所獲得的102組訓(xùn)練樣本進(jìn)行訓(xùn)練，在Matlab上構(gòu)建了12個Kriging模型，分別模擬不同射角下設(shè)計變量與彈丸出炮口時刻的炮口回轉(zhuǎn)角位移θy、高低角位移θz、回轉(zhuǎn)角速度ωy和高低角速度ωz之間的非線性映射關(guān)系。

經(jīng)過訓(xùn)練后的Kriging模型還需要檢驗其泛化能力和預(yù)測精度。在3種不同射角下，分別通過在設(shè)計變量的取值范圍內(nèi)均勻地選取10組數(shù)據(jù)，進(jìn)行相應(yīng)的動力學(xué)分析獲得炮口響應(yīng)量，從而得到30組檢測樣本。對于檢測樣本的數(shù)據(jù)擬合度，Kriging模型可采用復(fù)相關(guān)系數(shù)R2來檢驗，其表達(dá)式為

(5)

表4 檢驗結(jié)果

由表4可知，Kriging模型預(yù)測的復(fù)相關(guān)系數(shù)均接近1，表明構(gòu)建的Kriging模型具有較高的預(yù)測精度和較好的泛化能力。

2.4 優(yōu)化目標(biāo)函數(shù)

在單一射角工況下，對迫擊炮炮箍位置進(jìn)行結(jié)構(gòu)優(yōu)化時，由于優(yōu)化目標(biāo)函數(shù)同時包含彈丸出炮口時的炮口角位移和角速度，通過采用線性加權(quán)的方法構(gòu)建目標(biāo)函數(shù)，將多目標(biāo)優(yōu)化問題轉(zhuǎn)化為單目標(biāo)優(yōu)化問題[11]。在對目標(biāo)函數(shù)進(jìn)行歸一化處理后，所建立的優(yōu)化目標(biāo)函數(shù)為

(6)

其中：θy0、θz0、ωy0、ωz0分別為設(shè)計變量取初始值時彈丸出炮口瞬間的炮口回轉(zhuǎn)角位移、高低角位移、回轉(zhuǎn)角速度、高低角速度；α、β為加權(quán)系數(shù)，根據(jù)課題組前期靈敏度分析的經(jīng)驗，分別取為0.3、0.7。

本文研究的是多工況結(jié)構(gòu)優(yōu)化，結(jié)合上文45°、60°、80° 3種不同射角工況的權(quán)重系數(shù)，最終優(yōu)化目標(biāo)函數(shù)為

Oobj=aFa+bFb+cFc

(7)

其中：a、b、c分別為45°、60°、80° 3種射角所對應(yīng)的權(quán)重系數(shù)；Fa、Fb、Fc則分別為相應(yīng)工況下對應(yīng)式(6)的優(yōu)化目標(biāo)函數(shù)。

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

傳統(tǒng)的優(yōu)化方法無法針對基于Kriging模型建立的非線性函數(shù)關(guān)系尋優(yōu)，而作為一種不需要具體函數(shù)形式，可以通過模擬自然進(jìn)化過程搜索全局最優(yōu)解的方法，遺傳算法[12]可以較好地實現(xiàn)迫擊炮炮箍位置優(yōu)化。

以目標(biāo)函數(shù)Oobj最小化為目標(biāo)，運用遺傳算法在設(shè)計變量取值范圍內(nèi)優(yōu)化計算，得到結(jié)構(gòu)參數(shù)優(yōu)化結(jié)果為336 mm?；趦?yōu)化后的炮箍位置，對迫擊炮動力學(xué)模型做出相應(yīng)的修改并進(jìn)行有限元仿真分析，并將結(jié)果與優(yōu)化前進(jìn)行對比，如表5。

表5 優(yōu)化結(jié)果

由表5可以看出，采用Kriging模型和遺傳算法優(yōu)化后，最終目標(biāo)函數(shù)Oobj減少了46.02%。此外也可以看出，在80°射角工況下高低角速度ωz經(jīng)過優(yōu)化后增加了2.93%，但是該工況下的大部分?jǐn)?shù)值都得到了一定減少。同時，重要程度考慮最多的60°射角工況下的炮口擾動經(jīng)過優(yōu)化后，彈丸出炮口時刻的炮口回轉(zhuǎn)角位移θy、高低角位移θz、回轉(zhuǎn)角速度ωy和高低角速度ωz分別得到了37.56%、45.82%、58.98%、45.47%的顯著減少。除此之外，45°射角工況下的炮口擾動經(jīng)過優(yōu)化后也得到了的有效減小。優(yōu)化結(jié)果表明，在多工況條件下優(yōu)化后的迫擊炮炮箍位置，可以有效減少炮口擾動，實現(xiàn)迫擊炮炮口擾動優(yōu)化最小的目標(biāo)。

4 結(jié)論

本文主要研究了多工況下迫擊炮炮箍位置對炮口擾動情況的影響并進(jìn)行了優(yōu)化設(shè)計。研究結(jié)果表明優(yōu)化方法是行之有效的，可以為迫擊炮的炮箍位置設(shè)計提供理論參考。在接下來的研究工作中，進(jìn)一步增加對多工況問題的研究，以提高結(jié)構(gòu)優(yōu)化的普遍性。

[1] 唐治.迫擊炮設(shè)計[M].北京：兵器工業(yè)出版社,1994：106-108.

[2] 鄔欣明.某車載自動迫擊炮全炮動力學(xué)仿真[D].南京：南京理工大學(xué),2008.

[3] 陳杰,楊國來,葛建立.基于正交試驗的火炮上架多工況結(jié)構(gòu)優(yōu)化設(shè)計[J].彈道學(xué)報，2015,27(1)：92-96.

[4] 蘭鳳崇,賴番結(jié),陳吉清,等.考慮動態(tài)特性的多工況車身結(jié)構(gòu)拓?fù)鋬?yōu)化研究[J].機(jī)械工程學(xué)報,2014,50(20)：122-128.

[5] 李小剛,程錦,劉振宇,等.基于雙層更新Kriging模型的機(jī)械結(jié)構(gòu)動態(tài)特性穩(wěn)健優(yōu)化設(shè)計[J].機(jī)械工程學(xué)報,2014,50(3)：165-173.

[6] ATTY T L,VARGAS L G.Models,methods,concepts & applications of the analytic hierarchy process[M].2nd ed.New York：Springer，2012.

[7] SUBRAMANIAN N,RAMANATHAN R.A review of applications of analytic hierarchy process in operations management[J].International Journal of Production Economics,2012,138(2)：215-241.

[8] 劉克龍,姚衛(wèi)星,穆雪峰.基于Kriging代理模型的結(jié)構(gòu)形狀優(yōu)化方法研究[J].計算力學(xué)學(xué)報,2006,23(3)：344-347.

[9] 王紅濤,竺曉程,杜朝輝.自適應(yīng)Kriging近似模型及其在二維擴(kuò)壓器優(yōu)化設(shè)計中的應(yīng)用[J].計算力學(xué)學(xué)報,2011,28(1)：15-19.

[10]LOPHAVEN S N,NIELSEN H B,SONDERGAARD J.DACE—A Matlab Kriging Toolbox version 2.0[EB/OL].[2012-02-23] http://www2.imm.dtu.dk/～hbn/dace/dace.pdf.

[11]梁傳建,楊國來,王曉鋒.基于神經(jīng)網(wǎng)絡(luò)和遺傳算法的火炮結(jié)構(gòu)動力學(xué)優(yōu)化[J].兵工學(xué)報,2015,36(5)：789-794.

[12]周樂,楊國來,葛建立,等.基于遺傳算法的火炮反后坐裝置結(jié)構(gòu)多目標(biāo)優(yōu)化研究[J].兵工學(xué)報,2015,36(3)：433-436.

[13]徐志遠(yuǎn)，葛建立，楊國來.影響炮口擾動的火炮總體結(jié)構(gòu)參數(shù)靈敏度分析與優(yōu)化[J].兵器裝備工程學(xué)報，2016(6):45-48.

(責(zé)任編輯周江川)

Optimization Design for Gun Hoop Position of Mortar

LIN Zhi-wen1, GE Jian-li1, ZHANG Hong-hao2, WANG Xue-yan1, YANG Guo-lai1

(1.School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China;2.Beijing Institute of Special Electromechanical Technology, Beijing 100012, China)

Aiming at the problem of mortar muzzle disturbance, the structural dynamics model of a 120 mm mortar was established, and the dynamic optimization of mortar structure was carried out by using the combination of the finite element method, the analytic hierarchy process, the Kriging model and the genetic algorithm. Based on the test data for the training sample, a nonlinear mapping relationship between the gun hoop position of mortar and muzzle disturbance was established. By constructing the optimization objective function, the objective function was solved using genetic algorithm and the minimum quantity optimization of mortar muzzle disturbance was realized. Research shows that the muzzle disturbance could be effectively reduced at the gun hoop position which is optimized under multiple conditions, and it has certain reference significance for the study on mortar muzzle disturbance problem.

mortar; analytic hierarchy process; Kriging model; structural dynamics optimization

2016-09-04；

2016-09-28

中央高校基本科研業(yè)務(wù)費專項資金(30915118825)；國家重大科學(xué)儀器設(shè)備開發(fā)專項(2013YQ470765)

林智文(1992—)，男，碩士研究生，主要從事非線性有限元仿真研究。

葛建立(1980—)，男，博士，副教授，主要從事非線性有限元、虛擬樣機(jī)以及等幾何分析研究。

10.11809/scbgxb2017.01.019

林智文,葛建立,張鴻浩,等.迫擊炮炮箍位置優(yōu)化設(shè)計[J].兵器裝備工程學(xué)報，2017(1):79-82.

format:LIN Zhi-wen, GE Jian-li, ZHANG Hong-hao, et al.Optimization Design for Gun Hoop Position of Mortar[J].Journal of Ordnance Equipment Engineering,2017(1):79-82.

TJ31

A