一種自行火炮火力系統(tǒng)集成優(yōu)化設(shè)計方法

肖曉，徐誠，徐亞棟

(南京理工大學 機械工程學院，江蘇 南京 210094)

?

一種自行火炮火力系統(tǒng)集成優(yōu)化設(shè)計方法

肖曉，徐誠，徐亞棟

(南京理工大學 機械工程學院，江蘇 南京 210094)

摘要：為了提高某自行火炮的綜合性能，集成了對火炮射擊起直接影響作用的全彈道、炮身/反后坐裝置、調(diào)炮控制系統(tǒng)在內(nèi)的分析模型，建立了火力系統(tǒng)多學科協(xié)同仿真模型。確定了主要的優(yōu)化設(shè)計變量，以殺傷面積最大、反后坐阻力最小及火炮伺服系統(tǒng)調(diào)整時間最短作為總體的優(yōu)化目標，建立了火力系統(tǒng)集成優(yōu)化模型。采取多目標遺傳算法進行尋優(yōu)求解，獲得了系統(tǒng)的Parato最優(yōu)解集，根據(jù)實際需求，選取一組優(yōu)化解。優(yōu)化結(jié)果表明，模型的總體優(yōu)化目標得到了較大的改善，為自行火炮火力系統(tǒng)一體化設(shè)計和提高火炮綜合性能提供了有效的設(shè)計方法。

關(guān)鍵詞：火炮；全彈道；反后坐；多目標；多學科優(yōu)化

自行火炮的火力系統(tǒng)是火炮的重要組成部分之一，是決定火炮射擊性能的關(guān)鍵部件。火力系統(tǒng)主要包含彈藥、彈道、炮身/反后坐裝置及自動調(diào)炮等部分。彈藥與彈道是指內(nèi)彈道、外彈道以及彈藥終點效能，涉及到身管內(nèi)部結(jié)構(gòu)參數(shù)，彈丸結(jié)構(gòu)參數(shù)，裝藥參數(shù)，是決定火炮的射擊性能參數(shù)(如最大射程、射擊精度、彈丸飛行時間等)的直接因素；炮身/反后坐裝置是火炮上重要部件之一，在火炮射擊時提供制動力控制后坐部分的后坐運動，并使之復位。其性能直接影響著火炮射擊穩(wěn)定性。高低和方向伺服系統(tǒng)是火炮操控重要組成部分，決定了火炮的調(diào)炮時間和精度，對火炮的整個反應速度以及射擊精度起著決定性的影響。關(guān)于火炮優(yōu)化已有不少研究。李克婧等[1]對火炮全彈道過程進行了仿真分析；宗世增等[2]對火炮反后坐裝置進行了動力學耦合分析以及優(yōu)化；洪亞軍等[3]將火炮身管與反后坐裝置模型進行了集成優(yōu)化研究；李銀伢等[4]給出了一種基于粒子群優(yōu)化算法的火炮伺服系統(tǒng)控制器優(yōu)化設(shè)計方法。但是目前對火力系統(tǒng)集成優(yōu)化研究較少。本文將火力系統(tǒng)集成為一體，充分考慮了各部分之間的數(shù)據(jù)交換、耦合關(guān)系，建立了火力系統(tǒng)協(xié)同仿真模型和多學科優(yōu)化模型。在滿足約束條件下，以殺傷面積最大、反后坐阻力最小以及伺服系統(tǒng)穩(wěn)定時間最短為優(yōu)化目標，尋求系統(tǒng)的最優(yōu)解，為自行火炮火力系統(tǒng)一體化設(shè)計提供一種新的方法。

1火力系統(tǒng)多學科協(xié)同仿真模型

1.1　全彈道模型

a) 內(nèi)彈道計算模型

采用經(jīng)典內(nèi)彈道計算模型，模型如下：

式中，i=1,2,L,n，ψi為第i種火藥已燃百分數(shù)；χi、λi和μi為第i種火藥的藥型參數(shù)；Zi為第i種火藥的相對已燃厚度；μli為第i種火藥的燃速系數(shù)；ni為第i種火藥的燃速指數(shù)；p為火炮膛內(nèi)壓力；φ為次要功系數(shù)；m為彈丸質(zhì)量；v為彈丸速度；S為火炮身管橫截面積；l為彈丸行程長；fi為第i種火藥的火藥力;ωi為第i種火藥的質(zhì)量。

