第16章：非定常流动现象

当你站在风中，感受到阵阵吹来的风力时强时弱；当你看到旗帜在风中猎猎作响；当你听到电线在风中发出呜呜的哨音——这些都是非定常流动的表现。与我们之前讨论的定常流动不同，非定常流动中的流场参数随时间变化，这种时变特性带来了许多独特而重要的现象。从工程灾难（如塔科马海峡大桥的倒塌）到生物进化的奇迹（如昆虫的飞行），非定常流动塑造着我们的世界。本章将探讨这些动态现象背后的物理机制，以及如何在工程设计中考虑和利用它们。

16.1 涡激振动

16.1.1 涡脱落的物理机制

当流体绕过钝体（如圆柱）时，在一定的雷诺数范围内（约 $40 < Re < 10^5$），会发生周期性的涡脱落现象。这个过程可以这样理解：

流动方向 →
        ┌─────┐
        │     │  ← 圆柱
        │     │
        └─────┘
           ↓
    ╭──╮    ╭──╮
    │  │    │  │  ← 交替脱落的涡
    ╰──╯    ╰──╯
       ↘    ↙
        涡街

涡脱落的过程：

边界层分离：流体在圆柱背面分离，形成两个剪切层
剪切层卷起：上下剪切层交替卷起形成涡
涡的脱落：当涡增长到一定强度，它会从物体脱落
周期性重复：这个过程周期性地在上下两侧交替发生

16.1.2 斯特劳哈尔数

涡脱落频率由无量纲的斯特劳哈尔数（Strouhal number）描述：

$$St = \frac{fD}{V}$$ 其中：

$f$ 是涡脱落频率（Hz）
$D$ 是特征尺寸（如圆柱直径）
$V$ 是来流速度

对于圆柱绕流，在很宽的雷诺数范围内，$St \approx 0.2$。这个经验值非常有用：

实例估算：一根直径10cm的烟囱在10m/s的风中： $$f = \frac{St \cdot V}{D} = \frac{0.2 \times 10}{0.1} = 20 \text{ Hz}$$

16.1.3 锁频现象

当结构的固有频率接近涡脱落频率时，会发生"锁频"（lock-in）现象：

频率捕获：涡脱落频率被"拉"到结构振动频率
振幅放大：结构振动与流动形成正反馈
共振区间：锁频发生在 $0.8f_n < f_v < 1.2f_n$ 范围内

锁频时的振幅可以用下式估算： $$\frac{A}{D} \approx \frac{1}{2\zeta} \cdot \frac{\rho V^2 D^2}{\rho_s t^2 f_n^2 D^2}$$ 其中 $\zeta$ 是阻尼比，$\rho_s$ 是结构密度，$t$ 是壁厚。

16.1.4 工程实例

海洋立管：深海石油平台的立管在洋流作用下会发生涡激振动。设计时需要：

增加阻尼（如安装阻尼器）
改变涡脱落模式（螺旋列板）
避开共振区（调整张力改变固有频率）

高层建筑和烟囱：圆形截面的高层建筑特别容易发生涡激振动。常用对策：

螺旋形装饰条（破坏涡的相关性）
调谐质量阻尼器（TMD）
截面形状优化（如椭圆形或多边形）

经验法则：

当 $V_r = V/(f_n D) > 5$ 时，需要考虑涡激振动
螺旋列板的螺距约为 $5D$，高度约为 $0.1D$
增加5%的结构阻尼可减少50%的振动幅度

16.2 颤振与驰振

16.2.1 颤振的物理机理

颤振（Flutter）是一种由气动力、弹性力和惯性力相互作用引起的自激振动。以机翼为例：

初始扰动 → 机翼扭转 → 攻角变化 → 升力变化
    ↑                              ↓
    └──── 能量输入 ←─── 相位关系合适 ←┘

颤振的关键在于运动与气动力之间的相位关系。当系统从流动中吸收的能量大于结构阻尼耗散的能量时，振动会不断增大。

