消失模铸造充型过程数值模拟_图文

1、基于计算机模拟的消失模铸造充填特性研究张立强1,2,卢远志1,张楚惠3,刘金水1,2(1.湖南大学汽车车身先进设计与制造国家重点实验室,湖南长沙,4100822.湖南大学材料科学与工程学院,湖南长沙,4100823.湖南基特制造有限公司,湖南湘阴,410502摘要:通过计算机仿真软件ProCAST对某公司车身覆盖件修边冲孔模上模座的消失模铸造充型过程进行数值模拟,研究了消失模铸造的充填特性,并对比分析了消失模铸造和传统空腔铸造充型过程的差异。模拟结果表明:消失模铸造存在着与传统空腔铸造截然不同的充型过程,由于泡沫介质的存在,使金属液充型受阻,充型时间长,并且充型速度呈现周期性交替变化。关键词:

2、消失模铸造;数值模拟;充填特性;充型速度中图分类号:TG249.6 文献标识码:A 文章编号:Research on the Filling Characteristics of EPCbased on the Computer SimulationZHANG Li-qiang1,2, LU Yuan-zhi1,ZHANG Chu-hui3, LIU Jin-shui1,2(1. State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body (Hunan University, Changsha, 4100

3、82, Hunan, China;2.School of Materials Science Engineering, Hunan University, Changsha, 410082, Hunan, China;3.Hunan Jite Manufacturing Co., Ltd, Xiangyin, 410502, Hunan, ChinaAbstract: The filling process of EPC (Evaporated Pattern Casting are simulated in this paper by using the computer simulatio

4、n software to study the filling process of EPC. In order to study the filling characteristcs of EPC better, the comparation of the filling process of the traditional vacant casting have been also simulated at the same conditions. The simulation results show that the filling process of EPC is differe

5、nt entirely from the traditional vacant casting. The filling velocity becomes slow and varies periodically because of the existense of the foam in the EPC.Key words: EPC, numerical simulation, filling characteristics, filling velocity1 引言消失模铸造工艺由于其具有铸件尺寸精度高、机械加工余量少、生产适应性强等诸多优点, 被认为是一项很有发展前景的近净型加工技术,

6、并被铸造界誉为“21 世纪的铸造新技术”。目前,我国消失模铸造技术日趋成熟, 研究单位和应用厂家越来越多13。消失模铸造与其他铸造方法的显著差异在其充型过程的特殊性,金属液在充型过程中包含着复杂的物理化学现象,包括热量传输、充型流动、化学反应、冷却凝固等。这些现象耦合在一起,并相互影响,使得金属液的充型过程更加复杂,最终影响到铸件质量,同时,掌握消失模铸造型内液态金属的充填特性还关系到浇注系统的设计及充满或不充满等重要问题的确定45。在实际生产过程中,由于铸件结构上的差异和具体生产条件的限制,同时铸造过程又是一个在完全封闭的型腔内完成的生产过程,因此很难通过观察或物理试验来研究消失模铸造的充填

7、特性,而结合计算机模拟的消失模铸造过程研究是一种高效的分析方法,有助于揭示其充型过程的流动规律,发现充填中的不足,为改善铸造工艺设计,提高铸件质量提供一个可靠、快捷的分析手段69。基于此,本文作者结合某公司车身覆盖件修边冲孔模上模座铸件的实际生产,通过计算机仿真技术,对消失模铸造充型过程进行了数值模拟,并对比分析了消失模铸造和传统空腔铸造充型过程的差异,研究了消失模铸造的充填特性。第一作者简介:张立强,男,1978年出生,硕士研究生,湖南大学材料科学与工程学院,湖南长沙(410082,电话:013337315361,E-mail:zlq2k11032 充填过程的计算机模拟相片1是某公司车身覆盖

8、件修边冲孔模上模座的实际消失模铸件,材料为HT300,重1.09吨。铸件尺寸为长1870mm,宽750mm,高350mm,主要壁厚35mm。铸件要求内部组织致密,表面3-5mm以下无裂纹、气孔、缩松等铸造缺陷。本文根据实际生产的铸造工艺,对该铸件进行了数值模拟试验。计算机模拟包括三个基本步骤:首先利用CAD软件建立包括浇冒口系统的三维几何模型,再通过parasolid格式的中间文件将其导入计算机模拟软件;其次对三维几何模型进行必要的前处理建立有限元模型,包括几何模型简化、网格剖分、计算参数设置等处理;最后利用软件进行数值模拟计算和结果分析。 相片1 实际消失模铸造铸件Photo1 The pr

9、actical EPC casting2.1 铸件几何模型利用三维CAD软件在零件几何模型的基础上,根据车身覆盖件修边冲孔模上模座的结构特点、消失模铸造方法的特殊性、以及实际生产中的铸造工艺参数(见表1。添加了铸造工艺特征,包括铸型、浇铸位置、浇注系统等。为了保证后续网格剖分的质量和数值计算的精度,需要对上述三维几何模型做适当的简化处理。表1 铸造工艺参数Table 1 Technologic parameters of pouring直浇道截面积(cm2 横浇道截面积(cm2内浇口截面积(cm2暗冒口截面积/数量(cm2/个明冒口截面积/数量(cm2/个浇注方式28 36 21 10/6 1

10、2/8 底浇注2.2 有限元模型2.2.1 网格剖分根据数值计算原理和计算机性能要求,为了得到精确的模拟结果和减少计算量,有限元模型采用非均匀四面体网格对几何模型进行网格剖分,即浇道和铸件采用较小的网格尺寸进行剖分,而砂箱则采用较大的网格尺寸进行剖分。如图2所示计算所采用的有限元网格规模为:单元总数为1045847,节点总数为207186。 图2 铸件与铸型的网格剖分Fig.2 Mesh generation of casting and mould2.2.2 计算参数有限元模型中计算参数是根据生产实际进行设定的,主要包括材料参数和浇注工艺参数等。消失模铸造区别于其他铸造方法的是,它包括金属-

11、泡沫-铸型三类介质材料。铸件材料为HT300,泡沫材料选择EPS(聚苯乙烯,铸型材料为树脂砂,拥有很好的透气性。材料的热物性参数如表2所示。表中合金密度、潜热与热导率都是随温度变化的,比热数值指0100范围内。浇注工艺参数包括浇注温度1400和浇注压力1.1atm 。作为数值计算的边界条件主要为不同介质之间的热交换系数,包括泡沫和砂型间的热交换系数为100w/m 2.k ,金属和砂型间的热交换系数为500w/m 2.k 。表2 合金及泡沫材料的热物性参数密度/kg·m - 3 比热 /kJ·(kg·k - 1潜热 /kJ·kg - 1 热导率 /w

12、83;(m·k - 1 固相温度 /k 液相温度 /k 铸件(HT300 (6.37.4103 0.532601458 23.636.4 1406 1522 EPS (聚苯乙烯25 3.7 100 0.035 603 6232.3模拟结果数值模拟是在高性能微机工作站上利用专用的消失模铸造分析软件完成的,计算时间平均为58小时,消失模铸造的充填过程计算结果见图3。图3a图3e 是消失模铸造不同时刻的充型形态,从图3a 为15s 时的充型状态,可以看出金属液从最先流经的内浇道中进入,之后在内浇口与横浇道基本保持同步填充。从图3b 与c 为32s 时的充型状态,当金属液流经内浇口1与内浇口

13、2进入铸型时,充型速度明显增加,而横浇道内的金属液充型速度相对缓慢,并且进入铸型后的金属液在填充底层的同时液面逐渐上升,图3d 是也说明了这种充型形态。图3e 表示62s 时金属液完全充满型腔。图4为消失模铸造直浇道处金属液在不同时刻的充型速度矢量图,消失模铸造充型过程中,充型速度是不断变化的,图a,c,e 中的充型速度相对较大,而图b,d,f 中的充型速度则相对较小,呈现出周期性交替现象。 (a 15s (b 32s内浇口1内浇口2 (c 32s (截面 (d 49s(e 62s图3 浇注压力为1.1atm 的消失模铸造,浇注开始后不同时刻的充型形态Fig.3 The filling sta

14、te of different time of EPC when 1.1atm pouring pressure (a 4s (b 11s (c 17s(d 23s (e 25s (f 26s图4 浇注压力为1.1atm 的消失模铸造,不同时刻的充型速度变化情况Fig.4 The variation of filling velocity of different time of EPC when 1.1atm pouring pressure3 消失模铸造充填特性研究实际铸件的模拟结果表明,在同样的浇注工艺条件下,消失模铸造充型过程因为受到金属液-泡沫-砂型3个介质的相互影响与作用,与传统空

15、腔铸造金属液-砂型双介质充型顺序与充型速度差别很大。下面对其充填特性进行讨论分析。3.1 充型过程的理论基础与传统空腔铸造不同的是,消失模铸造的充填过程是一个包含泡沫气化过程的流场和温度场相互耦直浇道底部A 处节点号592合的数学物理过程,因此其本构关系包括分别包括流场、温度场以及气化模型三个部分,下面分别加以描述。3.1.1 流场本构关系在铸件充型过程数值模拟中,将液态金属看作不可压缩流体,其流动过程服从质量守恒和动量守恒,数学形式为连续性方程和Navier-Stokes 方程(简称N-S 方程。连续性方程0=+=z w y v x u D (1N-S 方程 w F z p z w w y

16、w v x w u tw v F y p z v w y v v x v u t v u F x p z u w y u v x u u tu z y x 222+-= +-= +-= + (2 式(1中,D 为散度;u ,v ,w 为速度矢量在x ,y ,z 方向上的分量;式(2中,为金属液的密度,t 为充型时间,p 为流体内部压力;为运动粘度;F x ,F y ,F z 为单位质量力在三个坐标轴x 、y 、z 的分量;2222222/z y x +=为拉普拉斯算子。3.1.2 温度场本构关系流体在运动过程中,温度场的计算遵循能量守恒定律,在直角坐标系下其本构方程为S T k z T w y

17、 T v xT u c t T c +=+2 (3 式(3中,T 为温度,t 为时间,S 为潜热项,k 为热传导率,2222222/z y x +=为拉普拉斯算子。3.2 气化模型3.2.1 气化过程的物理模型消失模铸造是金属液-泡沫-砂型3介质相互影响与作用的结果,由于存在着泡沫气化阻力的影响作用,其气化模型可用图5表示。根据这个模型可以看出,高温液态金属使泡沫迅速熔解气化从而在液态金属与泡沫模样之间形成气隙,气隙里存在着三个不断变化的相(泡沫固体颗粒、液体泡沫与泡沫气体,并且由于金属液与泡沫模样的热交换,气隙里同时发生着大量的传热与化学反应。 图5 泡沫气化过程示意图Fig.5 Schem

18、atic of the foam degradation process 113.2.2 气化过程的数学模型 图5中的气隙内部包括两部分气体，一部分是充型过程中退化泡沫气化的气体，一部分是气隙内部 分气体逃逸后剩余的气体，它们分别满足理想气体状态方程 PV = R GT ，故有 _ P1V = R G 1T _ （4） （5） 3 P2 V = R G 2 T _ _ V 实验测定为(44.5Vmetal Vmetal （ 式(4中, P1 为退化泡沫气化后形成的气体压力(Pa； 为气隙的体积(m ， 为每个时间步长所流入的金属液体积）10； R 为气体常数(m/k，为29.18m/k； G1

19、 为泡沫气化后气体 重量(N， 根据质量守恒定律， 它等于退化的泡沫重量；T 为气隙实际温度(k， 近似等于液态金属温度。 式(5中， P2 为气隙内部分气体逃逸后剩余的气体压力； G2 为剩余的气体重量。 为了求出式(5中的 G2 ，根据气体层流时的公式(6，需要计算逃逸的质量流量 q 。 q= p gap p 0 p 真 R r气 （6） 式(6中， p gap 为气隙处的绝对压力(Pa； p0 为大气压力(Pa； p真 为砂型中的真空度(Pa；R为气隙四 周涂层的气阻(N·s/m ； r气 气隙处气体的重度(N/m ， r气 Pgap / R Tgap 或 r气 P真 / R

20、T外 ，当压力较低的真空 状态时，用前面的式子计算 r气 ，当大气压或差压铸造时用后面的式子计算 r气 。 假设瞬时气体温度与体积恒定，根据道尔顿气体分压定律，有 5 3 Pgap = P1 + P2 （7） 3.2.3 气化区充型压力的计算 消失模铸造充型过程受到气隙压力 Pgap 的影响，所以气化区实际充型压力 P实 需要修正为 P实 = P充 Pgap 型速度，因此随着 Pgap 的变化，充型速度将不断的发生改变。 (8 式(8中，是 P充 金属液的充型压力； P实 是气化区实际充型压力。实际充型压力得到修正后，得到新的充 消失模铸造充填特性分析 3.3 消失模铸造充填特性分析 由消失模

21、铸造充型过程流场与温度场耦合的本构关系和气化模型可知， 消失模铸造充型过程中由于 受到气隙压力的影响，使得充型速度减慢，同时，充型速度不断发生变化。为了更好的研究消失模铸造 充型特性，对传统空腔铸造的充型过程进行了数值模拟，其模拟结果如图 6。从图 6a 与 b，传统空腔铸 造过程中， 金属液先充满横浇道后， 再经多个内浇口同时充填型腔， 6c,金属液从内浇道进入铸型后， 图 先填满底层，然后液面再逐渐上升。并且在同样的浇注工艺条件下，传统空腔铸造的充型速度更快，如 图 6e，在 10s 时，金属液基本充满铸型，而图 3e 消失模铸造需要 62s 的时间才可以充满型腔。 (a 0.4s (b

22、0.5s 6 (c1s (d 4s 直浇道底部 A 处 节点号 2438 (e 10s 图 6 浇注压力为 1.1atm 的空腔铸造，浇注开始后不同时刻的充型形态 Fig.6 The filling state of different time of vacant casting when 1.1atm pouring pressure 1.1 1.0 0.9 0.8 8 节点号 592 7 6 节点号 2438 Velocity(m/s 0.6 0.5 0.4 0.3 0.2 0.1 0.0 0 5 10 15 20 25 30 Time(s 35 40 45 50 55 60 65 70

23、 75 Velocity(m/s 0.7 5 4 3 2 1 0 0 1 2 3 4 Time(s 5 6 7 8 9 10 11 (a 消失模铸造充型速度与时间的关系 (b 传统空腔铸造充型速度与时间的关系 (a The relation between the filling velocity and (b The relation between the filling velocity and time for EPC time for vacant casting 图 7 直浇道底部充型速度与时间的关系 Fig.7 The relation between the filling v

24、elocity and time 图 7 为消失模铸造和传统空腔铸造直浇道底部充型速度的变化曲线， 从图 7a 可以看出,此零件消失 模铸造充型过程中可分为三个不同阶段，并且充型速度是呈现周期性交替变化。0 到 26s 为第一阶段， 充型的平均速度 0.12m/s；26s 到 42s 为第二阶段，充型的平均速度 0.46m/s；42s 到 62s 为第三阶段， 充型的平均速度 0.74m/s。结合图 3 中的可视化结果可知，第一阶段是浇注系统的充填阶段；第二阶段 属于金属液沿水平方向充填铸型阶段； 第三阶段是金属液沿垂直方向充填铸型阶段。 由传统空腔铸造的 模拟结果图 7b 可知，其充型速度相

25、对较大（0 到 10s 的平均速度为 5.17m/s）使得充型时间短，且充型 速度是基本稳定的。 充型形态差异的具体原因是金属液充型过程中， 泡沫会迅速的熔解、 气化并在金属液与泡沫之间形 成气隙，气隙会使金属液和泡沫的界面间压力上升，从而改变了液态金属在充型过程中的流动特性。同 7 时泡沫气化液化退让的空间让液态金属占领的置换过程是有代价的， 那就是液态金属必须付出足够的热 量才得以进行，泡沫模型的退让还要有个过程，就是泡沫吸收能量需要一定时间，导致消失模铸造的充 型速度较传统空腔铸造慢。总之，消失模铸造的气化过程造成了充型过程出现巨大阻力，也就是泡沫模 型消失过程的阻力制约着整个充型过程。

26、此外，在充型前期，金属液开始充填浇注系统，未进入铸型， 界面前沿排气面积小，泡沫气化的气体向外逃逸阻力大，则金属液前进的阻力大，充型的平均速度小； 而充型后期，金属液大面积填充铸型，界面前沿排气面积大，泡沫气化的气体逃逸变得更加容易，金属 液前进的阻力小，导致充型的平均速度增加。 此零件消失模铸造充型速度存在周期性交替变化的现象， 其可能原因是消失模铸造在充型的开始阶 段，在泡沫未气化之前，金属液未受到气隙阻力影响，流动速度较大，当泡沫首先气化时形成较大的气 隙，此时金属液的流动速度超过泡沫模样的热解速度，气隙间距减小，气体向外界逃逸的面积变小，透 气性变差，使得气隙阻力增加,金属液的充型速度

27、减小；而另一方面，气隙的减小,增加了金属液前沿向 泡沫模样换热， 泡沫模样的分解速度开始大于已经减小的金属液充型速度,所以气隙间距开始变大,， 排 气面积增大， 流动阻力减小,充型速度又开始变大。 随着充型的进行,金属液流动速度又开始大于泡沫的 热解速度,气隙间距又开始变小,同时金属液前进的阻力大增,充型速度又开始变慢，如此往复循环，直 到气隙厚度恰好能使消失模铸造的退让与液态金属的前进保持平衡为止。 因此， 充型过程中气隙间距的 不断变化，成了稳定消失模铸造的充型过程的重要因素。 由于存在上述消失模铸造充填特性， 因此在消失模铸造浇冒口以及浇注工艺设计时均需要采取与传 统空腔铸造不同的措施。

28、在设计消失模铸造浇冒口设计方面，为了设计出合理的浇注系统，通常采用底 浇方案。在设计浇注工艺方面，通常在浇注金属液的同时，点燃四周从铸型中逃逸出来的泡沫气体，这 样可以增加铸型负压度，加快气体逃逸，减少气隙压力变化辐度，从而使充型过程更加平稳，保证充填 的连续性与完整性，获得优质的铸件。 4 结论 通过上述分析与讨论， 说明基于计算机模拟的消失模铸造充填特性研究是一种有效的研究方法， 同 时利用计算机可视化技术可以较好的观察分析并揭示消失模铸造的充型规律，其主要充型特性是： (1 与传统空腔铸造金属液-砂型双介质不同，消失模铸造充型过程受着金属液-泡沫-砂型 3 个介 质的相互影响与作用,金属

29、液与泡沫间存在气隙， 形成充型阻力,使得消失模铸造充型过程和传统空腔铸 造差别很大,并且其充型速度较传统空腔铸造要慢的多。 (2 消失模铸造的充型过程中，因气隙间距和气体逃逸能力的不断变化，可能导致了充型速度呈现 周期性交替变化。 其它结构的零件消失模铸造充型过程中的充型速度是否也存在周期性交替变化， 还需 要进一步的实验验证。 (3 气化间隙是消失模铸造数值模拟的关键参数，同时也是影响其充填特性的关键参数。 参考文献： 参考文献： 1 俱英翠,李日,刘宏伟等. 模拟对比分析消失模铸造与普通铸造充型过程J.材料工艺,2006,6(2: 33-36. 2 陈 尧 剑 , 黄 天 佑 , 康 进 武 . 消 失 模 发