《自动飞行控制系统》课件04 05 06飞行器运动方程

单击此处编辑母版标题样式,,单击此处编辑母版文本样式,,第二级,,第三级,,第四级,,第五级,,,*,第二章 飞行器运动方程,飞行器运动方程组,,飞机的纵向运动,,飞机的横侧向运动,,2.1,飞行器运动方程组,建立飞行器运动方程时作出的假设条件（,Ma<3,）,:,,飞行器为刚体且质量是常数；,,地面坐标系为惯性坐标系,,,即假设地坐标为惯性坐标；,,忽略地球曲率,,,视地面为平面；,,重力加速度不随飞行高度而变化；,,假设机体坐标系的,OXZ,平面为飞行器的对称平面，飞行器不仅几何外形对称，而且内部质量分布也对称，即惯性积,I,XY,=I,ZY,=0,飞机运动的自由度,：对于飞机，若将其视为刚体，其在空间的运动需要六个自由度来描述质心的位移（,线运动,）：飞行速度的增减、升降和侧移运动；,,绕质心的转动（,角运动,）：俯仰角运动、偏航角运动以及滚转角运动纵向运动,（对称平面内运动）：速度的增减、质心的升降，绕,y,轴的俯仰角运动；,,横侧向运动,（非对称平面内运动）：质心的侧向移动、绕,z,轴的偏航角运动，饶,x,轴的滚转角运动2.1,飞行器运动方程组,飞机运动的特点,:,,飞机的基准运动为,等速直线平飞状态,，其小扰动线性化方程是常系数。

