CN113704908B - 一种曲轴连杆颈数据的处理方法 - Google Patents

一种曲轴连杆颈数据的处理方法 Download PDF

Info

Publication number
CN113704908B
CN113704908B CN202110923239.6A CN202110923239A CN113704908B CN 113704908 B CN113704908 B CN 113704908B CN 202110923239 A CN202110923239 A CN 202110923239A CN 113704908 B CN113704908 B CN 113704908B
Authority
CN
China
Prior art keywords
connecting rod
sampling
rod neck
data
crankshaft connecting
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
CN202110923239.6A
Other languages
English (en)
Other versions
CN113704908A (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202110923239.6A priority Critical patent/CN113704908B/zh
Publication of CN113704908A publication Critical patent/CN113704908A/zh
Application granted granted Critical
Publication of CN113704908B publication Critical patent/CN113704908B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开了一种曲轴连杆颈数据的处理方法,包括:根据曲轴连杆颈测量过程中实际采样角度分布,建立真实采样角度分布函数;采用比较法将被测轴颈与校准件进行比较,从而确定被测轴颈的圆轮廓数据;根据采样角度实际分布情况计算每个采样角对应的权函数;将每个采样角对应的权函数与该角度对应的采集数据进行时域离散循环卷积计算,所有的循环卷积再求和获得滤波后有效的圆轮廓数据;该处理方法通过回转轴线的运动角速度以及测头的采样间隔时间获得的采样角分布来逐点求取滤波权函数,解决了现有的数据处理方法与采样数据实际分布不相符的问题,精确实现连杆颈圆轮廓数据的滤波。

Description

一种曲轴连杆颈数据的处理方法
技术领域
本发明涉及精密测量领域,主要涉及一种基于非等间隔采样的曲轴连杆颈数据处理方法。
背景技术
利用随动式曲轴轴颈同步测量方法检测曲轴时,曲轴绕主轴颈轴线旋转模拟工况,主轴颈测头不动,被动跟随主轴颈一周采集数据;连杆颈测量采用随动式测量结构,测头相对连杆颈旋转一周采集数据,通过两测量结构的结合实现曲轴主轴颈与连杆颈的同步测量。主轴颈做回转运动的同时,连杆颈做绕回转中心的圆周运动,连杆颈的伸缩测量机构跟随连杆颈的运动而摆动,而顶端的弹性块测头在连杆颈表面运动一周完成截面径向数据的采集。由于设计结构的原因,曲轴围绕回转中心是匀速旋转,连杆颈的随动测头相对于回转中心并非匀速运动,所以连杆颈测头采集到的连杆颈表面圆周轮廓数据是非等间隔的采样样本。因此,基于等间隔样本假设和频域傅立叶变换的方法已不适用,需根据每个样本具体采样角度的实际分布进行操作,加以离散高斯加权函数进行处理。
发明内容
为了解决上述问题,本发明提供能够精确实现连杆颈圆轮廓数据的滤波的数据处理方法。
为了实现上述目的,本发明提供的技术方案是:一种曲轴连杆颈数据的处理方法,包括:
步骤1:根据曲轴连杆颈测量过程中实际采样角度分布,建立真实采样角度分布函数,为公式(1)所示:
式中,R为曲柄臂长度,L为测量导轨与主轴颈中心的距离,α为曲轴转过的角度,θ为曲轴连杆颈测头相对于回转中心的采样角度;
步骤2:采用比较法将被测曲轴连杆颈与校准件进行比较,从而确定被测曲轴连杆颈的圆轮廓数据;
步骤3:根据采样角度实际分布情况计算每个采样角对应的权函数,为公式(2)所示:
式中,n为采样点数,θi,θj表示不同的采样角度,Nc为截止波数,α为常数;exp表示以e为底的指数函数;
步骤4:将每个采样角对应的权函数与该角度对应的采集数据进行时域离散循环卷积计算,所有的循环卷积再求和获得滤波后有效曲轴连杆颈的圆轮廓数据。
作为优选的一种技术方案,所述步骤2具体通过公式(3)计算得出:
r=r0+Δd (3)
式中,r表示被测曲轴连杆颈的半径,r0是校准件曲轴连杆颈半径,Δd是实际测量曲轴连杆颈时测头相对于校准件的位移量。
作为优选的一种技术方案,所述步骤3进一步包括:
在等间隔采样情况下,高斯低通滤波器用于工件表面纹理分析的权函数由公式(4)给出:
式中x是离开权函数对称中心的距离,λc是高斯低通滤波器的截止波长,α是常数;
用截止波数Nc代替截止波长λc,再将权函数进行离散化从而获得每个采样角对应的权函数。
作为优选的一种技术方案,所述步骤4进一步包括:
利用公式(5)对所述权函数进行归一化处理:
在空间域中采用离散循环卷积和的运算进行高斯滤波,利用公式(6)得到非等间隔高斯滤波后的数据:
式中,r(j)是原始轮廓数据,y(i)是经过非等间隔权函数过滤后的数据。
作为优选的一种技术方案,所述曲轴连杆颈真实采样角度分布函数是在曲轴连杆颈做偏心运动实验中获取的。
本发明相对于现有技术的有益效果是:该处理方法通过回转轴线的运动角速度以及测头的采样间隔时间获得的采样角分布来逐点求取滤波权函数,解决了现有的数据处理方法与采样数据实际分布不相符的问题,精确实现连杆颈圆轮廓数据的滤波。
附图说明
图1是本发明提供的一种曲轴连杆颈数据的处理方法的流程图;
图2和图3是夹角θ随主轴颈运动角α的函数图像;
图4是连杆颈的采样轮廓显示和非等间隔滤波后的轮廓显示图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参照图1,图1本发明提出的一种曲轴连杆颈数据的处理方法的流程图,本实施例提供一种曲轴连杆颈数据的处理方法,具体包括如下步骤:
步骤1:根据曲轴连杆颈测量过程中实际采样角度分布,建立真实采样角度分布函数,为公式(1)所示:
式中,R为曲柄臂长度,L为测量导轨与主轴颈中心的距离,α为曲轴转过的角度,θ为曲轴连杆颈测头相对于回转中心的采样角度。
参照图2,图2是夹角θ随主轴颈运动角α的函数图像,其中逆时针方向为正,顺时针方向为负,主轴颈旋转一周,α由0°到360°变化过程中,θ相对于曲柄臂的角度由-90°变化到270°。由此可知在主轴颈绕回转中心旋转一周的过程中,测头在被测连杆颈圆周表面扫描一周并采集连杆颈径向尺寸变化量。
由于设计结构的原因,曲轴围绕回转中心是匀速旋转,连杆颈的随动测头相对于回转中心并非匀速运动。发现夹角θ与主轴颈运动角α之间并非线性关系,对夹角θ求导,得到其角速度的关系曲线,如图3所示,发现测头在连杆颈圆周表面的运动并非匀速运动,导致采样间隔不相等,即采样角的分布在整个运动过程中是非等间隔的。
步骤2:采用比较法将被测轴颈与校准件进行比较,从而确定被测轴颈的圆轮廓数据;
通过校准件的标量加上测头的位移量,获取采集数据为公式(2)所示:
r=r0+Δd (2)
式中,r表示被测曲轴连杆颈的半径,r0是校准件曲轴连杆颈半径,Δd是实际测量曲轴连杆颈时测头相对于校准件的位移量。
步骤3:等间隔采样情况下,高斯低通滤波器用于工件表面纹理分析的权函数可以由公式(3)给出:
式中x是离开权函数对称中心的距离,λc是高斯低通滤波器的截止波长,α是常数;exp表示以e为底的指数函数。
用截止波数Nc代替截止波长λc,再将权函数进行离散化。与等间隔采样相比,对于非等间隔的采样角度,高斯滤波的权函数由测量过程中连杆颈测头间隔的具体采样角决定,按照采样角度的实际分布,权函数是不同的,因此可以获得每个采样角对应的权函数。
根据采样角度实际分布情况计算每个采样角对应的权函数,为公式(4)所示:
式中,n为采样点数,θi,θj表示不同的采样角度,Nc为截止波数,α为常数。
步骤4:将每个采样角对应的权函数与该角度对应的采集数据进行时域离散循环卷积计算,所有的循环卷积再求和获得滤波后有效的圆轮廓数据。
对权函数进行归一化处理由公式(5)给出:
非等间隔情况下采集的样本使用角度信息计算滤波器的系数,还要与对应的幅度信息相乘,直接在空间域中采用离散循环卷积和的运算进行高斯滤波,得到非等间隔高斯滤波后的数据由公式(6)给出:
式中,r(j)是原始轮廓数据,y(i)是经过非等间隔权函数过滤后的数据。
图4所示为连杆颈的采样轮廓显示和非等间隔滤波后的轮廓显示。可以发现随着曲轴的回转运动,连杆颈上的采样点是先稀疏后密集再稀疏的过程。对非等间隔采样角度下测量到的数据进行重新排列,直接在空间域中采用不同高斯权函数进行高斯滤波,最大程度接近连杆颈圆轮廓,在保证轴颈同步测量的同时大大提高测量精度。
另外,本发明实施例还提供一种计算机可读存储介质,其中,该计算机可读存储介质可存储有程序,该程序执行时包括上述方法实施例中记载的曲曲轴连杆颈数据的处理方法的部分或全部步骤。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储器中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储器中,包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储器包括:U盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序可以存储于一计算机可读存储器中,存储器可以包括:闪存盘、只读存储器(英文:Read-Only Memory,简称:ROM)、随机存取器(英文:Random Access Memory,简称:RAM)、磁盘或光盘等。
以上参照附图描述了根据本发明的实施例的用于实现处理曲轴连杆颈数据的示例性流程图。应指出的是,以上描述中包括的大量细节仅是对本发明的示例性说明,而不是对本发明的限制。在本发明的其他实施例中,该方法可具有更多、更少或不同的步骤,且各步骤之间的顺序、包含、功能等关系可以与所描述和图示的不同。

Claims (6)

1.一种曲轴连杆颈数据的处理方法,其特征在于,包括:
步骤1:根据曲轴连杆颈测量过程中实际采样角度分布,建立真实采样角度分布函数,为公式(1)所示:
式中,R为曲柄臂长度,L为测量导轨与主轴颈中心的距离,α为曲轴转过的角度,θ为曲轴连杆颈测头相对于回转中心的采样角度;
步骤2:采用比较法将被曲轴连杆颈与校准件进行比较,从而确定被测曲轴连杆颈的圆轮廓数据;
步骤3:根据采样角度实际分布情况计算每个采样角对应的权函数,为公式(2)所示:
式中,n为采样点数,θi,θj表示不同的采样角度,Nc为截止波数,α为常数;exp表示以e为底的指数函数;
步骤4:将每个采样角对应的权函数与该角度对应的采集数据进行时域离散循环卷积计算,所有的循环卷积再求和获得滤波后有效的曲轴连杆颈的圆轮廓数据。
2.根据权利要求1所述的处理方法,其特征在于,所述步骤2具体通过公式(3)计算得出:
r=r0+Δd (3)
式中,r表示被测曲轴连杆颈的半径,r0是校准件曲轴连杆颈半径,Δd是实际测量曲轴连杆颈时测头相对于校准件的位移量。
3.根据权利要求1所述的处理方法,其特征在于,所述步骤3进一步包括:
在等间隔采样情况下,高斯低通滤波器用于工件表面纹理分析的权函数由公式(4)给出:
式中x是离开权函数对称中心的距离,λc是高斯低通滤波器的截止波长,α是常数;
用截止波数Nc代替截止波长λc,再将权函数进行离散化从而获得每个采样角对应的权函数。
4.根据权利要求1所述的处理方法,其特征在于,所述步骤4进一步包括:
利用公式(5)对所述权函数进行归一化处理:
在空间域中采用离散循环卷积和的运算进行高斯滤波,利用公式(6)得到非等间隔高斯滤波后的数据:
式中,r(j)是原始轮廓数据,y(i)是经过非等间隔权函数过滤后的数据。
5.根据权利要求1所述的处理方法,其特征在于:所述曲轴连杆颈真实采样角度分布函数是在曲轴连杆颈做偏心运动实验中获取的。
6.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至5任一项所述的一种曲轴连杆颈数据的处理方法的步骤。
CN202110923239.6A 2021-08-12 一种曲轴连杆颈数据的处理方法 Active CN113704908B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110923239.6A CN113704908B (zh) 2021-08-12 一种曲轴连杆颈数据的处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110923239.6A CN113704908B (zh) 2021-08-12 一种曲轴连杆颈数据的处理方法

Publications (2)

Publication Number Publication Date
CN113704908A CN113704908A (zh) 2021-11-26
CN113704908B true CN113704908B (zh) 2024-06-11

Family

ID=

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105841657A (zh) * 2016-03-25 2016-08-10 北京科技大学 一种金属表面形貌高斯滤波过程中的边缘效应抑制方法
CN105922133A (zh) * 2016-05-06 2016-09-07 上海交通大学 曲轴连杆颈随动磨削力动态测量方法及测量装置
CN108763803A (zh) * 2018-06-04 2018-11-06 东华大学 基于曲轴连杆颈随动磨削轮廓误差分解的磨床分析方法
CN109977352A (zh) * 2019-03-21 2019-07-05 哈尔滨工业大学 基于非等间隔滤波有效性处理的大型高速回转装备单级零部件跳动数据处理方法
CN110202418A (zh) * 2019-05-30 2019-09-06 中国人民解放军国防科技大学 轴类零件外圆表面的确定性砂带修形方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105841657A (zh) * 2016-03-25 2016-08-10 北京科技大学 一种金属表面形貌高斯滤波过程中的边缘效应抑制方法
CN105922133A (zh) * 2016-05-06 2016-09-07 上海交通大学 曲轴连杆颈随动磨削力动态测量方法及测量装置
CN108763803A (zh) * 2018-06-04 2018-11-06 东华大学 基于曲轴连杆颈随动磨削轮廓误差分解的磨床分析方法
CN109977352A (zh) * 2019-03-21 2019-07-05 哈尔滨工业大学 基于非等间隔滤波有效性处理的大型高速回转装备单级零部件跳动数据处理方法
CN110202418A (zh) * 2019-05-30 2019-09-06 中国人民解放军国防科技大学 轴类零件外圆表面的确定性砂带修形方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
高精度曲轴轴颈同步测量及误差分析研究;钱晓明;崔巍;楼佩煌;;仪器仪表学报;20200315(第03期);全文 *

Similar Documents

Publication Publication Date Title
JP5711145B2 (ja) 幾何学的タイヤ測定値のデータ品質を向上させるフィルタリング処理方法
US9322631B2 (en) Surface measurement instrument and calibration thereof
EP2541191B1 (en) Method of measuring a circle shape characteristic and circular shape characteristic measuring device and program
CN107917676B (zh) 一种基于条纹图像频谱分析的干涉测量方法
CN109716077A (zh) 使用叶尖定时(btt)监测涡轮机转子叶片的方法和系统
CN107196661B (zh) 用于转子系统的振动信号处理的信号重构方法
CN113447670A (zh) 一种单目线阵视觉追踪条纹编码标志的转速测量方法及系统
Huang et al. A pre-processing method for digital image correlation on rotating structures
CN113704908B (zh) 一种曲轴连杆颈数据的处理方法
US10444000B2 (en) Surface measurement instrument and method
CN113704908A (zh) 一种曲轴连杆颈数据的处理方法
CN109916361B (zh) 一种无需角度位置信息的圆度测量信号处理方法
CN116255930A (zh) 一种基于点云切片的截面提取及测量方法及系统
CN115930839A (zh) 一种基于经验模态分解的实时自适应动态角度测量方法和系统
CN110779471B (zh) 同轴度测量方法、装置、设备及存储介质
CN105133472B (zh) 路面平整度测量方法及惯性式检测仪
CN103471528B (zh) 一种干涉条纹倾斜角测量方法
Rossi A form of deviation-based method for coordinate measuring machine sampling optimization in an assessment of roundness
CN108534735B (zh) 基于v块的三点圆柱度测量装置及圆柱度误差分离方法
Gillich et al. Problem of detecting damage through natural frequency changes
CN110146025B (zh) 一种光栅尺的位移测量装置、方法及系统
CN112668125A (zh) 提高非完整小圆弧评价精度的方法、系统、介质及设备
CN112700387A (zh) 激光数据处理方法、装置及设备、存储介质
CN116701920B (zh) 一种提取oct色散失配系数的方法
CN114720110B (zh) 一种基于叶端定时传感器循环多尺度逆向优化布局的频率识别方法

Legal Events

Date Code Title Description
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant