CN109324079A - 一种基于超声的材料热膨胀系数的测量方法 - Google Patents

一种基于超声的材料热膨胀系数的测量方法 Download PDF

Info

Publication number
CN109324079A
CN109324079A CN201811236329.2A CN201811236329A CN109324079A CN 109324079 A CN109324079 A CN 109324079A CN 201811236329 A CN201811236329 A CN 201811236329A CN 109324079 A CN109324079 A CN 109324079A
Authority
CN
China
Prior art keywords
temperature
test specimen
thermal expansion
expansion coefficient
ultrasound
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
CN201811236329.2A
Other languages
English (en)
Other versions
CN109324079B (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.)
China Special Equipment Inspection and Research Institute
Original Assignee
China Special Equipment Inspection and Research Institute
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 China Special Equipment Inspection and Research Institute filed Critical China Special Equipment Inspection and Research Institute
Priority to CN201811236329.2A priority Critical patent/CN109324079B/zh
Publication of CN109324079A publication Critical patent/CN109324079A/zh
Application granted granted Critical
Publication of CN109324079B publication Critical patent/CN109324079B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N25/00Investigating or analyzing materials by the use of thermal means
    • G01N25/16Investigating or analyzing materials by the use of thermal means by investigating thermal coefficient of expansion

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)

Abstract

本发明提供了一种基于超声的材料热膨胀系数的测量方法,该方法基于声传播特性、介质温度、材料热膨胀系数三者之间的关联关系,采用超声回波法,根据超声传播测量时间和参数识别,实现了高温条件下材料热膨胀系数的无损、快速测量。本发明仅需进行升温测量一次,例如将被测试件加热面进行升温到预定温度值如500℃,即可获得室温至500℃不同温度下的材料热膨胀系数,具有测量速度快、周期短、成本低、操作简便等特点。

Description

一种基于超声的材料热膨胀系数的测量方法
技术领域
本发明属于超声检测技术领域,具体涉及一种基于超声的材料热膨胀系数的测量方法。
背景技术
热膨胀系数是材料的重要热学性能参数之一,也是结构热固耦合特性分析中的关键参数。目前,热膨胀系数的测量方法主要有相对比较法和绝对比较法两种,按测量原理又分为电测量法、光测量法和位移测量法。总体来说,国内外关于材料热膨胀系统的测量研究较为成熟,但多数实验测量存在测量设备价格昂贵、测量系统构成复杂、或者抗干扰或受人为因素制约等局限性。
本发明根据超声传播特性——介质温度——材料热膨胀系数三者之间的关联关系,采用热声理论,基于超声回波技术,可快速、简便、准确地测量材料随温度变化的热膨胀系数。相比现有的热膨胀系数测量方法,本发明一方面,具有测量装置简单,操作相对简便,适用于大尺寸试样测试以及试样内部处于非均匀温度场分布等条件;另一方面,由于该方法可以非接触式测量,提高了传感器耐高温性能,可有效提升测温范围。同时,本发明采用瞬态升温测量法,仅需进行升温测量一次,例如被测试件加热面进行升温到预定温度值如500℃,即可获得室温至500℃不同温度下的材料热膨胀系数,可极大地缩短测量周期。
发明内容
针对现有技术存在的以上问题,本发明提出一种基于超声的材料热膨胀系数的测量方法,适用于材料随温度变化的等效热膨胀系数的无损非接触式测量。其采用超声回波法,由瞬态传热条件下超声传播时间,基于热传导反问题的参数识别方法,可无损、快速地测量材料随温度变化的热膨胀系数常数。
为实现上述目的,本发明的技术方案具体如下:
包括以下步骤:
步骤一:加工一与被测物体材质相同的被测试件,被测试件为长方体或圆柱体,选定某一被测方向,设定其为x方向,L为被测试件沿被测方向的长度,L大于一个超声波波长;
步骤二:获取被测试件内超声波传播速度V与不同温度条件下温度T的相关关系,拟合出V与T的相关关系式V(T),V(T)为线性函数或二次函数,其中,T为被测试件整体温度达到一致的条件下的温度值;
步骤三:记录被测试件的初始温度T0,沿被测方向x方向从被测试件x=0处的一端进行加热,直至被测试件所述x=0处的加热端的温度高于测量的预定温度时,结束加热;
在所述加热初始至结束的过程中,进行以下步骤四至步骤五,随后,基于步骤四、步骤五的结果,进行步骤六;
步骤四:沿所选定的被测方向x方向,由超声波脉冲回波法,获得被测试件ti时刻的超声波传播测量时间ti,exp,即在步骤三所述加热初始至结束的过程中进行多次测量,下标i表示测量时间序数,i=1,2,……,n,共测量n次;当被测试件加热端的温度高于测量的预定温度时,进行第n次测量,获得tn,exp
步骤五:基于步骤二的相关关系式V(T),和步骤四的超声波传播时间ti,exp,建立被测试件内部热源的热传导反问题目标函数为:
式中,α(T)为待测量的温度T条件下的材料等效热膨胀系数;T(ti)为被测试件ti时刻加热端的温度值;ti,cal为超声波传播计算时间,是将ti时刻对应的被测试件加热端的温度值T(ti)代入步骤二所得的相关关系式V(T)得到的V,再根据L与V计算得到的;ti,exp为步骤四中实际测量得到的ti时刻的超声波传播测量时间;ΔT=T(x,ti)-T0为被测试件内沿被测方向x方向上坐标为x的点的温度的上升值;
所述目标函数的约束条件为:
T(x,t),t>0,x∈[0,L] (2b)
式中,k(T)为被测试件材料的导热系数,为温度T的函数;Cp和ρ分别为被测试件材料的比热容和密度;T(x,t)表示所述加热初始至结束的过程中所述试件内部形成的温度梯度场中,沿被测方向x方向上一点的温度随时间t和坐标x的分布;T(x,t)|x=0/x=L为边界条件,通过红外测温方法或热电偶测温方法获得;
步骤六:采用灵敏度法或梯度算法求解所述目标函数,获得不同温度条件下被测试件的等效热膨胀系数α(T),即得到了被测试件材料的热膨胀系数。
进一步地,所述步骤三的具体步骤为,当需要获得一预定温度下的材料热膨胀系数时,对被测试件加热的边界温度应高于所述预定温度,即步骤三所述加热过程结束时,加热端的温度高于预定温度。
进一步地,所述超声为电磁超声或激光超声。
本发明具有如下优点和积极效果:
1、所述方法为瞬态升温测量法,仅需进行升温测量一次,例如将被测试件加热面进行升温到预定温度值如500℃,即可获得室温至500℃不同温度下的等效热膨胀系数,具有测量速度快、周期短、成本低、操作简便等优势,适用于大尺寸试样测试以及试样内部处于非均匀温度场分布等条件。
2、基于电磁超声或激光超声进行非接触测量时,材料的等效热膨胀系数测量几乎不受传感器耐温性能的影响,具有测量范围大的优点;同时,基于超声测温原理,该方法对高温条件下的材料等效热膨胀系数具有更佳的测量精度。
附图说明
图1为材料热膨胀系数测量方法的流程。
图2为材料热膨胀系数随温度变化的测量结果。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的说明。
本发明所述基于超声的材料热膨胀系数的测量方法的具体流程如图1所示。
步骤一:加工一与被测物体材质相同的被测试件,被测试件为长方体或圆柱体,选定某一被测方向,设定其为x方向,L为被测试件沿被测方向的长度,L大于一个超声波波长;
步骤二:获取被测试件内超声波传播速度V与不同温度条件下温度T的相关关系,拟合出V与T的相关关系式V(T),V(T)为线性函数或二次函数,其中,T为被测试件整体温度达到一致的条件下的温度值;
步骤三:记录被测试件的初始温度T0,沿被测方向x方向从被测试件x=0处的一端进行加热,直至被测试件所述x=0处的加热端的温度高于测量的预定温度时,结束加热;
在所述加热初始至结束的过程中,进行以下步骤四至步骤五,随后,基于步骤四、步骤五的结果,进行步骤六;
当需要获得预定温度500℃下的材料热膨胀系数时,被测试件的底部即加热端由平板加热器加热至500℃以上,其余表面皆为绝热面,
步骤四:沿所选定的被测方向x方向,由超声波脉冲回波法,获得被测试件ti时刻的超声波传播测量时间ti,exp,即在步骤三所述加热初始至结束的过程中进行多次测量,下标i表示测量时间序数,i=1,2,……,n,共测量n次;当被测试件加热端的温度高于测量的预定温度时,进行第n次测量,获得tn,exp
沿所选定的被测方向,将超声探头置于被测试件与加热端即底部相对的上端面,采用垂直入射方式激发脉冲超声波,获得n个ti,exp,i=1,2,……,n;
步骤五:基于步骤二的相关关系式V(T),和步骤四的超声波传播时间ti,exp,建立被测试件内部热源的热传导反问题目标函数为:
式中,α(T)为待测量的温度T条件下的材料等效热膨胀系数;T(ti)为被测试件ti时刻加热端的温度值;ti,cal为超声波传播计算时间,是将ti时刻对应的被测试件加热端的温度值T(ti)代入步骤二所得的相关关系式V(T)得到的V,再根据L与V计算得到的;ti,exp为步骤四中实际测量得到的ti时刻的超声波传播测量时间;ΔT=T(x,ti)-T0为被测试件内沿被测方向x方向上坐标为x的点的温度的上升值;
所述目标函数的约束条件为:
T(x,t),t>0,x∈[0,L] (2b)
式中,k(T)为被测试件材料的导热系数,为温度T的函数;Cp和ρ分别为被测试件材料的比热容和密度;T(x,t)表示所述加热初始至结束的过程中所述试件内部形成的温度梯度场中,沿被测方向x方向上一点的温度随时间t和坐标x的分布;T(x,t)|x=0/x=L为边界条件,通过红外测温方法或热电偶测温方法获得;
步骤六:采用灵敏度法或梯度算法求解所述目标函数,获得不同温度条件下被测试件的等效热膨胀系数α(T),即得到了被测试件材料的热膨胀系数。
即经过步骤四和步骤五以基于测量底面回波传播时间的变化,通过步骤六求解所述目标函数,获得材料的等效热膨胀系数。
图2给出了根据本发明所述方法获得的材料随温度变化的热膨胀系数测量结果。图中实线部分是由实验数据预先拟合获得的已知被测材料室温至470℃范围内热膨胀系数真值为:
α(T)=4×10-14×T3-2×10-11×T2-1×10-8×T+1×10-5
其中T为温度。带方块标识的实线为本发明所述方法测量所获得的以6分段函数表征的测量结果,除初始状态外,其余温度条件下的等效热膨胀系数与真值基本吻合,平均误差小于1.2%。
上述实施方式为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合和/或简化,均应视为等效的置换方式,均应包含在本发明的保护范围之内。