b) 外彈道計算模型

針對某型自行火炮，采用質(zhì)點外彈道模型，編寫程序求解微分方程，求解出榴彈射擊時的落角、最大射程以及飛行時間等參量，計算模型如下:

式中，x、y是彈丸在每一時刻對應的坐標；v為彈丸的速度，vx、vy是彈丸在x、y方向的分速度；c是彈道系數(shù)；H(y)是氣重函數(shù)；G(v)是阻力函數(shù)；t是彈丸的飛行時間。

c) 終點彈道計算模型

在給定彈丸結(jié)構(gòu)參數(shù)(圖1)，一定的落角和落速時計算榴彈的殺傷面積。殺傷準則選用A-S殺傷準則，針對選用的炸藥類型和目標姿勢，對榴彈殺傷面積進行積分。

圖1　榴彈結(jié)構(gòu)參數(shù)圖

計算模型如下：

(1)

(2)

(3)

式中：Ni為第i組內(nèi)破片數(shù)目；Phk,i為該組內(nèi)破片的平均條件殺傷概率。

式(2)中的as(x,y)為殺傷破片的球面分布度，如式(4)所示：

(4)

1.2　反后坐裝置模型

本文所研究的火炮的反后坐裝置是由節(jié)制桿式駐退機和液體氣壓式復進機組成。炮身/反后坐裝置計算模型如下：

式中：Fpt為炮膛合力，F(xiàn)R為后坐阻力由液壓阻力Fφh、復進機力Ff、密封裝置摩擦力Fm、搖架導軌摩擦力FT及反后坐部分重力分量等構(gòu)成，mh為后坐部分質(zhì)量，K1、K2分別為主流和支流的液壓阻力系數(shù)；ρ為駐退液密度；A0為駐退機活塞工作面積；Ap為節(jié)制環(huán)孔面積；ax為節(jié)制桿任意截面的流液孔面積；Afj為復進節(jié)制器工作面積；A1為支流最小截面積；v為后坐速度；Af為復進機活塞工作面積；pf0為氣體初壓；Vf0為初始氣體體積；S為后坐行程。

以節(jié)制桿的尺寸參數(shù)、內(nèi)彈道的膛壓合力為輸入，求解微分方程，計算出火炮反后坐裝置的最大后坐阻力Frmax。

1.3　火炮高低和方向伺服系統(tǒng)模型

根據(jù)火炮的具體要求，高低機和方向機用交流伺服電機驅(qū)動，通過PID算法實現(xiàn)了三環(huán)控制，構(gòu)成了位置伺服系統(tǒng)的典型結(jié)構(gòu)，內(nèi)環(huán)為電流環(huán)和速度環(huán)，外環(huán)為位置環(huán)，實現(xiàn)了系統(tǒng)調(diào)炮過程定位的準確性和快速性， 子系統(tǒng)模型如圖2所示。

圖2　火炮伺服系統(tǒng)

對于本模型采用Z變換進行離散化，以C#語言對整個離散的系統(tǒng)傳遞函數(shù)進行編程，以電機參數(shù)及PID參數(shù)為輸入?yún)?shù)，通過調(diào)節(jié)PID參數(shù)，使系統(tǒng)達到穩(wěn)定，得到系統(tǒng)的最大超調(diào)率，精度跟穩(wěn)定時間等性能參數(shù)。

1.4　火力系統(tǒng)多學科協(xié)同仿真模型

通過軟件接口，建立了彈道仿真模型、反后坐仿真模型以及伺服系統(tǒng)模型的多學科協(xié)同仿真模型，實現(xiàn)了各模型之間的數(shù)據(jù)交流，模型如圖3所示。

圖3　火力系統(tǒng)多學科協(xié)同仿真模型

采用上述單學科分析模型和火力系統(tǒng)多學科協(xié)同仿真模型可以進行性能分析和設(shè)計參數(shù)靈敏度分析。

2多學科集成優(yōu)化模型

2.1　設(shè)計變量

通過對協(xié)同仿真模型分析，進行設(shè)計參數(shù)靈敏度分析，確定了包括彈道參數(shù)、彈丸結(jié)構(gòu)參數(shù)、身管參數(shù)和伺服控制環(huán)的控制參數(shù)在內(nèi)的15個設(shè)計變量，具體如下：

