CN113587992A - 固体材料预紧力和温度的超声波双波测量方法、用途及设备 - Google Patents

固体材料预紧力和温度的超声波双波测量方法、用途及设备 Download PDF

Info

Publication number
CN113587992A
CN113587992A CN202110864700.5A CN202110864700A CN113587992A CN 113587992 A CN113587992 A CN 113587992A CN 202110864700 A CN202110864700 A CN 202110864700A CN 113587992 A CN113587992 A CN 113587992A
Authority
CN
China
Prior art keywords
temperature
wave
solid material
ultrasonic
pretightening force
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
CN202110864700.5A
Other languages
English (en)
Other versions
CN113587992B (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.)
Lvyin Technology Wuhan Co Ltd
Original Assignee
Lvyin Technology Wuhan Co Ltd
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 Lvyin Technology Wuhan Co Ltd filed Critical Lvyin Technology Wuhan Co Ltd
Priority to CN202110864700.5A priority Critical patent/CN113587992B/zh
Publication of CN113587992A publication Critical patent/CN113587992A/zh
Priority to EP21951633.3A priority patent/EP4365557A1/en
Priority to PCT/CN2021/130376 priority patent/WO2023005061A1/zh
Application granted granted Critical
Publication of CN113587992B publication Critical patent/CN113587992B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D21/00Measuring or testing not otherwise provided for
    • G01D21/02Measuring two or more variables by means not covered by a single other subclass
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B1/00Measuring instruments characterised by the selection of material therefor
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L1/00Measuring force or stress, in general
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/72Wind turbines with rotation axis in wind direction

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Force Measurement Appropriate To Specific Purposes (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明公开了一种固体材料预紧力和温度的超声波双波测量方法,包括:获取固体材料在第一温度和应力状态下超声波双波的第一模态波声时、第二模态波声时;测量固体材料在未知状态下超声波双波的第一模态波时间和第二模态波时间;获取固体材料的第一模态波温度影响系数、第二模态波温度影响系数,以及相同规格和几何形状被测固体材料中的第一模态波预紧力影响系数、第二模态波预紧力影响系数;根据一阶或者二阶泰勒近似得到超声波双波声时相对变化的关系,联合解得测量预紧力、测量温度。本发明通过当前测量的双波声时及一些常数可同时得到预紧力与平均温度,计算方法简单,效率高,可以避免多种测量误差,具有广泛的应用。

Description

固体材料预紧力和温度的超声波双波测量方法、用途及设备
技术领域
本发明涉及超声波双波测量技术领域,具体地涉及一种固体材料预紧力和温度的超声波双波测量方法、用途及设备,尤其适用于螺栓轴力多次巡回检测。
背景技术
超声波监测检测固体材料(例如螺栓、金属等)预紧力方法是通过向固体材料发射超声波,测量超声的声时差来推导固体材料的预紧力变化。超声声时差的变化非常容易受到固体材料温度以及温度不均匀的影响。
以风电螺栓为例,风电螺栓的预紧力检测和监测是保证风电叶片和塔筒安全的重要手段。风电螺栓预紧力的监测和检测通常用超声波测量声时的方法配铂铑温度传感器Pt1000来测量螺栓上某点的温度,用这个测量的温度值来补偿超声波声速随温度的变化。由于现场温度的不均匀,导致螺栓本身温度的不均匀,螺栓的平均温度很难通过温度计单点准确测量得到,这样容易引起测量预紧力的温度补偿误差。例如对于长达数米的风机塔筒锚栓,其地面露出部分与地下部分的环境差异,都可能造成整个锚栓长度上的温度非常不均匀。又比如风机塔筒内部螺栓由于塔筒外阳光照射位置不同以及塔筒内部电器运行过程中的热辐射方向不均匀等原因,会发生整根螺栓温度分布不均匀的现象。
风电螺栓预紧力的超声波巡回检测通常是比螺栓预紧力超声在线监测成本更低的方法,但是这种方法要求每次测量螺栓预紧力时都记录螺栓的平均温度用于补偿测量螺栓预紧力的测量值,但是现有的方法中无法快速得到螺栓的平均温度,这在实际应用中非常不方便。
风电螺栓预紧力在线监测的超声波测量方法通常采用压电单波探头,每根螺栓配一个探头,每个探头需要配一只温度计。这增加了螺栓预紧力超声波方法监测设备的复杂性,而且配置的温度计通常只能测量螺栓上某一点的温度,无法获取螺栓整体的平均温度,而平均温度才是计算温度补偿需要用到的准确温度。一旦螺栓处于温度变化较快的环境,螺栓上的温度会变得非常不均匀,温度计的测量误差也会随之增大。
公告号为CN 111693190 A的发明专利公开了一种基于超声波的螺栓轴向应力测量装置,包括超声波测量装置,测量超声波在螺栓中的声时差;温度测量模块,测量螺栓温度;系数标定模块,通过所述超声波测量装置和温度测量模块至少测量三次超声波的声时差,根据矩阵公式确定系数;微处理器根据所述系数确定螺栓温度、轴向应力和声时差的函数关系,基于该函数关系确定所述螺栓的当前轴向应力。该方法使用超声波单波进行测量,并且至少测量三次超声波的声时差,需要多次测量螺栓温度,复杂度相对较高,成本也较高。
发明内容
针对上述存在的技术问题,本发明目的是:提供了一种固体材料预紧力和温度的超声波双波测量方法、用途及设备,通过直接使用超声波单一种类的方法同时测量了固体材料预紧力以及温度,排除了温度整体变化及温度不均匀变化对预紧力测量的影响。无需使用额外温度传感器测量固体材料温度来进行补偿,从而避免由于温度传感器测量误差、固体材料温度不均匀导致的温度测量误差引入预紧力测量误差,极大地提高了固体材料预紧力超声波巡检的速度与便捷性。同时,本发明也使用超声波这种单一的方法较高精度地代替了传统的温度计对固体温度的测量。
本发明的技术方案是:
一种固体材料预紧力和温度的超声波双波测量方法,包括以下步骤:
S01:获取固体材料在第一种状态下的温度T0及预紧力F0,以及在该状态下超声波双波的第一模态波声时TS0、第二模态波声时TL0
S02:测量固体材料在未知状态下超声波双波的第一模态波时间TS和第二模态波时间TL
S03:获取固体材料的第一模态波温度影响系数KTS、第二模态波温度影响系数KTL,以及相同规格和几何形状被测固体材料中的第一模态波预紧力影响系数KFS、第二模态波预紧力影响系数KFL
S04:根据一阶泰勒近似得到超声波双波声时相对变化的关系为:
(TS-TS0)/TS0=KTS·(T-T0)+KFS·(F-F0) (I)
(TL-TL0)/TL0=KTL·(T-T0)+KFL·(F-F0) (II)
其中,F为测量预紧力,T为测量温度;
联合公式(I)和(II)解得测量预紧力F:
Figure BDA0003187191410000021
解得测量温度T:
Figure BDA0003187191410000031
优选的技术方案中,所述步骤S04之后还包括:
S05:获取固体材料的第一模态波另一温度影响系数KTS2、第二模态波另一温度影响系数KTL2,以及相同规格和几何形状被测固体材料中的第一模态波另一预紧力影响系数KFS2、第二模态波另一预紧力影响系数KFL2
S06:根据二阶泰勒近似得到超声波双波声时相对变化的关系为:
(TS-TS0)/TS0=KTS·(T-T0)+KTS2·(T-T0)2+KFS·(F-F0)+KFS2·(F-F0)2 (III)
(TL-TL0)/TL0=KTL·(T-T0)+KTL2·(T-T0)2+KFL·(F-F0)++KFL2·(F-F0)2 (IV)
S07:联合公式(III)和(IV)解得测量预紧力F和测量温度T。
优选的技术方案中,所述第一模态波和所述第二模态波为横波、纵波、表面波、导波、爬波或以上任意两种模态的组合波。
优选的技术方案中,所述步骤S01中的数据通过测量得到或者查找记录数据得到。
优选的技术方案中,所述步骤S03和S05中的数据通过标定数据库查找得到。
将上述的固体材料预紧力和温度的超声波双波测量方法用于高强螺栓预紧力的定期超声波巡检,测量螺栓的预紧力和/或温度。
将上述的固体材料预紧力和温度的超声波双波测量方法用于测量自由拉力状态下金属在高温状态下的轴向平均温度。
将上述的固体材料预紧力和温度的超声波双波测量方法用于螺栓轴力的在线监测可以为每个在线监测超声波探头省掉一支铂铑温度计。
将所述的固体材料预紧力和温度的超声波双波测量方法用于测量压力容器外壳的轴向平均预紧力及轴向平均温度,进而推算出压力容器内部的温度。
本发明还公开了一种固体材料预紧力和温度的超声波双波测量设备,包括:
数据获取/测量模块,获取固体材料在第一状态下的温度T0及预紧力F0,及在该状态下超声波双波的第一模态波声时TS0、第二模态波声时TL0、;
测量模块,测量固体材料在未知状态下超声波双波的第一模态波时间TS和第二模态波时间TL
参数获取模块,获取固体材料的第一模态波温度影响系数KTS、第二模态波温度影响系数KTL,以及相同规格和几何形状被测固体材料中的第一模态波预紧力影响系数KFS、第二模态波预紧力影响系数KFL
建模计算模块,根据一阶泰勒近似得到超声波双波声时相对变化的关系为:
(TS-TS0)/TS0=KTS·(T-T0)+KFS·(F-F0) (I)
(TL-TL0)/TL0=KTL·(T-T0)+KFL·(F-F0) (II)
其中,F为测量预紧力,T为测量温度;
联合公式(I)和(II)解得测量预紧力F:
Figure BDA0003187191410000041
解得测量温度T:
Figure BDA0003187191410000042
优选的技术方案中,还包括精度提升模块,其包括:
第二参数获取模块,获取固体材料的第一模态波另一温度影响系数KTS2、第二模态波另一温度影响系数KTL2,以及相同规格和几何形状被测固体材料中的第一模态波另一预紧力影响系数KFS2、第二模态波另一预紧力影响系数KFL2
第二建模计算模块,根据二阶泰勒近似得到超声波双波声时相对变化的关系为:
(TS-TS0)/TS0=KTS·(T-T0)+KTS2·(T-T0)2+KFS·(F-F0)+KFS2·(F-F0)2 (III)
(TL-TL0)/TL0=KTL·(T-T0)+KTL2·(T-T0)2+KFL·(F-F0)+KFL2·(F-F0)2 (IV)
联合求解模块,联合公式(III)和(IV)解得测量预紧力F和测量温度T。
优选的技术方案中,所述第一模态波和所述第二模态波为横波、纵波、表面波、导波、爬波或以上任意两种模态的组合波。
优选的技术方案中,所述数据获取/测量模块中的数据通过测量得到或者查找记录数据得到。
与现有技术相比,本发明的优点是:
1、本发明的固体材料预紧力和温度的超声波双波测量方法,只需要通过采用当前状态测量的双波声时,对比状态时的声波声时以及一些常数就可以同时得到预紧力的变化与平均温度的变化,计算方法简单,效率高,具有广泛的应用。通过算法直接使用超声波单一种类的方法同时测量了螺栓预紧力以及温度,排除了温度整体变化及温度不均匀变化对预紧力测量的影响。无需使用额外温度传感器测量螺栓温度来进行补偿,从而避免由于温度传感器测量误差、螺栓温度不均匀导致的温度测量误差引入预紧力测量误差。
2、当温度的测量成为目标时,本发明也提供了一种抑制应力状态对温度测量影响的一种方法。这种方法可以推广到螺栓以外的物体或者材料仅仅用超声波单一种类方法来测量样品的温度。
3、当需要同时测量螺栓的轴向平均预紧力和平均温度时,本方法提供了可靠的超声波手段可以同时测量螺栓的轴向平均预紧力和平均温度。该方法也适用于其它可以传播超声波的物体,比如可以用本发明提供的超声波双波方法来同时测量压力容器外壳在某一时刻的平均温度和平均轴向预紧力。当采用适当的压力容器外壳轴向应力模型以及温度梯度模型时,本发明提供的超声波单一种类的方法可以以较高精度同时测量压力容器内部的温度和压力。
附图说明
下面结合附图及实施例对本发明作进一步描述:
图1是本发明固体材料预紧力和温度的超声波双波测量方法的流程图;
图2是本发明固体材料预紧力和温度的超声波双波测量设备的原理框图;
图3是本发明超声测量温度与热电偶温度计测量几百摄氏度大范围温度的对比图。
图4是本发明超声测量温度与铂铑温度计测量温度的对比图,测量的温度用于本发明超声波测量固体材料预紧力的温度补偿。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明了,下面结合具体实施方式并参照附图,对本发明进一步详细说明。应该理解,这些描述只是示例性的,而并非要限制本发明的范围。此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要地混淆本发明的概念。
如图1所示,一种固体材料预紧力和温度的超声波双波测量方法,包括以下步骤:
S01:获取固体材料在第一状态下温度T0及预紧力F0,以及在该状态下超声波双波的第一模态波声时TS0、第二模态波声时TL0
S02:测量固体材料在未知状态下超声波双波的第一模态波时间TS和第二模态波时间TL
S03:获取固体材料的第一模态波温度影响系数KTS、第二模态波温度影响系数KTL,以及相同规格和几何形状被测固体材料中的第一模态波预紧力影响系数KFS、第二模态波预紧力影响系数KFL
S04:根据一阶泰勒近似得到超声波双波声时相对变化的关系为:
(TS-TS0)/TS0=KTS·(T-T0)+KFS·(F-F0) (I)
(TL-TL0)/TL0=KTL·(T-T0)+KFL·(F-F0) (II)
其中,F为测量预紧力,T为测量温度;
联合公式(I)和(II)解得测量预紧力F:
Figure BDA0003187191410000061
解得测量温度T:
Figure BDA0003187191410000062
一较佳的实施例中,当需要进一步提高S04的超声波温度测量和预紧力测量精度时,所述步骤S04之后还包括:
S05:获取固体材料的第一模态波另一温度影响系数KTS2、第二模态波另一温度影响系数KTL2,以及相同规格和几何形状被测固体材料中的第一模态波另一预紧力影响系数KFS2、第二模态波另一预紧力影响系数KFL2
S06:根据二阶泰勒近似得到超声波双波声时相对变化的关系为:
(TS-TS0)/TS0=KTS·(T-T0)+KTS2·(T-T0)2+KFS·(F-F0)+KFS2·(F-F0)2 (III)
(TL-TL0)/TL0=KTL·(T-T0)+KTL2·(T-T0)2+KFL·(F-F0)+KFL2·(F-F0)2 (IV)
S07:选择联合公式(III)和(IV)解得测量预紧力F和测量温度T。
需要说明的是:步骤S07中在数学上有2组解,但是在物理上通常只有一组有实际物理意义。这组解比步骤S04中解出的测量预紧力和测量温度更加准确,当然计算也更加复杂些。当温度变化比较大(可以通过与阈值比较得到)或轴力变化范围比较大(可以通过与阈值比较得到)时,可以采用S05-S07的步骤来提高测量准确度。
这里的第一模态波和所述第二模态波为横波、纵波、表面波、导波、爬波或者以上两种模态波中任意两种组合的波。一般常用为横波和纵波,或者横波和横波纵波的组合波,或者纵波和横纵波的组合波。
这里的预紧力可以为轴力、载荷或者应力等等,当然应力在计算时需要乘以平均横截面积,以让公式中的系数保持不变。
步骤S01中的数据可以通过数据库查找得到,或者由数据记录文件得到,或者进行一次超声波双波测量得到。
步骤S03和S05中的数据通过标定数据库查找得到,KTS、KTS2、KFS、KFS2、KTL、KTL2、KFL、KFL2等都与具体材料和样品相关,是可以标定的比例常数。
当将上述方法用于高强螺栓预紧力的定期超声波巡检测量螺栓的预紧力和/或温度时,一般都可以通过数据库查到。以上公式说明巡检时螺栓的预紧力可以只通过当前测量的双波声时以及一些常数就可以得到。这样就免去了使用温度计逐点测量温度的麻烦,大大提高了螺栓预紧力超声波巡检的效率。
还可以将上述的固体材料预紧力和温度的超声波双波测量方法用于测量自由拉力状态下金属在高温状态下的轴向平均温度。
还可以将上述的固体材料预紧力和温度的超声波双波测量方法用于螺栓预紧力超声波在线监测过程中的温度测量,省去了通常的每颗螺栓都需要的铂铑温度计硬件,提高了本发明方法固体材料预紧力的温度补偿精度。
还可以将上述的固体材料预紧力和温度的超声波双波测量方法用于测量压力容器外壳的轴向平均预紧力及轴向平均温度,辅以适当的压力容器外壳轴向应力模型以及温度梯度模型,进而可以推算出压力容器内部的温度和压力。
如图2所示,本发明还公开了一种固体材料预紧力和温度的超声波双波测量设备,包括:
数据获取/测量模块10,获取固体材料在第一状态下的温度T0及预紧力F0,及该状态下超声波双波的第一模态波声时TS0、第二模态波声时TL0
测量模块20,测量固体材料在未知状态下超声波双波的第一模态波时间TS和第二模态波时间TL
参数获取模块30,获取固体材料的第一模态波温度影响系数KTS、第二模态波温度影响系数KTL,以及相同规格和几何形状被测固体材料中的第一模态波预紧力影响系数KFS、第二模态波预紧力影响系数KFL
建模计算模块40,根据一阶泰勒近似得到超声波双波声时相对变化的关系为:
(TS-TS0)/TS0=KTS·(T-T0)+KFS·(F-F0) (I)
(TL-TL0)/TL0=KTL·(T-T0)+KFL·(F-F0) (II)
其中,F为测量预紧力,T为测量温度;
联合公式(I)和(II)解得测量预紧力F:
Figure BDA0003187191410000081
解得测量温度T:
Figure BDA0003187191410000082
较佳的,当需要进一步提高由一阶泰勒近似得到的结果的精度时,还包括精度提升模块,其包括:
第二参数获取模块,获取固体材料的第一模态波另一温度影响系数KTS2、第二模态波另一温度影响系数KTL2,以及相同规格和几何形状被测固体材料中的第一模态波另一预紧力影响系数KFS2、第二模态波另一预紧力影响系数KFL2
第二建模计算模块,根据二阶泰勒近似得到超声波双波声时相对变化的关系为:
(TS-TS0)/TS0=KTS·(T-T0)+KTS2·(T-T0)2+KFS·(F-F0)+KFS2·(F-F0)2 (III)
(TL-TL0)/TL0=KTL·(T-T0)+KTL2·(T-T0)2+KFL·(F-F0)+KFL2·(F-F0)2 (IV)
联合求解模块,联合公式(III)和(IV)解得测量预紧力F和测量温度T。在数学上有2组解,但是在物理上通常只有一组有实际物理意义。
下标带0的数据可以通过数据库查找得到,或者由数据记录文件得到,或者进行一次超声波双波测量得到。其它常系数可以通过标定或者从数据库中得到。
下面以高强螺栓为例说明一下测量方法:
一种高强螺栓预紧力与温度的超声双波巡检方法,包括以下步骤:
S01:将需要巡检的螺栓编号。
S02:获取螺栓巡检前记录的每一颗螺栓以前任何一次双波测量预紧力时的横波声时TS0、纵波声时TL0、当时温度计测量出温度T0以及由当次测量得出的螺栓预紧力值F0;如果没有以前的记录,就用超声双波方法现场测量出这一组数据。
S03:测量螺栓在巡检时刻的超声横波回波时间TS和纵波回波时间TL。假设巡检时真实的温度为T,真实的预紧力为F。
S04:计算两次或者两个时刻的超声横波回波变化时间ΔTS=TS-TS0和超声纵波回波变化时间ΔTL=TL-TL0。两次温度变化记为ΔT=T-T0,两次预紧力值变化记为ΔF=F-F0
S05:在螺栓标定数据库中查找螺栓材料的横波温度影响系数KTS、纵波温度影响系数KTL,以及相同规格和几何形状被测螺栓中的横波预紧力影响系数KFS、纵波预紧力影响系数KFL。这些参数都是可以提前标定的材料学参数(也与螺栓的几何形状以及螺母的装配位置有关),在一阶泰勒的数学近似下比较完美地满足下面的关系式:
ΔTS/Ts0=KTS·ΔT (1)
ΔTL/TL0=KTL·ΔT (2)
ΔTS/TS0=KFS·ΔF (3)
ΔTL/TL0=KFL·ΔF (4)
其中Δ表示本次测量的量与测量前获得的对应的量之差。
S06:假设每一根具体的螺栓超声波双波声时的相对变化量只跟两个最显著的因素有关:温度和预紧力的变化。所以该具体的螺栓超声波双波声时相对变化在一阶泰勒近似下的数学表达式为:
ΔTS/TS0=KTS·ΔT+KFS·ΔF (5)
ΔTL/TL0=KTL·ΔT+KFL·ΔF (6)
求解上面2个表达式得出,
Figure BDA0003187191410000101
具体的分解算术表达式为,
Figure BDA0003187191410000102
Figure BDA0003187191410000103
螺栓巡检时当前绝对的预紧力和温度为:
Figure BDA0003187191410000104
Figure BDA0003187191410000105
S07:当需要进一步提高S06的超声波温度测量和预紧力测量精度时,根据二阶泰勒近似得到超声波双波声时相对变化的关系为:
(TS-TS0)/TS0=KTS·(T-T0)+KTS2·(T-T0)2+KFS·(F-F0)+KFS2·(F-F0)2 (12)
(TL-TL0)/TL0=KTL·(T-T0)+KTL2·(T-T0)2+KFL·(F-F0)++KFL2·(F-F0)2 (13)
其中,F为测量预紧力,T为测量温度。联合公式(12)和(13)解得测量预紧力F和测量温度T。在数学上有2组解,但是在物理上通常只有一组有实际物理意义。
下面以具体的实施例进行进一步的说明和验证。
实施例1
有某一厂家的高强度六角头螺栓M36X260,螺母安装好后,螺栓露出螺母12mm。已知螺栓材料为42CrMoA,这种材料标定的温度系数为KTS=1.35×10-4/℃,KTL=9.5×10-5/℃。对这种螺栓已按照安装后露出螺母12mm的夹持长度进行拉力标定,标定结果为KFS=6.49×10-6/KN,KFL=1.68×10-5/KN。
在恒定室温25℃环境下,对某一没有受轴力的螺栓进行第一次超声测量,测得超声横波回波时间TS0=175378ns和纵波回波时间TL0=95848ns。螺栓温度T0=25℃被认为等于室温。然后使用拉力实验机独立对螺栓进行轴力加载F’=350KN。
然后进行第二次超声测量,测得超声横波回波时间TS=175776ns和纵波回波时间TL=96410ns。所以ΔTS=TS-TS0=398ns,ΔTL=TL-TL0=562ns。根据公式(10)和(11)计算得出:
F=348.8KN
T=25.04℃
对比可知,计算出来的轴力值与拉力机上独立加载的轴力值几乎一致,计算出的螺栓平均温度值与温度计独立测量出的温度几乎一致。
实施例2
另一厂家的高强度六角头螺栓M42X240,螺母安装好后,螺栓露出螺母25mm。已知螺栓材料为42CrMoA,这种材料标定的KTS=1.35×10-4/℃,
KTL=9.5×10-5/℃。对这种螺栓已按照安装后露出螺母25mm的夹持长度进行拉力标定,标定结果为KFS=3.67×10-6/KN,KFL=9.93×10-6/KN。
在室温环境下,上述螺栓样品不受轴向拉力,即F0=0KN,对螺栓样品进行一发一收第一次超声测量,测得超声横波回波时间TS0=82520ns和纵波回波时间TL0=45092ns,通过温度传感器测量螺栓表面温度T0=21℃。
试验条件1:然后将螺栓样品放入-25℃的恒温箱内,此时螺栓实际上仍然不受轴向拉力,即我们用独立的方法获得F’=0KN。使用温度传感器独立地获得螺栓整体温度,T’=-25℃,这个独立的轴力值和独立的温度值用作稍后计算结果的参考对比。
此时进行第二次一发一收超声测量,测得超声横波回波时间TS=81990ns和纵波回波时间TL=44888ns。所以ΔTS=TS-TS0=-530ns,ΔTL=TL-TL0=-204ns。根据公式(10)和(11)计算出结果为:
F=-0.6KN
T=-26.6℃
计算得到的数值与独立获得的数值F’=0KN和T’=-25℃,相差不大,在允许范围内。
试验条件2:然后将螺栓样品放入60℃左右的恒温箱内放置一段时间,此时螺栓仍然不受轴向拉力,即我们用独立的方法获得F’=0KN,用温度计独立测量出螺栓的温度为T’=61℃。
此时进行第三次一发一收超声测量,测得超声横波回波时间TS=83000ns和纵波回波时间TL=45282ns。所以ΔTS=TS-TS0=480ns,ΔTL=TL-TL0=190ns。根据公式(10)和(11)计算出结果为:
F=16.4KN
T=63.6℃
上面计算得到的数值与独立方法获得的数值F’=0KN以及T’=61℃相差不大,在允许范围内。
实施例3
同实例2的高强度六角头螺栓M42X240,螺母安装好后,螺栓露出螺母25mm。已知螺栓材料为42CrMoA,这种材料标定的KTS=1.35×10-4/℃,KTL=9.5×10-5/℃。对这种螺栓按照安装后露出螺母25mm的夹持长度进行拉力标定,标定结果为KFS=3.67×10-6/KN,KFL=9.93×10-6/KN。
在室温21℃环境下,上述螺栓样品放置到拉伸台,已受静态轴向拉力F0=444KN,对该螺栓做第一次自发自收超声波测量测得超声横波回波时间TS0=165052ns和纵波回波时间TL0=90472ns,螺栓温度取环境温度T0=21℃。
然后对螺栓样品的螺母进行部分松卸,并同时使用热风枪加热螺栓中间局部区域,此时拉伸台独立的传感器显示螺栓的轴向拉力为F’=348KN,我们以此独立的螺栓轴力作为后面用超声波方法计算的轴力值的参考。此时,螺栓中间加热部位使用温度传感器独立测量其局部温度为T’=62.6℃。
此时进行第二次自发自收超声测量,测得超声横波回波时间TS=165172ns和纵波回波时间TL=90452ns。所以ΔTS=TS-TS0=-120ns,ΔTL=TL-TL0=20ns。利用公式(10)和(11)计算得到:
F=344.3KN
T=29.1℃
计算出的螺栓拉力值与独立数值F’=348KN相差不大,在允许范围内。此实例证明即使螺栓出现局部温度的非常不均匀和温度突变,通过本发明方法仍然可以测得准确的螺栓轴力变化。T和T’相差较大,是因为T为超声测量的螺栓平均温度,T’为温度传感器测量的螺栓加热部位的局部温度。
实施例4
又有另一厂家的高强度六角头螺栓M42X255,螺母安装好后,夹持长度196mm。已知螺栓材料为42CrMoA,这种材料标定的KTS=1.35×10-4/℃,KTL=9.5×10-5/℃。对这种螺栓已按照安装后夹持长度为196mm进行拉力标定,标定结果为KFS=3.86×10-6/KN,KFL=1.05×10-5/KN。
在室温环境下,上述螺栓样品不受轴向拉力,即F0=0KN,对螺栓样品进行第一次自发自收超声测量,测得超声横波回波时间TS0=173578ns和纵波回波时间TL0=94860ns,通过温度传感器测量螺栓表面温度T0=26.6℃。
然后将螺栓按照夹持长度为196mm在扭矩实验台上装好,使用液压扭矩扳手将螺栓拧紧。
扭矩实验台上有一个独立的力传感器用来显示螺栓当前的轴力,即我们用独立的方法获得F’=634.22KN。当前位于螺栓杆中部的温度传感器测得螺栓温度没变化,即T’=26.6℃。
此时进行第二次自发自收超声测量,测得超声横波回波时间TS=174028ns和纵波回波时间TL=95504ns。所以ΔTS=TS-TS0=450ns,ΔTL=TL-TL0=644ns。根据公式(10)和(11)计算出结果为:
F=637.8KN
T=27.55℃
计算得到的数值与独立获得的数值F’=634.22KN非常接近,在允许范围内。T和T’相差0.95℃的原因是螺栓扭矩法施工,螺母与螺栓摩擦产生热量,导致螺栓局部温度升高,0.95℃为螺栓整体温度变化量。
在上述扳手将螺栓拧紧过程中,同时使用超声单波方法实时测量螺栓拧紧过程的超声回波时间变化,从而达到实时测量螺栓轴力变化的目的。扳手停止时刻,横波变化时间和纵波变化时间ΔTS1=TS-TS0=450ns,ΔTL1=TL-TL0=644ns。根据公式(I)和(II),分别计算使用单横波和单纵波测量的螺栓轴力FS=671KN和FL=646KN。由此可得出结论,使用超声单波方法与扳手配合实时测量螺栓预紧过程的轴力值,结果会受到局部温度升高的影响。扳手紧固前后各进行一次双波测试,根据公式(10)计算得到的轴力不受局部温度变化影响。在与扳手配合中使用单波与双波相结合的方法,使得单波测量的实时性和双波测量的准确性相互结合。
实施例5
将一块初始温度为38.5℃的铁基高温合金,加热到350℃,然后慢慢冷却。用热电偶温度计和可耐高温的电磁超声波双波设备同时测量温度的变化,得到的热电偶温度数据和超声波温度测量一阶近似结果如图3所示。由于以上测量方法的步骤S06中的温度变化范围非常巨大,一阶泰勒近似出现明显偏差,为了提高超声波测量温度的精度,采用步骤S07的温度变化二阶泰勒近似,结果也见图3。图3中通过步骤S07的温度二阶泰勒近似以后,超声波测量的温度与热电偶温度计测量的温度在300摄氏度的范围内最大相差不过5摄氏度,这说明本发明超声波双波测量方法可以独立地大范围内测量自由状态下高温金属的温度。
实施例6
在一水电站采用超声双波设备监测一顶盖螺栓预紧力,相对处于风电塔筒环境的螺栓,水电站内螺栓的环境温度变化极其缓慢,正好可以作为普通温度计与本发明采用的超声波温度计做温度测量对比。图4是利用本发明方法得到的超声波双波温度计测量结果与现场Pt1000温度计测量的螺栓温度作的对比。由图可见,本发明提出的超声波双波方法测量的螺栓平均温度与用独立温度计测量的螺栓上某点的温度,或者说是环境温度一致性非常好,在环境温度变化10度的范围以内两种方法相互之间的测量误差没有超过0.5摄氏度。由于风电塔筒里面的螺栓温度一天之内变化较大较快,由通常采用的Pt1000温度计单点测量的温度难以真实地反映螺栓的平均温度,所以本发明为螺栓的超声在线监测也提供了可靠的、方法类别单一的以及操作方便的温度测量方法。
应当理解的是,本发明的上述具体实施方式仅仅用于示例性说明或解释本发明的原理,而不构成对本发明的限制。因此,在不偏离本发明的精神和范围的情况下所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。此外,本发明所附权利要求旨在涵盖落入所附权利要求范围和边界、或者这种范围和边界的等同形式内的全部变化和修改例。

Claims (10)

1.一种固体材料预紧力和温度的超声波双波测量方法,其特征在于,包括以下步骤:
S01:获取固体材料在第一状态下的温度T0和预紧力F0,以及在该状态下超声波双波的第一模态波声时TS0及第二模态波声时TL0
S02:测量固体材料在未知状态下超声波双波的第一模态波时间TS和第二模态波时间TL
S03:获取固体材料的第一模态波温度影响系数KTS、第二模态波温度影响系数KTL,以及相同规格和几何形状被测固体材料中的第一模态波预紧力影响系数KFS、第二模态波预紧力影响系数KFL
S04:根据一阶泰勒近似得到超声波双波声时相对变化的关系为:
(TS-TS0)/TS0=KTS·(T-T0)+KFS·(F-F0) (I)
(TL-TL0)/TL0=KTL·(T-T0)+KFL·(F-F0) (II)
其中,F为测量预紧力,T为测量温度;
联合公式(I)和(II)解得测量预紧力F:
Figure FDA0003187191400000011
解得测量温度T:
Figure FDA0003187191400000012
2.根据权利要求1所述的固体材料预紧力和温度的超声波双波测量方法,其特征在于,所述步骤S04之后还包括:
S05:获取固体材料的第一模态波另一温度影响系数KTS2、第二模态波另一温度影响系数KTL2,以及相同规格和几何形状被测固体材料中的第一模态波另一预紧力影响系数KFS2、第二模态波另一预紧力影响系数KFL2
S06:根据二阶泰勒近似得到超声波双波声时相对变化的关系为:
(TS-TS0)/TS0=KTS·(T-T0)+KTS2·(T-T0)2+KFS·(F-F0)+KFS2·(F-F0)2 (III)
(TL-TL0)/TL0=KTL·(T-T0)+KTL2·(T-T0)2+KFL·(F-F0)+KFL2·(F-F0)2 (IV)
S07:联合公式(III)和(IV)解得测量预紧力F和测量温度T。
3.根据权利要求1或2所述的固体材料预紧力和温度的超声波双波测量方法,其特征在于,所述第一模态波和所述第二模态波为横波、纵波、表面波、导波、爬波或以上任意两种模态的组合波。
4.根据权利要求1或2所述的固体材料预紧力和温度的超声波双波测量方法,其特征在于,所述步骤S01中的数据通过测量得到或者查找记录数据得到。
5.根据权利要求2所述的固体材料预紧力和温度的超声波双波测量方法,其特征在于,所述步骤S03和步骤S05中的数据通过标定数据库查找得到。
6.将权利要求1-5任一项所述的固体材料预紧力和温度的超声波双波测量方法用于高强螺栓预紧力的定期超声波巡检或者螺栓预紧力在线监测,同时测量螺栓的预紧力和温度。
7.将权利要求1-5任一项所述的固体材料预紧力和温度的超声波双波测量方法用于测量自由拉力状态下金属在高温状态下的轴向平均温度。
8.将权利要求1-4任一项所述的固体材料预紧力和温度的超声波双波测量方法用于测量压力容器外壳的轴向平均预紧力及轴向平均温度,进而推算出压力容器内部的温度。
9.一种固体材料预紧力和温度的超声波双波测量设备,其特征在于,包括:
数据获取/测量模块,获取固体材料在第一状态下温度T0及预紧力F0,以及该状态下超声波双波的第一模态波声时TS0及第二模态波声时TL0
测量模块,测量固体材料在未知状态下超声波双波的第一模态波时间TS和第二模态波时间TL
参数获取模块,获取固体材料的第一模态波温度影响系数KTS、第二模态波温度影响系数KTL,以及相同规格和几何形状被测固体材料中的第一模态波预紧力影响系数KFS、第二模态波预紧力影响系数KFL
建模计算模块,根据一阶泰勒近似得到超声波双波声时相对变化的关系为:
(TS-TS0)/TS0=KTS·(T-T0)+KFS·(F-F0) (I)
(TL-TL0)/TL0=KTL·(T-T0)+KFL·(F-F0) (II)
其中,F为测量预紧力,T为测量温度;
联合公式(I)和(II)解得测量预紧力F:
Figure FDA0003187191400000031
解得测量温度T:
Figure FDA0003187191400000032
10.根据权利要求9所述的固体材料预紧力和温度的超声波双波测量设备,其特征在于,还包括精度提升模块,其包括:
第二参数获取模块,获取固体材料的第一模态波另一温度影响系数KTS2、第二模态波另一温度影响系数KTL2,以及相同规格和几何形状被测固体材料中的第一模态波另一预紧力影响系数KFS2、第二模态波另一预紧力影响系数KFL2
第二建模计算模块,根据二阶泰勒近似得到超声波双波声时相对变化的关系为:
(TS-TS0)/TS0=KTS·(T-T0)+KTS2·(T-T0)2+KFS·(F-F0)+KFS2·(F-F0)2 (III)
(TL-TL0)/TL0=KTL·(T-T0)+KTL2·(T-T0)2+KFL·(F-F0)+KFL2·(F-F0)2 (IV)
联合求解模块,联合公式(III)和(IV)解得测量预紧力F和测量温度T。
CN202110864700.5A 2021-07-29 2021-07-29 固体材料预紧力和温度的超声波双波测量方法、用途及设备 Active CN113587992B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202110864700.5A CN113587992B (zh) 2021-07-29 2021-07-29 固体材料预紧力和温度的超声波双波测量方法、用途及设备
EP21951633.3A EP4365557A1 (en) 2021-07-29 2021-11-12 Method, use, and device for ultrasonic double-wave measurement of pre-tightening force and temperature of solid material
PCT/CN2021/130376 WO2023005061A1 (zh) 2021-07-29 2021-11-12 固体材料预紧力和温度的超声波双波测量方法、用途及设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110864700.5A CN113587992B (zh) 2021-07-29 2021-07-29 固体材料预紧力和温度的超声波双波测量方法、用途及设备

Publications (2)

Publication Number Publication Date
CN113587992A true CN113587992A (zh) 2021-11-02
CN113587992B CN113587992B (zh) 2022-07-12

Family

ID=78251882

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110864700.5A Active CN113587992B (zh) 2021-07-29 2021-07-29 固体材料预紧力和温度的超声波双波测量方法、用途及设备

Country Status (3)

Country Link
EP (1) EP4365557A1 (zh)
CN (1) CN113587992B (zh)
WO (1) WO2023005061A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023005061A1 (zh) * 2021-07-29 2023-02-02 律音科技(武汉)有限公司 固体材料预紧力和温度的超声波双波测量方法、用途及设备

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116358747B (zh) * 2023-04-23 2023-11-03 哈尔滨工业大学 一种基于双晶空耦换能器的高空间分辨率平面应力测量方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3027157B1 (ja) * 1999-06-07 2000-03-27 株式会社東日製作所 センサ
JP2001174343A (ja) * 1999-12-17 2001-06-29 Tadashi Onishi ボルト軸力の測定方法及び装置
JP2002277336A (ja) * 2001-03-15 2002-09-25 Tokyo Gas Co Ltd 低温下におけるボルトの軸力測定方法とその装置
JP2009031180A (ja) * 2007-07-30 2009-02-12 Central Res Inst Of Electric Power Ind 内部温度の測定方法および測定装置
CN102636307A (zh) * 2012-04-10 2012-08-15 交通运输部公路科学研究所 一种竖向预应力钢筋有效应力水平测试装置及其测试方法
CN108387338A (zh) * 2018-02-07 2018-08-10 大连理工大学 一种基于压电超声晶片的螺栓预紧力实时高精度检测方法及系统
CN110530571A (zh) * 2019-08-30 2019-12-03 航天精工股份有限公司 一种对带传感器外螺纹紧固件进行预紧力标定的方法
CN111537132A (zh) * 2020-04-20 2020-08-14 中物院成都科学技术发展中心 一种轴向预紧力双波测量方法
CN112577653A (zh) * 2020-12-11 2021-03-30 中铁大桥局集团有限公司 一种桥梁高强度螺栓紧固轴力的测量方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019153138A1 (zh) * 2018-02-07 2019-08-15 大连理工大学 一种基于压电超声晶片的螺栓预紧力实时高精度检测方法及系统
CN111693190A (zh) 2020-08-06 2020-09-22 浙江中自庆安新能源技术有限公司 一种基于超声波的螺栓轴向应力测量装置及方法
CN113587992B (zh) * 2021-07-29 2022-07-12 律音科技(武汉)有限公司 固体材料预紧力和温度的超声波双波测量方法、用途及设备

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3027157B1 (ja) * 1999-06-07 2000-03-27 株式会社東日製作所 センサ
JP2001174343A (ja) * 1999-12-17 2001-06-29 Tadashi Onishi ボルト軸力の測定方法及び装置
JP2002277336A (ja) * 2001-03-15 2002-09-25 Tokyo Gas Co Ltd 低温下におけるボルトの軸力測定方法とその装置
JP2009031180A (ja) * 2007-07-30 2009-02-12 Central Res Inst Of Electric Power Ind 内部温度の測定方法および測定装置
CN102636307A (zh) * 2012-04-10 2012-08-15 交通运输部公路科学研究所 一种竖向预应力钢筋有效应力水平测试装置及其测试方法
CN108387338A (zh) * 2018-02-07 2018-08-10 大连理工大学 一种基于压电超声晶片的螺栓预紧力实时高精度检测方法及系统
CN110530571A (zh) * 2019-08-30 2019-12-03 航天精工股份有限公司 一种对带传感器外螺纹紧固件进行预紧力标定的方法
CN111537132A (zh) * 2020-04-20 2020-08-14 中物院成都科学技术发展中心 一种轴向预紧力双波测量方法
CN112577653A (zh) * 2020-12-11 2021-03-30 中铁大桥局集团有限公司 一种桥梁高强度螺栓紧固轴力的测量方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023005061A1 (zh) * 2021-07-29 2023-02-02 律音科技(武汉)有限公司 固体材料预紧力和温度的超声波双波测量方法、用途及设备

Also Published As

Publication number Publication date
WO2023005061A1 (zh) 2023-02-02
EP4365557A1 (en) 2024-05-08
CN113587992B (zh) 2022-07-12

Similar Documents

Publication Publication Date Title
CN113587992B (zh) 固体材料预紧力和温度的超声波双波测量方法、用途及设备
US11131579B2 (en) Piezoelectric patch-based real-time and high-precision bolt preload detection method and system
EP3121586B1 (en) Fatigue tester
Zhang et al. Rapid determination of fatigue life based on temperature evolution
Schubel et al. Review of structural health and cure monitoring techniques for large wind turbine blades
US20090048788A1 (en) Rapid determination of fatigue failure based on temperature evolution
Meyendorf et al. Thermo-acoustic fatigue characterization
Gatowski et al. An experimental investigation of surface thermometry and heat flux
CN102840968B (zh) 一种航空发动机叶片大量程振动幅值检测装置及方法
CN111693190A (zh) 一种基于超声波的螺栓轴向应力测量装置及方法
CN104792435A (zh) 基于瞬态热边界反演的结构内部非均匀温度场的重建方法
CN112748007B (zh) 一种基于声发射的抗腐蚀疲劳性能测试装置
WO2022257409A1 (zh) 一种利用超声测量螺栓载荷的方法
CN111795765A (zh) 一种预紧件服役应力的无损监测装置和方法
Rokhlin et al. In situ ultrasonic measurement of crack closure
CN209946078U (zh) 一种超声波检测探头的标定装置
Beard et al. Direct shaft torque measurements in a transient turbine facility
CN212645940U (zh) 一种基于超声波的螺栓轴向应力测量装置
Tamura et al. Non-contact vibration measurement of the rotor blades that play a pivotal role in the reliability of gas turbines
CN105738407B (zh) 锚杯热膨胀系数检测方法和检测装置
Bin et al. Effect of temperature on solid ultrasonic propagation using finite element method and experiments
Li et al. Bolt axial stress detection and calibration by ultrasound shear and longitudinal wave
RU2261437C1 (ru) Способ теплового неразрушающего контроля многослойных объектов
Mohammed et al. Thermal Product of Fast Response Temperature Sensors for Transient Heat Transfer Applications with Numerically Determined Surface Heat Flux History
CN113984256B (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