飞机的操纵面有,升降舵,、,副翼,和,方向舵,；,,飞机的外形通常是,左右对称,而上下不对称的面对称形体，垂直尾翼安装在机身后上部，便于地面起降这种布局致使机体水平转弯的效率很低，所以飞机一般,采用倾斜转弯,飞机的偏航和滚转运动间的,交叉影响显著,2.1.1,动力学方程,飞机动力学方程可由牛顿第二定律导出,,,该定律的向量形式为,:,,,,,利用前面①和②假设,,,上式可写为,:,,根据理论力学,,,速度向量对时间的变化率为,:,,,,,可用机体坐标轴系上的分量表示,:,,i, j,和,k,分别表示沿机体坐标轴系,OX,OY,OZ,的单位向量,.,V,上的单位向量,速度标量,叉积,飞机相对于地面坐标轴系总角速度向量,.,,2.1.1,动力学方程,由此可得,:,,,,,,展开上式可得,:,,F,也可用分量表示为,:,,利用前面一系列式子可得线运动,(,重心运动,),方程,:,,2.1.1,动力学方程,下面推导角运动,(,绕重心的运动,),方程,.,利用②假设,,,可写为,:,,,H,代表旋转的角动量或动量矩,.,,单元质量,dm,因角速度,所引起的动量等于单元质量绕瞬时转动中心的切线速度,Vq,乘以,dm.Vq,又可表示成,,因此,,,切线速度所引起的动量增量为,:,,动量矩等于动量乘以旋转臂长,,,写成向量形式为,:,,,对飞机的全部质量进行积分,,,可得总的动量矩,:,,,式中,:,,表示瞬时转动中心到单元质量,dm,的距离向量,.,,2.1.1,动力学方程,对飞机的全部质量进行积分,,,可得总的动量矩,:,,,式中,:,所以带如上式后得,:,,,,,定义,:,为惯性矩,Ix:,绕,X,轴的转动惯量,;,,,为惯性积；其他积分定义依此类推,.,,依据第⑤假设,,,I,xy,=,I,zy,=0,,将上式的分量写为,:,,,,,可写成,:,,,2.1.1,动力学方程,依据第⑤假设,,,可写成,:,,,的分量是,:,,,,,,上式推导中假设飞机是质量刚体,,,内部质量不在机内移动，则惯性矩和惯性积对时间的变化率为零,.,而,:,,,,展开后得,:,,,2.1.1,动力学方程,外力矩,,L,的分量形式为,:,,利用前面的一系列式子可得角运动方程,:,,,,,,两个加框的方程组是描述飞机非定常运动的两组动力学方程组,.,,以上推导是研究了动坐标轴系,(,机体坐标轴系,),相对于静坐标轴系,(,地面坐标轴系,),的动力学问题,.,而各力和力矩项都于飞机的空间方位,(,,,,,,),和位置有关,,,上述两组方程显然是不够的,.,,在空间运动的飞机有,6,个自由度,,,每,1,个自由度用一个二阶微分方程描述,,,整个飞机的方程就有,12,阶,.,但是上述两个加框的方程组总起来只有,6,阶,,,另外,6,个一阶微分方程可由飞机的运动学方程来补充,.,,运动学方程描述飞机相对于地面坐标轴系的空间方位,.,,,2.1.1,动力学方程,(,本教材,),在惯性坐标系中应用牛顿第二定律,:,,飞机在外合力作用下的线运动方程为,:,,飞机在外合力矩作用下的角运动方程为,:,,选用机体坐标系作为动坐标系,,,将在地面坐标系中得到的运动速度,V,及动量矩,L,向机体坐标轴系上分解,,,假设机体坐标系相对于惯性坐标系的速度为,V,,角速度向量为,,,,则上式在动坐标系中表示为,:,,,,如果将总空气动力 和发动机推力,T,向动坐标系,(,机体坐标轴系,),内分解为,(FX,FY,FZ),,则上式可写成下列方程组,:,,,2.1.1,动力学方程,角运动方程,:,,,,,,整理上式可以得到下列方程组,:,,,,,,,其中,:,,,2.1.2,运动学方程,飞机相对于地面坐标系的位置,,,可由机体坐标轴系相对于地面坐标轴系的三个坐标以及这两个坐标轴系之间的三个夹角,(,俯仰角,,,滚转角,,,偏航角,),来确定,.,运动学方程建立了,V,x,,V,y,,V,z,,p,q,r,与,X,g,,Y,g,,Z,g,,,,,,之间的关系,.,,,机体坐标轴系与地面坐标轴系之间的关系,:,,根据机体坐标轴系,OXYZ,和地面坐标轴系,OX,g,Y,g,Z,g,之间的几何关系,,,可得方向余弦表,(1),机体坐标,,地面坐标,OX,OY,OZ,OX,g,cos,cos,,,cos,sin,sin,-sincos,,sin,,sin,+cossincos,,OY,g,sincos,,cos,,sin,+sinsinsin,,sinsincos-,cos,,sin,,OZ,g,-sin,cossin,,coscos,,上表说明机体坐标轴系,OX,上的单位向量在地面坐标轴系三个轴上的分量各自为,:,cos,cos,,,sincos,,和,-sin.,,2.1.2,运动学方程,其次还需建立三个姿态角变化率 与三个角速度分量,(,p,q,r,),间的几何关系,.,,三个姿态角变化率的方位如下,:,,:,沿水平面内与,OX,轴在水平面上的投影线相垂直,,,向右为正,.,,:,沿,OX,轴的向量,,,向前为正,.,,:,沿,OZ,轴的向量,,,向下为正,.,,为了得到三个姿态角变化率与绕机体轴三个角速度间的转换关系,,,将三个姿态角变化率向机体轴上投影,,,得,:,,,,,2.1.2,运动学方程,应该指出,:,在一般情况下并不是互相垂直的正交向量,,,但,(,p,q,r,),却互相正交,.,故,:,,上式表示飞机三个姿态角变化率或绕机体轴的三个角速度分量都能合成飞机总角速度向量,.,一般情况下 与 与 互相垂直,,,但 与 不互相垂直,.,只有,,=0,时,,,与 才互相垂直,.,由上式可解出 的表达式（运动方程组）,:,,2.1.2,运动学方程,由机体坐标系中的速度分量可得到速度矢量在地面坐标系中的各分量（导航方程组）,,,2.1.2,运动学方程,速度坐标,,地面坐标,OX,a,OY,a,OZ,a,OX,g,cos,,cos,,cos,,sin,,sin,-sincos,,cos,,sincos+,sin,,sin,,OY,g,sin,,cos,,sin,,sin,sin+cos,,cos,,sin,,sincos-,cos,sin,,OZ,g,-sin,cossin,,coscos,,上表与,(1),表形式完全相同,,,若将,,,,,分别换成,,,,,,,表,(1),就成为表,(2).,速度坐标轴系与地面坐标轴系之间的关系,:,,由于气动力,,,气动力矩都与,,,,有关,,,必然涉及速度坐标轴系,.,根据速度坐标轴系,OXaYaZa,和地面坐标轴系,OX,g,Y,g,Z,g,之间的几何关系,,,可得方向余弦表,(2),,2.1.2,运动学方程,速度坐标轴系与机体坐标轴系之间的关系,:,,根据速度坐标轴系,OXaYaZa,和机体坐标轴系,OXYZ,之间的几何关系,,,可得方向余弦表,(3),速度坐标,,机体坐标,OX,a,OY,a,OZ,a,OX,cos,,cos,,-,cos,,sin,,-sin,,OY,sin,cos,,0,OZ,sin,,cos,,-,sin,,sin,,cos,,,2.1.2,运动学方程,前面各方向余弦表可看作转换矩阵,,,通过它们可以从一个坐标轴系转换到另一坐标轴系,.,利用表,(1),转换矩阵可将机体坐标轴系的变量转换到地面坐标轴系上,,,利用,(3),转换矩阵可将速度坐标轴系的变量转换到机体坐标轴系上,.,若将表,(1),转换矩阵用表,(3),转换矩阵右乘,,,得到速度坐标轴系和地面坐标轴系的转换矩阵,,,由此可得下列有用的几何关系式,:,,2.1.3,飞机动力方程的讨论和线性化处理,力方程组,力矩方程组,运动方程组,导航方程组,,2.1.3,飞机动力方程的讨论和线性化处理,前面得出的,2,组飞机动力学方程共,6,个联立的非线性运动方程式，再加上那些复杂的结合关系，以及气动力、气动力矩等都是运动参数的非线性函数。

