羅中峰　徐　誠　徐亞棟

南京理工大學,南京,210094

輪式自行火炮底盤匹配與優(yōu)化

羅中峰徐誠徐亞棟

南京理工大學,南京,210094

為解決輪式自行火炮底盤匹配問題，以某8×8輪式自行火炮為例提供了一套匹配優(yōu)化方案。該方案采用多學科設計優(yōu)化的思想，首先依據(jù)匹配要求將輪式自行火炮底盤匹配模型分解，規(guī)定各個子模型相應的功能，并建立優(yōu)化模型;然后將單學科可行法優(yōu)化策略引入輔助設計變量，選用由第二代非劣排序遺傳算法和序列二次規(guī)劃法構(gòu)成的組合算法進行求解；最后通過計算實例驗證了方案的合理性。

輪式自行火炮;底盤;多學科優(yōu)化設計;機動性;發(fā)射動力學;平順性

0　引言

現(xiàn)代戰(zhàn)爭特別是高科技條件下的局部戰(zhàn)爭對輪式自行火炮的戰(zhàn)術(shù)性能提出了越來越高的要求，這使得輪式自行火炮底盤的匹配與優(yōu)化問題越來越突出。由于目前有關(guān)自行火炮的研究主要偏重于自行火炮某一性能[1-3]和火炮自身結(jié)構(gòu)[4],因此輪式自行火炮底盤匹配與優(yōu)化問題長期得不到解決。筆者基于多學科設計優(yōu)化分析復雜系統(tǒng)的思想[5],以某8×8輪式自行火炮為例進行了輪式自行火炮底盤匹配與優(yōu)化的研究工作。

1　學科分解與系統(tǒng)分析

1.1學科分解

本文以自行火炮作戰(zhàn)效能指標[6](圖1)為基礎(chǔ)進行輪式自行火炮底盤匹配與優(yōu)化的學科分解。與底盤相關(guān)的指標如下：機動性能指標、毀傷能力指標中的發(fā)射精度和攜彈量、防護能力指標、設備維護指標中的運輸適應能力。攜彈量和防護能力用底盤載重量來衡量，載重量越大，攜帶彈藥和防護裝備的數(shù)目就越多，但同時也會降低自行火炮的機動性能；發(fā)射精度可以通過基于底盤的發(fā)射動力學模型來考查;運輸適應能力通過基于底盤的行駛平順性模型來考查。根據(jù)上述各個性能指標之間的關(guān)系和考查方式的不同，可以將輪式自行火炮底盤匹配與優(yōu)化模型細分為以下3個基于底盤的子模型：機動性模型(考查機動性、載重量)、發(fā)射動力學模型、行駛平順性模型。

圖1　自行火炮作戰(zhàn)性能評估指標

1.2系統(tǒng)分析

圖2　自行火炮子系統(tǒng)耦合關(guān)系圖

自行火炮各子系統(tǒng)的耦合關(guān)系可用圖2表示。

1.2.1機動性模型

最大行駛速度(km/h)：

(1)

式中,ηt為行駛系統(tǒng)的總效率;PeN為動力裝置的最大輸出功率,kW;fv為地面變形阻力系數(shù);δb為質(zhì)量增加系數(shù);m為自行火炮的戰(zhàn)斗全質(zhì)量,t;g為重力加速度,m/s2。

單位質(zhì)量功率(kW/t)：

Pm=ηtPeN/m

(2)

作戰(zhàn)行程(km)：

(3)

式中,Q為油箱容積,L;Hu為燃油低燃值;ρ為燃油密度,kg/m3;ηs為動力裝置的熱效率。

縱向坡度角[7]：

(4)

式中,Me為動力裝置的最大輸出扭矩,N·m；ic為自行火炮行駛系統(tǒng)的最大傳動比;r為車輪半徑,m。

最小制動距離[8]:

(5)

式中,φ為地面附著系數(shù);CS為空氣阻力系數(shù);AV為自行火炮在行駛方向的正投影面積,m2;ρs為環(huán)境大氣密度,kg/m3;vr為行駛過程中自行火炮相對空氣的速度,m/s。

