发明名称 一种轧辊轴承座尺寸测量方法
摘要 本发明涉及轧辊轴承座尺寸测量方法,有效解决轧辊轴承座检测,对检测场地、设备和人员要求较高,测量受限制,不能满足轧辊轴承座尺寸实际测量需要的问题,对待检测轴承座进行摄影测量,获取轴承座表面的测量三维点数据,将摄影测量采集到的三维数据,综合轴承座的设计参数对轴承座的平面度,开档距和对称度进行分析计算,将分析结果利用报表形式输出;对分析计算得到的轴承座的平面度、对称度和开档距结果进行数据报表输出和图形化输出,本发明采用相机通过获取轧辊轴承座多幅数字图像,经图像处理、识别及图像匹配、空间交会等获取轴承座表面特征点的三维坐标,采用计算机快速计算出轴承座参数,实现轴承座几何参数快速、高精度检测。
申请公布号 CN102749047A 申请公布日期 2012.10.24
申请号 CN201210207804.X 申请日期 2012.06.21
申请人 郑州辰维科技股份有限公司;上海宝钢工业技术服务有限公司 发明人 王河伟;郭廷钧;李彦辉;张拥军;高钦强
分类号 G01B11/30(2006.01)I;G01B11/02(2006.01)I;G01B11/26(2006.01)I 主分类号 G01B11/30(2006.01)I
代理机构 郑州天阳专利事务所(普通合伙) 41113 代理人 宋金鼎
主权项 1.一种轧辊轴承座尺寸测量方法,其特征在于,由以下步骤实现:(1)、对待检测轴承座用相机进行数字近景摄影测量,获取轴承座表面的测量三维点数据,得到轴承座三维数据,方法是:A、在待检测轴承座两个衬板上布设摄影测量用圆形反光标志,轴承座表面由左右两个衬板平面、前后两个腔孔圆、以及上下面组成,在轴承座左右两个衬板表面按照满足覆盖全部衬板的要求布设25-50个圆形反光标志和8-16个编码标志,标志布设的原则是,摄影时,每张像片至少拍摄到4个编码标志,利用编码标志实现不同像片间的自动拼接;B、在轴承座的前后两个腔孔圆上布设6-9个工装点,用于测量腔孔圆,从而拟合圆计算,求得圆心,为计算轴承座的对称度提供基准;C、在轴承座前面或后面放置概略定向棒和基准尺,概略定向棒的作用是概略计算出像片的外方位元素,基准尺是已经法定计量单位精确标定过的,为整个测量工程提供长度基准;D、用相机对待测轴承座进行摄影,摄影的原则是从轴承座的上下左右各个方向和位置进行拍摄,由于一张像片不可能将所有标志点拍摄下来,故要确保相邻像片间要有4个以上的相同编码标志点;E、对所拍摄的像片用计算机及其软件进行图像处理、概略定向、同名像点自动匹配和光束法平差计算,得到标志点的三维坐标,即得到待测轴承座表面的测量三维点数据,从而得到轴承座三维数据;所述的标志点是指编码标志、圆形标志以及工装点;光束法平差:光束法平差是以每条空间光线为一单元,利用三点共线条件列出误差方程式;所谓三点共线条件,即是指在摄影时,测量点、相应像点和相机中心点是在同一条直线之上的;由共线条件出发,对每个像点可以列出下列两个关系公式:<maths num="0001"><![CDATA[<math><mrow><mfenced open='' close='}'><mtable><mtr><mtd><mi>x</mi><mo>=</mo><mo>-</mo><mi>f</mi><mfrac><mrow><msub><mi>a</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>-</mo><msub><mi>X</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>Y</mi><mo>-</mo><msub><mi>Y</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>Z</mi><mo>-</mo><msub><mi>Z</mi><mi>S</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>a</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>-</mo><msub><mi>X</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>Y</mi><mo>-</mo><msub><mi>Y</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>Z</mi><mo>-</mo><msub><mi>Z</mi><mi>S</mi></msub><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><mi>y</mi><mo>=</mo><mo>-</mo><mi>f</mi><mfrac><mrow><msub><mi>a</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>-</mo><msub><mi>X</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>Y</mi><mo>-</mo><msub><mi>Y</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>Z</mi><mo>-</mo><msub><mi>Z</mi><mi>S</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>a</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>-</mo><msub><mi>X</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>Y</mi><mo>-</mo><msub><mi>Y</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>Z</mi><mo>-</mo><msub><mi>Z</mi><mi>S</mi></msub><mo>)</mo></mrow></mrow></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mn>1</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>其中:x,y是像点坐标;f是摄影焦距;a<sub>1</sub>,a<sub>2</sub>,a<sub>3</sub>,b<sub>1</sub>,b<sub>2</sub>,b<sub>3</sub>,c<sub>1</sub>,c<sub>2</sub>,c<sub>3</sub>是旋转矩阵的9个参数;X<sub>S</sub>,Y<sub>S</sub>,Z<sub>S</sub>是摄站点的物方空间坐标;X,Y,Z是物方点的物方空间坐标;对上式进行线性化后即可得到光束法平差的误差方程式:光束法平差的误差方程式:V=CΔ+L    2);其中:V=[v<sub>x</sub> v<sub>y</sub>]<sup>T</sup>,V是指整体的改正数矩阵,v<sub>x</sub>,v<sub>y</sub>是像点坐标改正值,下标中的<sub>x</sub>,<sub>y</sub>是指像点坐标,<sup>T</sup>是指矩阵的转置;<img file="FDA00001794003200021.GIF" wi="872" he="278" />C是指误差方程系数矩阵,<img file="FDA00001794003200022.GIF" wi="94" he="127" />是像点坐标x对摄站点物方坐标X<sub>S</sub>的偏导值,<img file="FDA00001794003200023.GIF" wi="83" he="120" />是像点坐标x对摄站点物方坐标Y<sub>S</sub>的偏导值,<img file="FDA00001794003200024.GIF" wi="89" he="125" />是像点坐标x对摄站点物方坐标Z<sub>S</sub>的偏导值,<img file="FDA00001794003200025.GIF" wi="78" he="117" />是像点坐标x对旋转角ω的偏导值,<img file="FDA00001794003200026.GIF" wi="78" he="125" />是像点坐标x对旋转角<img file="FDA00001794003200027.GIF" wi="30" he="35" />的偏导值,<img file="FDA00001794003200028.GIF" wi="78" he="116" />是像点坐标x对旋转角κ的偏导值,<img file="FDA00001794003200029.GIF" wi="97" he="126" />是像点坐标y对摄站点物方坐标X<sub>S</sub>的偏导值,<img file="FDA000017940032000210.GIF" wi="85" he="120" />是像点坐标y对摄站点物方坐标Y<sub>S</sub>的偏导值,<img file="FDA000017940032000211.GIF" wi="92" he="128" />是像点坐标y对摄站点物方坐标Z<sub>S</sub>的偏导值,<img file="FDA00001794003200031.GIF" wi="79" he="117" />是像点坐标y对旋转角ω的偏导值,<img file="FDA00001794003200032.GIF" wi="79" he="126" />是像点坐标y对旋转角<img file="FDA00001794003200033.GIF" wi="30" he="35" />的偏导值,<img file="FDA00001794003200034.GIF" wi="79" he="116" />是像点坐标y对旋转角κ的偏导值,Δ=[ΔX<sub>S</sub> ΔY<sub>S</sub> ΔZ<sub>S</sub> Δω<img file="FDA00001794003200035.GIF" wi="64" he="46" />Δκ]<sup>T</sup>,Δ是指误差方程解的矩阵,ΔX<sub>S</sub>,ΔY<sub>S</sub>,ΔZ<sub>S</sub>是指摄站点物方坐标X<sub>S</sub>,Y<sub>S</sub>,Z<sub>S</sub>的改正数,Δω,<img file="FDA00001794003200036.GIF" wi="84" he="46" />Δκ是指旋转角ω,<img file="FDA00001794003200037.GIF" wi="51" he="35" />κ的改正数,<sup>T</sup>是指矩阵的转置;L=[x<sup>0</sup>-x y<sup>0</sup>-y]<sup>T</sup>L是指误差方程的常数项矩阵,x,y是像点坐标,x<sup>0</sup>,y<sup>0</sup>是像点坐标x,y的近似值,<sup>T</sup>是指矩阵的转置;(2)、将摄影测量采集到的轴承座三维数据,综合轴承座的设计参数对轴承座的平面度,开档距和对称度进行分析计算,并将分析结果利用报表形式输出;所述的轴承座的设计参数又称轴承座基本参数,包括轴承座的腔孔设计长度,衬板设计宽度,衬板面粘贴的编码点编号,方法是:A、利用轴承座的腔孔设计长度和衬板设计宽度来确定衬板的范围,从而在MPS/RBMS轧辊轴承座测量系统软件中自动剔除范围的杂点,利用衬板面粘贴的编码点编号在软件中分辨出轴承座的前后左右,方便自动分析出是轴承座哪个面的平面度;B、利用衬板上的测量三维点数据集,使用最小二乘平面拟合,得到一个拟合平面,并计算衬板上的测量三维点到拟合平面的最大距离,得轴承座的开档距;C、利用衬板上的测量三维点,测量到对面衬板上拟合平面的距离,得到当前位置点相对于对面的偏移程度;D、利用腔孔圆上的工装点拟合得到两个拟合腔孔圆,并计算得到圆心,从而以圆心轴线为基准,计算两衬板到此基准的对称度,从而得到衬板上各点相对于中心位置的对称性;(3)、数据报表及图形化输出:对分析计算得到的轴承座的平面度、对称度和开档距结果进行数据报表输出和图形化输出。
地址 450000 河南省郑州市郑州高新技术产业开发区冬青街7号国家科技园A2座12层