车载火炮射击过程动态响应仿真与参数匹配优化*

李 杰,贾长治,杜中华

(陆军工程大学石家庄校区,河北石家庄 05003)

近年来,车载火炮凭借其较好的战略战术机动能力、较大的火炮威力、较低的列装成本等优势获得各国青睐，并得到了迅速发展。但随着车载火炮的发展,车炮不匹配问题越来越突出,且显著影响了车载火炮的作战效能。车炮不匹配主要表现为射击稳定性较差,射弹散布较大。可见,对射击稳定性进行仿真与优化对车炮总体性能优化匹配具有重要意义。

对于火炮射击稳定性及参数匹配性问题,文献[1-4]分别从稳定性数学模型、刚柔耦合、土壤条件及轮胎与地面的接触情况等方面对某型现役火炮进行了射击稳定性仿真;文献[5-6]则基于层次分析法对坦克与车载火炮的系统匹配性问题进行了评估,但鲜见应用虚拟样机技术对“一种底盘,多种负载”的车炮匹配问题进行仿真研究。本文基于虚拟样机技术研究某A口径车载火炮底盘装载B火炮时的射击稳定性问题(B>A),应用ADAMS/Insight模块进行仿真试验并计算各结构参数的相对灵敏度,获取备选优化参数序列,进而引入序列二次规划法对备选优化参数进行优化,为车载火炮总体匹配性优化奠定基础。

1 虚拟样机建模

车载火炮在射击过程中,后坐部分与摇架通过平移副连接,在炮膛合力与后坐阻力的共同作用下沿炮膛轴线完成后坐与复进运动;起落部分与上架在耳轴处建立旋转副,并沿耳轴轴线上下振动;回转部分与车载底盘在座圈中心建立旋转副,并沿座圈轴线左右振动;悬挂系统液压驻锄及后支撑与底盘固连,随底盘运动;驻锄及后支撑与地面间建立平面接触,并将火炮系统的受力传递给地面并接受地面的约束反力。根据车载火炮在射击时的运动情况,本文建立了某车载火炮17刚体,22自由度的多体模型。车载火炮拓扑结构如图1所示。

图1 某车载火炮拓扑结构图

其中,H1为平移副,H2为旋转副,H3为固定副,H4为平面接触副。

火炮射击时,后坐部分受到的主动力为炮膛合力Ppt、后坐部分重力mhg、制退机力φ0(复进时提供液压阻力φ0f和复进节制器液压阻力φff)、复进机力Pf及相关摩擦力T等[7]。

(1)

(2)

(3)

(4)

(5)

由于火炮后坐部分的受力和运动情况较为复杂,模型通过Fortran语言编制用户自定义程序加载火炮射击时的炮膛合力与反后坐力。此外,模型中分别用线弹簧阻尼器模拟高低机、方向机中由于涡轮蜗杆间隙与齿轮侧隙等结构因素造成的约束反力。在火炮后坐的过程中,支撑机构会受到土壤的约束反力的作用,文中对土壤的结构和力学性能进行简化,以弹簧阻尼器模拟土壤的支持力与摩擦力。基于ADAMS建立的车载火炮虚拟样机模型如图2所示。

图2 某车载火炮虚拟样机模型

本文以高低角0°、方向角0°、全装药、常温环境为基本工况,对某现役A口径车载火炮虚拟样机进行射击仿真,并与相同条件下某基地试验结果进行对照。如图3～6分别为火炮后坐位移曲线、车尾跳高曲线、车尾侧移曲线及车尾侧倾角曲线。

图3 某车载火炮后坐位移曲线

图4 某车载炮车尾跳高曲线

图5 某车载火炮车尾侧移曲线

图6 某车载火炮车尾侧倾角曲线

依据某车载火炮试验数据对车载火炮虚拟样机模型进行可信性验证,试验与仿真对比如表1所示。

其中,各指标参数含义定义如下:

1) 最大后坐位移:炮身后坐时位移的最大值;

2) 车尾跳高:火炮射击时车尾中心参考点沿y轴的最大跳动;

表1 仿真结果与试验结果对比表

3) 车尾侧移:火炮射击时车尾中心参考点沿z轴的偏移量;

4) 车尾侧倾角:火炮射击时车尾中心参考点沿x轴的翻滚角。

2 仿真试验分析

2.1 仿真指标

根据某车载火炮射击过程中底盘系统的运动情况,选取车尾中心点为参考点,以车尾跳高、车尾侧移、车尾侧倾角为火炮射击静止性与射击稳定性衡量指标。

本文以高低射角0°、方向射角0°、全装药、常温装药为基本工况,加载力学模块链接库文件对虚拟样机进行仿真。图7～图10分别为某车载底盘分别装载A口径火炮与B口径火炮时的仿真指标对比曲线(B>A)。

图7 不同火力系统的后坐位移对比曲线

图8 不同火力系统的车尾跳高对比曲线

图9 不同火力系统的车尾侧移对比曲线

图10 不同火力系统的车尾侧倾角对比曲线

根据图8～图10仿真结果对比可知,当某现役A口径车载火炮底盘装载更大的B口径火炮时,车尾跳高与车尾侧移发生了明显变化,具有较大的优化必要性;而车尾侧倾角变化量较小,是否具有改进价值仍需结合参数灵敏度进行进一步分析。

2.2 相对灵敏度分析

参数灵敏度即系统目标函数对某一参数的偏导数,反映了输入参数对系统响应的影响情况。设系统目标函数为

F=f(b1,b2,b3,…,bn)

(6)

