CN111307662B - 一种玻璃粘度的测试方法 - Google Patents
一种玻璃粘度的测试方法 Download PDFInfo
- Publication number
- CN111307662B CN111307662B CN201911146009.2A CN201911146009A CN111307662B CN 111307662 B CN111307662 B CN 111307662B CN 201911146009 A CN201911146009 A CN 201911146009A CN 111307662 B CN111307662 B CN 111307662B
- Authority
- CN
- China
- Prior art keywords
- temperature
- glass
- viscosity
- thermal expansion
- change rate
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N11/00—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N11/00—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
- G01N2011/006—Determining flow properties indirectly by measuring other parameters of the system
- G01N2011/0093—Determining flow properties indirectly by measuring other parameters of the system thermal properties
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种玻璃粘度的测试方法，所述玻璃粘度η通过以下MYEGA公式计算获得：
其中，T为玻璃粘度所对应的温度；Tg是玻璃化转变温度，Tg取降温/升温速率q为10℃/分钟时得到的假想温度Tf值；mvis为经校正的粘度脆性系数，其确定过程为：(1)将n个预先获取的不同降温/升温速率q条件下的假想温度值Tf代入公式
计算得到线膨胀活化能Ea，n≥5，R是气体常数；(2)计算得到线性热膨胀脆性系数mDIL，
(3)最后校正得到mvis：mvis＝1.289(mDIL‑14.97)+14.97；所述玻璃不包括分相玻璃。本发明的一种玻璃粘度的测试方法，能测量玻璃10～1015Pa·s粘度范围，且测量精度高。
Description
技术领域
本发明属于玻璃粘度的测试技术领域，涉及一种玻璃粘度的测试方法。
背景技术
玻璃的现代定义非常广泛，凡是具有玻璃化转变过程的非晶材料都可以称之为玻璃，所以玻璃包括玻璃、玻璃纤维、部分高分子材料等。粘度是决定玻璃生产的关键性能，玻璃的粘度-温度曲线能直接指导产品生产及应用研究。
由于玻璃从液态到玻璃化转变区间通常会有约1015Pa·s的粘度变化，而目前为止还没有一种仪器能一次性直接测量玻璃10～1015Pa·s粘度范围的变化，通常需要至少三种以上设备去测量玻璃10～1015Pa·s粘度范围，其中高温旋转粘度仪通常能测量10～107Pa·s，平板粘度仪通常能测量104～1010Pa·s，三点弯曲粘度仪通常能测量1010～1015Pa·s。这三种仪器都是采用直接法，通过玻璃所受到的应力和应变率去计算玻璃的粘度。
此外，上述玻璃粘度测试的仪器价格昂贵，而且多数粘度仪在使用过程中还会出现铂金损耗，额外增加了使用成本。
因此，研究一种能一次性测量玻璃10～1015Pa·s粘度范围的测试方法具有十分重要的意义。
发明内容
本发明的目的是解决现有技术中至少需要三种以上设备去测量玻璃10～1015Pa·s粘度范围的问题，提供一种玻璃粘度的测试方法。
为达到上述目的，本发明采用的技术方案如下：
一种玻璃粘度的测试方法，先确定玻璃的粘度脆性系数mvis和玻璃化转变温度Tg，再利用MYEGA粘度模型计算得到玻璃在温度T时对应的粘度η，计算公式如下：
其中，T为玻璃粘度所对应的温度，Tg是玻璃化转变温度，mvis为经校正的粘度脆性系数；
mvis和Tg的确定过程为：
首先以n个不同的升降温速率q对玻璃样品进行热处理得到n条长度变化率-温度变化曲线，n≥5，一个q为10℃/分钟，其他q在1～30℃/分钟的范围内随机分布(“随机分布”是个数学概念)；
然后对n条长度变化率-温度变化曲线进行数学处理得到n个假想温度Tf，以q为10℃/分钟时得到的Tf作为Tg；
接着由n个q和n个Tf经公式
计算得到线膨胀活化能Ea，由Ea、气体常数R和Tg经公式
计算得到线性热膨胀脆性系数mDIL；
最后由mDIL经公式mvis＝1.289(mDIL-14.97)+14.97计算得到mvis；
式中，Tg、T和Tf的单位为℃，Eu的单位为kJ/mol，R的单位为J·K-1·mol-1。本发明的方法可直接测量玻璃10～1015Pa·s粘度范围，测量精度能达到±100.4Pa·s，且成本低，操作简单，并能额外获得玻璃的热膨胀系数。
作为优选的技术方案：
如上所述的一种玻璃粘度的测试方法，热处理前将玻璃样品加热至玻璃线性热膨胀软化点Td±10℃保温2～5分钟，以消除玻璃的热历史，热处理过程分为两速率相同的降温阶段和升温阶段(通常先降温后升温)。
如上所述的一种玻璃粘度的测试方法，降温阶段的终止温度为25±10℃，升温阶段的终止温度为Td±10℃。
如上所述的一种玻璃粘度的测试方法，数学处理的步骤如下：
(1)根据
对n条长度变化率-温度变化曲线进行处理得到n条热膨胀系数-温度变化曲线，其中，α为热膨胀系数，ΔL/L为长度变化率，ΔT为温度变化率；
(2)取温度区间[T1-ΔTi,T1]对应的长度变化率-温度变化曲线段，对其进行线性拟合得到粗实线a，取温度区间[T2-ΔTj,T2]对应的长度变化率-温度变化曲线段，对其进行线性拟合得到粗实线b，T1为热膨胀系数最大值对应的温度，单位为℃，ΔTi为5～10℃，T2为位于热膨胀系数-温度变化曲线上的热膨胀系数维持不变段且最靠近玻璃化转变区间的点对应的温度，单位为℃，ΔTj为50～200℃；
(3)对粗实线a和粗实线b进行延长至相交，交点对应的温度即为Tf。
如上所述的一种玻璃粘度的测试方法，线性拟合采用excel、origin或igor软件，本发明进行线性拟合的软件包括但不仅限于这几种，其他具有线性拟合功能的软件也可用于本发明。
有益效果：
(1)本发明的一种玻璃粘度的测试方法，能直接通过样品线性热膨胀性能测量玻璃10～1015Pa·s粘度范围，测量精度可以达到±100.4Pa·s，能与三点弯曲粘度仪达到相同的精度；
(2)本发明的一种玻璃粘度的测试方法，操作简单，成本较低；
(3)本发明的一种玻璃粘度的测试方法，还能同时获得测试样品的热膨胀系数，粘度和热膨胀系数的同时获取，对一些特殊的应用，比如电子产品等的封接玻璃具有重要意义；
(4)采用本发明的一种玻璃粘度的测试方法能直接获得玻璃所需的最关键参数，为实验研究和工业生产提供重要技术参数。
附图说明
图1为硼硅酸盐样品长度变化率(图中实线)和热膨胀系数(图中虚线)随温度变化的曲线；
图2为硼硅酸盐样品粘度随温度变化曲线，其中圆点代表利用三点弯曲粘度仪测得的粘度实验数据，虚线代表利用本发明的方法得到的粘度曲线；
图3为磷硅酸盐样品粘度随温度变化曲线，其中圆点代表利用高温旋转粘度仪测得的粘度实验数据，虚线代表利用本发明的方法得到的粘度曲线。
具体实施方式
下面结合具体实施方式，进一步阐述本发明。应理解，这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解，在阅读了本发明讲授的内容之后，本领域技术人员可以对本发明作各种改动或修改，这些等价形式同样落于本申请所附权利要求书所限定的范围。
一种玻璃粘度的测试方法，先确定玻璃的粘度脆性系数mvis和玻璃化转变温度Tg，再利用公式计算得到玻璃在温度T时对应的粘度η；具体如下：
(1)确定假想温度Tf；
(1.1)将玻璃样品加热至玻璃线性热膨胀软化点Td±10℃保温2～5分钟；
(1.2)以n个不同的降温/升温速率q对玻璃样品进行热处理得到n条长度变化率-温度变化曲线，q在1～30℃/分钟的范围内随机分布；其中，热处理过程分为两速率相同的降温阶段(终止温度为25±10℃)和升温阶段(终止温度为Td±10℃)；
(1.3)对n条长度变化率-温度变化曲线进行数学处理得到n个假想温度Tf；其中，数学处理的步骤如下：
(1.3.1)根据
对n条长度变化率-温度变化曲线进行处理得到n条热膨胀系数-温度变化曲线，其中，α为热膨胀系数，AL/L为长度变化率，ΔT为温度变化率；
(1.3.2)取温度区间[T1-ΔTi,T1]对应的长度变化率-温度变化曲线段，对其进行线性拟合(采用excel、origin或igor软件)得到粗实线a，取温度区间[T2-ΔTj,T2]对应的长度变化率-温度变化曲线段，对其进行线性拟合得到粗实线b，T1为热膨胀系数最大值对应的温度，单位为℃，ΔTi为5～10℃，T2为位于热膨胀系数-温度变化曲线上的热膨胀系数维持不变段且最靠近玻璃化转变区间的点对应的温度，单位为℃，ΔTj为50～200℃；
(1.3.3)对粗实线a和粗实线b进行延长至相交，交点对应的温度即为Tf，Tf的单位为℃；
(2)确定mvis(经校正的粘度脆性系数)；
(2.1)将n个预先获取的不同降温/升温速率q条件下的假想温度值Tf代入公式
计算得到线膨胀活化能Eu，Eu的单位为kJ/mol，R是气体常数，R的单位为J·K-1·mol-1；
(2.2)计算得到线性热膨胀脆性系数mDIL，且
式中，Tg是玻璃化转变温度，Tg取降温/升温速率q为10℃/分钟时得到的假想温度Tf值，单位为℃；
(2.3)校正得到mvis：mvis＝1.289(mDIL-14.97)+14.97；
(3)将上述得到的Tf和mvis代入下面的MYEGA公式计算得到玻璃在温度T时对应的粘度η；所述的MYEGA公式为：
其中，粘度η的单位为Pa·s，T为玻璃粘度所对应的温度，单位为℃。
将上述的测试方法分别用于测试硼硅酸盐样品(组分为70wt％SiO2、12wt％B2O3、10wt％Na2O、8wt％其他)和磷硅酸盐样品(组分为26wt％P4O10、25wt％CaO、21wt％B2O3、12wt％SiO2、16wt％其他)在温度T时对应的粘度η，其中，各参数的设定见下表1～2，硼硅酸盐样品长度变化率(图中实线)和热膨胀系数(图中虚线)随温度变化的曲线如图1所示，硼酸盐样品粘度随温度变化曲线如图2所示，磷硅酸盐样品粘度随温度变化曲线如图3所示，图2中，圆点代表利用三点弯曲粘度仪测得的粘度实验数据，图3中，圆点代表利用高温旋转粘度仪测得的粘度实验数据，图2～3中虚线代表利用本发明的方法得到的粘度曲线，将实验数据与本申请中的测试数据进行对比，可以看出，本申请中计算得到的粘度与实验数据接近，本申请的玻璃粘度的测试方法能直接通过样品线性热膨胀性能测量玻璃10～1015Pa·s粘度范围，测量精度可以达到±100.4Pa·s，能与三点弯曲粘度仪达到相同的精度。本申请的方法操作简单，成本较低；粘度和热膨胀系数的同时获取，对一些特殊的应用，比如电子产品等的封接玻璃具有重要意义。
表1硼硅酸盐样品的参数表
表2磷硅酸盐样品的参数表
Claims (4)
1.一种玻璃粘度的测试方法，其特征是：先确定玻璃的粘度脆性系数mvis和玻璃化转变温度Tg，再利用MYEGA粘度模型计算得到玻璃在温度T时对应的粘度η；
其中，mvis和Tg的确定过程为：
首先以n个不同的升降温速率q对玻璃样品进行热处理得到n条长度变化率-温度变化曲线，n≥5，一个q为10℃/分钟，其他q在1～30℃/分钟的范围内随机分布；
然后对n条长度变化率-温度变化曲线进行数学处理得到n个假想温度Tf，以q为10℃/分钟时得到的Tf作为Tg；
接着由n个q和n个Tf经公式
计算得到线膨胀活化能Ea，由Ea、气体常数R和Tg经公式
计算得到线性热膨胀脆性系数mDIL；
最后由mDIL经公式mvis＝1.289(mDIL-14.97)+14.97计算得到mvis；
式中，Tg、T和Tf的单位为℃，Ea的单位为kJ/mol，R的单位为J·K-1·mol-1；
数学处理的步骤如下：
(1)根据
对n条长度变化率-温度变化曲线进行处理得到n条热膨胀系数-温度变化曲线，其中，α为热膨胀系数，ΔL/L为长度变化率，ΔT为温度变化率；
(2)取温度区间[T1-ΔTi,T1]对应的长度变化率-温度变化曲线段，对其进行线性拟合得到粗实线a，取温度区间[T2-ΔTj,T2]对应的长度变化率-温度变化曲线段，对其进行线性拟合得到粗实线b，T1为热膨胀系数最大值对应的温度，单位为℃，ΔTi为5～10℃，T2为位于热膨胀系数-温度变化曲线上的热膨胀系数维持不变段且最靠近玻璃化转变区间的点对应的温度，单位为℃，ΔTj为50～200℃；
(3)对粗实线a和粗实线b进行延长至相交，交点对应的温度即为Tf。
2.根据权利要求1所述的一种玻璃粘度的测试方法，其特征在于，热处理前将玻璃样品加热至玻璃线性热膨胀软化点Td±10℃保温2～5分钟，热处理过程分为两速率相同的降温阶段和升温阶段。
3.根据权利要求2所述的一种玻璃粘度的测试方法，其特征在于，降温阶段的终止温度为25±10℃，升温阶段的终止温度为Td±10℃。
4.根据权利要求1所述的一种玻璃粘度的测试方法，其特征在于，线性拟合采用excel、origin或igor软件。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911146009.2A CN111307662B (zh) | 2019-11-21 | 2019-11-21 | 一种玻璃粘度的测试方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911146009.2A CN111307662B (zh) | 2019-11-21 | 2019-11-21 | 一种玻璃粘度的测试方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111307662A CN111307662A (zh) | 2020-06-19 |
| CN111307662B true CN111307662B (zh) | 2021-01-05 |
Family
ID=71144824
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911146009.2A Active CN111307662B (zh) | 2019-11-21 | 2019-11-21 | 一种玻璃粘度的测试方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111307662B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112763532A (zh) * | 2020-12-28 | 2021-05-07 | 天津大学 | 基于量热学测试的玻璃化转变脆性因子分析方法 |
| CN113340504B (zh) * | 2021-07-13 | 2022-03-01 | 中国工程物理研究院激光聚变研究中心 | 一种从熔石英假想温度分布获取残余应力分布的方法 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0416741A (ja) * | 1990-05-11 | 1992-01-21 | Mitsubishi Materials Corp | ペネトレーション試験用貫入ピン |
| DE19932601B4 (de) * | 1999-07-13 | 2004-09-02 | Schott Glas | Verfahren zur Bestimmung der Temperatur und/oder der Viskosität einer Glasschmelze |
| CN103558201B (zh) * | 2013-11-12 | 2015-07-22 | 中山大学 | 一种灵敏探测聚合物玻璃化转变的方法 |
| CN106950154A (zh) * | 2017-04-25 | 2017-07-14 | 武汉理工大学 | 一种玻璃形成的临界冷却速率评测方法 |
| CN207215651U (zh) * | 2017-09-27 | 2018-04-10 | 成都光明光电股份有限公司 | 玻璃粘度测量装置 |
| CN109101782B (zh) * | 2018-09-07 | 2021-12-03 | 山东大学 | 一种预测Cu-Zr非晶合金强脆转变参数f的方法 |
| CN209117521U (zh) * | 2018-11-23 | 2019-07-16 | 大连海事大学 | 一种快速全自动运动粘度测定装置 |
-
2019
- 2019-11-21 CN CN201911146009.2A patent/CN111307662B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN111307662A (zh) | 2020-06-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111307662B (zh) | 一种玻璃粘度的测试方法 | |
| Fulcher | Analysis of recent measurements of the viscosity of glasses | |
| Strawbridge et al. | The factors affecting the thickness of sol-gel derived silica coatings prepared by dipping | |
| Pascual et al. | Determination of the viscosity–temperature curve for glasses on the basis of fixed viscosity points determined by hot stage microscopy | |
| Simmons et al. | Non-Newtonian viscous flow in soda-lime-silica glass at forming and annealing temperatures | |
| Ray et al. | Simple method for theoretical estimation of viscosity of oxide melts using optical basicity | |
| CN107056249B (zh) | 高性能白瓷坯体及其制备方法 | |
| Mitang et al. | Effect of rare earths on viscosity and thermal expansion of soda-lime-silicate glass | |
| Mazurin et al. | Temperature dependences of viscosity of glass-forming substances at constant fictive temperatures | |
| Liu et al. | Influence of Y2O3 substitution for B2O3 on the structure and properties of alkali-free B2O3-Al2O3-SiO2 glasses containing alkaline-earth metal oxides | |
| Fulcher | Analysis of recent measurements of the viscosity of glasses | |
| O'Connell et al. | The non-linear viscoelastic response of polycarbonate in torsion: An investigation of time-temperature and time-strain superposition | |
| Nemilov | Physical ageing of silicate glasses at room temperature: general regularities as a basis for the theory and the possibility of a priori calculation of the ageing rate | |
| CN105776868A (zh) | 一种耐高温高绝缘封接玻璃及其制备方法 | |
| CN102494968A (zh) | 自动化测试凝胶化状态的测试方法 | |
| Robertson et al. | Refractive index: A probe for monitoring volume relaxation during physical aging of glassy polymers | |
| CN110749617B (zh) | 玻璃实用软化点温度的测量方法 | |
| Zheng et al. | Determining the liquidus viscosity of glass‐forming liquids through differential scanning calorimetry | |
| Ding et al. | Dilatometric fragility and prediction of the viscosity curve of glass‐forming liquids | |
| CN107365064B (zh) | 离子交换熔盐组合物和钢化盖板玻璃及其钢化方法 | |
| Lancelotti et al. | Structural relaxation dynamics of a silicate glass probed by refractive index and ionic conductivity | |
| He et al. | Viscosity and structure of lithium sodium borosilicate glasses | |
| Ding et al. | Pressure dependence of density and structural relaxation of glass near the glass transition region | |
| XIAO et al. | Effect of Al2O3/SiO2 ratio on the viscosity and workability of high-alumina soda-lime-silicate glasses | |
| CN107200479A (zh) | 一种耐碱耐热玻璃纤维及其制备方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |