CN107063145B - 超声测厚中的入射偏角自动辨识与误差补偿方法 - Google Patents

超声测厚中的入射偏角自动辨识与误差补偿方法 Download PDF

Info

Publication number
CN107063145B
CN107063145B CN201710047840.7A CN201710047840A CN107063145B CN 107063145 B CN107063145 B CN 107063145B CN 201710047840 A CN201710047840 A CN 201710047840A CN 107063145 B CN107063145 B CN 107063145B
Authority
CN
China
Prior art keywords
drift angle
thickness
incident
ultrasonic
calibration
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
CN201710047840.7A
Other languages
English (en)
Other versions
CN107063145A (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.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201710047840.7A priority Critical patent/CN107063145B/zh
Publication of CN107063145A publication Critical patent/CN107063145A/zh
Application granted granted Critical
Publication of CN107063145B publication Critical patent/CN107063145B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B17/00Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations
    • G01B17/02Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations for measuring thickness

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)

Abstract

本发明超声测厚中的入射偏角自动辨识与误差补偿方法属于超声检测技术领域,特别涉及一种超声在机测厚中的入射偏角自动辨识与误差补偿方法。该方法中,首先,进行超声入射偏角标定,提取不同标定入射偏角下的超声信号首次回波能量,建立标定入射偏角与首次回波能量间的关联关系。然后,采样获得被测件的实际超声测厚信号,利用所建立的标定入射偏角与首次回波能量间的关联关系自动辨识出入射偏角。最后,基于自相关分析方法计算当前测点回波声时差,利用建立的厚度误差补偿模型,完成入射偏角关联厚度误差补偿,获得精确测厚结果。本发明实现了超声在机测量中的入射偏角自动辨识及其关联厚度误差补偿,算法简单、可靠度高、实时性强。

Description

