DOI:10.19476/j.ysxb.1004.0609.2001.s1.009
铝合金铸件缩孔缩松预测新技术
徐宏 程军 钟雪友
华北工学院材料工程系
北京科技大学材料科学与工程学院 太原030051
北京100083
摘 要:
在铝合金铸件三维温度场计算的基础上 , 提出了采用动态多熔池等效液面收缩量法对缩孔进行预测 ;并综合质量守恒方程、动量方程和氢含量守恒方程 , 建立了数学模型对缩松进行预测。通过实验验证了数学模型的正确性
关键词:
温度场 ;缩孔 ;缩松 ;预测 ;
中图分类号: TG290
收稿日期: 2000-08-18
New technology of cavity and porosity predication for aluminum alloy casting
Abstract:
On the basis of 3 D temperature fields calculation of aluminum alloy casting, the liquid face shrinkage capacity method of varied several melt pool is put forward to predicate shrinkage cavity; also the mathematical model of porosity predication which consists of continuity equation, motion equation and conservation equation of gas content is established. The above models are proved correct by comparing the simulation and test results.
Keyword:
temperature field; shrinkage cavity; porosity; predication;
Received: 2000-08-18
一般认为缩孔是由凝固收缩引起液面下降而形成的, 采用等效液面收缩量法
[1 ]
定量预测缩孔形成时, 假设所有的补缩都来自铸件最高处的可流动单元, 显然这种处理方法没有考虑是否存在补缩通道。预测缩松时, Compell
[2 ]
指出, 在铸件的整个凝固过程中, 补缩由液体补缩、质量补缩、枝晶间补缩、爆炸充填和固体补缩5种方式组成, 形成缩松的主要原因是质量补缩和枝晶补缩产生困难。目前缩松预测判据很多, 但一般缩松预测只是基于凝固的热状况, 得出凝固时间、温度梯度、冷却速度等热参数来获得预测判据
[3 ,4 ]
, 这些判据对于窄结晶区合金缩松预测效果较好, 对于宽结晶区合金却不太适合
[5 ]
。缩松的形成应综合考虑铸型内热状况、枝晶间流动状况、金属液中过饱和气体析出、金属静压力的影响
[6 ]
。
本文作者在缩孔预测上提出了动态多熔池等效液面收缩量法, 缩松的预测综合运用了质量守恒方程、动量方程和气体含量守恒方程。
1 温度场模型
式中T为温度, c为比热, ρ为液相和固相的平均密度, L为结晶潜热, λ为导热系数, u, v, w为枝晶间液体流速, φL 为液相体积分数, 与温度有关。
等号左边第2项表示枝晶间液体流动的热对流, 由于结晶潜热释放, 两相区温度梯度和液相体积分数比较小, 可以忽略不计, 式 (1) 变为
其中φL 可采用式 (3)
[7 ]
计算:
式中Tm 为纯金属熔点, TL 为液相温度, k为溶质平衡分配系数。
2 动态多熔池等效液面收缩量法
在铸件的凝固过程中, 往往同时存在多个孤立的熔池, 凝固收缩时各个封闭熔池的最高液面金属液均向下补缩, 由此预测各孤立熔池的缩孔或单元致密度。
2.1 动态多熔池的判定
1) 在tn 时刻对所有铸件单元赋予熔池编号“-1”。
2) 判断全部单元, 定φSi >φc 的单元凝固状态号为“1”, 定φSi <φc 的单元凝固状态号为“0” (φSi 为i单元的固相体积分数, φc 为合金的临界固相体积分数, φSi >φc 的意义为此单元没有补缩能力) 。
3) 任选一个凝固状态编号为“0”且熔池编号为“-1”的单元进行孤立熔池判断。此单元与该单元相邻的“0”单元以及与这些“0”单元中任意一个相邻的单元必处于同一个孤立熔池给这些单元赋予熔池编号“1”。
4) 重复步骤3) , 依次赋予熔池编号“2, 3, …”直至tn 时刻全部单元的熔池编号不为“-1”。
5) tn+1 时刻重复步骤1) ~4) 。
用该方法能确定凝固过程中各个时刻的熔池个数及熔池状况。
2.2 孤立熔池的等效液面收缩量法
金属液刚充满型腔时所有单元处于一个熔池中, 随着凝固的进行, 铸件中有的部位补缩通道不畅而导致形成多个孤立熔池。将等效液面收缩量法用于各个孤立熔池, 并假设金属液补缩仅考虑液相向固相转变时的体积收缩, 不考虑热胀冷缩影响;当单元的固相体积分数达到临界值时, 此单元丧失补缩能力。
1) 单元从液相转变为固相时体积收缩率为
式中ρS 为固相密度, ρL 为液相密度。
2) i单元从 (t-Δt) 时刻到t时刻的体积收缩量为
3) 如果从 (t-Δt) 时刻到t时刻, P号熔池中有n个单元发生液固转变, P号熔池总的体积收缩量为
4) P号熔池最高液面 (假设为第n层) 单元的液相体积之和为Vn , 此时可能出现2种情况:
(1) 若Vn >ΔVP , 液面仍处在第n层, 对于φSi <φc 的各个单元, 其单元液量减少, 单元液相体积为
式中Nn 为第n层φsi <φc 的单元个数, Vi, n 为第n层i单元液相体积。
(2) 若Vn <ΔVP , 则第n层φSi <φc 的单元变成空单元, 液面下降至第 (n+1) 层, 第 (n+1) 层中φSi <φc 的单元其金属液量进一步减少, (n+1) 层单元液量变为
当然若第n层和第 (n+1) 层中可用于补缩的合金液不足, 则液面将下降至第 (n+2) 层, 处理方法类似。
5) 其余熔池处理方法重复步骤3) ~4) 。
3 缩松预测方法
假设:1) 枝晶骨架不能流动;2) 缩松形成起源于一次枝晶的根部, 这是由于一次枝晶间距比二次枝晶间距大几个数量级, 二次枝晶间距内形成气孔的表面阻力很大
[8 ]
;3) 达西定律描述的是枝晶域内液体流动的阻力
[9 ]
。在此基础上建立单元体积内缩松预测模型。
1) 质量守恒方程
式中第1项代表收缩量, 第2, 3, 4项代表枝晶流动的液体进入量, 第5项为缩松量。
2) 动量方程
在两相区内描述枝晶间流动的动量方程可用达西定律表示:
式中k为渗透率, k=φL 3 d2 /[180 (1-φL ) 2 ];d为晶胞直径, d=αΔti n ;α和n为常数, 其值分别为7.5×10-6 和0.39
[10 ]
;Δti 为单元i的凝固时间。
缩松中的气体压力为
式中p为金属液体压力, σLG 为液气表面能, r为气孔半径。
3) 气体含量守恒方程
在铝合金中溶解氢时, 氢含量满足下列守恒方程:
式中[H]0 为初始氢含量, [H]S , [H]L 分别为在固相和液相的氢含量。
式中KS 和KL 分别为固相和液相的气体平衡常数。
联立求解方程可对缩松产生的位置及数量进行预测。
4 实验研究及验证
铝合金成分如表1所示。实验条件如表2所示, 含氢量通过测氢仪测得。对实验铸件剖切后分析缩孔缩松分布发现每个孤立熔池都产生缩孔, 并且缩松的含量随氢含量的增加而增加, 如图1所示。缩松的含量随冷却速度的加快而减少, 如图2所示。
表1 实验用铝合金成分 Table 1 Composition of aluminium alloy
表2 实验条件 Table 2 Experimental parameters
图1 H含量对缩孔及缩松分布的影响 (砂型)
Fig.1 Influence of hydrogen content on distribution of shrinkage capacity and porosity (sand casting)
(a) —[H]=3×10-4 mL/g; (b) —[H]=1.6×10-3 mL/g
在三维温度场模拟的基础上, 对试样进行缩孔及缩松的模拟预测, 结果如图3~5所示。
图2 不同冷却速度对缩松分布的影响
Fig.2 Influence of cooling rate on porosity
(a) —Metal mould; (b) —Sand mould
图3 缩孔模拟结果 (砂型)
Fig.3 Simulation results of shrinkage capacity
(a) —Sample 1; (b) —Sample 2
图4 试件一缩松模拟结果 (砂型)
Fig.4 Simulation results of porosity in sample 1 (sand casting)
(a) —[H]=1.6×10-3 mL/g; (b) —[H]=3×10-4 mL/g
图5 试件二缩松模拟结果
Fig.5 Simulation results of porosoity in sample 2
(a) —Sand mould; (b) —Metal mould
从模拟结果可见每个孤立熔池顶部都形成了缩孔, 并且含氢量和冷却速度对缩松有较大影响。通过对比可以看出, 计算机模拟结果与实验结果一致。
5 结论
1) 采用动态多熔池等效液面收缩量法可正确定量预测缩孔。
2) 建立了缩松预测的综合数学模型。该模型综合了铸件形状、合金液枝晶间流动和含气量的影响, 提高了缩松预测的精度。
3) 通过实验研究及验证, 证明了数学模型的正确性。
参考文献
[1] ImafukuI , ChijiiwaKA .Mathematicalmodelforshrinkagecavitypredicationinsteelcastings[J].AFSTransactions, 1983, 10:527-540.
[2] CompellJ.Ontheoriginofporosityinlongfreezing rangealloy[J].BritishFoundryman, 1964 (4) :147-157.
[3] NiyamaE , UchidaT , MorikawaM , etal.Predicatingshrinkageinlargesteelcastingfromtemperaturegradi entcalculation[J].AFSInternationalCastMetalsJournal, 1981, 2:16-22.
[4] LeeYW , ChangE , ChieuCF .Modelingoffeederbe haviorofsolidifyingAl7Si0.3Mgalloy platecasting[J].MetTrans, 1990, 21B :715-722.
[5] SuriVK , PaulAJ, EI kaddahN , etal.Determinationofcorrelationfactorsforpredicationofshrinkageincast ing—partI :predicationofmicroporosityincastings, AGeneralizedCriterion[J].AFSTransactions, 1994, 102:861-867.
[6] ChiouIJ, TsaiHL .Modelingofporosityformationincastings[J].AFSTransactions, 1990, 98:823-830.
[7] FlemingsMC .SolidificationProcessing[M].McGrawillBookCo, NewYork, USA , 1974.
[8] PoirierDR , YeumKandMaplesAL .Athemodynamicpredicationformicroporosityformationinaluminum richAlCualloys[J].MetTrans, 1987, 18A :1979-1987.
[9] GanesanS , PoirierDR .Conservationofmassandmo mentumfortheflowofunterdendriticliquidduringso lidification[J].MetTrans, 1987, 18A :721-723.
[10] KimioK , RobertDP .Mathematicalmodelingofporos ityformationinsolidication[J].MetallurgicalTransac tionsB , 1985, 16B :359-365