1) 彈道參數(shù)：Lg,W0,Omg1,Omg2 ;

2) 彈丸參數(shù)：M,wh1,wh2,wh3,r2,ZYL;

3) 身管質(zhì)量:m;

4) 控制環(huán)參數(shù):kpn,Tin,kw,Tw,kd。

2.2　約束條件

約束條件包括設(shè)計變量和狀態(tài)變量約束，以下給出一些狀態(tài)變量約束：

Pmax≤340MPa(最大膛壓)

Pg≤43MPa(炮口膛壓)

Yt≥0.3(火藥能量利用率)

Lsc≥11000m(射程)

0.45(kg×dm-3)≤bz≤0.6(kg×dm-3) (火藥裝填密度)

0.25≤rsw≤0.7(火藥相對燃燒結(jié)束位置)

S≤0.91m(后坐長度)

dthmax≤0.09(伺服系統(tǒng)最大超調(diào)率)

Thtp≤1.0s(伺服系統(tǒng)峰值時間)

Wjd≤0.0015°(伺服系統(tǒng)控制精度)

2.3　目標

a) 集成優(yōu)化目標

1) 以殺傷面積最大為優(yōu)化目標 Max(area)

殺傷榴彈對目標的所造成的損害主要是通過榴彈爆炸時所產(chǎn)生的破片來實現(xiàn)對目標的殺傷作用，所以殺傷面積是衡量榴彈作用大小的一個重要指標，在保證相關(guān)要求下，應盡量追求榴彈殺傷面積最大化。

2) 以最大反后坐阻力最小為優(yōu)化目標 Min(Frmax)

射擊時傳遞到火炮架體上的載荷是后坐阻力Fr,因而最大后坐阻力Frmax是是設(shè)計和校核火炮架體強度和剛度的依據(jù)。后坐阻力Fr大，必然導致火炮架體結(jié)構(gòu)和構(gòu)件尺寸增大，增加火炮的質(zhì)量，因而應使最大后坐阻力Frmax盡量小以此減小火炮架體的受力。

3) 以火炮伺服系統(tǒng)調(diào)整時間最短為優(yōu)化目標 Min(ts)

火炮伺服系統(tǒng)一直追求準確性跟快速性，系統(tǒng)的穩(wěn)定時間是控制系統(tǒng)的一個重要指標，是衡量一個系統(tǒng)反應速度的重要指標。在滿足系統(tǒng)精度的基礎(chǔ)上，系統(tǒng)達到穩(wěn)定的時間越短，表明系統(tǒng)對外界指令響應越快。

b) Parato優(yōu)化解集

Parato解集的提出是由于多目標優(yōu)化問題中各個目標間是相互沖突的，優(yōu)化解不可能是單一的解，而是一個解集，由此出現(xiàn)了parato最優(yōu)解集。采用parato機制直接處理多個目標，不需要將多個目標轉(zhuǎn)化為單一目標，因此克服了歸一化方法的缺點。

本文采用NSGA-Ⅱ優(yōu)化算法，得出Parato最優(yōu)解集，根據(jù)實際需要，選取一組最優(yōu)解。

2.4　優(yōu)化流程圖

火力系統(tǒng)多學科優(yōu)化設(shè)計總體流程圖如圖4所示。

圖4　火力系統(tǒng)多學科優(yōu)化設(shè)計流程圖

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

采用NSGA-Ⅱ優(yōu)化算法，經(jīng)過1401次迭代，求得parato解集如圖5所示。

圖5　Parato解集

同時各目標隨著迭代次數(shù)的變化趨勢如圖6~圖8所示。

圖6　殺傷面積收斂歷程

圖7　最大反后坐阻力收斂歷程

圖8　伺服系統(tǒng)穩(wěn)定時間收斂歷程

根據(jù)實際需求，在優(yōu)化解集選取一組合理的解，與原來相比，殺傷面積由原值1198.57m2增加為1650.339 m2，增大37.6%；最大反后坐阻力由原值165.393kN下降為133.709kN，減小19.1%；伺服系統(tǒng)穩(wěn)定時間由原值3.808s下降為2.191s，減小42.4%。

