本文是 John D. Anderson《Fundamentals of Aerodynamics》（第六版）读书笔记系列的第四篇，覆盖教材第5章《翼型空气动力学》。阅读本文前需熟悉环量、库塔条件和库塔-儒科夫斯基定理。

第四篇：翼型空气动力学与薄翼理论

一、翼型的几何语言

1.1 概念引入

势流理论推导出了升力产生的数学本质——环量。但环量只适用于圆柱或简单的儒科夫斯基翼型。真实的机翼截面（翼型）有各种各样的形状，它们的几何参数直接影响气动性能。

Anderson 在教材中给出了翼型几何的完整语言体系。

1.2 翼型的几何要素

一个翼型的完整几何由以下要素定义：

弦线（Chord Line）：连接前缘（leading edge, LE）与后缘（trailing edge, TE）的直线。长度称为弦长（chord），记为 $c$ 。

中弧线（Mean Camber Line）：翼型上下表面之间中点的连线。如果中弧线是直线，称为对称翼型（symmetrical airfoil）；如果是曲线，称为弯翼型（cambered airfoil）。

厚度分布（Thickness Distribution）：垂直于中弧线的方向，上下表面距中弧线的距离。最大厚度与弦长之比称为最大相对厚度 $(t/c)_{\max}$ 。

前缘半径（Leading-edge Radius）：前缘处的曲率半径，决定了翼型在低迎角和大迎角下的性能。

后缘角（Trailing-edge Angle）：后缘处的夹角，影响库塔条件的数值实现。

1.3 NACA 翼型序列——标准化的几何体系

NACA（美国国家航空咨询委员会，NASA 的前身）在 1930-1940 年代建立了一套系统的翼型系列，用数字编码描述翼型的几何特征：

NACA 四位数字：NACA XXYY

第一位数字：最大弯度（% of chord）
第二位数字：最大弯度位置（tenths of chord，即十分之几弦长）
后两位数字：最大相对厚度（% of chord）

例：NACA 2412

最大弯度 2%（0.02c）
最大弯度位于 40% 弦长处（0.4c）
最大厚度 12%（0.12c）

NACA 五位数字：NACA 2XXYY

第一位数字：设计升力系数 $C_{l,\text{design}}$ 的 3/20 倍
第二、三位数字：最大弯度位置（% of chord 的一半）
后两位数字：最大厚度（% of chord）

例：NACA 23012

设计升力系数 $C_l = 0.3$ （ $2\times 0.15$ ）
最大弯度位置 15% chord（30/2）
最大厚度 12%

NACA 6 系列（层流翼型）：NACA 65-218

第一个数字表示系列
第二个数字表示最低压力点位置（tenths of chord）
减号后数字表示设计升力系数（tenths）
最后两位表示相对厚度

1.4 NACA 四位数翼型的数学公式

对于 NACA 四位数字翼型，几何可以用公式精确描述：

中弧线（ $\frac{dy_c}{dx}$ 表示斜率）：

y_c = \begin{cases} \frac{m}{p^2}\left(2p\frac{x}{c} - \left(\frac{x}{c}\right)^2\right), & 0 \le \frac{x}{c} \le p \\[6pt] \frac{m}{(1-p)^2}\left((1-2p)+2p\frac{x}{c} - \left(\frac{x}{c}\right)^2\right), & p \le \frac{x}{c} \le 1 \end{cases}

其中 $m$ 为最大弯度（%c）， $p$ 为最大弯度位置（%c）。

厚度分布（对称布置于中弧线两侧）：

y_t = \frac{t}{0.2}\left(0.2969\sqrt{\frac{x}{c}} - 0.1260\frac{x}{c} - 0.3516\left(\frac{x}{c}\right)^2 + 0.2843\left(\frac{x}{c}\right)^3 - 0.1015\left(\frac{x}{c}\right)^4\right)

其中 $t$ 为最大厚度。