Claims (3)

1.一种基于超声的材料热膨胀系数的测量方法,其特征在于,包括以下步骤:
步骤一:加工一与被测物体材质相同的被测试件,被测试件为长方体或圆柱体,选定某一被测方向,设定其为x方向,L为被测试件沿被测方向的长度,L大于一个超声波波长;
步骤二:获取被测试件内超声波传播速度V与不同温度条件下温度T的相关关系,拟合出V与T的相关关系式V(T),V(T)为线性函数或二次函数,其中,T为被测试件整体温度达到一致的条件下的温度值;
步骤三:记录被测试件的初始温度T0,沿被测方向x方向从被测试件x=0处的一端进行加热,直至被测试件所述x=0处的加热端的温度高于测量的预定温度时,结束加热;
在所述加热初始至结束的过程中,进行以下步骤四至步骤五,随后,基于步骤四、步骤五的结果,进行步骤六;
步骤四:沿所选定的被测方向x方向,由超声波脉冲回波法,获得被测试件ti时刻的超声波传播测量时间ti,exp,即在步骤三所述加热初始至结束的过程中进行多次测量,下标i表示测量时间序数,i=1,2,……,n,共测量n次;当被测试件加热端的温度高于预定温度时,进行第n次测量,获得tn,exp
步骤五:基于步骤二的相关关系式V(T),和步骤四的超声波传播时间ti,exp,建立被测试件内部热源的热传导反问题目标函数为:
式中,α(T)为待测量的温度T条件下的材料等效热膨胀系数;T(ti)为被测试件ti时刻加热端的温度值;ti,cal为超声波传播计算时间,是将ti时刻对应的被测试件加热端的温度值T(ti)代入步骤二所得的相关关系式V(T)得到的V,再根据L与V计算得到的;ti,exp为步骤四中实际测量得到的ti时刻的超声波传播测量时间;ΔT=T(x,ti)-T0为被测试件内沿被测方向x方向上坐标为x的点的温度的上升值;
所述目标函数的约束条件为:
T(x,t),t>0,x∈[0,L] (2b)
式中,k(T)为被测试件材料的导热系数,为温度T的函数;Cp和ρ分别为被测试件材料的比热容和密度;T(x,t)表示所述加热初始至结束的过程中所述试件内部形成的温度梯度场中,沿被测方向x方向上一点的温度随时间t和坐标x的分布;T(x,t)|x=0/x=L为边界条件,通过红外测温方法或热电偶测温方法获得;
步骤六:采用灵敏度法或梯度算法求解所述目标函数,获得不同温度条件下被测试件的等效热膨胀系数α(T),即得到了被测试件材料的热膨胀系数。
2.根据权利要求1所述的一种基于超声的材料热膨胀系数的测量方法,其特征在于,所述步骤三的具体步骤为,当需要获得一预定温度下的材料热膨胀系数时,对被测试件加热的边界温度高于所述预定温度,即步骤三所述加热过程结束时,加热端的温度高于预定温度。
3.根据权利要求1所述的一种基于超声的材料热膨胀系数的测量方法,其特征在于,所述超声为电磁超声或激光超声。
CN201811236329.2A 2018-10-23 2018-10-23 一种基于超声的材料热膨胀系数的测量方法 Active CN109324079B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811236329.2A CN109324079B (zh) 2018-10-23 2018-10-23 一种基于超声的材料热膨胀系数的测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811236329.2A CN109324079B (zh) 2018-10-23 2018-10-23 一种基于超声的材料热膨胀系数的测量方法