16.2.2 临界颤振速度

经典的颤振速度估算（对于机翼）： $$V_f \approx 2\pi f_\alpha \sqrt{\frac{I_\alpha}{\rho b^3 S_{C_L}}}$$ 其中：

$f_\alpha$ 是扭转固有频率
$I_\alpha$ 是扭转惯量
$b$ 是半弦长
$S_{C_L}$ 是升力线斜率

简化判据：

频率比 $\omega_h/\omega_\alpha \approx 0.5-0.8$ 时最危险
质量比 $\mu = m/(\rho b^2) > 20$ 相对安全
重心在气动中心之前有助于稳定

16.2.3 驰振的条件

驰振（Galloping）是另一种气动不稳定现象，主要发生在非流线型截面：

Glauert-Den Hartog判据： $$\frac{dC_L}{d\alpha} + C_D < 0$$ 这意味着当升力系数随攻角的变化率（负值）超过阻力系数时，会发生驰振。

典型易发生驰振的截面：

方形、矩形（冰覆盖的输电线）
D形截面（某些朝向）
三角形截面

16.2.4 工程实例

悬索桥设计：现代悬索桥采用扁平箱梁截面，具有良好的气动稳定性：

宽高比 $B/D > 10$
导流板和风嘴
中央开槽（减少涡激振动）

输电线防舞动：

防舞器（改变截面形状）
阻尼器（增加系统阻尼）
相间间隔棒（限制振幅）

经验数据：

输电线舞动多发生在风速 5-15 m/s
冰厚 5-20 mm 时最危险
振幅可达导线直径的 100-300 倍

16.3 拍动翼空气动力学

16.3.1 昆虫飞行的独特机制

昆虫飞行打破了传统定常空气动力学的限制。一只蜜蜂按照定常理论计算，其翅膀面积完全不足以产生支撑体重的升力。关键在于非定常机制：

拍动周期：
     ↑ 上拍
    ╱ ╲     前缘涡
   ╱   ╲    附着
  ╱     ╲
 ╱_______╲  翻转
     ↓ 下拍

16.3.2 前缘涡机制

拍动翼在大攻角下产生强烈的前缘涡（Leading Edge Vortex, LEV）：

涡的形成：急剧的攻角变化在前缘产生分离涡
涡的稳定：翼尖流动和科氏力稳定前缘涡
升力增强：前缘涡产生的低压区大幅增加升力

升力系数可达： $$C_L \approx 2-4$$ 远超定常流动的失速值（$C_L \approx 1.2$）

16.3.3 克拉普机制

"拍手-投掷"（Clap and Fling）机制，常见于小型昆虫：

拍手阶段：两翅在背部接触，挤出中间的流体
投掷阶段：两翅快速分开，形成强烈的环量

这种机制可以产生瞬时的高升力，特别适合悬停和起飞。

16.3.4 尺度效应

雷诺数对拍动翼性能影响巨大：

飞行者	翼展(mm)	频率(Hz)	Re数	主导机制
果蝇	3	200	100	粘性主导
蜜蜂	15	230	1000	前缘涡
蜂鸟	100	25	5000	混合机制
鸽子	300	5	50000	准定常

16.3.5 仿生应用

微型飞行器（MAV）：

翼展 < 15cm
采用拍动翼可提高效率30-50%
关键挑战：驱动机构和控制

设计准则：

最优频率：$f \propto L^{-1}$（L是特征长度）
拍动幅度：60°-120°
攻角变化：同步于位置变化，相位领先约90°

16.4 脉动流的特征

16.4.1 周期性流动参数

脉动流中，速度可表示为： $$V(t) = V_0 + V_1 \sin(\omega t)$$ 关键无量纲参数：

沃默斯利数（Womersley number）： $$Wo = R\sqrt{\frac{\omega}{\nu}}$$ 表征非定常惯性力与粘性力之比。

速度振幅比： $$\beta = \frac{V_1}{V_0}$$

