历史人物：古斯塔夫·埃菲尔与风洞实验
高级话题：风工程与极端天气设计
本章小结
练习题
常见陷阱与错误
钝体空气动力学是流体力学中极其重要yet经常被低估的分支。与追求流线型的航空器不同，我们日常生活中遇到的大部分物体——汽车、建筑物、桥梁、运动器材——都属于钝体范畴。理解钝体绕流不仅能帮助我们设计更安全、更高效的产品，还能解释许多看似反直觉的现象。
想象一下F1赛车尾翼产生的下压力、摩天大楼在强风中的摇摆、跳伞运动员稳定下降的姿态，这些都是钝体空气动力学的生动体现。本章将通过这些实例，建立对钝体绕流的物理直觉。
钝体最显著的特征是流动分离不可避免。不同于流线型物体能让流体平滑地绕过，钝体的尖锐边缘或大曲率表面必然导致流动分离：
流线型翼型 (Re = 10^6):
→→→→→→→→→
╱ ╲
→ →
╲ ╱
→→→→→→→→→
分离点在尾缘附近
钝体圆柱 (Re = 10^6):
→→→┐ ┌→→→
→→→ │ │ ←←←
→→→ │ │ ←←←
→→→ │ │ ←←←
→→→┘ └→→→
分离点在90°附近
形成卡门涡街
钝体阻力的构成与流线型体截然不同：
这个差异源于分离区的低压：
\[C_D = C_{D,pressure} + C_{D,friction} \approx C_{D,pressure}\]典型阻力系数对比：
钝体绕流对雷诺数的依赖性呈现非单调特征，最著名的例子是圆柱体的阻力危机：
C_D
↑
1.2 |-------- 亚临界区
| \
0.3 | \ 阻力危机
| \_____ 超临界区
|
+------------------→ log(Re)
10^4 10^5 10^6
这种突变源于边界层转捩位置的改变：
汽车空气动力学设计在过去一个世纪经历了巨大变革：
年代 车型示例 C_D值
1920s 箱型车 0.7-0.9
1930s 流线型概念车 0.5-0.6
1960s 大众甲壳虫 0.48
1980s 奥迪100 0.30
2000s 丰田普锐斯 0.25
2020s 特斯拉Model S 0.21
理论极限 水滴形 0.04
中国高铁CRH380A的头部设计是钝体减阻的杰作：
侧视图：
╱╲_______________
╱ \_____________
/ \___________
│ │__________
│ │
关键设计要素：
1. 长细比 L/D ≈ 3.5 (头部长度/车宽)
2. 双曲率过渡避免流动分离
3. 底部导流板减少地面效应
4. 转向架整流罩
设计目标的权衡：
经验公式（300km/h运行速度）： \(P_{aero} = \frac{1}{2}\rho V^3 C_D A \approx 2000 \text{ kW}\)
汽车尾部设计对总阻力的影响甚至超过头部：
不同尾部形状的流动模式：
方背式 (C_D ≈ 0.45)： 快背式 (C_D ≈ 0.35)：
┌─────┐ ╱─────╲
│ │←涡环 │ ╲←附着
│ │ │ ╲
└─────┘ └────────╲
阶背式 (C_D ≈ 0.32)：
┌─────╲
│ ╲←再附着点
│ ├─╲
└───────┘ ╲
主动流动控制：
横风作用下的侧向力和偏航力矩：
\[\begin{aligned} F_y &= \frac{1}{2}\rho V_{rel}^2 C_y A_{side} \\ M_z &= \frac{1}{2}\rho V_{rel}^2 C_{M_z} A_{side} L \end{aligned}\]其中相对风速：$V_{rel} = \sqrt{V_{vehicle}^2 + V_{wind}^2}$
经验准则：
建筑物表面的风压分布遵循伯努利原理的变形：
\[C_p = \frac{p - p_\infty}{\frac{1}{2}\rho V_\infty^2}\]典型矩形建筑的压力系数分布：
平面图（风从左吹来）：
+0.8 │ -0.5 -0.5 │ -0.3
←←← ┌──────────┐
←←← │ │ 尾流区
←←← │ │ C_p ≈ -0.2
←←← └──────────┘
+0.8 │ -0.5 -0.5 │ -0.3
迎风面：正压 C_p = +0.6 ~ +0.8
侧面：负压 C_p = -0.5 ~ -0.7
背风面：负压 C_p = -0.2 ~ -0.4
屋顶边缘：强负压 C_p = -1.5 ~ -2.0
高层建筑的涡脱落频率由斯特劳哈尔数确定：
\[St = \frac{f \cdot D}{V} \approx 0.2\]对于方形截面建筑：
共振条件：当涡脱落频率接近建筑自振频率时发生
实例：上海中心大厦（632m高）
建筑宽度 D = 40m
设计风速 V = 50m/s
涡脱落频率 f = 0.2 × 50/40 = 0.25 Hz
建筑一阶频率 f₁ ≈ 0.15 Hz
避免共振的措施：
1. 螺旋上升外形（扭转120°）
2. 破坏涡脱落的规律性
3. 增加结构阻尼（调谐质量阻尼器）
设计风载荷的计算：
\[F = C_f \cdot q \cdot A \cdot G\]其中：
50年一遇设计风速（中国标准）：
建筑风洞试验必须满足的相似条件：
边界层模拟：
大气边界层风速剖面：
高度z
↑
500m│ V(z) = V_ref(z/z_ref)^α
│ ╱
│ ╱ α = 0.16 (城市)
│ ╱ α = 0.28 (市中心)
│ ╱
│╱
0 └────────→ V
0 V_ref
降落伞通过最大化阻力系数实现减速：
圆形伞的阻力系数： \(C_D = 1.3 - 1.4 \text{ （完全张开）}\)
终端速度平衡： \(mg = \frac{1}{2}\rho V_t^2 C_D A\)
因此： \(V_t = \sqrt{\frac{2mg}{\rho C_D A}}\)
典型参数（运动跳伞）：
开伞过程的阶段分析：
时间轴 →
速度 ↓ 自由落体
55m/s │\
│ \ 开伞冲击
│ \ (8-10g)
│ \___
5m/s │ \___稳定下降
└────────────→ t
0 0.5s 2s 5s
开伞力峰值估算： \(F_{max} = C_{snatch} \cdot \frac{1}{2}\rho V_{deploy}^2 A\)
其中填充系数 $C_{snatch} ≈ 1.5-2.0$
降落伞的钟摆运动：
○ 伞盖
╱│╲
╱ │ ╲ 摆动周期：
╱ │ ╲ T ≈ 2π√(L/g)
╱ │ ╲ L = 伞绳长度
╱ │ ╲
●
跳伞者
稳定性因素：
不同用途的降落伞设计：
圆形伞（传统） 方形伞（滑翔伞） 带形伞（阻力伞）
___ _______ ══════
/ \ | | ║ ║
( ) | | ║ ║
\___/ |_______| ══════
C_D = 1.3 L/D = 3-4 C_D = 1.5-1.8
垂直降落 可控滑翔 飞机减速
高尔夫球表面的凹坑是空气动力学中最反直觉的设计之一：
飞行距离对比（初速70m/s，仰角15°）：
光滑球：约120米
凹坑球：约240米（翻倍！）
这个差异源于阻力危机的人为触发：
C_D
↑
0.5 |━━━━━ 光滑球
|
0.25 | ━━━━━ 凹坑球
|
+------------------→ V
20 40 60 m/s
凹坑的作用机理：
光滑球（Re = 10^5）： 凹坑球（Re = 10^5）：
→→→→╱╲←←←← →→→→╱╲←←←←
→→→╱ ╲←←← →→→╱ ╲←←←
→→→│ 大 │←←← →→→│ 小 │←←←
→→→╲尾流╱←←← →→→╲尾流╱←←←
→→→→╲╱←←←← →→→→╲╱←←←←
分离角：80° 分离角：120°
C_D ≈ 0.5 C_D ≈ 0.25
凹坑引发的微观涡流：
高尔夫球的后旋产生升力：
\[F_L = \frac{1}{2}\rho V^2 A C_L\]其中升力系数： \(C_L = \frac{4\pi r\omega}{V} \cdot f(\text{凹坑})\)
典型参数：
弹道差异：
高度(m)
↑
30│ ╱─╲ 后旋球
│ ╱ ╲ (凹坑+旋转)
20│ ╱ ╲
│╱ ╲
10│ ───────╲ 无旋球
│ ╱ ╲(仅凹坑)
0└──────────────────→ 距离(m)
0 50 100 150 200 250
高尔夫球规则限制：
设计参数的影响：
参数 影响 最优范围
─────────────────────────────────────
凹坑数量 覆盖率 350-450个
凹坑深度 涡强度 0.3-0.35mm
凹坑形状 流动分离 圆形/六边形
排列方式 对称性 正二十面体
覆盖率 有效面积 75-85%
经验关联式（凹坑球）： \(\frac{C_D(\text{dimpled})}{C_D(\text{smooth})} \approx 0.5 \cdot \left(1 + 0.1\log_{10}Re\right)\)
古斯塔夫·埃菲尔（Gustave Eiffel, 1832-1923）的科学贡献远超其著名的铁塔：
1889年：埃菲尔铁塔建成，高度300米
1909年：在铁塔上建立空气动力学实验室
1912年：建造首个大型风洞（Auteuil实验室）
埃菲尔风洞规格：
- 试验段直径：2m
- 最大风速：32 m/s
- 测力天平精度：0.1N
- 首次实现六分力测量
关键创新：
埃菲尔建立的经验公式（至今使用）：
风压公式： \(p = K \cdot V^2\) 其中K = 0.07（国际单位制）
阻力系数标定： 首次系统测量各种形状的$C_D$值：
他的名言：”风不是敌人，而是需要理解的朋友。”
极端风速下的非线性效应：
风速-压力关系的修正： \(p = \frac{1}{2}\rho V^2 \cdot (1 + I_u)^2\)
其中湍流强度$I_u = \sigma_u/\bar{V}$
台风眼墙特征：
风速剖面：
r (km)
0 10 20 30 40
│ │ │ │ │
V │ ╱╲
│ ╱ ╲ 最大风速半径
│ ╱ ╲ R_max ≈ 20km
│╱ ╲___
└───────────→
眼 眼墙 外围
龙卷风的独特载荷：
Enhanced Fujita Scale:
等级 风速(m/s) 典型破坏
EF0 29-38 树枝折断
EF1 38-49 屋顶瓦片脱落
EF2 49-60 大树连根拔起
EF3 60-74 列车脱轨
EF4 74-89 坚固房屋夷平
EF5 >89 钢筋混凝土严重破坏
现代风工程结合CFD与AI：
大涡模拟(LES)应用：
机器学习加速：
# 伪代码：风压预测
input: 建筑形状参数、风向、风速
model: CNN + LSTM混合网络
output: 表面压力分布时程
加速比: 10^4 ~ 10^5
钝体空气动力学揭示了一个重要真理：有时候，不完美反而更好。从高尔夫球的凹坑到建筑物的倒角，适度的”粗糙”往往能带来意想不到的性能提升。
| 应用场景 | 关键公式 | 典型参数 |
|---|---|---|
| 阻力计算 | $F_D = \frac{1}{2}\rho V^2 C_D A$ | $C_D$: 0.3-2.0 |
| 涡脱落频率 | $f = St \cdot V/D$ | $St ≈ 0.2$ |
| 风压系数 | $C_p = (p-p_\infty)/(\frac{1}{2}\rho V^2)$ | -2.0 ~ +1.0 |
| 终端速度 | $V_t = \sqrt{2mg/(\rho C_D A)}$ | - |
| Magnus升力 | $F_L = \frac{1}{2}\rho V^2 A C_L$ | $C_L$: 0.1-0.5 |
钝体空气动力学教会我们：
习题14.1 汽车阻力估算 一辆SUV迎风面积A=3m²，阻力系数$C_D=0.35$，在高速公路上以120km/h行驶。计算： (a) 空气阻力大小 (b) 克服空气阻力所需功率 (c) 若$C_D$降低到0.30，可节省多少功率？
提示：注意单位换算，空气密度取1.2 kg/m³
习题14.2 建筑物涡脱落 一座宽度D=40m的高层建筑，在V=20m/s的风速下发生涡脱落。 (a) 计算涡脱落频率 (b) 若建筑一阶自振频率为0.1Hz，是否存在共振风险？ (c) 共振风速是多少？
提示：使用斯特劳哈尔数St=0.2
习题14.3 降落伞设计 设计一个降落伞，使80kg的跳伞者以5m/s的安全速度着陆。 (a) 计算所需的伞面积（假设$C_D=1.3$） (b) 若采用圆形伞，直径应为多少？ (c) 从3000m高空跳下，大约需要多长时间着陆？
提示：海平面空气密度1.225 kg/m³
习题14.4 高尔夫球弹道分析 一个高尔夫球初速70m/s，发射角15°，后旋转速3000rpm。考虑Magnus效应： (a) 估算升力系数（球径42.7mm） (b) 无旋转时的飞行距离（仅考虑重力） (c) 定性描述旋转如何改变弹道
提示：Magnus升力系数可用$C_L ≈ 2\pi r\omega/V$估算
习题14.5 动车组会车压力波 两列CRH380动车组以350km/h相对速度会车，估算： (a) 车头处的动压 (b) 会车瞬间的压力脉动幅值 (c) 为什么高铁车窗必须密封？
提示：考虑流线型车头的压力系数约为0.8
习题14.6 风工程综合设计 设计一座300m高的超高层建筑，所在地50年一遇风速45m/s。要求： (a) 选择合适的截面形状，说明理由 (b) 估算顶部最大位移（假设悬臂梁模型） (c) 提出至少三种减振措施
提示：考虑涡激振动、舒适度、经济性
习题14.7 仿生设计思考题（开放题） 观察以下生物的形态，分析其空气动力学优势，并提出工程应用： (a) 翠鸟的流线型头部 (b) 鲨鱼皮的微小齿状结构 (c) 枫树种子的旋转下落
习题14.8 极端天气应对（综合分析题） 某沿海城市拟建设应急避难设施，需抵御17级超强台风（风速>60m/s）。 分析设计要点，包括： (a) 建筑外形选择 (b) 结构体系建议 (c) 附属设施防护
错误：直接用风洞小模型数据推算实际建筑 正确：考虑雷诺数效应，特别是在临界区附近
错误：所有物体都设计成流线型 正确：权衡空气动力学性能与功能需求
错误：仅考虑时均风载荷 正确：分析涡脱落、阵风响应等动态效应
错误：认为$C_D$是常数 正确：$C_D$随Re、表面粗糙度、来流湍流度变化
错误：仅追求最小阻力 正确：综合考虑升力、稳定性、噪声等多目标
错误：简单线性缩放 正确：考虑相似准则，特别是时间尺度