Publications (2)

Publication Number Publication Date
CN109324079A true CN109324079A (zh) 2019-02-12
CN109324079B CN109324079B (zh) 2020-11-06

Family

ID=65263123

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811236329.2A Active CN109324079B (zh) 2018-10-23 2018-10-23 一种基于超声的材料热膨胀系数的测量方法

Country Status (1)

Country Link
CN (1) CN109324079B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111366640A (zh) * 2020-04-30 2020-07-03 中国科学院光电技术研究所 一种超低膨胀玻璃热膨胀系数测量方法
CN115950916A (zh) * 2023-03-14 2023-04-11 中国空气动力研究与发展中心计算空气动力研究所 一种物体表面热流密度检测方法、装置以及设备
CN116660387A (zh) * 2023-05-16 2023-08-29 西南科大四川天府新区创新研究院 一种借助声速-温度模型快速计算热膨胀系数的方法
CN117452191A (zh) * 2023-12-22 2024-01-26 威科电子模块(深圳)有限公司 基于深度学习实现厚膜电路板的耐高温测试方法及系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110043787A1 (en) * 2009-08-20 2011-02-24 Carlos Duran Photoelastic method for absolute determination of zero cte crossover in low expansion silica-titania glass samples
WO2013116412A1 (en) * 2012-02-02 2013-08-08 Schlumberger Canada Limited Determination of thermodynamic properties of a fluid based on density and sound speed
CN103713009A (zh) * 2013-12-24 2014-04-09 中国科学院深圳先进技术研究院 热膨胀系数的测定方法
CN104596667A (zh) * 2015-01-05 2015-05-06 中国空气动力研究与发展中心计算空气动力研究所 超声波探测物体内部瞬态非均匀温度场的灵敏度方法
JP6264591B1 (ja) * 2017-03-28 2018-01-24 パルステック工業株式会社 熱膨張係数測定方法及びx線回折測定装置
CN108051472A (zh) * 2017-12-05 2018-05-18 中国空气动力研究与发展中心计算空气动力研究所 一种材料高温热物性参数的快速测量方法
KR20180077812A (ko) * 2016-12-29 2018-07-09 중앙대학교 산학협력단 탄성 멤브레인을 활용한 경화 중인 시멘트 복합체의 열팽창계수 측정방법

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110043787A1 (en) * 2009-08-20 2011-02-24 Carlos Duran Photoelastic method for absolute determination of zero cte crossover in low expansion silica-titania glass samples
WO2013116412A1 (en) * 2012-02-02 2013-08-08 Schlumberger Canada Limited Determination of thermodynamic properties of a fluid based on density and sound speed
CN103713009A (zh) * 2013-12-24 2014-04-09 中国科学院深圳先进技术研究院 热膨胀系数的测定方法
CN104596667A (zh) * 2015-01-05 2015-05-06 中国空气动力研究与发展中心计算空气动力研究所 超声波探测物体内部瞬态非均匀温度场的灵敏度方法
KR20180077812A (ko) * 2016-12-29 2018-07-09 중앙대학교 산학협력단 탄성 멤브레인을 활용한 경화 중인 시멘트 복합체의 열팽창계수 측정방법
JP6264591B1 (ja) * 2017-03-28 2018-01-24 パルステック工業株式会社 熱膨張係数測定方法及びx線回折測定装置
CN108051472A (zh) * 2017-12-05 2018-05-18 中国空气动力研究与发展中心计算空气动力研究所 一种材料高温热物性参数的快速测量方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
BEN-TING CHEN: "ultrasound thermal mapping based on a hybrid method combining physical and statistical models", 《ULTRASOUND IN MED.&BIOL.》 *
施超: "固体结构内部温度的超声测量", 《中国计量大学学报》 *
骆湘红: "基于想干探测的激光超声检测系统设计与实现", 《激光与红外》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111366640A (zh) * 2020-04-30 2020-07-03 中国科学院光电技术研究所 一种超低膨胀玻璃热膨胀系数测量方法
CN115950916A (zh) * 2023-03-14 2023-04-11 中国空气动力研究与发展中心计算空气动力研究所 一种物体表面热流密度检测方法、装置以及设备
CN115950916B (zh) * 2023-03-14 2023-05-26 中国空气动力研究与发展中心计算空气动力研究所 一种物体表面热流密度检测方法、装置以及设备
CN116660387A (zh) * 2023-05-16 2023-08-29 西南科大四川天府新区创新研究院 一种借助声速-温度模型快速计算热膨胀系数的方法
CN116660387B (zh) * 2023-05-16 2024-08-06 西南科大四川天府新区创新研究院 一种借助声速-温度模型快速计算热膨胀系数的方法
CN117452191A (zh) * 2023-12-22 2024-01-26 威科电子模块(深圳)有限公司 基于深度学习实现厚膜电路板的耐高温测试方法及系统
CN117452191B (zh) * 2023-12-22 2024-04-12 威科电子模块(深圳)有限公司 基于深度学习实现厚膜电路板的耐高温测试方法及系统