16.4.2 准定常假设的适用性

准定常假设认为瞬时流场等同于相同瞬时条件下的定常流场。

适用条件：

斯特劳哈尔数 $St < 0.1$
或减缩频率 $k = \omega L/V < 0.1$

当 $St > 1$ 时，流动完全由非定常效应主导。

16.4.3 相位滞后

在脉动流中，压力梯度与流量之间存在相位差： $$\Delta \phi = \arctan\left(\frac{Wo^2}{8}\right)$$

这种相位差在：

$Wo < 1$：粘性主导，相位差小
$Wo > 10$：惯性主导，相位差接近90°

16.4.4 生理流动实例

心血管系统：

主动脉：$Wo \approx 15$（惯性主导）
毛细血管：$Wo < 0.01$（粘性主导）
脉搏波速度：4-10 m/s

呼吸系统：

气管：$Wo \approx 2$（过渡区）
肺泡：$Wo < 0.1$（准定常）
潮气量：约500 mL

设计人工器官时必须匹配这些特征时间尺度。

16.5 历史人物：塔科马海峡大桥倒塌事件（1940）

16.5.1 "舞动格蒂"的故事

1940年11月7日，华盛顿州塔科马海峡大桥在仅19m/s（约68km/h）的风速下发生剧烈振动并最终倒塌。这座桥因其建成后就经常在风中起伏而被昵称为"舞动格蒂"（Galloping Gertie）。

大桥的关键参数：

主跨：853米（当时世界第三）
宽度：11.9米
高度：2.4米（非常薄）
宽高比：仅5:1（现代标准要求>10:1）

16.5.2 倒塌原因分析

长期以来，教科书将其归因于共振，但实际机理更复杂：

初期（建成-1940.11）：

竖向弯曲振动（频率约0.13Hz）
振幅1-1.5米
涡激振动主导

倒塌当天：

模态转换：从竖弯转为扭转（频率约0.2Hz）
扭转振幅达到45°
负气动阻尼导致颤振

关键因素：

截面形状：H型截面气动性能差
刚度不足：扭转刚度严重不足
气动耦合：扭转与竖弯模态耦合

16.5.3 工程教训

这次事故彻底改变了桥梁设计理念：

设计改进：

采用气动优化的箱型截面
增加扭转刚度
风洞试验成为必需
考虑气动弹性效应

现代标准：

检验风速等级：

- 施工阶段：10年一遇
- 运营阶段：100年一遇
- 颤振检验：1.2×100年一遇

16.5.4 科学遗产

塔科马大桥事故推动了多个学科的发展：

气动弹性理论：Theodorsen和Scanlan的开创性工作
风工程学科：成为独立的工程分支
计算方法：CFD-CSD耦合分析
实验技术：节段模型和全桥气弹模型试验

这个事件提醒我们：大自然的力量常常超出我们的直觉，谦逊和严谨的科学态度是工程安全的基石。

16.6 高级话题：流固耦合与气动弹性

16.6.1 流固耦合的分类

根据耦合强度和方式，流固耦合（FSI）可分为：

单向耦合：

流场影响结构，但结构变形不影响流场
适用于小变形（<5%特征尺寸）
计算简单，精度有限

双向耦合：

流场与结构相互影响
需要迭代求解
计算量大但精度高

强耦合vs弱耦合：

弱耦合：交替求解，显式时间推进
强耦合：联立求解，隐式格式
判据：$$\frac{\rho_f}{\rho_s} \cdot \frac{L}{h} > 1$$ 需要强耦合

16.6.2 气动弹性三角

气动弹性现象由三种力的相互作用决定：

        空气动力
           /\
          /  \
         /    \
        /  AE  \
       /________\
    弹性力    惯性力

不同区域的现象：

A-E边：静气动弹性（发散、反效）
A-I边：非定常空气动力
E-I边：结构动力学
中心：颤振、抖振等动气动弹性

16.6.3 降阶模型方法

对于复杂的FSI问题，降阶模型（ROM）提供了高效的解决方案：

