abaqus和Fluent的流固耦合模拟

abaqus和Fluent的流固耦合模拟
abaqus和Fluent的流固耦合模拟

耦合模拟

为耦合模拟ABAQUS需做如下工作:

l定义耦合步

l定义耦合区域

l定义耦合区域需要交换的物理量

以上每一步骤将在下面详细叙述

定义耦合步

ABAQUS耦合模拟界面是和存在的ABAQUS程序联合使用的。在你想定义的耦合步中,无论耦合情况如何,你必须先有效的载荷和边界条件。然后你再说明需要耦合的是这步,其中的一些量需要和三方软件进行数据交换。如下的一些过程ABAQUS是可以进行耦合分析的:

l准静态应力分析

l直接积分的隐式动态分析

l显式动态分析

l无耦合的热传导分析

l全积分热应力分析

与MPCCI server 数据交流始于耦合步,终于耦合步。

由于ABAQUS和其它三方软件在耦合分析过程中是实时的进行数据交换以及启动和终止三方程序,你可以在一个工作项目中只定义一个耦合步。

输入文件格式为:*CO-SIMULATION

定义接触区域

接触区域是系统之间的连接区域。这个表面对于ABAQUS而言必须是单元类型的面,任何对于MPCCI支持的单元类型均可以用于耦合步。而只有如下单元类型可以定义为接触区域,如表7.9.2-1

定义耦合区域的交换量

对于每个耦合区域你必须指定ABAQUS和其它三方软件进行交换的物理量,表7.9.2-2列出了可以用于交换和选择的物理量

输入输出的物理量的选择取决于分析的类型,如表7.9.2-3所示

输入文件的格式为:

*CO-SIMULA TION,IMPORT

surface_A,quantity_I1,quantity_I2,…

surface_B,quatity_I3

*CO-SIMULA TION,EXPORT

surface_A,quantity_E1

surface_B,quantity_E2

当前节点坐标和位移

因为在CFD代码中流体形状可以变化,不保持初始几何构型,所以在流固耦合(FSI)中选择当前节点坐标(COORD),而不是选择节点位移(U)。

不管是做小变形还是大变形,COORD的定义是当前节点坐标。这个定义和ABAQUS传统的定义有很大的区别,因为在传统的定义中再小变形情况下,初始坐标师部需要更新的。

位移再整体坐标系下总是交换的。如果对一个节点存在一个局部的变化,ABAQUS再将它们传给MPCCI server之前,先转化为整体坐标系下的位移。

集中力合法向压力

在流固耦合中,如果粘性剪切力需要导入,那么输入、输出的物理量就应该是集中节点力(CF)而不是法向压力(PRESS).在ABAQUS/Standard的耦合步中,从上一步到当前时间步,集中力和法向压力由一个斜坡的过渡。而在ABAQUS/Explicit的耦合步中,将保持常值,没有斜坡!

集中力总是在整体坐标系下进行数据交换的,如果一个节点存在一个局部坐标系,ABAQUS首先将集中力转化到局部坐标系下,然后再施加到结构上。

记中法向力可以在ABAQUS后处理中可视化。

热流和薄膜性质

对于流入一个表面热流分布可以用表面热流(HFL)这个概念。用薄膜性质(FILM)来模拟对流

这里q是进入表面的热流量,h是薄膜系数,为流体或者环境温度。薄膜系数可以通过热流和流体温度(FLUENT计算得到)以及壁温度(ABAQUS计算得到)来计算,如下

薄膜系数和流体温度传到ABAQUS内部,并保持为常值。当流体和壁的温度一致时,任意小的一个热传导系数付给ABAQUS。第一步耦合计算中,为了得到合理的薄膜常数,你应当保证在ABAQUS中壁温已经合适的初始化了,以及对初始的流体的温度场也有很好的估计。ABAQUS把初始化的壁温传给第三方软件。单位体系

ABAQUS对模型分析,没有特殊的单位要求。然而在耦合模拟中,ABAQUS运用的单位制必须和第三方

软件用的单位制相一致。MPCCI结构文件提供了转换单位的功能(详细见MPCCI用户手册)。

模型的维数和坐标系

向量在ABAQUS中可以转移,第一、第二、第三分别代表沿X,Y,Z方向的分量(对三维模型而言)。

需要注意的一点是,在ABAQUS中轴对称模型,对称轴式Y,而这个规则在三方软件中就未必是这么定义的。MPCCI提供了一个绘图工具来转换不同坐标系中计算结果。

限制

如果ABAQUS正在接受第三方软件的载荷值,那么耦合步就不能重新开算。这样载荷就不是作用在增量步的开始,这样会导致计算结果的发散。

对于壳单元这样的双面单元,必须分开指定。

abaqus与fluent流固耦合

基于MPCCI的流固耦合成功案例 基于MPCCI的流固耦合成功案例 (一)机翼气动弹性分析 1 问题陈述 机翼绕流问题是流固耦合中的经典问题。以前由于缺乏考虑流固耦合的软件,传统的分析方法是将机翼视为刚体,不考虑其弹性变形,通过CFD软件来计算机翼附近的流场。这个强硬的假设很难准确的描述流场的实际情况。更无法预测机翼的振动。MPCCI是基于代码耦合的并行计算接口,它可以同时调用结构和流体的软件来实现流固耦合。我们通过MPCCI,能很好的预测真实情况下的机翼绕流问题。采用ABAQUS结构分析软件来求解结构在流畅作用下的变形和应力分布,通过Fluent软件来计算由于固体运动和变形对整个流场的影响。 2 模拟过程分析顺序 MpCCI的图形用户界面可以方便的读入结构和流体的输入文件。后台调用ABAQUS和FLUENT。在MPCCI耦合面板中选择耦合面,然后选择在相应耦合面上流体和固体需要交换的量。启动MpCCI进行耦合。 3 边界条件设置

