CAE进阶之ANSYS Workbench 单双向流固耦合仿真计算快速入门

最新推荐文章于 2025-03-06 14:42:48 发布

CAE从入门到放弃

最新推荐文章于 2025-03-06 14:42:48 发布

阅读量3.7k

点赞数 69

文章标签：开发语言小程序 windows 经验分享程序人生

本文链接：https://blog.csdn.net/qq_34933927/article/details/145338755

版权

引言

输流管道的处于结构场和流场的共同作用下动态响应较为复杂，大多数工况下忽略流固耦合效应易造成误差。

管路流固耦合仿真计算分为两种，一种是单向流固耦合，单向流固耦合仿真过程中，先进行流场计算，得到对应边界条件下流场的载荷分布，随后将流场载荷（压力）施加于结构场上，单独进行结构场的动力学仿真；双向流固耦合则需要在流场和结构场中反复迭代，每完成一步流场计算，马上将所得载荷代入结构场得到该子步下结构场的变形，管路变形之后流域发生变化，流场也随之改变进行下一子步计算，并再次将流场载荷返回结构场，如此反复迭代。

商业软件中ANSYS提供了两种途径开展管路的流固耦合计算，主要是APDL和workbench，相对而言workbench对于初学者会有好一点。

流固耦合仿真计算过程：

问题类型		Workbench
问题类型		分析步	操作方法
单向流固耦合	稳态	Fluent+Static Structure	GUI
双向流固耦合（瞬态）		Fluent++Transient Structure+system couling	GUI

下面着重对Workbench的流固耦合步骤做一个讲解。

一、前期模型准备

建模：采用三维建模软件建立常规管路几何模型效率相对较低，在此处采用APDL建模的方式，便于管路的修改。管路建模需要建出流场域和管道的结构场。采用如下命令流，在APDL中运行即可。当然，初学者可以根据自己需要，在三维软件中建立所需要的模型，导出成STP格式或者IGES格式都可以，推荐stp。

FINISH  ! Make sure we are at BEGIN level   
/CLEAR,NOSTART  ! Clear model since no SAVE found   
/PREP7  
R1=0.005       !------------outer diameter-----------
R2=0.004       !-----------inner diameter------------
L1=1
L2=1
ALPHA=30*3.141592657/180
LR=0.4            !------------bench radius-----------------
K,1,0,0,0,  
K,2 ,0,0,L1,  
K,3 ,0,L2*SIN(ALPHA),L1,  
L,1,2
L,2,3
LFILLT,1,2,LR, ,   
LCOMB,ALL, ,0  
NUMCMP,LINE 
CYL4, , ,R2, ,R1  
CYL4, , ,R2 
VDRAG,2, , , , , ,1   
VDRAG,1, , , , , ,1 

NUMCMP,VOLU
IGESOUT,'PipeFluid','IGES',' ',1

二. 单向流固耦合

稳态计算

采用GUI的方式，在Workbench中分别拖入Geometry、FluidFlow（Fluent）和Static模块，并确定各模块相关步骤的衔接关系。

步骤1：在DM中导入利用APDL 建立的IGES格式的模型，并命名流体域和结构域的名称。

步骤2：在Fluid的mesh部分，需要禁用Pipe的结构部分，划分网格时采用四面体即可，并对inlet、outlet和Fsiwall进行命名，检查网格质量均值小于1，处于0.75附近，越高越好。

步骤3：Fluent中设置单位为毫米，模型为k-ε，边界条件设置时采用压力入口30MPa和压力出口0MPa，相应的水力直径需要修改为管道内径。随后进行初始化，有两种Initialization初始化方式均可以使用，最后Run Calculate，设置迭代步数。

步骤4：在Mechanical中禁用流场域，对结构场进行网格划分，对应尺寸与流场尺寸保持一致，在边界条件设置时，引入Support，可以使Elastic或者Fix Support。在管道内壁面引入流场计算所得FSIwall 的压力载荷（导入载荷过程较慢需要等待较长时间），并进行计算。

计算结果如下，可以输出流场流动特性云图和管路应力云图：

三. 双向流固耦合

瞬态计算

采用GUI的方式，在Workbench中分别拖入Geometry、FluidFlow（Fluent）和Transient Structure和System Coupling模块，并确定各模块相关步骤的衔接关系。

步骤1：在DM中导入利用APDL 建立的IGES格式的模型，并命名流体域和结构域的名称。

步骤3：Fluent中设置单位为毫米，模型为k-ε。将计算类型改为瞬态计算，设置重力加速度。边界条件设置时采用压力入口30MPa和压力出口0MPa，相应的水力直径需要修改为管道内径。

步骤4：需要进行动网格设置，勾选smooth和remeshing，其中smoothing中设置Diffusion参数为1.5，remeshing中设置参数默认并将size remeshing改为1，选择动网格区域为Fsiwall。随后进行初始化，有两种Initialization初始化方式均可以使用，最后Run Calculate，设置迭代步数。

步骤5：在Mechanical中禁用流场域，对结构场进行网格划分，对应尺寸与流场尺寸保持一致，在边界条件设置时，引入Support，可以使Elastic或者Fix Support。在analysis setting中设置计算时间最长1000s尽可能大一些，关闭自动步长设置，并设置按照子步进行计算，子步数为1。在管道内壁面设置FSI界面，并进行计算。

步骤6：在system coupling中分别进行analysis setting、data transfers和intermdiate restar data output。其中，analysis setting中主要设置计算总时长和时间步长。Datatransfers中勾选fsiwall和Fluid SolidInterface创建DataTransfers。输出根据子步数对应输出即可。