4. 空气动力学基础
分布式推进(DEP)是eVTOL区别于传统直升机的核心特征之一,其气动设计直接决定了整机的悬停效率、巡航升阻比、噪声水平与安全性。本章从螺旋桨/涵道的基本理论出发,系统梳理推力与功率的计算方法、桨尖速度的约束边界、载荷分布的优化逻辑以及多旋翼/多桨之间的气动干扰机理,为后续的电机选型、桨叶设计和总体布局提供量化依据。
上图展示了分布式推进气动设计的核心逻辑链:从一维动量理论到二维叶素理论,再到耦合的BEMT方法,最终落到桨尖速度、载荷分布与气动干扰等工程约束。下面逐一展开。
4.1 螺旋桨/涵道基本理论
4.1.1 动量理论:悬停与巡航的零阶模型
动量理论将螺旋桨视为一个无厚度的“致动盘”,气流通过桨盘时获得轴向动量增量。对于悬停状态(来流速度 V0=0),诱导速度 vi 与推力 T 的关系为:
vi = √(T / (2·ρ·A))
其中 ρ 为空气密度(海平面 1.225 kg/m³),A 为桨盘面积(A=πR²)。由此可得理想悬停功率:
P_ideal = T·vi = T^(3/2) / √(2·ρ·A)
该公式揭示了两个重要结论:推力一定时,桨盘面积越大,诱导功率越低;功率随推力的 3/2 次方增长,意味着大载重需求会急剧推高功率。对于巡航状态,动量理论给出:
P = T·(V0 + vi),其中 vi = (-V0 + √(V0² + 2T/(ρA))) / 2
当 V0 远大于 vi 时,功率趋近于 T·V0,即推进效率接近 100%。这正是固定翼巡航效率高于旋翼悬停的物理本质。
4.1.2 叶素理论:桨叶几何与气动力的映射
叶素理论将桨叶沿展向分成若干微段(叶素),每个叶素视为二维翼型,根据当地来流速度(合速度 W = √((ωr)² + (V0+vi)²))和攻角 α 计算升力与阻力:
dL = ½·ρ·W²·c·Cl(α)·dr
dD = ½·ρ·W²·c·Cd(α)·dr
将升力分解为轴向(推力)和切向(扭矩)分量,沿展向积分即可得到总推力 T 和扭矩 Q。叶素理论的关键输入是翼型升阻极曲线(Cl-α、Cd-α),以及桨叶的弦长分布 c(r) 和扭转分布 θ(r)。
4.1.3 叶素动量理论(BEMT):工程标准方法
BEMT将动量理论与叶素理论耦合,通过迭代求解轴向诱导速度 vi 和切向诱导速度 vt,同时引入 Prandtl 桨尖损失因子 F 修正桨尖区域的载荷衰减:
F = (2/π)·arccos(exp(-f)),其中 f = (N/(2r))·((R-r)/sinφ)
BEMT 的输出包括:沿展向的推力分布、功率分布、诱导效率、桨叶最佳弦长/扭转分布。这是目前eVTOL螺旋桨初步设计最常用的方法,计算量小且精度可接受(误差通常在 5% 以内)。
4.2 推力与功率计算
4.2.1 无量纲系数法
工程中常用无量纲系数进行推力与功率的快速估算:
| 系数 | 定义 | 典型值范围(eVTOL) | 用途 |
|---|---|---|---|
| 推力系数 CT | T / (ρ·n²·D⁴) | 0.05 ~ 0.15 | 推力估算 |
| 功率系数 CP | P / (ρ·n³·D⁵) | 0.04 ~ 0.12 | 功率估算 |
| 品质因子 FoM | P_ideal / P_actual | 0.65 ~ 0.80 | 悬停效率评价 |
| 桨盘载荷 DL | T / A | 30 ~ 120 N/m² | 悬停功率与噪声 |
其中 n 为转速(转/秒),D 为桨盘直径。品质因子 FoM 是衡量悬停效率的核心指标,理想值为 1.0,实际eVTOL螺旋桨通常在 0.7 左右。若 FoM 低于 0.6,说明桨叶设计或工作点选择存在明显问题。
4.2.2 功率分解与效率链
实际电机输入功率到螺旋桨气动功率的传递链为:
P_electric → P_shaft(机械损耗) → P_aero(气动功率) → P_induced(诱导功率) + P_profile(型阻功率) + P_parasite(废阻功率)
悬停时,诱导功率约占 70%~80%,型阻功率占 20%~30%。降低型阻功率的关键是减小桨叶面积(降低实度)和选用低阻翼型,但实度过低会导致桨叶载荷过高、失速提前。设计时需在两者间权衡。
P_hover ≈ (T^1.5) / (√(2·ρ·A) · FoM)
例如:T=5000 N,D=1.2 m(A=1.13 m²),FoM=0.72,则 P ≈ 5000^1.5 / (√(2×1.225×1.13) × 0.72) ≈ 353,553 / (1.664 × 0.72) ≈ 295 kW。
4.3 桨尖速度与马赫数约束
桨尖速度 Vtip = ω·R = π·n·D,是螺旋桨设计中最关键的转速约束。桨尖速度直接决定桨尖马赫数 M_tip = Vtip / a(a 为当地声速,海平面约 340 m/s)。当 M_tip 超过 0.7 时,桨尖区域出现显著压缩性效应:
- 激波诱导分离:桨尖翼型表面局部流速超过声速,产生激波,导致阻力急剧增加、升力下降。
- 噪声剧增:桨尖马赫数每增加 0.05,噪声级约增加 3~5 dB。
- 效率下降:桨尖损失增大,FoM 下降 0.02~0.05。
因此,eVTOL螺旋桨的桨尖速度通常限制在 Ma 0.55~0.65(对应海平面 Vtip ≈ 187~221 m/s)。对于分布式小直径桨(D=0.5~1.0 m),这意味着转速通常在 3000~6000 RPM 范围。若追求低噪声,可进一步降至 Ma 0.45~0.50,但需增大桨盘面积或增加桨叶数以补偿推力损失。
4.4 载荷分布与桨叶优化
4.4.1 最佳载荷分布理论
根据Betz最小诱导损失条件,螺旋桨在给定推力下达到最小诱导功率时,桨盘平面内的轴向诱导速度应为常数。由此可推导出最佳环量分布:
Γ(r) ∝ r / √(1 - (r/R)²)
该分布意味着桨尖附近环量趋于无穷大,但实际中受桨尖损失和结构强度限制,需在桨尖区域进行裁剪(tip unloading)。工程上常采用椭圆分布或带桨尖卸载的修正分布。
4.4.2 弦长与扭转分布设计
给定目标推力系数 CT 和桨尖速度,可通过BEMT反推最佳弦长分布 c(r) 和扭转分布 θ(r):
- 设定目标载荷分布(如椭圆分布或自定义卸载方案)。
- 迭代求解各叶素处的当地攻角 α(r),使升力系数 Cl(r) 接近翼型最高升阻比对应值(通常 Cl ≈ 0.6~0.8)。
- 计算弦长 c(r) = 2πr·F·vi / (N·W·Cl)。
- 计算扭转 θ(r) = φ(r) + α(r),其中 φ 为来流角。
典型分布式桨叶的弦长分布呈“根部宽、中部窄、桨尖略宽”的形态,扭转从根部约 30°~40° 线性减小至桨尖约 5°~10°。
4.4.3 载荷分布对噪声与结构的影响
载荷分布不仅影响效率,还直接影响噪声和结构强度:
- 桨尖载荷集中会增加桨尖涡强度,导致桨-涡干扰噪声(BPF及其谐波)显著增强。
- 根部载荷过大会增加桨毂弯矩,对电机轴和轴承寿命不利。
- 展向载荷波动会引起振动和疲劳问题,尤其在分布式多桨布局中,各桨载荷不一致会加剧整机振动。
4.5 气动干扰与分布式布局
4.5.1 桨间干扰
当多个螺旋桨近距离布置时,相邻桨叶的尾流相互卷吸,形成干扰。主要效应包括:
- 尾流收缩与加速:相邻桨的尾流相互吸引,导致有效桨盘面积减小,诱导功率增加。
- 桨尖涡相互作用:相邻桨的桨尖涡相互缠绕,加剧涡的耗散,但也会产生周期性载荷波动。
- 噪声叠加:多桨噪声在特定频率上相干叠加,可能形成指向性噪声增强。
定量研究表明,当桨间距(桨尖到桨尖)小于 0.3R 时,效率损失可达 5%~8%;间距大于 0.5R 时,损失降至 2% 以内。因此,分布式布局应尽量保持桨间距大于 0.5R。
4.5.2 机翼/机身干扰
分布式推进常将螺旋桨布置在机翼前缘或翼尖,产生有益的气动干扰:
- 螺旋桨滑流增升:桨后滑流速度增加,动压增大,使机翼局部升力增加,可降低起飞/着陆速度。
- 翼尖桨消涡:翼尖螺旋桨的旋转方向若与翼尖涡相反,可削弱翼尖涡强度,降低诱导阻力。
- 不利干扰:桨盘对机翼的遮挡效应会改变机翼表面压力分布,需通过CFD或风洞试验验证。
4.5.3 地面效应与机库效应
在起降阶段,地面效应会显著改变螺旋桨性能:
ΔP/P ≈ -0.1 ~ -0.15(当桨盘离地高度 h/R < 1.5 时)
地面效应使悬停功率降低 10%~15%,但也会带来“吸底”风险——当桨盘离地过近(h/R < 0.5)时,桨盘下方气流受地面阻塞,桨盘载荷分布畸变,可能导致推力脉动和噪声增大。
4.6 案例拆解:某四轴eVTOL螺旋桨初步设计
背景与约束:某城市空中交通(UAM)eVTOL,最大起飞重量 2200 kg,采用 4 个倾转旋翼(每个旋翼承担 550 kg 推力),巡航速度 240 km/h,悬停升限 1500 m(ρ≈1.058 kg/m³)。要求悬停噪声低于 75 dB(150 m 处),桨尖马赫数不超过 0.62。
初始参数设定:
- 单桨推力 T = 550×9.81 ≈ 5395 N。
- 选取桨盘载荷 DL = 60 N/m²,则桨盘面积 A = T/DL ≈ 90 m²,直径 D = √(4A/π) ≈ 3.38 m。但四轴布局下单个桨直径过大,结构重量和收纳困难。改为 8 桨布局(每个桨承担 275 kg),则单桨 T ≈ 2698 N,DL=60 时 A≈45 m²,D≈2.39 m。
- 桨尖速度限制:Vtip ≤ 0.62×340 ≈ 211 m/s,则最大转速 n = Vtip/(πD) ≈ 211/(π×2.39) ≈ 28.1 rps ≈ 1686 RPM。
BEMT迭代计算:采用 3 片桨叶,翼型选用 E387,目标 Cl=0.7,初始弦长分布按椭圆分布设定。迭代 15 次后收敛,得到:
- 悬停功率 P_hover ≈ 98 kW/桨,整机悬停功率 ≈ 784 kW。
- 品质因子 FoM ≈ 0.74。
- 桨尖载荷卸载至峰值的 65%,桨根载荷增加 15%。
气动干扰校核:8 桨布局中,相邻桨间距为 2.8 m(约 1.17R),桨间干扰损失约 2%。机翼滑流增升效应使巡航升阻比提升 8%,但巡航时桨叶倾转至水平位置,滑流对机翼的增升作用减弱,需通过倾转过渡段控制律补偿。
结果对比与复盘:初始 4 桨方案因桨盘过大导致结构超重(预估增加 120 kg),改为 8 桨后结构重量增加 60 kg,但悬停效率提升 5%,噪声降低 3 dB。最终选择 8 桨方案。复盘经验:分布式推进的桨数选择需在气动效率、结构重量、噪声和系统冗余之间综合权衡,不能单纯追求最小桨盘载荷。
4.7 总结与设计检查清单
本章核心要点可归纳为以下设计检查清单:
| 检查项 | 设计目标 | 验证方法 |
|---|---|---|
| 桨尖马赫数 | ≤0.62(海平面),≤0.58(高原) | 计算 Vtip/(a) |
| 品质因子 FoM | ≥0.70 | BEMT或CFD |
| 桨盘载荷 DL | 40~80 N/m² | T/A |
| 桨间距 | ≥0.5R | 几何检查 |
| 桨尖卸载 | 桨尖载荷≤峰值70% | BEMT输出 |
| 地面效应裕度 | 悬停功率裕度≥15% | 计入地效修正 |
在实际工程中,建议采用“动量理论初算 → BEMT详细设计 → CFD校核 → 风洞试验验证”的四级设计流程,每一级都需与电机选型、结构强度和噪声要求迭代闭环。分布式推进的气动设计不是孤立的螺旋桨设计,而是与整机布局、飞控策略深度耦合的系统工程。