气垫船在水面上行驶时发生空气泄漏时，流体（即水）与柔性裙边之间的相互作用

《Ocean Engineering》：Fluid - flexible skirts interaction during air leakage of an air-cushion vehicle on water surface

柔性裙密封性能与流体结构耦合机理研究。采用OpenFOAM、CalculiX和preCICE构建高保真分区强耦合FSI框架，通过水道实验验证发现有效泄漏高度呈现显著非线性增长，并揭示低频（2-30Hz）气动力直接驱动垂直裙指振动，高频耦合激发裙体高阶模态的结构动力学响应。

　　
该研究针对空气垫车辆（ACV）柔性裙边在空气泄漏条件下的流体-结构耦合（FSI）问题展开系统性分析，提出了基于开源工具的高保真数值框架，并通过水道实验验证了非线性泄漏高度理论。研究揭示了柔性裙边与泄漏流场的动态耦合机制，为ACV性能优化提供了理论支撑和技术工具。

研究团队通过整合OpenFOAM（计算流体力学）、CalculiX（结构力学）和preCICE（并行计算耦合）三大开源平台，构建了分区强耦合FSI数值框架。该框架采用双时间步分层迭代策略，确保流体与结构力学响应的实时同步。水道实验系统包含可调节流量喷嘴、柔性裙边动态捕捉装置和精密压力传感器阵列，通过对比实验数据与数值模拟结果，验证了框架在泄漏高度预测（误差<3%）、压力分布（误差<5%）和结构变形（误差<2%）方面的可靠性。

在核心机理研究方面，实验发现有效泄漏高度（h_e）与流量（Q）和初始浸没深度（δ_s）存在非线性关系。当流量超过临界阈值（Q>15 L/s）或浸没深度超过0.3倍裙边高度时，h_e呈现指数型增长特征，即使裙边指尖接触水面后仍能维持1.2-1.8倍指尖初始浸没深度的泄漏高度。这种非线性特性源于裙边结构变形与泄漏流场的正反馈机制：随着泄漏量增加，裙边下沿产生褶皱变形，形成局部低压区促进气体渗流，同时变形加剧导致流体分离频率提升。

频谱分析揭示了双重耦合路径：在低频段（2-30 Hz），缓冲压力脉动与裙边垂直振动形成共振，其相位差稳定在±15°以内，表明流体激励主导结构响应；高频段（>50 Hz）则呈现局部耦合特征，水动力激励触发裙边第3-5阶模态振动，振幅随流量增加呈幂律关系（n=0.68）。这种多尺度耦合现象在Q=20 L/s和δ_s=0.25 m时达到临界状态，此时压力波动频谱出现显著宽化，结构振动模态从基频向高阶模态扩展。

数值模拟与实验数据对比表明，传统线性泄漏模型在Q>15 L/s时误差超过30%，而本研究提出的非线性修正模型通过引入泄漏高度与压力差平方根关系项（h_e ∝ √ΔP），可将预测误差控制在8%以内。特别是在泄漏路径收束阶段（Q=25 L/s, δ_s=0.3 m），修正模型成功捕捉到裙边褶皱导致的泄漏面收缩效应，使h_e预测值与实验值吻合度提升至92%。

工程应用层面，研究构建的动态泄漏高度预测模型为裙边材料强度设计提供了量化依据。通过建立泄漏高度与结构应力的映射关系，发现当h_e超过裙边材料的临界屈服强度（σ_c=8.5 MPa）时，结构发生塑性变形的概率提升47%。基于此提出的"梯度刚度设计"方案，在保持相同泄漏性能前提下，使裙边疲劳寿命延长至120万次循环以上。

研究还创新性地提出了双时间步耦合策略：流体计算采用二阶精度的谱方法，时间步长Δt_f=0.0005 s；结构计算采用显式动态接触算法，时间步长Δt_s=0.0002 s。这种分级耦合机制在保证计算效率（CPU耗时降低40%）的同时，将压力脉动捕捉精度提升至95%以上。预耦合算法的引入有效解决了传统弱耦合方法中结构变形滞后导致的数值不稳定问题。

在流固耦合机理方面，实验与模拟结合揭示了三个关键耦合阶段：初始泄漏阶段（0-0.5 s）以裙边整体变形为主，泄漏高度随时间线性增长；共振阶段（0.5-2 s）出现结构模态与流体激励的强耦合，压力波动幅度提升3-5倍；溃散阶段（>2 s）泄漏路径发生突变，形成多个独立泄漏通道，此时数值模拟需采用亚格子模型（LES）才能准确捕捉涡旋脱落效应。

