CN109780965B - 基于同心圆网格遍历寻优的圆柱面多要素检测装置及方法 - Google Patents

基于同心圆网格遍历寻优的圆柱面多要素检测装置及方法 Download PDF

Info

Publication number
CN109780965B
CN109780965B CN201910085650.3A CN201910085650A CN109780965B CN 109780965 B CN109780965 B CN 109780965B CN 201910085650 A CN201910085650 A CN 201910085650A CN 109780965 B CN109780965 B CN 109780965B
Authority
CN
China
Prior art keywords
measured
section
circle
sensor
measuring
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.)
Expired - Fee Related
Application number
CN201910085650.3A
Other languages
English (en)
Other versions
CN109780965A (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 CN201910085650.3A priority Critical patent/CN109780965B/zh
Publication of CN109780965A publication Critical patent/CN109780965A/zh
Application granted granted Critical
Publication of CN109780965B publication Critical patent/CN109780965B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明公开了一种基于同心圆网格遍历寻优的圆柱面多要素检测装置及方法,采用动态旋转扫描测量,在行星架外圆上测量三个截面,每截面对径布置传感器;工件定位在精密回转平台上,对被测件进行高精度旋转定位;通过控制电机使左侧测量部由气缸导轨带动进入测量位,使传感器探头到达测量位置,实现在自动测量过程中的精准定位;工件在精密回转平台上缓慢旋转,每旋转一定角度,传感器采集一组数据,利用传感器测量得到的数据按照本文涉及的评定算法计算行星架外圆的外径、圆度、圆柱度。本发明提高了测量精度和测量效率,实现了在生产线上对行星架圆柱面的质量检测,保证了生产线的高效工作;实现对外圆面多要素的准确测量评定。

Description

基于同心圆网格遍历寻优的圆柱面多要素检测装置及方法
技术领域
本发明属于RV减速器行星架外圆圆柱面的加工精度检测技术领域,具体涉及一种基于同心圆网格遍历寻优的圆柱面多要素检测装置及方法。
背景技术
RV减速器由一个行星齿轮减速器的前级和一个摆线针轮减速器的后级组成,其中二级针摆传动部分由曲柄轴、摆线轮、针轮、针齿壳和行星架组件等组成。行星架轴承安装外圆形位误差关乎着减速器整体的装配精度,进而影响RV减速器的传动精度和传动效率。同时,较大的形位误差会使减速器运行过程中产生周期应力,影响减速器寿命。因此,对RV减速器的轴承安装外圆的形位误差检测至关重要。
行星架轴承安装外圆主要检测项目有外径尺寸、圆度、圆柱度、垂直度等。外径尺寸误差是轴承安装外圆的一个重要参数,直接影响轴承与行星架的配合精度,进而影响整机的传动精度;圆度误差和圆柱度误差影响减速器的工作性能,直接影响零件的互换性、配合精度、回转精度、摩擦性和使用寿命,进而影响RV减速器的传动精度和传动效率。
传统的测量方法以手动测量为主,测量仪器有专用测量仪、三坐标测量机测量。手动测量操作过程比较繁琐,基本都是离线测量,检测效率也比较低;专用测量仪测量精度高、性能稳定、操作及维护方便,但是由于技术难度较大,普及较少。三坐标测量机能精确测量出复杂结构形位误差,但是对测量环境有严格的要求,测量节拍较长,成本高,适合线下抽检,无法满足对生产线的在线测量。而面向生产线大批量多规格零件的多要素误差在线检测技术要待进一步的研究。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足,提供一种适应性强、准确度高、面向生产线的基于同心圆网格遍历寻优的圆柱面多要素检测装置及方法,应用于本文圆柱面多要素检测装置或辅助三坐标测量机,实现对RV减速器行星架轴承安装外圆面多几何要素的准确评定。
本发明采用以下技术方案:
一种基于同心圆网格遍历寻优的圆柱面多要素检测方法,包括以下步骤:
S1、以行星架基准底面为XOY平面,基准底面中心为原点,行星架理想外圆轴线为Z轴,建立坐标系;测量标准件或调用标准件传感器数据;
S2、动态旋转扫描测量时,精密回转平台转动,位移传感器采集数据,根据每个测量截面传感器的测量值,结合标准件数据,以采样角度间隔形成的离散点拟合出测量截面轮廓,得到各离散点的坐标信息,通过截面轮廓拟合测量截面圆及圆心,计算截面圆外径;利用最小区域法计算每个测量截面的圆度;
S3、根据步骤S2中得到的三个测量截面圆圆心,利用最小二乘法初步拟合被测件圆柱外圆轴线,计算最小二乘轴线两端点坐标;
S4、以步骤S3中得到的最小二乘法轴线两端点为基准圆心,在测量起始、终止截面内以合格件要求的圆柱度误差最大值为直径画圆,形成包络最小二乘轴线的斜柱体;以一定疏密程度对斜柱体上下底面进行同心圆网格划分,并计算同心圆网格节点坐标;
S5、两两连接斜柱体上下底面同心圆网格节点,形成直线群,通过带入同心圆网格节点坐标计算直线群中每条直线的方程,结合步骤S2中测量截面轮廓坐标,计算各测点距直线的距离及极差,在直线群中遍历搜索中最接近包容被测点的最小区域的理想轴线;利用该轴线计算圆柱面的圆柱度,判断圆柱度精度要求,若精度不满足要求,继续下面步骤;
S6、重新确定满足最小区域法的理想轴线的所在的网格范围,细化该部分网格,重复步骤S5,得到圆柱面的圆柱度。
具体的,步骤S1中,将标准件放置在精密回转平台,气缸带动测量部进入测量位,传感器探头与圆柱面接触,精密回转平台匀速转动,标准件旋转一周,每个传感器总共采集n组数据,对三个测量截面传感器进行编号,第i只传感器第j组示值记为Sbi[j],同一测量截面的两只传感器对径布置,当零件转过180°,零件同一位置被对径传感器重复测量;
同一截面两只传感器示值存在以下关系:
Figure GDA0002406616670000031
两传感器示值取平均值作为该截面传感器数据:
Figure GDA0002406616670000032
其中,j=1,2,…,n,当
Figure GDA0002406616670000033
时,取+号,当
Figure GDA0002406616670000034
时,取-号。
具体的,步骤S2中,将被测工件放置在精密回转平台,重复步骤S1中传感器采样步骤,得到被测件传感器数据,第i只传感器第j组示值记为Sci[j];同一截面两只传感器示值存在关系如下:
Figure GDA0002406616670000035
两传感器示值取平均值作为该截面传感器数据:
Figure GDA0002406616670000036
其中,j=1,2,…,n,当
Figure GDA0002406616670000037
时,取+号,当
Figure GDA0002406616670000038
时,取-号。
进一步的,根据标准件和被测工件所采集的传感器数据,计算每组传感器数值拟合的第k个测量截面圆半径ρk[j],拟合被测截面轮廓如下:
ρk[j]=Rb+(S'ci[j]-S'bi[j])
其中,Rb是标准件半径,根据得到的测量截面圆半径ρk[j]转换到直角坐标系,计算得到k截面轮廓采样点坐标(xkj,ykj,zk)如下:
Figure GDA0002406616670000041
利用改进最小二乘公式计算得到拟合圆圆心坐标Ok(xck,yck,zk)和半径Rck如下:
Figure GDA0002406616670000042
其中,f为最小二乘函数值,当f最小时,xck、yck、Rck的值即为拟合圆圆心坐标值和半径。
进一步的,以拟合圆圆心坐标Ok(xck,yck)为原点,画同心圆包容截面轮廓,测量截面的圆度fkR为:
fkR=Rmax-Rmin
其中,Rmax为包容截面轮廓的最大半径,Rmin为包容截面轮廓的最小半径。
具体的,步骤S3中,将空间直线标准方程简化为
Figure GDA0002406616670000043
计算拟合方程近似值与实际值之差:
Δx=xi-(k1·z+k2)
Δy=yi-(k3·z+k4)
根据最小二乘原理,对k1、k2、k3、k4进行求偏导,并令偏导等于零,代入三个测量截面圆心坐标O1(xc1,yc1,z1)、O2(xc2,yc2,z2)、O3(xc3,yc3,z3)求出方程中的参数k1、k2、k3、k4,根据被测圆柱面最小二乘轴线方程计算得到被测圆柱面最小二乘轴线端点坐标(x′c1,y′c1,z1)、(x′c3,y′c3,z3)。
具体的,步骤S4中,在测量起始、终止截面内以合格件要求的圆柱度误差最大值fC为直径画圆,形成包络最小二乘轴线的斜柱体;对斜柱体上下底面进行同心圆网格划分,设置同心圆个数m,以传感器测量旋转角度间隔设置轴向网格线,在斜柱体上下底面形成n·m个同心圆网格节点。
具体的,步骤S5中,在{dpqp'q'1max、dpqp'q'1min、dpqp'q'2max、dpqp'q'2min、dpqp'q'3max、dpqp'q'3min}中取最大值dpqp'q'max和最小值dpqp'q'min,为被测圆柱面的最小区域包容圆半径,以该直线为轴线的被测圆柱面圆柱度fpqp'q'为:
fpqp'q'=dpqp'q'max-dpqp'q'min
遍历斜柱体直线群中的(n·m)2条直线,得到(n·m)2个圆柱度组成的集合,取其中最小的fmin=min d'pqp'q',被测圆柱面圆柱度为fmin所对应的斜柱体下底面坐标为(x′p'q'1,y′p'q'1,z1),上底面坐标为(x′pq3,y′pq3,z3),此两点连线为对应的最优轴线。
进一步的,斜柱体上下底面每个同心圆网格节点坐标(xpq1,ypq1,z1)、(xpq3,ypq3,z3)计算如下:
Figure GDA0002406616670000051
Figure GDA0002406616670000052
其中,p为同心圆圈数,p=1,2,…,m,q为网格周向次序,q=1,2,…,n,fC为圆柱度误差最大值;
依次从斜柱体下底面上的网格点(x111,y111,z1)为起始点,遍历连接斜柱体上底面的网格点(xpq3,ypq3,z3),则构成的直线群中就会有(n·m)2条直线,每条直线的方程为:
Figure GDA0002406616670000061
计算斜柱体下底面截面轮廓采样点坐标(x1j,y1j,z1)到该直线的距离如下:
Figure GDA0002406616670000062
筛选出每个下底面截面轮廓采样点到直线距离的最大值dpqp'q'1max与最小值dpqp'q'1min
本发明的另一技术方案是,一种基于同心圆网格遍历寻优的圆柱面多要素检测装置,包括气缸固定座、精密导轨、气缸及滑块、测量部支撑座、限位挡块、传感器座、位移传感器和精密回转平台;气缸固定座设置在精密回转平台的一侧,测量部支撑座通过气缸及滑块固定在气缸固定座的精密导轨上,精密导轨的一端与精密回转平台之间设置有限位挡块,测量部支撑座上固定设置有传感器座,多个位移传感器设置在传感器座上,工件定位在精密回转平台上,精密回转平台上设置有圆光栅,能够对被测工件进行高精度旋转定位,采用动态旋转扫描测量,在工件外圆上测量三个截面,每截面对径布置2个位移传感器,2个位移传感器的测点在测量截面圆直径的延长线上。
与现有技术相比,本发明至少具有以下有益效果:
本发明一种基于同心圆网格遍历寻优的圆柱面多要素检测方法,在得到标准件信息后,只需将被测件放置该测量装置就能快速获得此被测圆柱面的外径、圆度、圆柱度等多参数信息;提高了测量精度和测量效率,实现了在生产线上对行星架圆柱面的质量检测,保证了生产线的高效工作;通过自行编写程序,也可辅助三坐标测量机,实现对外圆面多要素的准确测量评定。
进一步的,建立评定坐标系,采用比较测量法,通过调用或测量标准件得到标准件传感器数据,为后续计算被测件形位信息奠定基础。
进一步的,采用动态旋转扫描测量,快速采集得到被测件各测量截面传感器示值,结合标准件数据,拟合测量截面轮廓,通过计算能够准确的得到截面圆圆心、外径和圆度;同时为后续拟合圆柱外圆轴线奠定基础。
进一步的,通过得到的三个测量截面圆圆心,利用最小二乘法拟合得到圆柱外圆最小二乘轴线,计算得到该轴线两端点坐标,为后续同心圆网格划分奠定基础。
进一步的,以最小二乘法轴线两端点为基准圆心,以合格件要求的圆柱度误差最大值为直径画圆,形成包络最小二乘轴线的斜柱体;对斜柱体上下底面进行同心圆网格划分,计算得到同心圆网格节点坐标,为后续同心圆网格遍历寻优奠定基础。
进一步的,两两连接斜柱体上下底面同心圆网格节点,形成直线群,通过计算得到直线群中每条直线的方程,结合测量截面轮廓坐标,计算各测点距直线的距离及极差,在直线群遍历搜索最接近包容被测点的最小区域的理想轴线;利用该轴线即可得到满足最小条件原则的圆柱度。
本发明一种基于同心圆网格遍历寻优的圆柱面多要素检测装置,采用动态旋转扫描测量方式。被测工件定位在精密旋转轴系上,在被测外圆上测量三个截面,每截面对径布置2只传感器,测量时被测工件在精密回转平台上缓慢旋转,每旋转一定角度,传感器采集一组数据;根据每个测量截面传感器的测量值,结合标准件数据,拟合出测量截面轮廓和被测件圆柱外圆轴线,计算截面圆外径和圆度;通过同心圆网格划分,在包容拟合轴线的斜柱体区域内遍历搜索满足最小区域法的理想轴线,计算得到准确的圆柱度。同时,利用同心圆网格遍历寻优方法也可进行垂直度、同轴度的最优评定。
综上所述,本发明提供一种适应性强、准确度高、面向生产线的圆柱外圆面多要素在线检测装置及评定方法,在得到标准件信息后,只需将被测件放置该测量装置就能快速获得此被测圆柱面的外径、圆度、圆柱度等多参数信息;提高了测量精度和测量效率,实现了在生产线上对行星架圆柱面的质量检测,保证了生产线的高效工作;通过自行编写程序,也可辅助三坐标测量机,实现对外圆面多要素的准确测量评定。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为圆柱面多要素检测装置;
图2为行星架轴承安装外圆测量流程;
图3为被测行星架示意图;
图4为被测截面轮廓拟合图;
图5为被测截面圆度图;
图6为斜柱体同心圆网格图;
图7为网格节点坐标图。
其中:1.气缸固定座;2.精密导轨;3.气缸及滑块;4.支撑座;5.限位挡块;6.传感器座;7.位移传感器;8.工件;9.精密回转平台。
具体实施方式
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”、“一侧”、“一端”、“一边”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,在本发明的描述中,除非另有说明,“多个”的含义是两个或两个以上。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
请参阅图1,本发明提供了一种基于同心圆网格遍历寻优的圆柱面多要素检测装置,包括气缸固定座1、精密导轨2、气缸及滑块3、测量部支撑座4、限位挡块5、传感器座6、位移传感器7(三组,六个)和精密回转平台9,
气缸固定座1设置在精密回转平台9的一侧,测量部支撑座4通过气缸及滑块3固定在气缸固定座1的精密导轨2上,精密导轨2的一端与精密回转平台9之间设置有限位挡块5,测量部支撑座4上固定设置有传感器座6,多个位移传感器7设置在传感器座6上,工件8定位在精密回转平台9上,精密回转平台9上带有圆光栅,能够对被测件进行高精度旋转定位,采用动态旋转扫描测量,在工件8外圆上测量三个截面,每截面对径布置2个位移传感器7,2个位移传感器7的测点在测量截面圆直径的延长线上。
通过软件对电机进行控制,使左侧测量部由精密导轨2带动进入测量位,使位移传感器7的探头到达测量位置,实现在自动测量过程中的精准定位;工件8在精密回转平台9上缓慢旋转,每旋转一定角度,位移传感器7采集一组数据,利用位移传感器7测量得到的数据按照本发明测量方法采用的评定算法计算行星架外圆的外径、圆度、圆柱度等。
请参阅图2,本发明一种基于同心圆网格遍历寻优的圆柱面多要素检测方法,包括以下步骤:
S1、以行星架基准底面为XOY平面,基准底面中心为原点,行星架理想外圆轴线为Z轴,建立坐标系;测量标准件或调用标准件传感器数据;
S2、动态旋转扫描测量时,精密回转平台每转过一定角度,位移传感器采一组数据,根据每个测量截面传感器的测量值,结合标准件数据,以采样角度间隔形成的离散点拟合出测量截面轮廓,得到各离散点的坐标信息,通过截面轮廓拟合测量截面圆及圆心,计算截面圆外径;利用最小区域法计算每个测量截面的圆度;
S3、根据步骤S2中得到的三个测量截面圆圆心,利用最小二乘法初步拟合被测件圆柱外圆轴线,计算最小二乘轴线两端点坐标;
S4、以步骤S3中得到的最小二乘法轴线两端点为基准圆心,在测量起始、终止截面内以合格件要求的圆柱度误差最大值为直径画圆,形成包络最小二乘轴线的斜柱体;以一定疏密程度对斜柱体上下底面进行同心圆网格划分,并计算同心圆网格节点坐标;
S5、两两连接斜柱体上下底面同心圆网格节点,形成直线群,通过带入同心圆网格节点坐标计算直线群中每条直线的方程,结合步骤S2中测量截面轮廓坐标,计算各测点距直线的距离及极差,在直线群中遍历搜索最接近包容被测点的最小区域的理想轴线;利用该轴线计算圆柱面的圆柱度,判断圆柱度精度要求,若精度不满足要求,继续下面步骤;
S6、重新确定满足最小区域法的理想轴线的所在的网格范围,细化该部分网格,重复步骤S5,得到圆柱面的圆柱度。
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中的描述和所示的本发明实施例的组件可以通过各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图3,具体实施步骤如下:
1)根据零件设计图纸,以被测圆柱面下底面中心为坐标原点,圆柱面理想轴线为Z轴建立空间坐标系,根据传感器Z向布局得到各测量截面的Z轴坐标值zk
采用比较测量法,在对被测件进行检测之前,首先采集标准件信息,每批零件共用同一标准件信息,测量得到标准件信息后,将其信息保存。
具体采样方法如下:
将标准件放置在精密回转平台,气缸带动测量部进入测量位,使传感器探头与圆柱面接触。开启精密回转平台并以一定速度匀速转动,零件旋转一周,每个传感器总共采集n组数据,即每转过
Figure GDA0002406616670000115
传感器采一次样。对三个测量截面传感器进行编号,第i只传感器第j组示值记为Sbi[j],由于同一测量截面两只传感器对径布置,当零件转过180°,零件同一位置被对径传感器重复测量,所以同一截面两只传感器示值存在以下关系:
Figure GDA0002406616670000111
两传感器示值取平均值作为该截面传感器数据:
Figure GDA0002406616670000112
其中,j=1,2,…,n,当
Figure GDA0002406616670000113
时,取“+”号,当
Figure GDA0002406616670000114
时,取“-”号。
2)将被测工件放置在精密回转平台,重复(1)中传感器采样步骤,得到被测件传感器数据,第i只传感器第j组示值记为Sci[j];
同理,同一截面两只传感器示值存在以下关系:
Figure GDA0002406616670000121
两传感器示值取平均值作为该截面传感器数据:
Figure GDA0002406616670000122
其中,j=1,2,…,n,当
Figure GDA0002406616670000123
时,取“+”号,当
Figure GDA0002406616670000124
时,取“-”号。
3)根据标准件和被测件所采集的传感器数据,计算每组传感器数值拟合的第k个测量截面圆半径ρk[j],拟合被测截面轮廓,如图4所示。
ρk[j]=Rb+(S'ci[j]-S'bi[j])
其中,Rb是标准件半径。
4)将步骤(3)得到一系列测量截面圆半径ρk[j]转换到直角坐标系,得到k截面轮廓采样点坐标(xkj,ykj,zk):
Figure GDA0002406616670000125
利用最小二乘法原理,将每个截面轮廓点(xkj,ykj)代入最小二乘公式,求解当f最小时xck、yck、Rck的值,既拟合圆圆心坐标值和半径。
Figure GDA0002406616670000126
其中,f为最小二乘函数值,当f最小时,xck、yck、Rck的值即为拟合圆圆心坐标值和半径。
为得到解析解,在保持优化目标函数特征的前提上,对最小二乘公式进行改进,方便计算,利用下式计算得到拟合圆圆心坐标Ok(xck,yck,zk)和半径Rck,如图6所示。
Figure GDA0002406616670000131
5)以步骤(4)得到拟合圆圆心坐标Ok(xck,yck)为原点,画同心圆包容步骤(3)中的截面轮廓,如图5所示,该测量截面的圆度为:
fkR=Rmax-Rmin
其中,Rmax为包容截面轮廓的最大半径,Rmin为包容截面轮廓的最小半径。
6)利用最小二乘法将步骤(4)得到的三个测量截面圆心坐标O1(xc1,yc1,z1)、O2(xc2,yc2,z2)、O3(xc3,yc3,z3)拟合直线,得到被测圆柱面最小二乘轴线方程。
空间直线标准方程为:
Figure GDA0002406616670000132
简化可得
Figure GDA0002406616670000133
拟合方程近似值与实际值之差为
Δx=xi-(k1·z+k2)
Δy=yi-(k3·z+k4)
根据最小二乘原理,对k1、k2、k3、k4进行求偏导,并令偏导等于零,代入三个测量截面圆心坐标O1(xc1,yc1,z1)、O2(xc2,yc2,z2)、O3(xc3,yc3,z3)即可求出方程中的参数k1、k2、k3、k4,得到被测圆柱面最小二乘轴线方程。
根据步骤(1)可得,最小二乘轴线端点所在平面Z轴坐标为z1、z3,将其代入求得的最小二乘轴线方程,可得最小二乘轴线端点坐标(x′c1,y′c1,z1)、(x′c3,y′c3,z3)。
7)以最小二乘轴线两端点为圆心。在测量起始、终止截面内以合格件要求的圆柱度误差最大值fC为直径画圆,形成包络最小二乘轴线的斜柱体;以一定疏密程度对斜柱体上下底面进行同心圆网格划分,根据精度要求,设置同心圆个数m,以传感器测量旋转角度间隔设置轴向网格线(共点直线束,共
Figure GDA0002406616670000141
条过圆心的直线),在斜柱体上下底面形成n·m个同心圆网格节点,如图6所示。
8)两两连接斜柱体上下底面同心圆网格节点,形成直线群,根据圆柱度最小区域原理,该直线群中存在一条直线,以此为轴线可以形成空间直径差最小的两个同轴圆柱面包容被测实际圆柱面,在直线群搜索中最接近包容被测点的最小区域的理想轴线;
根据步骤(6),斜柱体上下底面分别有n·m个同心圆网格节点,如图7所示,计算斜柱体上下底面每个同心圆网格节点坐标(xpq1,ypq1,z1)、(xpq3,ypq3,z3)。
Figure GDA0002406616670000142
Figure GDA0002406616670000143
其中,p为同心圆圈数,p=1,2,…,m,q为网格周向次序,q=1,2,…,n,fC为圆柱度误差最大值。
依次从包容斜柱体下底面上的网格点(x111,y111,z1)为起始点,遍历连接斜柱体上底面的网格点(xpq3,ypq3,z3),则构成的直线群中就会有(n·m)2条直线。
每条直线的方程为:
Figure GDA0002406616670000144
计算斜柱体下底面采样点坐标(x1j,y1j,z1)到该直线的距离:
Figure GDA0002406616670000151
筛选出每个节点到采样点坐标到直线距离的最大值dpqp'q'1max和最小值dpqp'q'1min
同理可得中间采样截面采样点坐标(x2j,y2j,z2)到直线距离的最大值dpqp'q'2max和最小值dpqp'q'2min和斜柱体上底面采样点坐标(x3j,y3j,z3)到直线距离的最大值dpqp'q'3max和最小值dpqp'q'3min
在{dpqp'q'1max、dpqp'q'1min、dpqp'q'2max、dpqp'q'2min、dpqp'q'3max、dpqp'q'3min}中取最大值dpqp'q'max和最小值dpqp'q'min,即被测圆柱面的最小区域包容圆半径。所以以该直线为轴线的被测圆柱面圆柱度为:
fpqp'q'=dpqp'q'max-dpqp'q'min
遍历斜柱体直线群中的(n·m)2条直线。得到(n·m)2个圆柱度组成的集合,取其中最小的fmin=min d'pqp'q',所以被测圆柱面圆柱度为fmin,记该圆柱度所对应的斜柱体下底面坐标为(x'p'q'1,y'p'q'1,z1),上底面坐标为(x'pq3,y'pq3,z3),此两点连线即对应的最优轴线。
9)判断圆柱度精度是否达到要求,若不满足要求,重新确定满足最小区域法的理想轴线的所在的网格范围,以(8)中求得的轴线为两端点为原点,重复步骤(7)细化网格,重新进行同心圆网格遍历寻优,直至得到最优圆柱度。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。

Claims (4)

1.一种基于同心圆网格遍历寻优的圆柱面多要素检测方法,其特征在于,包括以下步骤:
S1、以行星架基准底面为XOY平面,基准底面中心为原点,行星架理想外圆轴线为Z轴,建立坐标系;测量标准件或调用标准件传感器数据,将标准件放置在精密回转平台,气缸带动测量部进入测量位,传感器探头与圆柱面接触,精密回转平台匀速转动,标准件旋转一周,每个传感器总共采集n组数据,对三个测量截面传感器进行编号,第i只传感器第j组示值记为Sbi[j],同一测量截面的两只传感器对径布置,当零件转过180°,零件同一位置被对径传感器重复测量;
同一截面两只传感器示值存在以下关系:
Figure FDA0002370290720000011
两传感器示值取平均值作为该截面传感器数据:
Figure FDA0002370290720000012
其中,j=1,2,…,n,当
Figure FDA0002370290720000013
时,取+号,当
Figure FDA0002370290720000014
时,取-号;
S2、动态旋转扫描测量时,精密回转平台转动,位移传感器采集数据,根据每个测量截面传感器的测量值,结合标准件数据,以采样角度间隔形成的离散点拟合出测量截面轮廓,得到各离散点的坐标信息,通过截面轮廓拟合测量截面圆及圆心,计算截面圆外径;利用最小区域法计算每个测量截面的圆度,将被测工件放置在精密回转平台,重复步骤S1中传感器采样步骤,得到被测件传感器数据,第i只传感器第j组示值记为Sci[j];同一截面两只传感器示值存在关系如下:
Figure FDA0002370290720000015
两传感器示值取平均值作为该截面传感器数据:
Figure FDA0002370290720000021
其中,j=1,2,…,n,当
Figure FDA0002370290720000022
时,取+号,当
Figure FDA0002370290720000023
时,取-号;
根据标准件和被测工件所采集的传感器数据,计算每组传感器数值拟合的第k个测量截面圆半径ρk[j],拟合被测截面轮廓如下:
ρk[j]=Rb+(S'ci[j]-S'bi[j])
其中,Rb是标准件半径,根据得到的测量截面圆半径ρk[j]转换到直角坐标系,计算得到k截面轮廓采样点坐标(xkj,ykj,zk)如下:
Figure FDA0002370290720000024
利用改进最小二乘公式计算得到拟合圆圆心坐标Ok(xck,yck,zk)和半径Rck如下:
Figure FDA0002370290720000025
其中,f为最小二乘函数值,当f最小时,xck、yck、Rck的值即为拟合圆圆心坐标值和半径;
S3、根据步骤S2中得到的三个测量截面圆圆心,利用最小二乘法初步拟合被测件圆柱外圆轴线,计算最小二乘轴线两端点坐标,将空间直线标准方程简化为
Figure FDA0002370290720000026
计算拟合方程近似值与实际值之差:
Δx=xi-(k1·z+k2)
Δy=yi-(k3·z+k4)
根据最小二乘原理,对k1、k2、k3、k4进行求偏导,并令偏导等于零,代入三个测量截面圆心坐标O1(xc1,yc1,z1)、O2(xc2,yc2,z2)、O3(xc3,yc3,z3)求出方程中的参数k1、k2、k3、k4,根据被测圆柱面最小二乘轴线方程计算得到被测圆柱面最小二乘轴线端点坐标(x′c1,y′c1,z1)、(x′c3,y′c3,z3);
S4、以步骤S3中得到的最小二乘法轴线两端点为基准圆心,在测量起始、终止截面内以合格件要求的圆柱度误差最大值为直径画圆,形成包络最小二乘轴线的斜柱体;以一定疏密程度对斜柱体上下底面进行同心圆网格划分,并计算同心圆网格节点坐标;
S5、两两连接斜柱体上下底面同心圆网格节点,形成直线群,通过带入同心圆网格节点坐标计算直线群中每条直线的方程,结合步骤S2中测量截面轮廓坐标,计算各测点距直线的距离及极差,在直线群中遍历搜索中最接近包容被测点的最小区域的理想轴线;利用该轴线计算圆柱面的圆柱度,判断圆柱度精度要求,若精度不满足要求,继续下面步骤;
在{dpqp'q'1max、dpqp'q'1min、dpqp'q'2max、dpqp'q'2min、dpqp'q'3max、dpqp'q'3min}中取最大值dpqp'q'max和最小值dpqp'q'min,为被测圆柱面的最小区域包容圆半径,以该直线为轴线的被测圆柱面圆柱度fpqp'q'为:
fpqp'q'=dpqp'q'max-dpqp'q'min
遍历斜柱体直线群中的(n·m)2条直线,得到(n·m)2个圆柱度组成的集合,取其中最小的fmin=min d'pqp'q',被测圆柱面圆柱度为fmin所对应的斜柱体下底面坐标为(x′p'q'1,y′p'q'1,z1),上底面坐标为(x′pq3,y′pq3,z3),此两点连线为对应的最优轴线;斜柱体上下底面每个同心圆网格节点坐标(xpq1,ypq1,z1)、(xpq3,ypq3,z3)计算如下:
Figure FDA0002370290720000041
Figure FDA0002370290720000042
其中,p为同心圆圈数,p=1,2,…,m,q为网格周向次序,q=1,2,…,n,fC为圆柱度误差最大值;
依次从斜柱体下底面上的网格点(x111,y111,z1)为起始点,遍历连接斜柱体上底面的网格点(xpq3,ypq3,z3),则构成的直线群中就会有(n·m)2条直线,每条直线的方程为:
Figure FDA0002370290720000043
计算斜柱体下底面截面轮廓采样点坐标(x1j,y1j,z1)到该直线的距离如下:
Figure FDA0002370290720000044
筛选出每个下底面截面轮廓采样点到直线距离的最大值dpqp'q'1max与最小值dpqp'q'1min
S6、重新确定满足最小区域法的理想轴线的所在的网格范围,细化该部分网格,重复步骤S5,得到圆柱面的圆柱度。
2.根据权利要求1所述的基于同心圆网格遍历寻优的圆柱面多要素检测方法,其特征在于,以拟合圆圆心坐标Ok(xck,yck)为原点,画同心圆包容截面轮廓,测量截面的圆度fkR为:
fkR=Rmax-Rmin
其中,Rmax为包容截面轮廓的最大半径,Rmin为包容截面轮廓的最小半径。
3.根据权利要求1所述的基于同心圆网格遍历寻优的圆柱面多要素检测方法,其特征在于,步骤S4中,在测量起始、终止截面内以合格件要求的圆柱度误差最大值fC为直径画圆,形成包络最小二乘轴线的斜柱体;对斜柱体上下底面进行同心圆网格划分,设置同心圆个数m,以传感器测量旋转角度间隔设置轴向网格线,在斜柱体上下底面形成n·m个同心圆网格节点。
4.一种利用权利要求1所述基于同心圆网格遍历寻优的圆柱面多要素检测方法的装置,其特征在于,包括气缸固定座(1)、精密导轨(2)、气缸及滑块(3)、测量部支撑座(4)、限位挡块(5)、传感器座(6)、位移传感器(7)和精密回转平台(9);气缸固定座(1)设置在精密回转平台(9)的一侧,测量部支撑座(4)通过气缸及滑块(3)固定在气缸固定座(1)的精密导轨(2)上,精密导轨(2)的一端与精密回转平台(9)之间设置有限位挡块(5),测量部支撑座(4)上固定设置有传感器座(6),多个位移传感器(7)设置在传感器座(6)上,工件(8)定位在精密回转平台(9)上,精密回转平台(9)上设置有圆光栅,能够对被测工件(8)进行高精度旋转定位,采用动态旋转扫描测量,在工件(8)外圆上测量三个截面,每截面对径布置2个位移传感器(7),2个位移传感器(7)的测点在测量截面圆直径的延长线上。
CN201910085650.3A 2019-01-29 2019-01-29 基于同心圆网格遍历寻优的圆柱面多要素检测装置及方法 Expired - Fee Related CN109780965B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910085650.3A CN109780965B (zh) 2019-01-29 2019-01-29 基于同心圆网格遍历寻优的圆柱面多要素检测装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910085650.3A CN109780965B (zh) 2019-01-29 2019-01-29 基于同心圆网格遍历寻优的圆柱面多要素检测装置及方法

