桥梁风荷载流固耦合分析
理论与物理
它们之间有什么区别呢？
控制方程
桥梁的风致响应用什么方程来描述呢？
Scanlan的颤振微分方程是基础。单位长度上的升力 $L$、阻力 $D$、力矩 $M$ 用自激气动力系数（颤振导数）表示。
其中 $K = B\omega/U$ 是折算频率，$H_i^*$, $A_i^*$ 是颤振导数，通过风洞试验或 CFD 求得。$B$ 是桥面宽度，$h$ 是竖向位移，$\alpha$ 是扭转角。
颤振导数不能从理论上求得吗？
对于薄翼理论可以从Theodorsen函数解析求得，但桥梁截面是钝体（bluff body），所以不得不依赖风洞试验或CFD。最近的方法是采用强迫振动法（forced oscillation）通过CFD来识别 $H_i^*$, $A_i^*$。
临界颤振风速
临界颤振风速怎么求呢？
求解二自由度（竖向位移 $h$、扭转 $\alpha$）耦合系统的特征值问题。
将颤振导数代入右边的自激气动力项，系统阻尼变为零时的风速即为临界颤振风速 $U_{cr}$。本州四国联络桥的明石海峡大桥设计要求是 $U_{cr} > 78$ m/s。
明石海峡大桥的“风洞试验一万小时”——桥梁FSI理论守护的最长悬索桥
全长3,911米、中央跨距1,991米的明石海峡大桥（1998年竣工）是世界最长的悬索桥。这座桥的设计中，“将颤振速度保持在台风时最大风速80m/s以上足够高”是最大的课题。设计团队使用1:100的缩尺模型进行了累计超过一万小时的风洞试验，对主塔间加劲梁的截面形状进行了反复试错优化。最终采用的纤细箱梁截面，在抑制卡门涡街产生的同时，确保了颤振速度达到设计风速的1.7倍以上。对于这种规模的桥梁FSI设计，目前CFD仍被用作风洞试验的“前期筛选”，将待试验的截面形状候选方案通过CFD计算从数十个缩减到数个，再进行试验，这是标准的工作流程。
各项的物理意义
- 结构-热耦合项：温度变化引起的热膨胀诱发结构变形，变形又影响温度场。$\sigma = D(\varepsilon - \alpha \Delta T)$。【日常例子】夏天铁轨伸长导致缝隙变窄——温度升高→热膨胀→产生应力的典型例子。电子基板在焊接后翘曲也是不同材料热膨胀系数差异导致的。发动机的汽缸体因高温部分和低温部分的温差产生热应力，最坏情况下会导致裂纹。
- 流体-结构耦合（FSI）项：流体压力·剪切力使结构变形，结构变形又改变流体域的双向相互作用。【日常例子】强风下悬索桥的缆索振动（涡激振动）——风力使结构摇晃，摇晃的结构改变风的流动，进而放大振动。心脏的血流与血管壁的弹性变形、飞机机翼的颤振（气动弹性不稳定性）也是典型的FSI问题。有时仅单向耦合即可，但变形较大时双向耦合是必须的。
- 电磁-热耦合项：焦耳发热 $Q = J^2/\sigma$ 引起温度升高，温度变化改变电阻的反馈循环。【日常例子】电暖炉的镍铬丝通电后发热（焦耳热）变红——温度升高电阻改变，电流分布也变化。IH电磁炉的涡流发热、输电线路温度升高导致的垂度增加也是这种耦合的例子。
- 数据传递项：通过插值解决不同物理场间网格不匹配的问题。【日常例子】天气预报中结合“气温数据”和“风的数据”计算体感温度时，如果各自的观测地点不同就需要插值——CAE的耦合分析中，结构网格和CFD网格通常也不一致，因此界面处的数据传递（插值）精度直接关系到结果的可信度。
假设条件与适用范围
- 弱耦合假设（单向耦合）：当一方物理场影响另一方而反向影响可忽略时有效
- 需要强耦合的情况：FSI中的大变形、电磁-热耦合中温度依赖性较强时
- 时间尺度的分离：各物理场特征时间差异较大时，可采用子循环提高效率
- 界面条件的协调性：需确认耦合界面处的能量·动量守恒在数值上得到满足
- 不适用的案例：三个以上物理场同时强耦合时，有时需要整体式方法
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 热膨胀系数 $\alpha$ | 1/K | 钢：约12×10⁻⁶、铝：约23×10⁻⁶ |
| 耦合界面力 | N/m²（压力）或N（集中力） | 确认流体侧与结构侧的力平衡 |
| 数据传递误差 | 无量纲（%） | 插值精度依赖于网格密度比。5%以下是目标 |
数值解法与实现
基于CFD的颤振导数识别
请告诉我用CFD求颤振导数的步骤。
强迫振动法是标准方法。
1. 创建2D截面模型（仅桥面截面）
2. 施加正弦波的竖向振动 $h(t) = h_0 \sin(\omega t)$ 并运行CFD
3. 从升力的时间历程中分离同相分量（刚度项）和反相分量（阻尼项）
4. 识别 $H_1^*$〜$H_4^*$、$A_1^*$〜$A_4^*$
5. 对扭转振动 $\alpha(t) = \alpha_0 \sin(\omega t)$ 同样实施
CFD使用2D RANS或LES。湍流模型在桥梁截面方面有实绩的是$k$-$\omega$ SST。
不需要3D全桥模型FSI吗？
实务主流是“条带理论”，即通过2D截面分析求得颤振导数，然后将其应用于模态法的全桥结构模型。但是，当3D效应（端部涡、跨度方向相位差）重要时，也会进行3D CFD-CSD耦合。
涡激振动的CFD分析
涡激振动如何进行模拟呢？
进行2D截面的CFD-CSD耦合。用LES（或DES）求解桥面截面周围的流动，结构用弹簧-质量模型表示并进行耦合。
锁定现象发生在涡脱落频率 $f_s = St \cdot U/D$（斯特劳哈尔数 $St \approx 0.1$〜$0.2$）接近结构固有频率时。要用CFD正确预测锁定范围和振幅，需要足够的时间积分（100个振动周期以上）。
抖振分析
自然风的随机性如何处理？
频域法效率较高。以Davenport的风速谱作为输入，结合气动导纳函数和结构的传递函数求得响应谱。
$S_u$ 是脉动风速谱，$\chi$ 是气动导纳，$H$ 是结构的频率响应函数。时域中也有生成人工脉动风场（von Karman型）并作为CFD入口条件的方法。
风洞实验与CFD的相互验证——桥梁设计的现场实情
在桥梁的抗风设计中，目前风洞实验仍是主角，但CFD和FSI的应用正在增加。风洞中使用1/100〜1/200比例的模型，因此存在雷诺数与实桥差异很大的问题。CFD在原理上可以进行实尺度计算，也能再现风洞中难以再现的3维效应和桥面交通引起的紊流。现场的标准做法是“用风洞把握大致趋势，用CFD进行细节推敲”。
整体式方法
将所有物理场作为一个联立方程组同时求解。对强耦合问题稳定，但实现复杂，内存消耗大。
分区法（分离迭代法）
各物理场独立求解，在界面交换数据。易于实现，可利用现有求解器。适用于弱耦合。
界面数据传递
最近邻法（最简单但精度低）、投影法（具有守恒性）、RBF插值（对网格不一致鲁棒性强）。守恒性与精度的平衡很重要。
子迭代
在每个耦合步内进行充分的迭代，确保界面条件的协调性。残差基准基于各物理场的典型值进行缩放。
Aitken松弛
自动调整耦合迭代的松弛系数。防止过度松弛导致发散，是加速收敛的自适应方法。
稳定性条件
注意附加质量效应（流体-结构耦合中结构密度≈流体密度时）。不稳定时可应用Robin型界面条件或IQN-ILS法。
Aitken松弛的比喻
Aitken松弛类似于“平衡跷跷板”。一方推得太用力，另一方就会弹起，其反作用力又导致推得过猛——为了抑制这种振荡，自动调整推力大小的就是Aitken松弛。当耦合迭代振荡不收敛时，根据上次的修正量自动调整下一次修正量的自适应方法。
实践指南
基于日本道路桥抗风设计手册的流程如下。
1. 截面形状气动特性评估: 通过风洞试验或CFD获取气动三分力系数（$C_D$, $C_L$, $C_M$）
2. 静力稳定性验算: 确认发散风速
3. 涡激振动验算: 预测锁定风速范围和响应振幅
4. 颤振验算: 临界颤振风速需为设计风速的1.2倍以上
5. 抖振响应: 设计风速下的位移·应力最大值
6. 减振措施探讨: TMD、襟翼、整流罩等
CFD和风洞试验如何区分使用？
实务中以风洞试验为主，CFD作为补充使用。不过最近CFD在截面形状的参数化研究中很活跃。风洞试验一个截面要花费数百万日元，但用CFD可以快速评估设计变更的效果。
CFD网格设计
桥梁截面CFD网格需要注意什么？