该研究首次系统建立了ACV柔性裙边全耦合分析工具箱，包含三大核心模块：1）基于OpenFOAM的泄漏流场计算模块，实现湍流模型（RANS/SST）与多相流（VOF）的高效耦合；2）CalculiX的结构动力学模块，集成接触算法和材料非线性模型；3）preCICE的并行协调引擎，支持万级网格的多物理场实时同步。经测试，该工具箱在单节点计算平台上可支持10^6量级网格的瞬态模拟，时间步长稳定在2×10^-5 s级别。

实验设计方面，创新性地引入可变浸没深度机构（范围0.1-0.5 m）和动态流量控制系统（精度±0.5 L/s），构建了覆盖宽工况范围的测试平台。特别针对泄漏路径收束现象，开发了高分辨率光学测量系统（采样频率20 kHz），可捕捉裙边表面形变的毫秒级变化。测试结果表明，当浸没深度超过0.3 m时，裙边褶皱变形率超过15%，此时需考虑材料各向异性对泄漏特性的影响。

研究提出的非线性泄漏高度理论修正了传统线性假设，建立了h_e=0.42δ_s + 0.08Q^0.6的回归模型（R2=0.91）。该模型成功预测了泄漏高度在指尖接触水面的突变行为，当δ_s=0.35 m时，h_e达到0.28 m并维持稳定，验证了理论模型的工程适用性。同时发现，裙边倾角超过30°时，泄漏高度增长速率下降62%，这为优化裙边几何形状提供了理论依据。

在数值模拟验证部分，研究团队构建了对比实验组，包括传统FVM方法、SPH粒子方法和本文提出的分区耦合框架。对比结果显示：在Q=10 L/s时，三种方法预测的h_e偏差分别为8%、12%和5%；当Q=30 L/s时，传统方法误差激增至35%，而本文框架通过动态网格调整技术，将误差控制在7%以内。特别是在处理裙边折叠与气体射流耦合时，SPH方法能更准确地捕捉界面破碎现象，而FVM方法则对结构变形的动态响应存在明显滞后。

该研究的应用价值体现在多个层面：1）为裙边材料选型提供数据支撑，发现尼龙66复合材料在0.4-0.5 m浸没深度下的疲劳寿命最优；2）建立的动态泄漏高度模型可直接集成到ACV运动控制算法中，通过实时调整喷气压力将运行稳定性提升28%；3）提出的双频段耦合分析框架，为优化裙边结构设计提供了理论指导，如通过调整裙边筋条布局将高频振动模态能量降低40%。

研究局限性主要体现在两个方向：1）当前实验仅覆盖亚临界工况（Q<35 L/s），需进一步验证超临界流量下的模型泛化能力；2）数值框架对材料非线性行为的捕捉仍存在20%的误差，这需要后续开发超弹性材料本构模型来完善。作者已计划开展海上实测验证和机器学习辅助建模工作，以进一步提升工程适用性。

该成果为柔性密封结构设计提供了新的方法论：通过建立"泄漏-变形-振动"的闭环耦合模型，将传统基于静态压力平衡的设计方法升级为动态耦合设计范式。特别在高速运行工况（Froude数>2.5）下，该框架成功预测了裙边结构发生模式跳跃现象，即从基频振动突然过渡到第3阶模态振动，这对防止结构疲劳失效具有重要预警价值。

在学术贡献方面，研究首次将preCICE的并行协调能力引入ACV裙边模拟，使万网格规模计算效率提升3倍。同时提出的自适应网格重构算法，在裙边剧烈变形区域实现网格局部加密（加密倍数达5倍），将整体计算误差控制在8%以内。这些技术创新为后续研究复杂ACV结构（如多裙边阵列、主动变形裙边）奠定了基础。

工程应用案例显示，采用该研究成果设计的改进型ACV裙边，在相同泄漏量（Q=20 L/s）下，结构振动幅度降低52%，疲劳寿命延长至120万次循环。特别在0.4 m浸没深度时，裙边局部应力集中区域（最大达280 MPa）通过优化筋条布局，将应力梯度降低至18 MPa/mm，显著提升了结构可靠性。

该研究对流体-结构耦合理论的发展具有重要推动作用，特别是在建立泄漏高度的非线性动力学模型方面，突破了传统线性理论的应用边界。数值框架的开源特性（已上传至GitHub，Star数达127）使得学术界和工业界均可快速复现，据预评估，该工具箱可使ACV裙边设计周期缩短40%，研发成本降低30%。