Publications (2)

Publication Number Publication Date
CN109780965A CN109780965A (zh) 2019-05-21
CN109780965B true CN109780965B (zh) 2020-05-22

Family

ID=66503472

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910085650.3A Expired - Fee Related CN109780965B (zh) 2019-01-29 2019-01-29 基于同心圆网格遍历寻优的圆柱面多要素检测装置及方法

Country Status (1)

Country Link
CN (1) CN109780965B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110567350A (zh) * 2019-09-26 2019-12-13 中国航空工业集团公司北京航空精密机械研究所 一种基于蓝牙的内螺纹锥度检测系统及检测方法
CN112025339B (zh) * 2020-08-28 2021-12-14 中钢集团西安重机有限公司 一种确定回转罐体中心孔位置的方法
CN112254642A (zh) * 2020-10-29 2021-01-22 赵芃芃 断路器灭弧室几何轴线的检测方法
CN114659482B (zh) * 2020-12-22 2023-06-13 核动力运行研究所 一种基于非均匀离散数据的高精度圆形轮廓尺寸测量算法
CN114018202B (zh) * 2021-11-08 2024-02-02 绍兴职业技术学院 一种快速评定圆度的新算法
CN114111526B (zh) * 2021-11-19 2024-03-19 国网四川省电力公司映秀湾水力发电总厂 一种高精度水轮机主轴密封测圆方法及系统
CN115256388A (zh) * 2022-07-29 2022-11-01 江苏航鼎智能装备有限公司 一种基于随机抽样一致算法的工业机器人圆运动评价方法
CN116294896B (zh) * 2023-03-21 2024-04-26 山东欧瑞安电气有限公司 一种电机转子外径圆柱度检测方法
CN116448046B (zh) * 2023-06-16 2023-09-22 山东科技大学 一种圆柱形零件端面垂直度误差的测量装置及方法
CN117058220B (zh) * 2023-10-11 2024-01-26 天津施格自动化科技有限公司 利用距离传感器测量金属棒直径及中心点拟合算法
CN117824487A (zh) * 2024-03-04 2024-04-05 南京师范大学 一种管道机器人差速器工装高精度智能检测方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU472246A1 (ru) * 1972-01-04 1975-05-30 Куйбышевский сельскохозяйственный институт Устройство дл измерени погрешности расположени двух цилиндрических поверхностей деталей
SE464828B (sv) * 1987-01-26 1991-06-17 Gustav Bernhard Rennerfelt Roerelseoeverfoerande anordning
CN1272599C (zh) * 2005-02-04 2006-08-30 哈尔滨工业大学 一种可分离圆柱度仪基准误差的圆柱度测量方法及装置
CN104316017B (zh) * 2014-11-14 2017-05-31 苏州爱德蒙得测控系统有限公司 一种圆柱度的三截面测量方法
CN106168471A (zh) * 2015-12-04 2016-11-30 日照瑞沃德精密科技有限公司 基于圆度圆柱度的精密测量装置及其测量方法
CN105956303B (zh) * 2016-05-10 2018-12-07 西安交通大学 一种以抵消变形误差为目标的机床装配结合面设计方法

