一种球墨铸铁铸锭石墨球尺寸数值预测的方法技术

本发明专利技术涉及一种球墨铸铁铸锭石墨球尺寸数值预测的方法，是为了解决现有技术中仅针对中心截面的二维计算由于没有考虑壁面效应，从而影响过冷度计算；在三维方向计算温度场和速度场会增加计算量、大大延长计算时间的缺点而提出的，包括：对铸造系统进行网格划分；对所有网格计算能量守恒方程，获得三维方向铸造系统温度场分布；对于铸锭中平行于重力方向的中心截面网格计算动量守恒方程，获得该计算网格中的金属液流动速度；对于铸锭网格和铸锭中平行于重力方向的中心截面网格，计算石墨球尺寸；重复上述步骤直至所有铸锭网格和铸锭中平行于重力方向的中心截面网格的温度均小于共晶线温度。本发明专利技术适用于各类尺寸的砂型和金属型中石墨球尺寸的预测。

A Method for Predicting the Size of Graphite Spheres in Ductile Iron Ingot

The present invention relates to a method for numerical prediction of graphite sphere size of nodular cast iron ingots, which aims to solve the problem that the calculation of supercooling degree is affected by the neglect of wall effect in the two-dimensional calculation of the central section in the prior art; the calculation of temperature field and velocity field in three-dimensional direction will increase the calculation amount and greatly prolong the calculation time. The disadvantages are as follows: meshing the casting system; calculating the energy conservation equation for all grids to obtain the temperature field distribution of the three-dimensional casting system; calculating the momentum conservation equation for the center section grid parallel to the gravity direction in the ingot to obtain the flow velocity of the liquid metal in the grid; and calculating the energy conservation equation for all grids to obtain the temperature field distribution of the three-dimensional casting system. Graphite sphere size is calculated by calculating the grid of center section parallel to gravity direction in ingot and ingot. Repeat the above steps until the temperature of all grid of ingot and grid of center section parallel to gravity direction in ingot is lower than that of eutectic line. The invention is suitable for predicting the size of graphite spheres in sand and metal moulds of various sizes.