其中,f(·)为广义函数;bj(j=1,2,…,n)为车载火炮结构参数。

根据有限差分法,使某一设计变量发生一个非常小的摄动Δbj,而其他变量保持不变,此时系统响应的变化量为ΔFj。因此,车载火炮后坐稳定性在bj处的灵敏度可表示为

(7)

当Δbj非常小的条件下,上式可表示为

(8)

若各参数的量纲不同,去除量纲不同对参数灵敏度的影响,则相对参数灵敏度可表示为

(9)

本文通过分析车载火炮的拓扑结构,设置后坐部分质量、底盘质量、方向机刚度、方向机阻尼等17个结构参数作为待研究变量,并应用ADAMS/Insight模块进行仿真试验,各结构参数设置如表2所示。

根据上述变量设计,应用ADAMS/Insight模块进行仿真试验设计并输出Html结果文件,得各参数的相对灵敏度如表3所示。

表2 结构参数变量设计表

(其中,由于炮身、摇架、上架与底盘的转动惯量Iyy与Izz取值较为接近,且考虑虚拟样机成功仿真需满足Ixx+Izz>

Iyy的必要条件,文中对同一部件的转动惯量在参数化过程中视为等比例变化,即Iyy=k·Izz。)

表3 各结构参数相对灵敏度

根据表3可知,各结构参数对车尾侧倾角的影响均较小,可见侧倾角的稳定性较好,可改进程度较低。因此,文中选择车尾跳高和车尾侧移作为稳定性匹配优化的研究目标。为减小参数优化过程计算规模,同时保证取得较好的参数优化结果,文中根据参数灵敏度分析结果剔除了部分影响程度相对较小的因素,将参数高低机刚度(DV-1)、高低机阻尼(DV-2)…后支撑半径(DV-17)等11个变量作为待优化参数,获得备选参数优化序列{DV-1、DV-2、DV-3、DV-5、DV-7、DV-11、DV-12、DV-14、DV-15、DV-16、DV-17}。

3 稳定性优化

3.1 序列二次规划法

序列二次规划法(SQP)是拟牛顿算法在约束优化方法上的推广应用,是目前公认的求解约束非线性优化问题的最有效的方法之一。其基本思想是在某一迭代点处构造一个二次规划子问题来近似替代原约束优化问题,并求解获得约束优化的一个改进迭代点,通过不断重复“构造-求解-构造”这一过程,直至得到满足一定收敛条件的近似最优解[8]。

minf(X),

s.t.hi(X)=0,i=1,2,…,l,

gj(X)≥0,j=1,…,m

(10)

X∈Rn

其中，X∈Rn为决策变量,f(X)为目标函数,hi(X)和gj(X)分别为等式约束函数和不等式约束函数。

当式(10)存在迭代点Xk时,运用泰勒展开在Xk处构造原问题的近似优化模型,化简后的二次规划问题如式(11)。

(11)

由于式(11)仅为原非线性约束最优化问题的近似,故其解未必是原方程的可行点。因此,设S=X-Xk,将上述二次规划问题转换为式(12)关于变量S的优化问题。

(12)

求解此二次规划问题,并将其最优解S*作为原问题的下一搜索方向Sk,则有Xk+1=Xk+Sk,将近似解Xk+1作为新的迭代点重复上述过程,直至满足终止条件即可求得二次规划问题的最优解。

3.2 多目标优化及结果分析

设置收敛精度为0.001,最大迭代次数100次,运用ADAMS/Insight模块进行优化求解。目标函数随迭代次数变化和目标函数优化前后对比分别如图11、图12所示,各结构参数及目标函数优化结果如表4所示。

图11 目标函数随迭代次数变化图

图12 目标函数优化前后对比图

表4 各结构参数及目标函数优化结果

根据表4所示结果,经过17次迭代之后目标函数值逐渐稳定于72.356 5 mm,与原104.796 3 mm相比降低了30.955 1%,优化效果较为明显。其中,车尾跳高由104.023 4 mm降低至72.168 7 mm,较初始值降低30.622 6%;车尾侧移由-10.588 8 mm降低至-5.209 1 mm,较初始值降低了50.805 6%。由此可见,在相同条件下车尾侧移较车尾跳高的改进空间略大。各子目标优化前后对比分别如图10、图11所示。

图13 车尾跳高优化前后对比图

图14 车尾侧移优化前后对比图

4 结束语

本文基于虚拟样机技术研究了某现役A口径车载火炮底盘换装更大的B口径火炮时的射击稳定性并进行了参数优化,得出三条结论。

1)通过两种火炮射击稳定性结果对比及参数相对灵敏度分析可知,当A口径车载火炮底盘装载B(B>A)口径火炮时,射击稳定性明显变差。其中,车尾跳高与车尾侧移变化较大,具有较大的改进空间;车尾侧倾角仅发生微小变化,且各参数影响不显著,可改进空间较小。

2)文中应用加权组合法将多目标优化问题转换为单目标优化问题,并引入序列二次规划法对稳定性目标函数进行了优化求解。经参数优化,车尾跳高与车尾侧移分别降低了30.622 6%与50.805 6%,优化结果为车载火炮总体性能优化匹配提供了参考。

3)射击静止性、射击稳定性和射击密集度是车载火炮车炮匹配性评价中功能匹配性的集中体现,而文中未能充分考虑车炮匹配性对射击密集度的影响。由于弹丸膛内运动时间较短,故同时监测炮口扰动情况与射击稳定性等指标并进行综合优化存在一定困难。因此,如何解决射击稳定性、静止性与炮口扰动的耦合问题并实现综合优化将是下一步的工作重点。