飞行动力学与控制大作业

1、飞行动力学与控制大作业报告院（系）航空科学与工程学院专业名称飞行器设计学 号学生姓名目录一 .飞机本体动态特性计算分析11.1 飞机本体模型数据1.1.2 模态分析1.1.3 传递函数2.1.4 升降舵阶跃输入响应 21.5 频率特性分析4.1.6 短周期飞行品质分析 52. 改善飞行品质的控制器设计 62.1 SAS空制率设计6.2.1.1 控制器参数选择 .7.2.1.2 数值仿真验证 112.2 CAS空制率设计 123. 基于现代控制理论的飞行控制设计方法 1 63.1 特征结构配置问题描述1.63.1.1 特征结构的可配置性 1.63.1.2 系统模型1.63.2 系统的特征结构配置

2、设计 1.73.2.1 设计过程1.73.2.2 具体的设计数据173.2.3 结果与分析184. 附录20一.飞机本体动态特性计算分析1.1飞机本体模型数据本文选取F16飞机进行动态特性分析及控制器设计，飞机的纵向状态方程形 式如下：状态变量为：x ux = Ax + Buy = CxTq(1.1)控制变量为：u e基准状态选择为V 120m/s, H 2000m的定直平飞。选取状态向量x uq T ,控制量为升降舵偏角，则在此基准状态下线化全量方程所得到的矩阵数据如下：-0.0312 -1.1095 -9.8066 -0.5083-0.0013 -0.6543A000-0.382800.9

3、18501.00000-0.6901(1.2)B -0.0167 -0.0014 0 -0.0956 T(1.3)C diag 1.0000 57.2958 57.2958 57.2958(1.4)1.2模态分析矩P$ A的特征值算出为：1,2 -0.6778 + 0.5926i3,4-0.0100 + 0.0769i对应的特征向量如下：0.98740.9874-1.0000-1.00000.1137 - 0.0053i 0.1137 + 0.0053i 0.0011 - 0.0000i 0.0011 + 0.0000i V0.0521 - 0.0629i 0.0521 + 0.0629i 0

4、.002 1 + 0.0078i 0.0021 - 0.0078i0.0019 + 0.0735i 0.0019 - 0.0735i-0.0006 + 0.0001i -0.0006 - 0.0001i由系统特征值可知，系统具有两对共腕复根，也即具有两种运动模态：长周期模态与短周期模态，其对应的模态频率及阻尼比如下:表一飞机长短周期模态特征特征值自然频率阻尼比长周期3,4 -0.0100 + 0.0769i0.07760.1284短周期1,2-0.6778 + 0.5926i0.90030.7529可以看出，在此飞行状态下，飞机纵向具有明显的长周期模态， 但不具备明 显的短周期的模态特征，模态

5、频率过低，需要使用纵向增稳系统，改善阻尼比和 自然频率。1.3传递函数飞机迎角与俯仰角速度对应于升降舵输入下的传递函数如下：32(1.5)(1.6)0.08021s 5.088s 0.1615s 0.06983s4 1.376 s3 0.8436 s2 0.02432s 0.00488Gq s5.477s3 3.724s2 0.1034s 9.536e 016 s4 1.376s3 0.8436s2 0.02432s 0.004881.4升降舵阶跃输入响应由上述传递函数可得迎角与俯仰角速度在升降舵单位阶跃输入下的响应分别如下:迎角对升降舵输入的阶跃响应图1升降舵单位阶跃输入迎角时域响应上面阶跃

6、响应的性能指标为：稳态值为- 14.3090,调节时间为332.0859s超调量是37.6120%,上升时间是40.94004图2升降舵单位阶跃输入俯仰角速度时域响应上面阶跃响应的性能指标为：稳态值为0,调节时间为558.8424s超调量是2.9663/0,上升时间是 42.4104s。Bode Diagram 5010-210-1100101102103Frequency (rad/s)00803590451 1 1 -LUS3V esanp40Bode Diagram20o o o 2 4 - -90-3-210IO1100Frequency (rad/s)i1010210图4俯仰角速度对

