基于有限质点法和光滑粒子流体动力学的流固耦合计算方法研究

RESEARCH ON COMPUTATIONAL METHOD OF FLUID-STRUCTURE INTERACTION BASED ON FINITE PARTICLE METHOD AND SMOOTHED PARTICLE HYDRODYNAMICS

  • 摘要: 流固耦合是结构工程研究热点与难点之一,涉及流体自由液面、结构强非线性响应及流固界面耦合等关键问题。本文提出了基于有限质点法(FPM)和光滑粒子流体动力学方法(SPH)的流固耦合计算方法。首先给出了结构模型计算的有限质点法,包括质点运动控制方程和平面梁单元内力计算方法。然后实现了流体模型计算的SPH方法,给出了流体控制方程的离散格式,以及拉伸稳定性控制方法和时间积分方案。在此基础上建立了FPM-SPH界面耦合方案,给出了耦合界面粒子布置方法,提出了适用于SPH粒子和结构质点的界面压力传递和界面位移协调方法,建立了界面处结构质点和流体粒子间的物理量对应关系,从而搭建了流固耦合分析框架。采用本文方法对溃坝水流冲击弹性梁、水压力作用弹性门与膜结构积水等3个算例进行了建模分析,得到的流体运动和结构变形过程与试验或文献结果接近,位移时程曲线也基本一致,验证了本文方法的有效性。

     

    Abstract: Fluid-structure interaction (FSI) is a significant area of research and a challenging topic in the field of structural engineering. FSI involves a range of crucial issues, including the behavior of free surfaces, the analysis of strong nonlinear structural responses, and the coupling at the fluid-structure interface. In this study, a novel computational method for fluid-structure interaction is proposed by combining the Finite Particle Method (FPM) and Smoothed Particle Hydrodynamics (SPH). Firstly, the FPM for the computation of structural model is presented, which includes the calculations of the motion equations of motion of particles and the internal forces within planar beam elements. Subsequently, the SPH method for the computation of fluid model is implemented, including the discrete formulation of the fluid governing equations, the control strategies for tension stability controlling, and the scheme for time integration. On these bases, an approach for the computation at FPM-SPH interfaces is put forward. Specifically, the arrangement of virtual particles at coupled interfaces is provided, and the algorithms for passing interface pressure and coordinating interface displacement are proposed. In these algorithms, the correspondence of physical quantities between structural particles and fluid particles at the interface is established. The computational framework of the method proposed is thusly constructed. Three examples, i.e., dam-break flows impacting an elastic beam, an elastic gate under water pressure, and water falling on a rectangular membrane structure, are modelled and analyzed using the method proposed. The numerical results of the fluid motion, of the structural deformation, and of the displacement time-histories are consistent with the experimental results or the results in the existing literatures, demonstrating the effectiveness of the method proposed.

     

/

返回文章
返回