CN102095395B - 大块金属玻璃形成厚度预测方法 - Google Patents
大块金属玻璃形成厚度预测方法 Download PDFInfo
- Publication number
- CN102095395B CN102095395B CN201010563783A CN201010563783A CN102095395B CN 102095395 B CN102095395 B CN 102095395B CN 201010563783 A CN201010563783 A CN 201010563783A CN 201010563783 A CN201010563783 A CN 201010563783A CN 102095395 B CN102095395 B CN 102095395B
- Authority
- CN
- China
- Prior art keywords
- temperature
- metallic glass
- partiald
- quartz ampoule
- lambda
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明涉及一种大块金属玻璃形成厚度预测方法,主要针对Zr-Ti-Cu-Ni-Be系列合金的金属玻璃形成厚度,建立了一个该合金系在水淬条件下的玻璃形成厚度的解析方法,进而确定出该合金系在不同冷却条件下的玻璃形成厚度。该方法具有预测成本低、预测方法简单、预测效率高,而且适用性广,不仅适用于Zr-Ti-Cu-Ni-Be系列合金,而且适用于其它系合金的玻璃形成厚度预测。
Description
一、技术领域:
本发明主要涉及一种大块金属玻璃形成厚度的预测方法,属于金属玻璃形成能力预测技术领域。
二、背景技术:
金属玻璃是一种具有短程有序的非稳态的凝聚态物质,它具有很多结晶态金属所不具备的性能而被广泛关注,如具有很高的强度、硬度、弹性、刚性,有高饱和磁感应、低铁损等性能。但大块金属玻璃的形成是与合金系、冷却条件分不开的,通常金属玻璃的形成是在超高冷却速度条件下形成的,而这大大限制了金属玻璃的形成尺寸。Zr-Ti-Cu-Ni-Be系合金是一种可成形的大尺寸金属玻璃材料,并已开始投入到实际生产,如最新的iphone手机外壳。但以往金属玻璃的形成厚度都是通过实验方法进行测定,实验方法测定具有成本较高,周期长的缺点,而且每次实验测定的结果受实验条件影响而产生波动。
三、发明内容:
1、发明目的:
本发明提出了一种大块金属玻璃形成厚度预测方法,解决了目前对金属玻璃的形成厚度实验测定方法成本高、周期长、实验结果易产生波动等缺陷,可以预测出不同的冷却条件下各合金系金属玻璃形成的厚度。
2、技术方案:
本发明是通过以下技术方案来实现的:
一种大块金属玻璃形成厚度预测方法,所述大块金属玻璃的形状能够看成半无限大固体,其特征在于:该方法步骤如下:
(1)、确定金属玻璃内部的传热与温度场:
在金属玻璃内部,金属玻璃的温度T与时间t之间存在着一定的关系,其关系如公式(1)所示
其中
其中T为温度(℃),t为时间(s),k1为合金导热率(W/m2K),ρ1为合金的密度(kg/m3),Cp1为合金比热容(J/kgK),α1为合金的热扩散系数(m2/s),r为金属玻璃厚度的一半;
(2)、确定石英管壳型内部的传热与温度场:
在石英管壳型内部,金属玻璃的温度T与时间t之间存在的关系如公式(3)所示
并且有
其中T为温度(℃),t为时间(s),k2为导热率(W/m2K),ρ2为石英管的密度(kg/m3),Cp2为石英管的比热容(J/kgK),α2为石英管热扩散系数(m2/s),δ为壳型的厚度;
(3)、确定金属玻璃与石英管壳型之间的热量传输:
在金属玻璃和石英管壳型之间存在下面的对流换热关系:
在石英管壳型和外部环境之间的对流换热关系式为:
其中λ1和λ2分别为内层金属玻璃和外层石英管壳型的导热率,h1和h2分别为内层金属玻璃与外层石英管壳型之间、外层石英管壳型与外部环境之间的界面换热系数,Tsuf为外部环境温度,h1,h2,α1,α2均为常数:T1(r,t),T2(r+dx,t),T3(r+δ)均为未知,其中dx为T1和T2之间的一个厚度,他们之间满足牛顿对流换热关系;
(4)、确定初始条件:
T(x,t=0)=T0 0<x<r+δ (7)
其中T0为液体金属的初始温度;
对于所述形状的金属玻璃,可以看成半无限大固体;该液体金属的初始温
其中高斯误差函数的定义为
求解得:
并且利用布尼兹法则对上式求微分可得:
求得:
即:
最后得:
由方程(15)便可以确定出在任一时刻,任意位置处的温度;
(5)、预测大块金属玻璃形成厚度
同时根据上述不同时刻计算的温度即可求出在某一时刻的温度变化率,即只要t1和t2时间间距足够小在0.1~0.01s之间,便可以认为是瞬时冷却速度;通过比较瞬时冷却速度和温度便可以确定该点是否达到玻璃转变条件,进而判断出金属玻璃形成厚度,完成大块金属玻璃形成厚度的预测。
所述方法用于长宽比大于10的液体金属在激冷条件下的玻璃形成厚度预测,所述液体金属的侧面和底部为激冷区,顶部为与空气对流换热区域。
所述方法用于Zr-Ti-Cu-Ni-Be系合金的大块金属玻璃形成厚度预测,而且可以应用到其它系合金的大块金属玻璃形成厚度预测。
3、优点及效果:
本发明提出的一种大块金属玻璃形成厚度预测方法,该方法具有如下优点:
由于采用解析方法求解,因此该方法不存在网格剖分或网格划分所导致的计算误差,即在任意时刻,任意位置的温度内均有解。因此本方法不仅适用于Zr-Ti-Cu-Ni-Be系合金,而且可以应用到其它系合金的大块金属玻璃形成厚度预测。预测成本低、预测方法简单、预测效率高,而且适用性广。
四、附图说明:
图1为本发明金属玻璃凝固形状的示意图;
图2为本发明B-B断面示意图。
附图标记说明:
1为液体金属;2为壳型;r为金属玻璃厚度的一半;δ为壳型的厚度;Tsuf=const表示环境温度为一定值。
五、具体实施方式:
本发明可以预测不同的冷却条件下Zr-Ti-Cu-Ni-Be系合金的玻璃形成厚度。这种方法对预测其它合金系的金属玻璃形成厚度也具有适用性。
下面结合附图对本发明做进一步的说明:
一种大块金属玻璃形成厚度预测方法,其特征在于:如图1和2中所示,该方法步骤如下:
(1)、确定金属玻璃内部的传热与温度场:
在金属玻璃内部,金属玻璃的温度T与时间t之间存在着一定的关系,其关系如公式(1)所示
其中
其中T为温度(℃),t为时间(s),k1为合金的导热率(W/m2K),ρ1为合金的材料的密度(kg/m3),Cp1为合金比热容(J/kgK),α1为合金的热扩散系数(m2/s),r为金属玻璃厚度的一半;
(2)、确定壳型内部的传热与温度场:
在壳型内部,金属玻璃的温度T与时间t之间存在的关系如公式(3)所示
并且有
其中T为温度(℃),t为时间(s),k2为石英管的导热率(W/m2K),ρ2为石英管材料的密度(kg/m3),Cp2为石英管的比热容(J/kgK),α2为石英管的热扩散系数(m2/s),δ为壳型的厚度;
(3)、确定金属玻璃与壳型之间的热量传输:
在金属玻璃和壳型之间存在下面的对流换热关系:
在壳型和外部环境之间的对流换热关系式为:
其中λ1和λ2分别为内层金属和外层石英管的导热率,h1和h2分别为金属玻璃与壳型之间、壳型与外部环境之间的界面换热系数,Tsuf为外部环境温度,h1,h2,α1,α2均为常数:T1(r,t),T2(r+dx,t),T3(r+δ)均为未知,其中dx为T1和T2之间的一个厚度,他们之间满足对牛顿流换热情况;
(4)、确定初始条件:
T(x,t=0)=T0 0<x<r+δ (7)
其中T0为液体金属的初始温度;
对于本发明形状的金属玻璃,可以看成半无限大固体;该液体金属的初始温度为T0(T=T0,t=0),在t>0后表面x=0处的温度保持为Tw(T=Tw,x=0),该温度场受热扩散方程控制,于是这个方程的解为
其中高斯误差函数的定义为
求解得:
并且利用布尼兹法则对上式求微分可得:
求得:
即:
最后得:
由方程(15)便可以确定出在任一时刻,任意位置处的温度;
(5)、预测大块金属玻璃形成厚度
同时根据上述不同时刻计算的温度即可求出在某一时刻的温度变化率,即只要t1和t2时间间距足够小在0.01~0.1s之间,便可以认为是瞬时冷却速度;通过比较瞬时冷却速度和温度便可以确定该点是否达到玻璃转变条件,进而判断出金属玻璃形成厚度,完成大块金属玻璃形成厚度的预测。
所述方法用于长宽比大于10的液体金属1在激冷条件下的玻璃形成厚度预测效果最佳,所述液体金属1的侧面和底部为激冷区,顶部为与空气对流换热区域。
所述方法可以用于Zr-Ti-Cu-Ni-Be系合金的大块金属玻璃形成厚度预测,而且还可以应用到其它系合金的大块金属玻璃形成厚度预测,适用范围广。
本发明提供的这种大块金属玻璃形成厚度预测方法,克服了以往实验测量方法的缺陷,降低了预测成本、预测方法简单、提高了预测效率,而且适用性广,不仅适用于Zr-Ti-Cu-Ni-Be系合金,而且可以应用到其它系合金的大块金属玻璃形成厚度预测,适于推广应用。
Claims (3)
1.一种大块金属玻璃形成厚度预测方法,所述大块金属玻璃的形状能够看成半无限大固体,其特征在于:该方法步骤如下:
(1)、确定金属玻璃内部的传热与温度场:
在金属玻璃内部,金属玻璃的温度T与时间t之间存在着一定的关系,其关系如公式(1)所示
其中
其中T为温度(℃),t为时间(s),k1为合金的导热率(W/m2K),ρ1为合金的密度(kg/m3),Cp1为合金的比热容(J/kgK),α1为合金的热扩散系数(m2/s),r为金属玻璃厚度的一半;
(2)、确定石英管壳型内部的传热与温度场:
在石英管壳型内部,金属玻璃的温度T与时间t之间存在的关系如公式(3)所示
并且有 (4)
其中T为温度(℃),t为时间(s),k2为石英管导热率(W/m2K),ρ2为石英管的密度(kg/m3),Cp2为石英管的比热容(J/kgK),α2为石英管的热扩散系数(m2/s),δ为壳型的厚度;
(3)、确定金属玻璃与石英管壳型之间的热量传输:
在金属玻璃和石英管壳型之间存在下面的对流换热关系:
在石英管壳型和外部环境之间的对流换热关系式为:
其中λ1和λ2分别为内层金属玻璃和外层石英管壳型的导热率,h1和h2分别为内层金属玻璃与外层石英管壳型之间、外层石英管壳型与外部环境之间的界面换热系数,Tsuf为外部环境温度,h1,h2,α1,α2均为常数:T1(r,t),T2(r+dx,t),T3(r+δ)均为未知,其中dx为T1和T2之间的一个厚度,它们之间满足牛顿对流换热关系;
(4)、确定初始条件:
T(x,t=0)=T0 0<x<r+δ (7)
其中T0为液体金属的初始温度;
对于所述形状的金属玻璃,可以看成半无限大固体;该液体金属的初始温
其中高斯误差函数的定义为
求解得:
并且利用布尼兹法则对上式求微分可得:
求得:
即:
最后得:
由方程(15)便确定出在任一时刻,任意位置处的温度;
(5)、预测大块金属玻璃形成厚度
2.根据权利要求1所述的大块金属玻璃形成厚度预测方法,其特征在于:所述方法用于长宽比大于10的液体金属(1)在激冷条件下的玻璃形成厚度预测,所述液体金属(1)的侧面和底部为激冷区,顶部为与空气对流换热区域。
3.根据权利要求1所述的大块金属玻璃形成厚度预测方法,其特征在于:所述方法用于Zr-Ti-Cu-Ni-Be系合金的大块金属玻璃形成厚度预测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010563783A CN102095395B (zh) | 2010-11-29 | 2010-11-29 | 大块金属玻璃形成厚度预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010563783A CN102095395B (zh) | 2010-11-29 | 2010-11-29 | 大块金属玻璃形成厚度预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102095395A CN102095395A (zh) | 2011-06-15 |
CN102095395B true CN102095395B (zh) | 2012-10-03 |
Family
ID=44128584
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201010563783A Expired - Fee Related CN102095395B (zh) | 2010-11-29 | 2010-11-29 | 大块金属玻璃形成厚度预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102095395B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6188079B1 (en) * | 1999-01-12 | 2001-02-13 | Owens-Brockway Glass Container Inc. | Measurement of hot container wall thickness |
CN1378240A (zh) * | 2001-03-28 | 2002-11-06 | 株式会社东芝 | 半导体器件的制造方法 |
CN2521584Y (zh) * | 2001-12-07 | 2002-11-20 | 中华人民共和国辽宁出入境检验检疫局 | 镀锡量测量仪 |
CN201000327Y (zh) * | 2006-09-06 | 2008-01-02 | 天津化工研究设计院 | 一种水结垢在线检测装置 |
CN101421583A (zh) * | 2006-01-19 | 2009-04-29 | S.A.朱塞佩克里斯蒂尼有限公司 | 测量材料层厚度的方法和装置 |
-
2010
- 2010-11-29 CN CN201010563783A patent/CN102095395B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6188079B1 (en) * | 1999-01-12 | 2001-02-13 | Owens-Brockway Glass Container Inc. | Measurement of hot container wall thickness |
CN1378240A (zh) * | 2001-03-28 | 2002-11-06 | 株式会社东芝 | 半导体器件的制造方法 |
CN2521584Y (zh) * | 2001-12-07 | 2002-11-20 | 中华人民共和国辽宁出入境检验检疫局 | 镀锡量测量仪 |
CN101421583A (zh) * | 2006-01-19 | 2009-04-29 | S.A.朱塞佩克里斯蒂尼有限公司 | 测量材料层厚度的方法和装置 |
CN201000327Y (zh) * | 2006-09-06 | 2008-01-02 | 天津化工研究设计院 | 一种水结垢在线检测装置 |
Also Published As
Publication number | Publication date |
---|---|
CN102095395A (zh) | 2011-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yang et al. | Numerical analysis on the thermal behavior of high temperature latent heat thermal energy storage system | |
CN103411996B (zh) | 固体材料导热系数测量装置及测量方法 | |
CN101279360B (zh) | 直径为Ф350~Ф400mm低合金钢连铸圆管坯的生产方法 | |
CN105352992A (zh) | 测定金属泡沫多孔介质导热热阻的方法 | |
Xie et al. | Numerical study of thin layer ring on improving the ice formation of building thermal storage system | |
CN102841983B (zh) | 一种工业窑炉热效率在线监测方法 | |
CN102095395B (zh) | 大块金属玻璃形成厚度预测方法 | |
Xing et al. | Magnetocaloric effect and critical behavior in melt‐extracted Gd60Co15Al25 microwires | |
CN202201949U (zh) | 定向凝固提纯高纯铝的熔炼炉 | |
CN202421105U (zh) | 瞬态热线法测量沥青混合料热物理参数的装置 | |
CN102851685A (zh) | 一种电解槽 | |
CN204718189U (zh) | 用于制冷设备的蒸发器组件和具有其的冷柜 | |
CN103285956A (zh) | 一种金属浴 | |
Hockaday et al. | The solar thermal treatment of manganese ore pellets using closed-loop forced convection of air | |
CN104043671B (zh) | 高精度超细薄壁铝合金盘管及其生产工艺 | |
Tao et al. | Impact of the usage of a slotted cathode carbon block on thermoelectric field in an aluminum reduction cell | |
CN201107121Y (zh) | 测量流槽中铝液温度的热电偶的固定装置 | |
CN101342581A (zh) | 一室双炉真空中频感应熔炼炉 | |
CN102156184B (zh) | 铝硅合金共晶组织片层间距预测方法 | |
CN111855739A (zh) | 加压凝固过程铸锭和铸型界面换热系数的确定方法及系统 | |
Tao et al. | Numerical simulation of current distribution in cathode carbon block of an aluminum reduction cell | |
Zhai et al. | Preparation of Multi nitrate molten salt and its properties tests | |
CN205192249U (zh) | 全纤维炉衬的轧钢加热炉 | |
CN102218515A (zh) | 连铸结晶器内钢液凝固传热过程计算方法 | |
CN114239315B (zh) | 一种用于工业硅造粒生产线热能系统的计算方法 |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121003 Termination date: 20131129 |