因此要直接用这些方程解算飞机的运动，一般不能用解析法,,,而只能用数值积分法求解，即只能利用计算机但是无论如何飞机运动方程都是一组复杂的非线性微分方程,,,在研究飞机的稳定性和操纵性时,,,常根据小扰动原理对这组方程进行线性化处理,,,以便采用较简便的求解方法,.,,飞机的飞行运动分为,基准运动,和,扰动运动,而稳定性的关键是扰动运动能否回到基准运动2.1.3,飞机动力方程的讨论和线性化处理,一般小扰动线性化是相对原点或某点进行的,.,这里的小扰动线性化则是相对于基准运动进行的,.,,基准运动,(,又称,未扰动运动,),:,是指在完全理想的条件下,,,飞机按照驾驶员或飞行控制系统的意图按预定规律进行的运动,.,,扰动运动,:,是指飞机在外干扰作用下偏离基准运动,,,一段时间内违背预定规律的运动,.,,外干扰,可能来自于大气的紊流,,,发动机工作情况的改变以及驾驶员的偶然操纵等,.,它可以是瞬时的，也可以是持续性的,.,,若扰动运动与基准运动之间差别甚小,,,则称为,小扰动,运动,.,由于是小扰动，因此，可将那些含有扰动运动参数与基准运动参数件差值得高于一阶的小量即所谓高阶小项略去，方程变为线性方程。

