CN1198971A - 在造型过程中预测湿砂的充填是否适当的方法 - Google Patents

在造型过程中预测湿砂的充填是否适当的方法 Download PDF

Info

Publication number
CN1198971A
CN1198971A CN98104153A CN98104153A CN1198971A CN 1198971 A CN1198971 A CN 1198971A CN 98104153 A CN98104153 A CN 98104153A CN 98104153 A CN98104153 A CN 98104153A CN 1198971 A CN1198971 A CN 1198971A
Authority
CN
China
Prior art keywords
sand
sand grains
force
damp
porosity
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
Application number
CN98104153A
Other languages
English (en)
Other versions
CN1108208C (zh
Inventor
泰育牧野
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sintokogio Ltd
Original Assignee
Sintokogio Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sintokogio Ltd filed Critical Sintokogio Ltd
Publication of CN1198971A publication Critical patent/CN1198971A/zh
Application granted granted Critical
Publication of CN1108208C publication Critical patent/CN1108208C/zh
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B22CASTING; POWDER METALLURGY
    • B22CFOUNDRY MOULDING
    • B22C19/00Components or accessories for moulding machines
    • B22C19/04Controlling devices specially designed for moulding machines

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Casting Devices For Molds (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Mold Materials And Core Materials (AREA)

Abstract

本发明披露了一种在造型过程中预测湿砂的充填是否适当的方法。该方法包括下列步骤:(a)分析湿砂的孔隙度;(b)分析作用在湿砂的砂粒之间的接触力;(c)分析所述砂粒周围的空气流体力;(d)根据作用在所述砂粒上的作用力计算所述砂粒的加速度,其中所述作用力包括接触力、流体力以及砂粒的重力;(e)根据所计算的加速度在砂粒运动很短的时间后分析多个运动方程式以得到所述砂粒的速率和位置;以及(f)重复上述步骤(a)、(b)、(c)、(d)和(e)直至所述砂粒停止运动。

Description

在造型过程中预测湿 砂的充填是否适当的方法
本发明涉及一种在湿砂造型过程中预测湿砂的充填是否适当的方法。
通常,在已制成砂型后才能测定湿砂的充填是否适当。因此,为了改变或者增加砂型的体积密度,人们不得不进行多次重复的造型实验,然后改进关于诸如造型工艺、造型条件以及湿砂性能的数据。因而,在某种程度上说,需要利用这些实验积累的数据才可能制出最佳的砂型。但是,这些实验积累的数据对于一种新的应用,例如需要铸造一种新的铸件(产品)或者使用一种新的造型工艺,或者使用具有新的性能的湿砂是无用的。所以,对于这样一种新应用来说,为了得到最佳条件,人们必须进行多次造型实验。这样就会花费很多时间。另外,在造型时,人们必须考虑膨润土或鲕状物(oolitiss)的影响,而这些因素的影响是不能从一般的粉料充填中预测出来。
本发明方法已经解决了这些问题。本发明的目的在于提供一种用于在诸如以提供压缩空气的方式进行造型、以喷砂的方式进行造型以及利用压实型砂的方法进行造型的造型方法中预测湿砂的充填是否适当的方法。
本发明所涉及的在湿砂造型中预测湿砂的充填是否适当的方法包括下列步骤:分析与湿砂充填程度相关的湿砂孔隙度;分析作用在湿砂的砂粒之间的接触力;分析所述砂粒周围的空气流体力;根据作用在所述砂粒上的作用力计算所述砂粒的加速度,其中所述作用力包括接触力、流体力和砂粒的重力;根据所计算的加速度在砂粒运动很短的时间后分析多个运动方程式以获得砂粒的速率和位置;重复上述的分析湿砂的孔隙度、接触力和流体力的步骤以及计算砂粒的加速度并且分析运动方程式的步骤,直至所述砂粒停止运动。
当在湿砂造型中使用空气流时,该方法还可包括这样一个步骤,即利用在上述分析湿砂孔隙度步骤中所得到的关于孔隙度的数据来分析所述空气流以得到所述空气流的速率。
在本发明中,“湿砂造型”一词通常指的是利用湿砂进行造型,其中将膨润土作为一种粘接剂。湿砂造型方法包括利用机械方式进行紧实的造型方法,例如振动造型或压实造型;还包括利用提供空气流进行造型的方法,诸如利用通入空气流、空气冲击或者吹砂的造型方法,或者将上述各种造型方法结合在一起。湿砂包括作为填充料的石英砂(或者其它砂)以及在填充料周围形成的鲕状物和膨润土的附加层。
在本发明中,“造型工艺”一词指的是根据产品的结构图制作铸件(产品)的工艺图。特别地,本发明涉及这样一种造型工艺,其中在制成砂型时可进行最佳的充填。“造型的条件”一词指的是造型过程中的条件,也就是说,例如在采用压缩空气造型的方法中的空气压力或者压实压力。“湿砂性能”一词通常包括湿砂中的含水量、湿砂的透气性以及紧实度。
图1是一个表示分析造型过程步骤的流程图。
图2表示一个能够得到砂粒接触力的砂粒模型。
图3表示可用于本发明以进行分析的一个金属砂箱和铸模的模型。
图4表示用于分析的在金属砂箱中自由下落并充填该砂箱的湿砂颗粒的一个示例。
图5示出了当从上方将空气流提供给湿砂颗粒后所述砂粒的状态。
现将参照附图描述本发明的优选实施例。图1是一个表示本发明方法步骤的流程图,该方法通过分析造型过程来预测湿砂的充填程度。现将根据该流程图描述本实施例。
在第一步中,输入关于造型方法、造型工艺、造型条件以及湿砂性能的数据。为了便于分析,将用于生产砂型的石英砂的体积分成一定数量的颗粒单元,每个颗粒单元的直径相同。根据所需要的分析的精确度来确定所述颗粒单元的数量。然后计算所述单元的直径。以同样的方式确定用于分析过程中的鲕状物和膨润土层的厚度。在本实施例中采用特殊单元方法(the distinct element method)。这种方法比其它方法具有更高的预测精度。
然后,建立用于分析孔隙度和空气流的网格(meshes)。“网格”一词表示一个计算所需要的网格(grid)。计算格点处的速率值和孔隙度值。这些网格也可用于分析空气流。
在第二步中,计算每个所述网格中的湿砂体积和每个网格的孔隙度。所述第一步和第二步共同构成了一个分析孔隙度的步骤。
在第三步中,如果使用的造型方法是提供压缩空气进行造型或吹砂的方法,并且所使用的气体是空气,那么将根据对一个考虑空气压力损失的方程式进行数学分析来得到所述空气流的速率。
第四步是一个分析接触力的步骤。这个分析步骤计算两个给定的颗粒i,j之间的距离,并且判断这两个颗粒是否彼此相互接触。如果它们相互接触,则确定两个向量。一个向量是正交向量,其方向从颗粒(i)的中心指向颗粒(j)的中心;另一个向量是正切向量,该向量与所述正交向量接逆时针方向成90°角。
如图2中所示,通过在正交和正切方向上给两个接触颗粒(特殊单元)设置虚拟弹簧和阻尼延迟器,可以得到颗粒(j)作用在颗粒(i)上的一个接触力。该接触力可通过将正交方向和正切方向的接触力合成来得到。
在第四步中,先得到正交方向的接触力。颗粒i,j在一个微小的时间段内的相对位移可通过方程式(1)来得出,在方程式(1)中使用一个与所述相对位移成正比的弹簧力增量和一个弹性弹簧系数(弹性系数)。
Δen=KnΔXn    (1)
其中,ΔXn:颗粒i,j在一个微小的时间段内的相对位移
      Δen:弹性力的增量
        Kn:与所述相对位移成正比的弹性系数(弹簧常数)。
此外,利用方程式(2)得出阻尼延迟器作用力,方程式(2)使用一个与所述相对位移速率成正比的粘性阻尼延迟器(粘度系数)。
Δdn=ηnΔXn/Δt   (2)
其中,Δdn=阻尼延迟器作用力
ηn:一个与所述相对位移的速率成正比的粘性阻尼延迟器(粘度系数)。
分别利用方程式(3)和(4)得出在一个给定的时间处作用在颗粒(i)上的颗粒(j)的正交弹簧力和阻尼延迟器作用力。
[en]t=[en]t-Δτ+Δen    (3)
[dn]t=Δdn               (4)
正切的接触力可利用方程式(5)得出。
[fn]t=[en]t+[dn]t        (5)
其中,[fn]t:一个正交接触力
因此,通过考虑所有来自其它颗粒的接触力可以得出在一个给定的时间(t)处作用在颗粒(i)上的接触力。
在第四步中,然后考虑鲕状物和膨润土的影响。换言之,因为湿砂是由诸如石英砂等的充填物以及鲕状物和膨润土的附加层构成的,所以应该按照下列公式根据与接触深度(相对位移)有关的鲕状物和膨润土层的厚度,有选择地使用相应的弹性系数值和粘度系数值:
        当δ<δh时  (6)
          kn=knh    (7)
          ηn=ηnh  (8)
其中δ:接触深度(相对位移)
    δh:鲕状物和膨润土层的厚度
      当δh<δ  时(9)
         kn=kns   (10)
         ηn=ηns (11)
其中,Knh:作用在鲕状物层和膨润土层之间的一个弹性常数
ηn:作用在鲕状物层和膨润土层之间的一个粘度系数
Kns:作用在鲕状物层和膨润土层以及石英砂砂粒之间的一个弹性常数
ηns:作用在鲕状物层和膨润土层以及石英砂砂粒之间的一个粘度系数
因为一个结合力作用在本发明所用的湿砂砂粒之间,所以必须考虑这个结合力或结合强度。当正交接触力等于或小于所述结合强度时,可以认为正交接触力为零。
在第四步中,第三是得到正切接触力。假设与正交换触力相似,正切接触力的弹性力与相对位移成正比,并且其阻尼延迟器作用力也与相对位移速率成正比。在这种情况下,利用方程式(12)可得出正切接触力。
[ft]t=[et]t+[dt]t    (12)
因为砂粒之间存在滑移或者它们在型壁上滑移,所以应该利用如下的Coulomb’s定律考虑所述的滑移量:
当  |[et]t|>μ0[en]t+fcoh    (13)时
[et]t=(μ0[en]t+fcoh)sign([en]t)    (14)
[dt]t=0    (15)
当|[et]t|<μ0[en]t+fcoh    (16)时
[et]t=[et]t-Δτ+Δet    (17)
[dt]t=Δdt      (18)
其中,μ0:摩擦系数
fcoh:结合强度
Sign(z):表示一个变量(Z)的正号或负号
在第五步骤中,得到作用在所述颗粒上的多个作用力。利用方程式(9)计算这些作用力。 f d = 1 2 ρ g C D A S U i 2 - - - ( 19 )
其中,ρg:气流密度
      CD:反应系数
      AS:喷射面积
      Ui:相对速度
当利用在第三步中由分析空气流所得到的数据计算诸如提供压缩空气进行造型或者喷砂造型的一种提供空气流进行造型的方法中的这些作用力时,计算气流和颗粒的相对速度。当造型工艺不是采用施加空气流的类型时,仅计算移动砂粒的速度。
在第六步中,利用方程式(20)得到由于颗粒间的碰撞或接触所导致的加速度,所述方程式(20)利用作用在颗粒上的作用力,即接触力、反应系数和重力。 r · · = 1 m ( f c + f d ) + g - - - ( 20 )
同样,当颗粒间发生倾斜碰撞(以一定的角度进行碰撞时),会出现转动。所述转动的角加速度利用方程式(21)给出。 ω · = T c I - - - ( 21 )
其中,
r:位置向量
m:颗粒质量
fc:接触力
fd:流体力
g:重力加速度
ω:角速度
Tc:碰撞所导致的扭矩
I:转动惯量
ω:W与时间的微分
根据从上述方程式以及公式(16)和(18)中所得到的加速度,可以得到颗粒在经过一个微小时间段后的速度和位置。 ν = ν 0 + r . . Δt - - - ( 22 ) r = r 0 + ν 0 Δt + 1 2 r . . Δ t 2 - - - ( 23 ) ω = ω 0 + ω · Δt - - - ( 24 )
其中,V:速度向量
      0:目前的数据
      Δt:一微小的时间段
在第七步中,重复上述计算过程直至颗粒停止运动。
现将根据图1中的流程图详细描述所述计算过程的一个实施例。
如图3中所示,在本实施例中使用一个金属砂箱和铸模。这里所采用的造型方法是将压缩空气提供给型砂的一种提供空气流进行造型的方法。表1中列出了湿砂的物理性能和金属砂箱以及铸模的尺寸。利用这两个尺寸对本实施例进行分析。分析中所用到的计算条件已列在表2中。
在本实施例中,如下面所描述的是一种空气流型湿砂造型方法。首先,利用数值计算得到砂粒的初始状态,所述砂粒自由地落入图3中所示的金属砂箱中。所得到的初始状态如图4所示。当从上方将空气流提供给处于初始状态的砂粒时,流体力作用在这些砂粒上。于是这些砂粒开始向下移动并且压实。
利用上述条件计算所述运动。计算结果表示在图5中。在所述铸模顶部的水平上,预测在铸模间的湿砂中有充填不充分的部分。这样,假设这些铸模不能成功地移开。因此,湿砂的性能、造型条件、造型工艺以及造型方法都将改变。利用类似的计算能够得到最佳的造型条件、造型工艺以及造型方法。虽然在本实施例中的计算是在二维分析的基础上得到的,但是也可以采用三维分析进行计算。
                      表1
填充物                    Flattery(商标)
压实性[%]                Volclay(商标)
砂粒的直径[m]             2.29×10-4
密度[kg/m3]              2500
结合强度[m/s2]           3.56×10-2
回弹系数                  0.028
砂粒的成形系数            0.861
金属砂箱的尺寸[mm]        250×110×110
每个铸模的尺寸[mm]        100×35×100
                    表2
单元数量                  1000
单元直径                  3.0×10-3
膨润土层的厚度[m]         3.0×10-4
石英砂的杨氏模量[MPa]     7.7
膨润土的杨氏模量[Mpa]     0.7
气罐的压力[Mpa]           0.5
时间间隔[S]               2.0×10-6

Claims (2)

1、一种在造型过程中预测湿砂的充填是否适当的方法,该方法包括下列步骤:
(a)分析与湿砂的充填程度相关的湿砂孔隙度;
(b)分析作用在湿砂的砂粒之间的接触力;
(c)分析所述砂粒周围的空气流体力;
(d)根据作用在所述砂粒上的作用力计算所述砂粒的加速度,其中所述作用力包括所述接触力、所述流体力以及所述砂粒的重力;
(e)根据所计算的加速度在砂粒运动很短的时间后分析多个运动方程式以得到所述砂粒的速度和位置;以及
(f)重复上述步骤(a)、(b)、(c)、(d)和(e)直至所述砂粒停止运动。
2、如权利要求1所述的方法,还包括一个分析作用在所述湿砂上的空气流的步骤,以利用在所述分析湿砂孔隙度的步骤中所得到的关于所述孔隙度的数据得到所述空气流的速率。
CN98104153A 1997-01-17 1998-01-16 在造型过程中预测湿砂的充填是否适当的方法 Expired - Lifetime CN1108208C (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP01977097A JP3346715B2 (ja) 1997-01-17 1997-01-17 生砂造型の充填不良の予測方法
JP19770/97 1997-01-17
JP19770/1997 1997-01-17

Publications (2)

Publication Number Publication Date
CN1198971A true CN1198971A (zh) 1998-11-18
CN1108208C CN1108208C (zh) 2003-05-14

Family

ID=12008578

Family Applications (1)

Application Number Title Priority Date Filing Date
CN98104153A Expired - Lifetime CN1108208C (zh) 1997-01-17 1998-01-16 在造型过程中预测湿砂的充填是否适当的方法

Country Status (6)

Country Link
US (1) US6021841A (zh)
EP (1) EP0853993B1 (zh)
JP (1) JP3346715B2 (zh)
KR (1) KR100503456B1 (zh)
CN (1) CN1108208C (zh)
DE (1) DE69826535T2 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101242919B (zh) * 2005-06-15 2010-12-22 新东工业株式会社 发泡性混合物的控制方法
CN103177194A (zh) * 2013-04-19 2013-06-26 重庆大学 一种细长型金属管药剂夯实状态的离散元分析方法
CN104968452A (zh) * 2013-02-26 2015-10-07 迪帕克·乔杜里 砂优化以减少铸造废品的计算机实施系统和方法

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3400356B2 (ja) * 1998-07-01 2003-04-28 新東工業株式会社 生型造型方法およびそのシステム
JPWO2006103996A1 (ja) * 2005-03-25 2008-09-04 北陸電力株式会社 数値解析装置および数値解析プログラム
JP4569629B2 (ja) * 2005-03-28 2010-10-27 新東工業株式会社 鋳型射出造型法
US20110202327A1 (en) * 2010-02-18 2011-08-18 Jiun-Der Yu Finite Difference Particulate Fluid Flow Algorithm Based on the Level Set Projection Framework
US9857493B2 (en) 2012-12-24 2018-01-02 Sintokogio, Ltd. Method for detecting powder and powder detection device

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CH423302A (de) * 1964-12-02 1966-10-31 Fischer Ag Georg Verfahren und Vorrichtung zum Prüfen von Formsand und ähnlichen Substanzen auf ihre Verwendungsfähigkeit für die Herstellung von Giessformen
NL131169C (zh) * 1966-07-20
JPH07239322A (ja) * 1993-04-21 1995-09-12 Mas Fab Gustav Eirich 型砂の成型特性を確定する方法と装置
CH687506A5 (de) * 1993-12-23 1996-12-31 Fischer Georg Giessereianlagen Einrichtung zur Messung von Formstoffeigenschaften.

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101242919B (zh) * 2005-06-15 2010-12-22 新东工业株式会社 发泡性混合物的控制方法
CN104968452A (zh) * 2013-02-26 2015-10-07 迪帕克·乔杜里 砂优化以减少铸造废品的计算机实施系统和方法
CN104968452B (zh) * 2013-02-26 2018-08-07 迪帕克·乔杜里 砂优化以减少铸造废品的计算机实施系统和方法
CN103177194A (zh) * 2013-04-19 2013-06-26 重庆大学 一种细长型金属管药剂夯实状态的离散元分析方法
CN103177194B (zh) * 2013-04-19 2015-10-21 重庆大学 一种细长型金属管药剂夯实状态的离散元分析方法

Also Published As

Publication number Publication date
JPH10202344A (ja) 1998-08-04
DE69826535D1 (de) 2004-11-04
DE69826535T2 (de) 2005-01-27
EP0853993B1 (en) 2004-09-29
EP0853993A1 (en) 1998-07-22
JP3346715B2 (ja) 2002-11-18
US6021841A (en) 2000-02-08
KR100503456B1 (ko) 2005-09-26
KR19980070551A (ko) 1998-10-26
CN1108208C (zh) 2003-05-14

Similar Documents

Publication Publication Date Title
CN1108209C (zh) 用于生砂铸模的方法和装置
Slonaker et al. General scaling relations for locomotion in granular media
CN1108208C (zh) 在造型过程中预测湿砂的充填是否适当的方法
Finn et al. Use of a coupled explicit—implicit solver for calculating spring-back in automotive body panels
CN103323575B (zh) 一种沉积岩石的重构方法及系统
Wu et al. Research on sieving performance of flip-flow screen using two-way particles-screen panels coupling strategy
Ghosh et al. Vertical oscillation of a bed of granular ma-terial.
Dent et al. Modeling of snow flow
CN110598168B (zh) 一种物体入水的数值仿真方法及系统
Rojek et al. Discrete element modelling and simulation of sand mould manufacture for the lost foam process
CN110083937B (zh) 一种基于Navier-Stokes方程的防波堤桨叶模态分析方法
CN112329168B (zh) 一种二维振动条件下颗粒阻尼器减振效果的3d数值仿真方法
Xiong et al. Development of an unresolved CFD-DEM method for interaction simulations between large particles and fluids
JP2008168350A (ja) 粒子媒体の詰固方法及びその装置
Enomoto et al. Effects of particle characteristics on the viscous properties of granular materials in shear
CN110849726B (zh) 一种获取压实作用下路基土压实度的方法
CN110929430B (zh) 柴油机颗粒捕集器中随机填充碳化硅颗粒的模量计算方法
CN109543337B (zh) 一种离散元模拟中的细观参数标定方法
Nishida et al. Discrete element method simulation of the restitutive characteristics of a steel spherical projectile from a particulate aggregation
Jiang et al. Simulation and experimental research on extra-squeeze forming method during gradient sand molding
JP2020015068A (ja) 横穴への鋳物砂の充填性評価方法および充填方法
Salehpour-Oskouyi et al. Three-dimensional modeling of system of vibrating spherical balls using discrete element method
Zhang et al. A Numerical Simulation of Motion of Particles under the Wafer in CMP
JP2005305540A (ja) 振動による砂型造型の充填不良を予測し充填性を高めた造型を行う方法
Rojek et al. Numerical simulation of sand mould manufacture for lost foam casting process

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CX01 Expiry of patent term
CX01 Expiry of patent term

Granted publication date: 20030514