机翼颤振、血管血流、轮胎充气……流固耦合仿真到底怎么做？

详细信息

作者优秀: 优秀教师/意见领袖/博士学历/特邀专家/独家讲师
平台推荐: 内容稀缺

1月前浏览246

点击文尾阅读原文查看

前言：当流体流过结构，结构在流体载荷下变形，变形又反过来改变流场——这种"你影响我、我影响你"的相互作用，就是流固耦合（FSI）。从飞机机翼的颤振、血管中血流的脉动，到换热器的强度校核、轮胎的充气变形，流固耦合几乎贯穿航空航天、能源动力、汽车、生物医疗等所有高端装备领域。

随着装备轻量化与服役工况日益严苛，流体与结构的相互作用愈发显著，单一的流体或结构仿真已无法满足设计需求。ANSYS Workbench 集成了从单向耦合、双向耦合到特种耦合的完整求解工具链，可覆盖常规稳态、动态强耦合、极端冲击等各类工况。

本文结合仿真理论与工程实践，系统梳理各类耦合方法的原理、流程与适用场景，帮助工程师快速完成方案选型与实操落地。

一、机理先行：流固耦合的本质与 ANSYS 全栈仿真体系

流固耦合的本质，是流体域与固体域通过交界面完成力学参数的交互传递。交界面需同时满足力平衡与运动学连续两类守恒条件：

σ_f · n = σ_s · n　　（1）

v_f = v_s，d_f = d_s　　（2）

其中：

σ_f、σ_s：流体/固体应力张量，单位 Pa；
n：交界面单位法向量，无量纲；
v：速度，m/s；
d：位移，m。

约束条件：仅适用于黏性流体无滑移、无穿透交界面；自由液面、滑移壁面不适用式（2）。

按数据传递方向，工程上将其分为单向耦合与双向耦合两大类（也常分别对应弱耦合、强耦合的说法），这是仿真方案选型的首要依据。

图 1　流固耦合交界面载荷传递联合云图：管内流体静压（Fluent）→ 管壁 von Mises 应力（Mechanical）

相较于单一的流体或结构仿真，流固耦合能够精准复现真实服役工况中"流体载荷导致结构变形、结构变形反向扰动流场"的动态效应，是解决装备振动、应力失效、气动弹性失稳等问题的关键技术。ANSYS Workbench 搭建了全覆盖的流固耦合仿真体系：Fluent 负责流体求解，Mechanical 负责结构求解，System Coupling 作为多物理场数据交互中枢，LS-DYNA 应对显式极端动力学问题，形成适配常规稳态、动态瞬态、极端大变形、封闭腔体承压等工况的完整求解梯队，兼顾仿真精度与计算效率，为各类流固耦合问题提供标准化、精细化的数值求解方案。

二、稳态工况首选：单向流固耦合的原理与三步实操

单向流固耦合是工程应用最广泛、计算成本最低、迭代收敛性最好的耦合方式，适用于结构刚度大、受力变形小、结构变形对流场影响可忽略的稳态工况。其核心逻辑是"单方向数据传递"，结构端求解退化为线性静力平衡方程：

[K]{u} = {F_f}（3），

[K]：结构刚度矩阵，N/m；{u}：节点位移向量，m；{F_f}：流体映射载荷向量，N。仅将 Fluent 求解得到的流体压力、壁面剪切力等载荷映射至固体结构，结构响应不再反馈流场。

依托 Workbench 的顺序求解模式，整个流程分为三步：

第一步，Fluent 流场求解——完成流体建模、网格划分、湍流模型与边界条件设置，迭代获得稳定流场；
第二步，交界面载荷映射——通过 Workbench 数据关联功能，将流固交界面上的压力与剪切载荷精准传递至结构网格，借助网格对齐或高精度插值算法避免应力失真；
第三步，Mechanical 结构分析——定义材料参数、约束边界与预应力条件，完成应力、变形乃至预应力模态分析。

图 2　变流速三通管稳态流场静压云图（入口 600 kPa，支管分流交汇处呈滞止特征）

实操中需重点把控两个细节：

一是 Fluent 端需正确配置壁面边界并确认载荷输出，保证流体数据完整传递；
二是结构端约束需贴合真实工况，避免过约束、欠约束引入求解误差。

该方法广泛应用于液压管路应力校核、换热器强度分析、稳态气动载荷评估等常规工业场景。

图 3　航空液压管路弯管段 von Mises 应力云图（单向耦合结构端：弯曲段内侧应力集中，峰值约 170 MPa）

三、动态工况攻坚：双向流固耦合的关键技术与调试要点

