Star ccm 与Abaqus 协同仿真（结构耦合）

最新推荐文章于 2025-02-24 11:59:40 发布

算法资料吧！

最新推荐文章于 2025-02-24 11:59:40 发布

阅读量2.9k

点赞数 31

文章标签： FSI starcm abaqus

本文链接：https://blog.csdn.net/jianqimingtian/article/details/144625644

版权

概述

本教程演示如何通过使用 Simcenter STAR-CCM+ 和 Abaqus 运行机械协同仿真案例。协同仿真涉及两个代码之间的强耦合。数据以频密间隔（称为耦合步）进行交换。利用求解器之间的这种水平的通信，您可以获得整个流体/固体交界面的完全求解。

在机械耦合中，Simcenter STAR-CCM+ 将拉力负载传递给 Abaqus（压力 + 壁面剪应力），而 Abaqus 将位移传递给 Simcenter STAR-CCM+。在 Abaqus 中，拉力负载作用于固体结构的表面。在 Simcenter STAR-CCM+ 中，位移作为网格变形器输入使用。

在协同仿真中，通过使用 SIMULIA 协同仿真引擎在 Simcenter STAR-CCM+ 和 Abaqus 之间自动交换数据。这种方法不同于基于文件的方法，这是因为两种代码之间的耦合相当强烈，而且 Simcenter STAR-CCM+ 和 Abaqus 同时运行。

要创建协同仿真案例，必须分别为流体和固体域准备单独的模型。在 Abaqus 创建固体模型，并在 Abaqus 输入文件中定义。在 Simcenter STAR-CCM+ 模拟中创建流体模型，并且在该模拟中定义适当的边界条件和物理属性。对于每个模型，必须指定一组协同仿真参数。这些参数确定构成流体结构交界面的边界/表面、由交界面导入/导出的场数据、耦合步骤的频率等。当设置完成后，使用 Simcenter STAR-CCM+ 来运行协同仿真。

在运行协同仿真之前，先分别运行流体模型和固定模型的单独案例以便确定这些模型各自收敛，这是很好的做法。在本案例中，可使用施加于平板表面的压力负载 *DSLOAD 和采用估计的恒定压力值，来单独运行 Abaqus 模型。Simcenter STAR-CCM+ 模拟可通过使用固定平面在不进行网格变形的情况下来运行。

本教程模拟固定在底部、处于 10 m/s 空气流率之下的的弹性舌片。平板的质量和刚度条件是，第一弯曲模式的频率接近 4 Hz，第一扭曲模式的频率接近 20 Hz。对于固定平板，10 m/s 入口速度下的涡流脱落接近 20 Hz，因此您期望涡流脱落和第一扭曲模式之间的实现共振。

在本教程中，通过运行准稳态模拟来初始化流体域。随后使用初始流体场作为协同仿真的初始条件。初始流体场建立后，将求解器设置更改为协同仿真所需要的值。运行协同仿真时，平板被松开，并在风中弯曲。在本教程中，运行协同仿真 0.125s，这大约是第一弯曲模式的基础时长的一半，以及大约 2.5 个涡流脱落时长。运行案例更长时间（大约 4 s）会发现，第一弯曲模式衰减，而第一扭曲模式主导平板运动。

前提条件

“Abaqus 机械耦合”教程中的说明适用于已熟悉 Simcenter STAR-CCM+ 的某些操作方法的用户。

本教程演示使用 Abaqus 2021 进行协同仿真。要查看所有支持的 Abaqus 版本，请参阅 Simcenter STAR-CCM+ 发行说明的外部软件包支持部分。

输入文件

将完成教程所需要的文件复制到工作目录。

如果运行协同仿真案例，为 Simcenter STAR-CCM+ 和 Abaqus 文件分别建立目录是很好的做法。

1.创建名称为 /mechanicalCoSimTutorial 的目录。

2.在此目录中，创建另外两个目录：/abaqus 和 /starccm+。

已为本教程准备了两个文件：Abaqus 输入文件和流体 Simcenter STAR-CCM+ 模拟。将这些文件复制到 /abaqus 和 /starccm+ 目录下：

3.导航至已下载教程文件的 couplingWithCaeCodes 文件夹。

4.将文件 plate.inp 复制到 /mechanicalCoSimTutorial/abaqus/ 目录。

5.将文件 plate-cosim.sim 复制到 /mechanicalCoSimTutorial/starccm+/ 目录。

查看 Abaqus 模型

查看 Abaqus 模型以及协同仿真设置的特征。

本教程的 Abaqus 模型是固定在底部的平板。网格和平板尺寸如下图所示。

网格	C3D8R 元素
材料特性	杨氏模量：38.4 MPa
	泊松系数：0.3
	密度：4096 kg/m3
假设/模型	弹性材料
	瑞利阻尼的 β 系数为 1e-6
	单个动态负载步骤
最大时间步	5.0e-4 s
总仿真时间	0.125 s

协同仿真设置

为本教程所提供的 Abaqus 输入文件 plate.inp 已包含协同仿真的必要命令。协同仿真在加载步骤中执行。协同仿真命令包含在步骤定义中。您可以自行选择一个文本编辑器打开 plate.inp 文件以查看协同仿真命令。

**
** STEP: Step-1
**

文件的该部分，上至 *结束步骤包含 Abaqus 作业的历史数据，包括求解器参数、边界条件、负荷和输出。

**
*Output, history, variable=PRESELECT
**
*CO-SIMULATION, NAME=Plate, PROGRAM=MULTIPHYSICS

在步骤定义中，加粗行指定在该步骤执行协同仿真，并为协同仿真任取名称 Plate。

**
*Output, history, variable=PRESELECT
*CO-SIMULATION, NAME=Plate, PROGRAM=MULTIPHYSICS
*CO-SIMULATION REGION, TYPE=SURFACE, EXPORT
ASSEMBLY_FSI_INTERFACE, U
*CO-SIMULATION REGION, TYPE=SURFACE, IMPORT
ASSEMBLY_FSI_INTERFACE, CF

加粗的一组命令确定交界平面以及在协同仿真过程中导入和导出的场。

第一行确定模型中的一个区域，在本例中确定的是表面，数据由此表面导出。在下一行，此表面通过下列方式来确认：输入定义该表面时所使用的组件的名称（即 Assembly），然后输入表面的名称 FSI_Interface。表面识别符的句法是 <组件名称>_<表面名称>。将导出场定义为节点位移 U。接下来的两行使用相同的格式来指定将拉力矢量 CF 导入 Abaqus。