⚡ electrical engineering

Time-Domain Two-Way Fluid–Structure Interaction Analysis of the HIRENASD Transonic Wing: Static Aeroelastic Response and Flutter Prediction

本研究通过针对 NASA 基准测试，验证了 ANSYS Workbench 中一种时域、松耦合的流固耦合框架，能够准确预测 HIRENASD 跨声速机翼的静态气动弹性响应和颤振边界，证明了其作为一种可复现且具备工业易用性的工具，在飞机机翼初步评估中的可行性。

原作者： Siddalingappa P Kodigaddi

发布于 2026-07-13

📖 1 分钟阅读☕ 轻松阅读

CC BY 4.0

原作者： Siddalingappa P Kodigaddi

原始论文采用 CC BY 4.0 许可（https://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象一下，一只巨大的、富有弹性的机翼在天空中翱翔，与风共舞。现在，想象这阵风不仅仅是吹过机翼；它实际上是在“推”和“拉”着机翼，而机翼也在随之弯曲和扭转，从而改变了气流的方式。这种看不见的舞蹈被称为流固耦合（Fluid–Structure Interaction, FSI）。如果这种现象失控，机翼就会开始剧烈抖动，直到断裂。这种可怕的抖动被称为颤振（flutter）。

在这项研究中，研究员 Siddalingappa P. Kodigalli 在计算机内部构建了一个超级精确的虚拟风洞，用以观察这种舞蹈如何在一种名为 HIRENASD 的特定机翼上发生。这可不是普通的机翼；它是一种为高速跨声速喷气机（飞行速度略低于音速）设计的尖端“超临界”机翼。该研究的目标是观察标准的商业化计算机程序（ANSYS）是否能够精确预测该机翼何时开始发生颤振，而无需使用极其昂贵的定制代码。

虚拟风洞测试

首先，研究员必须证明他的虚拟风洞是真实的。他设置了一个速度为 Mach 0.8（约 614 英里/小时）且空气密度极高（雷诺数为 7×10⁶）的模拟环境。他将计算机的结果与 NASA 的一个著名基准进行了对比。

结果： 计算机对升力、阻力和扭转力的预测极其精确，与 NASA 的数据误差不到 4%。这就像是在猜一只猫的重量，误差仅有几克。
机翼的自然节奏： 在风吹起之前，研究员检查了机翼在自身状态下的振动情况（就像拨动吉他弦）。现实中，机翼的主弯曲频率测量值为 26.0 Hz。计算机预测的结果为 26.504 Hz。这仅有 1.94% 的差异，证明了虚拟机翼在刚度和重量上与真实机翼高度一致。

静态拉伸（在抖动发生前）

接下来，他们通过改变不同的角度（从 -1.5° 到 +4.5°）来观察机翼在稳定压力下的弯曲程度。

弯曲： 在最高角度（4.5°）和速度（Mach 0.8）下，机翼尖端向下弯曲了 12 mm。
扭转： 由于机翼是后掠式的（像回力镖一样向后倾斜），这种弯曲实际上使机翼产生了轻微的扭转，这有助于减轻根部的载荷。这就像一根柔韧的树枝在风中弯曲以避免折断一样。计算机生成的这些弯曲数值与其它高端模拟结果几乎完全吻合。

重头戏：寻找颤振速度

这是最关键的时刻。研究员关闭了“阻尼”（通常会阻止抖动的内部摩擦），并逐渐增加模拟中的风速，以观察机翼何时会变得疯狂。他们测试了两种不同的“大气环境”：

高空（稀薄空气）： 速度为 Mach 0.8，密度为 0.2525 kg/m³。
低空（浓厚空气）： 速度为 Mach 0.7，密度为 0.50 kg/m³。

研究发现：

在 Mach 0.8（稀薄空气）下： 机翼一直保持平静，直到速度达到 266.56 m/s。就在这个精确的速度下，机翼开始以一种恒定的、稳定的节奏进行抖动（既没有增强也没有衰减）。这就是颤振速度。抖动发生的频率约为 22 Hz。
在 Mach 0.7（浓厚空气）下： 机翼一直保持平静，直到 230.46 m/s。此时，抖动频率约为 25.5 Hz。

“顿悟”时刻：
你可能会想：“等等，机翼在稀薄空气中的颤振速度（266.56 m/s）比在浓厚空气中（230.46 m/s）更高。这是否意味着在高空更安全？”
并不完全是这样。 论文解释说，虽然速度更高，但浓厚空气情况下的风力（动压）实际上要强得多。