当结构足够"柔"——大变形、高频振动、气动颤振等动态工况下，结构变形会显著改变流场边界、反向扰动流体分布，单向耦合的精度便不再够用，必须采用双向流固耦合。该方法由 Fluent、Mechanical、System Coupling 三大模块协同完成：System Coupling 作为数据交互中枢，在每个时间步内反复传递流体载荷与结构位移并进行多轮耦合迭代，直至多物理场残差收敛，从而精准还原动态耦合全过程。流体域求解不可压 Navier-Stokes 方程：

∇·v = 0　　（4a）

ρ_f (∂v/∂t + v·∇v) = −∇p + μ∇²v　　（4b）

结构域求解动力学方程：

M ü + C u̇ + K u = F_f(t)　　（5）

每个时间步内按下式松弛迭代并判定收敛：

u_(k+1) = u_k + ω (û_(k+1) − u_k)，‖û_(k+1) − u_k‖₂ ≤ ε　　（6）

其中：

ρ_f：流体密度，kg/m³；
μ：动力黏度，Pa·s；
M/C/K：质量（kg）、阻尼（N·s/m）、刚度（N/m）矩阵；
û：当前耦合步新解；
ω：松弛因子，无量纲，0 < ω ≤ 1（System Coupling 常用 Aitken 自适应取值）；
ε：耦合收敛容差，工程常取 1×10⁻⁴～1×10⁻³。

约束条件：式（4）要求连续介质（Kn < 0.01）、牛顿流体、马赫数 Ma < 0.3 方可按不可压处理；式（5）若含材料非线性须将 K 替换为切线刚度；瞬态时间步建议满足 Δt ≤ T_min/20（T_min 为关注模态的最小周期，s）。

图 4　复合材料机翼上表面压力分布云图（颤振分析，迎角 6°：前缘吸力峰，翼尖三维卸载）

相较单向耦合，双向耦合的难点集中在动网格配置与耦合参数调试。Fluent 端需启用动网格功能，根据变形特征选择算法——小幅弹性变形采用弹簧光顺法，大变形工况采用局部网格重构法，杜绝网格畸变与负体积问题；Mechanical 端需准确描述材料非线性、阻尼特性与初始状态，以适配柔性复合材料、金属簧片等特殊结构；System Coupling 端需合理设置数据映射方式、松弛因子与收敛准则，在计算精度与迭代效率之间取得平衡。

图 5　血管狭窄段血流速度云图（双向耦合收缩期峰值：喉部流速约 3.5 m/s，虚线为放大的弹性管壁变形）

该技术主要应用于机翼气动弹性（颤振）分析、血管血流耦合、安全阀动态启闭、风机叶片流致振动等高端复杂场景。

图 6　簧片阀开启过程流场速度云图（阀片最大开度时刻：阀口间隙高速射流及下游扩散）

四、极端工况破局：ALE 与静压流体单元两大特种技术

对于高速冲击、剧烈流体大变形、封闭腔体承压等非常规极端工况，常规双向耦合方法存在求解局限，ANSYS 提供两类专属特种耦合技术作为补充。

图 7　结构高速入水冲击压力云图（ALE 方法，v=30 m/s：砰击压力峰集中于楔形体舭部，自由液面飞溅抬升）

一类是 LS-DYNA 的 ALE（任意拉格朗日-欧拉）方法：固体采用拉格朗日描述精准追踪界面，流体采用欧拉描述适配大变形，二者融合可高效求解水下爆炸、结构高速入水、高速碰撞等强非线性瞬态耦合问题。仿真中需通过 K 文件关键字定义耦合约束、流体状态方程等核心参数，保障极端工况下的求解精度。

图 8　水下爆炸冲击波传播与结构反射压力云图（LS-DYNA：球面波阵面抵达壁面，反射区压力显著升高）

另一类是静压流体单元法，专门面向无宏观流动的封闭腔体场景：将腔内流体等效为均压的静压流体单元（Hydrostatic Fluid Element，Mechanical 中的 HSFLD 系列单元），腔内压力随结构变形引起的体积变化自动更新、均匀作用于腔壁，无需求解流场细节，大幅简化模型、提升效率，典型应用包括轮胎充气变形、密闭容器承压仿真等。封闭腔体压力-体积关系（等温假设）：

p·V = p₀·V₀　　（7）

p、V：当前腔压（Pa）与腔体容积（m³）；p₀、V₀：初始状态值。约束条件：腔体完全封闭、无泄漏；腔内无宏观流动、压力处处相等；等温假设（非等温时改用 pV/T = 常数）；不计流体惯性，故不适用于晃动类问题。

如果气体质量变化时（如充气增压过程），改用：

pV = mR_sT　　（8）

m：气体质量，kg；R_s：气体常数，J/(kg·K)；T：温度，K。