首先研究最常见的等速直线平飞状态的稳定性问题基准运动就选择为没有倾斜、没有侧滑的等速直线平飞运动2.1.3,飞机动力方程的讨论和线性化处理,小扰动原理,:,,设运动方程组中的某一方程为,:f(x,1,,x,2,,…,x,n,)=0,式中变量,(x,1,,x,2,,…,x,n,),可以是运动参数或其导数,.,变量,x,i,可表示为基准运动时参数,x,0,与偏差量,,x,i,之和,,,即,:,x,i,=x,0,+,x,i,,无论是基准运动还是扰动运动都应满足运动方程,f(x,1,,x,2,,…,x,n,)=0,,即,: f(x,10,,x,20,,…x,n0,)=0,,f(x,10,+,x,1,,,x,20,+,x,2,,…x,n0,+,,x,n,)=0,,将扰动方程式的左边展成泰勒级数,,,在小扰动假设下,,,二阶和二阶以上的小量可略,,,则得,:,,,从上式中减去基准方程得,:,,,,这就是线性化的小扰动方程,.,式中系数 是已知的,.,,2.1.3,飞机动力方程的讨论和线性化处理,飞机运动方程的线性化处理,,选取,定常直线无侧滑,飞行为,基准运动,,,得基准运动参数有,:,,,根据小扰动原理,,,扰动运动参数可用基准运动参数附加一小扰动量来表示,,,即,:,,,,,,同样,,,可将基准运动和扰动运动的外力和外力矩表示为,:,,,,2.1.3,飞机动力方程的讨论和线性化处理,将上述扰动各参数表示式以及扰动运动外力和外力矩表示式带入飞机的运动方程组,(,前,2,个加框的方程组,),,减去其对应的基准运动方程,,,并略去二阶及以上的小扰动量,,,以小扰动量为变量的线性化方程,:,,,2.1.3,飞机动力方程的讨论和线性化处理,运动方程的分组,,前页方程组是常系数线性微分方程,,,假设飞机外形和内部质量分布对称于,X,s,OZ,s,平面而且有基准运动的左右对称性,,,那么方程组还可以简化,.,,由于存在这种对称性,,,我们将运动参数,(,扰动量,),分成对称的和不对称的两类：,,前进的速度,u,,俯仰角速度,q,等运动参数变化时,,,并没有破坏绕飞机气流的对称性,,,是对称的参数,,,因而这些参数的变化引起的气动力和力矩始终处于飞机对称平面,(,纵向平面,),内,.,,另一类运动参数,(,,,p,r,,,等,),是不对称的,,,引起不对称的气动力和力矩,.,,对称的参数,不会引起不对称的的气动力和力矩,,,而不对称的运动参数除了引起不对称的气动力和力矩外，还对纵向平面的力和力矩,(X,Y,M,等,),有一定影响,.,,2.1.3,飞机动力方程的讨论和线性化处理,运动方程的分组,,因此,,,在基准运动对称的前提下,,,纵向平面的力和力矩在基准点对不对称运动参数的一阶导数必为零，即：,,,对前式应用上述结论就可将方程组分成互不相关的两组方程组,.,,不论在等式的左边还是右边都只含对称平面内的运动参数,(,,,q,u,等,),,称为纵向扰动运动方程,;,几何关系式对基准运动线性化,,,得,,2.1.3,飞机动力方程的讨论和线性化处理,运动方程的分组,,只含不对称的运动参数,(,,p,r,,,等,),,称为横侧向扰动运动方程,.,几何关系式对基准运动线性化,,,得,:,,2.2,飞机的纵向运动,2.2.1,纵向运动的传递函数,,飞机纵向运动只涉及纵向的运动参数和气动力,,,又由于习惯用速度坐标系来表示空气动力,,,所以用速度坐标系建立纵向运动一般方程,,,以此推导纵向小扰动线性运动方程,.,飞机纵向受力图,:,,2.2,飞机的纵向运动,2.2.1,纵向运动的传递函数,,发动机推力,T,，方向沿发动机轴线，与机身轴线形成发动机安装角,,T,。

一般情况下发动机推力线不一定通过飞机重心重心对推力线的垂距为,ZT,，当重心在推力线之上,ZT,为正值时，推力,T,对重心之矩为正；,,升力,L,，垂直于飞行速度,V,，向上为正；,,阻力,D,，平行于飞行速度,V,，向后为正；,,俯仰力矩,M,a,（仅指气动力矩），抬头为正2.2,飞机的纵向运动,-,纵向运动方程,沿重心轨迹的切向方程,:,沿重心轨迹的法向方程,:,绕,OY,轴转动俯仰力矩方程,:,,2.2,飞机的纵向运动,-2.2.1,纵向运动方程,实际上,:,飞机着陆下滑时飞行速度较小,,,迎角较大,.,一般情况如巡航飞行,,,则速度较快迎角较小,.,发动机安装角,,T,在一般飞机上是一个很小的角度,,,因此可近似认为,:,cos(+,,T,)=1;,推力远小于重力,,(,或升力,) ,,即,:,Tsin(+,T,)<