超声测厚中的入射偏角自动辨识与误差补偿方法
技术领域
本发明属于超声检测技术领域,特别涉及一种超声在机测厚中的入射偏角自动辨识与误差补偿方法。
背景技术
超声法因具有界面反射和材料透射特性,往往是生产现场中零部件厚度检测首选方法。将超声传感器集成于加工装备,实现在机自动测厚,可有效解决大型薄壁零件高效测厚难题。
目前,超声测厚主要采用脉冲反射原理,其要求传感器声束中心线与被测点处法矢重合,否则存在入射偏角。实际测量中,即使入射偏角在超声传感器许用域内,由于超声波在工件与耦合剂界面处的反射、折射状态发生变化,会致使超声传播路径改变,进而引起较大的厚度测量误差。然而,超声在机测厚中,由于被测面形复杂、跟踪响应时滞等因素影响,传感器声束中心线与测点法矢无法实时保持重合,不可避免地存在入射偏角。为此,如何在超声在机测量过程中快速、精确、自动辨识出超声入射偏角,补偿其引起的厚度测量误差、对于提高超声测厚精度,至关重要。
2014年,张晓琳等在发明专利CN201410810870.5中公开了一种大型壁板厚度在位测量方法。将超声波测厚仪的探头通过装卡机构安装在数控铣床的铣刀位置,机床控制超声波测厚仪探头运动,进行逐点接触测厚,通过数据传输处理模块分别将超声波测厚仪的存储器中的测量数据反馈给终端处理设备进行数据处理,完成壁板厚度的在位测量。2015年,康宜华等人在发明专利CN201510165390.2中公开了一种钢管超声测厚装置,包括机架、支撑座、钢管跟踪组件和压管限位组件等。支撑座安装于机架上,支撑座滑块与水槽固定连接带动水槽和超声测头装置靠近或远离压管限位组件,压管限位组件与水槽配合夹住钢管从而对钢管限位,钢管进入检测位置时水槽带动超声测头装置进入检测工位,实现钢管高速移动时的高速检测。
然而,上述超声自动测厚方法均未涉及超声入射偏角自动辨识及误差补偿方法。
发明内容
本发明主要解决的技术问题是克服现有方法的不足,针对超声在机自动测厚中,超声传感器声束中心线与测点法矢无法保持实时重合,导致入射偏角引起的厚度测量误差问题,发明了一种超声在机测厚中的入射偏角自动辨识与误差补偿方法。该方法基于超声波在被测材料表面的传播规律,建立厚度误差补偿模型,实现零部件厚度的超声准确测量,保证超声在机测量结果一致性。方法基于标定出的入射偏角-回波能量衰减的关联关系,通过计算实际超声测厚信号中的首次反射回波能量,自动辨识超声波入射偏角,方法可靠性高。
本发明采用的技术方案是一种超声测厚中的入射偏角自动辨识与误差补偿方法,其特征是:首先,进行超声入射偏角标定,提取不同标定入射偏角下的超声信号首次回波能量,建立标定入射偏角与首次回波能量间的关联关系;然后,采样获得被测件的实际超声测厚信号,利用所建立的标定入射偏角与首次回波能量间的关联关系自动辨识出入射偏角;最后,基于自相关分析方法计算当前测点回波声时差,利用建立的厚度误差补偿模型,完成入射偏角关联厚度误差补偿,获得精确测厚结果。方法的具体步骤如下:
第一步超声入射偏角在机标定
首先进行标定准备,将超声传感器1通过夹具5安装在数控机床主轴4上,转动平台8通过平台螺栓压板组件9固定在机床工作台10上,水槽6由水槽螺栓压板组件12固定放置在转动平台8上,标定件7放置于水槽6中央底部,然后向水槽6注入耦合剂2。主轴4带动超声传感器1,使其浸没于耦合剂2中,且超声传感器1与标定件7间保持一定的测量偏置距离。
按照固定转角步长增量调整转动平台8,在每个增量角度下均进行超声信号采集。设定转动平台8水平位置采集获得的超声回波信号作为基准信号s0。首先,对获得超声回波信号序列s(k),k∈[1N]进行小波分解,
式中,dj(k)为第j层小波分量,a(k)为剩余分量,r为分解层数,k为序列号。
根据传感器频率响应,用第p~q层小波分量表征滤波后超声回波信号g(k),
计算滤波后超声回波信号g(k)的首次回波能量E(θi,z0,Rq),
其中,θi为标定入射偏角,Rq为标定件表面粗糙度,m1为首次回波起始序号位置,m2为首次回波终止序号位置。
相对于基准信号s0,即入射无偏角时的超声回波信号,标定入射偏角θi下首次回波信号能量衰减为Ei,dBi,z0,Rq),
基于标定过程中获得的超声测厚信号入射偏角-首次回波能量衰减数据点(θi,Ei,dBi,z0,Rq)),i=1,2,…,M,共M组,采用多项式拟合法建立标定入射偏角与首次回波能量间的关联关系曲线。
第二步入射偏角自动辨识
将被测件3放置在水槽6并浸没于耦合剂2中,调节机床运动至被测件3待测量点上方,驱动主轴4向下运动,使得超声传感器1与被测件3间具有一定的测量偏置距离,采集获得当前测点的超声回波信号序列y(k)。通过公式(1)至公式(4),计算滤波后当前测点超声回波信号gy(k)以及首次回波能量衰减EdBL,z0,Rq)。根据入射偏角与首次回波能量间的关联关系曲线,自动辨识出被测件3当前测点入射偏角
第三步厚度误差补偿
对被测件3的当前测点处超声回波信号进行自相关分析,并利用回波峰值计算回波时间差。滤波后当前测点处超声回波信号gy(k)自相关函数r(n),
利用自相关函数r(n)曲线中最大峰值对应的序列号n,计算当前测点回波声时差△t,
△t=arg maxpeaks(r(n))/fs (6)
其中,fs为采集频率。
按如下公式计算被测件3当前测点的的厚度测量值d',
其中,cL2为超声波在被测件3中的传播速度。
考虑超声波斜入射时的界面折射,被测件3当前测点的厚度实际值d与厚度测量值d'存在如下关系,
d'=d/cosβL (8)
其中,βL为折射角度。
利用超声波在介质中的传播定律,折射角度βL满足如下关系,
其中,cL1为超声波在耦合剂2中的传播速度。
则,被测件3当前测点的厚度实际值d为,
将自动辨识出的当前测点入射偏角和回波声时差△t带入公式(9)和公式(10),计算获得补偿后的当前测点厚度。
本发明的有益效果是:在机标定了超声波入射偏角与信号首次回波能量衰减间的关联关系,根据标定曲线自动辨识当前测点入射偏角,运算速度快、反算精度高;入射偏角关联厚度误差补偿模型可实现不同入射角度下厚度的精确补偿,方法可靠度高,实时性强。
附图说明
附图1-入射偏角在机标定示意图,附图2-超声厚度测量原理图,其中:1-超声传感器,2-耦合剂,3-被测件,4-主轴,5-夹具,6-水槽,7-标定件,8-转动平台,9-平台螺栓压板组,10-机床工作台,11-信号线,12-水槽螺栓压板组,-当前测点入射偏角,βL-折射角,d'-厚度测量值,d-厚度实际值,P-当前测点位置。
附图3-标定入射角度1.8°时的超声回波信号波形。其中,横坐标-序号K,纵坐标-超声回波信号序列s(k)。
附图4-滤波后超声回波信号波形,其中:m1-首次回波起始序号位置,m2-首次回波终止序号位置。横坐标-序号K,纵坐标-超声回波信号g(k)。
附图5-标定入射偏角与首次回波能量间的关联关系曲线,其中:L1-标定数据,L2-拟合曲线,横坐标-标定入射角θi 0,纵坐标-回波能量衰减Ei,dB
具体实施方式
结合附图和技术方案详细说明本发明的实施方式,说明超声在机测厚中的入射偏角自动辨识与误差补偿方法。
测量参数:标定件7,材料45号钢、标准长方体量块、设计尺寸为30mm×30mm×5mm;被测件3,材料45号钢、弧形件、设计厚度8mm、半径48mm;水浸式超声平探头,脉冲宽度10MHz、晶片直径6.35mm、采集频率1.25GHz。
第一步,超声入射偏角标定准备。入射偏角标定示意图如图1所示。超声传感器1通过夹具5安装在主轴4上,转动平台8通过平台螺栓压板组9固定在机床工作台10上,水槽6由水槽螺拴压板组12固定在转动平台8上,标定件7放置在水槽6中央底部,水槽6注满耦合剂2。
第二步,在机标定。主轴4拖动超声传感器1浸没于耦合剂2中,且超声传感器1与标定件7间的测量偏置距离10mm。调整转动平台8至水平位置,则当前位置下的超声测量回波信号为基准信号s0。然后,控制转动平台8每次顺时针增量转动,转角步长0.2°,采集超声回波数据,直至转角达到2.8°,即传感器最大许用角度。标定入射角度1.8°时的超声回波信号波形如图3所示。对图3中的超声回波信心进行小波处理,经小波分解7次后采用第6、7层高频分量表征滤波后超声回波信号g(k)如图4所示。针对不同标定入射偏角下的超声回波信号,计算其首次回波能量衰减,进一步利用多项式拟合建立θ-EdB(θ,z0,Rq)标定关系曲线,如图5所示,其中,L1-标定数据,L2-拟合曲线。
第三步,入射偏角自动辨识。取下标定件7,并将被测件3放置在水槽6中央底部。机床控制超声传感器1至当前测点P,采集超声测量回波信号,通过公式(1)~(4)计算得到EdBL,z0,Rq),然后利用图5所示的标定曲线,辨识出当前测点P的入射偏角为1.8°。
第四步,入射偏角关联厚度误差补偿。通过对接收到的回波信号进行自相关运算,取相关函数峰值的出现时刻计算为回波信号声时差,当前测点P处的声时差△t=2.715μs。将△t=2.715μs、被测件中超声传播速度cL2=5950m/s代入到公式(7),计算当前测点P处的测量厚度d'=8.077mm。将标定入射偏角耦合剂中的超声传播速度cL1=1493.5m/s、被测件中超声传播速度cL2=5950m/s,代入公式(9),得到折射角度βL=7.75°。最后将d'=8.077mm与βL=7.75°代入公式(8),得到当前测点P补偿后的厚度实际值d=8.003mm。
本发明实现了超声在机测量入射偏角自动辨识与误差补偿,补偿结果准确,操作简单、方法可靠。