自行火炮載重量即戰(zhàn)斗全質(zhì)量減去火炮、懸掛和車輪的質(zhì)量(t)：

m1=m-m2-m4

(6)

式中,m2為火炮質(zhì)量(包括炮身、搖架、炮塔等),t;m4為懸掛和車輪質(zhì)量,t。

1.2.2發(fā)射動力學模型

(1)后坐阻力。反后坐裝置的引入,使得發(fā)射過程中火炮傳遞給自行火炮車體的力不再是炮膛合力,而是反后坐裝置傳遞的后坐阻力FR[9]:

(7)

式中,FR0為初始后坐阻力;FRg為后效時期后的坐阻力;tg為彈丸在膛內(nèi)運動時間;tr為從彈丸啟動到后效期結(jié)束的時間。

(2)發(fā)射動力學模型。自行火炮發(fā)射時,車體的振動是一個相當復雜的多自由度振動。為了便于優(yōu)化,本文將整個車體看成一個整體,并只考慮車體的上下振動、前后振動與俯仰振動。已知底盤的載重量為m1,繞車體質(zhì)心的轉(zhuǎn)動慣量為I1,火炮質(zhì)量為m2,繞車體質(zhì)心的轉(zhuǎn)動慣量為I2,其他變量如圖3所示。

圖3　某自行火炮發(fā)射動力學模型

取車體初始質(zhì)心位置C為坐標原點,以水平向左為X軸正向,鉛垂向上為Y軸正向建立圖3所示的坐標系,車體繞質(zhì)心的轉(zhuǎn)角為θ,逆時針方向為正。由拉格朗日方程,得到系統(tǒng)運動微分方程：

(8)

式中,FRy為后坐阻力在垂直方向的分力;FRx為后坐阻力在水平方向的分力;MR為后坐阻力作用于車體質(zhì)心的力矩;Li為第i排懸架與自行火炮的車體質(zhì)心的水平距離;x為車體質(zhì)心的水平坐標;xi為第i排懸掛的水平坐標,i=1，2，3，4。

為了完整評價車體跳動,以F的值作為車體跳動的評判(F的值越大跳動越大):

F=w1(maxy-miny)/y*+w2(maxx-

minx)/x*+w3(maxθ-minθ)/θ*

(9)

式中,w1、w2、w3為權(quán)重系數(shù);y*、x*、θ*為歸一化系數(shù)。

(3)后坐靜止性。自行火炮的后坐靜止性要求射擊時自行火炮車輪始終與地面接觸并沿水平方向不移動。若要達到上述要求，則必須滿足下式:

FRcosβ≤μ1mg

(10)

式中,β為自行火炮射擊高度角;μ1為路面的靜摩擦因數(shù)。

β=0時,水平后坐阻力最大,有FRcosβ=FR,保證此時自行火炮靜止應有

J1=FR-μ1mg≤0

(11)

(4)后坐穩(wěn)定性。自行火炮的后坐穩(wěn)定性要求火炮射擊時車輪始終與地面接觸。若要達到上述要求,則必須滿足下式:

W=MR-mgL4≥0

(12)

1.2.3行駛平順性模型

(1)路面激勵。本文采用濾波白噪聲作為路面輸入模型[10]，即

(13)

式中,q(t)為路面不平度位移輸入;Gq(n0)為路面不平度8級分類標準中規(guī)定的路面不平度系數(shù)幾何平均值;w0(t)為協(xié)方差為1 m2/s的單位白噪聲。

本文采用E級路面譜，測試距離為1000 m。

(2)平順性模型。輪式自行火炮行駛模型如圖4所示,其中,qi(t)為路面不平度位移輸入,a1、a2、a3為輪距,C為車體質(zhì)心。

圖4　某自行火炮行走動力學模型

由拉格朗日方程,得到系統(tǒng)運動微分方程：

(14)

(3)車輛平順性評價。車輛的平順性是評價車輛性能的重要指標,ISO2631/1給出了平順性評價的近似方法：用車體質(zhì)心垂直方向的加速度aw來評價車輛的平順性[11]。規(guī)定:aw≤0.315m/s2為舒適;0.3151.25m/s2為非常不適。

