2、过渡段动力学建模基础：坐标系定义与转换、刚体六自由度运动方程、气动与重力建模
各位同学，咱们今天聊聊过渡段动力学建模。说实话，这部分是整个倾转旋翼机控制设计的基石。你想想看，如果模型都搭错了，后面再花哨的控制算法也是白搭。我个人习惯，每次做新项目，第一件事就是把坐标系理清楚。
2.1 坐标系定义——别小看这一步
坐标系定义，说白了就是给飞行器找个「参照系」。我见过太多新手，上来就写方程，结果符号搞反了，仿真直接炸掉。嗯，这里要注意，我们至少需要三个坐标系：
- 地面坐标系（E系）：原点在地面某点，Z轴指向天。这个系用来描述飞行器的绝对位置和姿态。
- 机体坐标系（B系）：原点在飞行器质心，X轴指向机头，Y轴指向右翼，Z轴指向下。这个系用来描述飞行器自身的运动。
- 旋翼坐标系（R系）：原点在旋翼中心，随旋翼倾转。这个系是倾转过渡段特有的，用来描述旋翼拉力的方向变化。
核心要点：过渡段建模的难点，就在于旋翼坐标系相对于机体坐标系是时变的。倾转角从0°变到90°，所有力和力矩都在跟着变。
我曾经在一个项目中，因为搞混了旋翼坐标系和机体坐标系的旋转顺序，导致仿真结果和试飞数据差了30%。排查了整整两天，最后发现是欧拉角顺序写反了。所以，我建议大家在定义坐标系时，一定要把旋转矩阵写清楚。
2.2 坐标系转换——数学工具要趁手
坐标系转换，说白了就是「把同一个向量在不同坐标系下表示出来」。常用的转换矩阵有：
| 转换类型 | 转换矩阵 | 说明 |
|---|---|---|
| 地面→机体 | CEB = Rx(φ)Ry(θ)Rz(ψ) | 按Z-Y-X顺序旋转 |
| 机体→旋翼 | CBR = Ry(α) | α为倾转角，绕Y轴旋转 |
你可能会问，为什么旋翼转换只绕Y轴？因为在实际工程中，倾转机构通常只设计一个自由度。当然，也有双自由度倾转的构型，但那是进阶内容了。
个人经验：写代码时，我习惯把旋转矩阵封装成函数。比如 rotY(alpha)，这样可读性高，也不容易出错。千万别手算矩阵元素，容易算到怀疑人生。
2.3 刚体六自由度运动方程——飞行的物理本质
飞行器在空中的运动，说白了就是三个平动加三个转动。方程长这样：
平动方程：
m * (dV/dt + ω × V) = F_ext
转动方程：
I * (dω/dt) + ω × (I * ω) = M_ext这里V是速度向量，ω是角速度向量，F_ext和M_ext是合外力与合力矩。注意，这些方程是在机体坐标系下写的。为什么？因为惯性张量I在机体坐标系下是常数，省去了很多麻烦。
我记得刚入行时，总喜欢在地面坐标系下写方程，觉得直观。后来发现，每次都要对惯性张量做相似变换，计算量翻倍不说，还容易出错。所以，老老实实用机体坐标系吧。
2.4 气动建模——风洞数据是王道
气动建模，说白了就是描述「空气对飞行器的作用力」。对于倾转旋翼机，气动特性极其复杂：
- 机翼气动：升力、阻力、俯仰力矩。注意，过渡段机翼处于旋翼尾流中，气动数据会严重偏离常规值。
- 机身气动：通常简化为阻力项，因为机身产生的升力很小。
- 旋翼气动：拉力、反扭矩、挥舞力矩。这是最复杂的部分，因为旋翼既提供升力又提供推力。
避坑指南：我曾经在建模时，直接用了悬停状态下的旋翼气动数据来计算过渡段。结果仿真出来的爬升率完全不对。后来才发现，过渡段旋翼的入流比变化很大，必须用动态入流模型才能准确描述。
气动数据的获取，通常有两种方式：一是风洞实验，二是CFD计算。我个人更信赖风洞数据，因为CFD在复杂流场下误差较大。但风洞实验成本高，周期长，所以很多项目初期会用工程估算方法。
2.5 重力建模——简单但别忽略
重力建模，说白了就是「把重力向量从地面坐标系转换到机体坐标系」。公式很简单：
F_g = m * g * C_EB * [0, 0, 1]^T这里g是重力加速度，C_EB是地面到机体的转换矩阵。注意，重力方向始终指向地心，所以在机体坐标系下，它的分量会随着姿态变化而变化。
嗯，这里有个小细节：如果你在仿真中忽略了重力，飞行器会「飘」起来。我见过有人调试PID控制器时，发现积分项一直在增大，最后发现是忘了加重力补偿。所以，重力项虽然简单，但一定要加上。
2.6 知识体系总览
为了让大家更直观地理解本章的知识结构，我画了一张图：
这张图把本章的核心内容串起来了。你从「坐标系定义」出发，经过「运动方程」，最终落到「力与力矩建模」。三者缺一不可。
我的建议：初学者可以先从最简单的悬停状态开始建模，然后逐步加入倾转运动。别一上来就想搞定全包线模型，容易把自己绕晕。
好了，这一章的内容就到这里。记住，动力学建模是控制设计的地基，地基不稳，楼盖得再高也是危楼。下一章我们会讨论过渡段的线性化方法，到时候见。