本文是 John D. Anderson《Fundamentals of Aerodynamics》（第六版）读书笔记系列的第六篇，覆盖教材第7章《有限翼理论》（Prandtl’s Lifting-Line Theory）。这是Anderson全书中最具工程实用价值的一章之一。

第六篇：有限翼理论与升力线——从二维到三维的跨越

一、为什么二维理论不够？

1.1 概念引入

薄翼理论给出了一个完美的二维翼型：Cl=2παC_l = 2\pi\alpha，CdC_d 很小，唯一的力是升力。但现实中的机翼是有限展长的——有两个端部（翼尖），空气可以从下表面通过翼尖流到上表面。

这个”翼端泄漏”效应会带来三个重要后果：

  1. 升力降低：在翼尖处升力降为零
  2. 诱导阻力出现：额外的阻力，在二维中不存在
  3. 翼尖涡形成：真实机翼后面拖出两条集中的涡

Anderson 将普朗特（Prandtl）在 1918-1919 年左右提出的升力线理论描述为”空气动力学史上最优雅的理论之一”。

1.2 物理机制

升力沿着展向分布——翼根处大，翼尖处小。展向升力梯度意味着环量沿着展向变化。

根据亥姆霍兹涡量定理，涡量不能在流体内部终止。所以：

  • 升力变化的地方 → 涡面从机翼后缘拖出
  • 在翼尖处 → 集中卷成翼尖涡（wingtip vortices）

这些拖出的自由涡在机翼处诱导了一个向下的速度（下洗，downwash），改变了当地迎角，从而产生了诱导阻力。

二、升力线理论的数学框架

2.1 核心假设

普朗特升力线理论基于以下假设：

  1. 机翼用一根在展向位置的升力线（lifting line）替代，通常是 1/4 弦线
  2. 环量 Γ(y)\Gamma(y) 沿展向分布，在翼尖处 Γ(±b/2)=0\Gamma(\pm b/2)=0
  3. 后缘拖出的自由涡形成一个涡面（vortex sheet）
  4. 流动是不可压缩的
  5. 机翼是平面直机翼（无后掠，无扭转）

2.2 下洗与诱导迎角

下洗速度 w(y)w(y) 由展向涡量分布在 Biot-Savart 定律下产生：

w(y0)=14π∫−b/2b/2dΓ/dyy0−y dy w(y_0) = \frac{1}{4\pi}\int_{-b/2}^{b/2} \frac{d\Gamma/dy}{y_0 - y}\,dy

这个速度使来流迎角减小了 αi\alpha_i：

αi(y)=−w(y)V∞ \alpha_i(y) = -\frac{w(y)}{V_\infty}

诱导迎角（Induced Angle of Attack） 使得实际迎角 αeff\alpha_{\text{eff}} 小于几何迎角 α\alpha：

αeff(y)=α(y)−αi(y) \alpha_{\text{eff}}(y) = \alpha(y) - \alpha_i(y)

二维翼型的升力由有效迎角决定，而不是几何迎角。

2.3 升力线积分方程

升力线理论的核心——一个积分方程将环量分布 Γ(y)\Gamma(y) 与机翼的几何参数连接起来：

Γ(y)=12c(y)V∞a0(y)[α(y)−αL=0(y)−14πV∞∫−b/2b/2dΓ/dηy−η dη] \Gamma(y) = \frac{1}{2} c(y) V_\infty a_0(y) \left[\alpha(y) - \alpha_{L=0}(y) - \frac{1}{4\pi V_\infty}\int_{-b/2}^{b/2} \frac{d\Gamma/d\eta}{y-\eta}\,d\eta\right]

这个方程是普朗特的升力线方程。Mathematically formidable, physically beautiful.

2.4 傅里叶级数解

Glauert 引入了一个变换——将展向位置用角度表示：

y=−b2cos⁡θ,0≤θ≤π y = -\frac{b}{2}\cos\theta, \quad 0 \le \theta \le \pi

环量分布用傅里叶级数表示：

Γ(θ)=2bV∞∑n=1∞Ansin⁡(nθ) \Gamma(\theta) = 2b V_\infty \sum_{n=1}^\infty A_n \sin(n\theta)