图1 无人机模型和流体计算模型 结构部分单个机翼跨度在1.5m左右,厚度为0.1m左右。边界条件为机翼端部的固定,三个方向的位移完全固定,另一端完全自由。在固体中除了固定端的面外,其他三个面为耦合面。流体部分采用四面体网格,采用理想气体作为密度模型。流体的入口和出口以及对称性边界条件如下图所示。 图2 固体有限元模型 4 计算方法的选择 通过结合ABAQUS和FLUENT,使用MPCCI计算流固耦合。在本例中,固体在流场作用下产生很大的变形和运动。在耦合区域,固体结构部分计算耦合面上的节点位移,通过MPCCI传输给FLUENT的耦合界面,FLUENT 计算出耦合区域上的节点力载荷,然后通过MPCCI传给结构软件ABAQUS。在MPCCI的耦合面板中选择的耦合面如图所示,交换量为:节点位移、相对受力。采用ABAQUS中的STANDARD算法,时间增量步长为0.1毫秒。 5 计算结论 通过MPCCI结合ABAQUS和FLUENT,成功地计算在几何非线性条件下的气动弹性问题,得到了整个流体区域的流场分布以及结构的动态响应历程。

基于MpCCI的Abaqus和Fluent流固耦合案例1

CAE联盟论坛精品讲座系列 基于MpCCI的Abaqus和Fluent流固耦合案例 主讲人:mafuyin CAE联盟论坛总监 摘要:通过MpCCI流固耦合接口程序,对某薄壁管道流动中的传热过程进行了Abaqus和Fluent相结合的流固耦合仿真分析。信息介绍了从建模、设置到求解计算和后处理的全过程,对相关研究人员具有参考意义。 1 分析模型 用三维建模软件solidworks建立了一个管径为1m的弯管,结构尺寸如图1a所示,管的结构如图1b所示,流体的模型如图1c所示。值得注意的是,由于拓扑特征的原因,这样的管壁模型无法通过对圆环扫略直接生成,而需先通过对大圆的扫略生成实心的模型(类似于流体模型),然后进行抽壳得到管壁的模型。用同样的方法对大圆半径减去管壁厚度的圆进行扫略得到流体模型。 a. 尺寸关系 b. 管壁结构 c. 流体模型 图1. 几何模型示意图 图2. 流固耦合传热分析模型示意图 内壁面(耦合面) 速度入口 v=6m/s; T in=600K 外壁面 压力出口 P=0Pa;T out=300K

由于管壁结构和流体的热学行为不同,传热系数等都不一样,所以属于典型的流固耦合传热问题,热学模型如图2所示。即管的一端为流体速度入口,一端为压力出口,给定流体外壁面一个初始温度600K,流体入口速度为6m/s,温度为600K,出口相对大气压力为0Pa,出口温度为300K。需要求解流体和管壁的温度场分布情况。 2 流体模型 将图1c的流体模型以Step格式导入Fluent软件通常使用的前处理器Gambit中,如图3a所示。设置求解器为,然后划分体网格,网格尺寸为100mm,类型为六面体单元,一共生成4895个体单元,网格如图3b所示。 a. 导入Gambit软件中的流体模型 b. 流场的网格模型 图3. 流体模型及网格示意图 进行网格划分后,需定义边界条件,在Gambit软件中先分别定义速度入口(VELOCITY_INLET)、压力出口(PRESSURE_OUTLET)和壁面(Wall)三组边界条件,具体参数设置在Fluent软件中进行。然后定义流体属性,名称定义为air,类型为Fluid。这些定义的目的是能够在Fluent软件中识别出这些特征,具体类型和参数都可以在Fluent软件中进行设置和修改。定义完后点击【Export】,选择【Mesh】,选择路径和文件名称并进行输出。 打开Fluent6.3.26或以上的版本,选择3D求解器,点击【File】→【Read】→【Case】,然后选择Gambit中输出的msh文件,即可将网格文件读入Fluent 软件中。读入模型后,进行求解参数和条件的设置。

基于MpCCI的Abaqus和Fluent流固耦合案例

基于MpCCI 的Abaqus 和Fluent 流固耦合案例 mafuyin 摘要:通过MpCCI 流固耦合接口程序,对某薄壁管道流动中的传热过程进行了Abaqus 和Fluent 相结合的流固耦合仿真分析。信息介绍了从建模、设置到求解计算和后处理的全过程,对相关研究人员具有参考意义。 1 分析模型 用三维建模软件solidworks 建立了一个管径为1m 的弯管,结构尺寸如图1a 所示,管的结构如图1b 所示,流体的模型如图1c 所示。值得注意的是,由于拓扑特征的原因,这样的管壁模型无法通过对圆环扫略直接生成,而需先通过对大圆的扫略生成实心的模型(类似于流体模型),然后进行抽壳得到管壁的模型。用同样的方法对大圆半径减去管壁厚度的圆进行扫略得到流体模型。 a. 尺寸关系 b. 管壁结构 c. 流体模型 图1. 几何模型示意图 图2. 流固耦合传热分析模型示意图 内壁面(耦合面) 速度入口 v=6m/s; T in =600K 外壁面 压力出口 P=0Pa ;T out =300K

由于管壁结构和流体的热学行为不同,传热系数等都不一样,所以属于典型的流固耦合传热问题,热学模型如图2所示。即管的一端为流体速度入口,一端为压力出口,给定流体外壁面一个初始温度600K,流体入口速度为6m/s,温度为600K,出口相对大气压力为0Pa,出口温度为300K。需要求解流体和管壁的温度场分布情况。 2 流体模型 将图1c的流体模型以Step格式导入Fluent软件通常使用的前处理器Gambit 中,如图3a所示。设置求解器为,然后划分体网格,网格尺寸为100mm,类型为六面体单元,一共生成4895个体单元,网格如图3b所示。 a. 导入Gambit软件中的流体模型 b. 流场的网格模型 图3. 流体模型及网格示意图 进行网格划分后,需定义边界条件,在Gambit软件中先分别定义速度入口(VELOCITY_INLET)、压力出口(PRESSURE_OUTLET)和壁面(Wall)三组边界条件,具体参数设置在Fluent软件中进行。然后定义流体属性,名称定义为air,类型为Fluid。这些定义的目的是能够在Fluent软件中识别出这些特征,具体类型和参数都可以在Fluent软件中进行设置和修改。定义完后点击【Export】,选择【Mesh】,选择路径和文件名称并进行输出。 打开Fluent6.3.26或以上的版本,选择3D求解器,点击【File】→【Read】→【Case】,然后选择Gambit中输出的msh文件,即可将网格文件读入Fluent 软件中。读入模型后,进行求解参数和条件的设置。 (1)模型缩放:为了便于分析结果数据特征,统一采用国际单位制进行仿真,

