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

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

Info

Publication number
CN103822968B
CN103822968B CN201310551984.8A CN201310551984A CN103822968B CN 103822968 B CN103822968 B CN 103822968B CN 201310551984 A CN201310551984 A CN 201310551984A CN 103822968 B CN103822968 B CN 103822968B
Authority
CN
China
Prior art keywords
pressure
faying face
reflectance
scanning element
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.)
Active
Application number
CN201310551984.8A
Other languages
English (en)
Other versions
CN103822968A (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

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 (1)

1.一种面向结合面压强检测的压强-超声反射率曲线构建方法,其特征在于,包括以下步骤:
1)利用聚焦超声探头扫描结合面,计算不同压强下每个扫描点的超声反射率;然后利用结合面中心压强分布平坦区域扫描点的反射率均值R0及对应的压强P0,拟合求得初始压强-反射率关系,如式(1)所示:
P0=f(R0) (1)
反射率均值R0及对应的压强P0的获取方法具体如下:
(1)为检测螺栓结合面接触压强分布,按照与螺栓结合面相同的材料、表面粗糙度加工两个试件:上试件和下试件,将两试件安装到压力加载试验台上;
(2)获得参考信号:在结合面试件未接触的时候,利用奥林巴斯V319型15MHz水浸式聚焦超声波探头对上试件进行扫描,利用5800超声波脉冲发射接收器发射并接受超声脉冲信号,利用泰克MSO4104示波器进行接收并保存到电脑上,得到结合面上各个扫描点上的超声回波信号,该组信号为参考信号;
(3)保持上试件不动,让下试件与上试件接触,并施加不同的压力,其中压力的大小由压力传感器测得;按照第一次扫描的轨迹对结合面进行重新扫描,得到不同压力下结合面上各个扫描点的超声回波信号,该回波信号为结合面信号;
(4)分别对参考信号与结合面信号进行傅里叶变换,并计算每个扫描点在不同压力下的反射率;结合面所有点的反射率组成反射率矩阵,得到反射率均值R0
(5)根据压力计算结合面的平均压强,平均压强的80%即为中心区域压强,即反射率均值R0对应的压强P0
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重新拟合求得新的压强-反射率曲线。
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 CN103822968A (zh) 2014-05-28
CN103822968B true 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)

Families Citing this family (5)

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

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103017695A (zh) * 2012-12-17 2013-04-03 西安交通大学 一种主轴锥孔精度的准确检测装置及检测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2761928B2 (ja) * 1989-07-10 1998-06-04 日立建機株式会社 非破壊検査方法及び装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103017695A (zh) * 2012-12-17 2013-04-03 西安交通大学 一种主轴锥孔精度的准确检测装置及检测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《Evaluation of Nominal Contact Area and Contact Pressure Distribution in a Steel-Steel Interface by Means of Ultrasonic Techniques》;Francesco AYMERICH等;《JSME International Journal(Series C)》;20031231;第46卷(第1期);第297-305页 *
《Experimental procedure for mapping Hertzian contacts》;R.S.Dwyer-Joyce1等;《SCIENCE IN CHINA (Series A)》;20010831;第44卷;第423-430页 *
《结合面参数的超声检测方法研究》;杜飞等;《西安交通大学学报》;20130331;第47卷(第3期);第18-23页 *

Also Published As

Publication number Publication date
CN103822968A (zh) 2014-05-28

Similar Documents

Publication Publication Date Title
CN106989703B (zh) 一种全范围润滑膜厚度超声检测装置及方法
CN103822968B (zh) 面向结合面压强检测的压强-超声反射率曲线构建方法
CN102636303B (zh) 一种基于表面超声波测定薄镀层残余应力的方法
CN104142195A (zh) 基于超声波法的钢结构构件内部初始应力检测装置和方法
Hall et al. Minimum variance guided wave imaging in a quasi-isotropic composite plate
CN103292753B (zh) 采用超声波水浸聚焦技术测量热障涂层厚度的方法
CN102608212B (zh) 基于声压反射系数功率谱测量薄层声阻抗的方法
CA2608483A1 (en) Ultrasonic inspection method
CN103293225B (zh) 锅炉水冷壁管氢损伤超声波检测与诊断方法
IL281888B1 (en) An image reconstruction method based on non-linear mapping
CN106932277B (zh) 基于圆角平面接触理论的界面超声波反射率-压强关系曲线建立方法及加载试验台
CN104764803A (zh) 基于超声波波长变化的材料应力检测技术
CN106198739A (zh) 一种基于波型转换的tofd近表面盲区缺陷定位检测方法
CN103822971A (zh) 一种超声显微镜分辨力测试及校准方法
Ma et al. The reflection of guided waves from simple dents in pipes
CN103616436B (zh) 一种接触刚度的高精度超声检测方法
CN103543208B (zh) 基于频谱分析原理减小tofd检测近表面盲区的方法
CN104698089A (zh) 一种适用于倾斜裂纹定量和成像的超声相对时间传播技术
CN105954356A (zh) 一种基于有限幅度法的金属块闭合裂纹检测定位方法
CN114994175A (zh) 模态分解双谱分析的空耦超声应力检测装置及方法
JP4332531B2 (ja) 膜厚較正曲線の取得方法
EP3758608A1 (fr) Procédé d'élastographie hybride, sonde et dispositif pour élastographie hybride
CN111238702A (zh) 基于超声测量的螺栓轴向应力测试装置与测试方法
CN103075981A (zh) 一种超声波测厚方法
US10261054B2 (en) Method for the non-destructive ultrasonic testing of a part by echo analysis

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