这个公式有一条非常重要的特性：前缘附近的厚度变化与 $\sqrt{x}$ 成正比——这使得前缘是垂直的（曲率半径有限），从而准确描述了真实前缘的圆弧形状。

二、薄翼理论——用涡面解决升力问题

2.1 概念引入

库塔-儒科夫斯基定理给出了 $L' = \rho V_\infty \Gamma$ ，但其中环量 $\Gamma$ 是整个翼型的总环量。如果我想知道沿翼弦方向环量的分布——即各点的升力密度——呢？

薄翼理论（Thin Airfoil Theory）就是回答这个问题的理论。它的基本思想是：

用一个连续的涡面（vortex sheet）分布来替代翼型，通过求解这个涡分布的强度 $\gamma(s)$ 来得到翼型表面任意位置的压力和升力。

2.2 涡面概念

涡面（Vortex Sheet）：沿翼弦分布的一层连续的附着涡。每个位置 $s$ 处的环量密度为 $\gamma(s)$ ，单位是 $\text{m/s}$ 。

总环量：

\Gamma = \int_0^c \gamma(s)\,ds

根据库塔-儒科夫斯基定理，单位展长的升力为：

L' = \rho V_\infty \int_0^c \gamma(s)\,ds

而每个位置的局部升力密度为：

l'(x) = \rho V_\infty \gamma(x)

所以找到 $\gamma(x)$ 就找到了翼型上每点的升力分布。

2.3 涡面的边界条件

在涡面上，有一个关键的边界条件：涡面诱导的速度与来流叠加后，必须使得翼型表面是流线。

数学上，这意味着在翼型表面（近似在中弧线上）的法向速度为零：

V_{\infty,n} + w_n = 0

其中 $V_{\infty,n}$ 是来流的法向分量， $w_n$ 是涡面在翼型表面诱导的法向速度。

对于小迎角、薄翼型，边界条件可以线化：

\frac{1}{2\pi}\int_0^c \frac{\gamma(\xi)\,d\xi}{x-\xi} = V_\infty\left(\alpha - \frac{dy_c}{dx}\right)

这是薄翼理论的核心积分方程。左边是涡面诱导的法向速度（Biot-Savart 定律沿弦长积分），右边是来流的法向分量与中弧线斜率的组合。

2.4 薄翼理论的两个基本结果

Glauret 通过变换 $x = \frac{c}{2}(1-\cos\theta)$ ，将这个积分方程转化成了傅里叶级数的形式。解出 $\gamma(\theta)$ 后，得到两个极其重要的结果：

结果 1：升力系数

\boxed{C_l = 2\pi\left(\alpha - \alpha_{L=0}\right)}

其中 $\alpha_{L=0}$ 是零升力迎角——翼型不产生升力时的迎角。

对于对称翼型， $\alpha_{L=0} = 0$ ，所以 $C_l = 2\pi\alpha$ （与保角变换的结果一致）。

对于弯翼型， $\alpha_{L=0} < 0$ （通常是负数），这意味着即使在负迎角下，弯翼型也可以产生正向升力。

结果 2：力矩系数（关于前缘）

\boxed{C_{m,\text{LE}} = -\frac{C_l}{4} + \frac{\pi}{2}(\alpha_{L=0})_{\text{extra}}}

更重要的结论：关于 1/4 弦长点的力矩系数是常数：

C_{m,c/4} = \text{const}

这就是著名的气动中心（Aerodynamic Center）位于 1/4 弦长的结论。

2.5 手算例子：NACA 2412 翼型

以 NACA 2412 翼型为例（弯度 2%，最大弯度位置 40% chord，厚度 12%）。

求零升力迎角。

薄翼理论对 NACA 四位数字弯翼型的零升力迎角给出：

\alpha_{L=0} = -\frac{1}{\pi}\int_0^\pi \frac{dy_c}{dx} (1-\cos\theta)\,d\theta

代入 NACA 2412 的中弧线参数（ $m=0.02$ ， $p=0.4$ ）：