【技术实现步骤摘要】
一种球墨铸铁铸锭石墨球尺寸数值预测的方法
本专利技术涉及球墨铸铁铸锭仿真领域，具体涉及一种球墨铸铁铸锭石墨球尺寸数值预测的方法。
技术介绍
铸铁中球状石墨的形成可以有效提高铸铁的塑性、韧性和强度。含有大量石墨球的铸铁被称为球墨铸铁，基于其优异的力学性能，该材料被成功应用于一些受力复杂或对强度、韧性、耐磨性有较高要求的零件。球墨铸铁铸锭是一种母材，需经塑性加工得到各种用途。工业上可以用于球墨铸铁管、球墨铸铁井盖、船用气缸套的加工。为了提高球墨铸铁铸锭的综合力学性能，需要对石墨球尺寸进行控制。球化率高(石墨球个数多)、石墨球细小是提高球墨铸铁铸锭力学性能的关键。石墨球产生于铸锭凝固过程，凝固过程中的过冷度是石墨球形核和生长的驱动力。铸锭凝固是一个复杂的物理过程，涉及多个尺度物理现象的相互叠交：宏观尺度动量、热量、质量传输，介观尺度晶粒组织形成，微观尺度溶质扩散，不同尺度的物理现象之间相互作用、相互影响。采用实验方法研究铸锭凝固过程中石墨球尺寸不仅浪费了人力和财力，而且无法获得不同凝固条件下石墨球尺寸的变化特点，无法获得凝固过程中石墨球尺寸随时间变化曲线。更重要的是由于铸锭尺寸较大，每一次浇注成型都会耗费大量能源，因此大量的实验研究必然造成能源浪费，使环境受到破坏和污染。所以利用计算机模拟技术预测球墨铸铁铸锭凝固过程中石墨球尺寸变化，是控制和提升球墨铸铁铸锭产品质量的一个行之有效的手段。石墨球的形成由形核和长大两阶段组成。过冷度是形核驱动力，因为它反映了液相和固相自由能的相差程度。过冷度越大，自由能相差程度就越大，球状石墨更易形核。过冷度也是生长驱动力，过冷度增大表明固液界面前沿的碳元素成分梯度增大，大的成分梯度下碳元素从液相中向石墨球中的扩散通量增加，从而提高石墨球的生长速率。准确预测铸锭凝固过程温度场从而获得温度随时间变化曲线是预测石墨球尺寸的关键。目前针对铸锭的研究多采用二维计算或三维计算。其中仅针对中心截面的二维计算由于没有考虑壁面效应，会低估冷速，从而影响过冷度计算；在三维方向计算温度场和速度场会增加计算量、大大延长计算时间。
技术实现思路
本专利技术的目的是为了解决现有技术中仅针对中心截面的二维计算由于没有考虑壁面效应，会低估冷速，从而影响过冷度计算；在三维方向计算温度场和速度场会增加计算量、大大延长计算时间缺点，而提出一种球墨铸铁铸锭石墨球尺寸数值预测的方法，包括：步骤一、对铸造系统进行网格划分；步骤二、对所有网格计算能量守恒方程，获得三维方向铸造系统温度场分布；步骤三、对于铸锭中平行于重力方向的中心截面网格计算动量守恒方程，获得该计算网格中的金属液流动速度；步骤四、对于铸锭网格和铸锭中平行于重力方向的中心截面网格，计算石墨球尺寸；步骤五、重复步骤二至四，直至所有铸锭网格和铸锭中平行于重力方向的中心截面网格的温度Tin均小于共晶线温度TE。本专利技术具有以下有益效果：本专利技术的方法采用了瞬时形核二维石墨生长模型，考虑了三维温度场计算，更加准确地预测了凝固过程中石墨球尺寸随时间变化特点，解决了目前石墨球尺寸数值预测中不考虑三维能量传输、不考虑石墨连续生长特点的问题，为分析球墨铸铁中石墨球形成供了数据参考。附图说明图1为本专利技术一个实施例的铸造系统三维示意图；图2(a)为本专利技术一个实施例的铸锭三维示意图；图2(b)为本专利技术一个实施例的铸锭中心截面图，图中P1和P2是在中心界面上选取的两个相同高度的点；图3(a)为本专利技术一个实施例中铸锭温度场二维模拟、液体流动仅针对中心截面二维模拟时的P1点和P2点的冷却曲线；图3(b)为本专利技术一个实施例中铸锭温度场三维模拟、液体流动仅针对中心截面二维模拟时P1点和P2点的冷却曲线；图3(a)和图3(b)的横轴均为时间，纵轴均为温度；图4(a)为本专利技术一个实施例中铸锭温度场二维模拟、液体流动仅针对中心截面二维模拟时4805s时刻铸锭中心截面温度场分布图；图4(b)为本专利技术一个实施例中铸锭温度场二维模拟、液体流动仅针对中心截面二维模拟时和12005s时刻铸锭中心截面温度场分布图；图5(a)为本专利技术一个实施例中铸锭温度场三维模拟、液体流动仅针对中心截面二维模拟时4805s时刻铸锭中心截面温度场分布图；图5(b)为本专利技术一个实施例中铸锭温度场三维模拟、液体流动仅针对中心截面二维模拟时12005s时刻铸锭中心截面温度场分布图；图6(a)为铸锭温度场二维模拟、液体流动仅针对中心截面二维模拟时P1点和P2点处石墨球半径随时间变化曲线图；图6(b)为铸锭温度场三维模拟、液体流动仅针对中心截面二维模拟时P1点和P2点处石墨球半径随时间变化曲线图；图6(a)和图6(b)的横轴均为时间，纵轴均为石墨球半径；图7为本专利技术一个实施例的球墨铸铁铸锭石墨球尺寸数值预测的方法的流程图。具体实施方式具体实施方式一：本实施方式的球墨铸铁铸锭石墨球尺寸数值预测的方法，如图7所示，包括：步骤一、对铸造系统进行网格划分。其中X轴、Y轴、Z轴可以为相互正交的任意坐标轴。X轴、Y轴、Z轴的选取可以依照实际情况而定。坐标系的选择不同不会影响预测结果。步骤二、对所有网格计算能量守恒方程，获得三维方向铸造系统温度场分布。步骤三、对于铸锭中平行于重力方向的中心截面网格计算动量守恒方程，获得该计算网格中的金属液流动速度。步骤四、对于铸锭网格和铸锭中平行于重力方向的中心截面网格，计算石墨球尺寸。步骤五、重复步骤二至四，直至所有铸锭网格和铸锭中平行于重力方向的中心截面网格的温度Tin均小于共晶线温度TE。本专利技术所作出的改进是基于最新理论而提出的。最新的理论研究表明，石墨片可以通过二维搭接方式进行长大从而长成球状石墨。相比于基于螺旋位错石墨长大模型，二维搭接石墨长大模型可以再现凝固过程中球状石墨连续生长状态。因此本专利技术的球墨铸铁铸锭石墨球尺寸数值预测的方法，采用更符合物理实际的瞬时形核的二维搭接石墨长大模型，同时采用三维模型计算能量传输，二维模型计算动量传输，研究冷却速度与石墨球尺寸之间的关系，这在明晰不同铸造工艺下球墨铸铁铸锭石墨球尺寸变化特点方面具有重要意义。具体实施方式二：本实施方式与具体实施方式一不同的是：步骤一具体为：对X米×Y米×Z米的铸造系统进行宏观尺度网格剖分，X方向、Y方向和Z方向分别采用Δx、Δy、Δz作为网格剖分步长，Δx米＝Δy米＝Δz米，Δx、Δy和Δz的取值范围1×10-3米～4×10-3米，计算网格的标号为(i，j，k)char，其中i，j和k均为整数，i的取值范围是1～L，j的取值范围是1～M，k的取值范围是1～N，下角标char＝2表示铸型网格，下角标char＝0表示铸锭网格，下角标char＝21表示铸锭中平行于重力方向的中心截面网格，下角标char＝4、5、6和7分别表示内冷铁网格、冒口套网格、保温材料网格和绝热材料网格；铸造系统在X轴、Y轴、Z轴方向上的最小值分别为Xmin、Ymin、Zmin(米)，在X轴、Y轴、Z轴方向上的最大值分别为Xmax、Ymax、Zmax(米)。其它步骤及参数与具体实施方式一相同。具体实施方式三：本实施方式与具体实施方式一或二不同的是：步骤二中，针对所有计算网格(i，j，k)char只计算能量守恒方程，即计算温度场，获得三维方向铸造系统温度场分布。步骤二具体为：若计本文档来自技高网...

