CN103822968A - 面向结合面压强检测的压强-超声反射率曲线构建方法 - Google Patents
面向结合面压强检测的压强-超声反射率曲线构建方法 Download PDFInfo
- 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
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:
其中,ΔAi为单个扫描点代表的面积,即整个结合面的面积除以扫描点数;F为计算得到的结合面受到的压力;将上述F与利用压力传感器测量到的实际压力F0带入下面式(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为超声反射率幅值,Hi为结合面信号幅值,Ho为相同频率下参考信号幅值;结合面所有点的反射率组成反射率矩阵,参见图2。图2的反射率矩阵中,中间圆环区域压强分布均匀,因此计算该区域的平均反射率R0。
5)计算根据压力计算结合面的平均压强,平均压强的80%即为中心区域压强。利用中心区域压强值P0及该区域反射率均值R0,拟合求得初始压强-反射率关系,如式(1)所示:
P0=f(R0) (1)
6)根据初始压强-反射率关系及整个结合面的反射率矩阵,计算得整个结合面的压力,如式(2)所示:
其中,ΔAi为单个扫描点代表的面积,即整个结合面的面积除以扫描点数。F为计算得到的结合面受到的压力。而此时整个面的实际压力F0可以由压力传感器测得,利用下面式(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:
其中,ΔAi为单个扫描点代表的面积,即整个结合面的面积除以扫描点数;F为计算得到的结合面受到的压力;将上述F与利用压力传感器测量到的实际压力F0带入下面式(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。
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)
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)
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 | 西安交通大学 | 一种主轴锥孔精度的准确检测装置及检测方法 |
-
2013
- 2013-11-08 CN CN201310551984.8A patent/CN103822968B/zh active Active
Patent Citations (2)
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)
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)
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 |