\alpha_{L=0} \approx -2m\cdot\frac{1-p}{1+p} = -2\times 0.02 \times \frac{0.6}{1.4} \approx -0.0171\ \text{rad} \approx -0.98^\circ

验证：

C_l(\alpha=0^\circ) = 2\pi(0 - (-0.0171)) = 2\pi \times 0.0171 \approx 0.107

物理意义：即使 NACA 2412 在 0° 迎角下，也能产生约 0.107 的升力系数——这就是弯度带来的”免费升力”。

2.6 薄翼理论的七条结论

Anderson 在书中总结了薄翼理论的七条重要结论：

#	结论	物理意义
1	$C_l = 2\pi(\alpha - \alpha_{L=0})$	升力线斜率恒为 $2\pi$
2	$\alpha_{L=0}$ 仅由弯度决定	对称翼型 $\alpha_{L=0}=0$
3	气动中心位于 1/4 弦长	力矩系数在此处为常数
4	$C_{m,c/4} = \text{const}$	与迎角无关
5	压力中心随 $C_l$ 增大向前缘移动	直到接近 1/4 弦长
6	翼型厚度对升力系数影响很小	主要影响最大升力系数和阻力系数
7	薄翼理论在小迎角下非常准确	$\alpha \lesssim 8^\circ$

三、翼型失速——理论与现实的差距

3.1 什么是失速

薄翼理论预测升力系数随迎角线性增长——斜率为 $2\pi$ （每弧度 0.1096/度）。但现实世界中，当迎角超过某个临界值后，升力系数不再增加，反而急剧下降。这个现象称为失速（Stall）。

失速的物理原因是边界层分离：翼型上表面的气流无法再沿着壁面流动，从某个位置”脱体”。

3.2 三种失速类型

根据翼型的不同几何特征，失速有三种类型：

1. 薄翼失速（Thin Airfoil Stall）

特征：前缘分离
发生在：薄翼型（注意：这里”薄”指前缘半径很薄，不是相对厚度小）
升力下降陡峭
典型于 NACA 四位数薄翼型、对称翼型

2. 厚翼失速（Thick Airfoil Stall）

特征：后缘分离逐渐发展
发生在：厚翼型（ $t/c > 15\%$ ）
升力下降平缓
典型于大型运输机机翼

3. 前缘分离/薄翼混合失速

特征：前缘短气泡分离后重新附着，然后从后缘发展
发生在：中等厚度翼型（ $8\% < t/c < 15\%$ ）
升力变化在两者之间

3.3 最大升力系数

翼型最重要的性能指标之一就是最大升力系数 $C_{l,\max}$ ：

翼型类型	$C_{l,\max}$ （典型值）	失速迎角
对称翼型（NACA 0012）	1.2-1.4	12°-16°
常规弯翼型（NACA 2412）	1.4-1.6	14°-18°
高升力翼型	1.8-2.2	18°-25°
带襟翼的翼型	2.5-3.5	20°-30°

3.4 雷诺数对失速的影响

四旋翼桨叶的雷诺数约为 $1.4 \times 10^4$ ，而波音客机机翼的雷诺数约为 $10^8$ 。雷诺数对 $C_{l,\max}$ 影响巨大：

低雷诺数效应（ $Re < 10^5$ ）：

$C_{l,\max}$ 显著下降（可能只有高雷诺数的 60-70%）
失速迎角减小
滞后效应明显（增迎角和减迎角时的升力曲线不同）

这是四旋翼桨叶设计中一个经常被低估的问题——桨叶工作在远低于大型飞机的雷诺数下，最大升力系数大打折扣。

四、高升力装置——增加 $C_{l,\max}$ 的工程方法

为了让飞机能在更低的速度下起飞和着陆，工程上开发了高升力装置（High-Lift Devices）：

