📚 课程目录第4章螺旋桨气动建模

4. 螺旋桨气动建模

螺旋桨是多旋翼唯一的推力来源，其气动模型直接决定后续参数辨识与控制器设计的可信度。工程上常用两条互补路径：动量理论从宏观控制体出发，把桨盘视为一维流管，得到推力与诱导速度的解析关系；叶素理论从桨叶剖面出发，把桨叶沿展向切成若干微元，逐段积分升力与阻力。两者结合形成叶素-动量理论（BEMT），是理解拉力系数、扭矩系数与来流影响的基石。

一、动量理论：从控制体到诱导速度

动量理论假设桨盘是一个无限薄的作动盘，气流不可压、无粘、轴对称。设桨盘面积 A=πR²，远前方来流速度 V_∞（悬停时为 0），桨盘处诱导速度 v_i，远后方滑流速度为 2v_i。由动量方程与伯努利方程可得理想推力：

T = 2 ρ A v_i (V∞ + v_i)
悬停时 V∞ = 0：  T = 2 ρ A v_i²
诱导速度：        v_i = sqrt( T / (2 ρ A) )
理想功率：        P_ideal = T (V∞ + v_i)
悬停诱导功率：    P_i = T · v_i = T^1.5 / sqrt(2 ρ A)

该式说明：同样拉力下，桨盘面积越大，诱导速度越小，诱导功率越低。这解释了多旋翼为何倾向大直径、低转速桨。但动量理论无法给出桨叶几何、翼型与扭转的影响，也无法解释扭矩与型阻，因此必须引入叶素理论。

二、叶素理论：沿展向积分

叶素理论把桨叶沿半径方向切成 N 个微元，每个微元视为二维翼型。设半径 r 处弦长 c(r)、安装角 θ(r)、扭转分布已知，来流由旋转速度 Ωr 与轴向速度（V_∞+v_i）合成，得到入流角 φ 与攻角 α：

入流角：  φ = atan( (V∞ + v_i) / (Ω r) )
攻角：    α = θ(r) − φ
微元升力：dL = 0.5 ρ W² c C_l(α) dr
微元阻力：dD = 0.5 ρ W² c C_d(α) dr
W = sqrt( (Ω r)² + (V∞ + v_i)² )
dT = dL cosφ − dD sinφ
dQ = r (dL sinφ + dD cosφ)

沿展向积分即得总拉力 T 与总扭矩 Q。BEMT 的关键在于：叶素给出的 dT 必须与动量理论给出的 dT 相等，从而迭代求解每个剖面的诱导速度 v_i(r)。这一步通常用数值迭代完成，是螺旋桨性能计算的核心。

三、拉力系数与扭矩系数

为便于跨桨对比与参数辨识，工程上把拉力与扭矩无量纲化。定义拉力系数 C_T、扭矩系数 C_Q、功率系数 C_P 与前进比 J：

C_T = T / (ρ n² D⁴)
C_Q = Q / (ρ n² D⁵)
C_P = P / (ρ n³ D⁵) = 2π C_Q
J   = V∞ / (n D)
效率 η = (C_T / C_P) · J   （悬停时 J=0，用品质因数 FM 评价）

其中 n 为转速（转/秒），D 为直径。悬停时 C_T 与 C_Q 近似为常数，因此拉力与转速平方成正比、扭矩与转速平方成正比，这是多旋翼混控分配的基础。下表给出典型 10 英寸桨的静态系数范围，用于辨识初值设定：

参数	物理含义	典型范围（10 英寸桨）	辨识用途
C_T	拉力系数	0.09 ~ 0.13	由 T–n² 拟合斜率反推
C_Q	扭矩系数	0.008 ~ 0.012	由 Q–n² 拟合斜率反推
C_P	功率系数	0.05 ~ 0.075	电功率与气动功率换算
FM	悬停品质因数	0.55 ~ 0.72	评估桨效率与型阻损失
Re	雷诺数	5×10⁴ ~ 2×10⁵	判断是否处于低 Re 修正区

四、雷诺数影响：小桨低转速的“隐形修正”

多旋翼桨叶弦长通常只有 1~3 cm，转速 3000~8000 rpm，剖面雷诺数 Re = ρWc/μ 常落在 5×10⁴~2×10⁵，处于低雷诺数敏感区。此时翼型升力线斜率下降、最大升力系数降低、阻力系数显著上升，甚至出现层流分离泡。若直接用大 Re 的风洞数据，C_T 会被高估 10%~20%。工程处理方式：

使用低 Re 翼型数据表（如 XFOIL 在 Re=5×10⁴~1×10⁵ 下的极曲线）；
对 C_l、C_d 做 Re 插值，BEMT 迭代中随剖面 W 更新 Re；
在参数辨识中把 C_T、C_Q 视为随转速缓变的“有效系数”，而非绝对常数。

提示：若辨识数据只在单一转速点采集，C_T 与 C_Q 的 Re 依赖性会被掩盖。建议至少在 3 个转速档位（如 40%、60%、80% 油门）分别拟合，观察系数漂移量。

五、静态与动态特性

静态特性指悬停或匀速爬升下的稳态关系：T = C_Tρn²D⁴，Q = C_Qρn²D⁵，输入为转速、输出为拉力与扭矩，可用台架实验直接标定。动态特性则关注转速变化、来流突变、桨叶挥舞与入流动态。多旋翼中最重要的动态环节是入流动态与电机-桨转速动态：

入流一阶近似：  τ_i · dv_i/dt + v_i = v_i,ss
τ_i ≈ R / (2 v_i,ss)   （悬停时约 0.02~0.08 s）
转速一阶模型：  τ_m · dn/dt + n = K_m · u
τ_m ≈ 0.05~0.15 s（含电机电气与桨惯量）

动态辨识时，常用扫频或阶跃油门激励，记录拉力/扭矩与转速响应，再拟合 τ_i、τ_m。若忽略入流动态，在阵风或快速机动仿真中会出现推力响应偏快、相位超前等误差。

常见误区：把静态台架测得的 C_T 直接用于前飞或大机动场景。前飞时来流使入流角变化，C_T 随前进比 J 下降，同时桨叶前后行侧不对称产生挥舞与振动。正确做法是建立 C_T(J)、C_Q(J) 曲线，或在辨识中引入来流修正项。

六、案例：10 英寸桨的 BEMT 辨识流程

台架采集：固定桨，测量转速 n、拉力 T、扭矩 Q，覆盖 2000~8000 rpm；
计算 Re 与 W，按剖面更新 C_l、C_d；
BEMT 迭代求解 v_i(r)，积分得理论 T、Q；
与实测对比，拟合 C_T、C_Q 的修正因子；
扫频激励，辨识 τ_i、τ_m，完成静态+动态模型闭环。

完成本节建模后，螺旋桨已可用“系数 + 动态时间常数”描述，为后续电机-桨联合辨识与控制器设计提供可计算的气动输入。