Also Published As

Publication number Publication date
CN109324079B (zh) 2020-11-06

Similar Documents

Publication Publication Date Title
CN109324079A (zh) 一种基于超声的材料热膨胀系数的测量方法
CN101126729B (zh) 双热流计稳态法材料热导率测量方法
CN104792435B (zh) 基于瞬态热边界反演的结构内部非均匀温度场的重建方法
CN109813753A (zh) 双向热流法测定界面接触热阻的高精度方法
CN108051472B (zh) 一种材料高温热物性参数的快速测量方法
CN108008022B (zh) 一种随温度变化的超声波传播速度测量方法
CN105628790B (zh) 一种基于材料物性参数变化的结构内部温度场测量方法
CN108051475B (zh) 一种对流换热系数的快速测量方法
CN107748205A (zh) 一种随温度变化的弹性常数测量方法
CN107966472B (zh) 一种高温接触热阻的无损快速测量方法
CN109470772A (zh) 一种基于超声的内部热源强度大小和位置的无损测量方法
CN106525564A (zh) 热冲击‑力耦合加载与测试系统
Guo et al. Properties of Ceramic Substrate Materials for High‐Temperature Pressure Sensors for Operation above 1000° C
Zhang et al. An improved hot probe for measuring thermal conductivity of liquids
Hashimoto et al. Temperature wave analysis
CN207964949U (zh) 电介质材料微波性能测试装置
CN107356627B (zh) 一种基于虚拟热源采用四参数匹配测定材料异质含量的方法
RU2460063C1 (ru) Способ определения теплопроводности и температуропроводности твердого тела при нестационарном тепловом режиме
RU2551663C2 (ru) Способ определения теплопроводности твердого тела цилиндрической формы при стационарном тепловом режиме
CN109506807A (zh) 一种稳态条件下的高温结构内部温度及壁厚同时测量方法
Lisker et al. A transient technique for measuring the thermal conductivity of non-metals
Szczotka Non-invasive methods in diagnosis of wall dampness degree in sacral buildings
Sapozhnikov et al. Bismuth-based gradient heat-flux sensors in thermal experiment
CN106124556B (zh) 基于红外热像仪的热参数反演测量装置
Li et al. Temperature Field in Laser Line Scanning Thermography: Analytical Calculation and Experiment.

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