资源描述：
《计算材料学论文》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、“L”形AZ91铸造镁合金凝固过程数值模拟李宾112161366摘要：以“L”形镁合金铸件为例，对其凝固过程进行了合理的假设和简化,建立了凝固过程的二维非稳态温度场的计算数学模型；并运用有限差分方法对模型进行离散，得到显示差分格式。据此，利用Matlab编制了计算机程序，获得了铸件在凝固过程中的温度分布规律，数值模拟能够反映铸件冷却过程温度场的动态变化，表明凝固潜热的影响，也为提高铸件质量，确定凝固时间和优化工艺参数提供了参考。关键词：AZ91镁合金；铸件；潜热处理；凝固过程；温度场；有限差分；数值模拟铸造是获得机械产詁毛坯的主要方法之一，对国
2、民经济的发展起着十分重要的作用。据有关资料表明，铸件在机床、内燃机、重型机器中占70%-90%,农业机械屮ilf40%-70%o此外，铸造在冶金工业、航空、航天、船舶等重要产业部门都发挥着重要作用⑴。然而，传统的铸件生产不仅周期较长、浪费大，而且产品质量也难以保证。而在实际生产之前对铸件在浇注、凝固过程中可能产生的缺陷采用计算机数值模拟的方法进行有效的预测，从而在浇注前就采取有效对策，以确保铸件的质量。由于镁合金具冇较高的比强度、良好的减震性和切削加工性及尺寸的稳定性等，这些优良的特性使它成为非常重要的现代工业材料。目而镁合金压铸件被广泛应用于
3、汽车、航空航天和计算机制造业等各个领域以。而镁合金应用的关键的一个环节是在材料制备过程屮的金属凝固过程。因此可以通过数值模拟技术，凝固理论直接与镁合金铸造工艺联系起来，推动铸造工艺技术的发展，具冇着重耍的意义。本文以“L”形镁合金铸件为例，进行温度场的数值模拟，模拟结杲反映了铸件冷却过程温度场的动态变化。1模型的建立1.1研究对彖研究对象为L形镁合金铸件及其砂型，如图1(a)o砂型壁厚为120mm铸件厚度为50mm,横截面如图1(b)所示。本文研究边界温度为定值的情况，镁合金初始温度为其液相线温度596°C,砂型恒为25°Co(a)镁合金铸件及
4、其砂型(b)铸件横截而图图1镁合金铸件及其砂型示意图1.2铸件温度场的控制方程人们根据计算传热学、金属凝固理论以及能量守恒原理等理论，建立了铸件凝固过程三维温度场数值模拟的数学模型。这个模型的核心内容可以简单归纳为以下两个方而⑷：(1)三维不稳定导热偏微分方程的建立和求解；(2)单值条件的确定。铸件单元三维不稳定导热偏微分方程的直角坐标描述为：a?d2T+式中x,y,z——坐标轴方向的长度,mT温度,Kt时间,sp——密度,kg/n?2——热导率,W/(mK)Cp——比热,kJ/(kg・K)fs——固相率,％L——合金的结晶潜热，KJ/Kg考虑
5、到AZ91镁合金的结品温度范围较宽，本文采用等价比热容法对潜热作合理的处理⑴。其推导如Fo(2)固相率与温度密切相关，则器可表示为券畔，式⑴可表示曲(T•C—卩dt其中:式屮：7；为液相线温度；心为同相线温度；等效比热法是将潜热项折合成比热容,以等效比热容久"代替久o在实际计算中假设潜热均匀释放，则4=〃⑺一T》C；=Cp+L°=Cp+L/(Tl-TJ1.3基本假设和处理影响铸件凝固过程的因素非常多,在求解中若要把所有复杂因素都考虑进去是不现实的。因此在铸件凝固过程复杂的实际条件下，需做一些合理的简化和假设：(1)认为液体金属在瞬时充满铸型后开
6、始凝固。⑵忽略钢水在凝固过程中的相对流动,即只考虑凝固过程中钢液的对流传热。(1)由于铸模在钢液凝固过程屮的温度变化不是很大，因而将它的热物性参数看成是不随温度变化的常数。(2)不考虑合金的过冷，假定凝同是从给出的液相线温度开始,
7、古
8、相线温度结束。1.4热物性参数铸造模拟屮涉及到的最重要的参数就是热物性参数，因此准确的热物性参数是获得准确模拟结果的必要条件。在通常情况下，铸件与铸型的密度、比热、导热系数、换热系数等热物性值不是一恒定的值，I何随温度不同冇所变化，但数值变化不是很大，在模拟计算时均将系统的热物性值作为常数处理。因为这样不仅降低对
9、计算机内存的要求，大大提高模拟运算速度，而且对铸件凝固推移进程、最后凝固区域位置等相对值几乎没有影响。(1)铸件AZ91的热物性参数：表1AZ91镁合金热物性值网液相线温度596°C固相线温度468°C液态合金密度1650Kg/m3固态合金密度1750Kg/m3液态导热系数80W/(K-m)固态导热系数105W/(Km)液态比热1.35KJ/(Kg・K)固态比热1.20KJ/(Kg•K)凝固潜热370KJ/Kg(2)咲喃树脂砂型的热物性参数:⑸咲喃树脂砂的导热系数为1.01W/(mk),密度为1550kg/n?,比热为1.08kJ/(kg.K)
10、o木文屮为方便计算，假设砂型保持恒温。1.5初始条件及边界条件的确立⑺1)初始条件初始条件就是在过程开始的时刻，即t=0时整个温度场内部的温度分布。它