本征正交分解（POD）：将流场分解为基本模态： $$\mathbf{u}(\mathbf{x},t) = \bar{\mathbf{u}}(\mathbf{x}) + \sum_{i=1}^{N} a_i(t)\boldsymbol{\phi}_i(\mathbf{x})$$ 系统识别方法：

自回归模型（AR）
Volterra级数
神经网络方法

16.6.4 工程应用实例

飞机机翼设计：

静气动弹性：机翼扭转影响载荷分布
颤振边界：确定飞行包线
阵风响应：疲劳寿命评估

风力发电机叶片：

尺寸：长度可达100米
挑战：重力、离心力、气动力耦合
关键：避免失速颤振和古典颤振

高速列车受电弓：

速度：350km/h
问题：弓网耦合振动
解决：主动控制和优化设计

16.6.5 数值方法进展

分区耦合策略：

CFD求解器 ←→ 接口 ←→ CSD求解器
    ↓         ↓         ↓
  流场网格   数据传递   结构网格

时间推进格式：

松弛因子：$\alpha \in [0.3, 0.7]$
子迭代：3-5次通常足够
时间步长：$\Delta t < 0.1T_{min}$

网格变形技术：

弹簧类比法
径向基函数（RBF）
任意拉格朗日-欧拉（ALE）

16.7 本章小结

非定常流动现象展现了流体力学的动态美感和复杂性。本章的关键要点：

核心概念

涡激振动： - 斯特劳哈尔数 $St \approx 0.2$ 对圆柱 - 锁频发生在 $0.8f_n < f_v < 1.2f_n$ - 螺旋列板可有效抑制
颤振与驰振： - 颤振：自激振动，需要气动-弹性耦合 - 驰振：Glauert-Den Hartog判据 - 临界风速是设计关键
拍动翼机制： - 前缘涡提供高升力（$C_L = 2-4$） - 尺度效应决定主导机制 - 为MAV设计提供灵感
脉动流特征： - 沃默斯利数判断流态 - 准定常假设：$St < 0.1$ - 相位关系影响能量传递

关键无量纲数

| 参数 | 定义 | 物理意义 | 临界值 |

参数	定义	物理意义	临界值
$St$	$fD/V$	非定常性	0.1-1
$Wo$	$R\sqrt{\omega/\nu}$	脉动流特征	1, 10
$k$	$\omega L/V$	减缩频率	0.1
$V_r$	$V/(f_n D)$	约化速度	5

工程启示

非定常效应可能是破坏性的（塔科马大桥），也可能是有益的（昆虫飞行）
流固耦合在现代工程中无处不在
正确的时间尺度分析是理解非定常现象的关键
实验验证（特别是风洞试验）仍然不可替代

16.8 练习题

基础题

16.1 涡脱落频率计算 一根直径20cm的圆形烟囱在风速15m/s的环境中，计算涡脱落频率。如果烟囱的固有频率为18Hz，是否会发生锁频现象？

提示

使用斯特劳哈尔数 St = 0.2 计算涡脱落频率，然后检查是否在锁频范围内。

答案

涡脱落频率： $$f = \frac{St \cdot V}{D} = \frac{0.2 \times 15}{0.2} = 15 \text{ Hz}$$ 锁频范围：$0.8 \times 18 = 14.4$ Hz 到 $1.2 \times 18 = 21.6$ Hz

由于15 Hz在[14.4, 21.6]范围内，会发生锁频现象，振动会被放大。

16.2 昆虫翅膀雷诺数 一只蜜蜂翅膀长度15mm，拍动频率230Hz，飞行速度3m/s。估算其雷诺数并判断主导的流动机制。（运动粘度 $\nu = 1.5 \times 10^{-5}$ m²/s）

提示

特征速度应该考虑拍动速度而非飞行速度。拍动速度 = 2πfA，其中A是拍动幅度。

答案