2.2,飞机的纵向运动,,-,飞机纵向运动的特征（,1,）,采用系数冻结法将飞机纵向运动线性化方程组变为常系数线性微分方程组实际上各系数在基准运动中是变化的，若在飞机和飞行控制系统的过渡过程中，相对初始值的变化不超过,15-10%,，则可认为是常系数要了解纵向运动的特点，必须解出特征方程的根特征方程的根,,当外输入为零，纵向运动方程为：,,2.2,飞机的纵向运动,,-,飞机纵向运动的特征（,2,）,特征方程的根,,其特征方程为：,,,,,其中：,,,2.2,飞机的纵向运动,,-,飞机纵向运动的特征（,3,）,特征方程的根,,其特征方程为：,,,,,大量分析表明：升降舵在机翼之后的飞机的四个特征根中有两个大根、两个小根可近似写成两个二次因式之积：,,s,4,+a,1,s,3,+a,2,s,2,+a,3,s+a,4,=(s,2,+As+B)(s,2,+as+b),,A,B,a,b,均为待定系数可以用逐步近似法来确定其值2.2,飞机的纵向运动,,-,纵向扰动运动的典型示例、扰动运动的两种模态,(1),设：飞行高度,h=11000m,，以,M=0.9(V,0,=266m/s),作定常直线平飞,.,受扰动后,,,飞机偏离基准运动状态,.,现计算扰动因素消除后,,,飞机恢复到基准运动的过渡过程,.,,这完全靠飞机自身稳定性,,,驾驶员没有进行任何操纵,,,即,:,,,T,=0,,e,=0,,根据具体飞机结构参数和大导数表格,,,得到线性化方程的系数,,,并代入线性化方程得,:,,,,利用前述方法计算特征方程系数，可得：,,a,1,=1.4751 a,2,=8.9317 a,3,=0.1104 a,4,=0.01378,,2.2,飞机的纵向运动,,-,纵向扰动运动的典型示例、扰动运动的两种模态,(2),利用近似求解方法解出特征方程的根,,,由,,,s,2,+1.4752s+8.9317=0,,解得两个大根,(,近似解,):,,1,2,=-0.7375j2.896,,由,s,2,+0.01211s+0.001542=0,,解得两个小根,(,近似解,):,,3,4,=-0.00607j0.03895,,根据起始条件,t=0,时,,V==0,,0,=2,,解得近似解为,:,,,,由此可画出扰动运动的过渡过程曲线,,,见下页,.,,2.2,飞机的纵向运动,,-,纵向扰动运动的典型示例、扰动运动的两种模态,(3),t,s,,V,m/s,,,20,40,60,80,100,,,V,0,,-1,,-2,,2,,1,,0,-2,4,2,6,-4,-6,,V,, ,过渡过程曲线,可以看出,:,迎角,,在扰动运动的初期阶段变化剧烈,,,数秒后即平缓下来,.,速度,V,则缓慢增长,,,以后又缓慢减小,.,俯仰角,,兼有两者特点,,,开始阶段变化剧烈,,,以后又缓慢变化,.,,2.2,飞机的纵向运动,,-,纵向扰动运动的典型示例、扰动运动的两种模态,(4),两种扰动运动模态及其物理成因,:,,由前分析可以看出,:,扰动运动存在,两种模态,:,,短周期运动模态,(short period mode):,周期短,,,衰减快,.,其对应特征方程的一对大共轭复根,.,,长周期运动模态,(,phugoid,mode):,周期长,,,衰减慢,.,其对应特征方程的一对小共轭复根,.,,由此可以得出,结论,:,在外界瞬时扰动作用下,,,各运动参数随时间变化的规律正是这两种典型运动模态的迭加,.,,2.2,飞机的纵向运动,,-,纵向扰动运动的典型示例、扰动运动的两种模态,(5),两种扰动运动模态及其,物理成因,:,,飞机受到扰动后，出现不平衡的外力和外力矩，使飞机在受扰后的初始瞬间容易产生旋转运动,,,而其速度不易改变。