【技术保护点】
1.一种球墨铸铁铸锭石墨球尺寸数值预测的方法，其特征在于，包括：步骤一、对铸造系统进行网格划分；步骤二、对所有网格计算能量守恒方程，获得三维方向铸造系统温度场分布；步骤三、对于铸锭中平行于重力方向的中心截面网格计算动量守恒方程，获得该计算网格中的金属液流动速度；步骤四、对于铸锭网格和铸锭中平行于重力方向的中心截面网格，计算石墨球尺寸；步骤五、重复步骤二至四，直至所有铸锭网格和铸锭中平行于重力方向的中心截面网格的温度Tin均小于共晶线温度TE。

【技术特征摘要】
1.一种球墨铸铁铸锭石墨球尺寸数值预测的方法，其特征在于，包括：步骤一、对铸造系统进行网格划分；步骤二、对所有网格计算能量守恒方程，获得三维方向铸造系统温度场分布；步骤三、对于铸锭中平行于重力方向的中心截面网格计算动量守恒方程，获得该计算网格中的金属液流动速度；步骤四、对于铸锭网格和铸锭中平行于重力方向的中心截面网格，计算石墨球尺寸；步骤五、重复步骤二至四，直至所有铸锭网格和铸锭中平行于重力方向的中心截面网格的温度Tin均小于共晶线温度TE。2.根据权利要求1所述的球墨铸铁铸锭石墨球尺寸数值预测的方法，其特征在于，步骤一具体为：对X米×Y米×Z米的铸造系统进行宏观尺度网格剖分，X方向、Y方向和Z方向分别采用Δx、Δy、Δz作为网格剖分步长，Δx米＝Δy米＝Δz米，Δx、Δy和Δz的取值范围为1×10-3米～4×10-3米，计算网格的标号为(i，j，k)char，其中i，j和k均为整数，i的取值范围是1～L，j的取值范围是1～M，k的取值范围是1～N，下角标char＝2表示铸型网格，下角标char＝0表示铸锭网格，下角标char＝21表示铸锭中平行于重力方向的中心截面网格；铸造系统在X轴、Y轴、Z轴方向上的最小值分别为Xmin、Ymin、Zmin，在X轴、Y轴、Z轴方向上的最大值分别为Xmax、Ymax、Zmax。3.根据权利要求2所述的球墨铸铁铸锭石墨球尺寸数值预测的方法，其特征在于，步骤二具体为：若计算网格(i，j，k)char的下角标char＝0，则表明该计算网格为铸锭网格但不为中心截面网格，采用下式计算温度场：如果计算网格(i，j，k)char的下角标char＝21，则表明该计算网格为铸锭网格且为中心截面网格，采用下式计算温度场：如果计算网格(i，j，k)char的下角标char不为0且不为21，表明该网格不为铸锭网格，采用下式计算温度场：[H]m-char＝ρm-charcm-charTm-char其中下角标m-char表示非铸锭材质；cm-char为比热，ρm-char...

【专利技术属性】
技术研发人员：刘东戎，赵红晨，朱泓宇，赵思聪，郭二军，
申请(专利权)人：哈尔滨理工大学，
类型：发明
国别省市：黑龙江,23