血管流固耦合分析实例

Ansys14 workbench血管流固耦合实例 根据收集的一些资料,进行学习后,试着做了这个ansys14workbench的血管流固耦合模拟,感觉能够耦合上,仅是熟悉流固耦合分析过程,不一定正确,仅供参考,希望大家多讨论。谢谢! 1、先在proe5中建立血管与血液流体区的模型(两者装配起来),或者直接在workbench中建模。 图1 模型图 2、新建工程。在workbench中toolbox中选custom system,双击FSI: FluidFlow(fluent)->static structure. 图2 计算工程 3、修改engineering data,因为系统缺省材料是钢,需要构建血管材料,如图3所示。先复制steel,而后修改密度1150kg/m3,杨氏模量4.5e8Pa,泊松比0.3,重新命名,最后在主菜单中点击“update project”保存.

图3 修改工程材料 4、模型导入,进入gemetry模块,import外部模型文件。 图4 模型导入图 5、进入FLUENT网格划分。 在workbench工程视图中的Mesh上点击右键,选择Edit…,如图5所示,进入网格划分meshing界面,如图6所示。我们这里需要去掉血管部分,只保留血液几何。

图5 进入网格划分

图6 禁用血管模型 6、设置网格方法。 默认是采用ICEM CFD进行网格划分,设置方式如图7所示,截面圆弧边分为12份,纵截面的边均分为10份,网格结果如图8所示。另外在这个界面中要设置边界的几何面,如inlet、outlet、symmetry 图7 设置网格划分方式 图8 最终出网格

fluent单项流固耦合

流固耦合(Fluid-solid interaction,FSI)计算,通常用于考虑流体与固体间存在强烈的相互作用时,对流体流场与固体应力应变的考察。FSI计算按数据传递方式可分两类:单向耦合与双向耦合。所谓单向耦合,主要是指数据只从流体计算传递压力到固体,或者只从固体计算传递网格节点位移到流体。双向耦合则在每一时刻都同时向对方发送相应的物理量(流体计算发送压力数据,固体计算发送位移数据)。 ANSYS Workbench中可以利用Fluent与DS进行单向流固耦合计算。我们这里来举一个最简单的单向耦合例子:风吹挡板。我们假定挡板位移可忽略不计,固体变形对流场影响可以忽略,所考虑的是流体压力作用在固体上,固体的应力分布。当然这里的压力可以换成温度等其他物理量。 1、新建工程。注意是从Fluent –> Static Structure。连接图如1所示。 图1 计算工程关 系图2 进入DM建模 2、进入Fluent中的DM进行模型创建,如图2所示。 流固耦合计算中的几何模型与单纯的流体模型或固体模型不同,它要求同时具有流体和固体模型,而且流体计算中只能有流体模型,固体计算中只能有固体模型。建好后的模型如图3,4,5所示。由于固体模型需要从这里导入,所以我们保留固体与流体模型。

图3 实体模型 图4 固体模型

图5 流体模型 3、进入FLUENT网格设置。 在FLUENT工程视图中的Mesh上点击右键,选择Edit…,如图6所示,进入网格划分meshing界面,如图7所示。我们这里需要去掉固体部分,只保留流体几何。 图6 进入网格划 分图7 禁用固体模型

(整理)FLUENT14双向流固耦合案例.

说明:本例只应用于FLUENT14.0以上版本。 ANSYS 14.0是2011年底新推出的版本,在该版本中,加入了一个新的模块System Coupling,目前只能用于fluent与ansys mechanical的双向流固耦合计算。官方文档中有介绍说以后会逐渐添加对其它求解器的支持,不过这不重要,重要的是现在FLUENT终于可以不用借助第三方软件进行双向流固耦合计算了,个人认为这是新版本一个不小的改进。 模块及数据传递方式如下图所示。 一、几何准备 流固耦合计算的模型准备与单独的流体计算不同,它需要同时创建流体模型与固体模型。在geometry模块中同时创建流体模型与固体模型。到后面流体模型或固体模块中再进行模型禁用处理。 模型中的尺寸:v1:32mm,h2:120mm,h5:60mm,h3:3mm,v4:15mm。

由于流体计算中需要进行动网格设置,因此推荐使用四面体网格。当然如果挡板刚度很大网格变形很小时,可以使用六面体网格,划分六面体网格可以先将几何进行slice切割。这里对流体区域网格划分六面体网格,固体域同样划分六面体网格。 二、流体部分设置 1、网格划分 双击B3单元格,进入meshing模块进行网格划分。禁用固体部分几何。设定各相关部分的尺寸,由于固体区域几何较为整齐,因此在切割后只需设定一个全局尺寸即可划分全六面体网格。这里设定全局尺寸为1mm。划分网格后如下图所示。 2、进行边界命名,以方便在fluent中进行边界条件设置 设置左侧面为速度进口velocity inlet,右侧面为自由出流outflow,上侧面为壁面边界wall_top,正对的两侧面为壁面边界wall_side1与wall_side2(这两个边界在动网格设定中为变形域),设定与固体交界面为壁面边界(该边界在动网格中设定为system coupling类型)。 操作方式:选择对应的表面,点击右键,选择菜单create named selection,然后输入相应的边界名称。注意:FLUENT会自动检测输入的名称以使用对应的边界类型,当然用户也可以在fluent进行类型更改。完成后的树形菜单如下图所示。