其中 AnA_n 是由机翼几何确定的傅里叶系数。

三、椭圆翼——最优雅的解析解

3.1 椭圆环量分布

如果机翼是椭圆环量分布——即 Γ(y)=Γ01−(2y/b)2\Gamma(y) = \Gamma_0 \sqrt{1 - (2y/b)^2}，或者相应地，弦长 c(y)c(y) 沿展向呈椭圆分布——那么积分方程有精确的解析解。

椭圆翼的环量分布：

Γ(θ)=2bV∞A1sin⁡θ \Gamma(\theta) = 2b V_\infty A_1 \sin\theta

只有第一项傅里叶系数非零。这意味着：

  • 下洗速度 ww 沿展向是常数
  • 诱导迎角 αi\alpha_i 沿展向是常数

3.2 椭圆翼的关键公式

对于椭圆翼：

CL=a01+a0/(πAR)α C_L = \frac{a_0}{1 + a_0/(\pi AR)}\alpha

其中 a0a_0 是二维升力线斜率（2π2\pi 每弧度），AR=b2/SAR = b^2/S 是展弦比。

对于 a0=2πa_0 = 2\pi：

CL=2πARAR+2α C_L = \frac{2\pi AR}{AR + 2}\alpha

诱导阻力系数：

CD,i=CL2πAR \boxed{C_{D,i} = \frac{C_L^2}{\pi AR}}

这是升力线理论最著名的工程结果：诱导阻力与升力的平方成正比，与展弦比成反比。

3.3 手算例子

一个展弦比 AR=8AR = 8 的椭圆翼，在迎角 α=5∘\alpha = 5^\circ 下的性能：

二维升力系数：

Cl=2π×(5×π/180)=0.548 C_l = 2\pi \times (5\times\pi/180) = 0.548

三维升力系数（注意降级）：

CL=2π×88+2×(5×π/180)=50.2710×0.0873=0.439 C_L = \frac{2\pi \times 8}{8+2} \times (5\times\pi/180) = \frac{50.27}{10} \times 0.0873 = 0.439

三维效应使得 CLC_L 比 ClC_l 低了约 20%。

诱导阻力：

CD,i=0.4392π×8=0.19325.13≈0.0077 C_{D,i} = \frac{0.439^2}{\pi \times 8} = \frac{0.193}{25.13} \approx 0.0077

如果 AR=4AR = 4（典型的四旋翼臂/短翼），相同的 α=5∘\alpha=5^\circ：

CL=2π×44+2×0.0873=25.136×0.0873=0.366 C_L = \frac{2\pi \times 4}{4+2} \times 0.0873 = \frac{25.13}{6} \times 0.0873 = 0.366 CD,i=0.3662π×4=0.13412.57≈0.0107 C_{D,i} = \frac{0.366^2}{\pi \times 4} = \frac{0.134}{12.57} \approx 0.0107

低展弦比的机翼不仅升力更低，诱导阻力也显著增大。这就是为什么滑翔机有极长的展长（AR>20AR > 20），而战斗机的短翼（AR≈3−5AR \approx 3-5）虽然机动性好，但诱导阻力很大。

四、展弦比效应——最重要的工程参数

4.1 展弦比对升力线斜率的影响

升力线斜率随展弦比的变化：

ARAR aa（升力线斜率，1/rad） aa（1/deg） 相对二维值
∞\infty（二维） 2π=6.282\pi = 6.28 0.1096 100%
20 5.72 0.0998 91%
10 5.24 0.0914 83%
8 5.03 0.0877 80%
6 4.71 0.0822 75%
4 4.19 0.0731 67%
2 3.14 0.0548 50%

4.2 展弦比对诱导阻力的影响

诱导阻力与展弦比反比：

ARAR CL=0.5C_L=0.5 时 CD,iC_{D,i} CL=1.0C_L=1.0 时 CD,iC_{D,i}
∞\infty（二维） 0 0
20 0.0040 0.0159
10 0.0080 0.0318
8 0.0099 0.0398
6 0.0133 0.0531
4 0.0199 0.0796
2 0.0398 0.159

4.3 工程启示

