制造业信息化
仿真/建模/CAD/CAM/CAE/CAPP
MANUFACTURINGINFORMATIZATION基于ANSYSICEMCFD和CFX数值仿真技术许蕾,罗会信
)(武汉科技大学机械自动化学院,武汉430081
摘为了准确方便地计算和分析流体的传热和流动情况,文中以分析某浸入式水口浇钢温度场和流场为例介绍了要:
一种基于ANSYSICEMCFD和CFX数值仿真技术在连续铸钢过程中的应用。该方法是利用ANSYSICEMCFD中进行流体数值仿真的前处理操作,然后将生成的网格导入到CFX中施加边界条件进行求解计算,可以大大提高数值仿真效率,对CAE工程分析人员具有很大的借鉴意义。
CFD;CFX;数值仿真关键词:
中图分类号:TP391文献标识码:A文章编号:1002-2333(2008)12-0065-02
XULei,
LUOHui-xin
TheTechnologyofNumericalSimulationBasedonANSYSICEMCFDandCFXSoftware
(DepartmentofMachinery&Automation,WuhanUniversityofScience&Technology,Wuhan430081,China)
Abstract:Tocalculateandanalyzetheheattransferandflowtrulyandexpediently,takingtoanalyzethesubmergedentrynozzle'stemperaturefieldandflowfieldasanexample,thisarticleintroducesthetechnologyofnumericalsimulationinANSYSICEMCFDsoftwareapplyingtocontinuouscasting.DothepreprocessingoffluidnumericalsimulationinANSYSICEMCFDsoftwarefirst,andthenexportsomemeshestoCFXforsolution.Theengineerscanborrowsomeideasfromthismethod,whichcanenhancetheefficiencyofthenumericalsimulationgreatly.
Keywords:CFD;CFX;numericalsimulation
1前言ICEMCFD后的结果会有差异。模型复杂一点的IGES格式,导入后会有很多面的丢失,相对而言CATIA的文件格式比较方便。
在导入几何模型后,需要对模型进行清理。ICEMCFD的Geometry页面中的RepairGeometry菜单提供了大量的优化模型和清理模型的工具。对只导入面的模型自动生成了曲线和点,便于网格的划分。3.2
划分计算网格
ICEMCFD通过映射技术的六面体网格划分功能在拓扑空间进行网格划分,自动映射到物理空间,可在任意形状的模型中划分出六面体网格;也可以应用其自动划分网格功能来生成自由网格。要指出的是由于ICEMCFD自由划分网格的功能十分便捷,因此很多人为了方便直接采用tetra网格。然而对于流场分析,结构网格能够对关键区域进行局部划分,获得高质量的网格,进而得到更可靠的计算结果。如在处理近壁面附近层时,常常需要长的薄单元,由于四面体单元在某一特定方向拉长时,其数值求解特性会变差,计算误差大,方程不易收敛,而使用结构网格可以大大减小计算误差。因此学会建六面体网格将大大提高计算的收敛速度和准确度。除了这两种类型的网格外,ICEMCFD还可以建四面体网格、棱柱体网格、六面体混合网格及表面网格,适合不同边界条件的需要。
在使用ICEMCFD生成好网格后用Pre-MeshQuality命令对该网格进行质量检查并用Smooth命令对可以通过Move网格质量进行优化。若网格出现负体积,
CFD(ComputationalFluidDynamics)是计算流体力学的简称,其核心就是基于现代计算流体力学的有限单众所周知在计算流体力学中,其计算网格的元分析技术。
质量好坏,边界条件的施加以及求解器的求解方法等都将直接影响流体数值仿真的精确程度。目前市面上广泛使用的流体分析软件非常多,有的十分繁琐。本文提供了一种使用流体软件ANSYSCFX和其前处理软件-ANSYSICEMCFD相结合数值仿真方法,可以大大提高工程分析的效率。2
ANSYSICEMCFD和CFX简介
CFX是目前处于世界领先地位的CFD软件之一,广泛用于模拟各种流体流动、传热、燃烧和化学反应等问题。CFX的前处理模块是ICEMCFD,它提供广泛的CAD和CAE软件接口,可以快速生成多种形式的网格。它在生成网格时,可实现边界层网格自动加密、流场变化剧烈区域网格局部加密、分离流模拟等,这大大提高了数值仿真的精度。
下面以计算浸入式水口浇钢温度场和流场为例来说明在ANSYSICEMCFD和CFX中的数值仿真技术。33.1
运用ICEMCFD建立CAE模型几何模型的导入及清理
ICEMCFD提供了广泛的CAD软件接口,当几何实体模型建好以后可以通过ICEMCFD中的File/ImportGeometry直接导入CATIAV4、STEP、IGES、UG、Pro/E等格式的CAD几何模型文件。由于不同的文件格式导入
机械工程师2008年第12期
65
制造业信息化
MANUFACTURINGINFORMATIZATION仿真/建模/CAD/CAM/CAE/CAPP
Vertex命令拉动Block顶点,直到质量达到标准。对于湍流边界层,我们需要用Pre-MeshParams/MeshingParameters菜单来对边界
图1有限元模型
为水口上部定义速度为2m/s,入口温度恒为T=1580℃,湍流强度设为中等。出口处则设置为静压力为0个大气压边界。4.4
设置求解参数
在求解之前,可通过SolverControl设置对流格式和时间步长以及收敛精度。对于瞬态求解还需设置OutputControl。最后通过WriteSolverFile写出*.def文件进行求解计算。4.5
计算结果的后处理CFX可以用多种方式显示和输出计算结果,例如,显示速度矢量图、压力等值线等温线图、压力云图、
图、流线图、绘制XY散点图、参差图,生成流场变化的动画,进行粒子追踪等。图2、3分别表示水口浇钢后的温度分布云图和钢水速度矢量图。5
结论
在ANSYSICEMCFD中建立有限元模型十分方便快捷,能够满足不同流体边界条件的需要,而且其模型能够很好地导入到ANSYSCFX中施
图3钢水速度矢量图图2水口的温度分布云图
网格和某些特殊部位进行加密,这可以大大提高计算
图1是采用全六面体生成的该容器有限元模型。该精度。
71568个单元。模型共有68007个节点,3.3
导出ANSYSCFX格式的有限元模型
ANSYSICEMCFD本身附带了多种有限元求解软件的输入输出转换器。用户可以通过Output/SelectSolver选择求解器并通过Output/WriteInput来输出文输出文件格式默件。本文选择的求解器是ANSYSCFX,认为*.CFX5。44.1
运用ANSYSCFX求解计算导入有限元模型
CFX可以导入不同格式的网格,如CFXMesh(gtm)、Ansys及Nastran等,从ICEMCFD中导入网格直接选择ICEMCFD网格格式。4.2
设置流体或固体的物理参数
设置正确合理的材料参数是进行仿真分析的前提,在CFX中所有流体或固体参数均可通过Materials来创建、存储和管理。CFX中自建了理想气体、水和铝铜钢等常用流体和固体的物理参数。分析人员也可通过CreateMaterials来设置新的流体或固体的密粘度、导热系数和比热容等。由于ICEMCFD本身度、
不存在量纲,因此在CFX中设置流体或固体的物理参数时应该注意单位的协调。浸入式水口和钢水物理性质及如表1。
表1
参数名称水口耐火材料1580℃钢水加边界条件进行求解计算,在CFX后处理中也能够很好的反应计算结果,这种基于ANSYSICEMCFD和CFX的流体数值仿真技术能够大大提高工程流体分析的速度和精度。
[参考文献]
—CFD软件原理与应用[M].北[1]王福军.计算流体力学分析——
京:清华大学出版社,2004.
[2]NagashimaT,DevelopmentofaCAEsystembasedonthe
node-by-nodemeshlessmethod[J].Comput.MethodsAppl.Mech.2000,187:1-34.Engrg.,
[3]杜平安.有限元网格划分的基本原则[J].机械设计与制造,2000
(1):36-38.
[4]刘希云,等.流体力学中的有限元与边界方法[M].上海:上海交
通大学出版社,1993.
[5]PatankarSV.NumericalHeatTransferandFluidFlow[M].
WashingtonDC:Hemisphere,1980.
[6]罗会信,张学军.不同预热方式下浸入式水口的温度场[J].炼
钢,2002(1):42-45.
(编辑立
明)
浸入式水口和钢水物理参数及常用单位
m-3粘度/N·密度/kg·s·m-2热传导系数/W·m-1·K-1比热容/J·kg-1·K-1220069400.00510.12560.512008794.3施加边界条件
所谓边界条件,是指求解域的边界上所求解的变量
或其一阶导数随地点及时间变化的规律。所有CFD问题都需要有边界条件,对于瞬态问题还需要有初始条件。流场的解法不同,对边界条件和初始条件的处理方式也不一样。在CFD模拟时,基本边界条件包括:流动
[1]
给定压力边界、壁面边界、对称边界和周期进出口边界、
性边界等。不同的求解模型,边界条件的选择也不同。本文中浸入式水口浇钢温度场和流场采用稳态求解,边界条件如下:(1)固体表面
浸入式水口与外界接触的面
设为固体表面。与外界接触的表面则设置成对流换热,对流换热系数为12W/(m·K),设外界温度恒为300K。该压力使用无滑移边界条件。(2)进出口面对速度、
入口
的速度垂直于所在面,并视为稳定的一维流动。入口处
!!!!!!!!!!作者简介:许蕾(1983-),男,硕士研究生,研究方向为机械CAD/CAM
及机械系统虚拟现实。
收稿日期:2008-08-20
66
机械工程师2008年第12期