1. 流固耦合仿真概述与COMSOL工作流
1.1 什么是流固耦合
流固耦合（Fluid-Structure Interaction，FSI）研究流体流动与固体变形之间的双向作用：流体压力与黏性剪切力使固体产生位移和应力，固体的运动又反过来改变流场边界与流动状态。判断一个物理问题是否属于FSI，关键看两点：一是流体载荷是否足以引起可观测的结构变形或运动，二是结构变形是否显著改变流道几何或流动条件。若结构变形极小、对流场几乎无影响，通常可简化为“单向耦合”——先算流场，再把压力映射到固体上做静力分析；若变形量级与特征尺寸可比（例如柔性薄膜、血管壁、阀门开闭），则必须做“双向耦合”，每个时间步内交替求解流场与固体场并交换界面数据。
1.2 分类与工程应用场景
按耦合机制可分为三类：界面耦合（如管道振动、水弹性问题，流体与固体在交界面上交换力与位移）、体积耦合（如多孔介质中的渗流-变形，流体存在于固体骨架内部）、接触耦合（如活塞与缸套间的油膜润滑）。按求解策略又可分为整体式（monolithic）与分离式（partitioned）：整体式把流体与固体方程组装进同一矩阵求解，稳定性好但内存开销大；分离式分别求解两个物理场，通过界面插值传递数据，灵活但需处理“附加质量效应”带来的收敛困难。
典型工程场景包括：血管支架在脉动血流下的变形与疲劳评估、飞机机翼颤振与气动弹性分析、水轮机叶片在高速水流中的振动、柔性电子器件在风载下的响应、以及阀门开启过程中的瞬态流体力与阀芯运动。这些场景的共同点是：流体惯性、结构刚度与阻尼三者之间的匹配关系决定了系统是稳定还是发散。
1.3 COMSOL Multiphysics建模流程
COMSOL采用“模型开发器”树状结构组织仿真流程，核心步骤可归纳为：几何建模 → 材料与物理场选择 → 边界条件与耦合设置 → 网格划分 → 求解器配置 → 后处理。对于FSI问题，推荐使用“流固耦合”多物理场接口，它会自动在流体与固体交界面上生成“流固耦合边界”，并管理位移与力的传递。实际操作中，先建立流体域（如内流道）与固体域（如管壁），在“多物理场”节点下添加“流固耦合”，选择流体域与固体域，软件会自动识别交界面。若使用分离式求解，需在求解器中设置“分离式求解器”并指定流体与固体的求解顺序，通常先解流体再解固体，或采用“双向耦合”迭代直至界面残差收敛。
1.4 界面导航与关键设置
COMSOL界面左侧为“模型开发器”，包含全局定义、组件、网格、研究等节点；右侧为“设置”窗口，显示当前选中节点的参数。建模时需注意：单位制默认SI，但可在“全局定义”中修改；几何可通过参数化尺寸驱动，便于后续优化；材料可直接从内置库选取（如结构钢、水、空气），也可自定义；物理场需分别添加“层流”或“湍流”接口与“固体力学”接口，再通过“多物理场”节点耦合。求解前务必检查“研究”步骤中的“因变量”与“求解器配置”，对于瞬态FSI，时间步长应小于结构最小振动周期的1/20，且需开启“几何非线性”以捕捉大变形。
1.5 一个可复现的入门示例
以“二维轴对称管道内脉动流引起管壁振动”为例，说明关键设置。几何：流体域为半径5 mm、长100 mm的圆柱，固体域为壁厚1 mm的管壁。材料：流体为水（密度1000 kg/m³，动力黏度0.001 Pa·s），固体为结构钢（杨氏模量200 GPa，泊松比0.3，密度7850 kg/m³）。物理场：层流接口，入口设为“充分发展流动”，平均速度按正弦变化 u_in = 0.1·sin(2π·5t) m/s；出口设为压力0 Pa；固体力学接口，管壁两端固定约束。多物理场：添加“流固耦合”，自动识别内壁交界面。网格：流体域边界层第一层厚度0.01 mm，增长率1.2，共8层；固体域用扫掠网格，沿壁厚方向4个单元。研究：瞬态，时间范围0~0.5 s，步长0.005 s，求解器采用“分离式”，先解流体再解固体，每步迭代至残差小于1e-4。后处理：绘制管壁中点径向位移随时间变化曲线，若位移幅值在几个周期后趋于稳定，说明耦合收敛良好；若持续增大，则需减小时间步长或增加阻尼。