折返的Liung

我如溪水,静静流淌在山涧~~我如白云,俯瞰此生今生~~

飞行器相关基础知识

1、平均气动弦长

平均空气动力弦长是—个假想矩形机翼的弦长,这一假想机翼的面积S和实际机翼的面积相等,它的力矩特性和实际机翼也相同。

Wiki上的解释:

一般的飞行器,在机翼各个位置的翼弦长度不尽相同,因此便有其它的定义。在分析飞行器的性能时,最常用的是平均空气动力弦长(平均气动弦长,mean aerodynamic chord,MAC,或是用$\bar{c}$表示),其定义如下:

$$\bar{c} = \frac{1}{S}\int_{-\frac{b}{2}}^{\frac{b}{2}} c^2(y) dy$$

其中$y$座标的方向是沿着机翼翼展方向,而$c(y)$则是在该位置上翼剖面的弦长。

另外一种纯几何的定义为标准平均弦长(standard mean chord,SMC):

$$SMC = \frac{S}{b}$$

其中$S$为翼面积,$b$为翼展长度。

2、缩减频率(reduced frequency)

定义:

$$k = \frac{\omega l}{V} = 2\pi f \frac{l}{V}$$

其中:$\omega = 2\pi f $,圆频率。$V$ 表示风速。$l$ 表示特征长度(characteristic length),对于纵向和横航向,分别有如下定义:

$$l = \frac{\bar{c}}{2} \quad or \quad \frac{b}{2}$$

其中,$\bar{c}$为平均气动弦长。$b$为模型展长。

3、时间的无量纲化(nondimensional time)

用 $t'$ 来表示无量纲化的时间,则:

$$t' = \frac{t V}{l}$$

其中,$V$表示风速。$l$ 表示特征长度。

4、关于动导数面积积分法的理解

用能量交换法来分析模型运动的稳定性:当飞机模型做俯仰振荡运动时。瞬时法向气动力的做功为:

$$ \Delta E = F_m * \Delta l $$

而$\Delta l$可以看做一小段弧长。所以有:

$$\Delta l = rad(t)*R$$

这里,$R$ 代表飞机绕轴心做俯仰运动的半径。$rad(t)$ 代表瞬时的弧度。经过分析可以知道,飞行过程中,$R$ 可以看做是特征长度$b$。$rad(t)$ 可以认为有如下关系:

$$rad(t) = \omega (t) * \Delta t = \dot{\alpha} (t) *\Delta t$$

另外,$F_m$ 可以如下表示:

$$ F_m(t) = \frac{1}{2} \rho V_\infty^2 S C_m (t)= \bar{q} S C_m (t) $$

整理可以得到:

$$ \Delta E = \bar{q} S b C_m (t) \dot{\alpha} (t) \Delta t $$

沿时间路径积分可以得出:

$$ E = \bar{q} S b \int_{t_1}^{t_2} C_m (t) \dot{\alpha} (t) dt = \bar{q} S b \int_{C_{\alpha}} \!\!\!\!C_m (\alpha) d\alpha \overset{\underset{\mathrm{一个周期}}{}}{=====} \bar{q} S b \oint C_m (\alpha) d\alpha $$

5、求解一阶常微分方程的四阶龙格库塔方法(RK4)

数值分析中,龙格-库塔法(Runge-Kutta)是用于模拟常微分方程的解的重要的一类隐式或显式迭代法。这些技术由数学家卡尔·龙格和马丁·威尔海姆·库塔于1900年左右发明。龙格库塔法的家族中的一个成员如此常用,以至于经常被称为“RK4”或者就是“龙格库塔法”。该方法主要是在已知方程导数和初值信息,利用计算机仿真时应用,省去求解微分方程的复杂过程。

关于初值问题表示如下:

$$ y' = f(t,y)\quad,\quad y(t_0) = y_0 $$

则,对于该问题的RK4由如下方程给出

$$ y_{n+1} = y_n + \frac{h}{6}(k_1 + 2k_2 + 2k_3 + k4) $$

