一个基于fluent的柔性管道的流固双向耦合算例

一个基于fluent的柔性管道流固双向耦合算例

声明：

1. 本算例仅作为学术交流使用，未经作者允许不得转载。
2. 本算例使用的参数不准确，只是交待计算过程，计算思路及系统设置的注意事项，对于结果的精度不作探讨。

1. 物理模型

     物理模型为一根半径为5mm，厚度为0.1mm的管道，流体域在其内部，外壳为固体。使用ansys DesignModel对模型命名。进口压力为16000Pa，进口流速为周期性波动，其表达式为

如下图所示.

 

 2. 网格划分

网格划分分为两个阶段，一个为fluenti里面的网格划分，一个为结构仿真中的网格划分，这里分别使用静态结构分析和瞬态结构分析，其目的在于分析静态结构和瞬态结构的结果差异。出于对计算资源的考虑，管壁厚度选择一层网格，原则上需要3-5层网格进行网格划分，静态结构的网格划分如下图所示，使用四边形网格和面网格剖分。

 

Fluent的网格划分，使用fluent（带mesh网格划分功能）的功能模块，对流体域进行网格划分，考虑到进出口周期波动，故在添加局部尺寸一栏中，加入inlet和outlet两个，细化进出口两端的网格数。网格划分过程按照系统预设就可以，注意的是描述几何结构一栏中，选择几何结构同时有流体和固体区域和/或空区域组成，fluent网格划分算法在这些简单的图形上，问题不大。复杂的模型就需要使用ICEM作结构化网格划分，尤其是超声速的仿真中。

 

3. 工程数据的设定

双击静态结构中工程数据一栏，操作如下图所示；本算例泊松比为0.47，选择的是剪切模量为200Mpa，但其实系统会自动换算，只要选择合适自身的模量即可，操作完成后返回原界面

 

 

4.Fluent的设置

4.1 系统启动

系统使用双精度（Double Precision）和并行计算（Parallel），根据自身计算机的配置，选择适合的处理器个数和GPU个数。

 

4.2 通用界面的设定

进入通用界面，选择比例，调整适合的网格单位，这里是毫米。点击检查，系统通过自检，保证格网无负数后，选择瞬态运行。

 

4.3 材料的选择并更改属性

选择Fluent数据库内选择液态水后，点击复制，之后在创建与编辑材料的窗口内，以及自身材料的特性，修改密度和粘度，这里密度修改未1010，粘度修改为0.004.

4.4 边界条件的选择

在单元区域条件内，选择part-1，在流体窗口内选择液态水。

 

设置进口流速的表达式，如下图所示：

 

入口选择速度入口，这样计算收敛性较好，根据之前定下的边界条件，设定初始化的表压和速度

 4.5 动网格的设置

选择动网格进行计算，网格方法中选择光顺（不一定需要）和网格重构

 

在光顺法这里，选择扩散。在网格重构那里选择尺寸函数，并在最大长度尺寸一栏修改为10mm。避免过大的网格重构运算量。

 

 

在动网格一栏内点击创建/编辑，这是变形区域为流体区域。

设置系统耦合区域为流体的两个壁面，并设置单位高度，单位高度不好超过外型尺寸。

初始化系统，如果耦合成功，开始计算的按钮为灰色，设定计算步长，这里需要自行调试合适的步长，不是越小越好，步长越小有时候会导致系统的非线性项误差保留过多，容易放散。

5. 结构仿真设置

静态结构和瞬态结构均按照这个约束条件进行设定，流固耦合的边界面选择管道的内部。

 6. 耦合注意事项

拖动系统耦合的模块到workbench中，关联fluent和静态结构，关联fluent和瞬态结构，注意只能关联一次。

 

进入系统耦合模块后，修改仿真的结束时间和步长。

 

在数据传输一栏，右键增加数据传输类型，由于是双向耦合，故需要两次数据传输，一次为fluent到结构仿真，一次为结构仿真至fluent。完成之后开始仿真。

 

需要注意的是，系统只是一个计算器，如果设定不对，会产生非物理解，如下图所示。

 

对比静态结构和瞬态结构中，管道的平均总形变，平均等效弹性应变，和平均应力，如下图所示，可以在运行2秒内，系统的瞬态仿真结果平均总形变，平均等效弹性应变，和平均应力都叫静态仿真的要高。