7、升降舵响应传递函数Bode图迎角与俯仰角速度对应的传递函数的Bode图如下:B 0 n a -50图3迎角对升降舵响应传递函数Bode图1.6 短周期飞行品质分析飞机在当前状态下不具备短周期模态特征，短周期模态响应过大且频率过低, 操纵特性不符合飞行品质的要求，因此需要添加SAS 空制器来改善短周期模态阻 尼，提高短周期模态频率，使操纵品质满足要求。改善飞行品质的控制器设计2.1 SAS控制率设计增稳装置是在阻尼器的基础上发展而来的。阻尼器的作用主要是增加飞机的 俯仰阻尼，从而在一定程度上改善了飞机的短周期反应特性，但它不能改变飞机的纵向静稳定性，这时不能仅仅依靠阻尼器，必须借助于纵向增稳系统

8、（ SAS。纵向增稳装置除了俯仰角速度反馈回路之外， 还有对迎角或法向过载的反馈 回路，因而不仅能增加飞机俯仰阻尼，而且还能增加飞机的纵向静稳定性， 提高 飞机的短周期振荡频率，可以在更广阔的飞行范围内改善飞机的飞行品质。下图为纵向增稳装置的工作原理图。其中迎角变化是通过迎角传感器感受， 其信号输入经放大器放大后，再经舵机及助力器，推动舵面朝着减小迎角变化的 方向偏转。图5纵向增稳系统原理框图在本文中，对上述结构图做出如下简化：图6纵向增稳系统简化图2.1.1 控制器参数选择(2.1)在图6中，暂时忽略滤波器的作用，可得如下控制方程:eKqqK附加的气动导数增量为:Mq KqM e , M K

9、M q q e ,(2.2)通过俯仰角速度反从力学观点出发，通过迎角反馈，飞机的静稳定性增加, 馈，飞机的俯仰阻尼增加，从而可以改变飞机短周期的运动模态。为了使控制器设计更具实用价值，将迎角传感器与助力器的动力学特性分别20.2s 20.2表示为带宽10rad/s和20.2rad/s的一阶惯性环节,即:10G s , G ss 10首先考虑只有迎角反馈时的情况，即在原理图中kaWQ kq=0o取系统动态方程的状态变量为：x up e fT ,则系统的动态方程如下:Ax0 0 0 00 10 0 00 v0B 00 q20.20e010f20.2Ue(2.3)0 57.296 0000(2.4)

