流体绕过钝体时在尾流中交替脱落漩涡，产生周期性变化的脉动压力，激发结构在横向发生振动，这一现象称为涡激振动（VIV）。深海立管、跨海桥梁拉索、换热器管束与烟囱在风或海流作用下都可能发生涡激振动，严重的锁定（lock-in）会导致大幅振动与疲劳破坏。本项目的科学问题是：圆柱在均匀来流中何时发生涡激振动锁定？锁定状态下横向振幅与频率特征如何？振幅随约化速度如何变化？利用ANSYS Fluent与结构动力学耦合的流固耦合仿真可以再现涡激振动的发生与发展过程。
涡激振动的核心是漩涡脱落频率与结构固有频率的耦合。静止圆柱的漩涡脱落频率由斯特劳哈尔关系f_s=St·U/D描述，St约0.2。当约化速度Ur=U/(f_n·D)（f_n为结构固有频率）增大使脱落频率接近结构固有频率时，振动反馈调节漩涡脱落相位，使脱落频率被”锁定”在固有频率附近，振幅急剧增大。此时结构相当于一个非线性振子，其运动方程m·ÿ+c·ẏ+k·y=F_y(t)与流体Navier-Stokes方程通过动网格耦合求解。振幅随约化速度的响应呈典型的”驼峰”状：在锁定区（Ur约4至8）振幅大，离开锁定区后振幅回落，最大横向振幅通常可达0.5至0.6倍直径。
以二维圆柱为对象，直径D取0.1米，流体为水，来流速度按约化速度换算。在ANSYS Fluent中建立矩形计算域，入口距圆柱上游约10D、出口距下游约25D，上下边界距圆柱约10D。采用SST k-omega湍流模型或层流模型求解，圆柱表面为无滑移壁面。将圆柱建模为弹性支撑的刚性体，质量比取2、阻尼比取0.01，通过用户自定义函数或动网格模块将流体载荷传递给结构方程并更新网格，实现双向流固耦合。首先固定圆柱计算不同约化速度下的振幅响应曲线，然后在锁定约化速度下输出横向位移随时间的变化。
图17的左侧面板（a）给出锁定状态下圆柱横向位移y/D随时间的变化。位移在初始阶段较小，随后呈增长趋势并伴随明显的不规则调制（拍振），经过若干周期后振幅稳定在约0.5倍直径的平台，形成等幅正弦振动。拍振来源于脱落频率向固有频率靠近过程中的频率差，随着锁定建立，振动频率被吸入固有频率，拍振消失、振幅达到稳定。稳定振幅约0.5D与经典涡激振动实验结果一致，说明结构在该约化速度下处于完全锁定状态，漩涡脱落与结构振动之间形成了稳定的相位同步。
图17的右侧面板（b）给出最大横向振幅A/D随约化速度的响应曲线。曲线在低约化速度区振幅很小，接近零；当约化速度进入4.5至6.0的绿色锁定区间后振幅急剧上升，在约5.1处达到约0.55的峰值，随后随约化速度继续增大而回落，在高速区振幅重新减小。整个响应呈现典型的”驼峰”形态，锁定区间与峰值位置与St=0.2给出的共振条件吻合。振幅曲线的上分支与下分支（滞后）效应在低阻尼下可能出现双稳态，图中展示的是稳定分支。该响应曲线是涡激振动疲劳评估的核心输入，用于确定结构在流速变化时可能遭遇的最大振幅与循环次数。
两面板结合，左面板展示了锁定状态的时域振动特征，右面板给出了振幅随流速的频域响应包络，共同完成了Fluent流固耦合对圆柱涡激振动的定量分析。
本项目用ANSYS Fluent双向流固耦合模拟了弹性支撑圆柱的涡激振动，锁定约化速度约5.1，稳定横向振幅约0.5D，振幅响应呈驼峰形态。案例图以横向位移时间曲线与振幅-约化速度响应双面板呈现。后续可拓展：(1)考虑三维效应与端部条件，评估实际立管的涡激振动；(2)引入抑制装置（螺旋列板）的仿真，优化抑振设计；(3)统计锁定工况下的交变应力，开展疲劳寿命评估。