2　優(yōu)化模型

2.1優(yōu)化準則與目標函數(shù)

本文研究的目的是在火炮(包括炮身、搖架、炮塔等)的各種特性不變的情況下，選擇合適的底盤參數(shù)，使自行火炮的相關(guān)性能達到最優(yōu)。通過分析自行火炮的作戰(zhàn)效能指標和系統(tǒng)分解結(jié)果可以發(fā)現(xiàn),在同樣的條件下,底盤載重量m1是影響自行火炮機動性的最大因素,m1越小,機動性越強,作戰(zhàn)行程愈大,防護性能的加強、載彈量的提高都要求加大m1,并且m1也是影響自行火炮車體擾動、靜止性、穩(wěn)定性和行駛平順性的重要因素，所以選擇一個合適的m1相當重要。需要補充說明的是,研究過程中底盤的外形以及各懸掛的位置和輪胎都是不變的。綜合上面幾點,在保證機動性等要求的情況下,以底盤的載重量m1最大為目標函數(shù)。除此之外，還要考慮底盤參數(shù)對自行火炮發(fā)射速度和射擊精度的影響，因為這兩個指標是自行火炮射擊指標中與底盤參數(shù)直接相關(guān)的，通過控制底盤參數(shù)，可以控制發(fā)射時車體的擾動，降低其對發(fā)射速度和射擊精度的影響。針對這點,要求最大高低角發(fā)射時車體在垂直面內(nèi)的擾動最小。由上述分析得到如下目標函數(shù):輪式自行火炮底盤載重量最大即max(m1);高低角最大時發(fā)射后車體在垂直面內(nèi)的擾動最小即min(F)。

2.2設計變量

本模型有如下的設計變量:動力裝置的最大輸出功率PeN、動力裝置的最大輸出扭矩Me、自行火炮的戰(zhàn)斗全質(zhì)量m、自行火炮行駛系統(tǒng)的最大傳動比ic、自行火炮的油箱容積Q、第1排懸掛的阻尼c1和剛度k1、第2排懸掛的阻尼c2和剛度k2、第3排懸掛的阻尼c3和剛度k3、第4排懸掛的阻尼c4和剛度k4、自行火炮的車體質(zhì)心的水平坐標x、自行火炮載重量繞車體質(zhì)心的轉(zhuǎn)動慣量I1。

2.3設計約束

本模型有如下的設計約束:自行火炮最大行駛速度v≥vmin,最小單位質(zhì)量功率Pm≥Pmin,作戰(zhàn)行程s≥smin,縱向坡度角α1≥αmin,最小制動距離S1≤Smin,后坐時自行火炮的靜止性指標J1≤0,后坐時自行火炮的穩(wěn)定性指標W≤0,行駛時自行火炮車體質(zhì)心在垂直方向的加速度aw≤amin。

3　優(yōu)化策略及優(yōu)化算法

3.1優(yōu)化策略

優(yōu)化策略是在尋優(yōu)算法的基礎(chǔ)上提出的設計計算框架。選擇合適的策略能夠更好更快地找到最優(yōu)解。由上節(jié)的系統(tǒng)分析可知：各子系統(tǒng)間有較多的公共設計變量和耦合設計變量，即各子系統(tǒng)間存在較強的耦合。多級優(yōu)化策略和多學科可行法在處理強耦合系統(tǒng)時效率較低。同等情況下，多學科可行法比單學科可行法的計算復雜度高。考慮到上述幾點，本文采用單學科可行法優(yōu)化策略。單學科可行法是一種單級多學科優(yōu)化策略，它通過引入輔助設計變量，實現(xiàn)子模型之間的解耦。針對本文所述耦合模型,單學科可行法的具體處理過程如下：

對于耦合設計變量m,引入輔助設計變量m′,附加約束為

h1=m-m′=0

(15)

(16)

對于耦合設計變量v,引入輔助設計變量v′,附加約束為