装置类型	$C_{l,\max}$ 增量	原理
简单襟翼（Plain flap）	+0.5 ~ +0.8	增加有效弯度
开裂襟翼（Split flap）	+0.4 ~ +0.7	增大弯度 + 增加后缘压力
富勒襟翼（Fowler flap）	+0.8 ~ +1.5	增大弯度 + 增大有效面积
前缘缝翼（Slat）	+0.5 ~ +0.8	向后缘注入高能气流，延迟分离
双缝/三缝襟翼	+1.5 ~ +2.5	多重效应叠加

四旋翼没有传统的襟翼装置，但旋转桨叶的”前进/后退桨叶”效应在某种意义上类似——前进桨叶迎角大，后退桨叶迎角小，产生了不对称流动。

五、翼型力矩特性和气动中心

5.1 力矩特性

翼型不仅有升力，还会产生一个俯仰力矩（pitching moment），影响飞机的配平。

薄翼理论的扭矩结果：

对于对称翼型： $C_{m,c/4} = 0$

对于弯翼型： $C_{m,c/4} \neq 0$ ，且与迎角无关——这意味着气动中心在 $c/4$ 处。

5.2 压力中心

另一个重要的概念是压力中心（Center of Pressure, CP）——所有气动力的合力作用点。

对于对称翼型，压力中心在 $c/4$ 处（与迎角无关）。
对于弯翼型，压力中心随升力系数的增大向前缘移动。

六、完整概念地图

翼型几何
├── 弦线 c
├── 中弧线 (弯度)
├── 厚度分布
├── 前缘半径
└── 后缘角

NACA 翼型系列
├── 四位数 (弯度 + 位置 + 厚度)
├── 五位数 (设计 Cl + 弯度 + 厚度)
└── 6系列 (层流翼型)

薄翼理论
├── 涡面替代翼型
├── 核心积分方程
├── 傅里叶变换求解
├── Cl = 2π(α - αL=0)
├── 气动中心在 c/4
└── 七条重要结论

失速
├── 薄翼失速 (前缘分离)
├── 厚翼失速 (后缘分离)
└── 混合型

高升力装置
├── 襟翼 (Flap)
├── 缝翼 (Slat)
└── 多层组合

七、各概念在已发表文章中的出现

概念	出现文章	使用方式
翼型几何	《前进比完全解读》	螺旋桨桨叶翼型截面
$C_l$ 曲线	—	本文首次系统引入
失速	《前进比完全解读》《四旋翼飞行力学基础》	四旋翼高速飞行推力崩溃
高升力装置	—	本文首次引入
NACA 翼型	AirSim 源码深度解析	AirSim 中的翼型参数

八、核心公式速查卡

公式	含义	使用场景
$C_l = 2\pi(\alpha - \alpha_{L=0})$	薄翼升力系数	小迎角升力估算
$C_{m,c/4} = \text{const}$	气动中心力矩不变	配平计算
$\Gamma = \int_0^c \gamma(x)\,dx$	涡面总环量	升力分布
$\alpha_{L=0} \propto -\text{camber}$	零升力迎角与弯度成正比	弯翼型零升力

参考文献

Anderson, J. D. (2010). Fundamentals of Aerodynamics (6th ed., Chapter 5). McGraw-Hill.
Abbott, I. H., & Von Doenhoff, A. E. (1959). Theory of Wing Sections. Dover Publications.
Glauert, H. (1926). The Elements of Aerofoil and Airscrew Theory. Cambridge University Press.
Jacobs, E. N., Ward, K. E., & Pinkerton, R. M. (1933). “The Characteristics of 78 Related Airfoil Sections from Tests in the Variable-Density Wind Tunnel”. NACA Report No. 460.
Theodorsen, T. (1932). “Theory of Wing Sections of Arbitrary Shape”. NACA Report No. 411.

翼型特性与数据

从 NACA 风洞实测数据出发，深入分析各种翼型的 $C_l$ - $C_d$ 曲线、 $C_l$ - $C_m$ 曲线，比较不同类型翼型的性能差异（对称/弯翼型、层流翼型/常规翼型），给出工程选型的实用指南。