流固耦合现象在现实生活中随处可见,如何对其真实地模拟,一直是计算机图形工作者研究的热点和难点之一。流固耦合模拟的研究成果在自然灾害仿真、影视动画的特效制作和虚拟现实游戏的开发中有着重要的应用。目前,对流固耦合的模拟都是采用基于粒子的方法,尤属光滑粒子动力学(Smoothed Particle Hydrodynamics,SPH)方法应用最为广泛。虽然SPH方法在模拟大变形和自由表面流动等问题上有显著的优势,但是在模拟固体时存在着边界拉伸现象。另外,当前对流固耦合的模拟,大多都是将固体看作刚体进行研究,没有考虑流体与可变形固体的交互。针对这两个方面,本文利用SPH方法模拟流体域,对可变形固体域则采用有限元(Finite Element Method,FEM)方法模拟,通过耦合两种算法,从而完成流体与可变形固体的交互。FEM方法模拟固体时精度高且容易处理复杂的边界,而且FEM方法与SPH方法一样,都是基于拉格朗日格式的,便于在同一时间步长内进行求解。本文利用SPH-FEM耦合算法进行流固耦合的模拟,主要工作如下:(1)利用SPH方法对流体域进行求解和模拟。对影响算法精度和计算效率的邻域粒子搜索算法与核函数进行了重点分析和研究。利用蛙跳(Leap-Frog)方法显式求解时间积分,更新流体粒子的运动状态。(2)对交互场景中的可变形固体,采用有限元的方法进行分析和求解。利用应变张量和应变率建立变形模型,对可变形固体在受力时的变形进行控制。(3)对于SPH-FEM耦合算法边界问题,采用罚函数方法对边界处粒子和固体单元施加接触力,用以限制流体粒子的穿透。对边界粒子进行筛选,减少碰撞检测,加快计算效率。统一SPH方法与FEM方法的时间步长。(4)采用基于各向异性核的方法对流体表面进行重建,该方法与传统的表面重建方法相比,能够构造出更加光滑的流体表面。