其中:

\begin{align} k_1 &= f(t_n,y_n) \\ k_2 &= f(t_n + \frac{h}{2},y_n + \frac{h}{2}k_1) \\ k_3 &= f(t_n + \frac{h}{2},y_n + \frac{h}{2}k_2 \\ k_4 &= f(t_n + h,y_n + hk_3) \end{align}

这样,下一个值($y_n+1$)由现在的值($y_n$)加上时间间隔($h$)和一个估算的斜率的乘积决定。该斜率是以下斜率的加权平均:

  • k1是时间段开始时的斜率;
  • k2是时间段中点的斜率,通过欧拉法采用斜率k1来决定y在点tn + h/2的值;
  • k3也是中点的斜率,但是这次采用斜率k2决定y值;
  • k4是时间段终点的斜率,其y值用k3决定。
  • 当四个斜率取平均时,中点的斜率有更大的权值:

    $$\mbox{slope} = \frac{k_1 + 2k_2 + 2k_3 + k_4}{6}$$

    RK4法是四阶方法,也就是说每步的误差是h5,而总积累误差为h4阶。

    注意上述公式对于标量或者向量函数(y可以是向量)都适用。

    6.动导数

    1911年,Byran最早给出了动导数较为完整和严格的定义。飞行器受到的非定常气动力可以严格地表示为飞行状态变量的泛函。

    在线性的气动力系统中,气动力和力矩关于角位移的导数称为静导数,而气动力和力矩关于角速度的导数称为动导数

    动导数又称动稳定性导数(dynamic stability derivatives,DSD),用来描述飞行器进行机动飞行和受到扰动时的气动特性。在物理上一般表现为气动力对运动(状态变化)的阻尼作用。更一般地,动导数即飞行器所受气动力和力矩的各个分量,对飞行器各飞行状态参数的一阶时间变化率的偏导数。

    气动力和力矩分量:

    a.) 阻力D;侧向力Y;升力L

    b.) 滚转力矩$M_x$;俯仰力矩$M_y$;偏航力矩$M_z$;

    飞行器飞行状态变量:

    a.) 滚转角速度$\omega_x$;俯仰角速度$\omega_y$;偏航角速度$\omega_z$;

    b.) 迎角变化率$\dot{\alpha}$;偏航角变化率$\dot{\beta}$;飞行速度变化率$\frac{dV}{dt}$;

    共构成36个原始的动导数参数。

    实践表明,如下动导数是非常重要的:

    直接阻尼导数(direct derivatives):

    $$C_{lp} = \frac{\partial{M_x}}{\partial{p}},C_{nr} = \frac{\partial{M_y}}{\partial{r}},C_{mq} = \frac{\partial{M_z}}{\partial{q}}$$

    时差导数(plunging derivatives):

    $$C_{n \dot{\beta}} = \frac{\partial{M_y}}{\partial{\dot{\beta}}},C_{m \dot{\alpha}} = \frac{\partial{M_z}}{\partial{\dot{\alpha}}}$$

    交叉导数(cross derivatives):

    $$C_{lr} = \frac{\partial{M_x}}{\partial{\omega_y}},C_{np} = \frac{\partial{M_y}}{\partial{\omega_x}}$$

    交叉耦合导数(coupling cross derivatives):

    $$C_{lq} = \frac{\partial{M_x}}{\partial{\omega_z}},C_{nq} = \frac{\partial{M_y}}{\partial{\omega_z}},C_{mp} = \frac{\partial{M_z}}{\partial{\omega_x}},C_{mr} = \frac{\partial{M_z}}{\partial{\omega_y}}$$

    铰链阻尼力矩导数(hinge damping moment derivatives)(下标h代表铰链):

    $$C_{h \dot{\delta}_y} = \frac{\partial{M_h}}{\partial{\dot{\delta}_y}},C_{h \dot{\delta} _z}= \frac{\partial{M_h}}{\partial{\dot{\delta}_z}}$$

    Comments