由运动方程知,,t=0,时有,:,,2.2,飞机的纵向运动,,-2.2.3,纵向扰动运动的典型示例、扰动运动的两种模态,(6),两种扰动运动模态及其,物理成因,(1),:,可见,:,在扰动运动的初始阶段,,,飞机的角加速度变化得比飞机速度剧烈得多,.,一般飞机都是如此,.,因为一般飞机纵向静稳定力矩,M,,较大,,,起始迎角,0,可引起较大的恢复力矩,,,相比之下,,,飞机的转动惯量,I,y,并不大,,,因而在扰动运动初瞬产生较大的角加速度值,;,反向的静稳定恢复力矩又使飞机向相反方向转动,.,于是形成迎角和俯仰角的短周期振荡,.,另一方面,,,飞机的阻尼力矩,M,q,较大,,,因而飞机短周期振荡运动的衰减较快,.,一般情况下,,,在扰动的前几秒就基本结束,.,飞机的力矩也基本上恢复到原有的平衡状态,.,,2.2,飞机的纵向运动,,-,纵向扰动运动的典型示例、扰动运动的两种模态,(7),两种扰动运动模态及其,物理成因,(2),:,起始迎角,,0,所产生的阻力和俯仰角变化所产生重力分量,,,远远小于飞机质量,,,因而初期线加速度,dV/dt,很小,.,但是在力矩基本恢复平衡之后,,,作用于飞机上的外力仍然处于不平衡状态,.,飞机的航迹仍未恢复到原有水平直线飞行状态,,,而是,q0.,因此,,,当升力大于重力沿航迹的法向分量时,,,产生向上的法向加速度使航迹上弯,,,飞机高度逐渐增加,.,与此同时重力沿航迹的切向的分力使飞行速度不断减小,,,升力也就不断减小,.,当升力小于重力的法向分量时,,,出现向下的法向加速度,,,航迹便转为向下弯曲,,,高度逐渐降低,.,这时重力的切向分量使飞行速度不断增大,,,又使升力在下降过程中不断增大,,,航迹再次上弯,.,如此反复,,,就形成飞行速度,V,和航迹倾斜角,q,的振荡运动,.,一般来说,,,飞机质量较大,,,而起恢复作用的气动力,Z,V,V,和阻尼作用的力,X,V,V,较小,,,因此振荡周期较长,,,衰减较慢,,,形成长周期运动模态,.,此时飞机的重心时升时降,,,又称为浮沉运动,.,,2.2,飞机的纵向运动,,-,纵向扰动运动的典型示例、扰动运动的两种模态,(8),综上所述,:,,飞机纵向扰动运动可大致分为两个阶段,’,,初始阶段,:,是以迎角和俯仰角速度的变化为代表的短周期运动,,,飞行速度基本不变,.,,以后的阶段,:,是以飞行速度和航迹倾斜角的变化为代表的长周期运动,,,飞机的迎角基本不变,.,,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,现在研究驾驶员操纵飞机时,,,也就是操纵升降舵,,e,和油门杆,,T,时,,,飞机的纵向响应,.,,纵向运动传递函数,-,升降舵偏转为输入,,纵向运动方程令,,T,=0,,并认为各变量初始条件为零,,,经拉氏变换式得,:,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,纵向运动传递函数,-,升降舵偏转为输入,(1),,以,,e,为输入,,,V,为输出的传递函数可写为两个行列式之比,:,,依据克莱姆法则（,Cramer rule,）,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,纵向运动传递函数,-,升降舵偏转为输入,(2),,以,,e,为输入,,,,为输出的传递函数,:,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,纵向运动传递函数,-,升降舵偏转为输入,(3),,以,,e,为输入,,,,为输出的传递函数,:,,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,纵向运动的初始阶段,,,短周期运动占主导地位,,,其过渡过程时间很短,,,飞行速度变化不大,,,可以认为速度增量,,V=0.,这样,,,纵向运动方程式第一式,(,切向力方程,),可以删去,,,其它两式当,V=0,时,,,得,:,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,短周期近似传递函数为,(1):,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,短周期近似传递函数为,(2):,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,常规飞机的升降舵在距重心较远的平尾上,,,平尾上的舵面小偏转引起的法向力足以产生较大的纵向控制力矩,.,因此,,,从工程近似的意义上来说,,,可以认为,Z,e,=0,,传递函数近一步简化为,:,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,航迹倾斜角增量与,,和有以下关系,: =-,，故,:,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,长周期运动的近似传递函数,:,,纵向长周期模态主要是飞机质心的轨迹运动,.,与短周期相比,,,长周期运动响应的各参数变化缓慢得多,.,因此长周期运动期间,,,短周期动态过程已基本结束,,,对短周期模态起重要作用的微分方程组的第三式基本上处于静力矩平衡状态,.,在简化处理时,,,将第三式中惯性力矩项和阻尼力矩项忽略,,,即忽略力矩从不平衡到平衡的动态过程,,,则飞机微分方程组简化为,:,,,,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,长周期运动的近似传递函数,:,,显然,,,上式已简化为二阶系统,,,切认为,Z,e,=0,,由拉氏变换后求得传递函数为,:,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,长周期运动的近似传递函数,:,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,由长周期运动简化微分方程组第三式得,:,将前式代入上式,,,可得长周期运动的第三个传递函数,:,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,纵向运动传递函数,-,油门杆偏转的动力学影响,(1),,操纵油门杆即改变发动机推力,.,推力变化对纵向运动方程第一式有直接影响,.,若考虑一般飞机的发动机推力线都通过重心,,,或非常接近于通过通过重心,,,则,M,e,0,,那么推力变化对第二,,,第三式无影响,.,因此改变推力时,,,长周期模态的影响将占绝对优势,.,令,:,e,=0,,且认为,: M,T,0,,,并令,:P,2,=0,,以及,M,q,=M,*,=0,,可得操纵,油门杆,T,的长周期近似运动方程组,:,,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,操纵,油门杆,T,的长周期近似运动方程组,:,,,,将上式经拉氏变换后可得,:,,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,纵向运动传递函数,-,油门杆偏转的动力学影响,(1),,油门杆阶跃偏转的运动参数稳态值,,由终值定理知,:,某参数的稳态输出为,:,,,油门杆阶跃偏转的拉氏变换为,:,,将上式代入前式得,:,,,所以,:,,同理可得,:,,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,前面三个稳定值得出结论,:,,油门杆前推发动机推力加大,,,但速度和迎角的增量最终值为零,,,也就是回到了推油门前的状况,;,,只有,为正,,,即飞机抬头,.,由几何关系知,:,轨迹角,=-,得,:()=()-0=().,即飞机向上爬升,.,,实际过程是,:,增加推力后先是增加速度,,,随后动压加大使升力增加,,,因而轨迹上弯,.,待航迹倾斜角到达一定正值后,,,重力沿轨迹的分力又使速度减小,.,在长周期动态过程结束后航迹倾斜角,,达到某一稳态值,,(),,使,: T=T,T,,T0,=,Gsin,G,,,即,:,增加的推力完全用于平衡重力沿轨迹的分力,,,而速度回原值,.,此外,,,由于没有偏转升降舵,,,迎角只能回原值,.,,如果推油门杆的目的是为增加速度而不是向上爬生,,,那么就应该配合速度的增加逐渐推驾驶杆,,,使升降舵下偏以减小迎角,,,使,L=G,,这样才能达到加速平飞的目的,.,,2.2,飞机的纵向运动,,-2.2.1,纵向运动的传递函数,纵向操纵方面的重要结论,:,,单纯改变油门杆只能在过渡过程中改变速度,,,最终的稳态速度和迎角均不改变,,,但飞机轨迹上升,(,或下滑,).,,如果加大推力是为爬升而不是为加快速度,,,那么加大油门时最好相应地拉驾驶杆,(,升降舵上偏以增大迎角,),来加快轨迹向上弯曲,,,待达到一定的上升航迹倾斜角后推驾驶杆,,,使升降舵回到原位,.,若不动驾驶杆,,,虽然最终飞机还是要到大爬升状态,,,但是过渡时间就太长了,.,,反操纵问题,:,,长周期模态稳定才能实现上述结论,.,若长周期模态不稳定,,,虽然从数学上能解出稳态值,,,但实际上是不可能实现的,,,因为是不稳定平衡,.,,许多飞机在进入跨声速段后都有不同程度的速度不稳定现象,,,故有必要在控制系统中用,M,数配平系统来解决这个问题,.,,2.3,飞机的横侧向运动,飞机的横侧向运动包括滚转、偏航和侧移三个自由度的运动。