基于LSDYNA及FLUENT的板壳结构流固耦合分析

基于 LS-DYNA 及 FLUENT 的板壳结构流-固耦合分析
汪丽军 北京航空航天大学,交通科学与工程学院 100191
[摘 要]: 本文采用 ANSYS 显示动力分析模块 LS-DYNA 及流场分析模块 FLUENT,对水下的板壳 结构运动及其界面的流-固耦合现象进行了仿真分析。流场计算得到的界面压强数据以外载荷 的形式施加于结构表面,使其产生位移及变形;同时,结构的变化又进一步影响了流场的分 布。通过往复的双向耦合迭代,得到了板壳结构的动力学响应以及流场的分布情况。仿真结 果与试验结果的对比表明,此方法适用于解决兼有大位移及较大变形特征的流-固耦合问题。 [关键词]: 板壳结构 流-固耦合 有限元方法 ANSYS
Analysis of Fluid-Structure Interaction for Plate/Shell Structure Based on LS-DYNA and FLUENT
Wang Lijun School of Transportation Science & Engineering, Beihang University 100191
Abstract: In this paper,the movement of plate under water and the fluid-structure interaction(FSI) is simulated numerically by combining explicit dynamic solver LS-DYNA and computational fluid dynamics solver FLUENT in ANSYS. The pressure obtained from the calculation of flow field are applied as external loads on the surface of the plate, then the structural deformation and displacement can be calculated as well, which will affect the shape and pressure distribution of the flow field reversely. After sequential coupling iterations the dynamic response of the structure and flow field distribution are obtained consequently. By comparing numerical and experimental results it is proved that this proposed coupling method is suitable for solving such a kind of FSI problems considering both large displacement and comparatively large deformation. Keyword: Plate/shell structure, Fluid-Structure Interaction, Finite element method,ANSYS
1
前言
在自然界中,流-固耦合现象广泛存在于航空、航天、汽车、水利、石油、化工、海洋 以及生物等领域。很多实际问题中流体载荷对于结构的影响不可忽略;同时,结构的位移 和变形也会对流场的分布产生重要影响。例如各种水下运动机构都需要考虑这种现象。

(完整word版)ABAQUS实例分析

《现代机械设计方法》课程结业论文 ( 2011 级) 题目:ABAQUS实例分析 学生姓名 XXXX 学号 XXXXX 专业机械工程 学院名称机电工程与自动化学院 指导老师 XX 2013年 5 月8 日

目录 第一章Abaqus简介 (1) 一、Abaqus总体介绍 (1) 二、Abaqus基本使用方法 (2) 1.2.1 Abaqus分析步骤 (2) 1.2.2 Abaqus/CAE界面 (3) 1.2.3 Abaqus/CAE的功能模块 (3) 第二章基于Abaqus的通孔端盖分析实例 (4) 一、工作任务的明确 (6) 二、具体步骤 (6) 2.2.1 启动Abaqus/CAE (4) 2.2.2 导入零件 (5) 2.2.3 创建材料和截面属性 (6) 2.2.4 定义装配件 (7) 2.2.5 定义接触和绑定约束(tie) (10) 2.2.6 定义分析步 (14) 2.2.7 划分网格 (15) 2.2.8 施加载荷 (19) 2.2.9 定义边界条件 (20) 2.2.10 提交分析作业 (21) 2.2.11 后处理 (22) 第三章课程学习心得与作业体会 (23)

第一章: Abaqus简介 一、Abaqus总体介绍 Abaqus是功能强大的有限元分析软件,可以分析复杂的固体力学和结构力学系统,模拟非常庞大的模型,处理高度非线性问题。Abaqus不但可以做单一零件的力学和多物理场的分析,同时还可以完成系统级的分析和研究。 Abaqus使用起来十分简便,可以很容易的为复杂问题建立模型。Abaqus具备十分丰富的单元库,可以模拟任意几何形状,其丰富的材料模型库可以模拟大多数典型工程材料的性能,包括金属、橡胶、聚合物、复合材料、钢筋混泥土、可压缩的弹性泡沫以及地质材料(例如土壤、岩石)等。 Abaqus主要具有以下分析功能: 1.静态应力/位移分析 2.动态分析 3.非线性动态应力/位移分析 4.粘弹性/粘塑性响应分析 5.热传导分析 6.退火成形过程分析 7.质量扩散分析 8.准静态分析 9.耦合分析 10.海洋工程结构分析 11.瞬态温度/位移耦合分析 12.疲劳分析 13.水下冲击分析 14.设计灵敏度分析 二、Abaqus基本使用方法 1.2.1 Abaqus分析步骤 有限元分析包括以下三个步骤: 1.前处理(Abaqus/CAE):在前期处理阶段需要定义物理问题的模型,并生 成一个Abaqus输入文件。提交给Abaqus/Standard或 Abaqus/Explicit。 2.分析计算(Abaqus/Standard或Abaqus/Explicit):在分析计算阶段, 使用Abaqus/Standard或Abaqus/Explicit求解输入文件中所定义的

双向流固耦合实例Fluent与structure

