CN103822968A - 面向结合面压强检测的压强-超声反射率曲线构建方法 - Google Patents

面向结合面压强检测的压强-超声反射率曲线构建方法 Download PDF

Info

Publication number
CN103822968A
CN103822968A CN201310551984.8A CN201310551984A CN103822968A CN 103822968 A CN103822968 A CN 103822968A CN 201310551984 A CN201310551984 A CN 201310551984A CN 103822968 A CN103822968 A CN 103822968A
Authority
CN
China
Prior art keywords
pressure
faying face
reflectivity
pressure intensity
average
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
CN201310551984.8A
Other languages
English (en)
Other versions
CN103822968B (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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201310551984.8A priority Critical patent/CN103822968B/zh
Publication of CN103822968A publication Critical patent/CN103822968A/zh
Application granted granted Critical
Publication of CN103822968B publication Critical patent/CN103822968B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明公开了一种面向结合面压强检测的压强-超声反射率曲线构建方法,其首计算不同压强下每个扫描点的超声反射率;然后拟合求得初始压强-反射率关系;根据求得的初始压强-反射率关系,计算在不同压力下结合面所有扫描点压强,得到整个结合面的总的压力F:将F与利用压力传感器测量到的实际压力F0计算标定系数;然后重新计算每个扫描点的压强值Pi得到修正后的压强值Pi’,计算不同压力下结合面中心压强分布平坦区域扫描点的压强Pi’的均值Pα,利用所得的Pα与该区域对应的反射率均值R0重新拟合求得新的压强-反射率曲线。本发明采用迭代方式构建压强-反射率曲线,能够构建准确的压强-反射率关系,实现结合面接触压强的准确检测。

Description

面向结合面压强检测的压强-超声反射率曲线构建方法
技术领域
本发明属于超声检测技术领域,涉及一种压强-反射率曲线的构建,尤其是一种面向结合面压强检测的压强-超声反射率曲线构建方法。 
背景技术
机械系统中大量存在的结合面对系统性能有重要影响。随着高端装备制造、装配环节的日益精量化,实现对结合面接触压强分布的准确检测尤为迫切。而由于结合面的封闭性,给实验测量带来诸多困难。目前结合面接触压强分布的测量方法,主要有薄膜法和超声法。其中薄膜法是在结合面之间置入具有一定厚度的薄膜传感器,利用该传感器检测接触压强分布。但是薄膜法改变了工件的接触状况,进而会影响检测精度,且对于已经装配好的工件,薄膜法也不再适用。而超声检测能够实现对接触压强分布的直接无损检测。现有的超声检测方法,要首先制作与被测对象相同材料、表面形貌的试件,通过实验获取压强-反射率曲线,然后检测被测结合面的超声波反射率(即超声脉冲反射波与入射波的比值),将反射率带入上述曲线中即可得到结合面的压强分布。但是由于采用目前方法构建的压强-反射率曲线不够准确,导致超声法检测结果与理论模型计算结果差距较大,检测精度较低。 
发明内容
本发明的目的在于克服上述现有技术的缺点,提供一种面向结合面压强检测的压强-超声反射率曲线构建方法。 
本发明的目的是通过以下技术方案来实现的: 
这种面向结合面压强检测的压强-超声反射率曲线构建方法,包括以下步骤: 
1)利用聚焦超声探头扫描结合面,计算不同压强下每个扫描点的超声反射率;然后利用结合面中心压强分布平坦区域扫描点的反射率均值R0及对应的压强P0,拟合求得初始压强-反射率关系,如式(1)所示: 
P0=f(R0)    (1) 
2)根据式(1)求得的初始压强-反射率关系,计算在不同压力下结合面所有扫描点压强Pi,然后,将Pi带入下面式(2),得到整个结合面的总的压力F: 
F = Σ i = 1 n P i × ΔA i - - - ( 2 )
其中,ΔAi为单个扫描点代表的面积,即整个结合面的面积除以扫描点数;F为计算得到的结合面受到的压力;将上述F与利用压力传感器测量到的实际压力F0带入下面式(3): 
k = F F 0 - - - ( 3 )
其中,k为标定系数;每个扫描点的压强值Pi利用下面式(4)重新计算 
Pi'=kPi    (4) 
其中,Pi’为每个扫描点修正后的压强值; 
3)计算不同压力下结合面中心压强分布平坦区域扫描点的压强Pi’的均值Pα,利用所得的Pα与该区域对应的反射率均值R0重新拟合求得 新的压强-反射率曲线。 
进一步,以上步骤1)中,反射率均值R0及对应的压强P0的获取方法具体如下: 
(1)为检测螺栓结合面接触压强分布,按照与螺栓结合面相同的材料、表面粗糙度加工两个试件:上试件和下试件,将两试件安装到压力加载试验台上; 
(2)获得参考信号:在结合面试件未接触的时候,利用奥林巴斯V319型15MHz水浸式聚焦超声波探头对上试件进行扫描,利用5800超声波脉冲发射接收器发射并接受超声脉冲信号,利用泰克MSO4104示波器进行接收并保存到电脑上,得到结合面上各个扫描点上的超声回波信号,该组信号为参考信号; 
(3)保持上试件不动,让下试件与上试件接触,并施加不同的压力,其中压力的大小可以由压力传感器测得;按照第一次扫描的轨迹对结合面进行重新扫描,得到不同压力下结合面上各个扫描点的超声回波信号,该回波信号为结合面信号; 
(4)分别对参考信号与结合面信号进行傅里叶变换,并计算每个扫描点在不同压力下的反射率;结合面所有点的反射率组成反射率矩阵,得到反射率均值R0; 
(5)计算根据压力计算结合面的平均压强,平均压强的80%即为中心区域压强,即反射率均值R0对应的压强P0。 
本发明具有以下有益效果: 
本发明的面向结合面压强检测的压强-超声反射率曲线构建方法采用迭代方式构建压强-反射率曲线,与现有的方法相比,本发明能够构建准确的压强-反射率关系,实现结合面接触压强的准确检测;使超声法检测结果与理论模型计算结果差距小,检测精度高。 
附图说明
图1是本发明结合面试件压力加载试验台; 
图2是结合面试件的反射率矩阵; 
图3是单钉螺栓结合面在15000N预紧力下结合面中心线上压强分布检测结果与理论计算结果的对比。 
具体实施方式
下面结合附图对本发明做进一步详细描述: 
本发明面向结合面压强检测的压强-超声反射率曲线构建方法按照以下步骤实施: 
1)为检测螺栓结合面接触压强分布,按照与螺栓结合面相同的材料、表面粗糙度加工两个试件,将试件安装到压力加载试验台上,见图1,在图1中,1为超声探头,2为水,3为上试件,4为被测结合面,5为下试件,6为轴承,7为压力传感器,8为液压千斤顶。 
2)获得参考信号:在结合面试件未接触的时候,利用奥林巴斯V319型15MHz水浸式聚焦超声波探头对上试件(参见图1)进行扫描,利用5800超声波脉冲发射接收器发射并接受超声脉冲信号。利用泰克MSO4104示波器进行接收并保存到电脑上,得到结合面上各个扫描点上的超声回波信号,该组信号为参考信号。 
3)保持上试件不动,让下试件与上试件接触,并施加不同的压力,其中压力的大小可以由压力传感器测得。按照第一次扫描的轨迹对结合面进行重新扫描,得到不同压力下结合面上各个扫描点的超声回波信号,该回波信号为结合面信号。 
4)分别对参考信号与结合面信号进行傅里叶变换,并利用下面式子计算每个扫描点在不同压力下的反射率。 
R = H i H o
其中,R为超声反射率幅值,Hi为结合面信号幅值,Ho为相同频率下参考信号幅值;结合面所有点的反射率组成反射率矩阵,参见图2。图2的反射率矩阵中,中间圆环区域压强分布均匀,因此计算该区域的平均反射率R0。 
5)计算根据压力计算结合面的平均压强,平均压强的80%即为中心区域压强。利用中心区域压强值P0及该区域反射率均值R0,拟合求得初始压强-反射率关系,如式(1)所示: 
P0=f(R0)    (1) 
6)根据初始压强-反射率关系及整个结合面的反射率矩阵,计算得整个结合面的压力,如式(2)所示: 
F = Σ i = 1 n P i × ΔA i - - - ( 2 )
其中,ΔAi为单个扫描点代表的面积,即整个结合面的面积除以扫描点数。F为计算得到的结合面受到的压力。而此时整个面的实际压力F0可以由压力传感器测得,利用下面式(3)计算不同压力下的标定系数: 
k = F F 0 - - - ( 3 )
其中,k为标定系数。每个扫描点的压强值Pi利用下面式(4)重新计算 
Pi'=kPi    (4) 
其中,Pi’为每个扫描点修正后的压强值。 
7)利用修正后每个扫描点的压强值Pi’重新计算中心平坦区域的平均压强值Pα,将所得的Pα与之前超声测量的中心区域反射率R0重新拟合求得压强-反射率曲线。 
8)对单钉螺栓结合面施加15000N预紧力,并测量此时的螺栓结合面超声波反射率。利用所求的压强-反射率曲线计算结合面的压强分布。螺栓结合面中心线上压力分布检测值与理论计算结果的对比,结果参见图3。 

