基于SPH法的舱内液体晃荡数值模拟：理论、应用与优化.docxVIP

基于SPH法的舱内液体晃荡数值模拟：理论、应用与优化

一、引言

1.1研究背景与意义

在现代工程领域中，舱内液体晃荡问题广泛存在于航空航天、船舶、汽车运输等多个关键行业，其对系统的安全性、稳定性以及性能表现有着至关重要的影响，因而受到了学术界和工程界的高度关注。

在航空航天领域，飞行器在飞行过程中会经历各种复杂的飞行姿态变化，如加速、减速、转弯、颠簸等，这些飞行状态的改变会导致飞行器舱内的液体，如燃油、推进剂等产生晃荡现象。液体晃荡不仅会消耗飞行器的能量，降低其飞行效率，还可能引发一系列严重的问题。例如，液体晃荡产生的惯性力和冲击力可能会对飞行器的结构造成额外的载荷，长期作用下可能导致结构疲劳损坏，影响飞行器的使用寿命和可靠性。在极端情况下，液体晃荡还可能干扰飞行器的飞行姿态控制，使飞行器偏离预定的飞行轨道，甚至引发飞行事故，严重威胁飞行安全。

在船舶领域，液舱晃荡同样是一个不容忽视的问题。船舶在海上航行时，会受到风浪、水流等多种外部因素的作用，导致船体发生摇晃和颠簸。这些船体运动使得液舱内的液体产生晃荡，当晃荡较为剧烈时，会产生巨大的冲击力作用于液舱壁，可能导致液舱结构损坏，影响船舶的航行安全。历史上曾发生过多起因液舱晃荡导致船体结构受损甚至船舶沉没的事故，给生命财产带来了巨大损失。此外，液舱晃荡还会影响船舶的稳性和操纵性，增加船舶在航行过程中的风险。

为了深入研究舱内液体晃荡问题，数值模拟方法成为了一种重要的手段。在众多数值模拟方法中，光滑粒子流体动力学（SPH）方法凭借其独特的优势脱颖而出。SPH方法是一种无网格的拉格朗日粒子方法，它将连续介质离散为一系列具有质量、速度、密度等物理属性的粒子，通过粒子之间的相互作用来模拟流体的运动。与传统的网格方法相比，SPH方法具有诸多显著优点。首先，SPH方法无需预先划分网格，避免了网格生成过程中的复杂性和局限性，尤其适用于处理具有复杂边界形状和大变形的流体问题。在舱内液体晃荡问题中，液体自由表面在晃荡过程中会发生大变形，传统网格方法容易出现网格扭曲和畸变，导致计算精度下降甚至计算失败，而SPH方法能够自然地适应这种大变形，准确地捕捉自由液面的运动。其次，SPH方法在处理多相流问题时具有天然的优势。在实际工程中，舱内液体往往与周围的气体或其他介质共存，形成多相流系统，SPH方法可以方便地处理不同相之间的相互作用和界面追踪。此外，SPH方法还具有并行计算效率高的特点，能够充分利用现代计算机的多核处理器和分布式计算平台，大大缩短计算时间，提高计算效率。

由于SPH方法的这些优势，其在解决实际工程问题中具有极高的价值。通过基于SPH法的数值模拟，可以深入了解舱内液体晃荡的复杂流动特性，如液体的速度分布、压力分布、自由液面的波动规律等。这些模拟结果为飞行器、船舶等的设计和优化提供了重要的理论依据。设计人员可以根据模拟结果，合理调整液舱的形状、尺寸、内部结构以及装载方式，以减小液体晃荡的影响，提高系统的安全性和性能。此外，SPH方法还可以用于预测不同工况下舱内液体晃荡的情况，为工程实践中的操作和决策提供参考，有助于提前采取有效的措施来避免潜在的风险。

1.2国内外研究现状

自SPH方法提出以来，国内外学者针对其在舱内液体晃荡模拟方面开展了大量的研究工作，取得了一系列有价值的成果。

在国外，早期的研究主要集中在SPH方法的理论基础完善和算法改进上。Lucy和Gingold首次将SPH方法用于解决天体物理学问题，为该方法的发展奠定了基础。随后，Monaghan将SPH方法应用于任意三维可压缩流体运动的模拟，并进一步拓展到具有固体边界的不可压流体运动，推动了SPH方法在流体力学领域的应用。在舱内液体晃荡模拟方面，Colagrossi等通过SPH数值模拟和实验相结合的方式，研究了接近临界充液深度情况下流体晃动作用及在抨击作用下的二维流体晃动，分析了罐内充液高度、谐波激励的振幅和周期对晃动现象相关特征的影响，其研究结果得到了SPH理论分析和实验结果比较的支持。Landrini对二维矩形液舱的横荡运动进行了模拟，给出了波面图及波高随时间变化的曲线，直观地展示了液体晃荡的过程。Iglesias等采用SPH方法数值模拟了水槽中渔船的被动滚动阻尼器，研究了水槽的抗扭和晃动类型问题，考虑了不同形状水槽及三维效果的影响，并通过实验定量研究了相位滞后情况下流体的运动，验证了数值计算分析结果。

国内学者在SPH法模拟舱内液体晃荡方面也取得了丰硕的成果。张鹏、赵哲、李攀提出了一种基于SPH法的燃油转运舱内液体晃荡数值模拟方法，考虑了液体在飞行器里受到的初始压力和重力作用，以及空气流动的影响，通过数值模拟分析并查找了影响舱内液体晃荡的因素，