本文的主要作者为 Miguel R. Visbal 和 Caleb J. Barnes,分别隶属于 Air Force Research Laboratory(美国空军研究实验室),研究地点位于 Wright-Patterson 空军基地,俄亥俄州,美国。本研究发表于 AIAA SciTech 2025 Forum(AIAA 航空科学与技术论坛)上。该研究以题为“Effect of Impinging Shock Strength on the Response of a Flexible Panel”,主要探讨了冲击激波强度对柔性面板动态响应及稳定性的影响。
本研究的科学领域为流体-结构相互作用(fluid-structure interaction)及面板气动弹性稳定性,重点考察柔性面板在低超音速条件下受到冲击激波作用时的复杂动态响应。从学术背景来看,柔性面板的流动—声学—结构相互作用是关乎结构完整性、疲劳寿命和声辐射的重要课题。
面板抖振(panel flutter)作为一种著名的现象,近年来在不同马赫数范围内的稳定性问题受到广泛研究。Dowell 和 Mei 等专家对该领域的潜在流动理论进行了重要的综述性讨论。此外还有研究者发展出了求解非线性 Euler 方程的数值模型以研究无穷大长宽比条件下的面板行为。对于柔性面板与激波的作用,已有研究表明激波的强度显著影响面板的自激振荡行为,但此前的大多研究仅探讨了较为简单或者标准激波条件下的面板响应。
为填补这一空白,本文旨在进一步描述冲击激波强度对柔性面板稳定性影响的全貌,尤其是在涉及规则反射和马赫反射(Mach reflection)情形下的动态演化过程。此外,研究还探讨了冲击激波的作用位置和面板下腔体压力对其响应特性的影响。
流体-结构耦合方法
本文使用了经验证的耦合求解策略,通过将二维不可压 Euler 方程与 von Karman 型非线性板模型相结合,建立了适用于大幅板面变形的求解框架。面板变形被映射至可变形的流场网格中,采用子迭代(subiteration)方式实现气动力和结构动力学方程的隐式耦合,从而避免伪动力学结果的产生。
求解流程与数学方程
1. 流场方程
选用不可压二维 Euler 方程模拟动态变形表面附近的复杂流动。这些方程被引入强守恒格式(strong conservative form)并在非定常曲线坐标系中重写。 2. 结构动力学方程
柔性板的变形由非线性 von Karman 板模型控制。该模型假设板面材料为各向同性且厚度均匀,并允许大幅线性和非线性位移耦合作用。
网格生成与动态调整
网格采用拉伸型直角坐标系,沿面板和远场均布置密集点,动态网格通过代数方法根据面板形变重构,确保在变形过程中网格坐标保持精确。
研究主要分为以下实验情景: 1. 标准面板抖振(无激波施加)。 2. 不同强度下的规则斜激波反射(regular oblique shock reflection)。 3. 规则反射发展为非规则马赫反射的动力学过程。 4. 冲击激波作用位置变化的影响。 5. 面板下腔体压力变化的影响。
实验细节:
所有的计算假定马赫数 (M_\infty = 2.0),面板参数包括厚度比 (h/a\ =0.002)、质量比 (\mu_s=0.1)。通过多个仿真实验探索动态压力 (\lambda)、激波强度(由压力比 (p_3/p_1) 表征)以及其他参数对面板稳定性的作用。
面板在不同激波压力比下的动态响应
1. 无激波条件(标准面板抖振)
当 (p_3/p_1 = 1.0)(无激波)时,激振在动态压力 (\lambda\ > 625) 时开始出现,抖振以超临界分岔的形式发生,主要展现出了标准面板抖振下的一阶模态振型,且频率保持恒定。
弱激波(例如 (p_3/p_1 = 1.2))
激波引入后,面板停止抖振,变形呈现为稳定的二阶模态。这表明弱激波引发的静态压力差能够提升面板刚度,延迟抖振的发生。
中等激波(例如 (p_3/p_1 = 1.4))
此情况下出现了复合的面板响应:当动态压力超过临界值后,面板存在两种稳定状态:
强激波(例如 (p_3/p_1 = 1.8) 或 (2.2))
在更高压力比情形中,面板振荡以超临界分岔形式发生,而临界动态压力较低且振幅增加。
极强激波及马赫反射
当 (p_3/p_1 \geq 3.1) 时,面板停止抖振,进入发散(divergence)模式。马赫反射(Mach Stem)的形成使气动加载产生显著变化,但未导致新的动态稳定模式。
本文通过理论和数值模拟,首次系统性地探讨了冲击激波强度对柔性面板动态响应的作用规律,揭示了以下重要结论: 1. 激波强度对面板动态响应的显著影响,响应从标准抖振转变为发散、LCO,最终重新进入发散状态。 2. 激波强度增加可导致复杂的动力学分岔行为,包括亚临界至超临界分岔的切换过程。 3. 在面板未发生分离涡作用时,面板形变对马赫反射的结构变形具有显著影响,为实际工程提供了潜在优化路径。
本研究在科学价值和工程应用两方面提供了关键见解。理论上,它拓展了面板气动弹性问题的稳定性理论基础;实际中,这些结果对高超音速飞行器设计、航空航天器材料优化的关键性结构设计具有指导意义。
研究还探讨了冲击激波的作用位置和腔体压力对动态响应的影响,表明这些参数尽管未引发 LCO,却会对马赫反射的显示模式产生显著作用,为未来更复杂模型和非均匀压力场条件下的研究铺垫了基础。