跳到主要内容
Sandbox Physics

非定常空气动力学 · 扑翼风洞

Flapping Wing

一片三维翼段在透明风洞中做正弦升沉和俯仰。改变频率、振幅、相位和转轴位置，观察环量升力、附加质量、交替脱落的尾迹涡，以及整周期推力与功率代理怎样同步改变。

交互模型Flapping Wing
运动相位0∘0^\circ
约化频率 kk0.180.18
Strouhal 数 StSt0.080.08
瞬时升力 CLC_L0.0000.000
平均推力代理 CT∗‾\overline{C_T^{\ast}}0.0000.000
平均功率代理 CP∗‾\overline{C_P^{\ast}}0.0000.000
推进效率代理 η∗\eta^{\ast}0%0\%
Kelvin 环量账本漂移 ϵΓ\epsilon_\Gamma0.0×10−80.0\times10^{-8}

物理教程

扑翼风洞：相位怎样把振荡变成推进

背景Theodorsen 在 NACA 496 中给出了振荡薄翼的线性非定常空气动力学，Garrick 随后在 NACA 567 中把同一框架用于扑翼推进。这个 Lab 保留环量力、附加质量和尾迹记忆，再把每次束缚环量变化记入可见的后缘尾迹。

为什么值得理解升沉或俯仰本身只是往复运动；真正决定尾迹方向、升力相位和平均流向力的，是两种运动的幅值、频率和相位组合。

先抓重点

核心问题
为什么同一片翼用同样的振幅和频率，只改变俯仰相位，就能从制造阻力切换到产生推进力？
一句话直觉
升沉速度改变有效攻角，Theodorsen 函数让环量升力相对运动滞后，而 Kelvin 定理要求束缚环量的变化以相反环量进入尾迹；因此相位 ϕ\phi 同时重排力与尾迹。

核心数学模型

升沉、俯仰与有效攻角

h(t)=h0sin⁡ωt,θ(t)=θ0sin⁡(ωt+ϕ),αeff=θ−tan⁡−1 ⁣(h˙U∞)h(t)=h_0\sin\omega t,\qquad \theta(t)=\theta_0\sin(\omega t+\phi),\qquad \alpha_{\mathrm{eff}}=\theta-\tan^{-1}\!\left(\frac{\dot h}{U_\infty}\right)

翼段向上运动时，它看到的相对来流向下偏转；俯仰可以增强或抵消这个运动攻角。

Theodorsen 尾迹记忆

C(k)≈1−0.165 ik0.0455+ik−0.335 ik0.3+ik,k=ωc2U∞C(k)\approx1-\frac{0.165\,ik}{0.0455+ik}-\frac{0.335\,ik}{0.3+ik},\qquad k=\frac{\omega c}{2U_\infty}

Jones–Wagner 有理近似把尾迹造成的幅值衰减和相位滞后压缩成一个可实时计算的复响应。

环量与附加质量升力

L′=πρb2 ⁣(−U∞θ˙−h¨+abθ¨)+2πρU∞b C(k) ⁣[U∞θ−h˙+b(12−a)θ˙]L'=\pi\rho b^2\!\left(-U_\infty\dot\theta-\ddot h+ab\ddot\theta\right)+2\pi\rho U_\infty b\,C(k)\!\left[U_\infty\theta-\dot h+b\left(\frac12-a\right)\dot\theta\right]

第一项是不需要稳态环量的加速度力，第二项是带尾迹记忆的环量力；这里半弦长为 b，转轴位置由 a 表示。

Kelvin 环量账本

Γb+∑jΓw,j=constant,ΔΓw=−ΔΓb\Gamma_b+\sum_j\Gamma_{w,j}=\mathrm{constant},\qquad \Delta\Gamma_w=-\Delta\Gamma_b

束缚环量每次变化，后缘就释放相反的增量；场景中的双色尾迹线就是这本离散账本。

常见难点

反向涡街不等于完整推进预测

容易误解尾迹看起来像推进涡街，因此效率读数就是现实飞行器的效率。

正确理解这里的流向力由升力投影和固定型面阻力代理构成；缺少黏性尾迹、三维翼尖损失和动态失速，所以效率明确带星号，只用于比较相位。

线性模型不会自动产生失速

容易误解把俯仰振幅推到最大，就能用同一个模型研究动态失速涡。

正确理解Theodorsen 理论假设小扰动和附着势流。峰值有效攻角超过界限时，Lab 会继续算以暴露外推趋势，但不会把它称为可信失速结果。

动手实验

  1. 01

    拆开两种运动

    先比较“纯升沉”和“纯俯仰”，观察尾迹环量、瞬时升力与平均流向力。

    应该观察到： 两种运动都能制造交替环量，但当前升力投影模型中，只有带升沉速度的工况直接形成平均推进贡献。
  2. 02

    组成推进相位

    选择“推进相位”，保持振幅不变，观察后缘包络、推力代理与功率代理。

    应该观察到： 俯仰抵消过大的运动攻角，同时让升力与升沉速度在周期中以有利方式配对，平均流向力变为正值。
  3. 03

    只反转相位

    切换到“反向相位”，不要改变频率和振幅。

    应该观察到： 尾迹拓扑和力的相位一起反转，平均流向力从推进转为阻力，说明相位不是装饰参数。