由優(yōu)化結(jié)果(表1)可知，通過尋優(yōu)求解，彈丸的殺傷面積，反后坐裝置的最大反后坐阻力以及伺服控制系統(tǒng)的穩(wěn)定時間都得到了較大的改善，整體性能得到較大的提高。

表1　優(yōu)化結(jié)果比較

續(xù)表1

4結(jié)論

1) 建立了包括全彈道、炮身/反后坐裝置、火炮位置伺服系統(tǒng)在內(nèi)的火力系統(tǒng)多學科協(xié)同仿真模型和集成優(yōu)化模型。該模型以殺傷面積最大、最大反后坐阻力最小和系統(tǒng)調(diào)整時間最短為優(yōu)化目標，采用多目標遺傳算法進行尋優(yōu)。在滿足約束條件下，求得優(yōu)化解，為自行火炮火力系統(tǒng)一體化設(shè)計提供了一種新方法。

2) 進行了某自行火炮火力系統(tǒng)多學科集成優(yōu)化設(shè)計，彈丸殺傷面積、最大反后坐力及伺服控制系統(tǒng)調(diào)整時間都得到了大的改善。約束變量均在要求范圍內(nèi)，同時得到parato解集，為決策者提供了多組可行的方案，可按實際需要供決策者選擇，為該火炮綜合性能的提高提供了優(yōu)化的技術(shù)參數(shù)。

參考文獻:

[1] 錢林方，侯保林，徐亞棟，等. 火炮彈道學[M]. 北京：北京理工大學出版社，2009.

[2] 高躍飛. 火炮反后坐裝置設(shè)計[M]. 北京：國防工業(yè)出版社，2010.

[3] 李克婧，張小兵. 火炮全彈道過程仿真及應用[J]. 火炮發(fā)射與控制學報,2010,2:48-52.

[4] 宗士增，錢林方，徐亞棟. 火炮反后坐裝置動力學耦合分析與優(yōu)化[J]. 兵工學報,2007，28(3):272-275.

[5] 洪亞軍，曹巖楓，尹強，等. 火炮身管-反后坐裝置集成優(yōu)化設(shè)計方法研究[J]. 中國機械工程學報,2013,24(1):15-20.

[6] 李銀伢，盛安東，王遠鋼. 基于粒子群優(yōu)化的伺服系統(tǒng)比例積分微分控制器設(shè)計方法[J]. 兵工學報,2006,27(2):202-205.

[7] 張相炎，鄭建國，楊軍榮. 火炮設(shè)計理論[M]. 北京：北京理工大學出版社，2005.

[8] van MICHEL T，LARS K. Multi-disciplinary design of aircraft fuselage structures[C]. 45th AIAA Aerospace Sciences Meeting and Exhibit. Reno:AIAA,2007：1-13.

[9] 洪亞軍，趙彥峻，王亞平，等. 小口徑榴彈集成優(yōu)化模型與方法研究[J]. 彈道學報,2010,22(2):77-80.

歡迎投稿歡迎刊登廣告歡迎訂閱

A Self-propelled Artillery Firepower System Integration Optimization Design Method

XIAO Xiao, XU Cheng, XU Ya-dong

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

Abstract:In order to improve the comprehensive performance of a self-propelled gun, this paper integrates the trajectory model, gun recoil mechanism model and system control model which directly affect the effect to the gun shooting and establishes the collaborative simulation model of the fire system. In this paper, by determining the main optimization design variables and with the largest lethal area, least reverse recoil resistance and shortest adjusting time of artillery servo system as the general optimization goal, the integrated optimization model of the fire system is established. Then,a multi-objective genetic algorithm is used for optimization to obtain parato optimal solution set of this system and according to the actual demand, a set of optimized solution is chosen. Optimization results show that the model of overall optimization goal is improved greatly. This article provides the effective design method for self-propelled guns fire system integration design and enhancing the cannon comprehensive performance.

Keywords:artillery; trajectory model; reverse recoil device; multiple target; multidisciplinary optimization

中圖分類號：TJ3;TP391.9

文獻標志碼：B

文章編號：1671-5276(2015)02-0146-04

作者簡介：肖曉(1990-)，女，碩士研究生，主要研究方向：復雜機械建模、仿真與優(yōu)化。

基金項目：國防基礎(chǔ)科研項目