CN105598379A - 一种基于缩松缺陷预测的回转体薄壁铸件砂型铸造过程中冒口设计方法 - Google Patents
一种基于缩松缺陷预测的回转体薄壁铸件砂型铸造过程中冒口设计方法 Download PDFInfo
- Publication number
- CN105598379A CN105598379A CN201610169257.9A CN201610169257A CN105598379A CN 105598379 A CN105598379 A CN 105598379A CN 201610169257 A CN201610169257 A CN 201610169257A CN 105598379 A CN105598379 A CN 105598379A
- Authority
- CN
- China
- Prior art keywords
- sand
- casting
- pouring
- alloy
- sand mold
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B22—CASTING; POWDER METALLURGY
- B22C—FOUNDRY MOULDING
- B22C9/00—Moulds or cores; Moulding processes
- B22C9/08—Features with respect to supply of molten metal, e.g. ingates, circular gates, skim gates
- B22C9/088—Feeder heads
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Molds, Cores, And Manufacturing Methods Thereof (AREA)
Abstract
一种基于缩松缺陷预测的回转体薄壁铸件砂型铸造过程中冒口设计方法,涉及回转体薄壁铸件砂型铸造技术。本发明是要解决现有方法无法精确设计冒口尺寸的技术问题。本发明采用物理实验的方法确定合金/砂型界面换热系数,同时配合砂型浇注圆圈形状铸件的实验手段对界面换热系数的实用性进行实验验证,基于更为准确的缩松缺陷预测设计冒口。合金/砂型界面换热系数的精确选择和实验验证为缩松缺陷的准确预测以及合理的冒口尺寸设计奠定了坚实的基础,解决了目前缺陷预测和冒口设计过程中涉及大量试算且缺少必要实验验证的问题,提高了冒口设计的准确度,加速了铸造工艺优化进程,缩短产品开发周期,将铸件产品品质提高20~40%。
Description
技术领域
本发明涉及一种回转体薄壁铸件砂型铸造过程中冒口的设计方法,涉及回转体薄壁铸件砂型铸造技术。
背景技术
回转体薄壁铸件(例如,汽车发动机罩为缸体薄壁铸件、发电机轮为缸体叶片薄壁件)砂型铸造成型过程中,由于铸件厚度小(薄壁件的特点是最薄处小于50mm),则糊状同时凝固倾向大,铸件顶部明冒口的补缩通道极易被枝晶骨架所阻断,无法有效地对铸件下方进行补缩,因此凝固结束后在回转体铸件下方易产生缩松缺陷。缩松缺陷是一种重要的铸造缺陷。由于降低了承载载荷面的面积,它的存在会严重消弱铸件的强度和硬度。相对于厚壁铸件而言,薄壁铸件由于厚度小,缩松缺陷表面被氧化的倾向大,提高了后续焊合过程的难度,导致产品废品率增大。
回转体薄壁铸件砂型铸造成型过程中引入冒口是消除缩松缺陷的有效手段之一。合理的冒口设计取决于是否对缩松缺陷形成程度和位置有一个清晰明确的判断。采用实验试错法无法有效地对缩松缺陷形成进行预判。大量的实验不但延长了铸造工艺研发周期,同时也浪费了能源、人力、物力和财力,提高了铸件产品的附加值,削弱了中国铸造厂家在国际市场上的竞争力。随着数值计算技术的迅速发展,数值模拟已成为和理论研究和实验技术并行发展的第三种科学研究方法。采用数值模拟技术研究铸件凝固过程中温度场分布特点从而预测缩松缺陷形成已进入实用阶段。
国际上开发出很多铸造模拟软件(例如,ProCast,ViewCast,MagamaSoft),可以实现对铸件凝固过程中缩松缺陷形成进行预测。但是开展准确预测工作的前提是需要提供准确的合金/砂型界面换热系数以及对该界面换热系数的实用性进行实验验证。目前大多数数值预测工作中采用数值试算法确定合金/砂型界面换热系数且实验验证仅局限于冷却曲线的对比研究。数值试算法更注重结果而不注重物理过程的推演,且效率低具有盲目性。同时实验验证数据不应该仅来源于铸件局部(冷却曲线为局部数据),更应该考虑铸件整体状态。这就要求所开发的一种回转体薄壁铸件砂型铸造过程冒口设计方法从上述两个方面入手,第一采用实验所提供的热流密度随时间变化曲线直接确定合金/砂型界面换热系数,第二基于整体铸件对合金/砂型界面换热系数的实用性进行实验验证。基于准确的缩松缺陷预测进行冒口设计,加速铸造工艺优化进程,不但可以缩短铸件产品的试制周期,降低成本,节能减排,还符合“高科技引领绿色铸造过程”这样一个科学发展观。
发明内容
本发明是要解决现有方法无法精确设计冒口尺寸的技术问题,从而提供一种基于缩松缺陷预测的回转体薄壁铸件砂型铸造过程中冒口设计方法。
本发明为解决上述技术问题采取的技术方案是:
本发明的一种回转体薄壁铸件砂型铸造过程中冒口设计方法是按以下步骤进行:
步骤一、设计物理模拟实验,获得砂型中不同点的升温曲线。进行温度场数值模拟,采用实验所提供的热流密度随时间变化曲线直接确定合金/砂型界面换热系数。
步骤一(1)、物理模拟实验中铸件形状和尺寸特征:铸件为长方形板,板的长度Lplane为254mm,板的厚度δplane与回转体薄壁铸件最小壁厚δmincasting和最大壁厚δmaxcasting的平均值相等,即δplane=(δmincasting+δmaxcasting)/2,板的宽度Wplane和板的厚度δplane相等,即Wplane=δplane;单位均为mm;
步骤一(2)、物理模拟实验采用砂型铸造,砂型为长方形。尺寸特征:长度为Lsand=[10mm+Lplane/2+Lplane+(20%×dminsand)],宽度为Wsand=[Wplane+2×(20%×dminsand)],厚度为δsand=[δplane+2×(20%×dminsand)],其中dminsand为回转体薄壁铸件在砂箱中距离砂型外表面距离的最小值。铸件空腔处于砂型中,距左、右、前、后、上、下砂型外表面距离分别为(10mm+Lplane/2)、(20%×dminsand)、(20%×dminsand)、(20%×dminsand)、(20%×dminsand)、(20%×dminsand)。直浇道尺寸特征:纵向放置圆柱体,直径Dver为Lplane/5,高度Hver为[δplane+(20%×dminsand)]。横浇道尺寸特征:横向放置圆柱体,直径Dhori为δplane/4。单位均为mm;
步骤一(3)、物理模拟实验中在距离铸件空腔上表面不同位置处的砂型中放置测温热电偶(TC)。热电偶的总数量为NT,且3≤NT≤10。每个热电偶TCi距离铸件空腔上表面的距离为xi(i的取值范围1~NT)。
步骤一(4)、浇注长方形板物理模拟实验铸件,获得型砂中不同点的升温曲线(温度随时间变化曲线)。
步骤一(5)、进行铸件铸造过程三维传热过程的计算机仿真。物理模拟实验所得热流密度(单位为Jm-2s-1)作为温度场计算的边界条件,其中t为时间(单位s),ρsand为型砂密度(单位kgm-3),cpsand为型砂比热(单位Jkg-1K-1),由升温曲线提供,Δt(单位s)为升温曲线中相邻两点之间的时间间隔,ΔTi为升温曲线中t+Δt时刻温度与t时刻温度的差值(单位K),ΔTi=Ti(t+Δt)-Ti(t)。由合金/砂型(合金和砂型之间)界面处热流密度Qsimu(t)=halloy-sand(Talloy-Tsand)与Qexp(t)相等,推算平均合金/砂型界面换热系数halloy-sand(Wm-2K-1),其中Talloy为合金/砂型界面处合金液温度,Tsand为合金/砂型界面处型砂温度。
步骤二、采用砂型浇注圆圈形状铸件。测量不同浇注温度下铸件长度。针对圆圈形状铸件开展三维传热过程的计算机仿真,合金/砂型界面换热系数来自步骤一。模拟所得铸件长度与实验测量对比,验证步骤一中所得界面换热系数的实用性。
步骤二(1)、铸件型腔为圆圈形,共5圈。进行Npouring次实验,Npouring≥3,获得Npouring个实验铸件,每次实验的浇注温度Tpouring-i(1≤i≤Npouring)不同。测量实验铸件的长度,即金属液在圆圈形铸件型腔内流经的距离Lexppouring-i(1≤i≤Npouring),单位mm。绘制实验铸件长度Lexppouring随浇注温度Tpouring的变化曲线。
步骤二(2)、针对圆圈形状铸件开展三维传热过程的计算机仿真,合金/砂型界面换热系数halloy-sand来自步骤一。采用与实验相同的浇注温度Tpouring-i(1≤i≤Npouring),进行Npouring次模拟实验,测量模拟所得铸件的长度,绘制模拟结果中铸件长度Lsimupouring随浇注温度Tpouring的变化曲线。
步骤二(3)、模拟所得Lsimupouring-Tpouring变化曲线与实验所得Lexppouring-Tpouring曲线进行对比,如果最大差值的绝对值>100mm,则回到步骤一(3),增加热电偶个数,重复步骤一(4)和步骤一(5)。如果最大差值的绝对值≤100mm,则证明合金/砂型界面换热系数halloy-sand的实用性,进入步骤三。
步骤三、采用步骤一中所获得的合金/砂型界面换热系数halloy-sand,针对回转体薄壁铸件砂型铸造过程进行三维传热过程的计算机仿真。采用传统铸造工艺中制定的浇冒口系统,获得缩松缺陷分布特征。等温线闭合区域为缩松缺陷产生区,根据缩松率大于5%即为缩松缺陷的产生标准,显示缩松缺陷在铸件内的分布;
步骤四、基于步骤三所获得的传统铸造工艺下回转体薄壁铸件中缩松缺陷分布特点,进行冒口设计和三维传热过程的计算机仿真。在回转体薄壁件下部引入偶数个尺寸形同的腰形暗冒口。将缩松缺陷尽可能移置暗冒口内,则完成回转体薄壁铸件砂型铸造过程中冒口的设计。
步骤四(1)、采用圆环将缩松缺陷包裹,采用热节圆环法计算冒口尺寸。热节圆环的模数Mrejie-circle=2πrrejie-circle,其中,rrejie-circle为热节圆环半径,单位为mm;腰形冒口上限模数Mmaokou-circle-max=1.15Mrejie-circke;腰形冒口上限宽度为单位为mm;腰形冒口上限体积Vmaokou-circle-max=3.57(amaokou-circle-max)3,单位为mm3;
步骤四(2)、腰形冒口的个数为Nmaokou-an,2≤Nmaokou-an≤8,Nmaokou-an为偶数。每一个腰形冒口具有相同的体积和尺寸。其中Vmaokou-circle-i为每一个冒口的体积,单位为mm3,1≤i≤Nmaokou-an。每一个腰形暗冒口的宽度amaokou-circle-i=(Vmaokou-circle-i/3.57)1/3,高度为2amaokou-circle-i,长度为2amaokou-circle-i,圆弧半径为amaokou-circle-i/2,单位均为mm。
本发明包括以下有益效果:
本发明设计了一种回转体薄壁铸件砂型铸造过程中冒口设计方法,采用物理实验的方法确定合金/砂型界面换热系数,同时配合砂型浇注圆圈形状铸件的实验手段对界面换热系数的实用性进行实验验证,基于更为准确的缩松缺陷预测设计冒口。合金/砂型界面换热系数的精确选择和实验验证为缩松缺陷的准确预测以及合理的冒口尺寸设计奠定了坚实的基础,解决了目前缺陷预测和冒口设计过程中涉及大量试算且缺少必要实验验证的问题,提高了冒口设计的准确度,加速了铸造工艺优化进程,缩短产品开发周期,将铸件产品品质提高20~40%。
本发明适用于回转体薄壁铸件砂型铸造过程。利用本发明可以更为准确的预测缩松缺陷形成位置,为冒口设计提供准确的参数指导,从多个方面辅助铸造工艺开发和优化,具有巨大市场应用潜力,广泛采用本发明方法进行回转体薄壁铸件砂型铸造过程中冒口设计,能够制造高达数亿元的产值。
附图说明
图1为本发明涉及的缸体叶片铸件三维实体图,其中,a为主视图,b为俯视图,c为立体图;
图2为本发明实施例1中物理模拟实验铸造系统示意图,其中,a为三维实体图,b为主视图,c为俯视图;
图3为本发明实施例1中物理模拟实验所得铸件;
图4为本发明实施例1中物理模拟实验所得3个热电偶处的升温曲线,图中,横坐标表示时间,纵坐标表示温度;
图5为本发明实施例1中实验所得热流密度Qexp(t)随时间变化曲线;
图6为本发明实施例1中的各参数变化曲线图,图中:
a为合金/砂型(合金和砂型之间)界面处合金液温度Talloy随时间变化曲线、b为合金/砂型(合金和砂型之间)界面处型砂温度Tsand随时间变化曲线、c为合金/砂型(合金和砂型之间)界面换热系数为810Wm-2K-1时模拟所得热流密度Qsimu(t)随时间变化曲线;
图7为本发明实施例1中圆圈形状铸件三维实体图,图中:a为正视图,b为俯视图;
图8为本发明实施例1中4个不同浇注温度1486℃、1514℃、1555℃和1588℃下实验所得铸件,图中,a为浇注温度1486℃下实验所得铸件、b为浇注温度1514℃下实验所得铸件、c为浇注温度1555℃下实验所得铸件、d为浇注温度1588℃下实验所得铸件;
图9为本发明实施例1中4个不同浇注温度1486℃、1514℃、1555℃和1588℃下模拟所得铸件内温度场分布;
a为浇注温度1486℃下模拟所得铸件内温度场分布、b为浇注温度1514℃下模拟所得铸件内温度场分布、c为浇注温度1555℃下模拟所得铸件内温度场分布、d为浇注温度1588℃下模拟所得铸件内温度场分布;
图10为本发明实施例1中圆圈形状铸件长度模拟和实验测量对比;
图11为本发明实施例1中传统铸造工艺冒口和浇道位置,图中,a为主视图,b为俯视图;
图12为本发明实施例1中传统铸造工艺下缩松缺陷分布预测图;图中,浅灰色区域为无缩松缺陷形成区;深灰色区域为缩松缺陷形成区;
图13为本发明实施例1中腰形冒口尺寸示意图;
图14为本发明实施例1中优化工艺中腰形暗冒口尺寸示意图,a为主视图,b为俯视图;从图中可见共8个腰形暗冒口,标号分为1,2,3,4,5,6,7和8;
图15为本发明实施例1中添加腰形暗冒口后缩松缺陷分布预测图,浅灰色区域为无缩松缺陷形成区;深灰色区域为缩松缺陷形成区。
具体实施方式
具体实施方式一:本实施方式所述的基于缩松缺陷预测的回转体薄壁铸件砂型铸造过程中冒口设计方法,是按以下步骤进行:
步骤一、设计物理模拟实验,获得砂型中不同点的升温曲线:
进行温度场数值模拟,采用实验所提供的热流密度随时间变化曲线直接确定合金/砂型(合金和砂型之间)界面换热系数;
步骤一(1)、物理模拟实验中铸件形状和尺寸特征:铸件为长方形板,长方形板的长度用Lplane表示,板的厚度δplane与回转体薄壁铸件最小壁厚δmincasting和最大壁厚δmaxcasting的平均值相等,即δplane=(δmincasting+δmaxcasting)/2,板的宽度Wplane和板的厚度δplane相等,即Wplane=δplane;单位均为mm;
步骤一(2)、物理模拟实验采用砂型铸造,砂型为长方形,其尺寸特征:长度为Lsand=[10mm+Lplane/2+Lplane+(20%×dminsand)],宽度为Wsand=[Wplane+2×(20%×dminsand)],厚度为δsand=[δplane+2×(20%×dminsand)],其中dminsand为回转体薄壁铸件在砂箱中距离砂型外表面距离的最小值;铸件空腔处于砂型中,距左、右、前、后、上、下砂型外表面距离分别为(10mm+Lplane/2)、(20%×dminsand)、(20%×dminsand)、(20%×dminsand)、(20%×dminsand)、(20%×dminsand);直浇道尺寸特征:纵向放置圆柱体,直径Dver为Lplane/5,高度Hver为[δplane+(20%×dminsand)];横浇道尺寸特征:横向放置圆柱体,直径Dhori为δplane/4,单位均为mm;
步骤一(3)、物理模拟实验中在距离铸件空腔上表面不同位置处的砂型中放置测温热电偶TC;热电偶的总数量为NT,且3≤NT≤10;每个热电偶TCi距离铸件空腔上表面的距离为xi,i的取值范围1~NT;(这里的i也指TCi中i);
步骤一(4)、浇注长方形板物理模拟实验铸件,获得型砂中不同点的升温曲线,所述升温曲线为温度随时间变化曲线;
步骤一(5)、进行铸件铸造过程三维传热过程的计算机仿真:物理模拟实验所得热流密度作为温度场计算的边界条件,热流密度的单位为Jm-2s-1,式中t为时间,单位s;ρsand为型砂密度,单位kgm-3;cpsand为型砂比热,单位Jkg-1K-1;由升温曲线提供,Δt为升温曲线中相邻两点之间的时间间隔,单位为s;ΔTi为升温曲线中t+Δt时刻温度与t时刻温度的差值,单位K;ΔTi=Ti(t+Δt)-Ti(t);由合金/砂型(合金和砂型之间)界面处热流密度Qsimu(t)=halloy-sand(Talloy-Tsand)与Qexp(t)相等,推算平均合金/砂型界面换热系数halloy-sand(Wm-2K-1),其中Talloy为合金/砂型界面处合金液温度,Tsand为合金/砂型界面处型砂温度;
步骤二、采用砂型浇注圆圈形状铸件:测量不同浇注温度下铸件长度,针对圆圈形状铸件开展三维传热过程的计算机仿真,合金/砂型界面换热系数来自步骤一获得的平均合金/砂型界面换热系数halloy-sand(Wm-2K-1),模拟所得铸件长度与实验测量对比,验证步骤一中所得界面换热系数的实用性:
步骤二(1)、铸件型腔为圆圈形,共5圈,进行Npouring次实验,Npouring≥3,获得Npouring个实验铸件,每次实验的浇注温度Tpouring-i(1≤i≤Npouring)不同;测量实验铸件的长度,即金属液在圆圈形铸件型腔内流经的距离Lexppouring-i(1≤i≤Npouring),单位mm;绘制实验铸件长度Lexppouring随浇注温度Tpouring的变化曲线;
步骤二(2)、针对圆圈形状铸件开展三维传热过程的计算机仿真,合金/砂型界面换热系数halloy-sand来自步骤一;采用与实验相同的浇注温度Tpouring-i(1≤i≤Npouring),进行Npouring次模拟实验,测量模拟所得铸件的长度,绘制模拟结果中铸件长度Lsimupouring随浇注温度Tpouring的变化曲线;
步骤二(3)、模拟所得Lsimupouring-Tpouring变化曲线与实验所得Lexppouring-Tpouring曲线进行对比,如果最大差值的绝对值>100mm,则回到步骤一(3),增加热电偶个数,重复步骤一(4)和步骤一(5);如果最大差值的绝对值≤100mm,则证明合金/砂型界面换热系数halloy-sand具有实用性,进入步骤三;
步骤三、采用步骤一中所获得的合金/砂型界面换热系数halloy-sand,针对回转体薄壁铸件砂型铸造过程进行三维传热过程的计算机仿真;采用传统铸造工艺中制定的浇冒口系统,获得缩松缺陷分布特征;等温线闭合区域为缩松缺陷产生区,根据缩松率大于5%即为缩松缺陷的产生标准,显示缩松缺陷在铸件内的分布;
步骤四、基于步骤三所获得的传统铸造工艺下回转体薄壁铸件中缩松缺陷分布特点,进行冒口设计和三维传热过程的计算机仿真;在回转体薄壁件下部引入偶数个尺寸形同的腰形暗冒口;将缩松缺陷尽可能移置暗冒口内,则完成回转体薄壁铸件砂型铸造过程中冒口的设计:
步骤四(1)、采用圆环将缩松缺陷包裹,采用热节圆环法计算冒口尺寸。热节圆环的模数Mrejie-circle=2πrejie-circle,其中,rrejie-circle为热节圆环半径,单位为mm;腰形冒口上限模数Mmaokou-circle-max=1.15Mrejie-circke;腰形冒口上限宽度为单位为mm;腰形冒口上限体积Vmaokou-circle-max=3.57(amaokou-circle-max)3,单位为mm3;
步骤四(2)、腰形冒口的个数为Nmaokou-an,2≤Nmaokou-an≤8,Nmaokou-an为偶数;每一个腰形冒口具有相同的体积和尺寸;其中Vmaokou-circle-i为每一个冒口的体积,单位为mm3,1≤i≤Nmaokou-an;每一个腰形暗冒口的宽度amaokou-circle-i=(Vmaokou-circle-i/3.57)1/3,高度为2amaokou-circle-i,长度为2amaokou-circle-i,圆弧半径为amaokou-circle-i/2,单位均为mm。
具体实施方式二:本实施方式所述在步骤一(1)中,长方形板的长度Lplane为254mm。其它步骤与具体实施方式一相同。
具体实施方式三:本实施方式在步骤一中,
所述回转体薄壁铸件选择缸体叶片铸件,该缸体叶片铸件采用双相不锈钢制造,双相不锈钢材质为:C:0.02wt%,Si:0.58wt%,Mn:0.88wt%,S:0.04wt%,P:0.04wt%,Ni:5.00wt%,Mo:3.00wt%,Cu:0.60wt%,N:0.25wt%,Cr:22.00wt%,余量为Fe。其它步骤与具体实施方式一或二相同。
具体实施方式四:本实施方式在步骤一中,砂型材质为树脂砂,密度ρsand为1520.0kgm-3,比热cpsand为740Jkg-1K-1。其它步骤与具体实施方式一、二或三相同。
具体实施方式五:本实施方式所述缸体叶片铸件的最小壁厚δmincasting=35mm,最大壁厚δmaxcasting=85mm。其它步骤与具体实施方式一、二、三或四相同。
实施例1:
本试验所选回转体薄壁铸件为缸体叶片铸件,该缸体叶片铸件采用双相不锈钢制造,材质为C:0.02wt%,Si:0.58wt%,Mn:0.88wt%,S:0.04wt%,P:0.04wt%,Ni:5.00wt%,Mo:3.00wt%,Cu:0.60wt%,N:0.25wt%,Cr:22.00wt%。砂型材质为树脂砂,密度ρsand为1520.0kgm-3,比热cpsand为740Jkg-1K-1。
如图1为缸体叶片铸件三维实体。该铸件最小壁厚δmincasting=35mm,最大壁厚δmaxcasting=85mm。
如图2为物理模拟实验铸造系统示意图。物理模拟实验铸件材质为双相不锈钢:C:0.02wt%,Si:0.58wt%,Mn:0.88wt%,S:0.04wt%,P:0.04wt%,Ni:5.00wt%,Mo:3.00wt%,Cu:0.60wt%,N:0.25wt%,Cr:22.00wt%。砂型材质为树脂砂,密度ρsand为1520.0kgm-3,比热cpsand为740Jkg-1K-1。实验件为长方形板,厚度δplane与缸体叶片铸件最小壁厚δmincasting和最大壁厚δmaxcasting的平均值相等,即δplane=(δmincasting+δmaxcasting)/2=60mm。板的宽度Wplane和板的厚度δplane相等,即Wplane=δplane=60mm。缸体叶片铸件在砂箱中距离砂型外表面距离的最小值dminsand=300mm。物理模拟实验采用砂型铸造,砂型为长方形。尺寸特征:长度为Lsand=[10mm+Lplane/2+Lplane+(20%×dminsand)]=451mm,宽度为Wsand=[Wplane+2×(20%×dminsand)]=180mm,厚度为δsand=[δplane+2×(20%×dminsand)]=180mm。铸件空腔处于砂型中,距左、右、前、后、上、下砂型外表面距离分别为(10mm+Lplane/2)=137mm、(20%×dminsand)=60mm、(20%×dminsand)=60mm、(20%×dminsand)=60mm、(20%×dminsand)=60mm、(20%×dminsand)=60mm。直浇道尺寸特征:纵向放置圆柱体,直径Dver为Lplane/5约为51mm,高度Hver为[δplane+(20%×dminsand)]=120mm。横浇道尺寸特征:横向放置圆柱体,直径Dhori为δplane/4=15mm。
图3为物理模拟实验所得实际铸件。图2中物理模拟实验选取三个热电偶,TC1、TC2和TC3,热电偶距离合金/砂型界面的距离分别为x1=10mm、x2=5mm和x3=10mm。图4为实验中热电偶测温结果,温度曲线为升温曲线。
图5为物理模拟实验所得热流密度Qexp(t)随时间变化曲线。图6为合金/砂型界面换热系数为810Wm-2K-1时模拟所得热流密度随时间变化曲线。合金/砂型界面处合金液温度Talloy和合金/砂型界面处型砂温度Tsand由铸造模拟软件给出,根据Qsimu(t)=halloy-sand(Talloy-Tsand)=Qexp(t)推算平均合金/砂型界面换热系数为810Wm-2K-1。
为验证界面换热系数的实用性,采用砂型铸造圆圈形状铸件,图7为圆圈形状铸件三维实体图,该圆圈形状铸件共5圈。圆圈形状铸件材质为双相不锈钢:C:0.02wt%,Si:0.58wt%,Mn:0.88wt%,S:0.04wt%,P:0.04wt%,Ni:5.00wt%,Mo:3.00wt%,Cu:0.60wt%,N:0.25wt%,Cr:22.00wt%。砂型材质为树脂砂,密度ρsand为1520.0kgm-3,比热cpsand为740Jkg-1K-1。
图8为四个不同浇注温度下(1486℃、1514℃、1555℃和1588℃)实验所得铸件的实物图。变化趋势为:随着浇注温度的升高,所得铸件长度增加。
图9为四个不同浇注温度下(1486℃、1514℃、1555℃和1588℃)模拟所得铸件内温度场分布。灰色区域为铸件型腔没有被金属液填充,所以没有进行三维传热计算。变化趋势为:随着浇注温度的升高,金属液填充铸件型腔距离增加。
图10为模拟所得铸件长度和实验测量对比。模拟计算和实验测量的最大差值的绝对值为73mm,小于100mm。
图11为缸体叶片传统铸造工艺中冒口和浇道位置。图12为传统铸造工艺下缩松缺陷分布预测结果,缸体铸件中存在缩松缺陷。
对传统工艺进行优化,引入腰形冒口。图13为腰形冒口尺寸特征示意图。在缸体叶片底部引入8个腰形暗冒口。热节圆环的半径rrejie-circle由铸造模拟软件测出,其中rrejie-circle=9.8mm。热节圆环的模数Mrejie-circle=2πrrejie-circle=61.6。腰形暗冒口上限模数Mmaokou-circle-max=1.15Mrejie-circle=71.0mm,腰形暗冒口上限宽度为腰形暗冒口上限体积为Vmaokou-circle-max=3.57(amaokou-circle-max)3=9.6×107mm3。每一个腰形暗冒口具有相同的体积和尺寸,每一个腰形暗冒口的宽度amaokou-circle-i=(Vmaokou-circle-i/3.57)1/3=1.2×107/3.57)1/3=150mm,高度为2amaokou-circle-i=300mm,长度为2amaokou-circle-i=300mm,圆弧半径为amaokou-circle-i/2=75mm。图14为缸体叶片铸件中腰形暗冒口添加位置图。由图15可见,添加腰形暗冒口后缩松缺陷分布在浇道和冒口中,被移出铸件,确保了铸件质量。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (5)
1.一种基于缩松缺陷预测的回转体薄壁铸件砂型铸造过程中冒口设计方法,其特征在于,所述方法是按以下步骤进行:
步骤一、设计物理模拟实验,获得砂型中不同点的升温曲线:
进行温度场数值模拟,采用实验所提供的热流密度随时间变化曲线直接确定合金/砂型界面换热系数;
步骤一(1)、物理模拟实验中铸件形状和尺寸特征:铸件为长方形板,长方形板的长度用Lplane表示,板的厚度δplane与回转体薄壁铸件最小壁厚δmincasting和最大壁厚δmaxcasting的平均值相等,即δplane=(δmincasting+δmaxcasting)/2,板的宽度Wplane和板的厚度δplane相等,即Wplane=δplane;单位均为mm;
步骤一(2)、物理模拟实验采用砂型铸造,砂型为长方形,其尺寸特征:长度为Lsand[10mm+Lplane/2+Lplane+(20%×dminsand)],宽度为Wsand=[Wplane+2×(20%×dminsand)],厚度为δsand=[δplane+2×(20%×dminsand)],其中dminsand为回转体薄壁铸件在砂箱中距离砂型外表面距离的最小值;铸件空腔处于砂型中,距左、右、前、后、上、下砂型外表面距离分别为(10mm+Lplane/2)、(20%×dminsand)、(20%×dminsand)、(20%×dminsand)、(20%×dminsand)、(20%×dminsand);直浇道尺寸特征:纵向放置圆柱体,直径Dver为Lplane/5,高度Hver为[δplane+(20%×dminsand)];横浇道尺寸特征:横向放置圆柱体,直径Dhori为δplane/4,单位均为mm;
步骤一(3)、物理模拟实验中在距离铸件空腔上表面不同位置处的砂型中放置测温热电偶TC;热电偶的总数量为NT,且3≤NT≤10;每个热电偶TCi距离铸件空腔上表面的距离为xi,i的取值范围1~NT;
步骤一(4)、浇注长方形板物理模拟实验铸件,获得型砂中不同点的升温曲线,所述升温曲线为温度随时间变化曲线;
步骤一(5)、进行铸件铸造过程三维传热过程的计算机仿真:物理模拟实验所得热流密度作为温度场计算的边界条件,热流密度的单位为Jm- 2s-1,式中t为时间,单位s;ρsand为型砂密度,单位kgm-3;cpsand为型砂比热,单位Jkg-1K-1;由升温曲线提供,Δt为升温曲线中相邻两点之间的时间间隔,单位为s;ΔTi为升温曲线中t+Δt时刻温度与t时刻温度的差值,单位K;ΔTi=Ti(t+Δt)-Ti(t);由合金/砂型界面处热流密度Qsimu(t)=halloy-sand(Talloy-Tsand)与Qexp(t)相等,推算平均合金/砂型界面换热系数halloy-sand(Wm-2K-1),其中Talloy为合金/砂型界面处合金液温度,Tsand为合金/砂型界面处型砂温度;
步骤二、采用砂型浇注圆圈形状铸件:测量不同浇注温度下铸件长度,针对圆圈形状铸件开展三维传热过程的计算机仿真,合金/砂型界面换热系数来自步骤一获得的平均合金/砂型界面换热系数halloy-sand,单位为Wm-2K-1,模拟所得铸件长度与实验测量对比,验证步骤一中所得界面换热系数的实用性:
步骤二(1)、铸件型腔为圆圈形,共5圈,进行Npouring次实验,Npouring≥3,获得Npouring个实验铸件,每次实验的浇注温度Tpouring-i(1≤i≤Npouring)不同;测量实验铸件的长度,即金属液在圆圈形铸件型腔内流经的距离Lexppouring-i(1≤i≤Npouring),单位mm;绘制实验铸件长度Lexppouring随浇注温度Tpouring的变化曲线;
步骤二(2)、针对圆圈形状铸件开展三维传热过程的计算机仿真,合金/砂型界面换热系数halloy-sand来自步骤一;采用与实验相同的浇注温度Tpouring-i(1≤i≤Npouring),进行Npouring次模拟实验,测量模拟所得铸件的长度,绘制模拟结果中铸件长度Lsimupouring随浇注温度Tpouring的变化曲线;
步骤二(3)、模拟所得Lsimupouring-Tpouring变化曲线与实验所得Lexppouring-Tpouring曲线进行对比,如果最大差值的绝对值>100mm,则回到步骤一(3),增加热电偶个数,重复步骤一(4)和步骤一(5);如果最大差值的绝对值≤100mm,则证明合金/砂型界面换热系数halloy-sand具有实用性,进入步骤三;
步骤三、采用步骤一中所获得的合金/砂型界面换热系数halloy-sand,针对回转体薄壁铸件砂型铸造过程进行三维传热过程的计算机仿真;采用传统铸造工艺中制定的浇冒口系统,获得缩松缺陷分布特征;等温线闭合区域为缩松缺陷产生区,根据缩松率大于5%即为缩松缺陷的产生标准,显示缩松缺陷在铸件内的分布;
步骤四、基于步骤三所获得的传统铸造工艺下回转体薄壁铸件中缩松缺陷分布特点,进行冒口设计和三维传热过程的计算机仿真;在回转体薄壁件下部引入偶数个尺寸形同的腰形暗冒口;将缩松缺陷尽可能移置暗冒口内,则完成回转体薄壁铸件砂型铸造过程中冒口的设计:
步骤四(1)、采用圆环将缩松缺陷包裹,采用热节圆环法计算冒口尺寸。热节圆环的模数Mrejie-circle=2πrrejie-circle,其中,rrejie-circle为热节圆环半径,单位为mm;腰形冒口上限模数Mmaokou-circle-max=1.15Mrejie-circle;腰形冒口上限宽度为单位为mm;腰形冒口上限体积Vmaokou-circle-max=3.57(amaokou-circle-max)3,单位为mm3;
步骤四(2)、腰形冒口的个数为Nmaokou-an,2≤Nmaokou-an≤8,Nmaokou-an为偶数;每一个腰形冒口具有相同的体积和尺寸;其中Vmaokou-circle-i为每一个冒口的体积,单位为mm3,1≤i≤Nmaokou-an;每一个腰形暗冒口的宽度amaokou-circle-i=(Vmaokou-circle-i/3.57)1/3,高度为2amaokou-circle-i,长度为2amaokou-circle-i,圆弧半径为amaokou-circle-i/2,单位均为mm。
2.根据权利要求1所述的一种基于缩松缺陷预测的回转体薄壁铸件砂型铸造过程中冒口设计方法,其特征在于,在步骤一(1)中,长方形板的长度Lplane为254mm。
3.根据权利要求1或2所述的一种基于缩松缺陷预测的回转体薄壁铸件砂型铸造过程中冒口设计方法,步骤一中,
所述回转体薄壁铸件选择缸体叶片铸件,该缸体叶片铸件采用双相不锈钢制造,双相不锈钢材质为:C:0.02wt%,Si:0.58wt%,Mn:0.88wt%,S:0.04wt%,P:0.04wt%,Ni:5.00wt%,Mo:3.00wt%,Cu:0.60wt%,N:0.25wt%,Cr:22.00wt%,余量为Fe。
4.根据权利要求3所述的一种基于缩松缺陷预测的回转体薄壁铸件砂型铸造过程中冒口设计方法,步骤一中,砂型材质为树脂砂,密度ρsand为1520.0kgm-3,比热cpsand为740Jkg- 1K-1。
5.根据权利要求4所述的一种基于缩松缺陷预测的回转体薄壁铸件砂型铸造过程中冒口设计方法,所述缸体叶片铸件的最小壁厚δmincasting=35mm,最大壁厚δmaxcasting=85mm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610169257.9A CN105598379B (zh) | 2016-03-23 | 2016-03-23 | 一种基于缩松缺陷预测的回转体薄壁铸件砂型铸造过程中冒口设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610169257.9A CN105598379B (zh) | 2016-03-23 | 2016-03-23 | 一种基于缩松缺陷预测的回转体薄壁铸件砂型铸造过程中冒口设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105598379A true CN105598379A (zh) | 2016-05-25 |
CN105598379B CN105598379B (zh) | 2017-07-14 |
Family
ID=55979015
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610169257.9A Expired - Fee Related CN105598379B (zh) | 2016-03-23 | 2016-03-23 | 一种基于缩松缺陷预测的回转体薄壁铸件砂型铸造过程中冒口设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105598379B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106825422A (zh) * | 2017-03-06 | 2017-06-13 | 江苏理工学院 | 对铸钢件冒口补缩效率以及保温冒口质量的检验方法 |
CN107844852A (zh) * | 2017-10-20 | 2018-03-27 | 太原理工大学 | 一种模拟铸钢件砂型铸造过程的缩松缺陷预测方法 |
CN108344769A (zh) * | 2018-03-23 | 2018-07-31 | 中国航发哈尔滨东安发动机有限公司 | 一种测量铸造材料换热系数的装置 |
CN111872324A (zh) * | 2020-06-23 | 2020-11-03 | 上海交通大学 | 铸件凝固模拟用参数采集方法及浇冒系统网格化设计方法 |
CN112100819A (zh) * | 2020-08-24 | 2020-12-18 | 上海交通大学 | 一种基于集成计算与数据驱动的铸造冒口设计方法 |
CN114324449A (zh) * | 2021-12-29 | 2022-04-12 | 重庆大学 | 新型铸造模具试验装置及其对缩松的控制方法 |
CN117171895A (zh) * | 2023-11-03 | 2023-12-05 | 南通海泰科特精密材料有限公司 | 一种提高solidcast模拟准确性的方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE4028868A1 (de) * | 1990-09-12 | 1992-03-19 | W & W Schenk Gmbh & Co Kg Leic | Sandgussform zur herstellung von zylinderkoepfen aus aluminium |
US7270172B1 (en) * | 2004-11-05 | 2007-09-18 | Ashland Licensing And Intellectual Property Llc | Process for casting a metal |
CN101166595A (zh) * | 2005-04-26 | 2008-04-23 | 株式会社木村铸造所 | 铸铁铸件的缩孔的预测及防止方法 |
CN104439086A (zh) * | 2014-12-08 | 2015-03-25 | 哈尔滨理工大学 | 一种非轴对称复杂铸件砂型铸造过程中冒口的设计方法 |
CN104741535A (zh) * | 2013-12-25 | 2015-07-01 | 青岛玉光精铸厂 | 模数法计算钢铸件冒口的方法 |
CN105127372A (zh) * | 2015-08-31 | 2015-12-09 | 共享铸钢有限公司 | 一种大壁厚半环形内缸铸件的冒口设计方法 |
-
2016
- 2016-03-23 CN CN201610169257.9A patent/CN105598379B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE4028868A1 (de) * | 1990-09-12 | 1992-03-19 | W & W Schenk Gmbh & Co Kg Leic | Sandgussform zur herstellung von zylinderkoepfen aus aluminium |
US7270172B1 (en) * | 2004-11-05 | 2007-09-18 | Ashland Licensing And Intellectual Property Llc | Process for casting a metal |
CN101166595A (zh) * | 2005-04-26 | 2008-04-23 | 株式会社木村铸造所 | 铸铁铸件的缩孔的预测及防止方法 |
CN104741535A (zh) * | 2013-12-25 | 2015-07-01 | 青岛玉光精铸厂 | 模数法计算钢铸件冒口的方法 |
CN104439086A (zh) * | 2014-12-08 | 2015-03-25 | 哈尔滨理工大学 | 一种非轴对称复杂铸件砂型铸造过程中冒口的设计方法 |
CN105127372A (zh) * | 2015-08-31 | 2015-12-09 | 共享铸钢有限公司 | 一种大壁厚半环形内缸铸件的冒口设计方法 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106825422A (zh) * | 2017-03-06 | 2017-06-13 | 江苏理工学院 | 对铸钢件冒口补缩效率以及保温冒口质量的检验方法 |
CN107844852A (zh) * | 2017-10-20 | 2018-03-27 | 太原理工大学 | 一种模拟铸钢件砂型铸造过程的缩松缺陷预测方法 |
CN107844852B (zh) * | 2017-10-20 | 2021-10-01 | 太原理工大学 | 一种模拟铸钢件砂型铸造过程的缩松缺陷预测方法 |
CN108344769A (zh) * | 2018-03-23 | 2018-07-31 | 中国航发哈尔滨东安发动机有限公司 | 一种测量铸造材料换热系数的装置 |
CN111872324A (zh) * | 2020-06-23 | 2020-11-03 | 上海交通大学 | 铸件凝固模拟用参数采集方法及浇冒系统网格化设计方法 |
WO2021259386A1 (zh) * | 2020-06-23 | 2021-12-30 | 上海交通大学 | 铸件凝固模拟用参数采集方法及浇冒系统网格化设计方法 |
US11638953B2 (en) | 2020-06-23 | 2023-05-02 | Shanghai Jiao Tong University | Method for collecting parameters for casting solidification simulation and gridded design method for pouring and riser system |
CN112100819A (zh) * | 2020-08-24 | 2020-12-18 | 上海交通大学 | 一种基于集成计算与数据驱动的铸造冒口设计方法 |
CN112100819B (zh) * | 2020-08-24 | 2024-02-13 | 上海交通大学 | 一种基于集成计算与数据驱动的铸造冒口设计方法 |
CN114324449A (zh) * | 2021-12-29 | 2022-04-12 | 重庆大学 | 新型铸造模具试验装置及其对缩松的控制方法 |
CN114324449B (zh) * | 2021-12-29 | 2024-07-19 | 重庆大学 | 新型铸造模具试验装置及其对缩松的控制方法 |
CN117171895A (zh) * | 2023-11-03 | 2023-12-05 | 南通海泰科特精密材料有限公司 | 一种提高solidcast模拟准确性的方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105598379B (zh) | 2017-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105598379B (zh) | 一种基于缩松缺陷预测的回转体薄壁铸件砂型铸造过程中冒口设计方法 | |
Choudhari et al. | Casting design and simulation of cover plate using AutoCAST-X software for defect minimization with experimental validation | |
Khan et al. | A comparative study of simulation software for modelling metal casting processes | |
Ravi et al. | Casting solidification analysis by modulus vector method | |
CN104439086B (zh) | 一种非轴对称复杂铸件砂型铸造过程中冒口的设计方法 | |
Strohmandl | Use of simulation to reduction of faulty products | |
Dong et al. | Determination of wax pattern die profile for investment casting of turbine blades | |
CN107909189B (zh) | 一种模拟铝合金砂型铸造过程的缩孔缺陷预测方法 | |
Ingle et al. | A literature survey of methods to study and analyze the gating system design for its effect on casting quality | |
Sultana et al. | Solidification and filling related defects analysis using casting simulation technique with experimental validation | |
CN108108529A (zh) | 一种简便测量铸造界面换热系数的反算方法 | |
Rajkumar et al. | Influence of parameters on the smart productivity of modern metal casting process: an overview | |
Shehata et al. | Computer aided foundry die-design | |
Khan et al. | Mold design and casting of an impeller using MAGMASoft | |
Banerjee et al. | Analysis of Casting Defects in Investment Casting by Simulation | |
Sutaria et al. | Feeding system design and evaluation using temperature gradient (Feed Path) Maps | |
Namchanthra et al. | A CFD investigation into molten metal flow and its solidification under gravity sand moulding in plumbing components | |
CN108080604B (zh) | 螺旋桨铸件成形凝固顺序控制装置及其控制步骤 | |
Srinivasan et al. | Tribological Studies of Shrinkage Defect and Effective Yield Upgrade of Grey Cast-Iron Castings | |
Liu et al. | Numerical Simulation on Impeller Casting Defects and Optimization | |
Kamble et al. | Optimized Design of Risering System for Casted Component by Using Web Based Online Simulation E-Tool | |
Lagad | Air entrapment analysis of casting (turbine housing) for shell moulding process using simulation technique | |
Dmitriev et al. | Design of the gating system for production of casting blanks for space device cases | |
Guo et al. | Numerical simulation of the through process of aerospace titanium alloy casting filling, solidification, and hot isostatic pressing | |
Singh et al. | Applications of optimization and inverse modeling to alloy wheel casting |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170714 Termination date: 20190323 |
|
CF01 | Termination of patent right due to non-payment of annual fee |