展弦比大的机翼（AR>10AR > 10）：

  • 升力线斜率接近二维值
  • 诱导阻力低
  • 适合长续航飞行（滑翔机、侦察无人机）
  • 缺点：结构重、大过载时翼根弯矩大

展弦比小的机翼（AR<5AR < 5）：

  • 升力严重降级
  • 诱导阻力大
  • 适合大过载机动（战斗机、竞速机）
  • 结构紧凑，刚度好

四旋翼臂作为”机翼”的情况：四旋翼的机身/机臂通常展弦比为 1-3，诱导阻力效应非常显著。这也是为什么高速截击机会采用固定翼机身或者整机设计而非简单的四旋翼构型。

五、非椭圆翼和 Oswald 效率因子

5.1 非椭圆翼的修正

对于非椭圆环量分布的机翼，诱导阻力公式修正为：

CD,i=CL2πeAR \boxed{C_{D,i} = \frac{C_L^2}{\pi e AR}}

其中 ee 是Oswald 效率因子（0<e≤10 < e \le 1）。

  • 椭圆翼：e=1e = 1
  • 矩形翼：e≈0.7−0.8e \approx 0.7-0.8
  • 梯形翼（锥度比 λ≈0.3−0.5\lambda \approx 0.3-0.5）：e≈0.8−0.9e \approx 0.8-0.9

5.2 典型机翼的 Oswald 效率因子

机翼平面形状 锥度比 λ\lambda ee 备注
椭圆翼 — 1.0 理论最佳
矩形翼 1.0 0.7-0.8 制造简单
锥形翼 0.3-0.5 0.85-0.92 平衡设计
后掠翼 — 0.8-0.9 跨声速适用
三角翼 — 0.6-0.8 超声速适用

六、完整概念地图

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
二维翼型 (无限展长)
├── Cl = 2πα (无三维效应)
├── Cd 很小
└── 无诱导阻力

     ↓ 截断成有限展长

三维机翼 (有限展长)
├── 翼尖涡出现 → 下洗速度 w
├── 诱导迎角 αi = arctan(w/V∞)
├── 有效迎角 αeff = α - αi
├── 升力降低 (CL < Cl)
├── 诱导阻力出现 (CD,i > 0)
└── 展向环量分布 Γ(y)

普朗特升力线理论
├── 傅里叶级数求解 Γ(y)
├── 椭圆翼解析解 (e=1)
├── CD,i = CL²/(π e AR)
└── a = a0/(1+a0/(πAR))

工程结论
├── 大展弦比 → 低诱导阻力 → 长续航
├── 小展弦比 → 高诱导阻力 → 高机动
└── 椭圆翼 → 理论最优 (但制造复杂)

七、核心公式速查卡

公式 含义 使用场景
AR=b2/SAR = b^2/S 展弦比 三维机翼几何特性
CL=a01+a0/(πAR)αC_L = \frac{a_0}{1+a_0/(\pi AR)}\alpha 三维升力系数 下洗修正后的升力
CD,i=CL2/(πeAR)C_{D,i} = C_L^2/(\pi e AR) 诱导阻力 考虑展弦比和效率
ee Oswald 效率因子 <1< 1 表示非椭圆修正
αi(y)=CL/(πAR)\alpha_i(y) = C_L/(\pi AR) 平均诱导迎角 三维机翼迎角修正
w(y)=−V∞αi(y)w(y) = -V_\infty\alpha_i(y) 下洗速度 诱导速度场

参考文献

  1. Anderson, J. D. (2010). Fundamentals of Aerodynamics (6th ed., Chapter 7). McGraw-Hill.
  2. Prandtl, L. (1921). “Applications of Modern Hydrodynamics to Aeronautics”. NACA Report No. 116.
  3. Prandtl, L., & Tietjens, O. G. (1934). Applied Hydro- and Aeromechanics. Dover.
  4. Glauert, H. (1926). The Elements of Aerofoil and Airscrew Theory. Cambridge University Press.
  5. Anderson, J. D. (2016). Introduction to Flight (8th ed., Chapter 6). McGraw-Hill.

下一节：有限翼实验与修正

真实机翼的复杂性远超升力线理论——后掠翼、三角翼、边条翼、翼梢小翼和地面效应等。本篇将介绍这些工程修正和普朗特后缘涡面法的进阶应用。