h3=0.6v-v′=0

(17)

其中,v為自行火炮最大行駛速度,正常行駛速度取最大行駛速度的60%。

引入輔助設計變量之后,各個子系統(tǒng)就可以單獨建模。建模過程中，可以根據(jù)已引入輔助設計變量的情況，對子系統(tǒng)的輸入做適當?shù)恼{(diào)整。

3.2優(yōu)化算法

優(yōu)化算法應根據(jù)實際優(yōu)化問題的特點進行選擇。由于本文的優(yōu)化模型是一個多峰模型，為了獲得全局最優(yōu)解，本文選用了由第二代非劣排序遺傳算法和序列二次規(guī)劃法構(gòu)成的組合算法。在優(yōu)化過程中，首先應用第二代非劣排序遺傳算法在設計空間進行多目標全局搜索，求出Pareto解集，然后對Pareto解集中的各解進行如下計算：

f=1-F+(m1-m10)/m10=m1/m10-F

(18)

式中,m10為底盤載重量的初始值。

找出f值最大的X1(擾動下降的百分比和載重量上升的百分比之和最大的點)。接著應用序列二次規(guī)劃法以X1為起點進行局部精確尋優(yōu)。在精確尋優(yōu)過程中以f值最大為目標函數(shù)，設計變量的范圍以及約束同全局搜索時一樣。具體的尋優(yōu)過程如圖5所示。

圖5　優(yōu)化設計流程

4　計算實例

綜合圖6、圖7可以看出,優(yōu)化結(jié)果是收斂的。由表1可以看出，載重量增加了21.0%,擾動下降了18.5%,各項約束都滿足要求;同時附加約束的殘差也較小(h1=0.34,h2=0.53,h3=0.14)。綜合上述結(jié)果可知,優(yōu)化是成功的，優(yōu)化模型是合理的。雖然單位質(zhì)量功率Pm、最大行駛速度v、作戰(zhàn)行程s比優(yōu)化前減小，但還滿足設計約束,同時通過優(yōu)化也使先前不滿足約束的指標縱向坡度角α1和車體質(zhì)心垂直方向的加速度aw滿足了約束,滿足約束的指標后坐時自行火炮的靜止性指標J1和后坐時自行火炮的穩(wěn)定性指標W得到加強。總之,優(yōu)化提高了自行火炮的性能,使自行火炮總體資源得到了更加合理的利用。

圖6　載重量收斂歷程

圖7　車體跳動收斂歷程

參數(shù)優(yōu)化前優(yōu)化后約束優(yōu)化后比較Pm(kW/t)34.9932.1≥30滿足s(km)941.6866.8≥800滿足α1(°)26.433.6≥30滿足S1(m)3.83.36≤5.0滿足v(km/h)112.0103≥100滿足aw(m/s2)0.650.53≤0.63滿足J1(N)-9984-20000≤0滿足W(N·m)-13986-28450≤0滿足m1(t)12.815.5增加21.0%F1.00.815下降18.5%

5　結(jié)語

本文針對輪式自行火炮底盤匹配與優(yōu)化問題，提出了一種學科分解方案；以底盤載重量最大和最大射角時車體擾動最小為優(yōu)化目標，建立了一套輪式自行火炮底盤匹配優(yōu)化方案；應用單學科可行法優(yōu)化策略處理了該模型的耦合關(guān)系；選用由第二代非劣排序遺傳算法和序列二次規(guī)劃法序列二次規(guī)劃法構(gòu)成的組合算法進行求解;最后通過具體實例驗證了匹配優(yōu)化方案的合理性。

[1]向濤波,劉重發(fā),趙凌,等.某型自行火炮行駛平順性仿真研究[J].機械設計與制造,2011(10):220-222.

Xiang Taobo,Liu Chongfa,Zhao Ling,et al.Modeling and Simulation on Ride Comfort of a Self-propelled Gun[J].Machinery Design & Manufacture,2011(10)：220-222.

[2]吳大林,馬吉勝,杜中華,等.自行火炮扭力軸強化疲勞壽命計算及分析[J].計算機仿真,2012,29(10):30-33.

Wu Dalin,Ma Jisheng,Du Zhonghua,et al.Calculation and Analysis on Fatigue Life of Torsion Link for Self-propelled Gun on Enhancement Road[J].Computer Simulink,2012,29(10):30-33.

[3]李本朋,潘玉田,李永強.某自行火炮的彈藥架路面譜響應分析[J].機械管理開發(fā),2011(1):28-31.

Li Benpeng,Pan Yutian,Li Yongqiang.The Response Analysis of the Road about a Ammunition of the Self-propelled Artillery[J].Mechanical Management and Development,2011(1):28-31.

[4]宋焦，何永，趙威，等.某迫榴式火炮反后坐裝置的結(jié)構(gòu)優(yōu)化及分析[J].四川兵工報，2013，24(7)：22-24.

Song Jiao，He Yong，Zhao Wei，et al. Optimal Structural Design and Analysis on Recoil Mechanism of Howitzer Mortar[J]. Journal of Gun Launch & Control，2013，24 (7)：22-24.

[5]梅小寧，楊樹興.復雜系統(tǒng)的多學科設計優(yōu)化綜述[J].工程設計學報，2010，17(3)：173-180.

Mei Xiaoning，Yang Shuxing． Review of Multidisciplinary Design Optimization for Complex Systems[J]. Journal of Engineering Design，2010，17(3)：173-180.

[6]楊志強，易華輝，唐志躍.基于三角白化權(quán)函數(shù)的自行火炮作戰(zhàn)效能評估[J].四川兵工報，2011，32(10)：58-60.

Yang Zhiqiang，Yi Huahui，Tang Zhiyue. Operational Effect Evaluation of Self-propelled Guns Based on Triangle Whitening Function[J]. Journal of Gun Launch & Control，2011，32(10)：58-60.

[7]畢小平，王譜凱.坦克動力-傳動裝置性能匹配與優(yōu)化[M].北京：國防工業(yè)出版社，2004.

[8]王良曦，王紅巖.車輛動力學[M].北京：國防工業(yè)出版社,2008.

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

[10]吳志成，陳思忠，楊林，等.基于有理函數(shù)的路面不平度時域模型研究[J].北京理工大學報,2009,29(9):795-798.

Wu Zhicheng,Chen Sizhong,Yang Lin,et al.Model of Road Roughness in Time Domain Based on Rational Function[J].Transactions of Beijing Institute of Technology,2009,29(9):795-798.

[11]蔡樹新.自行火炮行駛試驗路面譜仿真研究[D].石家莊：軍械工程學院,2003.

(編輯張洋)

Matching and Optimization for Chassis System of Wheeled Self-propelled Guns

Luo ZhongfengXu ChengXu Yadong

Nanjing University of Science and Technology,Nanjing,210094

In order to deal with the optimum design for the chassis system(CS) of wheeled self-propelled guns,the way of MDO was introduced for analyzing the complicated system.Taking optimum design for a 8 by 8 wheeled self-propelled gun for instance,the optimization model was established,then was partitioned into several subsystems.In accordance with the requests for the chassis matching with the fire system the respective computation models were provided.Subsequently the coupling relationships among subsystems were analyzed in a manner of MDO.Finally,the instrumental variables were introduced to decouple subsystems from each other by the individual disciplinary feasible strategy.A combinatorial optimization algorithm of non-dominated sorting genetic algorithm Ⅱ and sequential quadratic programming was selected to optimize the resulting problem and the calculation example was verified the validity of this design pattern.

wheeled self-propelled gun;chassis;multidisciplinary design optimization(MDO);mobility;launch dynamics;riding comfort

2013-09-10

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

TJ303DOI：10.3969/j.issn.1004-132X.2015.02.004

羅中峰,男,1984年生。南京理工大學機械工程學院博士研究生。主要研究方向為多學科優(yōu)化。徐誠,男,1962年生。南京理工大學機械工程學院教授、博士研究生導師。徐亞棟,男,1978年生。南京理工大學機械工程學院副教授。