基于Hypermesh前处理与Fluent、Optistruct求解器的流固耦合分析（一）流场计算

一、概述        

        随着计算科学以及数值分析方法的不断发展，流固耦合或交互作用 （fluid structure coupling  或 fluid structure interaction）研究从 20 世纪 80 年代以来，受到了世界学术界和工业界的广泛 关注。流固耦合问题是流体力学（Computational Fluid Dynamics，CFD）与固体力学 （Computational  Solid Mechanics，CSM）交叉而生成的一门力学分支，同时也是多学科或多物理场研究的一个重要分支，它是研究可变形固体在流场作用下的各种行为以及固体变形对流场影响这二者相互作用的一门科学。了解流固耦合对于许多产品的设计至关重要。如果不考虑流体与固体之间的相互影响，则会导致产品性能被过高或过低估计。

        流固耦合一般分为单向耦合与双向耦合。如果结构变形非常小，并且可以认为结构的变形几乎不会对流场的各项参数产生影响，或产品本身不允许在流体的作用下发生较大的变形，这种情况下只需要先求解出流体与固体界面上的压强数据，并将压强数据传导到固体的表面进行结构力学计算。然而，如果结构发生大变形，流体的速度和压力场就会因此发生改变，此时我们需要将其作为双向耦合问题进行多物理场分析：流体流动和压力场会影响结构变形，而结构变形又反过来影响流体的流动和压力。实际工况中选择进行单向耦合分析还是双向耦合分析需要根据实际产品及作用工况进行判断。

        本文将执行一个单向流固耦合分析流程，先在Hypermesh前处理器进行流体域的建立和CFD网格划分，然后导入至Fluent求解器进行流场计算，得到流体与固体界面的压强信息，随后将Fluent中计算得到的压力信息映射至结构网格上，并使用Optistruct求解器进行结构力学分析。

二、问题描述

        在一个半圆形管道内部有一凸起物体，管道内水流速度为10m/s

        该凸起物体为空心结构，内部有加强筋，加强筋与外型面壁厚都为2mm，以下为凸起物体的内部结构示意

        将计算在恒定的水流下，该凸起结构受到水流冲击后的变形及应力情况。

三、流场计算

（1）流体域建模

        导入几何模型至Hypermesh

        提取管道内表面与凸起物体的外表面，并将管道两头封堵上，并修复拓补关系，形成一个封闭的流体域空间，将管道的一端作为流体的入口，另一端作为出口，如下图所示。  

        将各个面根据其边界条件类型的不同组织分布到不同的component里面，并按照类型进行命名，wall类型边界条件则将component的命名以wall为起始，入口类型边界条件以inflow为命名起始，出口类型边界条件以outflow为命名起始，流体网格则以fluid为命名起始，固体网格则以solid为起始，如下图。按照这种规则命名导入Fluent后边界条件类型可以被自动识别，只需调整其参数即可。

        用2d>automesh命令将边界面划分网格，将需要关注的凸起位置的网格细化，管道壁则可以用稍粗的网格，减少计算量。

        使用3d>CFD tetramesh命令进行流体网格划分,，选择需要做边界层的comp和不需要做边界层的comp，两种comp需要能构成完整封闭的空间才能生成实体流体域网格：

        流体网格生成完毕，可查看其内部实体网格：

        新建一个以fluid为起始命名的comp，将刚才生成的实体的四面体流体网格移动到该comp内，完成后的comp如下：

（2）导出网格        

        到这一步网格已经做好了，现在将该网格模型导入Fluent。先将模板切换到CFD：

        导出网格，导出过程中的提示直接确定。

        （3）导入Fluent并计算流场

        打开Fluent软件，切换到solution模块下：        

        导入CFD网格：

        注意单位问题，在Hypermesh内一般使用的是毫米单位，Fluent默认使用的米单位，因此需要将模型整体缩小1000倍

        设置入口边界流入速度：

创建水的流体材质：

将水的材质赋给流体域：

        求解：

计算结果-表面静压：

        到这一步便完成了从Hypermesh前处理划分网格到Fluent里面进行流场计算，得到关注位置的压强分布，下一篇博客将展示如何将流场计算结果单向耦合至结构网格上，进行结构力学计算。

        注：本文仅展示分析流程，对计算准确度暂无细纠。