Fluent数值模拟软件介绍及流体力学的数值处理方法.docxVIP

Fluent数值模拟软件介绍及流体力学的数值处理方法

3.2模型数值处理策赌

采用Flu旬t进行有限元数伯分析求解时，首先要创迁几何模型并合理的划分网格，定义材料热物理性销参数，如何确定时间步长的人小，以及处理相变潜热及固液移动界面等．陡否合理地处理好问题，对计算秸度和计箕效水以及结果的收敛性和和定性，甚至到数伯计算的成败的影响都是巨人的。

3.2.1求解区域的如何选择物理模型

根据模拟中会出现的物理过1用以及考虑到多方的作用影响，本文选择的是Flue11t中流动税型：非定倌流、湍流以及不可压缩流模型，选用的热传计税型是VOF多相流、副化凝固以及离散相模型。

3.2.2网情的划分

网格生成技术在数伯模拟中是关键的一环，计算的收敛性和计算结果的精度及计算效率都会受到网格生成技木的或多或少的影响。针对实际情况对工件进行合理划分，工件的三维杖型以及对其有限元网格的划分的具体悄况都在图 3-1中得以显现．祁规中越小的网格划分，计算精度则越高，但相对的对计算机的计算能力要求的也会和应的越禹。根据呏究说明，在数伯计纬网格生成时会有如下五方面的婓求：

(1)贴体性：在保证边界上的网格不被破碎的前提下，物理域上的边界与计釭网

格线相重合．

（2)光滑性：物理域上 的网格线必须是一条儿滑迕续的曲线：

(3)合理的叶t密分布．在温度场流场参数变化率较人的区域（如激光熔注熔池区域）采用较密的网格，

(4)正父1生 在T件表面上尽力保证网格线的正父性．井稚保边界上的计侍楷度·

(5)单伯性：要求物理域与计笢域1的点要一一对1芍，不能有网格线相交和霾叠。

图3-1 网格划分

Fig.3-1 Mesh

本义犊型的网格划分如图3-1所示，设激光诺移动的方向为X方向，垂直于激光移动方向的为Y方向，厚度方向为Z方向。由于热惯性，在整个工件上会存在较人的温

度梯度，靠近激光扫过的且线区墩尤为明显，必度梯度非常高。而扣反离热源越远，温度梯度则越小。为f保证计算精度及准度，在激光熔注域内采用密网格，血远离熔注区，则采用疏网格。中间采用自由网格划分进行过渡。这样，就可以在不卅加单元

和节点数炳的条件下提商计鸽栝度了。

2.3 材料喻透参数

本文在校拟中用的挂体材料为 10mm厚的 A231B镁合金，股初对为.?.OOtrunx50mt11xJ0mm。其主要元素含挝见表 3-1。密度为 1780kg/m气熔点为651°c,主要物理参数见表3-2。

表3-J AZ3l8镁令金的化，,为成分（匪径百分比/%)

化学成分

Al

Zn

Mn

Fe

Cu

s,

Mg

亚萤百分比

3.06

0.84

0.35

0.0024

0.00lS

0.012

Bal

农3-2 数值计算筷型中所用材料的物性参数

物涅虽

名称

AZ31B栈合金

p

密度1(kg／面）

1780

兀

熔点／k

924

丑

涕如（（．

1373

）＇

，

表面张力系数1(N的）

1.0

热膨胀系数/K·

IO·

。

Stcf画 Bo!比man常数/[W（I叫．K)]

52X10-$

g

垂力加速度／（m尼）

9.8

c

表面辐射系数

0.4

“

换热系数/{WI佃. K)]

25

Q

熔化措丹J[kJ/1电

83000

笃

环埮混反/(K)

98

其他物性参敌比如动力粘度、导热系数、比热容都是随温度变化为而变化的相关的函数，它们和温度的的关系曲线如图 3-2所示，函数关系如下·

对于AZ31B镁合金

Jt=

r-0.000017T+0.0220946

l

-0T+0.0360064

-0.0000019T+0.0025心3

0.000000426

[kg/(m·s)]

9.:.9TS!l09

1109Ts l.?.09

1209Ts1380

T1380

0.1071429T+43.607

｛入＝ －o283T+310459

0729r-471817

c?」:。T++7:3916

l10 5T+60835

1420

[W/m·K)

[J/(kg·K))

2935T673

673$;T773

T773

293SJ373

3