铸件充型凝固过程数值模拟

HEATS

1 概述
    欲获得健全的铸件，必先确定一套合理的工艺参数。数值模拟或称数值试验的目的，就是要通过对铸件充型凝固过程的数值计算，分析工艺参数对工艺实施结果的影响，便于技术人员对所设计的铸造工艺进行验证和优化，以及寻求工艺问题的尽快解决办法。
    铸件充型凝固过程数值计算以铸件和铸型为计算域，包括熔融金属流动和传热数值计算，主要用于液态金属充填铸型过程；铸件铸型传热过程数值计算，主要用于铸件凝固过程；应力应变数值计算，用于铸件凝固和冷却过程；晶体形核和生长数值计算，主要用于金属铸件显微组织形成过程和铸件机械性能预测；传热传质传动量数值计算，主要用于大型铸件或凝固时间较长的铸件的凝固过程。数值计算可预测的缺陷主要是铸件形成过程中易发生的冷隔、卷气、缩孔、缩松、裂纹、偏析、晶粒粗大等等，另外可以通过数值计算，提出合理的铸造工艺参数，包括浇注温度、铸型温度、铸件凝固时间、打箱时间、冷却条件等等。目前，用于液态金属充填铸型过程的熔融金属流动和传热数值计算以及用于铸件凝固过程的铸件铸型传热过程数值计算已经比较成熟，逐渐为铸造厂家在实际生产中采用，下面主要介绍这两种数值试验方法。
    1.1 数学模型
    熔融金属充型与凝固过程为高温流体于复杂几何型腔内作有阻碍和带有自由表面的流动及向铸型和空气的传热过程。该物理过程遵循质量守恒、动量守恒和能量守恒定律，假设液态金属为常密度不可压缩的粘性流体，并忽略湍流作用，则可以采用连续、动量、体积函数和能量方程组描述这一过程。
# D( @2 f% r" R z+ h, W, P( P" }$ M7 Q* z8 j4 v6 ~# @. g ], ~8 ]1 e" Z4 Q! N3 f) E7 Q5 @, S a: M9 o: e8 ]4 j' r: j; y: C! z$ Q# _1 w' W% H. }1 z2 N% A7 U' b# l z/ a( y7 r9 t. W: V" K1 j% p9 B% v: |0 }5 c% x8 r, P1 t) G |) l+ G: m& ~. f$ [2 c/ L* T( W( V5 J# X3 ^+ s5 z! w% ^

质量守恒方程

¶ u /¶ x+¶ v /¶ y+¶ w /¶ z= 0   (2-1)

动量守恒方程

¶（r u）/¶ t+u¶（r u）/¶ x+v¶（r u）/¶ y+w¶（r u）/¶ z              = -¶ p/¶ x+m (¶ 2u /¶ x2+¶ 2v /¶ y2+¶ 2w /¶ z2)+r gx   (2-2a)

6 S! P5 C, t4 r ¶（r v）/¶ t+u¶（r v）/¶ x+v¶（r v）/¶ y+w¶（r v）/¶ z 6 D6 V/ O- J- X" a" F       = -¶ p/¶ y+m (¶ 2u /¶ x2+¶ 2v /¶ y2+¶ 2w /¶ z2)+r gy              (2-2b)

¶（r w）/¶ t+u¶（r w）/¶ x+v¶（r w）/¶ y+w¶（r w）/¶ z 1 R% w6 { e# X1 X* M g           = -¶ p/¶ z+m (¶ 2u /¶ x2+¶ 2v /¶ y2+¶ 2w /¶ z2)+r gz             (2-2c)

体积函数方程

. l$ C( u# A& I3 m. T ¶F/¶ t+¶（Fu）/¶ x+¶（Fv）/¶ y +¶（Fw）/¶ z= 0             (2-3)

能量守恒方程

¶（rcpT）/¶ t+¶（r cpuT）/¶ x+¶（r cpvT）/¶ y +¶（r cpwT）/¶ z        = ¶（lT/¶ x）/¶ x+¶（lT/¶ y）/¶ y +¶（lT/¶ z）/¶ z +qv              (2-4)

/ u$ N1 B2 q+ L% n$ R     & c9 p% \, b- i [7 z4 N ||| ( A8 O( s" @) k! Z) c2 H, x4 `" k, c     式中  u,v,w —— x,y,z 方向速度分量(m/s)；r —— 金属液密度(kg/m3)；t —— 时间(s)；p —— 金属液体内压力(Pa)；μ —— 金属液分子动力粘度（Pa.s）；gx, gy, gz —— x,y,z 方向重力加速度(m/s2)；F —— 体积函数，0£F£1；cp —— 金属液比定压热容[J/(kg.K)]；T —— 金属液温度(K)；λ —— 金属液热导率[W/(m.K)]；q  —— 热源项[J/(m3.s)]。

    1.2 实体造型和网格剖分
    欲进行三维充型凝固过程数值模拟，首先需要铸件的几何信息，具体地说是要根据二维铸件图形成三维铸件实体，然后再对铸件实体进行三维网格划分以得到计算所需的网格单元几何信息。 
    利用市场上成熟的造型软件（如UG，ProE, Solid-Edge, AutoCAD等）进行铸件铸型实体造型，然后读取实体造型后产生的几何信息文件（如STL文件），编制程序对实体造型铸件进行自动划分，这种方法可以大大缩短几何条件准备时间。剖分后的网格信息包括单元尺寸和单元材质标识。

    1.3 数值计算方法
    用于铸件充型凝固过程数值计算的方法主要有3种：有限差分法、控制容积法（又称有限体积法）和有限元法，后两种方法采用的较少，目前在铸造市场上推广的一些数值模拟软件大部分采用的是有限差分法。以有限差分法为例，方程（2—1）的离散采用中心差分方法，方程（2—2）和（2—4）的离散采用上风方案和中心差分方案相结合的方法。
    充型过程中液态金属自由表面是不断变化的，每个时间步长对应的计算域均不相同，新的计算域的确定是通过求解方程(2-3)得到的。普通的数值方法在离散方程(2—3)时将造成很大的假扩散问题，计算结果将出现界面模糊（Smearing）现象，在 F=1 与 F=0 之间存在大量自由表面单元。为了得到清晰的自由表面，美国的科研人员发展了一种VOF（Volume of Fluid）方法，较好地处理了流体流动过程的自由表面计算问题。目前在计算流体力学领域已经在VOF方法的基础上开发了一些更准确的方法，可以获得更精确的流体流动过程自由表面变化。
    充型凝固过程数值计算步骤如下: 
    1) 将铸件和铸型作为计算域，进行实体造型、剖分和单元标识。
    2) 给出初始条件、边界条件和金属、铸型的物性参数。
    3) 求解体积函数方程得到新时刻流体流动计算域。
    4) 求解连续性方程和动量方程，得到新时刻计算域内流体速度场和压力场。
    5) 求解能量方程，得到铸件和铸型的温度场及液态金属固相分数场。
    6) 增加一个时间步长，重复3)～6)步至充型完毕。
    7) 计算域内流体流动速度置零，调整时间步长。
    8) 将充型完毕时计算得到的铸件和铸型温度场作为初始温度条件，求解能量方程至铸件凝固完毕。
    9) 计算结果后处理，进行铸造工艺分析、铸件缺陷预报和工艺参数优化工作。 
    1.4 应用实例
    下面给出一个利用ZCAST软件对低压铸造铝合金轮毂铸件的铸造工艺参数进行优化的实例。 
    在轮毂铸件低压铸造过程中，自动控制模具温度对防止产生缩孔缩松缺陷、组织粗大以及生产周期延长很重要。在金属型模具设计中模拟计算和水冷控制器可以作为一种重要的工具。通过模拟循环过程中轮毂铸件的流动和凝固设计合理的冷却系统。通过使用水冷控制器控制模具的热量散失或积聚。在铸造厂家的低压铸造轮毂铸件生产线上通过减少循环时间和铸造缺陷可以提高生产率。
   

