3、飞行动力学建模：六自由度刚体方程、坐标系定义、运动学与动力学方程
倾转旋翼机的飞行动力学建模是整个控制律设计与仿真验证的数学基础。与常规固定翼飞机或直升机不同，倾转旋翼机在飞行过程中会经历直升机模式、过渡模式和固定翼模式三种截然不同的飞行阶段，其旋翼倾转角度从0°变化到90°，导致气动构型、惯性分布、推力矢量方向以及操纵效能发生剧烈变化。因此，建立准确、完整且适用于控制律设计的六自由度刚体动力学模型，是后续进行线性化、控制器设计和实时仿真验证的前提。本节将系统阐述建模所需的坐标系定义、运动学方程、动力学方程以及面向倾转旋翼机特殊性的建模要点。
上图展示了倾转旋翼机飞行动力学建模的整体框架，从坐标系定义出发，经过运动学与动力学方程，结合气动、旋翼、重力和操纵输入，最终形成完整的六自由度刚体方程。下面将逐一展开各部分的详细建模方法。
3.1 坐标系定义与转换关系
倾转旋翼机建模涉及多个坐标系，各坐标系之间的转换关系是建立运动方程的前提。建模中常用的坐标系包括：地面惯性坐标系（E系，原点取地面某点，X轴指向北，Z轴向下）、机体坐标系（B系，原点位于飞机重心，X轴指向机头，Y轴指向右翼，Z轴向下）、速度坐标系（W系，X轴沿空速方向）、以及针对倾转旋翼机特有的旋翼倾转坐标系（N系）和桨毂坐标系（H系）。
坐标系之间的转换通常采用欧拉角（滚转角φ、俯仰角θ、偏航角ψ）描述。从地面系到机体系的转换矩阵 LBE 由三个基本旋转矩阵按 Z-Y-X 顺序相乘得到：
L_BE = R_x(φ) · R_y(θ) · R_z(ψ)其中 R_x、R_y、R_z 分别为绕 X、Y、Z 轴的基本旋转矩阵。对于倾转旋翼机，旋翼倾转坐标系 N 系相对于 B 系绕 Y 轴旋转倾转角 δ_t，其转换矩阵为：
L_NB = R_y(δ_t)这一转换关系将旋翼产生的拉力和力矩从桨毂坐标系转换到机体坐标系，是倾转旋翼机建模中区别于常规飞行器的关键环节。
3.2 六自由度运动学方程
运动学方程描述飞机的位置和姿态随时间的变化关系，不涉及力的作用。位置方程将机体系下的速度矢量转换到地面系：
Ṗ_E = L_BE^T · V_B其中 V_B = [u, v, w]T 为机体坐标系下的速度分量。姿态方程描述欧拉角速率与机体角速度之间的关系：
Θ̇ = W(Θ) · ω_B展开形式为：
| φ̇ | | 1 sinφ·tanθ cosφ·tanθ | | p |
| θ̇ | = | 0 cosφ -sinφ | | q |
| ψ̇ | | 0 sinφ/cosθ cosφ/cosθ | | r |当 θ 接近 ±90° 时，矩阵 W(Θ) 出现奇异，此时应切换为四元数表示。四元数运动学方程为：
q̇ = 0.5 · Ω(ω_B) · q其中 Ω(ω_B) 是由机体角速度构成的反对称矩阵。四元数法避免了三角函数的计算，且不存在奇异问题，适合倾转旋翼机全包线飞行的仿真需求。
3.3 六自由度动力学方程
动力学方程基于牛顿-欧拉法建立，描述飞机在外力作用下的线运动和角运动。力方程在机体坐标系下表示为：
m · V̇_B + ω_B × (m · V_B) = F_aero + F_grav + F_prop其中 F_aero 为气动力，F_grav 为重力在机体系下的分量，F_prop 为旋翼产生的推力矢量。展开为标量形式：
m( u̇ + qw - rv ) = X_aero + X_grav + X_prop
m( v̇ + ru - pw ) = Y_aero + Y_grav + Y_prop
m( ẇ + pv - qu ) = Z_aero + Z_grav + Z_prop力矩方程在机体坐标系下表示为：
I · ω̇_B + ω_B × (I · ω_B) = M_aero + M_prop + M_gyro其中 I 为惯性张量矩阵，M_gyro 为旋翼陀螺力矩。惯性张量在倾转旋翼机中不是常数，它随旋翼倾转角 δ_t 变化：
I(δ_t) = I_fus + I_wing + I_tail + I_nacelle(δ_t) + I_rotor(δ_t)旋翼陀螺力矩是倾转旋翼机特有的重要力矩项，其表达式为：
M_gyro = Σ (I_rotor · Ω_rotor × ω_B)其中 I_rotor 为旋翼转动惯量，Ω_rotor 为旋翼转速矢量。当旋翼倾转时，该力矩对飞机的横航向动态特性有显著影响。
3.4 倾转旋翼机的特殊建模考虑
倾转旋翼机在建模中需要特别处理以下几个问题：
| 特殊项 | 影响机理 | 建模方法 | 典型量级 |
|---|---|---|---|
| 旋翼滑流干扰 | 旋翼下洗/滑流改变机翼和尾翼的有效动压和迎角 | 动量理论+涡流理论，建立滑流速度增量模型 | 滑流速度增量可达来流的30%-50% |
| 惯性张量时变 | 旋翼倾转导致质量分布变化 | 分段线性插值或多项式拟合 I(δ_t) | I_yy 变化可达15%-25% |
| 操纵冗余 | 旋翼总距、周期变距、倾转角、气动舵面同时作用 | 控制分配矩阵，加权伪逆法 | 4-6个操纵通道 |
| 气动非线性 | 大迎角、大侧滑角下气动系数非线性 | 风洞数据插值+CFD修正 | 失速迎角约15°-20° |
在旋翼滑流干扰建模中，需要计算机翼和尾翼处于旋翼滑流区内的比例，并修正当地动压和迎角。对于 V-22 这类倾转旋翼机，过渡飞行时滑流对机翼的干扰是影响纵向稳定性的关键因素。
3.5 案例拆解：过渡飞行段的建模与仿真
以某型倾转旋翼机从直升机模式（δ_t=0°）过渡到固定翼模式（δ_t=90°）的仿真为例，说明六自由度建模的实际应用。
背景与目标：该机型总重 15000kg，翼展 14m，旋翼直径 11.6m，两台发动机分别驱动两侧旋翼。要求在 30 秒内完成倾转过渡，验证过渡段姿态稳定性和高度保持能力。
关键输入参数：
初始状态: 悬停, 高度 500m, 空速 0m/s
倾转规律: δ_t(t) = 3°/s (0-30s)
旋翼转速: 589 rpm (直升机模式) → 412 rpm (固定翼模式)
发动机功率: 2×4580 kW建模与仿真过程：
- 建立六自由度模型，将旋翼拉力分解为沿机体轴的分量，并考虑倾转角变化引起的推力矢量变化。
- 引入旋翼陀螺力矩模型，计算倾转过程中旋翼角动量变化对机体的影响。
- 建立滑流干扰模型，计算机翼和尾翼在滑流区内的有效动压增量。
- 采用四元数法避免大姿态角下的奇异问题。
- 设置仿真步长 0.01s，采用四阶龙格-库塔积分算法。
仿真结果分析：
| 时间(s) | 倾转角(°) | 空速(m/s) | 俯仰角(°) | 高度变化(m) |
|---|---|---|---|---|
| 0 | 0 | 0 | 0 | 0 |
| 10 | 30 | 25 | 3.2 | -8 |
| 20 | 60 | 55 | 5.1 | -15 |
| 30 | 90 | 85 | 2.8 | -12 |
仿真结果显示，在倾转过程中飞机出现约 5° 的俯仰姿态变化和 15m 的高度损失。通过分析发现，高度损失主要来源于旋翼拉力垂直分量减小和机翼升力尚未完全建立之间的“升力缺口”。这一结果指导控制律设计时需要在过渡段增加总距补偿和俯仰姿态保持策略。
3.6 模型验证与确认方法
建立六自由度模型后，必须进行系统的验证与确认。验证方法包括：
- 静态验证：检查配平状态下力和力矩是否平衡，各状态量是否稳定。
- 动态验证：对比阶跃响应、频率响应与理论分析或风洞试验数据。
- 能量一致性检查：验证仿真过程中总能量（动能+势能）的变化与外力做功一致。
- 极限状态测试：在包线边界（大迎角、高速、大倾转角）下检查模型数值稳定性。
对于倾转旋翼机，特别需要验证过渡段的模型连续性。建议在 δ_t = 0°、30°、60°、90° 四个典型状态进行配平计算，检查状态量是否连续变化，避免出现数值跳变。
3.7 总结与建模流程复盘
六自由度飞行动力学建模是倾转旋翼机控制律设计的基石。本节核心要点可归纳为：
- 坐标系定义是建模的基础，倾转旋翼机需要额外定义旋翼倾转坐标系，并处理时变转换矩阵。
- 运动学方程描述几何关系，推荐使用四元数法避免奇异问题。
- 动力学方程基于牛顿-欧拉法，需特别处理时变惯性张量和旋翼陀螺力矩。
- 滑流干扰、操纵冗余、气动非线性是倾转旋翼机建模的三大难点。
- 模型验证必须贯穿建模全过程，不能等到建模完成后再进行。
完整的建模流程为：定义坐标系 → 建立运动学方程 → 建立动力学方程 → 建立力和力矩模型 → 集成六自由度方程 → 模型验证与修正 → 输出用于控制律设计的线性化模型。每一步都需要与物理概念和实际飞行数据相互印证，确保模型的准确性和可靠性。