Claims (1)

1.一种超声测厚中的入射偏角自动辨识与误差补偿方法,其特征是:首先,进行超声入射偏角标定,提取不同标定入射偏角下的超声信号首次回波能量,建立标定入射偏角与首次回波能量间的关联关系;然后,采样获得被测件的实际超声测厚信号,利用所建立的标定入射偏角与首次回波能量间的关联关系自动辨识出入射偏角;最后,基于自相关分析方法计算当前测点回波声时差,利用建立的厚度误差补偿模型,完成入射偏角关联厚度误差补偿,获得精确测厚结果;方法具体步骤如下:
第一步超声入射偏角在机标定
首先进行标定准备;将超声传感器(1)通过夹具(5)安装在数控机床主轴(4)上,转动平台(8)通过平台螺栓压板组件(9)固定在机床工作台(10)上,水槽(6)由水槽螺栓压板组件(12)固定放置在转动平台(8)上,标定件(7)放置于水槽(6)中央底部,然后向水槽(6)注入耦合剂(2);主轴(4)拖动超声传感器(1),使其浸没于耦合剂(2)中,且超声传感器(1)与标定件(7)间保持一定的测量偏置距离;
按照固定转角步长增量调整转动平台(8),在每个增量角度下均进行超声信号采集;设定转动平台(8)水平位置采集获得的超声回波信号作为基准信号s0;首先,对获得超声回波信号序列s(k),k∈[1N]进行小波分解,
式中,dj(k)为第j层小波分量,a(k)为剩余分量,r为分解层数,k为序列号;
根据传感器频率响应,用第p~q层小波分量表征滤波后超声回波信号g(k),
计算滤波后超声回波信号g(k)的首次回波能量E(θi,z0,Rq),
其中,θi为标定入射偏角,Rq为标定件表面粗糙度,m1为首次回波起始序号位置,m2为首次回波终止序号位置;
相对于基准信号s0,即入射无偏角时的超声回波信号,标定入射偏角θi下首次回波信号能量衰减为Ei,dBi,z0,Rq),
基于标定过程中获得的入射偏角-首次回波能量衰减超声测厚信号数据点(θi,Ei,dBi,z0,Rq)),i=1,2,…,M,共M组,采用多项式拟合法建立标定入射偏角与首次回波能量间的关联关系曲线;
第二步入射偏角自动辨识
将被测件(3)放置在水槽(6)并浸没于耦合剂(2)中,调节机床运动至被测件(3)待测量点上方,驱动主轴(4)向下运动,使得超声传感器(1)与被测件(3)间具有一定的测量偏置距离,采集获得当前测点的超声回波信号序列y(k);通过公式(1)至公式(4),计算滤波后当前测点超声回波信号gy(k)以及首次回波能量衰减EdBL,z0,Rq);根据入射偏角与首次回波能量间的关联关系曲线,自动辨识出被测件(3)当前测点入射偏角
第三步厚度误差补偿
对被测件(3)的当前测点处超声回波信号进行自相关分析,并利用回波峰值计算回波时间差;滤波后当前测点处超声回波信号gy(k)自相关函数r(n),
利用自相关函数r(n)曲线中最大峰值对应的序列号n,计算当前测点回波声时差Δt,
Δt=argmaxpeaks(r(n))/fs (6)
其中,fs为采集频率;
按如下公式计算被测件(3)当前测点的的厚度测量值d',
其中,cL2为超声波在被测件(3)中的传播速度;
考虑超声波斜入射时的界面折射,被测件(3)当前测点的厚度实际值d与厚度测量值d'存在如下关系,
d'=d/cosβL (8)
其中,βL为折射角度;
利用超声波在介质中的传播定律,折射角度βL满足如下关系:
其中,cL1为超声波在耦合剂(2)中的传播速度;
被测件(3)当前测点的厚度实际值d为:
将自动辨识出的当前测点入射偏角和回波声时差Δt带入公式(9)和公式(10),计算获得补偿后的当前测点厚度。
CN201710047840.7A 2017-01-26 2017-01-26 超声测厚中的入射偏角自动辨识与误差补偿方法 Active CN107063145B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710047840.7A CN107063145B (zh) 2017-01-26 2017-01-26 超声测厚中的入射偏角自动辨识与误差补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710047840.7A CN107063145B (zh) 2017-01-26 2017-01-26 超声测厚中的入射偏角自动辨识与误差补偿方法