Also Published As

Publication number Publication date
CN109780965A (zh) 2019-05-21

Similar Documents

Publication Publication Date Title
CN109780965B (zh) 基于同心圆网格遍历寻优的圆柱面多要素检测装置及方法
CN102538700B (zh) 螺旋转子型面轮廓误差测量仪
CN111993159A (zh) 一种轴类工件在位非接触检测方法
CN107560583B (zh) 圆柱工件的轴心线校正方法及其分段截面的直径测量方法
CN102749041B (zh) 一种螺旋桨型面轮廓误差测量仪及测量方法
CN202350753U (zh) 螺旋转子型面轮廓误差测量仪
CN110470242B (zh) 一种大型零件内孔圆度在位测量装置及方法
CN201014954Y (zh) 内径自动检测仪
CN109059766B (zh) 一种深沟球轴承内圈沟位置的检测方法
US11754387B2 (en) Noncontact sensor calibration using single axis movement
CN115077452A (zh) 一种超高精度rv减速器曲柄轴的检测方法
CN110160454A (zh) 一种轴承装配游隙在线测量机及其游隙测量方法
CN105333828B (zh) 一种四轴联动跟踪式曲轴测量方法
CN117260389A (zh) 多传感器融合驱动的大型深孔零件形状误差在位测量系统
CN111829479A (zh) 一种用于测量零件深孔内表面形状误差的装置及测量方法
CN117146706A (zh) 一种轴类工件在位检测方法及系统
CN115752344B (zh) 一种rv减速器曲柄轴相位夹角的检测方法
CN107907071A (zh) 一种从动式轴承沟曲率半径测量装置及方法
CN201014953Y (zh) 全自动内径检测仪
CN211085044U (zh) 一种管壁结构件内外径测量装置
CN109827903B (zh) 一种活塞检测缺陷检测装置及方法
CN209961177U (zh) 用于测量内孔圆柱度的装置
CN112697056B (zh) 一种曲面薄壁件厚度的检测方法
CN113701643B (zh) 非平行面尺寸精确测量和圆周角度孔定位的测量方法
CN116448046B (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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200522