Claims (2)

1.一种面向结合面压强检测的压强-超声反射率曲线构建方法,其特征在于,包括以下步骤:
1)利用聚焦超声探头扫描结合面,计算不同压强下每个扫描点的超声反射率;然后利用结合面中心压强分布平坦区域扫描点的反射率均值R0及对应的压强P0,拟合求得初始压强-反射率关系,如式(1)所示:
P0=f(R0)  (1)
2)根据式(1)求得的初始压强-反射率关系,计算在不同压力下结合面所有扫描点压强Pi,然后,将Pi带入下面式(2),得到整个结合面的总的压力F:
F = Σ i = 1 n P i × Δ A i - - - ( 2 )
其中,ΔAi为单个扫描点代表的面积,即整个结合面的面积除以扫描点数;F为计算得到的结合面受到的压力;将上述F与利用压力传感器测量到的实际压力F0带入下面式(3):
k = F F 0 - - - ( 3 )
其中,k为标定系数;每个扫描点的压强值Pi利用下面式(4)重新计算
Pi'=kPi  (4)
其中,Pi’为每个扫描点修正后的压强值;
3)计算不同压力下结合面中心压强分布平坦区域扫描点的压强Pi’的均值Pα,利用所得的Pα与该区域对应的反射率均值R0重新拟合求得新的压强-反射率曲线。
2.根据权利要求1所述的面向结合面压强检测的压强-超声反射率曲线构建方法,其特征在于,步骤1)中,反射率均值R0及对应的压强P0的获取方法具体如下:
(1)为检测螺栓结合面接触压强分布,按照与螺栓结合面相同的材料、表面粗糙度加工两个试件:上试件和下试件,将两试件安装到压力加载试验台上;
(2)获得参考信号:在结合面试件未接触的时候,利用奥林巴斯V319型15MHz水浸式聚焦超声波探头对上试件进行扫描,利用5800超声波脉冲发射接收器发射并接受超声脉冲信号,利用泰克MSO4104示波器进行接收并保存到电脑上,得到结合面上各个扫描点上的超声回波信号,该组信号为参考信号;
(3)保持上试件不动,让下试件与上试件接触,并施加不同的压力,其中压力的大小可以由压力传感器测得;按照第一次扫描的轨迹对结合面进行重新扫描,得到不同压力下结合面上各个扫描点的超声回波信号,该回波信号为结合面信号;
(4)分别对参考信号与结合面信号进行傅里叶变换,并计算每个扫描点在不同压力下的反射率;结合面所有点的反射率组成反射率矩阵,得到反射率均值R0
(5)根据压力计算结合面的平均压强,平均压强的80%即为中心区域压强,即反射率均值R0对应的压强P0
CN201310551984.8A 2013-11-08 2013-11-08 面向结合面压强检测的压强-超声反射率曲线构建方法 Active CN103822968B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310551984.8A CN103822968B (zh) 2013-11-08 2013-11-08 面向结合面压强检测的压强-超声反射率曲线构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310551984.8A CN103822968B (zh) 2013-11-08 2013-11-08 面向结合面压强检测的压强-超声反射率曲线构建方法