Publications (2)

Publication Number Publication Date
CN107063145A CN107063145A (zh) 2017-08-18
CN107063145B true CN107063145B (zh) 2018-07-24

Family

ID=59598382

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710047840.7A Active CN107063145B (zh) 2017-01-26 2017-01-26 超声测厚中的入射偏角自动辨识与误差补偿方法

Country Status (1)

Country Link
CN (1) CN107063145B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP4160142A1 (en) * 2021-09-29 2023-04-05 Renishaw PLC Ultrasound method for inspecting a part

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107702668A (zh) * 2017-11-16 2018-02-16 哈尔滨工程大学 一种超声波测薄油膜用标定测试台
CN108692684A (zh) * 2018-05-30 2018-10-23 昆山睿力得软件技术有限公司 一种自动探测刹车片厚度的系统及方法
CN109059788B (zh) * 2018-10-26 2021-04-06 中国特种设备检测研究院 厚度测量方法及设备
CN112731405B (zh) * 2021-01-04 2023-12-29 中国航空制造技术研究院 一种飞机壁板的在线超声测量方法及加工方法
CN113063593B (zh) * 2021-03-15 2022-02-22 西安交通大学 基于超声回波能量系数的滚动轴承信号特征信息获取方法
CN114858059B (zh) * 2022-07-06 2022-09-20 北京航空航天大学 降低沿视线积分测量位置系统误差的装置、方法和系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3425675B2 (ja) * 1993-09-09 2003-07-14 エフニック株式会社 真珠の巻きの厚み測定方法
CN101710102B (zh) * 2009-12-16 2011-02-16 内蒙古北方重工业集团有限公司 变壁、厚壁管材超声波检测方法
CN103033153B (zh) * 2012-12-11 2015-05-20 浙江大学 扫描超声波显微镜同时测量薄层材料机械特性参数的方法
CN104764423B (zh) * 2015-04-09 2017-06-30 武汉华宇一目检测装备有限公司 一种钢管超声测厚装置
CN105021142B (zh) * 2015-07-15 2017-09-29 中国科学院金属研究所 一种激光搭接焊缝宽度的测量方法和所用装置

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP4160142A1 (en) * 2021-09-29 2023-04-05 Renishaw PLC Ultrasound method for inspecting a part
WO2023052756A1 (en) * 2021-09-29 2023-04-06 Renishaw Plc Ultrasound method for inspecting a part

