非线性动力学建模：倾转过渡段的六自由度动力学方程
倾转旋翼机在过渡段飞行中，气动构型随短舱倾转角连续变化，旋翼滑流与机翼/机身/尾翼的干扰耦合显著，导致动力学模型呈现强非线性、时变性和多体耦合特征。建立准确的六自由度（6-DOF）非线性动力学方程，是后续进行非线性控制设计、仿真验证和飞行品质分析的基础。本节将系统推导倾转过渡段的刚体动力学方程，明确气动力、力矩和重力项的表达形式，并讨论建模中的关键处理技巧。
3.1 坐标系定义与转换关系
倾转过渡段动力学建模首先需要明确坐标系定义。通常采用地面坐标系（惯性参考系，NED）、机体坐标系（原点位于重心，xb指向机头，yb指向右翼，zb指向下）和速度坐标系（用于定义迎角α与侧滑角β）。短舱倾转角δn是过渡段建模中的关键变量，它描述了旋翼拉力方向相对于机体坐标系的偏转角度。在δn=0°时对应固定翼模式，δn=90°时对应直升机模式，过渡段中δn在0°~90°之间连续变化。
坐标系之间的转换通过欧拉角（滚转φ、俯仰θ、偏航ψ）实现。地面坐标系到机体坐标系的转换矩阵Rbe采用Z-Y-X顺序旋转，其表达式为：
R_be = [cosθ·cosψ, cosθ·sinψ, −sinθ;
sinφ·sinθ·cosψ − cosφ·sinψ, sinφ·sinθ·sinψ + cosφ·cosψ, sinφ·cosθ;
cosφ·sinθ·cosψ + sinφ·sinψ, cosφ·sinθ·sinψ − sinφ·cosψ, cosφ·cosθ]在过渡段建模中，由于短舱倾转导致旋翼拉力方向变化，还需要引入短舱坐标系或旋翼轴坐标系。旋翼拉力在机体坐标系中的分量可表示为：
X_P = T·cos(δ_n)
Z_P = −T·sin(δ_n) （z轴向下为正，拉力向上为负）其中T为旋翼总拉力，δn为短舱倾转角。这种坐标变换关系是后续建立力与力矩方程的基础。
3.2 六自由度非线性动力学方程推导
基于牛顿-欧拉方法，倾转过渡段的六自由度刚体动力学方程可写为如下形式。力方程（平动动力学）：
m(u̇ + qw − rv) = X_A + X_G + X_P
m(v̇ + ru − pw) = Y_A + Y_G + Y_P
m(ẇ + pv − qu) = Z_A + Z_G + Z_P力矩方程（转动动力学）：
I_x·ṗ − (I_y − I_z)·qr = L_A + L_G + L_P
I_y·q̇ − (I_z − I_x)·rp = M_A + M_G + M_P
I_z·ṙ − (I_x − I_y)·pq = N_A + N_G + N_P其中，u、v、w为机体坐标系下的线速度分量，p、q、r为角速度分量，Ix、Iy、Iz为绕机体轴的转动惯量（假设惯性积为零或已通过坐标变换消除）。下标A表示气动力/力矩，G表示重力项，P表示旋翼/动力项。
值得注意的是，在过渡段中，由于短舱倾转导致质量分布变化，转动惯量Ix、Iy、Iz实际上是δn的函数。严格建模时应考虑转动惯量随短舱倾转角的变化率，即增加惯性时变项。但在工程实践中，若短舱倾转速度较慢（通常小于10°/s），可忽略惯性时变项，采用准定常假设。
3.3 气动力与力矩建模
过渡段的气动力建模是整节内容的难点。由于旋翼滑流对机翼、尾翼的干扰，气动导数不再是常数，而是迎角α、侧滑角β、短舱倾转角δn以及动压比的函数。气动力可分解为三部分：
- 机体气动力：升力Lwb、阻力Dwb、侧力Ywb，基于风洞实验或CFD数据插值获得。过渡段中机翼部分浸没在旋翼滑流中，需采用滑流/自由流分区积分方法。
- 尾翼气动力：平尾升力Lht、垂尾侧力Lvt，需考虑旋翼下洗对尾翼的干扰，下洗角ε是α和δn的函数。
- 旋翼滑流干扰项：旋翼滑流使机翼局部动压增大，产生附加升力和阻力增量ΔCL、ΔCD。常用动量理论或涡流理论估算滑流速度增量。
气动力矩包括俯仰力矩MA、滚转力矩LA和偏航力矩NA，其表达式为：
M_A = q̄·S·c̄·[C_m0(α,δ_n) + C_mα·α + C_mq·(q·c̄/2V) + C_mδe·δe]
L_A = q̄·S·b·[C_lβ·β + C_lp·(p·b/2V) + C_lδa·δa + C_lδr·δr]
N_A = q̄·S·b·[C_nβ·β + C_nr·(r·b/2V) + C_nδa·δa + C_nδr·δr]其中q̄为动压，S为参考面积，c̄为平均气动弦长，b为翼展。过渡段中，Cm0、Cmα等导数随δn变化显著，需建立多维插值表。
3.4 重力项与旋翼力/力矩项
重力项在机体坐标系中的表达式为：
X_G = −mg·sinθ
Y_G = mg·sinφ·cosθ
Z_G = mg·cosφ·cosθ重力项仅与欧拉角有关，与速度、角速度无关，在过渡段中形式不变。旋翼力/力矩项则较为复杂，包括：
- 旋翼拉力：T = C_T·ρ·A·(ΩR)²，其中C_T为拉力系数，是旋翼总距和前进比的函数。过渡段中旋翼来流速度随飞行速度变化，需使用叶素理论或动量-叶素组合理论计算。
- 旋翼反扭矩：Q = C_Q·ρ·A·(ΩR)²·R，反扭矩方向与旋翼旋转方向相反，在过渡段中需要通过尾桨或气动舵面配平。
- 陀螺力矩：旋翼高速旋转时，短舱倾转会产生陀螺力矩Mgyro = Irotor·Ω×δ̇n，该力矩在过渡段中不可忽略。
- 旋翼拉力在机体轴上的分量：如3.1节所述，XP = T·cos(δn)，ZP = −T·sin(δn)，同时拉力作用点与重心不重合时还会产生附加力矩。
3.5 建模中的关键处理与工程简化
在实际工程建模中，完全精确的非线性模型往往过于复杂，需要根据控制设计目的进行合理简化。下表总结了不同应用场景下的建模策略：
| 应用场景 | 模型复杂度 | 关键保留项 | 可忽略项 |
|---|---|---|---|
| 控制律初步设计 | 中等 | 重力项、旋翼拉力分量、基本气动导数 | 陀螺力矩、惯性时变项、高阶气动耦合 |
| 飞行品质评估 | 较高 | 动态导数、滑流干扰、短舱倾转速率效应 | 弹性变形、高阶非定常气动 |
| 飞行仿真器 | 高 | 全量非线性模型、气动插值表、旋翼动态 | 极小时间尺度效应 |
| 实时控制验证 | 低-中 | 线性化模型+增益调度 | 非线性耦合项 |
在过渡段建模中，一个重要的工程简化是采用准定常假设：认为气动力和力矩瞬时响应于当前飞行状态，忽略非定常气动效应。该假设在短舱倾转速率较低（<10°/s）时成立。另一个常用技巧是将旋翼滑流对机翼的干扰等效为动压比kw = q̄w/q̄，其中q̄w为机翼处有效动压，kw是δn和前进比μ的函数。
3.6 案例拆解：某倾转旋翼机过渡段建模实践
以某型倾转旋翼无人机为例，其质量m=1200kg，翼展b=10.5m，机翼面积S=8.2m²，平均气动弦长c̄=0.85m。短舱倾转角范围0°~95°，最大倾转速率8°/s。建模目标是为过渡段控制器设计提供仿真模型。
建模步骤与关键数据：
- 气动数据获取：通过CFD计算获得δn=0°、15°、30°、45°、60°、75°、90°七个状态下的气动导数表。每个状态计算α∈[−10°, 20°]、β∈[−15°, 15°]范围内的气动系数。关键导数Cm0随δn的变化呈现明显非线性：δn=0°时Cm0=0.02，δn=45°时Cm0=0.15，δn=90°时Cm0=−0.08。
- 旋翼模型：旋翼直径D=3.2m，转速Ω=680rpm（恒定），拉力系数CT通过叶素理论计算，考虑前进比μ的影响。在过渡段中，μ从0（悬停）变化到0.25（巡航），CT变化约15%。
- 滑流干扰修正：通过风洞实验获得动压比kw随δn和μ的变化规律。当δn=90°时，机翼完全处于滑流中，kw≈2.8；当δn=0°时，kw≈1.0。
- 模型集成与验证：将上述子模型集成到6-DOF方程中，在MATLAB/Simulink中搭建仿真环境。使用悬停配平点（δn=90°，V=0）和巡航配平点（δn=0°，V=60m/s）进行开环验证，仿真结果与飞行试验数据对比，俯仰角响应误差小于2°。
建模中发现的关键问题：在δn=30°~50°区间，俯仰力矩系数Cm0出现非单调变化，导致配平点不唯一。进一步分析发现，这是旋翼滑流在机翼上表面分离与再附着的非线性效应所致。该现象在控制律设计中必须考虑，否则可能导致过渡段中段出现俯仰振荡。
3.7 总结与建模复盘
本节系统建立了倾转过渡段的六自由度非线性动力学方程，涵盖坐标系定义、力/力矩方程推导、气动模型构建和旋翼/重力项处理。核心要点可归纳为：
- 坐标系与转换：明确地面系、机体系、速度系和短舱坐标系之间的转换关系，特别是旋翼拉力在机体轴上的投影随δn变化。
- 方程结构：6-DOF方程包含力方程（3个）和力矩方程（3个），各项来源清晰：气动、重力、旋翼/动力。
- 气动建模：过渡段气动导数随α、β、δn、μ变化，需建立多维插值表；滑流干扰是过渡段气动建模的核心难点。
- 旋翼建模：拉力、反扭矩、陀螺力矩均需考虑，且与δn及其速率相关。
- 工程简化：根据应用场景选择合适的模型复杂度，准定常假设和动压比修正是常用简化手段。