双向流固耦合实例( Fluent 与 structure) 说明:本例只应用于FLUENT14.0 以上版本。 ANSYS 14.0 是2011 年底新推出的版本,在该版本中,加入了一个新的模块System Coupling ,目前只能用于fluent 与ansys mechanical 的双向流固耦合计算。官方文档中有介绍说以后会逐渐添加对其它求解器的支持,不过这不重要,重要的是现在FLUENT 终于可以不用借助第三方软件进行双向流固耦合计算了,个人认为这是新版本一个不小的改进。 模块及数据传递方式如下图所示。 一、几何准备流固耦合计算的模型准备与单独的流体计算不同,它需要同时创建流体模型与固体模型。在geometry 模块中同时创建流体模型与固体模型。到后面流体模型或固体模块中再进行模型禁用处理。 模型中的尺寸:v1:32mm ,h2:120mm ,h5:60mm ,h3:3mm ,v4:15mm 。由于流体计算中需要进行动网格设置,因此推荐使用四面体网格。当然如果挡板刚度很大网格变形很小时,可以使用六面体网格,划分六面体网格可以先将几何进行slice 切割。这里对流体区域网格划分六面体网格,固体域同样划分六面体网格。 二、流体部分设置 1、网格划分 双击B3 单元格,进入meshing 模块进行网格划分。禁用固体部分几何。设定各相关部分的尺寸,由于固体区域几何较为整齐,因此在切割后只需设定一个全局尺寸即可划分全六面体网格。这里设定全局尺寸为1mm 。划分网格后如下图所示。

2 、进行边界命名,以方便在 fluent 中进行边界条件设置 变形域),设定与固体交界面为壁面边界(该边界在动网格中设定为 操作方式:选择对应的表面,点击右键,选择菜单 create named selection ,然后输入相 应的边界名称。 注意: FLUENT 会自动检测输入的名称以使用对应的边界类型, 在 fluent 进行类型更改。完成后的树形菜单如下图所示。 本部分操作完毕后,关闭 meshing 模块。返回工程面板。 3 、进入 fluent 设置 FLUENT 主要进行动网格设置。其它设置与单独进行 FLUENT 仿真完全一致。 设置使用瞬态计算,使用 K-Epsilon 湍流模型。 这里的动网格主要使用弹簧光顺处理(由于使用的是六面体网格且运动不规律),需要使 用 TUI 命令打开光顺对六面体网格的支持。使用命令 /define/dynamic-mesh/controls/smoothing-parameters 。 动态层技术与网格重构方法在六面体网格中失效。因此,建议使用四面体网格。我们这里 由于变形小,所以只使用光顺方法即可满足要求。 点击 Dynamic mesh 进入动网格设置面板。如下图所示,激活动网格模型。 设置左侧面为速度进口 velocity inlet ,右侧面为自由出流 outflow ,上侧面为壁面边界 wall_top ,正对的两侧面为壁面边界 wall_side1 与 wall_side2 这两个边界在动网格设定中为 system coupling 类型) 当然用户也可以

ABAQUS顺序热力耦合分析实例

ABAQUS顺序热力耦合分析实例此实例中需要确定一个冷却栅管的温度场分布。温度场的求解采用稳态热分析,在此之后还将进行热应力分析来求出冷却栅管在温度作用下产生的位移和应力分布。由于冷却栅管比较长,并且是轴对称结构,根据上述特点,可以简化有限元分析模型。此实例中使用国际单位制。 1、part中创建轴对称可变形壳体,大致尺寸为1,通过creat line创建一个封闭曲线(0.127,0) (0.304,0)(0.304,0.006)(0.152,0.006)(0.152,0.031)(0.127,0.031)(0.127,0) 使用creat Fillet功能对模型倒角处设置0.005的倒圆角。倒角后,模型并未改变,需要在模型树中,part下的Features右键,Regenerate,最终模型如下图所示。 2、在材料模块中定义密度7800,弹性模量1.93E11,泊松比0.3。所不同的是,热分析还需 要指定热传导系数以及比热。在Thermal里输入参数,热铲刀系数25.96,比热451。 3、创建截面属性以及装备部件,和普通的静力分析设置一样。 4、Step有所不同,分析类型仍为通用分析步,下面要更改为Heat Transfer。在Edit Step窗 口中,使用默认的瞬态分析(Transient),时长设置为3s。切换到Incrementatin进行相应的设置,如下图。

5、Load模块中,设置左边温度为100度,右边及上边温度为20度。Creat BC,类型选择 Other>Temperature。在纯粹的热传导分析方程中,没有位移项,因此不会发生刚体位移,这里也就不需要设置位移边界条件。 6、接下来划分网格,种子尺寸给0.005,单元类型需要在单元族中选择专门用来热分析的 Heat Transfer,查看下面确保使用的单元为DCAX4。使用结构化的全四边形网格划分方法。 7、到此,热分析的设置已经完成,可以提交计算,完成后,查看变量NT11即为节点温度。

ABAQUS实例分析论文

目录 第一章Abaqus简介 (1) 一、Abaqus总体介绍 (1) 二、Abaqus基本使用方法 (2) 1.2.1 Abaqus分析步骤 (2) 1.2.2 Abaqus/CAE界面 (3) 1.2.3 Abaqus/CAE的功能模块 (3) 第二章基于Abaqus的通孔端盖分析实例 (4) 一、工作任务的明确 (5) 二、具体步骤 (5) 2.2.1 启动Abaqus/CAE (4) 2.2.2 导入零件 (5) 2.2.3 创建材料和截面属性 (6) 2.2.4 定义装配件 (7) 2.2.5 定义接触和绑定约束(tie) (10) 2.2.6 定义分析步 (14) 2.2.7 划分网格 (15) 2.2.8 施加载荷 (19) 2.2.9 定义边界条件 (20) 2.2.10 提交分析作业 (21) 2.2.11 后处理 (22) 第三章课程学习心得与作业体会 (22)

第一章: Abaqus简介 一、Abaqus总体介绍 Abaqus是功能强大的有限元分析软件,可以分析复杂的固体力学和结构力学系统,模拟非常庞大的模型,处理高度非线性问题。Abaqus不但可以做单一零件的力学和多物理场的分析,同时还可以完成系统级的分析和研究。 Abaqus使用起来十分简便,可以很容易的为复杂问题建立模型。Abaqus具备十分丰富的单元库,可以模拟任意几何形状,其丰富的材料模型库可以模拟大多数典型工程材料的性能,包括金属、橡胶、聚合物、复合材料、钢筋混泥土、可压缩的弹性泡沫以及地质材料(例如土壤、岩石)等。 Abaqus主要具有以下分析功能: 1.静态应力/位移分析 2.动态分析 3.非线性动态应力/位移分析 4.粘弹性/粘塑性响应分析 5.热传导分析 6.退火成形过程分析 7.质量扩散分析 8.准静态分析 9.耦合分析 10.海洋工程结构分析 11.瞬态温度/位移耦合分析 12.疲劳分析 13.水下冲击分析 14.设计灵敏度分析 二、Abaqus基本使用方法 1.2.1 Abaqus分析步骤 有限元分析包括以下三个步骤: 1.前处理(Abaqus/CAE):在前期处理阶段需要定义物理问题的模型,并生

