改进的SPH方法及其在结构...水流固耦合问题中的应用研究.pdfVIP

改进的SPH方法及其在结构入水流固耦合问题中的应用研究

摘要

结构物入水砰击问题是一种复杂的流固耦合问题，涉及到自由液面大变形、瞬态

砰击载荷和多相界面非线性演化等诸多复杂力学现象。传统基于欧拉网格或拉格朗日

网格的数值方法在求解自由液面大变形和动边界问题时，面临网格畸变和物质界面追

踪重构等难题，需要额外实施较为复杂的数值技术，存在一定的困难和局限性。作为

一种拉格朗日粒子法，光滑粒子流体动力学方法（SmoothedParticleHydrodynamics,

SPH），在求解流体大变形和流固交界面问题具有较强优势，在自由表面流动、多相流

和流固耦合等领域得到了广泛的应用。

本文首先对结构入水砰击问题的研究现状进行了综述，主要包括理论、数值和实

验三部分，详细梳理了国内外研究进展，而后分析了SPH方法发展现状及现阶段SPH

方法求解结构入水问题时存在的局限性，确定了本文的研究内容及范围，为全文奠定

了基础。本文在基础SPH方法的架构上，瞄准固壁边界处理、大规模多相计算和流固

耦合分析难点，重点发展了单层粒子边界计算算法、轴对称多相SPH算法和多分辨率

SPH-FEM耦合算法，并基于发展的方法实现了对多类结构入水流固耦合问题尤其是救

生艇投放和飞机水上迫降等问题的模拟分析。

针对SPH方法计算入水问题时结构固壁边界离散精度和计算精度难以兼顾的难题，

发展了单层粒子边界处理方法，补齐了核函数支持域的缺失，避免了数值穿透现象，

实现了对复杂外形结构砰击载荷准确计算。通过简单几何结构入水和复杂曲面结构入

水等算例，充分验证了本文提出的改进的单层粒子边界计算方法和建立的SPH入水计

算模型，准确分析了结构入水过程中的砰击载荷、流场复杂流动、结构运动特征等。

开展了救生艇投放数值模拟，探明救生艇在入水过程中的自由液面演化和结构运动状

态特征，为救生艇的成功释放提供技术支撑。

针对SPH方法模拟大密度比的多相流问题时由多相时间步长很小引起的计算效率

较低问题，基于坐标系变换发展了轴对称SPH多相模型，能够实现对三维轴对称结构

多相入水的快速求解。同时采用体积自适应方法保证计算过程中粒子体积一致性，消

除了对称轴处数值奇异性，进而保证SPH计算精度。在充分校核的基础上，开展多相

流条件下的圆锥、小球和圆盘入水数值模拟研究，对比单向流结果，揭示了空气在带

底升角结构和平板结构入水过程中的影响。

针对SPH方法在模拟弹性体砰击入水这类结构和流体变形较大的问题时存在拉伸

I

哈尔滨工程大学博士学位论文

不稳定等、精度低等问题，发展了一套SPH-FEM流固耦合算法，两种方法间采用分

区松耦合的方式进行信息交互，既发挥了FEM方法在求解结构动力问题时计算精度好、

效率高的优点，又发挥了SPH无网格算法在模拟和表征复杂自由表面流动的天然优势。

此外，针对这种耦合求解器，还进一步提出了时空多分辨率耦合算法，采用额外积分

点等手段来保证任意流体粒子/结构单元尺寸比下的高效模拟。

针对目前结构入水问题中弹塑性响应对入水过程影响研究不足的问题，在发展的

时空多分辨率SPH-FEM流固耦合算法的基础上，引入板壳结构的延性损伤本构模型，

实现了对板壳结构及简易结构入水的弹塑性变形准确模拟，讨论了结构弹塑性对砰击

载荷和结构入水运动姿态的影响。此外开展了考虑结构弹塑性的三维客机水上迫降模

拟，分析了飞机迫降过程中的载荷和结构响应特征。通过对比刚性飞机入水结果，充

分说明了结构弹塑性在飞机迫降过程中的降载作用，为飞机的迫降问题的数值预报提

供了技术支撑。

关键词：光滑粒子流体动力学；结构入水；固壁边界条件；轴对称SPH模型；流

固耦合

II

改进的SPH方法及其在结构入水流固耦合问题中的应用研究

Abstract

Inmarineandoceanengineering,theproblemofwaterentryisacomplex