低速（浓厚空气）下的“颤振压力”为 13.27 kPa。
高速（稀薄空气）下的“颤振压力”仅为 8.97 kPa。
这意味着，尽管在浓厚空气中移动速度较慢，但机翼承受的压力却增加了 47.9%。论文证实了一个经验法则：危险程度与密度乘以速度的平方成比例。因此，即使在爬升过程中保持相同速度，一个在低空安全的机翼也可能变得危险，因为空气变薄了，气动“刚度”下降了。

这是什么样的抖动？

研究员仔细观察了抖动过程，以了解其发生的机制。

单模态摆动： 机翼并没有以复杂的结状方式同时进行扭转和弯曲。相反，它是一种简单的上下弯曲运动（就像跳水板一样）。
原因： 机翼的扭转频率（279.13 Hz）远高于其弯曲频率（26.504 Hz）。它们之间的差距非常大，以至于永远不会发生“合并”（coalesce）来产生复杂的致命颤振。这是一种简单的、单模态的弯曲。
偏移： 在较高的速度（Mach 0.8）下，气压实际上让机翼感觉更“软”，使其抖动频率从 26.504 Hz 下降到了 22 Hz。在较低的速度（Mach 0.7）下，空气的“软化”作用较小，因此频率仅下降到 25.5 Hz。

这对未来意味着什么

这篇论文并不声称已经永久“解决”了所有飞机的颤振问题。相反，它表明商业化的现成软件（ANSYS Fluent 和 ANSYS Mechanical）可以用于高置信度地预测这些危险速度，而无需使用秘密的、政府专用的代码。

信心： 这些结果是基于经过现实世界数据和其他超级计算机基准严格校验的模拟。作者估计其方法在 5–8% 的误差范围内是准确的。
局限性： 模拟并未考虑风洞中的极端低温（这会使金属稍微变软）或复杂的气温变化，但其主要弯曲模态的准确性足以信任其颤振预测。

简而言之，这项研究构建了一个高科技机翼的数字孪生，证明了它的行为与真实情况一致，并利用它找到了机翼开始不受控制地“起舞”的精确速度。它表明，虽然在稀薄空气中机翼可以飞得更快而不发生抖动，但风的力量才是真正的危险所在，而这种新的计算机方法可以帮助工程师设计出更安全的机翼，造福未来。

技术摘要：HIRENASD 跨声速机翼时域双向流固耦合分析

问题陈述
运输机机翼在跨声速状态下的结构完整性受复杂的气动弹性耦合（尤其是颤振）控制，这可能导致灾难性的失效。线性频域方法（如双层涡格法 DLM）往往缺乏捕捉跨声速流中固有关键非线性现象（包括移动激波、激波诱导边界层分离以及强流固耦合）的保真度。虽然高雷诺数气动弹性动力学（HIRENASD）项目为超临界机翼（BAC3-11 翼型）提供了全面的实验数据库，但在使用商用通用型 CFD-CSD 求解器验证这些构型方面，公开文献中仍存在显著空白。此外，针对该基准模型，在多个马赫数下进行系统性的时域颤振边界识别，以及定量分析自由流密度和马赫数如何共同影响颤振动压，尚未得到充分解决。

方法论
本研究采用在 ANSYS Workbench 中实现的基于时间的、松耦合的双向流固耦合（F least-coupled, two-way FSI）框架。该方法将求解过程划分为：

CFD： 使用 ANSYS Fluent（基于密度的求解器），利用 Reynolds 平均 Navier-Stokes (RANS) 方程和 SST $k-\omega$ 湍流模型。流体域采用非结构化四面体网格（约 600,000 个单元），并包含棱柱层以实现 $y^+ \approx 1$ 。
CSD： 使用 ANSYS Mechanical，采用 HIRENASD 机翼结构（镍马氏体钢）的六面体有限元模型（41,923 个单元）。模态分析使用 Block-Lanczos 算法，瞬态积分采用 Newmark- $\beta$ 方法。
耦合： ANSYS System Coupling 通过通用网格接口（GGI）管理力守恒的数据交换，并通过智能桶（Smart Bucket）插值进行位移映射。网格变形通过弹簧模拟平滑处理及局部重构网格来完成。

