要不是那个大桥塌了，谁在乎流固耦合
说起来挺讽刺的。很多重大工程进步，都是被灾难逼出来的。1940年塔科马海峡大桥的颤振倒塌视频，学土木和机械的应该都看过——风速才19m/s，桥就像纸片一样扭麻花，最后垮了。那时候人们才意识到，风不是静载荷，它和结构会“互动”的。这就是流固耦合（FSI）的经典案例。不过咱们工程师讲究实用，今天不追悼历史，聊聊我自己趟过的一些坑，和心得。
FSI这玩意，说穿了就是流体压力和结构变形互相影响。你吹一口气，树叶摇一摇，那就是最简单的流固耦合。但要命的是，工程上的FSI往往是非线性的、多物理场的，还带湍流。比如海上立管的涡激振动、飞机的颤振、心脏瓣膜的血液作用——别怕，我主要聊机械领域的，像阀门开闭时的水锤啊、换热器管束的流致振动啊，嗯，一堆麻烦。
强耦合弱耦合，分清楚才有活路
刚入行那会，我傻傻分不清。上司丢给我一个泵的流固耦合项目，我直接把Fluent算的压力场导出来，当作恒载加到结构上——结果自然是不准的！这就是单向耦合与双向耦合的差别。其实更准确的说法是弱耦合和强耦合。弱耦合就是交替求解，流体算一步，结构算一步，数据传一次；强耦合则是流体和结构方程同时求解，什么雅可比矩阵全耦合进去——计算量恐怖。
大多数工程问题，弱耦合够了。比如风力发电机叶片，变形小，对气流场影响不大，单向足矣。但碰上大变形，像血管支架展开，就必须双向，不然变形压扁流道，流场早变了。这里有个判据：特征时间比！如果结构响应周期远大于流体变化周期，可以解耦；反之必须考虑耦合。举个数值：韦伯数We>10^-3？或者无量纲数，自己翻文献去。我习惯直接试算——先跑个单向，看变形是否超过了边界层厚度的5%，超过了，重来吧，双向走起。
网格！网格！网格！重要的事说三遍
做FSI，最让人崩溃的不是算法，是网格。动网格技术，ALE（任意拉格朗日-欧拉方法）用得最多。但ALE有个坏毛病：网格变形大了，单元畸变，计算发散。你辛辛苦苦建了个模型，跑了三天三夜，结果出来个负体积——想砸电脑。我有次做生物瓣膜仿真，瓣叶开启，网格扭曲得跟麻花似的，最后乖乖用浸入边界法，虽然精度差点吧，好歹算出来了。
还有一个坑：数据映射。流体网格和结构网格往往不一致。流体密集，结构稀疏。压力插值不当，会在界面产生虚假应力。径向基函数（RBF）映射比保形映射稳，但计算量大。我的经验是，结构网格别太粗，至少流体边界层内要对应得上。另外，时间步长也得讲究，算FSI你得满足库朗数条件吧？流体那边CFL<1，结构这边也得跟上，最好取流体时间步的1/10到1/100？看问题。做VIV（涡激振动）的，必须抓得住涡脱落频率。
软件，是伙伴也是噩梦
市面上工具不少。Ansys Workbench里Fluent+Mechanical是最常见的组合，有系统耦合器，界面还算友好——就是许可证死贵。OpenFOAM+CalculiX是开源党的爱，但调试起来真要命。COMSOL强耦合直接解，对于小规模问题简直是神器，前处理爽啊，可三维大规模一跑就内存溢出。Abaqus有CEL（耦合欧拉-拉格朗日），做极端变形，比如鸟撞、液体晃荡，特别好用。也踩过坑：当年用某个版本（不提了），两个求解器时间不同步，界面数据丢了，结果完全不准——所以验证很重要！总之一句话：别迷信软件，它只是个工具。
参数调得怀疑人生？来点土办法
FSI收敛性是个玄学。你设置亚松弛因子，流体0.3，结构0.8，算得慢但稳。有时候压力-速度耦合换个算法，SIMPLE换成PISO就能过。最大的痛苦是：你不知道是物理上振荡，还是数值振荡。我的土办法：先跑流体（结构固定），看出不出奇异的压力脉动；再跑结构（加正弦载荷），看响应对不对；最后合起来。监控点位移和力的时间历程，搞不好就得重来。还有个野路子：如果只关心稳态，用好伪瞬态或者阻尼，骗过求解器也是本事。
这些年，ML（机器学习）也开始抢饭碗了。用降阶模型代替部分FSI，实时仿真不是梦。不过，暂时还替代不了老工程师的判断——哪个地方容易坏，你瞄一眼就心中有数，对吧？
好了，该去洗瓶子了。明天还有一根立管的作业等着呢。流固耦合，痛并快乐着。