Publications (2)

Publication Number Publication Date
CN103822968A true CN103822968A (zh) 2014-05-28
CN103822968B CN103822968B (zh) 2016-08-17

Family

ID=50758132

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310551984.8A Active CN103822968B (zh) 2013-11-08 2013-11-08 面向结合面压强检测的压强-超声反射率曲线构建方法

Country Status (1)

Country Link
CN (1) CN103822968B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105277149A (zh) * 2015-11-18 2016-01-27 西安建筑科技大学 结合面真实接触面积测量装置与测量方法
CN106932277A (zh) * 2017-04-28 2017-07-07 大连理工大学 基于圆角平面接触理论的界面超声波反射率‑压强关系曲线建立方法及加载试验台
WO2018196148A1 (zh) * 2017-04-28 2018-11-01 大连理工大学 基于圆角平面接触理论的界面超声波反射率-压强关系曲线建立方法及加载试验台
CN108908378A (zh) * 2018-05-14 2018-11-30 莱茵检测认证服务(中国)有限公司 一种人与机器人接触过程中接触压强的测试方法
CN108981989A (zh) * 2018-10-10 2018-12-11 郑州大学 一种检测轴套过盈配合结合面压强的装置及其使用方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0342567A (ja) * 1989-07-10 1991-02-22 Hitachi Constr Mach Co Ltd 非破壊検査方法及び装置
CN103017695A (zh) * 2012-12-17 2013-04-03 西安交通大学 一种主轴锥孔精度的准确检测装置及检测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0342567A (ja) * 1989-07-10 1991-02-22 Hitachi Constr Mach Co Ltd 非破壊検査方法及び装置
CN103017695A (zh) * 2012-12-17 2013-04-03 西安交通大学 一种主轴锥孔精度的准确检测装置及检测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
FRANCESCO AYMERICH等: "《Evaluation of Nominal Contact Area and Contact Pressure Distribution in a Steel-Steel Interface by Means of Ultrasonic Techniques》", 《JSME INTERNATIONAL JOURNAL(SERIES C)》, vol. 46, no. 1, 31 December 2003 (2003-12-31), pages 297 - 305 *
R.S.DWYER-JOYCE1等: "《Experimental procedure for mapping Hertzian contacts》", 《SCIENCE IN CHINA (SERIES A)》, vol. 44, 31 August 2001 (2001-08-31), pages 423 - 430 *
杜飞等: "《结合面参数的超声检测方法研究》", 《西安交通大学学报》, vol. 47, no. 3, 31 March 2013 (2013-03-31), pages 18 - 23 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105277149A (zh) * 2015-11-18 2016-01-27 西安建筑科技大学 结合面真实接触面积测量装置与测量方法
CN105277149B (zh) * 2015-11-18 2017-11-10 西安建筑科技大学 结合面真实接触面积测量装置与测量方法
CN106932277A (zh) * 2017-04-28 2017-07-07 大连理工大学 基于圆角平面接触理论的界面超声波反射率‑压强关系曲线建立方法及加载试验台
WO2018196148A1 (zh) * 2017-04-28 2018-11-01 大连理工大学 基于圆角平面接触理论的界面超声波反射率-压强关系曲线建立方法及加载试验台
US11204291B2 (en) 2017-04-28 2021-12-21 Dalian University Of Technology Interface ultrasonic reflectivity-pressure relation curve establishment method and loading testbed
CN106932277B (zh) * 2017-04-28 2023-05-12 大连理工大学 基于圆角平面接触理论的界面超声波反射率-压强关系曲线建立方法及加载试验台
CN108908378A (zh) * 2018-05-14 2018-11-30 莱茵检测认证服务(中国)有限公司 一种人与机器人接触过程中接触压强的测试方法
CN108981989A (zh) * 2018-10-10 2018-12-11 郑州大学 一种检测轴套过盈配合结合面压强的装置及其使用方法
CN108981989B (zh) * 2018-10-10 2024-01-26 郑州大学 一种检测轴套过盈配合结合面压强的装置及其使用方法