本研究通过与 NASA FUN3D 空气动力学基准数据以及 HIRENASD 实验结构频率进行对比来验证该框架。研究在马赫数 ( $M$ ) 为 0.7 和 0.8、攻角从 $-1.5^\circ$ 到 $+4.5^\circ$ 的条件下进行静态气动弹性分析。动态颤振分析通过系统地改变自由流速度来识别颤振边界，旨在表征两种不同飞行条件下亚临界（阻尼）、临界（中性）和超临界（发散）响应机制的转变： $M=0.8$ ( $\rho = 0.2525$ kg/m³) 和 $M=0.7$ ( $\rho = 0.50$ kg/m³)。

核心贡献

框架验证： 本研究成功建立了一个仅使用商业工具（ANSYS Fluent, Mechanical, System Coupling）即可实现的、可复现且具有工业相关性的跨声速气动弹性 FSI 工作流。
静态气动弹性基准测试： 该框架重现了稳态空气动力学系数（ $C_L, C_D, C_M$ ），其误差在 NASA FUN3D 参考数据的 1.87% 至 3.52% 之间，并能以 1.94% 的误差预测基本弯曲频率（26.504 Hz）。
动态颤振特征分析： 研究提供了系统性的时域颤更边界确定方法，明确表征了从阻尼振荡到发散振荡的转变。研究表明，在两个测试的马赫数范围内，颤振均由单一模式的弯曲机制控制。
物理缩放洞察： 研究量化了马赫数、自由流密度与颤振动压之间的关系，证明了尽管在低密度（高高度）下颤振速度可能更高，但由于遵循缩放律 $q_f \propto \rho V_f^2$ ，颤振动压在较高密度下显著更高。

结果

静态响应： 在 $M=0.8$ 时，升力线斜率约为每度 0.107。在 $4.5^\circ$ 攻角下，最大翼尖挠度达到 12 mm（无量纲值 $\delta/b = 0.0093$ ）。 $M=0.8$ 下的阻力极曲线表现出剧烈的激波诱导阻力上升，而 $M=0.7$ 的极曲线则更接近抛物线形状，并由诱导阻力主导。
模态分析： 第一阶面外弯曲模态预测为 26.504 Hz（与实验值的误差为 1.94%）。高阶模态呈现出系统性的高估（9–13%），这归因于有限元模型中的结构简化（例如缺少紧固件的柔性）。
颤振边界：
- 在 $M=0.8$ ( $\rho = 0.2525$ kg/m³) 条件下： 颤振速度 ( $V_f$ ) 为 266.56 m/s，颤振动压 ( $q_f$ ) 为 8.97 kPa。颤振频率约为 22 Hz，由于跨声速气动刚度修正，较真空状态下的结构频率降低了 16.9%。
- 在 $M=0.7$ ( $\rho = 0.50$ kg/m³) 条件下： 颤振速度为 230.46 m/s，颤振动压为 13.27 kPa。颤振频率约为 25.5 Hz，显示出较小的频率降幅（3.8%），表明跨声速非线性效应较弱。
缩放律： 尽管 $M=0.7$ 时的颤振速度较低，但由于自由流密度高出 98%，其颤振动压比 $M=0.8$ 时高出 47.9%。这证实了维度缩放关系 $q_f \propto \rho V_f^2$ 。

意义与主张
论文主张，经过验证的 ANSYS Workbench 框架为运输类飞机机翼的初步颤振边界评估提供了一条实用的、具有工业相关性的路径，无需依赖定制的内部代码或专有求解器。研究强调，仅使用商业工具使得该流程能够立即转移到工业设计环境中。

作者断言，对颤振动压随密度变化的定量确认，为将风洞模型结果外推至全尺寸飞行条件（特别是关于高度相关的颤振风险）提供了科学严谨的基础。研究指出，即使在恒定马赫数下，由于有效气动刚度的降低，机翼在低高度（高密度）下可能处于无颤振状态，但在巡航高度（低密度）下可能会接近其颤振边界。研究结论认为，基于与 AGARD 445.6 机翼已验证框架的对比，所确定的颤振速度具有约 5–8% 的方法论不确定度，且观察到的单模态弯曲机制与该特定构型中弯曲与扭转模态之间较大的频率间隔是一致的。

虚拟风洞测试

静态拉伸（在抖动发生前）

重头戏：寻找颤振速度

这是什么样的抖动？

这对未来意味着什么

类似论文