CN109470772B - 一种基于超声的内部热源强度大小和位置的无损测量方法 - Google Patents

一种基于超声的内部热源强度大小和位置的无损测量方法 Download PDF

Info

Publication number
CN109470772B
CN109470772B CN201811236326.9A CN201811236326A CN109470772B CN 109470772 B CN109470772 B CN 109470772B CN 201811236326 A CN201811236326 A CN 201811236326A CN 109470772 B CN109470772 B CN 109470772B
Authority
CN
China
Prior art keywords
temperature
heat source
tested
intensity
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
CN201811236326.9A
Other languages
English (en)
Other versions
CN109470772A (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 CN201811236326.9A priority Critical patent/CN109470772B/zh
Publication of CN109470772A publication Critical patent/CN109470772A/zh
Application granted granted Critical
Publication of CN109470772B publication Critical patent/CN109470772B/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
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04Analysing solids
    • G01N29/07Analysing solids by measuring propagation velocity or propagation time of acoustic waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/01Indexing codes associated with the measuring variable
    • G01N2291/011Velocity or travel time

Landscapes

  • Physics & Mathematics (AREA)
  • Acoustics & Sound (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 Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明提供了一种基于超声的内部热源强度大小和位置的无损测量方法,其特征在于:同时考虑超声在结构内部传播过程中的时延特性和边界上的温度变化,结合热传导反问题求解,用于无损、快速地测量结构内部随时间变化的热源的强度和位置。本发明不需要获取结构内部点的温度信息,简化了测量装置,操作简便;相比内部热源的红外检测技术,能够更有效地探测结构内部热源的强度和位置的状况,具有测量速度快、测量范围大等优点。

Description

一种基于超声的内部热源强度大小和位置的无损测量方法
技术领域
本发明属于超声检测技术领域,具体涉及一种基于超声的内部热源强度大小和位置的无损测量方法。
背景技术
结构内部热源的测量在无损检测、微生物发酵、弹药和食物存储等领域都有重要的应用。例如,医学上的肿瘤诊断和治疗等问题本质上可归结为对物体内部热源的位置识别的问题。内部热源的识别问题,特别是强度和位置同时识别,属于高度的不适定(ill-posed)问题,即解的惟一性和解的稳定性都不满足。尽管内部热源识别问题的研究已有大量报道,但研究多集中在通过数值算法求解热源逆问题,求解的稳定性和准确度不易达到预期,并且多数研究方法也并没有进行实验验证,同时一些方法的局限性在于仍需要探测少数内部点的温度信息。
本发明以超声在结构内部传播过程中的时延特性和边界温度变化为共同出发点,根据超声在材料介质中传播特性随温度变化的特性反映了内部热源引发的结构内部温度变化,同时边界温度进一步提供了优化的边界约束条件,从而可快速、简便、准确地测量内部热源的强度和位置。相比现有方法,本发明不需要获取内部点的温度信息,这大大简化了测量装置,操作相对简便;与内部热源的红外探测技术相比,本发明提供了一种可进行内部热源探测的新途径,将超声传播的时间变为识别内部热源强度和位置的参数,能够更有效地探测结构内部热源参数的分布状况。此外,由于声学设备可以实现非接触式测量,其提高了传感器耐高温性能和温度效应的抗干扰能力,有效地提升了测温范围和测温精度。
发明内容
针对现有技术存在的以上问题,本发明提出一种基于超声的内部热源强度大小和位置的无损测量方法,适用于单个或多个内部热源的强度和位置的无损测量。该方法根据超声在金属、复合材料等介质中随温度变化而出现的传播速度变化特性,采用超声回波法,由瞬态传热条件下超声传播时间,基于热传导反问题的参数识别方法,可无损、快速地测量结构内部随时间变化的热源的强度和位置。
为实现上述目的,本发明的技术方案具体如下:
包括以下步骤:
步骤一:获取被测构件内超声传播速度V与不同温度条件下温度T的相关关系,拟合出V与T的关系式V(T),其中,T为被测试件整体温度达到一致的条件下的温度值;
步骤二:对被测构件进行加热,选定某一被测方向,设定其为x方向,L为被测试件沿被测方向的长度,L大于一个超声波波长;沿被测方向x方向从被测构件x=0处的一端进行加热,直至被测试件所述x=0处的加热端的温度高于测量的预定温度时,结束加热;
在所述加热初始至结束的过程中,进行以下步骤三至步骤五,随后,基于步骤三至步骤五的结果,进行步骤六;
步骤三:由超声波脉冲回波法,获得被测试件ti时刻的超声波传播测量时间ti,exp,即在步骤二所述加热初始至结束的过程中进行多次测量,下标i表示测量时间序数,i=1,2,……,n,共测量n次;
步骤四:基于步骤一的关系式V(T),和步骤三的超声传播测量时间ti,exp,建立被测构件内部热源的热传导反问题目标函数为:
Figure BDA0001838226640000031
式(1)中,x表示点热源的位置,即沿被测方向x方向被测点热源的坐标值;Q(x,t)为点热源的强度,其随时间t而变化;ti,cal为超声波传播计算时间,是将ti时刻对应的被测试件加热端的温度值T(ti)代入步骤二所得的相关关系式V(T)得到的V,再根据L与V计算得到的;ti,exp为步骤三中实际测量得到的ti时刻的超声波传播测量时间;
所述热传导反问题目标函数的约束条件为:
Figure BDA0001838226640000032
T(x,t),t>0,x∈[0,L] (3)
式中,k(T)为被测构件材料的导热系数,为温度T的函数;Cp和ρ分别为被测构件材料的比热容和密度;T(x,t)表示所述加热初始至结束的过程中所述试件内部形成的温度梯度场中,沿被测方向x方向上一点的温度随时间t和坐标x的分布;T(x,t)|x=0/x=L为边界条件;
步骤五:由红外测温方法或热电偶测温方法,获得T(x,t)|x=0/x=L的边界温度;
步骤六:根据步骤五所得的温度边界条件T(x,t)|x=0/x=L,求解上述公式(1)和(2),获得x和t,从而测出构件内部热源沿被测方向的强度和位置;
进一步地,该方法能够测量被测构件内部已知或未知位置的单个或多个热源的强度。
进一步地,该方法适用于金属或复合材料的被测构件。
进一步地,该方法步骤一所述V与T的关系式V(T)是线性关系式或二次函数关系式。
本发明具有如下优点和积极效果:
1、所述方法不需要获取内部点的温度信息,大大简化了测量装置,操作相对简便;相比内部热源的红外探测技术,能够更有效地探测结构内部热源的强度和位置。
2、基于超声进行非接触测量时,内部热源强度和位置的测量几乎不受传感器耐温性能的影响,具有测量范围大的优点。
附图说明
图1为基于超声的构件内部热源强度和位置的无损测量方法的流程;
图2为随时间变化的内部热源测量。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的说明。
实施例1:
被测构件为长方体,长度为0.1m,中心位置设置单个点热源Q=1e+6W/m3,即理论上近似为一维含内部热源的传热问题。
步骤一:获取被测构件内超声传播速度V与不同温度条件下温度T的相关关系,拟合出V与T的关系式V(T),其中,T为被测试件整体温度达到一致的条件下的温度值;
步骤二:选定长度方向为被测方向,设定其为x方向,L为被测试件沿被测方向的长度,L=0.1m;沿长度方向除上、下表面外,加装隔热层,沿被测方向x方向从被测构件x=0处的一端进行加热,直至被测试件所述x=0处的加热端的温度高于测量的预定温度时,结束加热;
预定温度是为加热端加热的加热器所事先设定的停止加热工作的温度值;采用平板加热器,加热器与被测构件的加热端紧密接触;
在所述加热初始至结束的过程中,进行以下步骤三至步骤五,随后,基于步骤三至步骤五的结果,进行步骤六;
步骤三:将超声探头置于x=L的一端的端面上,采用垂直入射方式激发脉冲超声,由超声波脉冲回波法,获得被测试件ti时刻的超声波传播测量时间ti,exp,即在步骤二所述加热初始至结束的过程中进行多次测量,下标i表示测量时间序数,i=1,2,……,n,共测量n次;
步骤四:基于步骤一的关系式V(T),和步骤三的超声传播测量时间ti,exp,建立被测构件内部热源的热传导反问题目标函数为:
Figure BDA0001838226640000051
式(1)中,x表示点热源的位置,即沿被测方向x方向被测点热源的坐标值;Q(x,t)为点热源的强度,其随时间t而变化;ti,cal为超声波传播计算时间,是将ti时刻对应的被测试件加热端的温度值T(ti)代入步骤二所得的相关关系式V(T)得到的V,再根据L与V计算得到的;ti,exp为步骤三中实际测量得到的ti时刻的超声波传播测量时间;
所述热传导反问题目标函数的约束条件为:
Figure BDA0001838226640000052
T(x,t),t>0,x∈[0,L] (3)
式中,k(T)为被测构件材料的导热系数,为温度T的函数;Cp和ρ分别为被测构件材料的比热容和密度;T(x,t)表示所述加热初始至结束的过程中所述试件内部形成的温度梯度场中,沿被测方向x方向上一点的温度随时间t和坐标x的分布;T(x,t)|x=0/x=L为边界条件;
步骤五:由红外测温方法或热电偶测温方法,获得T(x,t)|x=0/x=L的边界温度;
步骤六:根据步骤五所得的温度边界条件T(x,t)|x=0/x=L,求解上述公式(1)和(2),获得x和t,从而测出构件内部热源沿被测方向的强度和位置;点热源强度和位置的测量误差均小于0.01%。
实施例2:被测构件为长方体,长度为1m,宽度为0.1m,设置两个点热源,设长度方向为x方向,其中,长方体的中心位置即x=0.5m,y=0.05m处设置点热源Q=1e+6W/m3,距离x=0的端面0.3m位置即x=0.3m,y=0.05m处设置点热源Q=1.2e+6W/m3。设定宽度方向即y方向为检测方向,从被测构件y=0处的一端进行加热,将超声探头置于y=1m的端面上,采用垂直入射方式激发脉冲超声,沿x方向移动超声探头,进行步进超声扫查检测,依次按照前述步骤一至步骤六,基于测量底面回波传播时间的变化,通过求解热源的热传导反问题,获得两个内部热源的强度和位置。内部热源位置测量结果分别为x=0.29m和x=0.5m,而两个点热源强度误差分别为3.75%和5.21%。
实施例3:同实施例1的被测构件,中心位置设置单个点热源Q(t)=133928.57t2+1848214.29t+9250000。将超声探头置于x=L的一端的端面上,采用垂直入射方式激发脉冲超声,依次按照前述步骤一至步骤六,基于测量底面回波传播时间的变化,通过求解热源的热传导反问题,获得随时间变化内部热源的强度如图2所示。在6个测量点处将直值与测量值比较,平均误差小于1.39%。
上述实施方式为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合和/或简化,均应视为等效的置换方式,均应包含在本发明的保护范围之内。

Claims (4)

1.一种基于超声的内部热源强度大小和位置的无损测量方法,其特征在于,具体步骤包括:
步骤一:获取被测构件内超声传播速度V与不同温度条件下温度T的相关关系,拟合出V与T的关系式V(T),其中,T为被测试件整体温度达到一致的条件下的温度值;
步骤二:对被测构件进行加热,选定某一被测方向,设定其为x方向,L为被测试件沿被测方向的长度,L大于一个超声波波长;沿被测方向x方向从被测构件x=0处的一端进行加热,直至被测试件所述x=0处的加热端的温度高于测量的预定温度时,结束加热;
在所述加热初始至结束的过程中,进行以下步骤三至步骤五,随后,基于步骤三至步骤五的结果,进行步骤六;
步骤三:由超声波脉冲回波法,获得被测试件ti时刻的超声波传播测量时间ti,exp,即在步骤二所述加热初始至结束的过程中进行多次测量,下标i表示测量时间序数,i=1,2,……,n,共测量n次;
步骤四:基于步骤一的关系式V(T),和步骤三的超声传播测量时间ti,exp,建立被测构件内部热源的热传导反问题目标函数为:
Figure FDA0001838226630000011
式(1)中,x表示点热源的位置,即沿被测方向x方向被测点热源的坐标值;Q(x,t)为点热源的强度,其随时间t而变化;ti,cal为超声波传播计算时间,是将ti时刻对应的被测试件加热端的温度值T(ti)代入步骤二所得的相关关系式V(T)得到的V,再根据L与V计算得到的;ti,exp为步骤三中实际测量得到的ti时刻的超声波传播测量时间;
所述热传导反问题目标函数的约束条件为:
Figure FDA0001838226630000021
T(x,t),t>0,x∈[0,L] (3)
式中,k(T)为被测构件材料的导热系数,为温度T的函数;Cp和ρ分别为被测构件材料的比热容和密度;T(x,t)表示所述加热初始至结束的过程中所述试件内部形成的温度梯度场中,沿被测方向x方向上一点的温度随时间t和坐标x的分布;T(x,t)|x=0/x=L为边界条件;
步骤五:由红外测温方法或热电偶测温方法,获得T(x,t)|x=0/x=L的边界温度;
步骤六:根据步骤五所得的温度边界条件T(x,t)|x=0/x=L,求解上述公式(1)和(2),获得x和t,从而测出构件内部热源沿被测方向的强度和位置。
2.根据权利要求1所述的基于超声的内部热源强度和位置的无损测量方法,其特征在于,该方法能够测量被测构件内部已知或未知位置的单个或多个热源的强度。
3.根据权利要求1所述的基于超声的内部热源强度和位置的无损测量方法,其特征在于,该方法适用于金属或复合材料的被测构件。
4.根据权利要求1所述的基于超声的内部热源强度和位置的无损测量方法,其特征在于,步骤一所述V与T的关系式V(T)是线性关系式或二次函数关系式。
CN201811236326.9A 2018-10-23 2018-10-23 一种基于超声的内部热源强度大小和位置的无损测量方法 Active CN109470772B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811236326.9A CN109470772B (zh) 2018-10-23 2018-10-23 一种基于超声的内部热源强度大小和位置的无损测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811236326.9A CN109470772B (zh) 2018-10-23 2018-10-23 一种基于超声的内部热源强度大小和位置的无损测量方法

Publications (2)

Publication Number Publication Date
CN109470772A CN109470772A (zh) 2019-03-15
CN109470772B true CN109470772B (zh) 2021-05-07

Family

ID=65664020

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811236326.9A Active CN109470772B (zh) 2018-10-23 2018-10-23 一种基于超声的内部热源强度大小和位置的无损测量方法

Country Status (1)

Country Link
CN (1) CN109470772B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110133101B (zh) * 2019-05-27 2021-05-28 东北大学 一种纤维增强复合材料板高温动力学性能退化分析方法
CN113639911B (zh) * 2021-07-28 2022-06-10 北京理工大学 一种用于测量含能药柱表层周向压应力的超声无损检测方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101825497B (zh) * 2010-03-30 2011-10-26 华南师范大学 一种基于热声效应的温度实时测量系统与方法
JP5957425B2 (ja) * 2013-08-06 2016-07-27 株式会社Ihi検査計測 内部付着物の厚さ計測装置と方法
CN104596667B (zh) * 2015-01-05 2017-12-01 中国空气动力研究与发展中心计算空气动力研究所 超声波探测物体内部瞬态非均匀温度场的灵敏度方法
CN106768466B (zh) * 2016-11-16 2019-03-29 石友安 一种基于超声波的瞬态转捩热流的无损探测方法
CN107462637A (zh) * 2017-08-31 2017-12-12 北京金风慧能技术有限公司 超声波回波检测方法及装置
CN107505395A (zh) * 2017-08-31 2017-12-22 北京金风慧能技术有限公司 工件内部损伤检测方法及装置
CN108051472B (zh) * 2017-12-05 2020-08-14 中国空气动力研究与发展中心计算空气动力研究所 一种材料高温热物性参数的快速测量方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Reconstruction of internal temperature distributions in heat materials";Dong Wei.et al;《Applied Thermal Engineering》;20161013;全文 *
"固体结构内部瞬态非均匀温度场的重建方法研究";石友安;《兵工学报》;20161231;全文 *

Also Published As

Publication number Publication date
CN109470772A (zh) 2019-03-15

Similar Documents

Publication Publication Date Title
WO2012167403A1 (zh) 脉冲红外热波技术测厚方法
WO2016101903A1 (zh) 一种导热系数测定装置
WO2019184076A1 (zh) 一种基于虚拟热源原理测量有限厚度材料异质含量的方法
Kruczek et al. In situ measurement of thermal diffusivity in anisotropic media
CN109470772B (zh) 一种基于超声的内部热源强度大小和位置的无损测量方法
CN104040327A (zh) 用于测量热导率的方法
CN109324079B (zh) 一种基于超声的材料热膨胀系数的测量方法
CN108008022B (zh) 一种随温度变化的超声波传播速度测量方法
Yamada et al. Noncontact monitoring of surface temperature distribution by laser ultrasound scanning
CN108344390A (zh) 基于k范围温度变化斜率曲线的涂层厚度检测方法
Schlichting et al. Defect sizing by local excitation thermography
CN106124078A (zh) 一种采用双热电偶测量强瞬变流体温度的方法
WO2011065877A1 (en) Measurement of thermal properties
JP3670869B2 (ja) コーティング層熱抵抗測定法
Ihara et al. New ultrasonic thermometry and its applications to temperature profiling of heated materials
RU2659617C1 (ru) Термографический способ контроля объектов и устройство для его осуществления
CN107966472B (zh) 一种高温接触热阻的无损快速测量方法
Silva et al. Analyzing the influence of thermal NDT parameters on test performance
RU2415408C1 (ru) Способ определения теплофизических характеристик теплозащитных покрытий
Kubicar et al. Transient methods for the measurement of thermophysical properties: The pulse transient method
Krapez et al. Time-resolved pulsed stimulated infrared thermography applied to carbon-epoxy non destructive evaluation
CN109613054A (zh) 一种直接通电纵向导热系数测试方法
CN107356627B (zh) 一种基于虚拟热源采用四参数匹配测定材料异质含量的方法
RU2478939C1 (ru) Способ измерения коэффициента температуропроводности теплоизоляционных материалов методом регулярного режима третьего рода
RU2462703C2 (ru) Способ определения тепловой активности материалов и устройство для его осуществления

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