泡沫夹层梁在黏性水流下的稳定性分析：非线性流固耦合模型的构建与应用
本研究由伊朗塔比阿特莫达勒斯大学机械工程系的Ehsan Taati*（通讯作者）主导完成，并于2025年9月在线发表于期刊《Journal of Ocean Engineering and Science》。该研究聚焦于海洋工程领域的关键结构问题，旨在探究在黏性水流环境下，泡沫金属夹层梁的屈曲与后屈曲力学行为。
一、 研究背景与目标
夹层结构（由两层面板和轻质芯层构成）因其优异的抗冲击性、高比刚度及卓越的承载能力，在民用、海洋及近海工程中展现出巨大应用潜力，例如用于浮动平台、船体和海底设施。其中，泡沫金属作为一种典型的胞状结构材料，具有轻质、良好能量吸收和低热导率等特性，成为高性能夹层结构的理想芯材。尽管已有大量研究关注泡沫夹层梁在静态或空气中的力学行为（如屈曲、自由振动和非线性振动），但实际工程中的梁结构常处于流体环境中，流体-结构相互作用（Fluid-Structure Interaction, FSI）对其稳定性的影响不容忽视。
现有文献中，关于均质梁或夹层梁在流致振动方面的研究较多，但针对流体（特别是黏性流体）环绕下梁的屈曲与后屈曲问题的研究非常稀少。更重要的是，尚未有研究系统地考察周围黏性水流对泡沫夹层梁稳定性行为的影响。此外，现有的线性FSI模型由于忽略了接触表面的几何非线性效应，可能不足以准确预测结构的屈曲行为。
因此，本研究旨在填补这一空白。其核心目标是：首次建立一个非线性水弹性流固耦合（hydroelastic FSI）模型，用于分析黏性水流中泡沫夹层梁的屈曲与后屈曲行为。具体研究目标包括：1）推导考虑中等大变形的非线性流体压力分布表达式；2）基于等效单层理论和Timoshenko梁理论，结合von Kármán几何非线性，建立梁的力学控制方程；3）通过Galerkin方法求解，得到静态平衡路径、屈曲临界载荷和临界上游流速的闭合形式解；4）系统研究几何参数、边界条件、孔隙率分布及流动参数对夹层梁稳定性的影响。
二、 详细研究流程
本研究主要包含三个核心且环环相扣的理论建模与分析流程：1）泡沫夹层梁的几何非线性力学建模；2）无界黏性水流的非线性水动力压力建模；3）基于Galerkin方法的屈曲与后屈曲求解与分析。
第一流程：泡沫夹层梁的几何非线性力学建模 研究的对象为一个长L、宽b、总厚h的泡沫金属夹层梁，该梁由两片相同的正交各向异性面板和一片泡沫金属芯层对称粘合而成，芯层孔隙率沿厚度方向呈余弦分布。梁受到轴向压缩载荷Na和由黏性水流引起的横向分布载荷q(x)的作用。 研究采用等效单层（Equivalent Single-Layer, ESL）Timoshenko梁理论进行建模，同时引入了von Kármán几何非线性以考虑结构在屈曲前后可能发生的中等大变形。位移场由轴向位移u(x)、横截面转角φ(x)和横向挠度w(x)描述。通过应变-位移关系（包含了非线性项(1⁄2)(dw/dx)²）和最小总势能原理，推导出梁的控制微分方程组。 梁的本构关系考虑了面板和泡沫芯的不同材料属性。对于泡沫芯，其杨氏模量Ec和剪切模量Gc随厚度z变化，遵循公式：Ec = E_ref [1 - e0 cos(πz/tc)]，其中e0为无量纲孔隙率系数，tc为芯层厚度，E_ref为基体固体模量。材料属性的对称分布使得拉伸-弯曲耦合系数B11为零。 最终，通过数学推导，将控制方程化简为一个关于横向挠度w(x)的四阶非线性微分方程，这是进行屈曲与后屈曲分析的基础方程。同时，轴向位移u(x)也可以由挠度w(x)的平方梯度（即(dw/dx)²）的积分表达，体现了非线性变形导致的轴力与弯曲的耦合。
第二流程：无界黏性水流的非线性水动力压力建模 这是本研究的核心创新点之一。为了精确描述流体对结构的反馈作用，研究首次推导了无界、充分发展层流、不可压缩黏性水流沿梁长度方向的非线性压力分布的闭合形式表达式。 建模基于势流理论的基本思想，但针对黏性流进行了关键扩展。首先，假设由梁运动扰动的流体垂向速度分量vz随离梁表面距离呈指数衰减，即vz = Φ(x) e^{αz}，其中α为衰减率，其物理意义是受结构运动影响的流体体积，需通过实验确定。 最关键的一步是引入了非线性不可渗透条件。传统线性FSI模型假设流体速度垂向分量等于结构表面法向速度的线性近似。而本研究采用了Taati等人此前发展的关系式，在梁与流体的接触面（z = -h/2）上，考虑挠度梯度的影响：vz = U_∞ (dw/dx) [1 - 1⁄2 (dw/dx)²]。这个条件将流体运动与梁的非线性变形直接耦合起来。 随后，结合不可压缩流体的连续性方程和Navier-Stokes方程（N-S方程），进行一系列数学推导。利用近梁前缘处流速均匀（vx ≈ U_∞）的边界条件，以及从梁表面（z = -h/2）到无穷远处（z → -∞）压力积分（假设无穷远处压力为p_∞），最终得到了流体压力pf沿梁分布的显式表达式。该表达式由线性项(pl_f)和非线性项(pnl_f1, pnl_f2)组成，具体形式复杂，包含了挠度w(x)及其一阶、二阶、三阶导数与流体参数（密度ρ_∞、动力粘度μ、上游流速U_∞、衰减率α）的组合。这一推导过程成功地将非线性几何效应（通过不可渗透条件）和流体黏性（通过N-S方程）同时纳入了水动力载荷模型。
第三流程：基于Galerkin方法的屈曲与后屈曲求解与分析 此流程将前两个流程的成果结合，并进行求解。将第二流程得到的非线性压力分布q(x)=b*pf(x)代入第一流程得到的梁控制方程中，得到一个包含复杂非线性项的、以w(x)为未知函数的四阶微分方程。 为了求解这个方程并获得静态平衡路径（即载荷-挠度关系），研究采用了半解析的Galerkin方法。该方法假设横向挠度w(x)可以表示为一系列已知形状函数wm(x)的线性叠加，即w(x)= Σ w_m * w_m(x)，其中w_m为待求的未知系数。形状函数wm(x)根据梁的边界条件选取为欧拉-伯努利梁的屈曲模态形状，并针对简支（SS）、两端固支（CC）和悬臂（CF）三种典型边界条件给出了具体的表达式和特征根λ_m。 将假设的w(x)表达式代入控制方程，利用Galerkin法的加权残差原理（即要求残差与各形状函数的内积为零），将复杂的偏微分方程转化为一组关于未知系数w_m的非线性代数方程。经过详细的积分运算，最终得到了一个关于单个模态系数w_m（以第一模态m=1为例）的闭合形式方程。该方程清晰地展示了轴向载荷Na、上游流速U_∞、静水压力p_∞与模态挠度系数w_m之间的非线性关系，即静态平衡路径的数学描述。 从这个方程出发，可以通过两种方式获取临界值：1）忽略所有非线性项，得到一个线性特征值问题，从而解析地求出屈曲临界压缩载荷Na,cr与临界上游流速U_∞,cr的表达式，它们与衰减率α、材料刚度、边界条件等参数直接相关。2）直接求解完整的非线性方程，通过追踪载荷参数变化时挠度解的分岔点，来确定屈曲临界状态，并进一步获得屈曲后的平衡路径（后屈曲行为）。
三、 主要研究结果
本研究通过理论推导和数值算例，得到了一系列揭示黏性水流下泡沫夹层梁稳定性规律的重要结果。
1. 黏性水流显著降低梁的稳定性承载力。 数值结果表明，在给定的水流条件下（U_∞ = 50 m/s，α = 20 m^{-1}），与干燥（无流）情况相比，简支泡沫夹层梁的屈曲临界压缩载荷降低了约25%。通过与商业有限元软件ANSYS的模拟结果对比，验证了本解析模型的有效性，两者差异在可接受范围内（对于不同长细比梁，差异在1.3%至8.7%之间），且趋势一致。
2. 临界参数对流体衰减率高度敏感。 衰减率α是模型中的一个关键参数，表征了结构运动对流体扰动的影响范围。研究显示，临界上游流速U_∞,cr随α的增大而显著变化。例如，当α从10 m^{-1}增大到30 m^{-1}时，对于简支、固支和悬臂梁，U_∞,cr的最大变化幅度可达约73%。这强调了对衰减率进行精确确定对于评估流体中结构稳定性的至关重要。
3. 孔隙率和芯层厚度的影响。 参数研究发现，水流对夹层梁屈曲行为的影响程度，与泡沫芯层的属性密切相关。对于高孔隙率（e0大）和较薄芯层（tc小）的梁，水流引起的临界载荷降幅更为显著。例如，当e0增大时，干/湿工况临界载荷的百分比差δNa,cr(%)增加；而当芯层厚度增至三倍时，δNa,cr(%)可减小约12.8倍。这说明，芯层的刚度贡献越弱，流体耦合效应的影响就越突出。
4. 几何非线性与流体非线性效应的必要性。 研究通过对比分析，得出了两个关于模型选择的重要结论： * 线性FSI模型的不足：如果将压力分布中的非线性项全部忽略（即使用线性压力），虽然在屈曲前路径上与非线性模型结果吻合较好，但无法准确预测屈曲分岔点和后屈曲平衡路径。如图8所示，在线性压力模型下预测的后屈曲行为与非线性模型存在明显差异。 * 上游静水压力的影响：线性屈曲分析（公式35b）无法体现上游静水压力p_∞对临界流速的影响。然而，非线性分析（图7b）表明，p_∞的增大会导致屈曲前挠度显著增加，从而使失稳在更低的U_∞下发生。这进一步证实了考虑几何非线性对于准确分析流载梁稳定性是必需的。
5. 复杂载荷下的失稳模式。 对于已发生屈曲的梁（U_∞超过临界值），进一步增加上游静水压力p_∞会引发突弹跳变（snap-through）行为，即结构从一个平衡位形突然跳跃到另一个稳定的、挠度更大的向下位形。这揭示了在流体与轴向载荷共同作用下，夹层梁可能表现出丰富的非线性失稳现象。
四、 研究结论与意义
本研究成功建立了一个用于分析黏性水流中泡沫夹层梁屈曲与后屈曲行为的非线性水弹性流固耦合模型。其主要结论是：黏性水流会显著降低泡沫夹层梁的稳定性承载力；结构的临界参数（如屈曲载荷、临界流速）对流体衰减率极为敏感；泡沫芯层的高孔隙率和较小厚度会放大水流的不利影响；而线性FSI模型由于忽略了接触表面的几何非线性以及压力分布的非线性部分，不足以准确预测此类结构的稳定性行为，尤其是在后屈曲 regime 和考虑静水压力效应时。
该研究的科学价值在于：1）首次系统地将黏性流体力学与具有几何非线性的夹层结构力学相结合，为流固耦合稳定性分析提供了一个创新的理论框架。2）推导出的非线性水动力压力闭合解是该领域的理论贡献，为后续相关研究奠定了基础。3）采用的等效单层梁理论与Galerkin半解析法相结合的策略，在保证一定精度的同时，提供了极高的计算效率，非常适合于海洋结构的概念设计和参数化研究。
其工程应用价值直接面向海洋工程领域：为在复杂流体环境中（如海水绕流）工作的泡沫夹层梁、板、壳等轻质高性能结构的设计与安全性评估，提供了关键的理论依据和分析工具。设计人员可以利用本研究给出的闭合解或方法，快速评估不同设计参数（材料、几何、边界条件）和流动条件对结构稳定性的影响，从而优化设计方案，避免流致失稳事故。
五、 研究亮点与创新
六、 其他有价值内容
研究附录中详细列出了Galerkin方法求解过程中产生的所有系数（a1至a9）的积分定义式，确保了研究过程的透明性和结果的可复现性。此外，引言部分对等效单层（ESL）理论与层合（Layerwise）理论进行了简要对比和评述，为读者理解本研究建模选择的背景和适用范围提供了有用的信息。文末引用了大量相关文献，涵盖了夹层结构力学、泡沫金属材料、流固耦合振动等多个方面，为感兴趣的读者提供了深入探索的路线图。