CN111272088B - 滚动直线导轨副滑块型面中径的测量算法 - Google Patents

滚动直线导轨副滑块型面中径的测量算法 Download PDF

Info

Publication number
CN111272088B
CN111272088B CN202010157674.8A CN202010157674A CN111272088B CN 111272088 B CN111272088 B CN 111272088B CN 202010157674 A CN202010157674 A CN 202010157674A CN 111272088 B CN111272088 B CN 111272088B
Authority
CN
China
Prior art keywords
displacement sensor
laser displacement
axis
coordinate system
rolling
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
CN202010157674.8A
Other languages
English (en)
Other versions
CN111272088A (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 Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN202010157674.8A priority Critical patent/CN111272088B/zh
Publication of CN111272088A publication Critical patent/CN111272088A/zh
Application granted granted Critical
Publication of CN111272088B publication Critical patent/CN111272088B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/08Measuring arrangements characterised by the use of optical techniques for measuring diameters

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种滚动直线导轨副滑块型面中径的测量算法,将传感器置于滑块外侧两端,对滑块内滚道倾斜扫掠测量;对于滑块的大基准面与侧基准面,分别通过两个激光位移传感器测出其平面度,通过倾斜扫掠滑块内滚道得到半个椭圆轮廓数据点,经过最小二乘法进行椭圆拟合,求解出拟合椭圆的形心位置,改变滑块位置,扫掠出滑块不同截面的半椭圆轮廓,求解出对应的拟合椭圆形心位置,将同一滚道的拟合椭圆形心位置坐标点通过数学分析,计算出滚道中心轴线的直线方程,计算滑块左右滚道中心轴线的空间距离即为滚动直线导轨副滑块中径,进而测定滑块大基准面与侧基准面的平行度。为滚动直线导轨副的中径测量提出了一套完整的算法。

Description

滚动直线导轨副滑块型面中径的测量算法
技术领域
本发明属于滚动直线导轨副滑块检测领域,具体涉及一种滚动直线导轨副滑块型面中径的测量算法。
背景技术
滑块是滚动直线导轨副的一个重要构件,滚动直线导轨副的滑块中径尺寸及位置精度要求高,直接影响了钢球和滚道的接触角,进而影响直线导轨副预紧力,最终改变滚动直线导轨副的使用性能。从提高国内滚动直线导轨副的产品出发,很有必要对滚动直线导轨副滑块型面中径进行检测。
滚动直线导轨副滑块中径的测量主要分为接触式测量与非接触式测量。接触式测量,这种方法精度高,但是操作繁琐受认为因素影响大,且无法实现自动化测量。非接触式测量方法,通过激光位移传感器扫掠滚道,在计算机中拟合出滑块中径。
由于滚动直线导轨副滑块的跨距小,激光位移传感器尺寸大,二者容易产生干涉,无法将传感器置于滑块凹槽内部,直接对滑块中径进行测量。目前采用非接触式方法对滚动直线导轨副滑块中径检测仍处于一片空白。
发明内容
本发明的目的在于提供一种滚动直线导轨副滑块型面中径的测量算法。
本发明所采用的技术方案是:
一种滚动直线导轨副滑块型面中径的测量算法,包括以下步骤:
步骤1、构建坐标系:构建夹具空间直角坐标系、待测滑块的空间直角坐标系、激光位移传感器***组成的空间直角坐标系、第一激光位移传感器和第二激光位移传感器的空间直角坐标系、第一激光位移传感器和第二激光位移传感器的空间斜坐标系;
构建夹具的空间直角坐标系,具体为夹具上固连的标定块的空间直角坐标系o0-x0y0z0,其中x0轴和夹具台面垂直,y0轴沿着夹具台的长度方向,z0轴沿着夹具台的宽度方向,x0轴、y0轴、z0轴遵循右手法则;
构建待测滑块的空间直角坐标系o1-x1y1z1,其中x1轴与待测滑块的侧基准面垂直,y1轴方向与待测滑块导向方向相同,z1轴垂直于待测滑块的大基准面,x1轴、y1轴、z1轴遵循右手法则;
构建激光位移传感器***组成的空间直角坐标系o2-x2y2z2,其中y2轴方向与y1轴方向一致,x2轴在竖直面内与y2轴垂直,指向向上,z2轴与x2轴、y2轴遵循右手法则;
构建第一激光位移传感器的空间直角坐标系o3-x3y3z3,其中z3轴为第一激光位移传感器运动方向,x3轴垂直于z3轴,且位于z3轴、激光射出线组成的面内,y3轴与z3轴、x3轴遵循右手法则;
构建第一激光位移传感器的空间斜坐标系o4-x4y4z4,实际采集数据过程中,数据点在第一激光位移传感器的空间斜坐标系o4-x4y4z4中采集,其中z4轴即z3轴,y4轴即y3轴,x4轴代表传感器光线射出方向,其与x3轴之间的夹角为安装倾斜角β1
构建第二激光位移传感器的空间直角坐标系o5-x5y5z5,其中z5轴为第二激光位移传感器运动方向,x5轴垂直于z5轴,且位于z5轴、激光射出线组成的面内,y5轴与z5轴、x5轴遵循右手法则;
构建第二激光位移传感器的空间斜坐标系o6-x6y6z6,实际采集数据过程中,数据点在第二激光位移传感器的空间斜坐标系o6-x6y6z6中采集,z6轴即z5轴,y6轴即y5轴,x6代表传感器光线射出方向,其与x5轴之间的夹角为安装倾斜角β2
所述激光位移传感器***中的第一激光位移传感器、第二激光位移传感器沿z2轴方向测量待测滑块的两个滚道,滑块托架上固连标定圆柱、待测滑块以及标定块,测量过程中,滑块托架带动标定圆柱、待测滑块以及标定块沿y2轴方向移动,第三激光位移传感器、第四激光位移传感器、第五激光位移传感器、第六激光位移传感器固连在实验台上,其中第三激光位移传感器、第四激光位移传感器用于采集待测滑块和标定块侧基准面数据,第五激光位移传感器、第六激光位移传感器用于采集待测滑块以及标定块大基准面的数据;
步骤2、利用标定块得到夹具空间直角坐标系下的***误差曲线:侧基准面的第三激光位移传感器、第四激光位移传感器与大基准面的第五激光位移传感器、第六激光位移传感器测量标定块的两个与夹具配合的面,得到大基准面与侧基准面各两条直线数据,作为夹具***的误差曲线,用于后续补偿;
步骤3、确定待测滑块大基准面与侧基准面的平面度:第三激光位移传感器、第四激光位移传感器、第五激光位移传感器、第六激光位移传感器测量待测滑块的大基准面与侧基准面,得到待测滑块大基准面与侧基准面各两条直线数据,并将步骤2中得到的夹具***误差曲线作为补偿曲线,得到待测滑块大基准面、侧基准面在夹具空间直角坐标系下的坐标值,利用平面的最小二乘拟合方程,求出待测滑块的大基准面和侧基准面的平面度;
步骤4、第一激光位移传感器和第二激光位移传感器扫掠标定圆柱,利用椭圆的最小二乘法拟合,求解出标定椭圆一般方程;
步骤5、构建步骤4中所得椭圆一般方程与第一激光位移传感器和第二激光位移传感器安装倾斜角β1、β2并求解出所述安装倾斜角β1、β2:所述安装倾斜角β1是在第一激光位移传感器射出光线与第一激光位移传感器运动方向构成的平面内第一激光位移传感器运动方向的垂直线与第一激光位移传感器射出光线之间的夹角;安装倾斜角β2是第二激光位移传感器射出光线与第二激光位移传感器运动方向构成的平面内第二激光位移传感器运动方向的垂直线与第二激光位移传感器射出光线之间的夹角;
步骤6、确定第一激光位移传感器和第二激光位移传感器绕自身运动方向轴旋转的安装偏转角α1、α2并求解出所述安装偏转角α1、α2;所述安装偏转角α1是x3轴方向与滑块托架运动方向y2轴之间的夹角;安装偏转角α2是x5轴方向与滑块托架运动方向y2轴之间的夹角;
步骤7、计算待测滑块内滚道中心点坐标:第一激光位移传感器和第二激光位移传感器扫掠待测滑块内滚道,采集分别位于第一激光位移传感器和第二激光位移传感器空间斜坐标系下的待测滑块滚道数据点,采用椭圆的最小二乘拟合,计算出待测滑块滚道Roll1与滚道Roll2在空间斜坐标系下的滚道中心点的坐标(x4Roll1,y4Roll1,z4Roll1)、(x6Roll2,y6Roll2,z6Roll2),、将空间斜坐标系下的滚道Roll1、滚道Roll2中心点的坐标(x4Roll1,y4Roll1,z4Roll1)、(x6Roll2,y6Roll2,z6Roll2)转化为位移传感器***组成的空间直角坐标系下的滚道Roll1、滚道Roll2中心点的坐标(x2Roll1,y2Roll1,z2Roll1)、(x2Roll2,y2Roll2,z2Roll2);
步骤8、对待测滑块的不同截面进行测量,求解出不同截面的滚道Roll1、Roll2中心点在位移传感器***组成的空间直角坐标系下的坐标,拟合出滚道Roll1、Roll2中心轴线的方向向量;
步骤9、根据步骤8拟合出的滚道Roll1、Roll2中心轴线的方向向量,求解出滚道Roll1中心轴线与滚道Roll2中心轴线的距离,即为滚动直线导轨副滑块滚道中径。
本发明的有益效果是:
(1)本发明将激光位移传感器置于滑块外侧两端,并倾斜安装进行测量,避免了激光位移传感器检测期间与滑块干涉的问题;
(2)本发明巧妙运用了圆柱的界面形状为椭圆的原理,圆柱中心轴线上的点即为对应截面椭圆的圆心,有效地将滑块内滚道圆孔位置信息提取出来,观点新颖;
(3)本发明巧妙地利用标定圆柱作为检测基准来反求激光位移传感器地安装倾斜角,标定所用圆柱体容易加工到高精度以满足使用要求;
(4)本发明巧妙地将滑块中径测量问题转化为异面直线地距离,将不同坐标系下采集的数据通过几何关系与坐标变换转化到同一坐标系下,通过求解异面直线距离来求得滑块中径,形象直观便于计算。
除了上面所描述的目的、特征和优点之外,本发明还有其它的目的、特征和优点。下面将参照图,对本发明作进一步详细的说明。
附图说明
图1为本发明传感器的布局图。
图2为夹具空间直角坐标系、待测滑块空间直角坐标系、位移传感器***组成的空间直角坐标系、第一激光位移传感器和第二激光位移传感器的空间直角坐标系、第一激光位移传感器和第二激光位移传感器的空间斜坐标系图。
图3为坐标系o3-x3y3z3与坐标系o4-x4y4z4之间的关系图。
图4为坐标系o5-x5y5z5与坐标系o5-x5y5z5之间的关系图。
图5为本发明扫掠标定圆柱块所得出的半个椭圆的轮廓图。
图6为本发明根据扫掠出标定圆柱块所得半个椭圆拟合出的整个椭圆轮廓图。
图7为本发明扫掠滑块内滚道所得的两个半椭圆轮廓图。
图8为本发明根据扫掠滑块内滚道所得的两个半椭圆拟合成的两个椭圆轮廓图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明的一种滚动直线导轨副滑块型面中径的测量算法,包括以下步骤:
步骤1、构建坐标系,具体是构建夹具空间直角坐标系、待测滑块3空间直角坐标系、位移传感器***组成的空间直角坐标系、第一激光位移传感器2和第二激光位移传感器4的空间直角坐标系、第一激光位移传感器2和第二激光位移传感器4的空间斜坐标系;
构建夹具的空间直角坐标系,具体为夹具上固连的标定块10的空间直角坐标系o0-x0y0z0,其中x0轴和夹具台面垂直,y0轴沿着夹具台的长度方向,z0轴沿着夹具台的宽度方向,x0轴、y0轴、z0轴遵循右手法则;
构建待测滑块3的空间直角坐标系o1-x1y1z1,其中x1轴与待测滑块3的侧基准面垂直,y1轴方向与待测滑块3导向方向相同,z1轴垂直于待测滑块3的大基准面,x1轴、y1轴、z1轴遵循右手法则;
构建激光位移传感器***组成的空间直角坐标系o2-x2y2z2,其中y2轴方向与与y1轴方向一致,x2轴在竖直面内与y2轴垂直,指向向上,z2轴与x2轴、y2轴遵循右手法则;
构建第一激光位移传感器2的空间直角坐标系o3-x3y3z3,其中z3轴为第一激光位移传感器2运动方向,x3轴垂直于z3轴,且位于z3轴、激光射出线组成的面内,y3轴与z3轴、x3轴遵循右手法则;
构建第一激光位移传感器2的空间斜坐标系o4-x4y4z4,实际采集数据过程中,数据点在第一激光位移传感器2的空间斜坐标系o4-x4y4z4中采集,其中z4轴即z3轴,y4轴即y3轴,x4代表传感器光线射出方向,其与x3轴之间的夹角为安装偏转角β1
构建第二激光位移传感器4的空间直角坐标系o5-x5y5z5,其中z5轴为第二激光位移传感器4运动方向,x5轴垂直于z5轴,且位于z5轴、激光射出线组成的面内,y5轴与z5轴、x5轴遵循右手法则;
构建第二激光位移传感器4的空间斜坐标系o6-x6y6z6,实际采集数据过程中,数据点在第二激光位移传感器4的空间斜坐标系o6-x6y6z6中采集,z6轴即z5轴,y6轴即y5轴,x6代表传感器光线射出方向,其与x5轴之间的夹角为安装偏转角β2
所述激光位移传感器***中的第一激光位移传感器2、第二激光位移传感器4沿z2轴方向测量待测滑块3的两个滚道,滑块托架9上固连标定圆柱1、待测滑块3以及标定块10,测量过程中,滑块托架9带动标定圆柱1、待测滑块3以及标定块10沿y2轴方向移动,第三激光位移传感器5、第四激光位移传感器6、第五激光位移传感器7、第六激光位移传感器8固连在实验台上,其中第三激光位移传感器5、第四激光位移传感器6用于采集待测滑块3和标定块10侧基准面数据,第五激光位移传感器7、第六激光位移传感器8用于采集待测滑块3以及标定块10大基准面的数据;
通过构建七个空间坐标系,方便后续说明与计算。
步骤2、利用标定块10得到夹具空间直角坐标系下的***误差曲线:第三激光位移传感器5与第四激光位移传感器6采集标定块10的侧基准面数据,第五激光位移传感器7与第六激光位移传感器8采集标定块10的大基准面数据,滑块托架9沿y2轴运动的距离由光栅采集。分别以光栅采集的数据为横轴,以第三激光位移传感器5、第四激光位移传感器6采集的数据为纵轴,得到标定块10侧基准面的两条曲线数据;分别以光栅采集的数据为横轴,以第五激光位移传感器7、第六激光位移传感器8采集的数据为纵轴,得到标定块10大基准面的两条曲线数据,获取的四条曲线作为夹具***的误差曲线,用于后续补偿;
滑块托架9沿着y2轴方向运动过程中,第三激光位移传感器5、第四激光位移传感器6、第五激光位移传感器7、第六激光位移传感器8以及光栅测量出一系列坐标值(0x2,0y2,0z2),对于第五激光位移传感器7、第六激光位移传感器8,0x2值由测距仪测得,实验过程中第五激光位移传感器7、第六激光位移传感器8与标定块10在x2轴方向上没有数值变化,为定值,0y2值由光栅尺读取,0z2值由激光位移传感器实时采集。对于第三激光位移传感器5、第四激光位移传感器6,0z2值由测距仪测得,实验过程中第三激光位移传感器5、第四激光位移传感器6与标定块10在z2轴方向上没有数值变化,为定值,0y2值由光栅尺读取,0x2值由激光位移传感器实时采集。通过测量标定块10的两个与夹具配合的面,通过测量标定块两个面上具有代表性的数据来反映对应平面特性的原理,获取夹具空间直角坐标系下的***误差曲线,用于后续补偿。
步骤3、确定待测滑块3大基准面与侧基准面的平面度:第三激光位移传感器5、第四激光位移传感器6采集待测滑块3的侧基准面数据,第五激光位移传感器7、第六激光位移传感器8采集待测滑块3的大基准面数据,滑块托架9沿y2轴运动的距离由光栅采集。分别以光栅采集的数据为横轴,以第三激光位移传感器5、第四激光位移传感器6采集的数据为纵轴,得到待测滑块3侧基准面的两条曲线数据;分别以光栅采集的数据为横轴,以第五激光位移传感器7、第六激光位移传感器8采集的数据为纵轴,得到待测滑块3大基准面的两条曲线数据,并将步骤2中得到的夹具***误差曲线作为补偿曲线,得到待测滑块3大基准面、侧基准面在夹具空间直角坐标系下的坐标值,利用平面的最小二乘拟合方程,求出待测滑块3的大基准面和侧基准面的平面度;
同步骤2,第三激光位移传感器5、第四激光位移传感器6、第五激光位移传感器7、第六激光位移传感器8在沿着y2轴方向运动过程中,测量待测滑块3的大基准面、侧基准面,得到一系列坐标值(1x2,1y2,1z2),对于第三激光位移传感器5、第四激光位移传感器6,1z2值与步骤2中测得的0z2值相等,为定值,1y2由光栅尺读取,1x2值由第三激光位移传感器5、第四激光位移传感器6实时采集;对于第五激光位移传感器7、第六激光位移传感器8,1x2值与步骤2中测得的0x2值相等,为定值,1y2由光栅尺读取,1z2值由第五激光位移传感器7、第六激光位移传感器8实时采集。将步骤2中得到夹具空间直角坐标系下的***误差曲线用于补偿,得到待测滑块3的大基准面、侧基准面相对夹具空间直角坐标系的坐标:
(1x0,1y0,1z0)=((1x2-0x2),(1y2-0y2),(1z2-0z2)) (1)
对激光位移传感器7、8所测量的数据。建立大基准面的最小二乘拟合方程为:
Figure GDA0003211546310000081
式中(1x0,1y0,1z0)由第五激光位移传感器7、第六激光位移传感器8所测得,i=1,2...n,n为采集的点的数量;
对方程(2)系数A1、B1、C1求偏导得到方程组为:
Figure GDA0003211546310000082
求解式(3)中的方程组,结果为:
Figure GDA0003211546310000091
由式(4)求得方程系数A1,B1,C1,得到大基准面平面度最小二乘法评定结果t1为:
Figure GDA0003211546310000092
第三激光位移传感器5、第四激光位移传感器6所测量的面,建立侧基准面的最小二乘拟合方程为:
Figure GDA0003211546310000093
式中(1x0,1y0,1z0)由第三激光位移传感器5、第四激光位移传感器6所测得,i=1,2...n,n为采集的点的数量;
同理通过式(3)、式(4)、式(5)得侧基准面平面度的最小二乘法评定结果t2为:
Figure GDA0003211546310000094
步骤4、通过第一激光位移传感器2、第二激光位移传感器4扫掠标定圆柱1,利用椭圆的最小二乘法拟合,求解出椭圆一般方程的系数。
当第一激光位移传感器2测量标定圆柱1时,可得一系列轮廓数据(3x4,3y4,3z4)构建椭圆的最小二乘拟合方程为:
Figure GDA0003211546310000095
其中(3x4,3y4,3z4)为第一激光位移传感器2以及光栅,测出的标定圆柱1的轮廓数据,i=1,2...n,n为采集的点的数量;其中3x4i值由第一激光位移传感器2所采集,3z4i值由光栅所采集,3y4值由测距仪测得,测量过程中第一激光位移传感器2与标定圆柱1在y4方向上没有数值变化,为定值。
对式(8)中方程系数A0,B0,C0,D0,E0求偏导数为:
Figure GDA0003211546310000101
展开得:
Figure GDA0003211546310000102
对式(10)中的方程组求解,得到标定圆柱1所拟合的椭圆方程系数A0,B0,C0,D0,E0
步骤5、构建步骤4中所得椭圆方程与第一激光位移传感器2、第二激光位移传感器4安装倾斜角β1、β2。安装倾斜角β1是在第一激光位移传感器2射出光线与第一激光位移传感器2运动方向构成的平面内,第一激光位移传感器2运动方向的垂直线与第一激光位移传感器2射出光线之间的夹角;安装倾斜角β2是第二激光位移传感器4射出光线与第二激光位移传感器4运动方向构成的平面内,第二激光位移传感器4运动方向的垂直线与第二激光位移传感器4射出光线之间的夹角。求解出倾斜安装的第一激光位移传感器2和第二激光位移传感器4实际的安装倾斜角β1、β2
对于第一激光位移传感器2,空间斜坐标系(x4,y4,z4)下实际坐标点与空间直角坐标系下的理论坐标点(x3,y3,z3)存在如下关系:
Figure GDA0003211546310000103
对于任意位置的椭圆,可以根据5个独立的参数描述:中心点(x0,z0)、长半轴a、短半轴b(假设a>b)、倾角θ,在空间直角坐标系o4-x4y4z4中任意位置的椭圆方程式表达为:
Figure GDA0003211546310000111
由于存在β1角产生的偏差,将斜坐标系o3-x3y3z3与直角坐标系o4-x4y4z4之间坐标转化关系带入,得到实际椭圆方程式:
Figure GDA0003211546310000112
式(13)展开是一个五元四次非线性方程,用下述变量代换方法将此复杂的非线性方程变换乘线性方程,即A1x4 2+B1x4z4+C1z4 2+D1x4+E1y4+F1=0,其中
A1=cos2θ/a2+sin2θ/b2
B1=2cosθ(cosβ1sinθ-cosθsinβ1)/a2-2sinθ(sinβ1sinθ+cosβ1cosθ)/b2
C1=(cosβ1sinθ-cosθsinβ1)2/a2+(sinβ1sinθ+cosβ1cosθ)2/b2
D1=(2sinθ(y0cosθ-x0sinθ)/b2-2cosθ(x0cosθ+y0sinθ)/a2
E1=-2(cosβ1sinθ-cosθsinβ1)(x0cosθ+y0sinθ)/a2-2(sinβ1sinθ+cosβ1cosθ)(y0cosθ-x0sinθ)/b2
F1=(x0cosθ+y0sinθ)2/a2+(y0cosθ-x0sinθ)2/b2-1
(14)
Figure GDA0003211546310000113
则在斜坐标系下的椭圆方程转化为:
A4x4 2+B4x4z4+C4z4 2+D4x4+E4y4+1=0 (16)
令式(8)构建椭圆的最小二乘拟合方程系数与式(16)斜坐标系下椭圆方程系数相等,组成方程组:
Figure GDA0003211546310000121
方程组(17)中,椭圆长轴2a、中心点(x0,z0)、倾角θ以及安装倾斜角β1为未知参数,短轴2b为标定圆柱1直径为已知,五个未知数联立五个方程组,计算求出第一激光位移传感器2的安装倾斜角β1的值;
对于第二激光位移传感器4,同理根据式(11)至式(17)可以求出第二激光位移传感器4的安装倾斜角β2的值。
步骤6、确定第一激光位移传感器2、第二激光位移传感器4绕自身运动方向轴旋转的安装偏转角α1、α2。安装偏转角α1是x3轴方向与滑块托架9运动方向y2轴之间的夹角;安装偏转角α2是x5轴方向与滑块托架9运动方向y2轴之间的夹角。求解出第一激光位移传感器2和第二激光位移传感器4实际的安装偏转角α1、α2
滑块托架9沿y2方向运动ΔL,距离由光栅尺读出。第一激光位移传感器2前后数值变化量为Δd1,第二激光位移传感器4前后数值变化量为Δd2。第一激光位移传感器2、第二激光位移传感器4实际的安装偏转角为
Figure GDA0003211546310000122
步骤7、计算待测滑块3内滚道中心点坐标,通过第一激光位移传感器2和第二激光位移传感器4扫掠待测滑块3内滚道,采集分别位于第一激光位移传感器2和第二激光位移传感器4空间斜坐标系下的待测滑块3滚道数据点,采用椭圆的最小二乘拟合,计算出待测滑块3滚道Roll1与滚道Roll2在空间斜坐标系下的滚道Roll1中心点坐标、滚道Roll2中心点坐标(x4Roll1,y4Roll1,z4Roll1)、(x6Roll2,y6Roll2,z6Roll2),、将空间斜坐标系下的滚道Roll1中心点坐标、滚道Roll2中心点坐标(x4Roll1,y4Roll1,z4Roll1)、(x6Roll2,y6Roll2,z6Roll2)转化为位移传感器***组成的空间直角坐标系下的椭圆中心点的坐标(x2Roll1,y2Roll1,z2Roll1)、(x2Roll2,y2Roll2,z2Roll2);
当第一激光位移传感器2沿着z4轴方向测量待测滑块3滚道内轮廓时,测量方向上扫掠所得的内滚道轮廓为一半椭圆,可得一系列内滚道数据(1x4,1y4,1z4)。
构建椭圆的最小二乘拟合方程为:
Figure GDA0003211546310000131
其中(1x4i,1y4i,1z4i)为第一激光位移传感器2以及光栅测出的滚道Roll1的轮廓数据,i=1,2...n,n为采集的点的数量。
对式(19)中系数A,B,C,D,E求偏导数为:
Figure GDA0003211546310000132
式(20)展开得:
Figure GDA0003211546310000133
对式(21)求解,得到A,B,C,D,E;
对一般椭圆方程:
x4Roll1 2+Ax4Roll1z4Roll1+Bz4Roll1 2+Cx4Roll1+Dz4Roll1+E=0 (22)
公式(22)一般椭圆方程所表示椭圆中心点的坐标为第一激光位移传感器2空间斜坐标系下的滚道Roll1中心点的坐标,表示为:
Figure GDA0003211546310000134
将第一激光位移传感器2空间斜坐标系下的滚道Roll1中心点的坐标(x4Roll1,y4Roll1,z4Roll1)转化为第一激光位移传感器2空间直角坐标系下的滚道Roll1中心点的坐标(x3Roll1,y3Roll1,z3Roll1),二者存在如下关系:
Figure GDA0003211546310000141
将第一激光位移传感器2空间直角坐标系下的滚道Roll1中心点的坐标(x3Roll1,y3Roll1,z3Roll1)转化为位移传感器***组成的空间直角坐标系下的滚道Roll1中心点的坐标(x2Roll1,y2Roll1,z2Roll1),二者存在如下关系:
Figure GDA0003211546310000142
其中,(x20Roll1,y20Roll1,z20Roll1)为第一激光位移传感器2直角坐标系原点在位移传感器***组成空间直角坐标系下的坐标值。x20Roll1由电机编码器读取,y20Roll1、z20Roll1由光栅尺读取。
至此,完成拟合出的滚道Roll1椭圆中心点坐标从第一激光位移传感器2空间斜坐标系(x4Roll1,y4Roll1,z4Roll1)到激光位移传感器***组成空间直角坐标系(x2Roll1,y2Roll1,z2Roll1)的转换。
对于第二激光位移传感器4,同理运用式(19)至式(25),可以得到待测滑块3另一滚道Roll2某一截面中心点在位移传感器***组成空间直角坐标系下的滚道Roll2中心的坐标点(x2Roll2,y2Roll2,z2Roll2)。
步骤8、根据步骤7,对待测滑块3的不同截面进行测量,求解出不同截面的滚道Roll1,Roll2中心点在位移传感器***组成的空间直角坐标系下的坐标。拟合出滚道Roll1,Roll2中心轴线的直线方程。
对于滚道Roll1,待测滑块3固定在滑块托架9上,移动滑块托架9,当第一激光位移传感器2刚好能够采集到待测滑块3内滚道数据时记为位置1,通过倾斜安装的激光位移传感器2扫掠滑块内滚道。将扫掠到的待测滑块3内滚道数据点拟合出中心点坐标,并将中心点坐标从第一激光位移传感器2空间斜坐标系o4-x4y4z4变换到位移传感器***组成的空间直角坐标系o2-x2y2z2,得到所述位置1下,滚道Roll1截面的拟合椭圆中心点在位移传感器***组成空间直角坐标系o2-x2y2z2下的坐标(1x2Roll1,1y2Roll1,1z2Roll1);
移动滑块托架使之处于位置k,联系步骤8-1,可以测出移动滑块托架处于不同位置时,每个滑块托架位置之间的距离为1cm,距离由伺服电机控制。测出移动滑块托架9处于不同位置时,滚道Roll1截面的拟合椭圆中心点在位移传感器***组成空间直角坐标系o2-x2y2z2下的坐标(kx2Roll1,ky2Roll1,kz2Roll1),k=1,2...10。
滑块托架9处于不同位置所得坐标点,拟合出滚道Roll1中心轴线的直线方程。由于滚道Roll1中心轴线穿过x2o2z2面,则空间直线方程可以简化为:
Figure GDA0003211546310000151
其中(x01,y01,z01)为滚道Roll1中心轴线上的任意一点,(mRoll1,1,nRoll1)为滚道Roll1中心轴线的法向向量,滚道Roll1径向截面为椭圆弧形状,滚道Roll1中心轴线表示由滚道Roll1不同径向截面的椭圆弧中心组成的直线,y01=0x01,z01,mRoll1,nRoll1为要求参数。
直线方程可以简化成:
Figure GDA0003211546310000152
化为矩阵形式:
Figure GDA0003211546310000153
第k个点方程为:
Figure GDA0003211546310000154
并联10个点的方程为:
Figure GDA0003211546310000161
最小二乘拟合:
Figure GDA0003211546310000162
化简成:
Figure GDA0003211546310000163
求解矩阵即可得到mRoll1,nRoll1,x01,y01,为直线方程的各个参数。得到滚道Roll1中心轴线的方向向量为
Figure GDA0003211546310000164
对于滚道Roll2中心轴线,同理运动式(26)至式(32),可以得到滚道Roll2中心轴线的方向向量为
Figure GDA0003211546310000165
步骤9、根据步骤8拟合出的滚道Roll1,Roll2中心轴线的方向向量,求解出滚道Roll1中心轴线与滚道Roll2中心轴线的距离,即为滚动直线导轨副滑块滚道中径。
根据步骤8拟合出的滚道Roll1,Roll2中心轴线的方向向量,求解滚道Roll1与滚道Roll2中心轴线的距离,即为滑块滚道中径,两中心轴线距离为:
Figure GDA0003211546310000166
其中
Figure GDA0003211546310000167
为滚道Roll1与滚道Roll2中心轴线的方向向量,
Figure GDA0003211546310000168
Figure GDA0003211546310000169
是两条直线上任意一点的连线方向,
Figure GDA00032115463100001610
本发明通过将高精度激光位移传感器置于滑块端部两侧,并对滑块滚道进行倾斜扫掠,扫掠的结果数据点为半个椭圆形,对采集到的数据采用最小二乘法进行拟合。求解出拟合椭圆的形心坐标,拟合椭圆的形心坐标即为检测滑块内滚道的形心。通过测量已知尺寸和安装位置坐标的标定圆柱,可以将激光位移传感器的实际安装倾斜角标定出来;通过测量已知安装位置的标定块,可以将激光位移传感器的实际安装偏转角标定出来;在考虑激光位移传感器倾斜角与偏转角的基础上,通过检测拟合滑块不同截面的内滚道形心点坐标,将同一滚道的一系列形心点坐标拟合出中心轴的直线方程,求解滑块两侧对称的内滚道中心直线方程的异面直线距离即为滑块的中径;通过另外两对激光位移传感器将滑块的大基准面、侧基准面进行检测,提出了检测滚动直线导轨副滑块型面中径的算法。
下面结合实施例对本发明做进一步详细的描述。
实施例1:
滚动直线导轨副滑块的型面中径检测结合说明书附图中的图1、图2、图3、图4来进行说明。标定块10与标定圆柱块1固连在滑块托架9上面,所用标定圆柱的标准尺寸为9.996mm,圆度为0.003mm,直线度为0.003mm。第三激光位移传感器5、第四激光位移传感器6、第五激光位移传感器7、第六激光位移传感器8组成的滑块基准面测量***固定在试验台上,滑块托架9及固连在滑块托架9上的标定块10、标定圆柱1、待测滑块3同时沿着y2轴运动,第一激光位移传感器2、第二激光位移传感器4组成的滑块内滚道测量***能够沿着z2轴与x2轴运动,即第一激光位移传感器2、第二激光位移传感器4组成的滑块内滚道测量***能够相对第三激光位移传感器5、第四激光位移传感器6、第五激光位移传感器7、第六激光位移传感器8组成的滑块基准面测量***沿着z2轴运动。
首先滑块托架9带着固定在上面的标定块10,沿着y2轴运动,第三激光位移传感器5、第四激光位移传感器6、第五激光位移传感器7、第六激光位移传感器8组成的滑块基准面测量***测量标定块10与滑块托架9配合的面作为补偿曲线。接着第三激光位移传感器5、第四激光位移传感器6、第五激光位移传感器7、第六激光位移传感器8组成的滑块基准面测量***测量待测滑块3与滑块托架9配合的面,将上一步采集到的误差曲线用于滑块大基准面与侧基准面的数据补偿,可以得到滑块大基准面与侧基准面在夹具坐标系o0-x0y0z0下的坐标点。运用式(1)至式(7)可以求出待测滑块大基准面与侧基准面的平行度。
第一激光位移传感器2、第二激光位移传感器4组合的滑块内滚道测量***沿着z2轴运动,采集到标定圆柱1的外轮廓数据点坐标,如图5所示;将采集到的标定圆柱1数据点拟合成椭圆,如图6所示;标定圆柱1的实际直径为9.996mm,联系式(8)至式(11)可以求解出第一激光位移传感器2、第二激光位移传感器4的安装倾斜角为:
β1=0.1321°,β2=0.1569°
滑块托架9沿着z2轴方向移动,第一激光位移传感器2、第二激光位移传感器4采集到标定块10之间的距离,采用式(18)可以求解出第一激光位移传感器2、第二激光位移传感器4的安装偏转角为:
α1=119.561°,α2=60.964°
第一激光位移传感器2、第二激光位移传感器4组成的滑块内滚道测量***沿着z2轴运动,采集出待测滑块3的内滚道坐标顶,如图7所示;运用式(19)至式(23)拟合出被截斜面轮廓形心在激光位移传感器2、4空间斜坐标系下的坐标值为:
(x4Roll1,y4Roll1,z4Roll1)=(5.135,0,5.896),(x6Roll2,y6Roll2,z6Roll2)=(5.231,0,5.896),其中Roll1,Roll2表示滚道Roll1与滚道Roll2。图8为内滚道拟合全椭圆图。
改变滑块托架9位置,令第一激光位移传感器2、第二激光位移传感器4沿z2轴方向扫掠滑块内滚道,将所得到数据拟合出内滚道形心坐标点。如此往复10次,然后运用式(24)、式(25)将得到的在第一激光位移传感器2、第二激光位移传感器4空间斜坐标系下的形心坐标点转化到位移传感器***组成的空间直角坐标系下的中心点坐标(x2Roll1,y2Roll1,z2Roll1),(x2Roll2,y2Roll2,z2Roll2)。
运用式(26)至式(32)拟合出内滚道Roll1与内滚道Roll2的中心轴线方程。运用式(33)求解出滑块中径为:d=45.324mm。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (9)

1.一种滚动直线导轨副滑块型面中径的测量算法,其特征在于,包括以下步骤:
步骤1、构建坐标系:构建夹具空间直角坐标系、待测滑块(3)的空间直角坐标系、激光位移传感器***组成的空间直角坐标系、第一激光位移传感器(2)和第二激光位移传感器(4)的空间直角坐标系、第一激光位移传感器(2)和第二激光位移传感器(4)的空间斜坐标系;
构建夹具的空间直角坐标系,具体为夹具上固连的标定块(10)的空间直角坐标系o0-x0y0z0,其中x0轴和夹具台面垂直,y0轴沿着夹具台的长度方向,z0轴沿着夹具台的宽度方向,x0轴、y0轴、z0轴遵循右手法则;
构建待测滑块(3)的空间直角坐标系o1-x1y1z1,其中x1轴与待测滑块(3)的侧基准面垂直,y1轴方向与待测滑块(3)导向方向相同,z1轴垂直于待测滑块(3)的大基准面,x1轴、y1轴、z1轴遵循右手法则;
构建激光位移传感器***组成的空间直角坐标系o2-x2y2z2,其中y2轴方向与y1轴方向一致,x2轴在竖直面内与y2轴垂直,指向向上,z2轴与x2轴、y2轴遵循右手法则;
构建第一激光位移传感器(2)的空间直角坐标系o3-x3y3z3,其中z3轴为第一激光位移传感器(2)运动方向,x3轴垂直于z3轴,且位于z3轴、激光射出线组成的面内,y3轴与z3轴、x3轴遵循右手法则;
构建第一激光位移传感器(2)的空间斜坐标系o4-x4y4z4,实际采集数据过程中,数据点在第一激光位移传感器(2)的空间斜坐标系o4-x4y4z4中采集,其中z4轴即z3轴,y4轴即y3轴,x4轴代表传感器光线射出方向,其与x3轴之间的夹角为安装倾斜角β1
构建第二激光位移传感器(4)的空间直角坐标系o5-x5y5z5,其中z5轴为第二激光位移传感器(4)运动方向,x5轴垂直于z5轴,且位于z5轴、激光射出线组成的面内,y5轴与z5轴、x5轴遵循右手法则;
构建第二激光位移传感器(4)的空间斜坐标系o6-x6y6z6,实际采集数据过程中,数据点在第二激光位移传感器(4)的空间斜坐标系o6-x6y6z6中采集,z6轴即z5轴,y6轴即y5轴,x6代表传感器光线射出方向,其与x5轴之间的夹角为安装倾斜角β2
所述激光位移传感器***中的第一激光位移传感器(2)、第二激光位移传感器(4)沿z2轴方向测量待测滑块(3)的两个滚道,滑块托架(9)上固连标定圆柱(1)、待测滑块(3)以及标定块(10),测量过程中,滑块托架(9)带动标定圆柱(1)、待测滑块(3)以及标定块(10)沿y2轴方向移动,第三激光位移传感器(5)、第四激光位移传感器(6)、第五激光位移传感器(7)、第六激光位移传感器(8)固连在实验台上,其中第三激光位移传感器(5)、第四激光位移传感器(6)用于采集待测滑块(3)和标定块(10)侧基准面数据,第五激光位移传感器(7)、第六激光位移传感器(8)用于采集待测滑块(3)以及标定块(10)大基准面的数据;
步骤2、利用标定块(10)得到夹具空间直角坐标系下的***误差曲线:侧基准面的第三激光位移传感器(5)、第四激光位移传感器(6)与大基准面的第五激光位移传感器(7)、第六激光位移传感器(8)测量标定块(10)的两个与夹具配合的面,得到大基准面与侧基准面各两条直线数据,作为夹具***的误差曲线,用于后续补偿;
步骤3、确定待测滑块(3)大基准面与侧基准面的平面度:第三激光位移传感器(5)、第四激光位移传感器(6)、第五激光位移传感器(7)、第六激光位移传感器(8)测量待测滑块(3)的大基准面与侧基准面,得到待测滑块(3)大基准面与侧基准面各两条直线数据,并将步骤2中得到的夹具***误差曲线作为补偿曲线,得到待测滑块(3)大基准面、侧基准面在夹具空间直角坐标系下的坐标值,利用平面的最小二乘拟合方程,求出待测滑块(3)的大基准面和侧基准面的平面度;
步骤4、第一激光位移传感器(2)和第二激光位移传感器(4)扫掠标定圆柱(1),利用椭圆的最小二乘法拟合,求解出标定椭圆一般方程;
步骤5、构建步骤4中所得椭圆一般方程与第一激光位移传感器(2)和第二激光位移传感器(4)安装倾斜角β1、β2并求解出所述安装倾斜角β1、β2:所述安装倾斜角β1是在第一激光位移传感器(2)射出光线与第一激光位移传感器(2)运动方向构成的平面内第一激光位移传感器(2)运动方向的垂直线与第一激光位移传感器(2)射出光线之间的夹角;安装倾斜角β2是第二激光位移传感器(4)射出光线与第二激光位移传感器(4)运动方向构成的平面内第二激光位移传感器(4)运动方向的垂直线与第二激光位移传感器(4)射出光线之间的夹角;
步骤6、确定第一激光位移传感器(2)和第二激光位移传感器(4)绕自身运动方向轴旋转的安装偏转角α1、α2并求解出所述安装偏转角α1、α2;所述安装偏转角α1是x3轴方向与滑块托架(9)运动方向y2轴之间的夹角;安装偏转角α2是x5轴方向与滑块托架(9)运动方向y2轴之间的夹角;
步骤7、计算待测滑块(3)内滚道中心点坐标:第一激光位移传感器(2)和第二激光位移传感器(4)扫掠待测滑块(3)内滚道,采集分别位于第一激光位移传感器(2)和第二激光位移传感器(4)空间斜坐标系下的待测滑块(3)滚道数据点,采用椭圆的最小二乘拟合,计算出待测滑块(3)滚道ROLL1与滚道ROLL2在空间斜坐标系下的滚道中心点的坐标(x4Roll1,y4Roll1,z4Roll1)、(x6Roll2,y6Roll2,z6Roll2),、将空间斜坐标系下的滚道ROLL1、滚道ROLL2中心点的坐标(x4Roll1,y4Roll1,z4Roll1)、(x6Roll2,y6Roll2,z6Roll2)转化为位移传感器***组成的空间直角坐标系下的滚道ROLL1、滚道ROLL2中心点的坐标(x2Roll1,y2Roll1,z2Roll1)、(x2Roll2,y2Roll2,z2Roll2);
步骤8、对待测滑块(3)的不同截面进行测量,求解出不同截面的滚道ROLL1、ROLL2中心点在位移传感器***组成的空间直角坐标系下的坐标,拟合出滚道ROLL1、ROLL2中心轴线的方向向量;
步骤9、根据步骤8拟合出的滚道ROLL1、ROLL2中心轴线的方向向量,求解出滚道ROLL1中心轴线与滚道ROLL2中心轴线的距离,即为滚动直线导轨副滑块滚道中径。
2.根据权利要求1所述的滚动直线导轨副滑块型面中径测量算法,其特征在于,所述步骤2获取夹具空间直角坐标系的***误差曲线具体包括:滑块托盘(9)沿着y2轴方向运动过程中,第三激光位移传感器(5)、第四激光位移传感器(6)、第五激光位移传感器(7)、第六激光位移传感器(8)以及光栅测量出一系列坐标值(0x2,0y2,0z2),对于第五激光位移传感器(7)、第六激光位移传感器(8),0x2值在实验前由测距仪测得,为定值,0y2值由光栅尺读取,0z2值由第五激光位移传感器(7)、第六激光位移传感器(8)实时采集,对于第三激光位移传感器(5)、第四激光位移传感器(6),0z2值由测距仪测得,为定值,0y2值由光栅尺读取,0x2值由第三激光位移传感器(5)、第四激光位移传感器(6)实时采集,分别以光栅采集的数据为横轴,以第三激光位移传感器(5)、第四激光位移传感器(6)采集的数据为纵轴,得到标定块(10)侧基准面的两条曲线数据,分别以光栅采集的数据为横轴,以第五激光位移传感器(7)、第六激光位移传感器(8)采集的数据为纵轴,得到标定块(10)大基准面的两条曲线数据,获取的四条曲线作为夹具***的误差曲线,用于后续补偿。
3.根据权利要求2所述的滚动直线导轨副滑块型面中径测量算法,其特征在于,步骤3确定待测滑块(3)大基准面、侧基准面的平面度具体包括:
步骤3-1:滑块托盘(9)在沿着y2轴方向运动过程中,第三激光位移传感器(5)、第四激光位移传感器(6)、第五激光位移传感器(7)、第六激光位移传感器(8)以及光栅测量待测滑块(3)的大基准面、侧基准面,得到一系列坐标值(1x2,1y2,1z2),对于第三激光位移传感器(5)、第四激光位移传感器(6),1z2值与0z2值相等,为定值,1y2值由光栅尺读取,1x2值由第三激光位移传感器(5)、第四激光位移传感器(6)实时采集;对于第五激光位移传感器(7)、第六激光位移传感器(8),1x2值与0x2值相等,为定值,1y2由光栅尺读取,1z2值由第五激光位移传感器(7)、第六激光位移传感器(8)实时采集,分别以光栅采集的数据为横轴,以第三激光位移传感器(5)、第四激光位移传感器(6)采集的数据为纵轴,得到待测滑块(3)侧基准面的两条曲线数据,分别以光栅采集的数据为横轴,以第五激光位移传感器(7)、第六激光位移传感器(8)采集的数据为纵轴,得到待测滑块(3)大基准面的两条曲线数据,将步骤2中得到夹具空间直角坐标系下的***误差曲线用于补偿,得到待测滑块(3)的大基准面、侧基准面相对夹具空间直角坐标系的坐标:
(1x0, 1y0,1z0)=((1x2-0x2),(1y2-0y2),(1z2-0z2)) (1);
步骤3-2:对第五激光位移传感器(7)、第六激光位移传感器(8)所测量的数据,建立大基准面的最小二乘拟合方程为:
Figure FDA0003252026630000041
上式中(1x0,1y0,1z0)由步骤3-1误差补偿后的数据点,i=1,2...n,N为采集的点的数量;
步骤3-3、联系步骤3-2中的方程,对方程系数A1、B1、C1求偏导得到方程组为:
Figure FDA0003252026630000042
步骤3-4:求解步骤3-3中的方程组,结果为:
Figure FDA0003252026630000043
步骤3-5:由步骤3-4求得方程系数A1,B1,C1,得到大基准面平面度最小二乘法评定结果t1为:
Figure FDA0003252026630000051
步骤3-6:第三激光位移传感器(5)、第四激光位移传感器(6)所测量的面,建立侧基准面的最小二乘拟合方程为:
Figure FDA0003252026630000052
同理通过步骤3-3、3-4得侧基准面平面度的最小二乘法评定结果t2为:
Figure FDA0003252026630000053
4.根据权利要求3所述的滚动直线导轨副滑块型面中径测量算法,其特征在于,所述步骤4具体包括:
步骤4-1:当第一激光位移传感器(2)测量标定圆柱(1)时,得到一系列轮廓数据(3x4,3y4,3z4);
步骤4-2:构建标定椭圆的最小二乘拟合方程为:
Figure FDA0003252026630000054
其中(3x4,3y4,3z4)为第一激光位移传感器(2)以及光栅在o4-x4y4z4坐标系下测出的标定圆柱(1)的轮廓数据,i=1,2...n,N为采集的点的数量,其中3x4i值由第一激光位移传感器(2)所采集,3z4i值由光栅所采集,3y4值由测距仪测得,为定值;
步骤4-3:对步骤4-2中方程求偏导数为:
Figure FDA0003252026630000055
展开得:
Figure FDA0003252026630000061
步骤4-4:对步骤4-3中的方程组求解,得到标定圆柱(1)所拟合的椭圆方程系数A0,B0,C0,D0,E0
5.根据权利要求4所述的滚动直线导轨副滑块型面中径测量算法,其特征在于,所述步骤5具体包括:
步骤5-1:对于第一激光位移传感器(2),空间斜坐标系(x4,y4,z4)下实际坐标点与空间直角坐标系下的理论坐标点(x3,y3,z3)存在如下关系:
x3=x4-z4sinβ1
y3=y4
z3=z4cosβ1 (11)
步骤5-2:对于任意位置的椭圆,根据5个独立的参数:中心点(x0,z0)、长半轴A、短半轴B、倾角θ,在空间直角坐标系o4-x4y4z4中任意位置的椭圆方程式表达为:
Figure FDA0003252026630000062
由于存在β1角产生的偏差,将空间斜坐标系o3-x3y3z3与空间直角坐标系o4-x4y4z4之间坐标转化关系带入,得到实际椭圆方程式:
Figure FDA0003252026630000063
上式(13)展开是一个五元四次非线性方程,用下述变量代换方法将此复杂的非线性方程变换乘线性方程,即A1x4 2+B1x4z4+C1z4 2+D1x4+E1y4+F1=0,其中
A1=cos2θ/a2+sin2θ/b2
B1=2cosθ(cosβ1sinθ-cosθsinβ1)/a2-2sinθ(sinβ1sinθ+cosβ1cosθ)/b2
C1=(cosβ1sinθ-cosθsinβ1)2/a2+(sinβ1sinθ+cosβ1cosθ)2/b2
D1=(2sinθ(y0cosθ-x0sinθ)/b2-2cosθ(x0cosθ+y0sinθ)/a2
E1=-2(cosβ1sinθ-cosθsinβ1)(x0cosθ+y0sinθ)/a2-2(sinβ1sinθ+cosβ1cosθ)(y0cosθ-x0sinθ)/b2
F1=(x0cosθ+y0sinθ)2/a2+(y0cosθ-x0sinθ)2/b2-1 (14)
A4=A1/F1
B4=B1/F1
C4=C1/F1
D4=D1/F1
E4=E1/F1 (15)
则在空间斜坐标系下的椭圆方程转化为:
A4x4 2+B4x4z4+C4z4 2+D4x4+E4y4+1=0 (16);
步骤5-3:令构建椭圆的最小二乘拟合方程系数与空间斜坐标系下椭圆方程系数相等,组成方程组:
A4=A0
B4=B0
C4=C0
D4=D0
E4=E0 (17)
方程组(17)中,椭圆长轴2a、中心点(x0,z0)、倾角θ以及安装倾斜角β1为未知参数,短轴2b为标定圆柱(1)直径为已知,五个未知数联立五个方程组,计算求出第一激光位移传感器(2)的安装倾斜角β1的值;
步骤5-4:对于第二激光位移传感器(4),同样采用步骤5-1、5-2、5-3求出第二激光位移传感器(4)的安装倾斜角β2的值。
6.根据权利要求5所述的滚动直线导轨副滑块型面中径测量算法,其特征在于,所述步骤6具体包括:
步骤6-1:滑块托架(9)沿y2方向运动ΔL,距离由光栅尺读出,第一激光位移传感器(2)前后数值变化量为Δd1,第二激光位移传感器(4)前后数值变化量为Δd2
步骤6-2:第一激光位移传感器(2)、第二激光位移传感器(4)实际的安装偏转角为:
Figure FDA0003252026630000081
Figure FDA0003252026630000082
7.根据权利要求6所述的滚动直线导轨副滑块型面中径测量算法,其特征在于,所述步骤7具体包括:
步骤7-1:当第一激光位移传感器(2)沿着z4轴方向测量待测滑块(3)滚道内轮廓时,测量方向上扫掠所得的内滚道轮廓为一半椭圆,得到一系列内滚道数据(1x4,1y4,1z4);
步骤7-2:构建椭圆的最小二乘拟合方程为:
Figure FDA0003252026630000083
其中(1x4i,1y4i,1z4i)为第一激光位移传感器(2)以及光栅测出的滚道ROLL1的轮廓数据,i=1,2...n,N为采集的点的数量;
步骤7-3:对式(19)中的系数A,B,C,D,E求偏导数为:
Figure FDA0003252026630000084
式(20)展开得:
Figure FDA0003252026630000085
步骤7-4:对式(21)求解,得到A,B,C,D,E;
步骤7-5:对一般椭圆方程:x4Roll1 2+Ax4Roll1z4Roll1+Bz4Roll1 2+Cx4Roll1+Dz4Roll1+E=0(22),
公式(22)所表示椭圆中心点的坐标为第一激光位移传感器(2)空间斜坐标系下的滚道ROLL1中心点的坐标,表示为:
Figure FDA0003252026630000091
步骤7-6:将第一激光位移传感器(2)空间斜坐标系下的滚道ROLL1中心点坐标(x4Roll1,y4Roll1,z4Roll1)转化为第一激光位移传感器(2)空间直角坐标系下的滚道ROLL2中心点的坐标(x3Roll1,y3Roll1,z3Roll1),二者存在如下关系:
x3Roll1=x4Roll1-z4Roll1sinβ1
y3Roll1=y4Roll1=0
z3Roll1=z4Roll1cosβ1 (24)
步骤7-7:将第一激光位移传感器(2)空间直角坐标系下的滚道ROLL1中心点坐标(x3Roll1,y3Roll1,z3Roll1)转化为位移传感器***组成的空间直角坐标系下的滚道ROLL1中心点坐标(x2Roll1,y2Roll1,z2Roll1),二者存在如下关系:
x2Roll1=x20Roll1+x3Roll1sinα1
y2Roll1=y20Roll1+y3Roll1cosα1
z2Roll1=z20Roll1+z3Roll1 (25)
其中,(x20Roll1,y20Roll1,z20Roll1)为第一激光位移传感器(2)空间直角坐标系原点在位移传感器***组成空间直角坐标系下的坐标值,x20Roll1由电机编码器读取,y20Roll1、z20Roll1由光栅尺读取;至此,完成拟合出的滚道ROLL1的滚道中心点坐标从第一激光位移传感器(2)空间斜坐标系(x4Roll1,y4Roll1,z4Roll1)到位移传感器***组成空间直角坐标系(x2Roll1,y2Roll1,z2Roll1)的转换;
步骤7-8:对于第二激光位移传感器(4),同样采用步骤7-1至步骤7-7.得到待测滑块(3)另一滚道ROLL2在位移传感器***组成空间直角坐标系下的滚道ROLL2的滚道中心点坐标(x2Roll2,y2Roll2,z2Roll2)。
8.根据权利要求7所述的滚动直线导轨副滑块型面中径测量算法,其特征在于,所述步骤8具体包括:
步骤8-1:对于滚道ROLL1,待测滑块(3)固定在滑块托架(9)上,移动滑块托架(9),当第一激光位移传感器(2)刚好能够采集到待测滑块(3)内滚道数据时记为位置1,第一激光位移传感器(2)扫掠滑块内滚道,将扫掠到的待测滑块(3)滚道ROLL1数据点拟合出滚道中心点坐标,并将滚道中心点坐标从第一激光位移传感器(2)空间斜坐标系o4-x4y4z4变换到位移传感器***组成的空间直角坐标系o2-x2y2z2,得到所述位置1下,滚道ROLL1截面的拟合滚道中心点在位移传感器***组成空间直角坐标系o2-x2y2z2下的坐标(1x2Roll1,1y2Roll1,1z2Roll1);
步骤8-2:移动滑块托架(9)使之处于不同位置K,相邻位置之间的距离为1CM,根据步骤8-1,滚道ROLL1在不同位置下,截面的拟合滚道中心点在位移传感器***组成空间直角坐标系o2-x2y2z2下的坐标(kx2Roll1,ky2Roll1,kz2Roll1),k=1,2...10;
步骤8-3:根据步骤8-2所得数据点,拟合出滚道ROLL1中心轴线的直线方程,由于滚道ROLL1中心轴线穿过x2o2z2面,则空间直线方程能够简化为:
Figure FDA0003252026630000101
其中(x01,y01,z01)为滚道ROLL1中心轴线上的任意一点,
Figure FDA0003252026630000102
为滚道ROLL1中心轴线的法向向量,滚道ROLL1径向截面为椭圆弧形状,滚道ROLL1中心轴线表示由滚道ROLL1不同径向截面的椭圆弧中心组成的直线,y01=0,x01,z01,mRoll1,nRoll1为要求参数;
直线方程简化成:
Figure FDA0003252026630000103
化为矩阵形式:
Figure FDA0003252026630000104
第K个点方程为:
Figure FDA0003252026630000105
并联10个点的方程为:
Figure FDA0003252026630000106
最小二乘拟合:
Figure FDA0003252026630000107
化简成:
Figure FDA0003252026630000111
求解矩阵即可得到mRoll1,nRoll1,x01,y01,为直线方程的各个参数,得到滚道ROLL1中心轴线的方向向量为
Figure FDA0003252026630000112
步骤8-4:对于滚道ROLL2中心轴线,同样运用步骤8-1至步骤8-3,得到滚道ROLL2中心轴线的方向向量为
Figure FDA0003252026630000113
9.根据权利要求8所述的滚动直线导轨副滑块型面中径测量算法,其特征在于,步骤9:根据步骤8拟合出的滚道ROLL1,ROLL2中心轴线的方向向量,求解滚道ROLL1与滚道ROLL2中心轴线的距离,即为滑块滚道中径,两中心轴线距离为:
Figure FDA0003252026630000114
其中
Figure FDA0003252026630000115
分别为滚道ROLL1与滚道ROLL2中心轴线的方向向量,
Figure FDA0003252026630000116
Figure FDA0003252026630000117
是两条直线上任意一点的连线方向,
Figure FDA0003252026630000118
CN202010157674.8A 2020-03-09 2020-03-09 滚动直线导轨副滑块型面中径的测量算法 Active CN111272088B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010157674.8A CN111272088B (zh) 2020-03-09 2020-03-09 滚动直线导轨副滑块型面中径的测量算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010157674.8A CN111272088B (zh) 2020-03-09 2020-03-09 滚动直线导轨副滑块型面中径的测量算法

Publications (2)

Publication Number Publication Date
CN111272088A CN111272088A (zh) 2020-06-12
CN111272088B true CN111272088B (zh) 2021-10-26

Family

ID=70999371

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010157674.8A Active CN111272088B (zh) 2020-03-09 2020-03-09 滚动直线导轨副滑块型面中径的测量算法

Country Status (1)

Country Link
CN (1) CN111272088B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112504146B (zh) * 2020-07-20 2022-10-28 南京理工大学 一种滚动直线导轨副导轨滚道中径的检测方法
CN117146720B (zh) * 2023-11-01 2024-02-09 中机生产力促进中心有限公司 一种滚动直线导轨副导轨型面检测方法及检测平台

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09323319A (ja) * 1996-06-07 1997-12-16 Tokyo Seimitsu Co Ltd ワイヤソー用溝付ローラの溝形状測定方法及び装置
CN202485598U (zh) * 2012-03-23 2012-10-10 上海飞华汽车部件有限公司 一种测量齿轮中径跳动的自动找准检具
CN104197881B (zh) * 2014-08-28 2017-01-11 中北大学 圆柱外螺纹中径直接测量仪
CN105136097B (zh) * 2015-08-19 2017-08-29 山东建筑大学 非等适应比滚道的滚珠丝杠副弹性变形接触角的确定方法
CN106441168B (zh) * 2016-08-30 2019-01-04 南京理工大学 滚动直线导轨副滑块型面精度的测量方法
CN108871229B (zh) * 2018-06-11 2020-04-17 南京理工大学 一种滚珠螺母螺旋内滚道曲面及外径的测量方法

Also Published As

Publication number Publication date
CN111272088A (zh) 2020-06-12

Similar Documents

Publication Publication Date Title
CN111993159B (zh) 一种轴类工件在位非接触检测方法
CN109341546B (zh) 一种点激光位移传感器在任意安装位姿下的光束标定方法
CN108871229B (zh) 一种滚珠螺母螺旋内滚道曲面及外径的测量方法
CN106441168B (zh) 滚动直线导轨副滑块型面精度的测量方法
CN111272088B (zh) 滚动直线导轨副滑块型面中径的测量算法
CN106767558A (zh) 一种导轨基面直线度误差的解耦辨识方法
Zhang A study on the Abbe principle and Abbe error
JP6431157B1 (ja) 鉄道における建築限界の測定装置及び測定方法
CN109341471A (zh) 基于球列实现三轴机床几何误差检测的鉴定方法
CN111854587B (zh) 一种导轨五自由度运动误差在线测量装置及方法
CN100549614C (zh) 用于检测可在坐标轴上移动的滑架的空间位置的装置
CN114623784A (zh) 一种用于大长度导轨直线度测量的方法和***
CN106052568A (zh) 丝杠型面检测装置位移传感器***安装误差的补偿方法
CN109520417A (zh) 机床几何误差及旋转台转角定位误差检定装置和方法
CN112504146B (zh) 一种滚动直线导轨副导轨滚道中径的检测方法
CN111060056B (zh) 一种平行轮廓精确重构的重构装置及重构方法
CN111708321B (zh) 数控机床刀轴方向动态误差检测装置及方法
CN113091653A (zh) 基于五棱镜测量直线导轨角自由度误差的装置及方法
CN108801187B (zh) 基于坐标变换的导轨滑台运动的几何误差辨识方法
CN216846033U (zh) 基于深矢高工件的内壁测量***
Zha et al. A strategy to evaluate and minimize parallelism errors of a rotor system in a precision rotary table
Xin et al. Measurement techniques for complex surface based on a non-contact measuring machine
CN106989661A (zh) 一种测试机床静压导轨表面形状误差的方法
Ouyang et al. Ball array calibration on a coordinate measuring machine using a gage block
CN110455188B (zh) 单轴平移台与结构光3d传感器联合测量标定方法

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