操纵面是副翼和方向舵，它是飞机横侧向运动动力学环节的两个输入量滚转运动和偏航运动紧密地交联在一起,,,因而侧向运动问题是多变量,(,多输入，多输出,),系统的问题,.,,横侧向运动线性化方程,,因转动惯量是相对机体轴系求得的，故按机体轴系与地面坐标系来描述侧向运动方程,.,侧向力和力矩均与侧滑角,有关,,,应将它引入运动方程,.,,基准运动为等速直线平飞状态的横侧小扰动线性化方程如前面,P37,画框的方程组所示,,,现改写成,:,,2.3,飞机的横侧向运动,上式是对稳定轴系建立的，由于我们习惯用侧滑角,,,,而不习惯用侧移分速度,,,两者关系如下,:=/V,0,,(V,0,为基准运动的飞行速度,),,由于基准运动是对称直线平飞,,,因此有,:,,令,:,,及,,前式可改写为,:,,几何关系,:,,2.3,飞机的横侧向运动,改写成,:,,,,,,上式中各大导数的意义及计算公式列于表,.,,已知对机体轴系的各气动参数,,,转动惯量和惯性积等参数,,,可通过下列各式转换到稳定轴上,:,,2.3,飞机的横侧向运动,飞机横侧运动的特征,,侧向奇次方程的根,,令下式飞机侧向小扰动方程组右边为零,,,则奇次方程为,:,,2.3,飞机的横侧向运动,飞机横侧运动的特征,,侧向奇次方程的根,,令特征行列式为,:,,,,展开行列式,,,并令其等于零,,,得特征方程,:P,4,+b,1,p,3,+b,2,p,2,+b,3,p+b,4,=0 (b1b2b3b4,可从上展开行列式得出,.),,一般情况下正常式飞机侧向运动的,特征根,具有以下,特点,:,,具有一个绝对值较大的负实根,;,,具有一个绝对值很小的负实根,,,或是有一个很小的正实根,;,,具有一对共轭复根,,,其实部的绝对值介于两个实根的绝对值之间,.,,2.3,飞机的横侧向运动,飞机横侧运动的特征,,侧向运动实例,:,,侧向运动特征方程为,:,,P,4,+2.484P,3,+5.024P,2,+8.666P-0.0236=0,,解得特征方程,0,(s,）,=0,的根为,:,,模态,1 ,1,=-2.0768 (,大根,),,模态,2 ,2,=0.00272 (,小根,),,模态,3 ,3,4,=-0.205j2.0341 (,共轭复根,),,由此可见,:,,飞机侧向扰动通常由两个非周期运动模态和一个振荡模态组成,.,大的负实根对应于,滚转,运动模态,,,衰减很快,.,小实根,(,可正可负,),对应于,螺旋,运动模态,,,这种模态运动非常缓慢,.,如果具有小的正实根,,,则对应于缓慢发散的不稳定模态,.,共轭复根所对应的是,振荡,运动模态,.,,2.3,飞机的横侧向运动,横侧扰动运动的三种模态,:,,滚转模态,,荷兰滚模态,(,振荡模态,),,螺旋模态,,横侧运动的传递函数,(,略,),,横侧运动,(,三自由度,),传递函数,,二自由度传递函数,,一自由度传递函数,,,。