最新fluent流固耦合传热设置问题

FLUENT流固耦合传热设置问题 看到很多网友对于fluent里模拟流固耦合传热(同时有对流和导热)有很多疑问,下面说说我的解决方法。 1,首先要分清你的问题是否是流固耦合传热。 (1)如果你的传热问题只是流体与固体壁面的传热,不涉及到固体壁面内部的导热,那么这就是一个对流传热问题,不是流固耦合传热问题, 这时候你只需要设置壁面的对流换热系数即可。如下图 注意右边这几个参数的含义:从上往下依次为:壁面外部的对流传热系数;外部流体温度;壁面厚度;壁面单位体积发热率。 这里没有内部流体的对流传热设置,因为fluent会根据流体温度以及壁面温度,利用能量守恒,自动计算内壁流体与壁面的对流换热情况。 (2)流固耦合传热问题。在建模的时候你应该定义两个区域,流体区域和固体区域,并且在切割区域的时候,你应该选中connect,如下图所 示 边界条件设置:交界面为wall。在导入fluent以后,fluent就会自动生成wall-shadow。这样在流固交界面上就生成了一对耦合的面,如下图所示,

。 2,耦合传热设置问题 (1)首先就是求解器的设置问题,应该选择耦合求解器,虽然计算速度会慢一些,但是这更符合实际情况,更容易收敛,误差更小。如果是非 稳态过程还应选择unsteady。如下图所示 (2)交界面设置问题,这个是关键。不用过多的设置只需要选择coupled。 这样fluent就会自动计算耦合面的传热问题。如下图所示

(3)当然还要选择能量方程。其他诸如湍流模型、材料设置、进出口条件等等,需要你根据实际情况设定,这里不再雷述。1.在国际单位制中,电荷的单位是 A. 伏特 B. 安培 C. 库仑 D.瓦特 2.小明家装修房屋需要购买导线,关于导线种类的选择,最恰当的是: A.强度大的铁丝B.细小价格较便宜的铝丝 C.粗一点的铜丝D.性能稳定的镍铬合金丝 3.小明在研究通过导体的电流时,根据测量数据绘制出如图 所示的I-U图像。对此作出的判断中,错误 ..的是: A.通过R1的电流与它两端所加电压成正比 B.通过R2的电流与它两端所加电压不成正比 C.将它们串联接入到同一电路中时,通过R1的电流较小 D.将它们并联连接到两端电压为1.5V的电路中时,通过 干路的电流大约是0.46A 4.小灯泡L上标有“2.5V”字样,它的电阻随它两端电压变化的图像如图甲所示。将小灯泡L和电阻R0接入图乙所示的电路中,电源电压为6V,且保持不变。当开 关S闭合时,小灯泡L恰好能正常发光。 下列说法正确的是: A.开关S断开时,小灯泡L的电阻为0Ω B.开关S闭合时,小灯泡L的电阻为8Ω C.小灯泡L的额定功率为0.5W D.电阻R0的阻值为14Ω 5.假设导体没有电阻,当用电器通电时,下列说法正确的是() A.白炽灯仍然能发光B.电动机仍然能转动 C.电饭锅仍然能煮饭D.电熨斗仍然能熨衣服 6.在图8所示电路中,闭合开关S后,在滑片P 向右滑动过程中,各电表示数变化正确的 是() A.A1、A3示数不变,A2、V示数变小 B.A1、V 示数不变,A2、A 3示数变大R1 R2

ABAQUS流固耦合

当只进行渗流计算时: 1.由于Abaqus中缺乏非耦合的孔压单元,这时可采用耦合单元,但要约束住所有位移的 自由度。 2.渗流材料参数选择。在CAE中都是在(Material-creat-other-pore fluid)选项中。 1)Gel:定义凝胶微粒吸湿膨胀的发育过程,这在一般的岩土分析中应用不多。 2)Moisture swelling:定义由于吸湿饱和所引起的固体骨架体积膨胀(或负吸力引起的 骨架收缩)。 3)(3)Permeability:定义饱和介质的渗透系数,该渗透系数可以在type选项中定义 为各向同性、正交各向异性和各向异性,并且可以根据Void Ratio定义为孔隙比的 函数。在Suboptions中选择Saturation Dependent 参数来指定与饱和度相关性系数 ks(s),缺省设置为ks=s3,而非饱和介质渗透系数k’=ksk 选择Velocity dependence 参数可以激活Forchheimer定律,缺省的是Darcy定律 4)Pore Fluid Expansion:定义固体颗粒与流体体积热变化效应。 5)Porous Bulk Moduli:定义固体颗粒与流体体积模量。 6)Sorption:定义负孔隙压力与饱和度之间的相关性。当type=Absorption时,定义吸 湿曲线,type=Exsorption时定义排水曲线。 3、载荷及边界条件 1)通过(Load-creat-step-fluid-surface pore fluid)选项定义沿着单元表面的外法线方向 的渗流速度vn,当考虑降雨影响时可采用此载荷 2)边界条件(Boundary condition-creat-other-pore pressure)选项定义孔压边界条件,此 时要先假定浸润面的位置,然后定义浸润面上的孔压为零,Abaqus会在后续的分 析计算中自动计算出浸润面的位置。Abaqus默认的是不透水边界。 3)当渗流自由面遇到临空的自由排水面时,需要定义一个特殊的边界条件。此时可 以通过在inp文件中加入*Flow或*Sflow来定义 4)初始条件的定义。初始条件中一般要定义以下几种:*initial condition,type=saturation 初始饱和度initial condition,type=pore pressure 初 始孔压initial condition,type=ratio 初始孔隙比当进行耦合分析时,基本步骤同 上,但要去掉除边界条件之外的约束,同时还要在边界上加上流体压力

