CN108620955B - Measurement and Identification Method of Rotary Axis Error of Machine Tool Based on Monocular Vision - Google Patents

Measurement and Identification Method of Rotary Axis Error of Machine Tool Based on Monocular Vision Download PDF

Info

Publication number
CN108620955B
CN108620955B CN201810372989.7A CN201810372989A CN108620955B CN 108620955 B CN108620955 B CN 108620955B CN 201810372989 A CN201810372989 A CN 201810372989A CN 108620955 B CN108620955 B CN 108620955B
Authority
CN
China
Prior art keywords
error
lathe
point
machine tool
axis
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
CN201810372989.7A
Other languages
Chinese (zh)
Other versions
CN108620955A (en
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.)
Dalian University of Technology
Original Assignee
Dalian University of 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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201810372989.7A priority Critical patent/CN108620955B/en
Publication of CN108620955A publication Critical patent/CN108620955A/en
Application granted granted Critical
Publication of CN108620955B publication Critical patent/CN108620955B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B23MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23QDETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
    • B23Q17/00Arrangements for observing, indicating or measuring on machine tools
    • B23Q17/24Arrangements for observing, indicating or measuring on machine tools using optics or electromagnetic waves
    • B23Q17/2452Arrangements for observing, indicating or measuring on machine tools using optics or electromagnetic waves for measuring features or for detecting a condition of machine parts, tools or workpieces
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B23MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23QDETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
    • B23Q17/00Arrangements for observing, indicating or measuring on machine tools
    • B23Q17/24Arrangements for observing, indicating or measuring on machine tools using optics or electromagnetic waves
    • B23Q17/2409Arrangements for indirect observation of the working space using image recording means, e.g. a camera
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30108Industrial image inspection
    • G06T2207/30164Workpiece; Machine component

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Mechanical Engineering (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Machine Tool Sensing Apparatuses (AREA)

Abstract

本发明基于单目视觉的机床回转轴误差测量与辨识方法属于机床几何误差测量领域,涉及一种利用单目视觉测量机床回转轴几何误差测量与辨识的方法。该方法采用单目相机在回转轴间歇、连续两次回转的过程中,采集带有编码标记点的定位靶标图像,经过对间歇旋转的图像处理,重建编码标记点在机床坐标系下的坐标,利用理论点与实际点坐标在误差矩阵模型中的对应关系,求解机床独立位置误差,再经过对连续旋转拍摄得到的的图像处理求得机床综合误差,最后在机床综合误差中消除机床独立位置误差得到机床依赖位置误差。该方法有效解决了测量过程繁琐以及辨识模型复杂的问题,采用单目视觉实现非接触测量,设备成本低及设备鲁棒性高。

The invention relates to a monocular vision-based method for measuring and identifying errors of machine tool rotary axes, belonging to the field of machine tool geometric error measurement, and relates to a method for measuring and identifying geometric errors of machine tool rotary axes using monocular vision. This method uses a monocular camera to collect positioning target images with coded marking points during the intermittent and continuous two rotations of the rotary shaft, and reconstructs the coordinates of the coded marking points in the machine tool coordinate system after image processing for intermittent rotation. Using the corresponding relationship between the theoretical point and the actual point coordinates in the error matrix model, the independent position error of the machine tool is solved, and then the comprehensive error of the machine tool is obtained through image processing obtained by continuous rotation shooting, and finally the independent position error of the machine tool is eliminated in the comprehensive error of the machine tool Get the machine tool dependent position error. This method effectively solves the problems of cumbersome measurement process and complex identification model, and uses monocular vision to realize non-contact measurement, with low equipment cost and high equipment robustness.

Description

基于单目视觉的机床回转轴误差测量与辨识方法Measurement and Identification Method of Rotary Axis Error of Machine Tool Based on Monocular Vision

技术领域technical field

本发明属于机床几何误差测量领域,涉及一种利用单目视觉测量机床回转轴几何误差测量与辨识的方法。The invention belongs to the field of geometric error measurement of machine tools, and relates to a method for measuring and identifying the geometric error of a rotary axis of a machine tool by using monocular vision.

背景技术Background technique

随着航空、航天、核电、冶金等高技术领域的高速发展和加工水平的快速提高,对高端数控机床及装备的需求迅速增多,也对数控机床的加工精度提出更高要求。五轴联动数控加工作为空间连续、平滑、和复杂曲面的主要加工手段,较传统的三轴数控加工在加工精度、质量及效率上均有较大提升,而回转轴作为五轴数控机床的重要组成部分,其几何误差在加工中通过多轴耦合形成复杂的空间误差,严重影响机床的加工精度,所以研究五轴联动数控机床中回转轴误差的测量与辨识具有现实意义。目前主要的机床误差检测方法有:标准件测量、激光干涉仪、球杆仪、平面光栅、R-test等,其中视觉测量方法成本低、结构简便、测量精度高等优势,可实现机床误差的非接触测量,因此,发明一种操作方便、测量成本低,测量精度高的回转轴误差检测与辨识技术尤为重要。With the rapid development of aviation, aerospace, nuclear power, metallurgy and other high-tech fields and the rapid improvement of processing levels, the demand for high-end CNC machine tools and equipment has increased rapidly, and higher requirements have been placed on the processing accuracy of CNC machine tools. Five-axis linkage CNC machining is the main processing method for continuous, smooth, and complex surfaces. Compared with traditional three-axis CNC machining, the machining accuracy, quality and efficiency have been greatly improved, and the rotary axis is an important part of five-axis CNC machine tools. Its geometric error forms a complex spatial error through multi-axis coupling during machining, which seriously affects the machining accuracy of the machine tool. Therefore, it is of practical significance to study the measurement and identification of the rotary axis error in the five-axis linkage CNC machine tool. At present, the main machine tool error detection methods include: standard part measurement, laser interferometer, ballbar, plane grating, R-test, etc. Among them, the visual measurement method has the advantages of low cost, simple structure, and high measurement accuracy, which can realize non-destructive detection of machine tool errors. Contact measurement, therefore, it is particularly important to invent a rotary axis error detection and identification technology that is easy to operate, low in measurement cost, and high in measurement accuracy.

华中科技大学的周向东、唐小琦等人发明的专利号为CN 105371793 A的“一种五轴机床旋转轴几何误差一次装卡测量方法”基于机构的空间误差模型,可通过一次装卡测量,计算出五轴数控机床两个旋转轴的十二项几何误差参数,但该方法属于接触式测量方法,且要求多测头逐一触碰正方体基准试件的不同测量面,操作过程较为繁琐。西安瑞特快速制造工程研究有限公司的林京、王琇峰等人发明的专利号为CN 102744648 A的“一种数控机床回转工作台误差测量与分离的方法”通过电涡流传感器使用测量在机床工作台面中心的金属圆环,结合误差解耦方法可完成机床几何误差的非接触测量,该方法结构简单、操作方便,但较适用于直线轴的误差检测,且成本较高,并不适用于回转轴误差的测量与辨识。The patent No. CN 105371793 A invented by Zhou Xiangdong, Tang Xiaoqi and others from Huazhong University of Science and Technology is "a method for measuring the geometric error of the rotation axis of a five-axis machine tool by one-time clamping". Based on the spatial error model of the mechanism, it can be calculated by one-time clamping measurement. Twelve geometric error parameters of the two rotation axes of the five-axis CNC machine tool, but this method is a contact measurement method, and requires multiple probes to touch different measurement surfaces of the cube reference specimen one by one, and the operation process is relatively cumbersome. The patent No. CN 102744648 A invented by Lin Jing and Wang Xiufeng of Xi'an Ruite Rapid Manufacturing Engineering Research Co., Ltd. is "A Method for Measuring and Separating Errors of Rotary Tables of CNC Machine Tools". The metal ring in the center, combined with the error decoupling method, can complete the non-contact measurement of the geometric error of the machine tool. This method is simple in structure and easy to operate, but it is more suitable for the error detection of the linear axis, and the cost is high, and it is not suitable for the rotary axis. Error measurement and identification.

发明内容Contents of the invention

本发明要解决的技术难题是克服现有技术的问题,发明一种基于单目视觉的机床回转轴误差测量与辨识方法。该方法采用单目视觉系统结合定位靶标检测五轴数控机床的回转轴误差,实现了回转轴的单次测量多项误差分离。首先在机床工作台上安装含有编码标记点的定位靶标并对回转轴进行分度,在机床每隔一定分度值间歇旋转一周的过程中采集标记点图像,利用摄像机成像模型重建编码标记点机床坐标系下的坐标,将测量点坐标带入到误差矩阵模型中求解机床独立位置误差;之后设定机床连续回转一周,控制相机每隔相同分度值采集标记点图像,同理可求得机床综合误差,消除机床独立位置误差即为机床依赖位置误差,该方法利用编码标记点结合单目视觉实现数控机床回转轴误差的检测与辨识,有效解决了测量过程繁琐以及辨识模型复杂问题。辨识模型简单且为非接触测量,操作简易,设备成本低、鲁棒性高。The technical problem to be solved by the present invention is to overcome the problems of the prior art and to invent a method for measuring and identifying the error of the rotary axis of the machine tool based on monocular vision. This method uses a monocular vision system combined with a positioning target to detect the error of the rotary axis of a five-axis CNC machine tool, and realizes the separation of multiple errors in a single measurement of the rotary axis. First, install a positioning target containing coded marking points on the machine tool workbench and index the rotary axis, collect the images of the marked points during the intermittent rotation of the machine tool every certain index value, and use the camera imaging model to reconstruct the coded marking point machine tool Coordinates in the coordinate system, bring the coordinates of the measurement points into the error matrix model to solve the independent position error of the machine tool; then set the machine tool to rotate continuously for one circle, control the camera to collect the image of the marker point at the same division value, similarly, the machine tool can be obtained Comprehensive errors, elimination of independent position errors of machine tools is machine-dependent position errors. This method uses coded marking points combined with monocular vision to realize the detection and identification of CNC machine tool rotary axis errors, which effectively solves the problems of cumbersome measurement process and complex identification models. The identification model is simple and non-contact measurement, easy to operate, low equipment cost and high robustness.

本发明所采用的技术方案是一种基于单目视觉的机床回转轴误差测量与辨识方法,其特征为,该方法采用单目相机在回转轴间歇、连续两次回转的过程中,采集带有编码标记点的定位靶标图像,经过对间歇旋转的图像处理,重建编码标记点在机床坐标系下的坐标,利用理论点与实际点坐标在误差矩阵模型中的对应关系,求解机床独立位置误差,再经过对连续旋转拍摄得到的的图像处理可求得机床综合误差,最后在机床综合误差中消除机床独立位置误差得到机床依赖位置误差;方法的具体步骤如下:The technical solution adopted in the present invention is a method for measuring and identifying errors of machine tool rotary axes based on monocular vision, which is characterized in that the method uses a monocular camera to collect information with The positioning target image of the coded mark point, after image processing of intermittent rotation, reconstructs the coordinates of the coded mark point in the machine tool coordinate system, and uses the correspondence between the theoretical point and the actual point coordinates in the error matrix model to solve the independent position error of the machine tool, Then, the comprehensive error of the machine tool can be obtained through image processing obtained by continuous rotation shooting, and finally the independent position error of the machine tool is eliminated in the comprehensive error of the machine tool to obtain the dependent position error of the machine tool; the specific steps of the method are as follows:

第一步,安装实验装置The first step is to install the experimental device

实验装置由单目相机1、定位靶标2、机床工作台3、系统安装架4,夹具5组成,单目相机1垂直机床工作台3安装在系统安装架4上,先将棋盘格标定板装夹于机床工作台3上准备标定单目相机;夹具5安装在机床工作台3上,以方便定位靶标2与棋盘格标定板的更换装夹,将机床待测的回转轴C轴归零;The experimental device consists of a monocular camera 1, a positioning target 2, a machine tool table 3, a system mounting frame 4, and a fixture 5. The monocular camera 1 is installed on the system mounting frame 4 vertically to the machine tool table 3. Clamp it on the machine tool workbench 3 to prepare for calibrating the monocular camera; Fixture 5 is installed on the machine tool workbench 3 to facilitate the replacement and clamping of the positioning target 2 and the checkerboard calibration plate, and reset the C-axis of the machine tool’s rotary axis to be tested to zero ;

定位靶标2为光刻玻璃板,上表面光刻有呈菱形分布的第一、第二、第三、第四圆形编码点69、71、73、75,及位于图像中心的中心编码点119,并以中心编码点119为原点建立机床坐标系;The positioning target 2 is a photolithographic glass plate, and the upper surface is photolithographically engraved with the first, second, third and fourth circular coding points 69, 71, 73, 75 distributed in a diamond shape, and the central coding point 119 located in the center of the image , and establish the machine tool coordinate system with the center code point 119 as the origin;

第二步,相机标定The second step, camera calibration

本发明依据张氏标定法结合高精度棋盘格标定板来标定单目相机的内参数及畸变参数,选取空间一点坐标为(Xp,Yp,Zp),其在像平面上投影点坐标为(x,y),综合考虑实际成像的畸变现象,使用摄像机非线性透视投影模型表达式如下:The present invention calibrates the internal parameters and distortion parameters of the monocular camera based on the Zhang’s calibration method combined with a high-precision checkerboard calibration board, and selects a point coordinate in space as (X p , Y p , Z p ), which projects the point coordinates on the image plane is (x, y), considering the distortion phenomenon of the actual imaging comprehensively, the expression of the nonlinear perspective projection model of the camera is as follows:

其中,f为单目相机的焦距,αx=f/dx与αy=f/dy分别定义为x、y两轴上的归一化焦距,(x0,y0)为像素坐标系原点坐标,由αx、αy、u0、v0四个参数构成的M0为相机内参数矩阵,R为3×3的单位旋转正交阵,t为平移向量,组成的M1为相机外参数矩阵;即由内参数αx、αy、u0、v0及外参数R、t确定机床坐标系与像素坐标系的关系;点(x,y)为投影点的理想位置,(x′,y′)为考虑畸变的实际坐标,δx、δy为横轴与纵轴的非线性畸变值,r为像素坐标系下投影点与原点的距离,k1、k2与k3分别为一阶、二阶与三阶径向畸变系数,p1、p2是一阶、二阶离心畸变系数;通过标定单目相机可得内参数αx、αy、u0、v0及畸变参数k1、k2、k3、p1、p2,进而确定相机成像模型;单目相机1标定好后,将棋盘格标定板更换成定位靶标2,并固定在机床工作台3上;Among them, f is the focal length of the monocular camera, α x =f/dx and α y =f/dy are defined as the normalized focal lengths on the x and y axes respectively, and (x 0 , y 0 ) is the origin of the pixel coordinate system Coordinates, M 0 composed of four parameters α x , α y , u 0 , v 0 is the camera internal parameter matrix, R is a 3×3 unit rotation orthogonal matrix, t is a translation vector, and the M 1 formed is the camera External parameter matrix; that is, the relationship between the machine tool coordinate system and the pixel coordinate system is determined by the internal parameters α x , α y , u 0 , v 0 and external parameters R, t; the point (x, y) is the ideal position of the projected point, ( x′, y′) are the actual coordinates considering the distortion, δ x and δ y are the nonlinear distortion values of the horizontal and vertical axes, r is the distance between the projected point and the origin in the pixel coordinate system, k 1 , k 2 and k 3 are the first-order, second-order and third-order radial distortion coefficients respectively, p 1 and p 2 are the first-order and second-order centrifugal distortion coefficients; the internal parameters α x , α y , u 0 , v can be obtained by calibrating the monocular camera 0 and distortion parameters k 1 , k 2 , k 3 , p 1 , p 2 , and then determine the camera imaging model; after the monocular camera 1 is calibrated, replace the checkerboard calibration plate with a positioning target 2, and fix it on the machine tool table 3 up;

第三步,标记点识别The third step, mark point recognition

对回转轴进行分度,在机床每隔一定分度值间歇旋转一周的过程中采集标记点图像,利用摄像机成像模型重建编码标记点机床坐标系下的坐标;本发明采用圆环编码标记点的定位靶标2来表征机床运动位置,圆环编码中心为黑色圆标记点,外圈为同心圆环区域,用于表征圆环编码的身份信息,称为编码带;该圆环区域按照角度平均分为10段,每段对应一个二进制位;背景色为白色与黑色,相对应的二进制编码为“0”、“1”;从标记点圆心出发,按照一定方向扫描编码带,扫描到黑色码带记为1,白色码带记为0;如果没有扫描到编码带,则从中心开始重新扫描;扫描一周后,整个编码点的码值序列即被全部读出,形成一个二进制序列,最终转化为十进制整数,从而获得编码点的码值;解码后,根据不同编码标记点的码值,再利用标定好的摄像机成像模型内外参数,即可重建各个对应标记点的三维坐标,并经过坐标系基准转换,最终可得机床坐标系下的编码标记点的三维坐标为(Xw,Yw,Zw);Index the rotary axis, collect the image of the mark point during the intermittent rotation of the machine tool every certain index value, and use the camera imaging model to reconstruct the coordinates of the code mark point in the machine tool coordinate system; the present invention uses the circular code mark point Locate the target 2 to represent the movement position of the machine tool. The center of the ring code is a black circle mark point, and the outer ring is a concentric ring area, which is used to represent the identity information of the ring code, which is called a code belt; the ring area is divided according to the angle. There are 10 segments, each segment corresponds to a binary bit; the background color is white and black, and the corresponding binary codes are "0" and "1"; starting from the center of the marked point, scan the coding belt in a certain direction until the black code belt is scanned It is recorded as 1, and the white code band is marked as 0; if the code band is not scanned, start scanning again from the center; after scanning for one week, the code value sequence of the entire code point is read out, forming a binary sequence, and finally converted into Decimal integers, so as to obtain the code value of the code point; after decoding, according to the code value of different code mark points, and then use the internal and external parameters of the calibrated camera imaging model, the three-dimensional coordinates of each corresponding mark point can be reconstructed, and through the coordinate system benchmark Conversion, and finally the three-dimensional coordinates of the coding mark points in the machine tool coordinate system can be obtained as (X w , Y w , Z w );

第四步,回转轴误差测量与辨识The fourth step, rotary axis error measurement and identification

数控机床回转轴误差主要有两种误差,分别是机床独立位置误差(PIGE)与机床依赖位置误差(PDGE),前者主要为机床的装配误差且与机床命令位置无关,所以为常量,后者主要来源于数控机床零、部件的制造缺陷造成的几何误差,与机床命令位置相关,所以为变量且回转角度的一个函数,两者结合可将机床回转轴的综合误差表示为:The rotary axis error of CNC machine tools mainly has two kinds of errors, which are machine tool independent position error (PIGE) and machine tool dependent position error (PDGE). The former is mainly the assembly error of the machine tool and has nothing to do with the machine tool command position, so it is a constant, and the latter is mainly The geometric error caused by the manufacturing defects of the parts and components of the CNC machine tool is related to the command position of the machine tool, so it is a function of the variable and the rotation angle. The combination of the two can express the comprehensive error of the machine tool rotary axis as:

e=ePIGE+ePDGE(θ) (2)e=e PIGE +e PDGE (θ) (2)

首先辨识机床回转轴机床独立位置误差,当控制机床间歇运动到某一位置时,此时回转轴的综合误差项只含有ePIGE,又连接误差共有4项,包括2项线性位置误差,2项转角误差,则其误差矩阵可表示为First, identify the independent position error of the rotary axis of the machine tool. When the machine tool is controlled to move intermittently to a certain position, the comprehensive error term of the rotary axis only contains e PIGE , and there are 4 connection errors, including 2 linear position errors and 2 corner error, then its error matrix can be expressed as

其中,εxc、εyc为C轴实际中心与理想中心在X方向和Y方向的线性位置误差;δxc、δyc为C轴实际轴线与理想轴线绕X轴与Y轴的转角误差。由于误差矩阵T中只含有4项未知量,即通过对回转轴的分度,在机床坐标系的0°,90°,180°和270°的4个角度下,根据单目相机采集编码标记点的图像,利用像素坐标系下的坐标及摄像机成像模型,分别重建编码标记点机床坐标系下的坐标为利用公式(4)得出的方程即可求得数控机床回转轴的机床独立位置误差。设定机床坐标系下4个角度的理论点坐标为求解公式为:Among them, ε xc and ε yc are the linear position errors between the actual center and the ideal center of the C-axis in the X and Y directions; δ xc and δ yc are the rotation angle errors between the actual axis and the ideal axis of the C-axis around the X-axis and the Y-axis. Since the error matrix T contains only 4 unknown quantities, that is, through the graduation of the rotary axis, at the 4 angles of 0°, 90°, 180° and 270° of the machine tool coordinate system, the encoding mark is collected according to the monocular camera The image of the point, using the coordinates in the pixel coordinate system and the camera imaging model, respectively reconstructs the coordinates of the coding mark point in the machine tool coordinate system as The machine tool independent position error of the rotary axis of the CNC machine tool can be obtained by using the equation obtained from formula (4). Set the theoretical point coordinates of the four angles in the machine tool coordinate system as The solution formula is:

之后再辨识机床回转轴机床依赖位置误差,控制机床以一定速度运行,此时回转轴的综合误差项为ePIGE+ePDGE(θ),则其误差矩阵可表示为:Then identify the position error of the rotary axis of the machine tool, and control the machine tool to run at a certain speed. At this time, the comprehensive error term of the rotary axis is e PIGE + e PDGE (θ), and its error matrix can be expressed as:

其中,Δεxc、Δεyc为C轴实际中心与理想中心在X方向和Y方向的线性位置误差;Δδxc、Δδyc为C轴实际轴线与理想轴线绕X轴与Y轴的转角误差,上述4项误差即为所求机床回转轴机床依赖位置误差;之前在静态条件下已经辨识机床回转轴机床依赖位置误差,即εxc、εyc、δxc、δyc为已知量,则T中仍然含有4项未知量,同样在机床坐标系中0°,90°,180°和270°的4个角度下触发相机拍摄并分别重建编码标记点坐标,记为求解公式为:Among them, Δε xc and Δε yc are the linear position errors between the actual center and the ideal center of the C-axis in the X and Y directions; Δδ xc and Δδ yc are the rotation angle errors between the actual axis of the C-axis and the ideal axis around the X-axis and the Y-axis. The four errors are the machine tool-dependent position errors of the rotary axis of the machine tool; the machine-tool rotary axis machine-dependent position errors have been identified under static conditions before, that is, ε xc , ε yc , δ xc , and δ yc are known quantities, then T It still contains 4 unknown quantities, and also triggers the camera to shoot at the 4 angles of 0°, 90°, 180° and 270° in the machine tool coordinate system and reconstructs the coordinates of the coded marker points respectively, which are recorded as The solution formula is:

即在准静态工作下求得数控机床回转轴的综合误差,减去已知的机床依赖位置误差,即为回转轴的机床依赖位置误差。That is, the comprehensive error of the rotary axis of the CNC machine tool is obtained under quasi-static work, and the known machine tool-dependent position error is subtracted, which is the machine tool-dependent position error of the rotary axis.

本发明的有益之处在于利用编码标记点结合单目视觉实现数控机床回转轴误差的检测与辨识,有效解决了测量过程繁琐以及辨识模型复杂问题,采用单目视觉实现非接触测量,具有结构简单,操作简便,设备成本低及设备鲁棒性高。该方法辨识模型简单且为非接触测量,利用编码标记点结合单目视觉实现数控机床回转轴误差的检测与辨识,设备成本低、鲁棒性高。The advantage of the present invention is that the detection and identification of the rotary axis error of the CNC machine tool is realized by using the coding mark points combined with monocular vision, which effectively solves the problems of cumbersome measurement process and complex identification model, and realizes non-contact measurement by monocular vision, which has a simple structure , easy to operate, low equipment cost and high equipment robustness. The identification model of this method is simple and non-contact measurement. The detection and identification of the rotary axis error of the CNC machine tool is realized by using the coded marking points combined with the monocular vision. The equipment cost is low and the robustness is high.

附图说明Description of drawings

图1为基于单目视觉的机床回转轴误差测量与辨识模型图。其中,1-单目相机,2-定位靶标,3-机床工作台,4-系统安装架,5-夹具。Figure 1 is a model diagram of the error measurement and identification of the rotary axis of the machine tool based on monocular vision. Among them, 1-monocular camera, 2-positioning target, 3-machine tool table, 4-system mounting frame, 5-fixture.

图2为定位靶标上编码标记点的图像。其中,69、71、73、75分别为第一、第二、第三、第四圆形编码点,119为图像中心的中心编码点,X、Y为机床坐标系的坐标轴。Figure 2 is an image of the coded markers on the positioning target. Among them, 69, 71, 73, and 75 are the first, second, third, and fourth circular coding points respectively, 119 is the central coding point of the image center, and X and Y are the coordinate axes of the machine tool coordinate system.

图3为机床回转轴误差测量与辨识原理图。其中,εxc、εyc分别为C轴实际中心与理想中心在X方向和Y方向的线性位置误差,δxc、δyc分别为C轴实际轴线与理想轴线绕X轴与Y轴的转角误差,C轴为机床理论回转轴,C’轴为机床实际回转轴。Figure 3 is a schematic diagram of the error measurement and identification of the rotary axis of the machine tool. Among them, ε xc , ε yc are the linear position errors between the actual center and the ideal center of the C-axis in the X direction and the Y direction, respectively, and δ xc , δ yc are the rotation angle errors of the actual axis and the ideal axis of the C-axis around the X-axis and the Y-axis respectively , C axis is the theoretical rotary axis of the machine tool, and C' axis is the actual rotary axis of the machine tool.

图4为机床回转轴误差测量与辨识流程图。Figure 4 is a flow chart of the error measurement and identification of the rotary axis of the machine tool.

具体实施方式Detailed ways

以下结合技术方案及附图详细叙述本发明的具体实施方式。The specific embodiments of the present invention will be described in detail below in conjunction with the technical solutions and accompanying drawings.

图1为基于单目视觉的机床回转轴误差测量与辨识模型图。本发明采用定位靶标2为厚度3mm的光刻玻璃板,上表面光刻有四个呈菱形分布的圆形编码点及一个位于中心的编码点,测量时将其装夹在机床工作台3上,采用单目相机1在回转轴间歇、连续两次回转的过程中,采集带有编码标记点的定位靶标2图像,采集带有编码标记点的定位靶标图像,经过对间歇旋转的图像处理,重建编码标记点在机床坐标系下的坐标,利用理论点与实际点坐标在误差矩阵模型中的对应关系,求解机床独立位置误差,再经过对连续旋转采集的图像处理求得机床综合误差,最后在机床综合误差中消除机床独立位置误差可求得机床依赖位置误差。方法的流程如图4所示,具体步骤如下:Figure 1 is a model diagram of the error measurement and identification of the rotary axis of the machine tool based on monocular vision. In the present invention, the positioning target 2 is a photolithographic glass plate with a thickness of 3mm, and the upper surface is photoengraved with four circular code points distributed in a diamond shape and one code point located in the center, which is clamped on the machine tool workbench 3 during measurement. , the monocular camera 1 is used to collect the image of the positioning target 2 with coded marking points in the process of intermittent and continuous two rotations of the rotary shaft, collect the image of the positioning target with coded marking points, and after the image processing of the intermittent rotation, Reconstruct the coordinates of the coding mark points in the machine tool coordinate system, use the correspondence between the theoretical point and the actual point coordinates in the error matrix model to solve the independent position error of the machine tool, and then obtain the comprehensive error of the machine tool through image processing of continuous rotation acquisition, and finally The machine tool dependent position error can be obtained by eliminating the machine tool independent position error in the machine tool comprehensive error. The flow of the method is shown in Figure 4, and the specific steps are as follows:

第一步,实验装置的安装The first step, the installation of the experimental device

实验装置平台由单目相机1、特征靶标2、机床工作台3、系统安装架4和夹具5注组成,通过夹具5将棋盘格标定板安装到机床工作台3上,再将单目相机1垂直机床工作台3安装在系统安装架4上,将机床待测的回转轴C轴归零。The experimental device platform is composed of a monocular camera 1, a characteristic target 2, a machine tool workbench 3, a system mounting frame 4 and a fixture 5. The checkerboard calibration plate is installed on the machine tool workbench 3 through the fixture 5, and then the monocular camera 1 The vertical machine tool table 3 is installed on the system mounting frame 4, and the C-axis of the rotary axis of the machine tool to be tested is reset to zero.

第二步,相机标定The second step, camera calibration

本方法依据张氏标定法结合高精度棋盘格标定板来标定单目相机的内参数及畸变参数,标定后可确定内参数αx、αy、u0、v0,外参数R、t,及畸变参数k1、k2、k3、p1、p2,进而确定相机成像模型;单目相机标定后将棋盘格标定板更换为定位靶标2,并固定在机床工作台3上。This method calibrates the internal parameters and distortion parameters of the monocular camera based on Zhang’s calibration method combined with a high-precision checkerboard calibration board. After calibration, the internal parameters α x , α y , u 0 , v 0 , external parameters R, t, and distortion parameters k 1 , k 2 , k 3 , p 1 , p 2 , and then determine the camera imaging model;

第三步,标记点特征识别The third step, mark point feature recognition

为准确表征机床运动位置,使用带有圆环编码标记点的定位靶标,如图2所示。将圆环编码带区域按照角度平均分为10段,每段对应一个二进制位,图像处理时按照一定方向扫描编码带,扫描到黑色码带记为1,白色码带记为0,扫描一周后,得到整个编码点的二进制序列,转化为十进制整数即为码值,利用不同编码标记点的码值及摄像机成像模型即可重建编码标记点在机床坐标系下的三维坐标为(Xw,Yw,Zw);In order to accurately characterize the movement position of the machine tool, a positioning target with a ring-coded marking point is used, as shown in Figure 2. Divide the area of the circular coding belt into 10 segments on average according to the angle, and each segment corresponds to a binary bit. Scan the coding belt in a certain direction during image processing, and record the black code belt as 1 and the white code belt as 0. After scanning for one week , to obtain the binary sequence of the entire coding point, which is converted into a decimal integer to be the code value. Using the code values of different coding points and the camera imaging model, the three-dimensional coordinates of the coding point in the machine tool coordinate system can be reconstructed as (X w , Y w , Z w );

第四步,回转轴误差测量与辨识The fourth step, rotary axis error measurement and identification

图3为机床回转轴误差测量与辨识原理图,本方法中,依据误差矩阵模型中的未知量数量,设定机床C轴的分度值为90°,即在机床间歇回转一周的过程中,分别位于0°,90°,180°和270°四个角度下采集一次图像,并重建含有PIGE的编码点三维坐标同时在以第一幅图像的中心点为原点建立的机床坐标系中,根据高精度定位靶标上编码点的位置关系求取理论点三维坐标利用公式(3)、(4),依据理论点与实际点坐标在误差矩阵模型中的对应关系,求解机床独立位置误差(PIGE);Figure 3 is a schematic diagram of the error measurement and identification of the rotary axis of the machine tool. In this method, according to the number of unknown quantities in the error matrix model, the indexing value of the C-axis of the machine tool is set to 90°, that is, during the intermittent rotation of the machine tool, Collect an image at four angles of 0°, 90°, 180° and 270°, and reconstruct the 3D coordinates of the coded points containing PIGE At the same time, in the machine tool coordinate system established with the center point of the first image as the origin, the three-dimensional coordinates of the theoretical point are obtained according to the positional relationship of the coded points on the high-precision positioning target Using formulas (3) and (4), according to the corresponding relationship between the theoretical point and the actual point coordinates in the error matrix model, the machine tool independent position error (PIGE) is solved;

将机床C轴归零后,再控制机床连续回转一周,同样每隔分度值90°触发相机采集图像,经过后续图像处理还原编码点三维坐标经过公式(5)、(6),可求解含有PIGE与PDGE的机床综合误差,此时PIGE为上一步求解的已知量,最后在机床综合误差中消除机床独立位置误差(PIGE)可得到机床依赖位置误差(PDGE)。After returning the C-axis of the machine tool to zero, control the machine tool to rotate continuously for one round, and also trigger the camera to collect images every 90° of the index value, and restore the three-dimensional coordinates of the coding point after subsequent image processing Through the formulas (5) and (6), the comprehensive error of the machine tool including PIGE and PDGE can be solved. At this time, PIGE is the known quantity solved in the previous step. Finally, the independent position error (PIGE) of the machine tool can be obtained by eliminating the comprehensive error of the machine tool (PIGE). Dependent Position Error (PDGE).

本发明根据单目视觉非接触测量的特点,使用简单的误差辨识模型,解决了测量过程繁琐及辨识模型复杂问题,具有结构简单、操作方便的优点,利用编码标记点与单目视觉结合实现数控机床回转轴误差的检测与辨识,具有设备经济性良好、测量稳定性高等特点。According to the characteristics of non-contact measurement of monocular vision, the present invention uses a simple error identification model to solve the problems of cumbersome measurement process and complex identification model, and has the advantages of simple structure and convenient operation. The detection and identification of machine tool rotary axis errors has the characteristics of good equipment economy and high measurement stability.

Claims (1)

1. a kind of measurement of machine tool rotary axis error and discrimination method based on monocular vision, it is characterized in that, this method uses monocular Camera rotating shaft interval, twice in succession turn round during, acquisition have coded markings point positioning target image, by pair The image procossing of intermittent rotary rebuilds coordinate of the coded markings point under lathe coordinate system, utilizes mathematical point and actual point coordinate Corresponding relationship in error matrix model solves lathe independent position error, using what is shot to continuous rotation Image procossing can acquire lathe composition error, and lathe independent position error is finally eliminated in lathe composition error can be obtained lathe Rely on location error;Specific step is as follows for method:
The first step installs experimental provision
Experimental provision is by monocular camera (1), positioning target (2), platen (3), system mounting rack (4) and fixture (5) group At;The vertical platen (3) of monocular camera (1) is mounted on system mounting rack (4), then gridiron pattern scaling board is clamped in Platen is prepared on (3) for camera calibration;Fixture (5) is mounted on platen (3), positioning target is facilitated (2) with the replacement clamping of gridiron pattern scaling board;
Positioning target (2) is photoetching glass plate, and upper surface photoetching has the first circle codification point (69) of the distribution that assumes diamond in shape, the second circle Shape encoded point (71), third circle codification point (73), the 4th circle codification point (75) and centrally located centre code point It (119), is and with centre code point (119) that origin establishes lathe coordinate system;Positioning target (2) is photoetching glass plate, upper surface Photoetching has the first circle codification point (69), the second circle codification point (71), third circle codification point (73), for the distribution that assumes diamond in shape 4 circle codification points (75) and centre code point (119) positioned at picture centre, and be that origin is established with centre code point (119) Lathe coordinate system, four additional circle codification point are used for the motion information of the accurate transmission lathe;Gridiron pattern is demarcated after calibration is good Plate is replaced with positioning target (2), and fixed on platen (3), and lathe rotating shaft C axis to be measured is zeroed;
Second step, camera calibration
The present invention demarcates the intrinsic parameter and distortion ginseng of monocular camera according to Zhang Shi standardization combined high precision gridiron pattern scaling board Number, selection some coordinates in space are (Xp, Yp, Zp), subpoint coordinate is (x, y) in picture plane, comprehensively considers actual imaging Distortion phenomenon uses the non-linear perspective projection model expression of video camera are as follows:
Wherein, f is the focal length of monocular camera, αx=f/dx and αy=f/dy is respectively defined as the normalization focal length on two axis of x, y, (x0, y0) it is pixel coordinate system origin, by αx、αy、u0、v0The M that four parameters are constituted0For camera Intrinsic Matrix, R 3 × 3 unit rotating orthogonal battle array, t are translation vector, the M of composition1For Camera extrinsic matrix number;I.e. by intrinsic parameter αx、αy、u0、v0 And outer parameter R, t determines the relationship of lathe coordinate system and pixel coordinate system;Point (x, y) is the ideal position of subpoint, (x ', y ') For the actual coordinate for considering distortion, δx、δyFor the nonlinear distortion variate of horizontal axis and the longitudinal axis, r is subpoint and original under pixel coordinate system The distance of point, k1、k2With k3Respectively single order, second order and three rank coefficient of radial distortion, p1、p2It is single order, second order centrifugal distortion system Number;Intrinsic parameter α can be obtained by demarcating monocular camerax、αy、u0、v0And distortion parameter k1、k2、k3、p1、p2, and then determine camera at As model;After monocular camera (1) calibration is good, gridiron pattern scaling board is replaced with positioning target (2), and in platen (3) Upper fixation;
Third step, reference point identifying
Rotating shaft is indexed, acquisition label point image during lathe was every certain scale division value intermittent rotary one week, The coordinate under coded markings point lathe coordinate system is rebuild using camera imaging model;The present invention is using annulus coded markings point Positioning target characterizes machine tool motion position, and annulus encoding centre is dark circles mark point, outer ring is concentric annular regions, is used for Characterize the identity information of annulus coding, referred to as coding-belt;The circle ring area is equally divided into 10 sections, every section correspondence one according to angle Binary digit;Background colour is white and black, and corresponding binary coding is " 0 ", " 1 ";From the mark point center of circle, according to Certain orientation scanning encoding band, scanning are denoted as 1 to black code band, and white code band is denoted as 0;If not scanning coding-belt, It is rescaned since center;After run-down, the code value sequence of entire encoded point is all read, and forms a binary system Sequence is eventually converted into decimal integer, to obtain the code value of encoded point;After decoding, according to the code of different coding mark point Value recycles the camera imaging model inside and outside parameter demarcated, can rebuild the three-dimensional coordinate of each correspondence markings point, and pass through Coordinate system Reference Transforming is crossed, the three-dimensional coordinate that can finally obtain the coded markings point under lathe coordinate system is (Xw, Yw, Zw);
4th step, error of rotary axle measurement and identification
Rotating shaft of numerical control machine error is that lathe independent position error and lathe rely on position and miss respectively there are mainly two types of error Difference, the former is mainly the rigging error of lathe and unrelated with lathe command position, so being constant, the latter is mainly derived from numerical control Geometric error caused by the manufacturing defect of lathe parts and components, it is related to lathe command position, so being variable and angle of revolution One function, the two combine and can indicate the composition error of machine tool rotary axis are as follows: e=ePIGE+ePDGE(θ) (2)
Lathe rotating shaft lathe independent position error is recognized first, when controlling lathe intermittent movement to a certain position, is returned at this time The composition error item of shaft contains only ePIGE, and connect error and share 4, including 2 linear position errors, 2 angular errors, Then its error matrix is represented by
Wherein, εxc、εycFor C axis practical center and the linear position error of desired center in the x direction and the y direction;δxc、δycFor C axis The angular errors of practical axis and ideal axis around X-axis and Y-axis;Due to containing only 4 unknown quantitys in error matrix T, by right The indexing of rotating shaft is acquired according to monocular camera and is encoded under 0 ° of lathe coordinate system, 90 °, 180 ° and 270 ° of 4 angles The image of mark point is rebuild coded markings point lathe respectively and is sat using the coordinate and camera imaging model under pixel coordinate system Mark system under coordinate beIt can be acquired using the equation that formula (4) obtains The lathe independent position error of rotating shaft of numerical control machine;Set the mathematical point coordinates of lower 4 angles of lathe coordinate system asSolution formula are as follows:
Lathe rotating shaft lathe is recognized again and relies on location error, and control lathe is run with certain speed, at this time the synthesis of rotating shaft Error term is ePIGE+ePDGE(θ), then its error matrix indicates are as follows:
Wherein, Δ εxc、ΔεycFor C axis practical center and the linear position error of desired center in the x direction and the y direction;Δδxc、 ΔδycIt is the practical axis of C axis and ideal axis around the angular errors of X-axis and Y-axis, above-mentioned 4 errors are that required lathe relies on position Set error;Lathe rotating shaft lathe independent position error, i.e. ε are recognized in a static condition beforexc、εyc、δxc、δycFor The amount of knowing then still contains 4 unknown quantitys, touches under 0 °, 90 °, 180 ° and 270 ° of 4 angles equally in lathe coordinate system in T It sends out camera to shoot and rebuild coded markings point coordinate respectively, be denoted asIt solves Formula are as follows:
The composition error that rotating shaft of numerical control machine is acquired under quasi-static work subtracts known lathe independent position error i.e. Location error is relied on for the lathe of rotating shaft.
CN201810372989.7A 2018-04-18 2018-04-18 Measurement and Identification Method of Rotary Axis Error of Machine Tool Based on Monocular Vision Active CN108620955B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810372989.7A CN108620955B (en) 2018-04-18 2018-04-18 Measurement and Identification Method of Rotary Axis Error of Machine Tool Based on Monocular Vision

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810372989.7A CN108620955B (en) 2018-04-18 2018-04-18 Measurement and Identification Method of Rotary Axis Error of Machine Tool Based on Monocular Vision

Publications (2)

Publication Number Publication Date
CN108620955A CN108620955A (en) 2018-10-09
CN108620955B true CN108620955B (en) 2019-11-19

Family

ID=63694296

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810372989.7A Active CN108620955B (en) 2018-04-18 2018-04-18 Measurement and Identification Method of Rotary Axis Error of Machine Tool Based on Monocular Vision

Country Status (1)

Country Link
CN (1) CN108620955B (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110136204B (en) * 2019-03-19 2021-08-03 浙江大学山东工业技术研究院 Sound-diaphragm dome assembly system based on machine tool position calibration of double-sided telecentric lens camera
CN111604598B (en) * 2020-05-06 2021-04-20 大连理工大学 A tool setting method for a robotic arm feeding laser etching system
CN111678471B (en) * 2020-06-09 2021-07-30 无锡身为度信息技术有限公司 Error identification and compensation method for rotary table of cylindrical coordinate measuring machine
CN112611318B (en) * 2020-12-03 2023-01-06 深圳数马电子技术有限公司 Method and device for measuring motion axis error
CN113280734B (en) * 2021-05-24 2022-04-22 南京大学 Method for correcting errors of manual displacement table by using visual calibration
CN114102256B (en) * 2021-11-04 2022-09-06 清华大学 Method, device and storage medium for recognizing geometric error of rotary axis of machine tool
CN114178907A (en) * 2021-12-10 2022-03-15 上海诺倬力机电科技有限公司 Visual positioning system and method for numerical control machine tool
CN114427833B (en) * 2021-12-17 2024-02-06 哈尔滨理工大学 On-machine detection device for aviation blade based on structured light measurement and point cloud acquisition method
CN119973725B (en) * 2025-04-17 2025-07-25 浙江大学 A dynamic optimization motion variable measuring device for machine tools

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN86102519A (en) * 1986-09-26 1988-04-06 上海市机床研究所 The elimination discrimination technology and the device of reference sphere setting-up eccentricity signal
CN103499416A (en) * 2013-09-28 2014-01-08 大连理工大学 Quick monocular rotational inertia measuring method and quick monocular rotational inertia measuring device
CN105043259A (en) * 2015-08-25 2015-11-11 大连理工大学 Numerical control machine tool rotating shaft error detection method based on binocular vision
JP2016040531A (en) * 2014-08-12 2016-03-24 コニカミノルタ株式会社 Working tool measuring method and measuring device

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN86102519A (en) * 1986-09-26 1988-04-06 上海市机床研究所 The elimination discrimination technology and the device of reference sphere setting-up eccentricity signal
CN103499416A (en) * 2013-09-28 2014-01-08 大连理工大学 Quick monocular rotational inertia measuring method and quick monocular rotational inertia measuring device
JP2016040531A (en) * 2014-08-12 2016-03-24 コニカミノルタ株式会社 Working tool measuring method and measuring device
CN105043259A (en) * 2015-08-25 2015-11-11 大连理工大学 Numerical control machine tool rotating shaft error detection method based on binocular vision

Also Published As

Publication number Publication date
CN108620955A (en) 2018-10-09

Similar Documents

Publication Publication Date Title
CN108620955B (en) Measurement and Identification Method of Rotary Axis Error of Machine Tool Based on Monocular Vision
CN109458958B (en) Calibration method for center position of rotary table in four-axis vision measurement device
CN109520420B (en) Method for determining space coordinates of rotary center of rotary table
CN108340211B (en) Three-dimensional measurement method of CNC machine tool contour error based on monocular vision
CN109596073A (en) A kind of origin position scaling method of the revolving platform central axis of view-based access control model measurement
CN108765494B (en) Multi-camera calibration method based on cylinder calibration object
CN105798704B (en) A kind of lathe face profile error monocular measuring method
CN105043259B (en) Digit Control Machine Tool rotary shaft error detection method based on binocular vision
CN103615980B (en) Method and system for measuring parameters of round holes in plate
JP2013503380A (en) Calibration method for machine tools
CN105382631A (en) Equipment and method for detecting error of rotating shaft of five-axis numerical control machine tool
CN102699733A (en) Method and device for measuring movement locus of automatic tool changing mechanical arm
CN109015110A (en) A kind of machine tool rotary shaft position error modeling and discrimination method
CN112013797B (en) A method and its application for locating the axis of rotation in space based on cylinder and line structure cursors
CN110695982A (en) Mechanical arm hand-eye calibration method and device based on three-dimensional vision
CN113119129A (en) Monocular distance measurement positioning method based on standard ball
CN205192431U (en) Case valve barrel parameter automatic checkout device
CN111338290A (en) Multi-vision-based five-axis numerical control machine tool multifunctional detection method
CN113984279B (en) Batch full-automatic calibration method for pointer type pressure gauges
CN111256612A (en) Machine vision-based method for measuring straight tooth involute small-modulus gear
CN111649670B (en) Method for calibrating center coordinate of rotation axis of swing table
CN113989368B (en) A method and system for high-precision positioning of object surface
CN205588066U (en) Automatic aligning device of machining center
CN216012589U (en) Pointer type pressure gauge batch automatic calibration system
CN112179291B (en) A method for calibrating a self-rotating scanning line structured light three-dimensional measuring device

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