基于有限元剛?cè)狁詈蠀?shù)化建模的火炮射擊穩(wěn)定性研究

趙躍躍，顧克秋，焦 偉

(南京理工大學(xué) 機械工程學(xué)院， 南京 210094)

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

基于有限元剛?cè)狁詈蠀?shù)化建模的火炮射擊穩(wěn)定性研究

趙躍躍，顧克秋，焦 偉

(南京理工大學(xué) 機械工程學(xué)院， 南京 210094)

將對射擊穩(wěn)定性影響較大的部件考慮為柔性體，建立了全炮剛?cè)狁詈戏蔷€性有限元參數(shù)化模型。通過對射擊穩(wěn)定性的總體結(jié)構(gòu)參數(shù)靈敏度分析，研究了重要結(jié)構(gòu)參數(shù)對射擊穩(wěn)定性的影響規(guī)律。同時利用響應(yīng)面(RSM)近似模型方法進行優(yōu)化，有效降低了前大架跳高。該研究對火炮總體結(jié)構(gòu)設(shè)計與優(yōu)化具有借鑒意義。

火炮；射擊穩(wěn)定性；剛?cè)狁詈希混`敏度分析；RSM近似模型

火炮系統(tǒng)輕量化符合現(xiàn)代戰(zhàn)爭發(fā)展的需要，射擊穩(wěn)定性是火炮系統(tǒng)輕量化必須考慮的核心性能指標(biāo)。射擊穩(wěn)定性差已成為制約超輕型火炮發(fā)展的主要因素之一。目前，對火炮射擊穩(wěn)定性問題的研究方法主要有測試、多體動力學(xué)理論分析和有限元分析三種[1],很多學(xué)者在這些方面做了大量的研究工作[2-4]。但大多數(shù)的研究還是采用多剛體動力學(xué)模型[5-7]，忽略了關(guān)鍵部件彈性變形對射穩(wěn)定性的影響。即使建立了剛?cè)狁詈蟿恿W(xué)模型進行優(yōu)化研究，也只是對總體參數(shù)中剛性部件結(jié)構(gòu)參數(shù)進行靈敏度分析與優(yōu)化[8-10]。主要原因在于復(fù)雜柔性體建模難度大不易實現(xiàn)，對柔性體結(jié)構(gòu)參數(shù)的優(yōu)化還沒找到一種行之有效的方法。

為得到更接近實際的某輕型火炮射擊穩(wěn)定性分析結(jié)果，把對射擊穩(wěn)定性影響較大的前后大架、身管、搖架和活動上架用柔性體模擬，建立全炮剛?cè)狁詈戏蔷€性有限元參數(shù)化模型，實現(xiàn)自動參數(shù)化建模，對包括柔性體結(jié)構(gòu)參數(shù)在內(nèi)的火炮總體結(jié)構(gòu)參數(shù)進行靈敏度分析，確定對射擊穩(wěn)定性有重要影響的結(jié)構(gòu)參數(shù)，探索研究這些結(jié)構(gòu)參數(shù)對射擊穩(wěn)定性的影響規(guī)律，然后利用響應(yīng)面近似模型(RSM)方法構(gòu)建近似模型，采用粒子群算法(PSO)對射擊穩(wěn)定性進行優(yōu)化，從而為提高火炮系統(tǒng)射擊穩(wěn)定性提供參考。

1 某輕型炮剛?cè)狁詈戏蔷€性有限元建模

1.1 建模說明

根據(jù)火炮實際發(fā)射過程，建模說明如下：

1) 將全炮劃分為炮口制退器、身管、炮尾、左右前大架、左右后大架、下架、活動上架、上架底板、搖架、搖架襯瓦、左右高平機共14個部分，考慮身管、搖架、活動上架和前后大架彈性變形對射擊穩(wěn)定性的影響，采用柔性體建模，其余部件均為剛體。

2) 考慮身管與搖架襯瓦接觸間隙的影響；

3) 炮尾滑塊與搖架導(dǎo)軌在后坐過程中接觸/碰撞；

4) 火炮放置于水平地面，在0°方向角0°射角工況下發(fā)射，下部架體與土壤之間的相互作用采用質(zhì)量-彈簧-阻尼單元的集中參數(shù)模型模擬。

1.2 火炮下部架體-土壤相互作用模型

某輕型火炮在發(fā)射過程中，其下架和左右前大架與土壤貼合，后大架埋入土壤中。由于下部架體與土壤之間的相互作用十分復(fù)雜，建立土壤三維實體模型計算量巨大，這會增加全炮有限元計算代價甚至導(dǎo)致不收斂，可以采用質(zhì)量-彈簧-阻尼單元的集中參數(shù)模型來模擬。結(jié)構(gòu)-土壤作用原理如圖1所示，對于長為L寬為b的一個矩形基礎(chǔ)，可將其轉(zhuǎn)換為等效半徑為r0的圓形基礎(chǔ)進行求解，彈性半無限空間的剛度和阻尼的計算公式如表1所示[11]。

利用集中參數(shù)模型，前大架、下架與土壤之間只考慮豎直方向上的剛度與阻尼，后大架與土壤之間考慮豎直、前后、水平三個方向的剛度與阻尼，不考慮扭轉(zhuǎn)剛度和轉(zhuǎn)動剛度。根據(jù)表1中的計算方法，分別計算出下架、前大架和后大架與土壤相互作用的剛度和阻尼值，如表2所示。

圖1 結(jié)構(gòu)-土壤作用原理

振動類型剛度阻尼豎直振動4Gdr01-μ3．4r201-μρGd水平振動32(1-μ)7-8μGdr018．4(1-μ)r207-8μρGd

表1中：Gd=Ed/2(1+μ)為動剪切模量，Ed為土壤動彈性模量；r0為等效圓形的半徑，r0=(bL/π)1/2；ρ為土壤密度，μ為泊松比。

表2 集中參數(shù)模型的剛度值和阻尼值

1.3剛?cè)狁詈戏蔷€性有限元模型

建模時將關(guān)鍵部件身管、搖架、活動上架采用柔性體建模，同時考慮前后大架對全炮支撐和射擊穩(wěn)定性的影響，用簡化的梁單元模擬，其余均為剛體，在其質(zhì)心處賦予相應(yīng)的質(zhì)量和轉(zhuǎn)動慣量。身管與搖架襯瓦、炮尾滑塊與搖架導(dǎo)軌建立接觸；搖架與活動上架、活動上架與上架底板、上下座圈之間的連接用施加等效剛度和阻尼的Cartesion-Cardan復(fù)合連接器模型；簡化為下部架體與土壤之間的相互作用簡化為集中參數(shù)模型；高平機簡化為施加等效剛度和阻尼的軸向連接器；復(fù)進機、制退機和緩沖器的作用力通過對有限元軟件進行二次開發(fā)實現(xiàn)。

模型求解分兩步進行，第一步靜力分析求得全炮在重力作用下的初始狀態(tài)；第二步模擬火炮發(fā)射過程，采用隱式求解，不僅延續(xù)重力作用，同時身管施加炮膛合力。全炮剛?cè)狁詈戏蔷€性有限元模型如圖2所示。

1.4 有限元模型參數(shù)化

選取一組基于總體結(jié)構(gòu)的可變參量，通過對參數(shù)值的改變自動實現(xiàn)相應(yīng)設(shè)計方案，這種方法叫作總體結(jié)構(gòu)參數(shù)化?；谟邢拊幊坦ぞ逷ython編寫腳本程序，實現(xiàn)全炮剛?cè)狁詈戏蔷€性有限元模型的參數(shù)化。

圖2 全炮剛?cè)狁詈戏蔷€性有限元模型

2 射擊穩(wěn)定性總體參數(shù)靈敏度分析

2.1 設(shè)計變量的選取和目標(biāo)函數(shù)的確定

為探索各結(jié)構(gòu)參數(shù)對射擊穩(wěn)定性的影響，選取前大架長度(前大架與下架連接位置不變)、后大架長度(后大架駐鋤支點位置不變)和后大架等效剛度、下耳軸前后位置和高度(相對全炮坐標(biāo)系原點)、高平機下支點距活動上架下耳軸前后距離和豎直距離、炮尾質(zhì)量及質(zhì)心位置、下架質(zhì)量及質(zhì)心位置和高平機等效剛度共16個結(jié)構(gòu)參數(shù)為設(shè)計變量。各設(shè)計變量的初值及取值范圍如表3所示。

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

根據(jù)某輕型火炮實際發(fā)射動作特征，以前大架跳高即前大架跳離地面的距離作為衡量射擊穩(wěn)定性好壞的依據(jù)，故目標(biāo)函數(shù)為：

f=Hqdj

(1)

式中，Hqdj即為前大架跳高。

2.2 結(jié)構(gòu)參數(shù)靈敏度分析

由于最優(yōu)拉丁超立方設(shè)計會使所有的試驗點盡量均勻分布在設(shè)計空間，具有非常好的空間填充性和均衡性，使得因子和響應(yīng)的擬合更加精確真實，因而采用Isight中的最優(yōu)拉丁超立方試驗設(shè)計對火炮總體結(jié)構(gòu)參數(shù)進行靈敏度分析。試驗設(shè)計共選取300個試驗點，通過試驗設(shè)計得到對目標(biāo)函數(shù)貢獻程度前10位的設(shè)計變量如圖3所示。

圖3 各設(shè)計變量對目標(biāo)函數(shù)貢獻程度

2.3主要結(jié)構(gòu)參數(shù)對射擊穩(wěn)定性的影響規(guī)律

通過靈敏度分析可知，下耳軸高度、后大架等效剛度、下耳軸前后方向位置、炮尾質(zhì)量以及炮尾質(zhì)心Y方向坐標(biāo)、后大架長度、下架質(zhì)量對射擊穩(wěn)定性影響最為顯著，其對射擊穩(wěn)定性的影響規(guī)律如圖4所示，由圖可知：

1)Hqdj隨著下耳軸高度XEZ2的增大而增大，這是因為下耳軸高度增大使得翻轉(zhuǎn)力臂隨之增大，不利于提高射擊穩(wěn)定性；Hqdj隨著下耳軸前后位置XEZ1的增大而減小，這是因為下耳軸前移使得全炮重心也隨之前移，增大了穩(wěn)定力矩，有利于提高射擊穩(wěn)定性。

2)Hqdj隨著后大架等效剛度K1的增大而減小，這說明后大架材料剛度的增大有利于提高射擊穩(wěn)定性；Hqdj隨著后大架長度L2的增大呈線性增大的趨勢，即在保持后大架駐鋤位置不變的情況下，后大架不宜過長。

3) 炮尾質(zhì)心Y方向坐標(biāo)HZT3和炮尾質(zhì)量HZT1對射擊穩(wěn)定性的影響規(guī)律相同，Hqdj都隨它們的增大而減小。炮尾質(zhì)心Y方向坐標(biāo)值HZT3增大有利于提高射擊穩(wěn)定性，這是由于HZT3的改變直接影響動力偶臂的長短。炮尾質(zhì)量HZT1增大使得后坐體質(zhì)心后移，后坐體質(zhì)心后移也有利于提高射擊穩(wěn)定性。但Hqdj隨著HZT3和HZT1的增大下降緩慢，說明小范圍的改變炮尾質(zhì)量和質(zhì)心Y方向坐標(biāo)對射擊穩(wěn)定性的提高是有限的。

4) 增大下部架體質(zhì)量XBJT1在一定程度上可以提高射擊穩(wěn)定性。

圖4 主要結(jié)構(gòu)參數(shù)對目標(biāo)函數(shù)的主效應(yīng)圖

3 射擊穩(wěn)定性優(yōu)化

3.1 構(gòu)造近似模型

近似模型方法是通過數(shù)學(xué)模型的方法逼近一組輸入變量(獨立變量)與輸出變量(響應(yīng)變量)的方法。在上述試驗設(shè)計的基礎(chǔ)上構(gòu)造出代替原始全炮剛?cè)狁詈蠀?shù)化模型的近似模型，可以避免成千甚至上萬次的尋優(yōu)計算。在近似模型中輸入變量與輸出響應(yīng)之間關(guān)系為：

(2)

利用試驗設(shè)計分析結(jié)果，選取靈敏度分析篩選出的影響程度較大的K1、XEZ1、L2、HZT1、HZT3和XJ1為輸入變量，以前大架跳高Hqdj為輸出響應(yīng)。因XEZ2對射擊穩(wěn)定性的影響最為敏感，且其值由于實際情況制約無需調(diào)整，為避免在優(yōu)化時對其他設(shè)計變量最優(yōu)值的影響，故不作為輸入變量。同時，為了得到最為精確的近似模型，分別利用RSM(響應(yīng)面)、RBF(神經(jīng)網(wǎng)絡(luò))和Kriging(克里格)三種方法構(gòu)造近似模型，通過計算，提取誤差分析R2值分別為0.996 48、0.985 47和0.968 78，均大于0.9。由于R2值越接近1可信度越高，故選用RSM近似模型。

3.2 優(yōu)化的實施

選取近似模型中的輸入變量作為優(yōu)化設(shè)計變量，目標(biāo)函數(shù)為前大架跳高最小。同時，為了滿足全炮靜止穩(wěn)定性，不能使得全炮在靜止條件下向前傾覆，還需設(shè)置約束條件。故優(yōu)化數(shù)學(xué)模型為：

(3)

式中，Hqdj(xi)為目標(biāo)函數(shù)；Hzc為全炮在重力作用下后大架駐鋤豎直方向的位移，當(dāng)Hzc≥0時，說明全炮向前傾覆，不滿足靜態(tài)穩(wěn)定性；xi為設(shè)計變量，xLi和xUi分別為設(shè)計變量的上下限。

優(yōu)化采用粒子群優(yōu)化算法，該算法具有較強的魯棒性和全局搜索能力，非常適用于求解復(fù)雜優(yōu)化問題，是一種高效的并行搜索優(yōu)化算法。算法主要參數(shù)設(shè)置如表4所示。

表4 粒子群(PSO)算法參數(shù)

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

優(yōu)化后的結(jié)果如表5所示，優(yōu)化前后的前大架跳高隨時間變化曲線如圖5所示。

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

圖5 前大架跳高

從圖5可以看出，優(yōu)化后的前大架跳高明顯小于優(yōu)化前，跳高峰值減小了33.26%，大大提高了射擊穩(wěn)定性，達到了優(yōu)化的目的。值得注意的是，由于是單目標(biāo)優(yōu)化，優(yōu)化后的結(jié)構(gòu)參數(shù)雖然使得射擊穩(wěn)定性得到了提高，但同時也會帶來其他問題。例如后大架長度的改變使得下架與后大架連接角度的改變，從而導(dǎo)致后大架受力的變化，而后大架受力較為惡劣，容易出現(xiàn)局部應(yīng)力過大。因此，為使火炮綜合性能最優(yōu)，還需進一步研究。

4 結(jié)論

下耳軸位置、后大架長度和后大架等效剛度、炮尾質(zhì)量及質(zhì)心位置以及下架質(zhì)量對射擊穩(wěn)定性有重要影響，需合理取值。

本文的研究過程對某輕型火炮總體結(jié)構(gòu)優(yōu)化有一定的參考價值，這種研究方法對其他工程實際問題也有一定的借鑒意義。

[1] 劉露陽，楊國來，肖輝.某輕型牽引榴彈炮射擊穩(wěn)定性分析[J].動力學(xué)與控制學(xué)報，2015，13(1)：37-41.

[2] WEBB D W，SOENCKSEN GUIDOS B J，GARNER J W.Flatfire jump performance of a 155mm-M198 howitzer[J].ADA369710，1999：1-15.

[3] 趙紀(jì)華.火炮最優(yōu)射擊穩(wěn)定性設(shè)計研究[J].振動與沖擊，2010，29(11)：91-94.

[4] 彭迪.超輕型火炮結(jié)構(gòu)動力學(xué)有限元分析與總體優(yōu)化匹配技術(shù)研究[D].南京：南京理工大學(xué)，2012.

[5] CHEN M.Research on influence of structure parameter on firing stability of machine gun based on the virtual prototype[J].ActaArmamentarii，2008，29(10)：1167-1171.

[6] 萬李，徐亞棟.某輕型火炮結(jié)構(gòu)參數(shù)靈敏度分析及穩(wěn)定性優(yōu)化[J].四川兵工學(xué)報，2012，33(8)：34-36.

[7] 李三群，劉海平，鄭智剛等.火炮射擊穩(wěn)定性的序列二次規(guī)劃算法與虛擬樣機融合優(yōu)化[J].武器裝備自動化，2005，24(3)：1-5.

[8] 王飛，楊國來，葛建立等.某火炮射擊穩(wěn)定性的總體參數(shù)靈敏度分析與優(yōu)化[J].火炮發(fā)射與控制學(xué)報，2014，35(4)：58-61.

[9] 周樂.大口徑火炮反后坐裝置多目標(biāo)優(yōu)化設(shè)計研究[D].南京：南京理工大學(xué)，2015.

[10] 葛建立，楊國來，曾晉春.某自行火炮總體結(jié)構(gòu)參數(shù)靈敏度分析與優(yōu)化[J].火炮發(fā)射與控制學(xué)報，2007(1)：16-19.

[11] 常士驃，張?zhí)K民.工程地質(zhì)手冊第四版[M].北京：中國建筑工業(yè)出版社，2008.

ResearchonGunFiringStabilityBasedonFiniteElementRigid-FlexibleCouplingParametricModeling

ZHAO Yueyue， GU Keqiu， JIAO Wei

(Mechanical Engineering College， Nanjing University of Science and Technology， Nanjing 210094， China)

The components that have a greater impact on firing stability were considered as flexible body, and the nonlinear finite element parametric model of the whole gun was established. Through the analysis of the sensitivity of the overall structural parameters of the firing stability, the structural parameters which have influence rule on the firing stability were studied. At the same time, the RSM approximation model was used to optimize the firing stability, which effectively reduced the high jump of the front trails. This study was of great significance to the design and optimization of the overall structure of gun.

gun；firing stability；rigid-flexible coupling；sensitivity analysis；RSM approximation model

2017-06-01；

2017-06-29

武器裝備重點預(yù)先研究項目(301070504)

趙躍躍(1991—)，男，碩士研究生，主要從事動力學(xué)仿真與優(yōu)化技術(shù)研究。

10.11809/scbgxb2017.10.015

本文引用格式：趙躍躍，顧克秋，焦偉.基于有限元剛?cè)狁詈蠀?shù)化建模的火炮射擊穩(wěn)定性研究[J].兵器裝備工程學(xué)報，2017(10):69-73.

formatZHAO Yueyue，GU Keqiu，JIAO Wei.Research on Gun Firing Stability Based on Finite Element Rigid-Flexible Coupling Parametric Modeling[J].Journal of Ordnance Equipment Engineering,2017(10):69-73.

TJ303

A

2096-2304(2017)10-0069-05

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