[发明专利]一种铝合金重力金属型铸造工艺优化方法

摘要

本发明涉及一种铝合金重力金属型铸造工艺优化方法，是基于光滑粒子流体动力学计算方法对平板件重力金属型铸造渣团和气体随金属液流动进行模拟，可以有效的模拟渣团及气体的运动轨迹，在铸造之前进行预测，根据计算结果进行浇注系统及排气系统的优化设计，能有效的避免渣团进入铸件内部及气体无法及时排出，有利于在实际铸造中预测、减少铸造缺陷，通过制备铝合金平板件，采用光滑粒子流体动力学方法建立数学计算模型，以VS2010为开发平台编写程序，进行计算机运行，得出预测结果，显示铝合金平板件重力金属型铸造中渣团及气体的运动轨迹和分布情况；此预测方法使用设备少，计算方法通用、合理，计算速度快，模拟结果准确。

主权项

一种铝合金重力金属型铸造工艺优化方法，其特征在于：是基于光滑粒子法对金属液中渣团、气体运动轨迹的仿真方法应用于铝合金重力金属型铸造工艺优化，具体的步骤如下：(1)预制铝合金平板件样件①制备铝合金平板件用开合式模具，模具材料为工具钢，浇注时，模具的预热温度为300℃±5℃；②熔炼获得铝合金熔液称取铝合金4kg±0.1kg，置于熔炼坩埚中，加热至720℃±5℃，并对熔液进行搅拌，采用六氯乙烷除气，随后除渣，静置后铝合金熔液温度降至690℃±5℃，待用；③重力金属型铸造过程将铝合金熔液注入开合式模具浇口，在重力作用下金属液进入模具型腔并充满型腔；④铸件冷却及出型待铸件凝固，打开开合式模具，取出平板铸件；⑤清理铸件表面用金属刷清理铸件表面、用机械切除铸件余头、用砂纸打磨铸件表面，平板铸件成型；(2)建立粒子运动轨迹预测模型①对粒子进行初始化设计首先对粒子进行初始配置，光滑粒子法是一种基于粒子的计算方法，要对模具固壁、铝合金熔液、渣团及型腔中的气体进行粒子化，确定粒子的物性参数，具体步骤如下：1)三维实体的建立及粒子化用三维建模软件建立模具固壁、铝合金熔液、渣团及型腔中气体的三维模型，采用计算机程序进行粒子化；2)确定粒子的物性参数按照粒子类型，进行质量、密度、初始速度、粘度的设置；②建立相互作用粒子的搜索方法采用链表搜索方法，确定相互作用粒子，并进行配对，为下一步计算建立基础，先把计算区域划分为多个相同大小的单元格，单元格的尺寸为0.0706mm，所有粒子均分布在各个单元格内；分别对每个粒子的支持域内与其相互作用的粒子进行搜索配对，对于三维问题，在搜索过程中，单元格按照从左至右，从下至上，从前至后的搜索顺序，每次搜索十三个单元格，在每个时间步长计算完成后，重新进行粒子配对；③计算粒子下一时刻的位置及速度金属液粒子、气体粒子和渣团粒子下一时刻的位置及速度的计算步骤如下：1)计算金属液粒子或气体粒子的加速度值；粒子i所受的力为压力、重力、粘性力、表面张力，计算公式如下：fi＝fip+fiv+fig+fitfip=-Σjmjpi+pjρj▿iwij]]>fiv=Σjmjρj(4ηiηjrij▿iwij(ηi+ηj)rij2)vij]]>fig＝g∑jmjwij对于粒子i表面张力的求解：当相互作用的粒子i和粒子j为金属液粒子与气体粒子或者为气体粒子与金属液粒子时，考虑表面张力，计算公式如下：fit=β(▿n^i)ni]]>则粒子i的动量方程表达式如下：dvidt=-Σjmjpi+pjρiρj▿iwij+Σjmjρiρj(4ηiηjrij▿iwij(ηi+ηj)rij2)vij+g+1ρiβ(▿n^i)ni]]>其中，fi表示粒子i所受的合力，fip表示粒子i受到的压力，fiv表示粒子i受到的粘性力，fig表示粒子i受到的重力，fit表示粒子i受到的表面张力，∑j表示对粒子i支持域内与其相互作用的粒子加权求和，mj表示粒子j的质量，pi pj分别表示粒子i和j的压力值，通过求解状态方程获得，ρi表示粒子i的密度，ρj表示粒子j的密度，通过连续性方程获得，rij表示粒子间的位置差，rij2表示粒子间距离的平方，wij表示光滑函数，g表示重力加速，表示光滑函数梯度，ηi表示粒子i的动力粘性系数，ηj表示粒子j的动力粘性系数，vij表示粒子i和粒子j的速度差，β表示表面张力系数，表示粒子i的单位法向量的散度，ni表示法向量，表示粒子i的加速度；粒子在力f的作用下位置及运动状态发生变化，为完成上述的计算需要求解压力，粒子i的压力求解方程如下：Pi=c02ρ0γ[(ρiρ0)γ-1]]]>其中c0表示声速，γ是一个常数，ρ0表示初始密度，Pi表示粒子i的压力，ρi表示i粒子的密度；求得ρi值，需要求出密度变化率，求解表达式如下：dρidt=ρiΣjvij▿iwijmjρj]]>其中ρi表示粒子i的密度，ρj表示粒子j的密度，vij表示粒子i与粒子j的速度差，wij表示光滑函数，表示光滑函数梯度，表示粒子i密度的变化率；粒子在力f的作用下，可以计算出粒子加速度；2)计算渣团粒子的加速度值；渣团粒子i所受的力为压力、重力、粘性力和渣团粒子间的相互作用力，计算公式如下：fi＝fip+fiv+fig+fisfip=-Σjmjpi+pjρj▿iwij]]>fiv=Σjmjρj(4ηiηjrij▿iwij(ηi+ηj)rij2)vij]]>fig＝g∑jmjwij则粒子i的动量方程表达式如下：dvidt=-Σjmjpi+pjρiρj▿iwij+Σjmjρiρj(4ηiηjrij▿iwij(ηi+ηj)rij2)vij+g+fis/ρi]]>其中，fi表示粒子i所受的合力，fip表示粒子i受到的压力，fiv粒子i受到的粘性力，fig粒子i受到的重力，fis表示渣团粒子间的相互作用力(包括渣团间的阻尼力和排斥力)，∑j表示对粒子i支持域内与其相互作用的粒子加权求和，mj表示粒子j的质量，pi pj分别表示粒子i和j的压力值，通过求解状态方程获得，ρi表示粒子i的密度，ρj表示粒子j的密度，通过求解连续性方程获得，rij表示粒子间的位置差，rij2表示粒子间距离的平方，wij表示光滑函数，g表示重力加速，表示光滑函数梯度，ηi表示粒子i的动力粘性系数，ηj表示粒子j的动力粘性系数，vij表示粒子i和粒子j的速度差，表示粒子i的加速度；粒子在力f的作用下位置及运动状态发生变化，为完成上述的计算需要求解压力，粒子i的压力求解方程如下：Pi=c02ρ0γ[(ρiρ0)γ-1]]]>其中c0表示声速，γ是一个常数，ρ0表示初始密度，Pi表示粒子i的压力，ρi表示i粒子的密度；求得ρi值，需要求出密度变化率，求解表达式如下：dρidt=ρiΣjvij▿iwijmjρj]]>其中ρi表示粒子i的密度，ρj表示粒子j的密度，vij表示粒子i与粒子j的速度差，wij表示光滑函数，表示光滑函数梯度，表示粒子i密度的变化率；粒子在力f的作用下，可以计算出粒子加速度；3)经过一个时间步长后，对粒子的速度、位置进行更新在经过一个时间步长后，对金属液粒子、渣团粒子及气体粒子进行速度、位置更新，具体步骤如下：首先，获得作用于各个金属液粒子、渣团粒子和气体粒子的合力，进而求得加速度；其次，金属液粒子、渣团粒子和气体粒子在n+Δt时刻的速度值等于n‑Δt时刻的速度值加上加速度乘以两倍的时间步长；金属液粒子、渣团粒子和气体粒子在n+Δt时刻的位置值等于n时刻的位置值加上其在n时刻的速度值乘以时间步长再加上n时刻的加速度乘以时间步长的平方乘以二分之一，其中n表示当前计算时刻，Δt表示时间步长；一个时间步长计算完成后，重新进行粒子的搜索配对和粒子速度、位置的计算，直至充满型腔，并得到渣团粒子、气体粒子的运动轨迹；(3)预测结果数值模拟结果表明有渣团进入铸件内部，且铸件出现浇不足、气体难以排除的现象，根据模拟结果，修改浇注系统及排气系统的设计，经过再次模拟计算，渣团集中到集渣包中，未进入到铸件中，并且气体能够及时排出，没有浇不足现象。