abaqus和Fluent的流固耦合模拟

耦合模拟 为耦合模拟ABAQUS需做如下工作: l定义耦合步 l定义耦合区域 l定义耦合区域需要交换的物理量 以上每一步骤将在下面详细叙述 定义耦合步 ABAQUS耦合模拟界面是和存在的ABAQUS程序联合使用的。在你想定义的耦合步中,无论耦合情况如何,你必须先有效的载荷和边界条件。然后你再说明需要耦合的是这步,其中的一些量需要和三方软件进行数据交换。如下的一些过程ABAQUS是可以进行耦合分析的: l准静态应力分析 l直接积分的隐式动态分析 l显式动态分析 l无耦合的热传导分析 l全积分热应力分析 与MPCCI server 数据交流始于耦合步,终于耦合步。 由于ABAQUS和其它三方软件在耦合分析过程中是实时的进行数据交换以及启动和终止三方程序,你可以在一个工作项目中只定义一个耦合步。 输入文件格式为:*CO-SIMULATION 定义接触区域 接触区域是系统之间的连接区域。这个表面对于ABAQUS而言必须是单元类型的面,任何对于MPCCI支持的单元类型均可以用于耦合步。而只有如下单元类型可以定义为接触区域,如表7.9.2-1 定义耦合区域的交换量 对于每个耦合区域你必须指定ABAQUS和其它三方软件进行交换的物理量,表7.9.2-2列出了可以用于交换和选择的物理量

输入输出的物理量的选择取决于分析的类型,如表7.9.2-3所示 输入文件的格式为: *CO-SIMULA TION,IMPORT surface_A,quantity_I1,quantity_I2,… surface_B,quatity_I3 *CO-SIMULA TION,EXPORT surface_A,quantity_E1 surface_B,quantity_E2 当前节点坐标和位移 因为在CFD代码中流体形状可以变化,不保持初始几何构型,所以在流固耦合(FSI)中选择当前节点坐标(COORD),而不是选择节点位移(U)。 不管是做小变形还是大变形,COORD的定义是当前节点坐标。这个定义和ABAQUS传统的定义有很大的区别,因为在传统的定义中再小变形情况下,初始坐标师部需要更新的。 位移再整体坐标系下总是交换的。如果对一个节点存在一个局部的变化,ABAQUS再将它们传给MPCCI server之前,先转化为整体坐标系下的位移。 集中力合法向压力 在流固耦合中,如果粘性剪切力需要导入,那么输入、输出的物理量就应该是集中节点力(CF)而不是法向压力(PRESS).在ABAQUS/Standard的耦合步中,从上一步到当前时间步,集中力和法向压力由一个斜坡的过渡。而在ABAQUS/Explicit的耦合步中,将保持常值,没有斜坡! 集中力总是在整体坐标系下进行数据交换的,如果一个节点存在一个局部坐标系,ABAQUS首先将集中力转化到局部坐标系下,然后再施加到结构上。 记中法向力可以在ABAQUS后处理中可视化。 热流和薄膜性质 对于流入一个表面热流分布可以用表面热流(HFL)这个概念。用薄膜性质(FILM)来模拟对流 这里q是进入表面的热流量,h是薄膜系数,为流体或者环境温度。薄膜系数可以通过热流和流体温度(FLUENT计算得到)以及壁温度(ABAQUS计算得到)来计算,如下 薄膜系数和流体温度传到ABAQUS内部,并保持为常值。当流体和壁的温度一致时,任意小的一个热传导系数付给ABAQUS。第一步耦合计算中,为了得到合理的薄膜常数,你应当保证在ABAQUS中壁温已经合适的初始化了,以及对初始的流体的温度场也有很好的估计。ABAQUS把初始化的壁温传给第三方软件。单位体系 ABAQUS对模型分析,没有特殊的单位要求。然而在耦合模拟中,ABAQUS运用的单位制必须和第三方

双向流固耦合实例(Fluent与structure)

双向流固耦合实例(Fluent与structure) 说明:本例只应用于FLUENT14.0以上版本。 ANSYS 14.0是2011年底新推出的版本,在该版本中,加入了一个新的模块System Coupling,目前只能用于fluent与ansys mechanical的双向流固耦合计算。官方文档中有介绍说以后会逐渐添加对其它求解器的支持,不过这不重要,重要的是现在FLUENT终于可以不用借助第三方软件进行双向流固耦合计算了,个人认为这是新版本一个不小的改进。 模块及数据传递方式如下图所示。 一、几何准备 流固耦合计算的模型准备与单独的流体计算不同,它需要同时创建流体模型与固体模型。在geometry模块中同时创建流体模型与固体模型。到后面流体模型或固体模块中再进行模型禁用处理。 模型中的尺寸:v1:32mm,h2:120mm,h5:60mm,h3:3mm,v4:15mm。 由于流体计算中需要进行动网格设置,因此推荐使用四面体网格。当然如果挡板刚度很大网格变形很小时,可以使用六面体网格,划分六面体网格可以先将几何进行slice切割。这里对流体区域网格划分六面体网格,固体域同样划分六面体网格。 二、流体部分设置 1、网格划分 双击B3单元格,进入meshing模块进行网格划分。禁用固体部分几何。设定各相关部分的尺寸,由于固体区域几何较为整齐,因此在切割后只需设定一个全局尺寸即可划分全六面体网格。这里设定全局尺寸为1mm。划分网格后如下图所示。