10、0 57.296 00x000 57.296迎角反馈回路的开环传递函数如下:(2.5)216.2033(s 63.39)(s2 0.03153s 0.01373)Z 72-"""" 2""-(s 20.2)(s 10)(s0.01993s 0.00602)(s 1.356s 0.8106)迎角反馈回路闭环后的根轨迹如下图所示:10内环迎角反馈回路的根轨迹sanoceocepYA aa m0.870.87-100.780.667.52.50.780.520.36*0.180.660.520.360.18-20-18-16-14-12-1

11、0-8-6-4-2Real Axis (seconds -1)图7迎角反馈回路的根轨迹因此将短周由于相较于助力器与迎角传感器，短周期模态对应的极点较小, 期极点局部放大，局部放大后的根轨迹图如下:内环迎角反馈回路的根轨迹0 12- - sanocect epxA 03 ml- 3- 40.70.820.910.9750.580.440.28 0.14System: sys_1 Gain: 1.62Pole: -0.102 + 2.82iDamping: 0.0361Overshoot (%): 89.3Frequency (rad/s): 2.820.9750.910.820.70.58-4-

12、3-2-10123Real Axis (seconds -1)图8迎角反馈回路的根轨迹局部放大图由图7可见，迎角传感器和助力器极点在左半平面离原点很远处，对其他模态影响较小。由图8可见，迎角反馈使短周期模态的频率逐渐增加但是阻尼比逐渐降低， 其中，频率的增加为主要变化。对短周期的影响基本可以忽略。因此，可以通过 选取适当的反馈增益使短周期的模态频率达到理想的值。选取k1.62,此时，短周期的频率为 2.82rad/s,阻尼比为0.036,因止匕，当前的短周期阻尼比偏低，需要引入俯仰角速度反馈增加短周期模态的阻尼比。当k 1.62时，以反馈迎角后的系统作为新的被控对象，此时俯仰角速度反 馈回路的

13、开环传递函数如下:Gqf S(2.6)110.6455 s (s+10) (s+0.6509) (s+0.02901) (s 19.89)(s 11.35)(s2 0.02944s 0.01307)(s2 0.3088s 8.078)从开环传递函数中可知迎角反馈对助力器与迎角传感器的极点几乎没有影 响，在俯仰角反馈时也有相同的结论， 此时，俯仰角速度反馈回路闭环后的根轨 迹如下所示:sdnocesT sx A 03(图9俯仰角速度反馈回路的根轨迹0.7440.58模态特征根自然频率阻尼比模态特性短周期-2.5431 ± 2.5427i3.59620.7072振荡收敛长周期-0.012

14、9 ± 0.0984i0.09920.1298振荡收敛通过上述分析可以看出，迎角反馈增益主要影响短周期模态的自然频率， 俯仰角速度反馈增益主要影响短周期的阻尼比。 选取适合的迎角反馈增益和俯仰角速度反馈增益进行组合，可以得到满意的自然频率和阻尼比， 进而在很大程度上改善飞机的纵向短周期模态特征图10俯仰角速度反馈回路的根轨迹局部放大图由上图可知，当选取kq 0.824时，短周期模态具有最佳阻尼比0.707,此时3205-1-2-3-5模态频率为3.59rad/s。可见，在k 1.62, kq 0.824时，短周期模态得到了很大的改善，具有良好的阻尼比和自然频率，满足良好的操纵性能要求

15、。采用 SAS控制器后，长短周期模态的特征值，自然频率及阻尼比如下:表二采用SAS控制器后飞行器长短周期模态特性1-4.5-4-3.5-3-2.5-2-1.5-1-0.50-1、Real Axis (seconds )俯仰角速率反馈回路的根轨迹0.30.140.9750.97540.720.440.31 &0.920.920.820.823System: sys_2Gain: 0.824Pole: -2.54 + 2.54iDamping: 0.707Overshoot (%): 4.31Frequency (rad/s): 3.592.1.2数值仿真验证当上述SAS控制器参数选为k1

16、.62, kq 0.824时，迎角与俯仰角速度对升降舵单位阶跃输入的响应如下所示:图11改善纵向稳定性后a的阶跃响应上面阶跃响应的性能指标为：稳态值为 0.5918,调节时间为196.7070s超调量是18.9655%,上升时间是30.8512so俯仰角速率对升降舵的单位阶跃响应50100150200250t(s) (seconds)300350400图12改善纵向稳定性后 q的阶跃响应上面阶跃响应的性能指标为：稳态值为0,调节时间为357.1331s超调量是 0.2/0,上升时间是31.1391s由图可以看出，对飞机本体进行SAS空制律设计后，短周期模态得到了很好 的改善，受到扰动后的震荡能

17、够迅速收敛且响应适中，飞行品质较好。2.2 CAS控制率设计无论阻尼器还是增稳系统，其目的都是改善飞机的模态特性（即稳定性），但经常还会导致静稳定性的下降。为解决稳定性和静操纵性之间的矛盾，对于以 机械式操纵为基础的飞机，在增稳系统的基础上增加前馈，即增加杆力（或杆位 移）传感器和指令模型，将驾驶员的操纵指令与飞机的响应构成闭环控制，形成所谓控制增稳系统，与机械操纵系统并联工作。因此可以采用CAS在SAS的基础 上改善飞机的激动性，以便满足操纵性要求以及多种响应类型的需要。CAS空制器的原理图如下:Kv伺服机构三卷 传卷罪一.工.闿叫图13CAS纵向控制原理框图若用公式表示如下:KxX KxK

18、ZyX KZyny KzZ z由上述原理图及公式可知，CAS控制器是在SAS空制器的基础上通过输入端 增加了一个PI控制器，即在反馈的同时增强控制作用从而提高飞机的机动性。下图为俯仰角速度的CASS制律原理图。图14俯仰角速度CAS控制律原理框图增加前向PI控制器，引入-2的零点，内回路选着k 0.02。迎角回路俯仰角速度反馈回路的根轨迹如下图所示：图15俯仰角速度反馈回路的根轨迹System: sys Gain: 0.53Pole: -1.95 + 1.95iDamping: 0.706Overshoot (%): 4.35Frequency (rad/s): 2.76-sanocawwt&

19、#171;pxA wya«poan.在上述根轨迹图中，选取kp控制图中选择虚线路线代替 者效果有所不同。俯仰角速度在0.53，此时短周期模态，con=2.76rad/s,己=706PI控制器的前向通道可以去除零点的影响，两CAS下加零点和不加零点对比图如图14所示:阶跃响应对比图1.41.210.80.6p m A0.40.2Time (seconds)图16俯仰角速度的阶跃响应由上图可知，在俯仰角速率单位阶跃输入的条件下，移除 PI控制器的零点 将会有效的减小系统的超调量，且与保留零点时系统的调节时间几乎相同，具有 更好的操纵特性。三.基于现代控制理论的飞行控制设计方法3.1 特征

20、结构配置问题描述从前面的分析中可以看出，系统的短周期模态并不明显，同时各个模态作用 量之间也存在着相互偶合的现象。在进行控制设计时不仅要考虑稳定性、操纵性 设计，同时也要对飞机各模态间进行解耦设计。 在现代控制理论中，特征结构配 置设计不仅可以配制出理想的特征模态同时可以完成对各个方向的解耦设计。3.1.1 特征结构的可配置性利用状态或输出反馈任意配置闭环零极点的充分必要条件是被控系统可控。设被控系统为 A,B,C,D ，当被控系统可控时，通过输出反馈和前馈校正可以进行系统的特征结构配置，其状态空间结构图如下:Uc增控部分系统配置状态空间结构图图173.1.2 系统模型系统的状态方程如下:(3

21、.1)x = Ax + Bu y = Cx状态变量为：x u控制变量为：U eT3.2 系统的特征结构配置设计3.2.1 设计过程下面简略介绍一下特征结构配置的设计过程。1）根据具体的增稳和解耦需求给出期望的飞机机动力学模型d d(3.2)x = A x + B uC y = Cx + Du C这里Uc （Vzc，Wxc，Wyc，Wzc）是角速度和垂向速度指令。2）用MATLAB等相关软件计算期望模型状态矩阵 Ad的特征值和特征向量 获得期望的特征结构。3）根据特征结构配置的实现过程计算出输出反馈矩阵K和前馈校正矩阵P （具体过程参考有关文献，附录给出了 MATLA瞅计程序）。4）根据图16所

22、示的系统状态空间结构图得出新的状态空间模型。(3.3)x = (A + BKC)x + BPu y = Cx5）对比配制出的特征结构是否与期望的特征结构基本一致，如果有一定出入首先检验期望模型是否合理，然后可对反馈矩阵K和前馈校正矩阵P 进行局部微调。3.2.2 具体的设计数据参考上文SAS设计结果，在保证飞机的飞行品质下给出期望模型的特征结构:表三期望模型的特征结构特征值-0.01 ± 0.1i-2.5 ±2.5i特征向量u100111q00长周期短周期设计过程中的输出反馈矩阵K和前馈校正矩阵P-1.7834 + 0.0000i 1.8584 - 0.0000i -0.0

23、139 - 0.0000i 0.6378 - 0.0000i K-0.0010 - 0.0000i 0.0023 - 0.0000i 0.0054 - 0.0000i 0.0007 - 0.0000i由于参考了上文SAS设计已经获得的设计结果，没有重新设计1）中的期望模型,所以P在此设计为单位矩阵。实际获得的特征结构：表四 实际获得的特征结构特征值-0.01 ± 0.1i-2.5 ±2.5i特征向量u0.99510.0004 ±0.0010i0.0121 ±0.0008i-0.1492 ±0.2209i0.9274-0.1855 ±

24、0.1855iq0.0097 ± 0.0968i0.0098 + 0.0000i长周期短周期3.2.3 结果与分析经特征结构配置后系统对升降舵的阶跃响应如下特征结构配置后迎角对升降舵输入的阶跃响应上面阶跃响应的性能指标为：稳态值为-0.0395,调节时间为535.5539s超调 量是715.4832%,上升时间是35.0219so图19配置后q的阶跃响应上面阶跃响应的性能指标为：稳态值为 0,调节时间为463.6558s超调量是0.5/0,上升时间是34.3331s由上面表三和表四的数据对比可以看出： 用特征结构的配置方法，特征值是 可以全部配置出的，也就是可以配制出期望的长短周期模

25、态， 但对应的特征向量 并不完全按照期望的配出，因为特征向量的配置还受到输入的影响，基本上是有 几个输入就可以精确地配置特征向量中的几个值。对比配置后得到阶跃响应图和 SAS得到的阶跃响应图可以看出迎角稳态值 有所减小，其他性能基本一致，从而印证了现代控制理论中的特征结构配置法在 改善飞机的操稳特性中是可行的。四.附录1)获取响应性能指标的MATLAB®序 function RPI(sys) y,t=step(sys);Y,k=max(y);% 峰值timetopeak=t(k-1)%到达峰值时间C=dcgain(sys)班值percentovershoot=100*(Y-C)/C%超调量i=length(t);while (y(i)>0.98*C)&(y(i)<1.02*C) i=i-1;endsettingtime=t(i)%调节时间if Y>Cn=1;while y(n)<