第三章 非定常气动力基础：Theodorsen理论、Wagner函数、Küssner函数
各位同行，今天我们来聊聊非定常气动力。说实话，这是气动弹性力学里最绕、也最核心的一块。你想想看，叶片在转，气流在变，力也在变——这可不是简单的“升力=动压×面积×升力系数”能搞定的。
我当年刚接触这个领域时，总觉得定常气动力已经够用了。直到有一次做某型风扇叶片的颤振分析，算出来稳定裕度挺高，结果试验时叶片在某个转速下直接拍断了。嗯，从那以后，我再也不敢小看非定常效应。
3.1 为什么需要非定常气动力？
先问个问题：叶片在振动时，周围的气流会怎么反应？
答案很简单——气流不会瞬间跟上叶片的运动。这种“滞后”效应，就是非定常气动力的根源。说白了，叶片动一下，气流要过一会儿才反应过来。这个“一会儿”里，力的幅值和相位都变了。
我个人的习惯是，把非定常气动力拆成两部分看：
- 准定常部分：叶片当前姿态对应的瞬时力
- 记忆效应部分：过去运动历史对当前力的影响
这两部分加在一起，才是真实的非定常气动力。你想想看，如果只算准定常部分，那颤振边界估计能偏出去30%以上。
3.2 Theodorsen理论：翼型振荡的经典解
Theodorsen理论是1935年提出的，距今快90年了，但至今仍是工程界的基石。它解决的是：一个薄翼型在不可压流中做简谐振荡时，升力和力矩怎么算。
我记得第一次推导Theodorsen公式时，被那些Bessel函数搞得头大。但后来发现，其实核心思想就两条：
- 翼型运动分解为沉浮和俯仰两个自由度
- 尾涡的诱导速度用Theodorsen函数C(k)来修正
这里k是减缩频率，k=ωb/U，b是半弦长，U是来流速度。k越大，非定常效应越强。
核心公式（升力系数）：
C_L = 2π [ (ḣ/U + α) · C(k) + (b/2U)·α̇ · C(k) + (b/2U²)·ḧ ]其中：h为沉浮位移，α为俯仰角，上点表示对时间求导
这里有个坑，我踩过。Theodorsen理论假设尾涡是平面的，且翼型做小振幅简谐运动。如果你算的是大攻角或非简谐运动，这个理论就不太准了。我曾经在一个项目中直接用Theodorsen算非简谐响应，结果跟CFD差了40%。
工程经验：Theodorsen理论最适合做颤振边界的初步扫参。我一般用它先筛一遍参数空间，再用CFD验证关键点。这样效率高，也不容易漏掉危险工况。
3.3 Wagner函数：阶跃响应的时域描述
Wagner函数描述的是：翼型突然改变攻角后，升力随时间的变化过程。说白了，就是给翼型一个阶跃输入，看它怎么响应。
你可能会问：为什么要研究阶跃响应？因为任何复杂的运动都可以分解成一系列阶跃的叠加。这就是Duhamel积分的思路。
Wagner函数Φ(s)的表达式长这样：
Φ(s) = 1 - 0.165·exp(-0.0455s) - 0.335·exp(-0.3s)其中s=2Ut/b，是无量纲时间。你看，这个函数从0开始，逐渐趋近于1。物理意义是：攻角突变后，升力不会立刻达到最终值，而是慢慢爬上去。
我个人的理解是：Wagner函数就像给系统加了个低通滤波器。高频成分被衰减，低频成分保留。这在做时域仿真时特别有用。
避坑指南：Wagner函数只适用于不可压流。如果你算的是高马赫数（Ma>0.3），必须用可压缩修正。我曾经在Ma=0.6时直接套用Wagner函数，结果颤振速度算出来偏低了15%。后来改用可压缩的Wagner函数才把结果拉回来。
3.4 Küssner函数：突风响应的核心工具
Küssner函数和Wagner函数是“孪生兄弟”。区别在于：Wagner是翼型自己动，Küssner是来流在动（比如突风）。
Küssner函数Ψ(s)的表达式：
Ψ(s) = 1 - 0.5·exp(-0.13s) - 0.5·exp(-s)你看，形式跟Wagner很像，但系数不同。物理上，突风引起的升力响应比攻角变化更快一些。为什么？因为突风直接改变了翼型表面的速度分布，而攻角变化还要通过环量调整来传递。
我在做发动机进气畸变分析时，经常用Küssner函数来评估突风对叶片载荷的影响。比如，进气道出口有总压畸变，相当于叶片周期性遭遇“突风”。用Küssner函数可以快速估算出叶片振动响应的幅值。
三个函数的对比：
| 函数 | 输入类型 | 应用场景 | 频率范围 |
|---|---|---|---|
| Theodorsen C(k) | 简谐振荡 | 颤振分析（频域） | 全频段 |
| Wagner Φ(s) | 阶跃攻角 | 时域响应、机动载荷 | 低频为主 |
| Küssner Ψ(s) | 阶跃突风 | 突风响应、畸变载荷 | 中低频 |
3.5 三个函数的内在联系
这三个函数不是孤立的。它们通过傅里叶变换和拉普拉斯变换联系在一起。简单说：
- Theodorsen函数是频域表达
- Wagner和Küssner是时域表达
- 三者可以通过积分变换互相转换
我画了一张图，帮你理清它们的关系：
你看这张图，三个函数就像三兄弟，各有专长但又血脉相连。实际工程中，我通常这样用：
- 做频域颤振分析时，用Theodorsen函数
- 做时域响应仿真时，用Wagner函数
- 做突风或畸变载荷时，用Küssner函数
但不管用哪个，都要记住它们的适用范围：不可压、无粘、小扰动。超出这个范围，就要考虑修正或换用更高级的方法。
个人建议：刚开始学的时候，别急着钻公式推导。先把三个函数的物理意义搞清楚：谁在动？响应是什么？频率/时间怎么转换？这些想明白了，公式自然就记住了。
好了，这一章的内容就到这里。非定常气动力是气动弹性分析的“心脏”，搞懂了它，后面的颤振分析、响应计算都会顺畅很多。下一章我们会把这些理论用到实际的叶片颤振分析中，到时候再细聊。
公众号:蓝海资料掘金营，微信deep3321