飞行器六自由度仿真

1 引言

现在的战争已不是过去大刀长矛的时代，他早已成为国家综合实力的体现，这很大程度取决于军事高科技。这其中导弹作为精确打击的利器关乎国家的战略安全。而研究其包括导弹在内的飞行器精确制导与控制便显得十分的重要。

飞行器最优制导规律研究是进行武器系统总体方案论证和提高制导性能及精度的关键技术之一。而要进行制导规律最优性研究一方面需要研究合适的制导规律,另一方面需要进行接近实际情况的全面的大量的仿真研究。

仿真验证包括建模与仿真两个方面。在大型工程的方案论证阶段甚至包括实际研制的各个阶段，都要进行仿真检验以论证可行性、合理性和最优性。仿真技术在工程应用特别在高端武器系统总体设计和方案论证中具有极为重要的作用。对制导问题的研究在国外倍受重视。在公开发表的文献中,专门讨论制导规律方面的研究论文很多，可见制导规律的研究是非常重要的。但是仅有理论研究是远远不够的，因为设计的所谓最优制导规律大都是把实际系统进行了大量简化情况下推导出来的，因而与实际情况差别较大。也就是说理论上是最优的制导规律或参数在实际系统中不一定是最优的。因此,必须建立接近实际状态的数学模型和仿真软件。通过仿真计算确定出制导系统的最优参数及制导规律的控制效果，才能最终确定制导规律的最优性。

目前国内外这类问题研究主要存在下列三个问题:其一是模型被简化。从众多公开发表的文献资料看,大都是把控制系统简化为一阶、二阶、或三阶等根模型来推导制导规律,并据此模型进行仿真计算。其二是把飞行器的六自由度空间运动状态简化为平面运动状态进行仿真研究，以这种把飞行器运动限制在平面范围内进行仿真计算是有局限性的。其三是在全弹道数字仿真中仅选取几个特征点参数来代表全弹道的气动力参数，这些参数要表征全弹道动态过程是比较片面的，因而仿真结论的可信度是比较差的。

若把飞行器看成一个刚体，则它在空间的运动，可以看做是质心的移动和绕质心的转动的合成运动。质心的移动取决于作用在飞行器上的力，绕质心的转动则取决于作用在飞行器上相对于质心的力矩。在飞行中，作用在飞行器上的力主要有：总空气动力、发动机的推力和重力等。作用在飞行器上的力矩有：空气动力引起的空气动力矩，由发动机推力(若推力作用线不通过飞行器质心时)引起的推力矩等。作用在导弹上的空气动

力和空气动力矩在其他条件相同的情况下，取决于飞行器的气动外形。

飞行器的推力，是由发动机内的燃气流以高速喷出而产生的反作用力等组成。推力是飞行器飞行的动力。火箭发动机推力只与火箭的飞行高度有关，而与火箭的其他运动参数无关，它的大小主要取决于发动机的性能参数。

根据万有引力定律，所有物体之间都存在着相互作用力。飞行器在空间飞行就要受到地球、太阳、月球等等的引力。对于贴近地球表面的大气层内飞行的飞行器，所以只计地球对其引力。在考虑地球自转的情况下，除了受地心的引力外，还要受到因地球自转所产生的离心惯性力。

然而，高速飞行的飞行器一般是薄翼细长体弹性结构，因此有可能产生气动力和结构弹性的相互作用，造成弹体外形的弹性或塑性变形；操纵机构(如空气动力舵面)的不时偏转也相应改变飞行器的外形。同时运动着的飞行器也不是常质量的，对于装有火箭发动机的导弹，工作着的火箭发动机不断地以高速喷出燃料燃烧后的产物，使导弹的质量不断发生变化；对于装有空气喷气发动机的导弹来说，一方面使用空气作为氧化剂，空气源源不断地进入发动机内部，另一方面燃烧后的燃气与空气的混合气又连续地往外喷出。由此可见，每一瞬时工作着的反作用式发动机内部的组成不断地发生变化，即装有反作用式发动机的导弹是一个变组成系统。由于导弹的质量、外形都随时间变化，因此研究导弹运动不能直接应用经典动力学理论，而采用变质量力学来研究，这比研究刚体运动要繁杂得多。研究导弹的运动规律时，为使问题易于解决，可以把导弹质量与喷射出的燃气质量合在一起考虑，转换成为一个常质量系，即采用所谓“固化原理”，指在任意研究瞬时，设把变质量系的导弹视为虚拟刚体，把该瞬时在导弹所包围的“容积”内的质点“固化”在虚拟的刚体上作为它的组成。同时，把影响飞行器运动的一些次要因素通常也略去，如弹体结构变形对运动的影响等。这时，在这个虚拟的刚体上作用有如下诸力：对该物体的外力(如气动力、重力等)、反作用力(推力)、哥氏惯性力(液体发动机内流动的液体由于导弹的转动而产生的一种惯性力)、变分力(由火箭发动机内流体的非定态运动引起的)等。其中后两种力较小，也常被略去。

采用了“固化原理”，可把所研究瞬时的变质量系的导弹的动力学基本方程可写成常质量刚体的形式，这时，要把反作用力作为外力来看待，把每研究瞬时的质量取代原来的常质量。研究飞行器绕质心转动运动也可以用同样方式来处理。

为了验证这些问题，本文对六自由度飞行器进行了建模与仿真研究。

2 六自由度飞行器对象建模

2.1 飞行器质心动力学模型仿真

做如下假设：认为飞行器不仅是刚体，且质量是常数；假设地面为惯性参考系，即假设地坐标为惯性坐标；忽略地球曲率，把地球看成平面；假设重力加速度不随飞行高度而变化；假设机体X轴和Y轴处于飞行器对称面内，因此惯性积Jxz和Jxy等于0。根据设计要求，我们来建立仿真飞行器质心运动动力学和飞行器绕质心转动的动力学方程的数学模型。

2.1.1 飞行器作用力仿真模块

由于质心的运动动力方程是建立在弹道坐标系上的，根据

{m dV

dt

=p cosαcosβ−X−mg cosθ

mV dθ

dt

=P(sinαcosγv+cosαsinβsinγv)+Y cosγv−Zsinγv−mgcosθ

mV cosθdψv

dt

=P(sinαsinγv−cosαsinβcosγv)+Y sinγv+Z cosγv

(2.1)

将2.1方程三项作为一个向量来考虑，将右边作为力的分量F x，F y，F z，则为：

[dV dt

dθdt

dψv dt ]

=

[

F x

m

⁄

F y

mV

⁄

−F z mV cosθ

⁄]

=[

100

01V⁄0

00−1V cosθ

⁄

][

F x

F y

−F z

]/m

(2.2)

对上式2.2进行积分，可以算出来V，θ，φv，具体建立simulink仿真框图如下：