单向流固耦合是流固耦合分析的简化形式,其核心特征为流场对固体场施加单向载荷作用,固体变形不反向影响流场计算,仅满足弱耦合条件。在控制方程层面,流体域通过壁面剪切应力与压力将载荷传递至固体域,固体域求解应力、应变与位移,但不更新流体域网格与边界条件,流体控制方程与固体控制方程独立求解,仅通过载荷传递实现耦合。
该方法适用于固体变形极小、对流场扰动可忽略的场景,如刚性结构受力、小变形构件强度校核、稳态流场下结构响应分析等。相较于双向耦合,单向FSI计算成本低、收敛稳定、求解效率高,无需迭代更新流固交界面,可大幅缩短仿真周期。其本质是以工程精度为前提的近似耦合,在满足小变形假设时,可在计算效率与结果可靠性间实现良好平衡,为结构强度校核与流场特性分析提供高效数值手段。
举个最直观的例子:风吹过桥梁,风(流体)会对桥梁(固体)产生压力和推力,让桥梁轻微晃动,但桥梁的这种小晃动,并不会改变风的流动方向和速度——这就是典型的单向流固耦合。简单说,就是“流体影响固体,固体不影响流体”,是一种“单向奔赴”的耦合关系。
今天,我们做一个单向流固耦合的案例。
我们考虑如下的管道模型:假设流体从下端入口流入,从右上方倾斜角度流出,入口端和弯头上游端均固定。在这种情况下,受流体冲击影响,倾斜管段将发生如箭头所示的变形。
对于流固耦合问题,建议在Workbench平台中开展。虽然FLUENT内置了结构力学模块,但功能相对有限,不建议用于专业细致的模拟。在Workbench中创建如下工作流程,本案例采用Fluent Meshing划分流体网格(多面体)。
需要注意的是,结构静力学模块的几何模型来源于流体计算模块,无法单独进行修改,需要在上游模块进行处理。我们在SCDM中进行模型前置处理,必须同时包含流体域和固体域,且两者交界面需进行共享(Share),最终效果等同于DM中的多体部件。为减小计算量,本案例仅取管道的一半进行计算。
强烈建议对模型的各个部分进行Part命名(流体域是必须的),便于后续边界条件的设定。
对流体域划分多面体网格。注意,在Update Regions中将固体部分类型设置为Dead,以此抑制固体域,避免划分固体网格。需要特别指出的是,抑制固体域并非必须操作,需根据实际情况而定——若还需考虑固体的导热问题,则固体域仍需保留。
接下来进入Fluent进行流体部分的计算。介质为空气,将入口设置为30m/s的速度入口边界,出口设置为表压为0的压力出口条件,详细过程从略。
壁面的压力分布以及入口的静压结果如下图所示。需要特别指出的是,由于壁面压力需要传递给静力学分析模块作为边界条件,因此在CFD分析时,要尽可能获得准确的壁面压力结果。
接下来是静力学计算部分。为夸张地模拟变形效果,这里我们先进行了“取巧”操作——将钢材的弹性模量设置为1.5×106Pa(原值为2×10¹¹Pa),这样就可以用线性模型获得较大的变形量。但在实际应用中,这种方法并不可取:能发生很大变形的结构,极大概率涉及非线性过程,需要使用非线性模型(如超弹性模型等),不过这些模型的计算量相当巨大。
进入静力学模块,首先要抑制流体域。
划分如下结构化网格。该模型较为简单,可采用Sweep方式获得高质量的结构化网格。
设置如下力学边界条件:入口端面和弯头上游边缘线均设置为固定边界,内壁面的压力来自CFD计算结果,中间剖面设置为无摩擦支撑边界(等同于对称边界)。
壁面压力边界的导入方式如下。注意定义传递面的CFD Surface要选择正确,前文强调的Part命名作用就在于此。
至此,力学边界条件设置完毕。这里需要着重说明压力边界:此处的压力边界完全从CFD传递而来,其代表的实际意义必须与CFD中的定义完全一致,同时也要与力学边界对应。
例如:若CFD中定义的压力为表压,则此处也为表压,因此不需要对外壁面进行力学边界设置,相当于默认对内外壁施加了压差;若CFD中定义的是绝对压力,则此处也为绝对压力,这种情况下需要在外壁面施加参考压力才能正确对应内外压差。
管道的变形结果如下,达到了预期效果。
通过本案例,我们完整演示了单向流固耦合分析的流程:从几何建模、网格划分,到流体计算、压力数据传递,再到结构静力学分析。单向流固耦合作为一种高效的工程近似方法,在小变形假设成立的情况下,能够在保证结果可靠性的同时大幅提升计算效率,是工程实践中常用的分析手段。