拍动幅度约为翼长：A ≈ 15mm = 0.015m 拍动速度：$V_{flap} = 2\pi f A = 2\pi \times 230 \times 0.015 \approx 21.7$ m/s

这远大于飞行速度，所以用拍动速度计算： $$Re = \frac{V_{flap} \cdot L}{\nu} = \frac{21.7 \times 0.015}{1.5 \times 10^{-5}} \approx 2170$$ 在Re ≈ 1000-5000范围，前缘涡机制主导。

16.3 脉动流相位差 在直径2cm的血管中，血流脉动频率为1.2Hz（心率72次/分），血液运动粘度 $\nu = 4 \times 10^{-6}$ m²/s。计算沃默斯利数并估算压力梯度与流速的相位差。

提示

先计算角频率ω，然后用Wo数公式，最后估算相位差。

答案

角频率：$\omega = 2\pi f = 2\pi \times 1.2 = 7.54$ rad/s 半径：R = 0.01 m

沃默斯利数： $$Wo = R\sqrt{\frac{\omega}{\nu}} = 0.01 \sqrt{\frac{7.54}{4 \times 10^{-6}}} = 13.7$$ 相位差： $$\Delta \phi = \arctan\left(\frac{Wo^2}{8}\right) = \arctan\left(\frac{13.7^2}{8}\right) = \arctan(23.5) \approx 87°$$ 压力梯度领先流速约87°，接近90°，说明惯性效应主导。

挑战题

16.4 悬索桥颤振分析 某悬索桥主梁宽30m，高3m，质量线密度8000 kg/m，扭转惯量 $I_\alpha = 1.2 \times 10^7$ kg·m²/m，扭转频率0.3Hz。估算颤振临界风速。（空气密度1.2 kg/m³）

提示

使用简化的颤振速度公式，考虑质量比和频率的影响。对于扁平箱梁，升力线斜率约为2π。

答案

半宽度 b = 15m 质量比：$\mu = \frac{m}{\rho b^2} = \frac{8000}{1.2 \times 15^2} = 29.6$

颤振速度估算： $$V_f \approx 2\pi f_\alpha b \sqrt{\frac{\mu}{2\pi}} = 2\pi \times 0.3 \times 15 \times \sqrt{\frac{29.6}{2\pi}}$$ $$V_f \approx 28.3 \times 2.17 = 61.4 \text{ m/s}$$

这相当于221 km/h的风速，属于超强台风级别。

16.5 涡激振动抑制设计 某海洋立管直径0.5m，长度100m，在1m/s的洋流中工作。设计螺旋列板参数来抑制涡激振动。分析列板对阻力的影响。

提示

螺旋列板的标准设计：螺距5D，高度0.1D，三股螺旋。阻力会增加20-30%。

答案

螺旋列板设计：

螺距：5 × 0.5 = 2.5 m
列板高度：0.1 × 0.5 = 0.05 m
螺旋数：3股，相位差120°
螺旋角：$\theta = \arctan(2.5/(\pi \times 0.5)) \approx 58°$

效果分析：

破坏涡的相关长度，从100m降至约2.5m
振幅减少80-90%
阻力系数从1.2增至约1.5（增加25%）
年疲劳损伤降低100倍以上

权衡：阻力增加的代价远小于疲劳寿命延长的收益。

16.6 拍动翼MAV设计 设计一个翼展10cm的微型飞行器，重量5g。基于昆虫飞行原理，确定拍动频率、拍动幅度和所需功率。

提示

使用频率-尺度关系 f ∝ L^(-1)，参考蜜蜂的数据进行缩放。功率密度约为30-50 W/kg。

答案

基于尺度律（参考蜜蜂）：

蜜蜂：翼展30mm，频率230Hz
MAV缩放：$f = 230 \times (30/100) = 69$ Hz

设计参数：

拍动频率：70 Hz
拍动幅度：±60°（约±50mm行程）
翼面积：约25 cm²（双翼）
展弦比：4

升力估算：

需要升力：5g × 9.8 = 0.049 N
升力系数（时均）：$\bar{C_L} \approx 1.5$
验证速度：$V = \sqrt{\frac{2mg}{\rho S C_L}} \approx 2.5$ m/s

功率需求：

比功率：40 W/kg
总功率：40 × 0.005 = 0.2 W
考虑效率（30%）：需要0.7 W输入功率

这在现有电池和电机技术范围内可实现。

16.7 开放思考题：流固耦合的未来 随着计算能力的提升和新材料的发展，讨论以下问题：

智能材料（如形状记忆合金、压电材料）如何改变流动控制策略？
机器学习在预测和控制非定常流动中的潜力？
仿生设计中，哪些生物的流固耦合机制最值得深入研究？

思考方向

智能材料应用： - 自适应翼型：根据流动状态改变形状 - 主动涡发生器：按需激活 - 智能蒙皮：分布式传感和致动
机器学习潜力： - 降阶模型的自动构建 - 非线性系统识别 - 实时流动状态预测和控制
仿生研究重点： - 蝙蝠翼膜的被动适应性 - 鱼类侧线系统的流动感知 - 鸟类羽毛的多尺度流动控制

16.9 常见陷阱与错误

概念误区

共振 ≠ 颤振 - 错误：塔科马大桥是共振导致的倒塌 - 正确：颤振是自激振动，能量来自流动
涡脱落频率不是常数 - 错误：固定St数适用于所有情况 - 正确：St数随Re变化，特别在临界Re附近
准定常假设的滥用 - 错误：拍动翼可以用准定常理论分析 - 正确：St > 0.1时必须考虑非定常效应

计算错误

特征速度选择 - 错误：总是使用来流速度 - 正确：拍动问题用拍动速度，可能远大于飞行速度
相位关系忽略 - 错误：脉动流中压力和流速同相 - 正确：高Wo数时相位差接近90°
尺度效应忽视 - 错误：大尺度经验直接用于微尺度 - 正确：雷诺数相似才能保证机理相同

设计陷阱

单一模态考虑 - 错误：只检查基频共振 - 正确：高阶模态和模态耦合同样危险
阻尼的过度依赖 - 错误：增加阻尼总是有益的 - 正确：过大阻尼可能激发其他不稳定性
风洞试验的局限 - 错误：缩尺模型完全代表实际 - 正确：需要考虑雷诺数效应和弹性相似

16.10 最佳实践检查清单

设计阶段

[ ] 识别所有可能的非定常源
涡脱落（钝体）
流动分离（失速）
激波振荡（跨声速）
尾流干扰（多体）
[ ] 计算关键无量纲数
斯特劳哈尔数
约化速度
质量比
阻尼比
[ ] 频率分析
列出所有固有频率
计算激励频率范围
识别潜在共振区
考虑频率的环境变化（温度等）

分析验证

[ ] 稳定性边界确定
颤振速度计算
驰振判据检查
涡激振动评估
参数敏感性分析
[ ] 试验计划
风洞试验（刚性模型+弹性模型）
相似准则满足
测量点布置合理
数据采样率足够（>10倍最高频率）
[ ] 数值仿真策略
选择合适的湍流模型
时间步长满足CFL条件
网格独立性验证
FSI耦合方式选择

控制措施

[ ] 被动控制方案
气动外形优化
质量/刚度调整
阻尼器设计
涡流发生器布置
[ ] 主动控制考虑
传感器布置
执行器响应时间
控制律鲁棒性
失效安全设计
[ ] 监测与维护
振动监测系统
疲劳寿命评估
定期检查计划
应急响应预案

文档记录

[ ] 设计依据完整
所用标准和规范
关键假设说明
安全系数选择理由
与类似项目对比
[ ] 风险评估报告
失效模式分析
概率风险评估
缓解措施效果
残余风险说明

记住：非定常流动现象往往是灾难性失效的根源，但也可能是创新设计的灵感。理解其物理本质，合理利用和控制，是流体工程师的核心能力。