首页/文章/ 详情

FLUENT单向流固耦合分析案例

5月前浏览808

1 前言

单向流固耦合是流固耦合分析的简化形式,其核心特征为流场对固体场施加单向载荷作用,固体变形不反向影响流场计算,仅满足弱耦合条件。在控制方程层面,流体域通过壁面剪切应力与压力将载荷传递至固体域,固体域求解应力、应变与位移,但不更新流体域网格与边界条件,流体控制方程与固体控制方程独立求解,仅通过载荷传递实现耦合。

该方法适用于固体变形极小、对流场扰动可忽略的场景,如刚性结构受力、小变形构件强度校核、稳态流场下结构响应分析等。相较于双向耦合,单向FSI计算成本低、收敛稳定、求解效率高,无需迭代更新流固交界面,可大幅缩短仿真周期。其本质是以工程精度为前提的近似耦合,在满足小变形假设时,可在计算效率与结果可靠性间实现良好平衡,为结构强度校核与流场特性分析提供高效数值手段。

举个最直观的例子:风吹过桥梁,风(流体)会对桥梁(固体)产生压力和推力,让桥梁轻微晃动,但桥梁的这种小晃动,并不会改变风的流动方向和速度——这就是典型的单向流固耦合。简单说,就是“流体影响固体,固体不影响流体”,是一种“单向奔赴”的耦合关系。

今天,我们做一个单向流固耦合的案例。

2 问题描述

我们考虑如下的管道模型:假设流体从下端入口流入,从右上方倾斜角度流出,入口端和弯头上游端均固定。在这种情况下,受流体冲击影响,倾斜管段将发生如箭头所示的变形。

3 建模、网格与求解

对于流固耦合问题,建议在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中定义的是绝对压力,则此处也为绝对压力,这种情况下需要在外壁面施加参考压力才能正确对应内外压差。

管道的变形结果如下,达到了预期效果。

结语

通过本案例,我们完整演示了单向流固耦合分析的流程:从几何建模、网格划分,到流体计算、压力数据传递,再到结构静力学分析。单向流固耦合作为一种高效的工程近似方法,在小变形假设成立的情况下,能够在保证结果可靠性的同时大幅提升计算效率,是工程实践中常用的分析手段。

来源:仿真与工程
MeshingFluent MeshingFluentWorkbench静力学非线性控制管道
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-03-12
最近编辑:5月前
余花生
签名征集中
获赞 259粉丝 606文章 391课程 0
点赞
收藏
作者推荐

讨论一个等效热阻问题

正文共:927字4图预计阅读时间:3分钟今天,我们讨论一个等效热阻问题。有如下一个固体结构,由三个长方体组合而成,外部大长方体为材质1,中间两个小长方体为材质2,材质1和材质2具有不同的导热系数。假设大小长方体完美接触,因此忽略接触热阻。如果我们要计算固体结构在z向的导热能力,严谨的做法是把三个长方体都创建计算域。但是如果为了简化计算,则可以考虑只创建大的长方体,此时需要获得一个可以表征这个固体结构整体导热能力的等效导热系数。我们知道,多孔介质导热系数的计算采用体积加权平均,这并非一个普适的物理定律,而是一种在特定条件和理想假设下非常有效且常用的近似模型。体积加权平均在物理上对应的是“并联模型”。当实际情况偏离并联模型的假设时,体积加权平均会失效或不精确。将多孔介质想象成由固体骨架和孔隙流体并列排列组成的复合材料块。热量可以同时、平行地通过固体路径和流体路径传导。其等效导热系数计算式为:与并联模型对应的是串联模型,等效导热系数如下,它使用调和平均,对应热量必须依次穿过固体层和流体层,串联模型最常见的场景是复合保温层。我们思考,其实多孔介质模型思想完全可以用于纯固体导热,比如本案例开头所述的固体结构。根据串并联模型,我们可以在z向将整个固体域分为A、B、C三个区域,这三个区域在z向为串联,而区域C又是并联模型(材质1和2),故此区域可以考虑多孔介质模型。这里我们假定材质1和2的导热系数分别为202.4W/m℃和20W/m℃,根据模型尺寸可得区域B的材质2体积为0.00036m3,材质1的体积为0.00064m3,即体积分数为0.36和0.64。于是,根据多孔介质模型,区域C在z向的等效导热系数为136.736W/m℃。区域A和B材质相同,整合为一体,则对于整个固体域区域C而言,A+B和C的体积分数均为0.5。根据串联计算式,整个固体域在Z向的等效导热系数为163.211W/m℃。我们在FLUENT分别采用简化方法(等效导热系数163.211W/m℃)和非简化方法计算上述固体结构的导热能力。计算过程太简单了,不做描述了。在z向顶部和底部施加60℃温差,简化方法和非简化方法的传热功率分别为489.65W和479.95W,也就是说简化模型和真实模型的误差约为2.02%。来源:仿真与工程

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