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

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

Info

Publication number
CN109324079B
CN109324079B CN201811236329.2A CN201811236329A CN109324079B CN 109324079 B CN109324079 B CN 109324079B CN 201811236329 A CN201811236329 A CN 201811236329A CN 109324079 B CN109324079 B CN 109324079B
Authority
CN
China
Prior art keywords
temperature
tested piece
thermal expansion
expansion coefficient
measurement
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
Application number
CN201811236329.2A
Other languages
English (en)
Other versions
CN109324079A (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

Images

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,建立被测试件内部热源的热传导反问题目标函数为:
Figure BDA0001838226930000031
式中,α(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的点的温度的上升值;
所述目标函数的约束条件为:
Figure BDA0001838226930000032
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,建立被测试件内部热源的热传导反问题目标函数为:
Figure BDA0001838226930000051
式中,α(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的点的温度的上升值;
所述目标函数的约束条件为:
Figure BDA0001838226930000061
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,建立被测试件内部热源的热传导反问题目标函数为:
Figure FDA0002618459040000011
式中,α(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的点的温度的上升值;
所述目标函数的约束条件为:
Figure FDA0002618459040000021
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 CN109324079A (zh) 2019-02-12
CN109324079B true 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)

Families Citing this family (3)

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

Family Cites Families (7)

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

Also Published As

Publication number Publication date
CN109324079A (zh) 2019-02-12

Similar Documents

Publication Publication Date Title
CN109324079B (zh) 一种基于超声的材料热膨胀系数的测量方法
WO2016101903A1 (zh) 一种导热系数测定装置
WO2012167403A1 (zh) 脉冲红外热波技术测厚方法
CN108008022B (zh) 一种随温度变化的超声波传播速度测量方法
CN101126729A (zh) 双热流计稳态法材料热导率测量方法
WO2019184076A1 (zh) 一种基于虚拟热源原理测量有限厚度材料异质含量的方法
CN108051472B (zh) 一种材料高温热物性参数的快速测量方法
Kosugi et al. Accuracy evaluation of surface temperature profiling by a laser ultrasonic method
CN110220940B (zh) 一种不规则样品导热性能测试方法
CN109813753A (zh) 双向热流法测定界面接触热阻的高精度方法
CN113092523A (zh) 薄膜材料导热性能测试装置及测试方法
CN108051475B (zh) 一种对流换热系数的快速测量方法
CN109470772B (zh) 一种基于超声的内部热源强度大小和位置的无损测量方法
CN107966472B (zh) 一种高温接触热阻的无损快速测量方法
Ihara et al. New ultrasonic thermometry and its applications to temperature profiling of heated materials
CN108918580B (zh) 一种无损稳态导热率测量方法
CN114113207B (zh) 一种光学材料热扩散系数的测量方法
CN109506807B (zh) 一种稳态条件下的高温结构内部温度及壁厚同时测量方法
CN109613054A (zh) 一种直接通电纵向导热系数测试方法
Kubicar et al. Transient methods for the measurement of thermophysical properties: The pulse transient method
CN109506806B (zh) 一种瞬态条件下高温结构内部温度及厚度的同时测量方法
Liu et al. Health monitoring of RTV silicone rubber coating on insulators based on multimode features of active infrared thermography
Lisker et al. A transient technique for measuring the thermal conductivity of non-metals
Brady et al. Thermal image analysis for the in-situ NDE of composites
Halloua et al. Nondestructive inverse approach for determining thermal and geometrical properties of internal defects in CFRP composites by lock-in thermography

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