Also Published As

Publication number Publication date
CN107063145A (zh) 2017-08-18

Similar Documents

Publication Publication Date Title
CN107063145B (zh) 超声测厚中的入射偏角自动辨识与误差补偿方法
KR101804484B1 (ko) 센서 장치 및 상기 장치를 사용한 잔류 응력 검측 시스템
CN104076089B (zh) 环形锻件自动超声波c扫描检测系统
US20150177194A1 (en) Dual Robot Detection Apparatus For Non-Damage Detection
CN102494645B (zh) 一种基于超声的内腔尺寸精密测量装置及其方法
CN203275373U (zh) 一种非金属超声检测仪校准装置
CN210621786U (zh) 一种基桩超声波透射法导向轮结构改进
CN109540043A (zh) 一种使用激光测距原理测量板坯锥度仪锥度的方法
CN109521092A (zh) 一种表面波和板波声速的非接触式测量装置及其方法
JP2007212410A (ja) 膜厚較正曲線の取得方法
CN205749405U (zh) 水浸法超声波检测探头角度调节尺
CN103033566A (zh) 一种超声波探头扩散角自动测定装置
CN103822968B (zh) 面向结合面压强检测的压强-超声反射率曲线构建方法
CN116818057B (zh) 一种流量计现场计量系统及方法
CN108594212A (zh) 一种全方位测距装置
CN107941179B (zh) 一种超声波测气缸活塞环油膜厚度测试装置
CN1011544B (zh) 孔口测绘及表面时间测量系统
CN111693611A (zh) 一种利用激光超声检测金属亚表面缺陷的方法及系统
CN216246171U (zh) 一种基于超声波的涂层厚度测量装置
CN110360946A (zh) 一种测量物体深度的装置及方法
CN111982034B (zh) 一种用于测量液体静压导轨油膜厚度的超声波检测装置
CN112484836B (zh) 一种超声探头装置及工件声速测量方法
CN209471081U (zh) 一种表面波和板波声速的非接触式测量装置
CN104568079B (zh) 超声波液位计的模拟测试方法
CN208255408U (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