|||

" G4 Q( A+ N8 G) ^9 Z

    模拟计算步骤如图2-1所示。首先进行前处理，前处理的主要工作是计算域内铸件铸型的3维实体造型，然后在ZCAST软件中导入实体造型文件，输入计算边界条件,包括入流边界、初始温度、热电偶位置设置等。

9 q, B: F4 m7 N; B% K

图2-1 数值计算过程

    前处理完毕之后开始数值模拟计算，包括一个循环当中的合型和开型过程，合型过程考虑充型和凝固，开型过程只考虑铸型型腔的传热。数值模拟计算部分可以考虑多个循环过程，以便于观察生产过程是否稳定。

4 [1 T. j0 m+ j5 b8 X5 E) D4 L6 c# \

    模拟计算完成后是后处理过程，主要是对模拟计算结果进行可视化分析，对铸造工艺参数进行调整。

8 M# v E( V1 ^ M

- S' H* G% I0 R# }" A

图2-2 试验设备，TC为热电偶位置

/ m6 B, ?+ e! g, E5 k2 _( F9 f: ], ^

    图2-2为试验设备，包括冷却水供给系统、多通道温度采集系统（图左）、低压铸造设备和铸型（图右）。

; _& h" m$ h* Y

（a）                                  （b）

) b3 L- f6 o- T: G) b3 Y

    图2-3 实体造型和剖分，实体造型中包括铸件、浇注系统、铸型和冷却通道。TC为热电偶位置。

4 g% P S! p6 X8 n" q7 [8 L

    图2-3（a）为计算域内实体造型结果，图2-3（b）为计算域内网格剖分结果。表2-1列出了模拟计算条件。界面换热系数（cal/cm²×sec×°C）：铸件/铸型为0.03，铸型/空气为0.004，液态金属动力粘度为0.1Pa×s。

7 R, L) X% Z- c2 }

表2-1 模拟计算中用到的物理参数

+ D" x: z- V5 B' H

$ I/ G% G. d1 z! z5 v2 U) ^& }( `5 X8 b! e% F) r6 M5 z o7 i; q! r% d% L* _& D3 m$ c; g) P, w. b) B$ v& {3 S/ F- h4 c, _- m* O6 m( `! P4 ?9 W$ o, k' z/ M/ U {7 F$ ~% H1 L) C* H. S, D( x# W& _# l4 m( y! e3 g s& O( \: J/ J9 q4 E3 k7 r# P. |( i. l- ?6 V6 G0 {! \ j, |; w; A; ?

材料	密度 g/cm3	比热容 cal/g×°C	热导率 W/m×°C	潜热 cal/g	$ q0 |3 g/ D! P3 K 液相线 °C	; J8 t) k. F0 y& M4 l+ O, m" z+ C/ _* a 固相线 °C	初始温度 °C	环境温度 °C
金属	% K% r+ @$ K+ B: [3 x6 n) x 2.69	% @: M9 j; C* M r" T5 d 0.23	0.34	93	' W' Y6 u, I' T8 m6 h9 R 620	577	: A. V5 [( q! [* U 740	; }- b* X/ h1 j3 U -
$ L9 ?7 X- B0 F5 n6 k! q 铸型	4 ^; W4 \6 E4 c" X * L# K5 c9 C( u$ d6 L3 y