2、进行边界命名,以方便在fluent中进行边界条件设置 设置左侧面为速度进口velocity inlet,右侧面为自由出流outflow,上侧面为壁面边界wall_top,正对的两侧面为壁面边界wall_side1与wall_side2(这两个边界在动网格设定中为变形域),设定与固体交界面为壁面边界(该边界在动网格中设定为system coupling类型)。 操作方式:选择对应的表面,点击右键,选择菜单create named selection,然后输入相应的边界名称。注意:FLUENT会自动检测输入的名称以使用对应的边界类型,当然用户也可以在fluent进行类型更改。完成后的树形菜单如下图所示。 本部分操作完毕后,关闭meshing模块。返回工程面板。 3、进入fluent设置 FLUENT主要进行动网格设置。其它设置与单独进行FLUENT仿真完全一致。 设置使用瞬态计算,使用K-Epsilon湍流模型。 这里的动网格主要使用弹簧光顺处理(由于使用的是六面体网格且运动不规律),需要使用TUI命令打开光顺对六面体网格的支持。使用命令 /define/dynamic-mesh/controls/smoothing-parameters。 动态层技术与网格重构方法在六面体网格中失效。因此,建议使用四面体网格。我们这里由于变形小,所以只使用光顺方法即可满足要求。 点击Dynamic mesh进入动网格设置面板。如下图所示,激活动网格模型。

ABAQUS金属切削实例

CAE联盟论坛精品讲座系列【二】 ABAQUS金属切削实例 主讲人:fuyun123CAE联盟论坛—ABAQUS版主 背景介绍: 切削过程是一个很复杂的工艺过程,它不但涉及到弹性力学、塑性力学、断裂力学,还有热力学、摩擦学等。同时切削质量受到刀具形状、切屑流动、温度分布、热流和刀具磨损等影响,切削表面的残余应力和残余应变严重影响了工件的精度和疲劳寿命。利用传统的解析方法,很难对切削机理进行定量的分析和研究。计算机技术的飞速发展使得利用有限元仿真方法来研究切削加工过程以及各种参数之间的关系成为可能。近年来,有限元方法在切削工艺中的应用表明,切削工艺和切屑形成的有限元模拟对了解切削机理,提高切削质量是很有帮助的。这种有限元仿真方法适合于分析弹塑性大变形问题,包括分析与温度相关的材料性能参数和很大的应变速率问题。ABAQUS作为有限元的通用软件,在处理这种高度非线性问题上体现了它独到的优势,目前国际上对切削问题的研究大都采用此软件,因此,下面针对ABAQUS的切削做一个入门的例子,希望初学者能够尽快入门,当然要把切削做好,不单单是一个例子能够解决问题的,随着深入的研究,你会发现有很多因素影响切削的仿真的顺利进行,这个需要自己去不断探索,在此本人权当抛砖引玉,希望各位切削的大神们能够积极探讨起来,让我们在切削仿真的探索上更加精确,更加完善。 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 切削参数:切削速度300m/min,切削厚度0.1mm,切削宽度1mm 尺寸参数:本例作为入门例子,为了简化问题,假定刀具为解析刚体,因为在切削过程中,一般我们更注重工件最终的切削质量,如应力场,温度场等,尤其是残余应力场,而如果是要进行刀具磨损或者涂层刀具失效的分析的话,那就要考虑建立刀具为变形体来进行分析了。工件就假定为一个长方形,刀具设置前角10°,后角6°,具体尺寸见INP文件。 下面将切削过程按照ABAQUS的模块分别进行叙述,并对注意的问题作出相应的解释。 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 建模:建模过程其实没有什么好注意的,对于复杂的模型,我一般用其他三维软件导入进来,注意导入的时候尽量将格式转化为IGES格式,同时要把一些不必要的东西去掉,比如一些尖角,圆角之类的,如果不是分析那个部位的应力集中的话就没必要导入它,如果导入,还要进行一些细化,大大降低了计算的效率。我一般做的是二维切削,模型相对比较简单,所以一般都是直接在ABAQUS中进行建模。由于此处为刚体,要在part里面建立刚体参考点,而且注意不要在装配模块建立参考点,因为有时候ABAQUS找不到装配模块相应的参考点。 1、工件

fluent流固耦合传热设置问题

f l u e n t流固耦合传热设置 问题 Prepared on 21 November 2021

FLUENT流固耦合传热设置问题 看到很多网友对于fluent里模拟流固耦合传热(同时有对流和导热)有很多疑问,下面说说我的解决方法。 1,首先要分清你的问题是否是流固耦合传热。 (1)如果你的传热问题只是流体与固体壁面的传热,不涉及到固体壁面内部的导热,那么这就是一个对流传热问题,不是流固耦合传热问 题,这时候你只需要设置壁面的对流换热系数即可。如下图 注意右边这几个参数的含义:从上往下依次为:壁面外部的对流传热系数;外部流体温度;壁面厚度;壁面单位体积发热率。 这里没有内部流体的对流传热设置,因为fluent会根据流体温度以及壁面温度,利用能量守恒,自动计算内壁流体与壁面的对流换热情况。 (2)流固耦合传热问题。在建模的时候你应该定义两个区域,流体区域和固体区域,并且在切割区域的时候,你应该选中connect,如下图 所示 边界条件设置:交界面为wall。在导入fluent以后,fluent就会自动生成wall-shadow。这样在流固交界面上就生成了一对耦合的面,如下图所示, 。 2,耦合传热设置问题 (1)首先就是求解器的设置问题,应该选择耦合求解器,虽然计算速度会慢一些,但是这更符合实际情况,更容易收敛,误差更小。如果 是非稳态过程还应选择unsteady。如下图所示 (2)交界面设置问题,这个是关键。不用过多的设置只需要选择coupled。这样fluent就会自动计算耦合面的传热问题。如下图所 示 (3)当然还要选择能量方程。其他诸如湍流模型、材料设置、进出口条件等等,需要你根据实际情况设定,这里不再雷述。

相关文档
最新文档