Also Published As

Publication number Publication date
CN103822968B (zh) 2016-08-17

Similar Documents

Publication Publication Date Title
US11193838B2 (en) Method for determining plane stresses on in-service steel structure member based on phase spectrum of ultrasonic transverse wave
CN102636303B (zh) 一种基于表面超声波测定薄镀层残余应力的方法
CN103822968A (zh) 面向结合面压强检测的压强-超声反射率曲线构建方法
CN102608212B (zh) 基于声压反射系数功率谱测量薄层声阻抗的方法
CN104034287B (zh) 一种弹性各向异性金属基体热障涂层厚度超声测量方法
CN103543208B (zh) 基于频谱分析原理减小tofd检测近表面盲区的方法
CN111855809B (zh) 一种基于复合模式全聚焦的裂纹形貌重建方法
CN106198739A (zh) 一种基于波型转换的tofd近表面盲区缺陷定位检测方法
CN104698089A (zh) 一种适用于倾斜裂纹定量和成像的超声相对时间传播技术
CN104764803A (zh) 基于超声波波长变化的材料应力检测技术
CN105954356A (zh) 一种基于有限幅度法的金属块闭合裂纹检测定位方法
CN109387564A (zh) 一种基于锂枝晶生长的锂离子电池在线检测方法及装置
CN103616436B (zh) 一种接触刚度的高精度超声检测方法
CN107024535A (zh) 一种基于表面波的垂直缺陷的多系数深度检测方法
CN103075981A (zh) 一种超声波测厚方法
CN111678988A (zh) 混凝土材料表面损伤的非线性超声评价装置及方法
CN114910565A (zh) 一种非线性超声检测中相对非线性系数的修正方法
CN110646119A (zh) 一种超声波测量轧制金属材料表面应力张量的方法
CN105738017A (zh) 元素含量影响超声波评价金属材料表层应力的修正方法
US10261054B2 (en) Method for the non-destructive ultrasonic testing of a part by echo analysis
CN211783988U (zh) 基于超声测量的螺栓轴向应力测试装置
CN203117167U (zh) 公路钢桥用超声波检测装置
CN112484836A (zh) 一种超声探头装置及工件声速测量方法
CN114964581B (zh) 一种基于超声相控阵聚焦原理的应力检测方法
CN104458915A (zh) 一种风塔焊缝非线性超声检测方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant