说说铸造气孔物理模型研究进展.doc

1、说说铸造气孔物理模型研究进展时间:2009-10-29 09:57 作者:铸造模型 点击: 57次成型铸造是生产高质量铝合金铸件的经济实惠的方法,在铸造过程所形成的缺陷中,气孔缺陷是长期以来困扰业界的一个难题.1 背景介绍成型铸造是生产高质量铝合金铸件的经济实惠的方法，在铸造过程所形成的缺陷中，气孔缺陷是长期以来困扰业界的一个难题。与气泡形成有关的主要现象有 ：液相主体的流体流动与质量传递；热传导与凝固；固相晶粒的形核与长大，以及它与热和溶质浓度的相互关系；枝晶液体的流动约束及以后的局部压力改变；偏析与溶解的氢的求解；气泡的形核与长大；气泡对晶粒长大的影响。由于各种因素对气泡形成与长大的促进作

2、用，使其成为不可避免的缺陷。而且，气孔被认为是对抗疲劳与抗断裂性能最有害的缺陷之一 。人们就此进行了广泛的研究，Gabb等人研究了铸造合金的疲劳寿命，推断出裂纹多数发生在靠近表面的微孔处 。由于缺陷伴随整个构件，铸件的生产过程控制要求很高，然而，昂贵的质量控制仅使缺陷的大小控制在允许范围之内。因此，研究铸造气孔的形成原因是非常重要的 。铝合金中气孔形成有两大主要原因：（1）凝固体积收缩，枝晶间补缩失败；（2）氢从合金中析出，氢气泡形成；或者是两个原因的共同作用 。体积收缩是由于固液相合金的密度差异。当凝固继续进行时，体积减少，而且周围的金属液流动补缩。依据固相的数量与分布，补缩流动被阻碍或者完

3、全阻塞。当补缩不充分，就形成空隙。氢在固/液相中不同的溶解度导致气泡的生成。氢在固相铝合金中的溶解度远远低于其在液相中的溶解度。凝固时，多余的氢原子将从新形成的固相中逃逸到周围的液相中 。当氢含量达到了金属液中的临界溶解度，氢分子气泡形成，依据局部氢浓度与扩散率，它或长大或减小。最终气孔的数量、大小、和分布将决定于许多因素，如晶粒提纯剂、变质剂、化学成分与冷却率等 。人们为了提高这些材料的力学性能，通过控制工艺参数和冶金的化学过程，实现晶粒微结构的精炼 。如铝硅合金在熔体中加入少量晶粒提纯剂来控制铸件的晶粒微结构 ，变质剂把铝合金中共晶硅的形状从针状变为纤维状的亚共晶合金，改变了它的力学性能

4、。不同的方法精炼不同的微结构，但是，它们对多数的气缩孔缺陷有相同的影响 。文献 考虑了冷却率、氢浓度、晶粒提纯剂，共晶变质剂等因素，建立了这些变量对气孔形成影响的二维模型。所建立的二维模型将形核的数学模型与晶粒结构的预测模型进行了耦合，在给定的位置预测微孔的大小、形态和分布。人们在熔炼过程和铸造过程进行了大量的研究来减少气孔缺陷，发现即使是真空冶炼，仍然会产生有害的气孔。为了满足铸造工厂的需求，冶金研究者作了大量的努力来探索铸造铝合金潜在的物理原理，试图通过实验分析与数值模拟相结合的方法来对气孔的形成进行预测。这些研究成果对实际生产起到了很好的指导作用，并且成效显著。然而，影响气泡的形核与长大

5、的因素，人们还不能够完全理解，而且由于缺乏一致的、可靠的物理属性与参数，使建模过程更为复杂。理想的气孔模型应能够预测所有气泡的位置、大小和形状，因为即使单个气孔也会对铸造合金的力学性能产生决定性的不利影响。目前为止，尚不存在这样的模型，如果这种模型能够实现，将非常复杂。2 模型的现状2.1 模型的发展一些早期的商业软件致力于凝固过程的模拟，可以预测宏观的缺陷（如自由表面上的缩孔以及封闭区出现的气孔）。但是这些软件都是基于纯粹的热物理数据，没有耦合微孔的预测模型。而且封闭区域与未封闭区域之间的转变并没有考虑。人们做了大量的工作来对气泡的形成进行建模。最多的是铝合金中气孔模型，其次是镍合金，钢。但

6、是先前的模型在应用上存在一定的局限性 ，它们基于经验的判据函数，和达西定律与质量守恒的耦合求解模型 ，大多数模型只能预测气孔的百分比。许多研究者试图建立一个综合的模型，把所预测的微孔的大小和位置作为过程变量的函数。已有的研究表明，缺陷面积与缺陷到自由表面的距离是决定缺陷对铸件的疲劳寿命影响程度的重要因素 ，位于铸件中心的大的气泡，可能对力学性能和疲劳强度影响并不大，但是接近于表面的非常小的气泡，却能产生显著的影响。设计者想要了解零件的几何模型与生产过程是如何对气泡的形成产生影响。Vo等人 提出了二维轴对称模型，并用ABAQUS进行了定性的预测。Lee等人在1995年首次提出了二维连续随机模型

7、，之后在1999年到2000年，发展了三维模型 。研究者也曾用达西定律的半解析解来预测糊状区的压力降低，并将其与气泡百分比相联系 。考虑到很小的形核曲率半径，气泡内有非常大的压力。通常情况下，影响机械性能的控制因素是最大气泡的半径大小，而不是气泡的百分率。文献24里的模型超越了所有早期的模型，解算了气泡的百分率和大小、形状。这些模型通过联系过程参数和最终的微观结构的关系，可用于工艺过程的优化。Bounds等人 基于求解自由表面流、热传导和凝固建立了宏观缺陷的三维模型（宏观气缩孔）。对于质量填充和枝晶间填充，他们采用了不同的流动方程。比较计算结果与由实验浇注得到的结果，达到了很好的一致，但是仅限

8、于定性的预测。他们引入了自由表面的流动、热传导与凝固的相互作用。但是求解Navier-Stokes方程比较复杂，形核潜在的物理因素、特别是边界条件（包括模具中气体的影响）没有明确指出，模型的求解范围也不容易取得。Sabau 与Viswanathan和Pequet等人考虑了气体溶解度的影响。他们都认为气体在熔体中的局部扩散速度特别快。通过大量的实验与理论分析，Lee等人证明，铝合金中气泡长大时，氢从过饱和金属液向气泡中扩散是一个速率控制因子。Sabau与Viswanathan 提出了铝合金中预测气缩孔的三维模型，它基于Navier-Stokes方程，氢质量守恒方程，气泡形核与长大模型的求解。试图

9、求解整个液相区域和糊状区的压力控制方程，导致了内在的模型数值求解的困难：液相区的压力降低（通常低于1Pa）,明显低于糊状区的压力降低（通常是几千帕）。他们计算了不同结构的A356铝合金盘形铸件的凝固，并比较了数值结果与试验结果，总体趋势一致。铸造合金构件的最终质量，取决于在时间和空间尺度上发生的许多现象。许多研究者通过运用耦合微宏观尺度参数 或计算晶粒的演变来预测凝固的微观结构，以及预测它们对宏观偏析的影响 。在宏观的尺度上，熔体合金以相对低的温度充型，充型与冷却的时候会形成氧化物类的缺陷。为了模拟这个过程，对Navier-Stokes方程与能量方程进行耦合求解，计算液相表面的位置。许多方法以

10、前曾被应用，如Welch与Harlow提出的MAC技术 ，Hirt与Nichols提出的VOF方法等 。Kothe等人对不同的方法进行了总结归纳 。充型过程要用数十秒的时间，而且模具以米为尺度。尽管截面宽度与特征长度是以毫米为尺度的。一旦充型，必须对热传导进行解算。通常用连续性方法（如FEM）执行。需要计算几秒到几千秒的时间。铸件凝固时，晶粒形核，长大成枝晶。晶粒的枝晶长大取决于溶质。它以几十微米级进行扩散，产生微观偏析。缺陷的形成，如气泡的形核与长大，与宏观的热传导、动量传递、微观扩散相互作用。文献32运用多尺度模型来模拟代表不同宏观现象的物理现象。这些宏观现象控制了凝固过程中微米级的偏析到

11、米级的热传导与流体流动中气泡的产生。首先提出了模型的理论，然后描述了实验来验证计算。用模型计算的结果与实验结果比较，并阐述了计算最大气泡半径的方法。预测的精度也是一个重要问题，有限元网格好于有限差分网格的一个方面是：需要更少的单元与节点描述几何模型和复杂零件的凝固，更少的糊状区节点用于计算压力下降和微孔形成，而有限差分网格计算节点数量庞大而且耗时多。为保持有限元网格的优点，并提高精度。Pequet等人提出了糊状区提纯技术。当糊状区移动到粗有限元网格时，提纯的网格单元被激活，并且控制方程仅在这些单元上求解。只求解糊状区，所以要在糊状区限定边界条件。另外，许多数学模型试图通过现象去耦或假设单一因素

12、控制物理过程，来简化复杂的过程，这对于测试不同过程对整体的重要性是一个有用的办法。2.2 模型的分类Lee等人对Al-Si合金的气孔进行了广泛的研究，回顾了铝合金微孔形成的所有模型 ，气孔的预测模型可以被分为五类：（1）解析模型；（2）判据函数模型；（3）气体扩散模型；（4）连续随机模型；（5）达西流动模型。这些模型的数学复杂性不同：简单的类型，解析算法，需要最少的计算；高度复杂的计算模型，随机连续扩散模型。没有哪个模型可以预测所有生产条件下的微孔的形成。从总体上讲，综合适用性和预测质量是与计算成本成比例的。例如：解析算法模型，是收缩驱动气泡的增长模型，仅适用于高度物理约束的情况，这种情况在工

13、厂很少见。随机连续模型，可预测在一定温度梯度和扩散梯度下的气泡形核与长大，并预测凝固结构的演变。但是仅在很小的体积范围内（体积控制在很小的范围，是因为计算要求高。因此不能适用于工厂铸造生产的大尺度要求）。解析解与判据函数受它们应用的限制 ，达西定律与质量守恒求解复杂 。通过耦合已存在的微观气孔模型与热传导模型而形成多尺度模型，它可以优化模型的预测能力，降低计算成本。验证了这种耦合算法：三相系统中气泡长大扩散控制模型耦合微观枝晶臂距气泡模型和商业软件中的宏观模型。然后把该模型应用于一系列不同形状和铸造条件的L169铝合金中，并用光学图像分析仪分析铸件的微结构。用X射线照相术对两个样本进行检测，三

14、维数据给出了气泡的形状，并可与预测结果比较。阐明了收缩与气泡的关系。2.2.1 解析模型Walther等人 建立了流体流过一个长圆管，内部呈辐射状凝固的过程中形成气泡的解析模型。他们认为充型金属将以持续减小的凝固半径流过有限长度的圆管。最后，圆管的半径变的非常小，以致压力显著降低，在流管中成穴，形成气泡。T.S.Piwonka与M.C.Flemings 在Walther等人的解析模型基础之上提出了纯金属或接近纯金属的凝固过程中气泡预测的解析模型，研究了铝合金中气孔的形成。Piwonka与Flemings认为，圆管可被看作是一束有n个半径为r的圆筒，它们包含在长度为L半径为R的圆管内。n等于初晶

15、臂距的平方。在这种假设的前提下，弯曲度系数 来补偿弯曲的流体通道，得到是周围的压力， 是末端的压力（圆管末端气孔形成处）， 是常数， 是热扩散系数， 是动力粘度。上式是Hagen-Poiseuille方程在流体流过一束直的、平行的毛细管，将收缩率取代速度常数项的情况 。采用此模型，气泡的大小可被预测。为了使预测结果更接近实际，他们引入了液-固表面张力 和溶解气体对压力的影响。应用杠杆定律，把气相与液固相分割，计算了气泡内的压力为 ， 其中是液体内压力。预测了大气泡，与实验结果吻合较好。Piwonka与Flemings把杠杆定律作为铝合金中的气体守恒方程。这个模型，采用了两个假设，假设温度梯度和

16、凝固速度均为常数。在这两个假设的前提下，才能应用解析模型控制流体流过长流管的辐射状凝固过程。为了满足这些条件，需要满足充型速度为常数，这是非常不现实的。因此，解析模型不能用于复杂的工业铸造需要。2.2.2 判据函数模型最常用的判据函数模型时是新山判据 。新山判据认为气缩孔的形成主要是由于凝固收缩，糊状区静止单向的压力降低可进行计算：是温度梯度，冷却率 ， 为压力的降低。温度梯度减小，冷却率增大，都使气泡形成的可能性增大。基于判据函数（Niyama, KC , FRN ）的气缩孔模型，其应用受定性预测的限制。新山判据与其它的判据函数，使依赖于局部的温度参数，而实际上，收缩同时依赖于其它部分的凝固。新山判据用低碳钢铸造中比较精确，但是Spittle等人 最近的研究工作表明，新山判据在Al-Si合金中并不适用。2.2.3 气体扩散模型通过考虑凝固收缩与气体析出两种因素对气缩孔形成的驱动作用，Lee与Hunt 得出结论，对于小的糊状区，气孔的形成主要受气体扩散与限制