飞行动力学 - 第23节-运动方程的线化 之 基础点摘要

1. 运动方程的线化

飞机的全量运动方程耦合且非线性,很难获得用于理论设计和评估的解析解,但可以基于小扰动理论对其进行线化和解耦。

小扰动假设: x = x 0 + Δ x x = x_0 + \Delta x x=x0+Δx

飞机受扰动后,所有参数 = 基准状态量 + 小扰动量

  • 上世纪四十年代出现的小扰动理论已被证明能在实际使用中得到很好地结果,对于工程应用来说具有足够的精度
  • 不适合急剧机动等大扰动问题

2. 小扰动符号

在这里插入图片描述

3. 运动方程

在这里插入图片描述

4. 气动力模型

假设条件:

  1. 作用在飞机上的瞬时气动力/力矩仅与当时瞬时的运动参数有关
  2. 气动力/力矩随运动参数线性变化
  3. 纵向气动力与力矩(F_x, F_z, M)仅与纵向变量(u, a, q)相关
  4. 横航向气动力与力矩(F_y, L, N)仅与横航向变量(\beta, p, r)相关

这些假设仅适用于小 α / β \alpha / \beta α/β,当 α / β \alpha / \beta α/β较高时,气动力/力矩随 α / β \alpha / \beta α/β非线性变化。

5. 一阶状态方程

纵向与横航向方程均为一阶状态方程形式。

x ˙ = A x + B u \dot{x} = A x + B u x˙=Ax+Bu

在这里插入图片描述

6. 参考资料

飞行动力学-第23节-运动方程的线化

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
飞机纵向运动方程可以表示为: m * q_dot_dot = -m * g * sin(theta) - C_L * 0.5 * rho * V^2 * S - C_D * 0.5 * rho * V^2 * S * sign(V) + T 其中,m是飞机的质量,q_dot_dot是飞机的俯仰角加速度,g是重力加速度,theta是飞机的俯仰角,C_L是升力系数,rho是空气密度,V是飞机的速度,S是机翼面积,C_D是阻力系数,sign(V)是速度的符号函数(即当V>0时为1,当V<0时为-1),T是推力。 为了进行方程的线性化,我们需要对方程进行泰勒展开,并保留一阶项(即舍去高阶项)。假设当前状态为q1, q2, q3, q4, q5, q6,对应的状态变量分别为: q1: 飞机的俯仰角 q2: 飞机的俯仰角速度 q3: 飞机的速度 q4: 飞机的高度 q5: 发动机推力 q6: 大气密度 则可以将飞机纵向运动方程表示为: q1_dot = q2 q2_dot = (-m*g*cos(q1) - 0.5*rho*q3^2*S*C_L - 0.5*rho*q3^2*S*C_D*sign(q3))/m q3_dot = -q3*1.0/q4*(0.5*rho*q3^2*S*(C_L*cos(q1)-C_D*sin(q1))+q5-m*g*sin(q1)) q4_dot = -q3*sin(q1) q5_dot = 0 q6_dot = 0 其中,q1_dot表示飞机的俯仰角速度,q2_dot表示飞机的俯仰角加速度,q3_dot表示飞机的速度变化率,q4_dot表示飞机高度的变化率,q5_dot表示发动机推力的变化率,q6_dot表示大气密度的变化率。 将上述方程用matlab代码表示为: function [dqdt] = linearized_equations(t,q) % 参数定义 m = 10000; % 飞机质量 S = 50; % 机翼面积 g = 9.8; % 重力加速度 rho = 1.225; % 大气密度 C_L = 1.2; % 升力系数 C_D = 0.04; % 阻力系数 q5 = 50000; % 发动机推力 % 状态变量 q1 = q(1); % 飞机的俯仰角 q2 = q(2); % 飞机的俯仰角速度 q3 = q(3); % 飞机的速度 q4 = q(4); % 飞机的高度 q6 = q(6); % 大气密度 % 方程 dqdt = zeros(6,1); dqdt(1) = q2; dqdt(2) = (-m*g*cos(q1) - 0.5*rho*q3^2*S*C_L - 0.5*rho*q3^2*S*C_D*sign(q3))/m; dqdt(3) = -q3*1.0/q4*(0.5*rho*q3^2*S*(C_L*cos(q1)-C_D*sin(q1))+q5-m*g*sin(q1)); dqdt(4) = -q3*sin(q1); dqdt(5) = 0; dqdt(6) = 0; end 其中,t是时间变量,q为状态变量,dqdt为状态变量的导数。可以使用ode45函数对该方程进行求解。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值