阅读说明：本技术 基于向量方向测量的数控机床平动轴几何误差辨识方法 (Method for identifying geometric error of translational axis of numerical control machine tool based on vector direction measurement ) 是由 王金栋 唐雷雨 郭浩然 赖洋 谢宇鸿 闫府 于 2021-08-26 设计创作，主要内容包括：本发明公开了基于向量方向测量原理的数控机床平动轴几何误差辨识方法,在机床运动部件上设置不同测量点,并将这些测量点依次相连构成空间向量,利用在机床运动过程中构成空间向量方向的改变仅与机床轴的角位移误差有关,而与线位移误差无关这一特性,利用激光跟踪仪采用多站分时测量原理对机床平动轴运动过程中测量点构成向量方向变化的测量,从而依次辨识出机床平动轴的角位移误差和线位移误差。建立基于空间向量方向测量原理的数控机床平动轴几何误差测量的数学模型,依次推导出了基于精密数控转台测量的激光跟踪仪基站位置自标定算法、测量点确定算法、机床平动轴单项误差分离算法。该方法具有精度高、效率快、成本较低、算法较简单的优点。(The invention discloses a method for identifying geometric errors of a translational shaft of a numerical control machine tool based on a vector direction measurement principle, wherein different measurement points are arranged on a moving part of the machine tool, the measurement points are sequentially connected to form a space vector, the characteristic that the change of the direction of the space vector formed in the moving process of the machine tool is only related to the angular displacement error of the machine tool shaft but is not related to the linear displacement error is utilized, and the change of the direction of the vector formed by the measurement points in the moving process of the translational shaft of the machine tool is measured by a laser tracker by adopting a multi-station time-sharing measurement principle, so that the angular displacement error and the linear displacement error of the translational shaft of the machine tool are sequentially identified. A mathematical model for measuring the geometric errors of the translational axis of the numerical control machine tool based on the space vector direction measurement principle is established, and a laser tracker base station position self-calibration algorithm, a measuring point determination algorithm and a machine tool translational axis single error separation algorithm based on precise numerical control turntable measurement are deduced in sequence. The method has the advantages of high precision, high efficiency, low cost and simple algorithm.)

1. The method for identifying the geometric error of the translational axis of the numerical control machine tool based on vector direction measurement is characterized by comprising the following steps of:

s1, mounting the precision numerical control turntable near the main shaft of the machine tool, and translating along with the main shaft or the workbench; the cat eye is arranged on the precise numerical control rotary table and can rotate along with the precise numerical control rotary table;

s2, installing the laser tracker near the front of the machine tool in a time-sharing manner, and defining the position of the laser tracker as the position of a base station;

s3, controlling the precise numerical control rotary table to rotate by a certain angle theta every time, recording the distance measurement data of the current laser tracker, constructing a least square model based on the ohm distance of the measurement point by using the distance data of the measurement point obtained by the measurement of the laser tracker, and determining the homogeneous transformation matrix between the instrument coordinate system of the laser tracker and the coordinate system of the precise numerical control rotary table, thereby determining the initial position coordinates of each base station of the laser tracker on the rotary table;

s4, controlling a machine tool translation component to feed along the axis of the machine tool according to a preset path, measuring once when the machine tool moves for a certain distance, sequentially rotating a precision turntable arranged on a main shaft of the machine tool for 120 degrees, wherein measuring points are respectively named as A, B, C, and a laser tracker sequentially measures at four base station positions;

s5, establishing a nonlinear redundancy equation set by using the distance data between the measuring point and each base station measured by the laser tracker, and solving according to the minimum two-component principle to determine the space coordinate of the measuring point;

s6, using the space coordinates of a series of measuring points when the translational component moves along the axis obtained in the step S5, and 3 vectors formed by two A, B, C three pointsThe direction of the vector will be changed in the movement process of the machine tool, the change of the direction of the vector is only related to the angular displacement error in the movement of the machine tool, a series of vector establishment equation sets are formed according to a homogeneous transformation matrix and measuring points caused by the angular displacement error, and the angular displacement errors can be identified by solving the equation sets;

s7, determining a series of new points by taking the space coordinates of a series of measuring points when the translational component moves along the axis obtained in the step S5 and the angular displacement errors of each item obtained in the step S6, taking the midpoint of A, B as T, the midpoint of B, C as R and the midpoint of C, A as S, respectively establishing a motion error equation for T, R, S to form a redundancy equation set, and substituting the redundancy equation set into the solution of each angular displacement error, so that each linear displacement error can be identified, and the complete identification of the geometric errors of the numerical control machine tool is completed.

2. The method of claim 1, wherein in step S2, at least four base station positions are selected, and any three base station positions are not collinear.

3. The method for identifying geometric errors of a translational axis of a numerically controlled machine tool based on vector direction measurement according to claim 1, wherein the step S3 further comprises the steps of:

s31, determining an initial value of the distance from the centre of the cat eye to the rotation center of the rotary table;

s32, determining an initial value of the position of the base station where the laser tracker is located;

and S33, calibrating the base station.

4. The method for identifying geometric errors of a translational axis of a numerically controlled machine tool based on vector direction measurement as claimed in claim 1, wherein in step S3, when the precision numerically controlled turntable rotates by a certain angle θ, the rotating part is controlled to stop moving, the measurement result of the laser tracker at the current position is recorded, then the laser tracker is moved to the next base station position, and the above measurement process is repeated until the measurement is completed at all the base station positions.

5. The method for identifying geometric errors of a translational axis of a numerical control machine tool based on vector direction measurement according to claim 1, wherein in the step S4, at least four different cat eye positions are measured at each base station position to obtain data redundancy, thereby improving reliability and measurement accuracy of the measurement system.

6. The method for identifying geometric errors of translational axes of numerical control machine tools according to claim 1, wherein in step S4, when the translational member moves to each measurement point, the translational member is controlled to stop moving; when the rotating part of the precise numerical control rotary table moves to a measuring position, the rotary table is controlled to stop moving.

7. The method for identifying the geometric errors of the translational axis of the numerical control machine tool based on the vector direction measurement according to claim 1, wherein in the step S1, the precise numerical control rotary table is installed at a proper position near the spindle of the turning and milling center and translates along the axial direction of the numerical control machine tool along with the machine tool workbench, the translational part of the numerical control machine tool is controlled to feed along the axial direction of the machine tool along a preset path, and in a motion area, the laser tracker tracks and measures the motion of the translational part in real time.

8. The method for identifying the geometric errors of the translational shaft of the numerical control machine tool based on the vector direction measurement of the claim 7 is characterized in that when the translational part of the numerical control machine tool is fed along the axis direction of the machine tool, three measurement points are respectively measured every certain distance of movement.

Technical Field

The invention belongs to the technical field of laser precision measurement, and particularly relates to a method for identifying geometric errors of a translational shaft of a numerical control machine tool based on vector direction measurement.

Background

With the continuous development of modern manufacturing technology, the demand of high-precision numerical control machine tools is increasing. The precision is an important performance index of the numerical control machine tool and directly influences the processing quality of a workpiece. How to economically improve the machining precision of the machine tool is a key problem, and scholars at home and abroad carry out extensive and intensive research in this respect. Compared with a machine tool error hardware compensation technology, error measurement and compensation are used as a soft technology for improving the machine tool precision, the method has the advantages of wide applicability, low cost and the like, and is widely applied to the field of machine tool precision compensation. In the process of processing a workpiece by a machine tool, a plurality of error factors exist, such as geometric errors, thermal deformation errors, cutting force errors, clamping errors and the like, and the final processing precision of the workpiece is influenced by the errors. Among the error factors, the geometric error is a main error factor, and meanwhile, the geometric error is relatively stable, has repeatability in a long working time of a machine tool and is easy to measure and compensate. The measurement and compensation of the geometric error of the machine tool are one of the effective ways for improving the machining precision of the machine tool at present.

At present, there are many methods for measuring the geometric error of the translational axis of a machine tool, and the methods mainly comprise the following steps: a physical reference method, a cue-meter measuring method, an orthogonal grating measuring method, a laser interferometry method, and the like. In recent years, a laser tracker is applied to error detection of machine tools and three-dimensional measuring machines as a new type of three-dimensional measuring instrument. Laser trackers typically employ a spherical coordinate measurement method, and for this mode of measurement, angular measurement errors are a major source of coordinate measurement errors. The laser tracker angle error is assumed to be 1 radian per second, resulting in an error of about 5 μm per meter. At present, the requirement on the detection precision of a machine tool is higher and higher, and the existing measuring method based on the spherical coordinates is difficult to meet the detection requirement on the precision of the machine tool. In order to improve the measurement precision and consider the measurement cost, a multi-station time-sharing measurement method is adopted. The method adopts the GPS measurement principle, only relates to distance measurement in the measurement process, and effectively avoids the influence of angle measurement error on the measurement result. However, this method also has certain disadvantages: 1) during measurement, a machine tool is required to feed in a 3D space, and the measurement track is complex. When the motion trail of the machine tool is a 2D plane and a 1D straight line, the machine tool error cannot be identified; 2) in the base station calibration process, the calibration precision of the base station can be influenced by the self error of the machine tool, so that the measurement precision of the method is further improved to a certain extent. How to simplify the measurement process and the error identification algorithm and overcome the defect that the base station calibration precision is influenced by the error of the machine tool, which has positive effect on further improving the overall machining precision of the multi-axis machine tool.

Disclosure of Invention

The invention aims to solve the problems and provides a method for identifying the geometric error of a translational shaft of a numerical control machine tool based on vector direction measurement, which can quickly and accurately identify the geometric error of the numerical control machine tool.

In order to solve the technical problems, the technical scheme of the invention is as follows: the method for identifying the geometric error of the translational axis of the numerical control machine tool based on vector direction measurement comprises the following steps:

s1, mounting the precision numerical control turntable near the main shaft of the machine tool, and translating along with the main shaft or the workbench; the cat eye is arranged on the precise numerical control rotary table and can rotate along with the precise numerical control rotary table;

s2, installing the laser tracker near the front of the machine tool in a time-sharing manner, and defining the position of the laser tracker as the position of a base station;

s3, controlling the precise numerical control rotary table to rotate by a certain angle theta every time, recording the distance measurement data of the current laser tracker, constructing a least square model based on the ohm distance of the measurement point by using the distance data of the measurement point obtained by the measurement of the laser tracker, and determining the homogeneous transformation matrix between the instrument coordinate system of the laser tracker and the coordinate system of the precise numerical control rotary table, thereby determining the initial position coordinates of each base station of the laser tracker on the rotary table;

s4, controlling a machine tool translation component to feed along the axis of the machine tool according to a preset path, measuring once when the machine tool moves for a certain distance, sequentially rotating a precision turntable arranged on a main shaft of the machine tool for 120 degrees, wherein measuring points are respectively named as A, B, C, and a laser tracker sequentially measures at four base station positions;

s5, establishing a nonlinear redundancy equation set by using the distance data between the measuring point and each base station measured by the laser tracker, and solving according to the minimum two-component principle to determine the space coordinate of the measuring point;

s6, using the space coordinates of a series of measuring points when the translational component moves along the axis obtained in the step S5, and 3 vectors formed by two A, B, C three pointsThe direction of the vector will be changed in the movement process of the machine tool, the change of the direction of the vector is only related to the angular displacement error in the movement of the machine tool, a series of vector establishment equation sets are formed according to a homogeneous transformation matrix and measuring points caused by the angular displacement error, and the angular displacement errors can be identified by solving the equation sets;

s7, determining a series of new points by taking the space coordinates of a series of measuring points when the translational component moves along the axis obtained in the step S5 and the angular displacement errors of each item obtained in the step S6, taking the midpoint of A, B as T, the midpoint of B, C as R and the midpoint of C, A as S, respectively establishing a motion error equation for T, R, S to form a redundancy equation set, and substituting the redundancy equation set into the solution of each angular displacement error, so that each linear displacement error can be identified, and the complete identification of the geometric errors of the numerical control machine tool is completed.

Further, in step S2, at least four base station positions are selected, and any three of the base station positions are not collinear.

Further, the step S3 further includes the following steps:

s31, determining an initial value of the distance from the centre of the cat eye to the rotation center of the rotary table;

s32, determining an initial value of the position of the base station where the laser tracker is located;

and S33, calibrating the base station.

Further, in step S3, when the precision numerical control turntable rotates by a certain angle θ, the rotating part is controlled to stop moving, the measurement result of the laser tracker at the current position is recorded, then the laser tracker is moved to the next base station position, and the above measurement process is repeated until the measurement is completed at all base station positions.

Further, in step S4, at least four different cat eye positions are measured at each base station position to obtain data redundancy, thereby improving the reliability and the measurement accuracy of the measurement system.

Further, in step S4, when the translation component moves to each measurement point position, controlling the translation component to stop moving; when the rotating part of the precise numerical control rotary table moves to a measuring position, the rotary table is controlled to stop moving.

Further, in step S1, the precision numerical control rotary table is installed at a suitable position near the spindle of the turn-milling machining center, and translates along the axial direction of the numerical control machine along with the machine tool table, so as to control the translational component of the numerical control machine to feed along the axial direction of the machine tool along a preset path, and in the motion region, the laser tracker tracks and measures the motion of the translational component in real time.

Further, when the numerical control machine tool translation component feeds along the axis direction of the machine tool, three measuring points are measured respectively when the numerical control machine tool translation component moves for a certain distance.

The invention has the beneficial effects that: according to the method for identifying the geometric errors of the translational shaft of the numerical control machine tool based on the vector direction measurement, the laser tracker tracks and measures the points with the same translational component of the multi-shaft machine tool at least four base station positions on the rotary table of the multi-shaft machine tool in sequence, and the space coordinates of a series of measuring points on the precise numerical control rotary table are respectively determined by using the distance data of the measuring points obtained by measurement. Through the obtained space coordinates of a series of measuring points, every two measuring points obtained from the same position of a translation part of the machine tool form space vectors, and the angular displacement error is identified firstly only related to the angular displacement error by utilizing the change of the directions of the vectors; and respectively fitting out a vector in the direction of the plane axis of the machine tool, and introducing the angular displacement error to identify the linear displacement error, thereby realizing complete identification of the geometric error of the numerical control machine tool. The method solves the problem that the geometric error of the numerical control machine tool is measured by using the laser tracker at present. Meanwhile, the method is based on the GPS principle, only the distance quantity is measured in the measuring process, and the influence of the angle measurement error of the laser tracker on the overall measuring precision is effectively avoided, so that the on-site overall measuring precision of the laser tracker is greatly improved. The advantages of rapid measurement of the laser tracker and high precision of the precise numerical control rotary table are utilized to realize high-efficiency high-precision measurement of the geometric error of the numerical control machine. The method has the advantages of rapidness, high precision, simple error identification and the like, and can be used for rapidly and accurately identifying the geometric error of the numerical control machine tool. And a foundation is laid for further improving the integral processing precision of the numerical control machine tool.

Drawings

FIG. 1 is a schematic diagram of the position of a base station, a precision numerically controlled turntable and a cat eye relative to a numerically controlled machine tool in the present invention;

FIG. 2 is a schematic diagram of a method for identifying angular displacement errors of a numerically controlled machine tool according to the present invention;

FIG. 3 is a schematic view of a turret coordinate system of the present invention;

FIG. 4 is a schematic diagram of a variation of the coordinate system of the present invention;

FIG. 5 is a diagram of a base station calibration coordinate system according to the present invention.

Detailed Description

The invention is further described with reference to the following figures and specific embodiments:

as shown in fig. 1 to 5, the method for identifying geometric errors of a translational axis of a numerical control machine based on vector direction measurement provided by the invention comprises the following steps:

s1, mounting the precision numerical control turntable near the main shaft of the machine tool, and translating along with the main shaft or the workbench; the cat eye is installed on the accurate numerical control revolving stage to can follow accurate numerical control revolving stage and rotate together.

In the step, the precise numerical control rotary table is arranged at a proper position near a main shaft of a turn-milling machining center, and translates along the axial direction of the numerical control machine along with a machine tool workbench, so that the translation component of the numerical control machine is controlled to feed along the axial direction of the machine tool along a preset path, and in a motion area, a laser tracker tracks and measures the motion of the translation component in real time.

And S2, installing the laser tracker near the front of the machine tool in a time-sharing mode, namely, at different positions, and defining the position of the laser tracker as the position of the base station.

In this step, at least four base station positions are selected, and any three are not collinear.

S3, controlling the precise numerical control rotary table to rotate by a certain angle theta every time, recording the distance measurement data of the current laser tracker, constructing a least square model based on the ohm distance of the measurement point by using the distance data of the measurement point obtained by the measurement of the laser tracker, and determining the homogeneous transformation matrix between the instrument coordinate system of the laser tracker and the coordinate system of the precise numerical control rotary table, thereby determining the initial position coordinates of each base station of the laser tracker on the rotary table.

In the step, when the precise numerical control rotary table rotates by a certain angle theta, the rotating part is controlled to stop moving, the measurement result of the laser tracker at the current position is recorded, then the laser tracker is moved to the next base station position, and the measurement process is repeated until the measurement is completed at all the base station positions.

In the base station calibration process, a nonlinear redundancy equation set for base station calibration can be established by using a large amount of measured data. When solving the system of nonlinear redundancy equations, it should be linearized first and converted into a linear redundancy equation to solve. The main approach of the present linearization is to perform taylor expansion at a certain point position of the function and omit the terms after the first-order partial derivative, which involves selecting which position to perform taylor expansion, i.e. determining the initial value of the taylor expansion position. Whether the initial value position is selected accurately or not directly influences the calculation precision and efficiency. When the selected initial value is far from the true value, it may cause the iterative computation to be non-convergent, and thus a solution result cannot be obtained, so how to determine the initial value of the base station location is a key problem. And constructing a laser tracker instrument coordinate system based on the measured point ohm distance and a four-axis machine tool self machine tool coordinate system transformation least square model by utilizing a large number of measured point space coordinates obtained by measuring by the laser tracker, thereby determining the initial position of each base station of the laser tracker on the rotary table.

Step S3 specifically includes the following steps:

and S31, determining the initial distance value from the centre of the cat eye to the rotation centre of the rotary table.

During measurement, the coordinates of a series of measurement points on the turntable measured by the laser tracker in the instrument coordinate system are A'_i(x′_i,y′_i,z′_i) As shown in FIG. 5, the series of measurement points are K, G, J, I, M … …, A'_iDenotes these points, x'_i、y′_i、z′_iThe space plane equation of the measuring points in the instrument coordinate system is fitted through a space plane, and the space plane equation where the measuring points are located is assumed to be z ' ═ Ax ' + By ' + C, and A, B, C are respectively constant. During calibration, the rotary table rotates a certain angle for one time to measure, and N measuring points are shared. According to the least squares fitting principle, the objective function is:

according to the extreme principle, if F (a, B, C) is to be minimized, then:

whereinThe partial derivatives are calculated according to the extreme value principle,for the sign of the partial derivative, the coordinates of each point are measured as known quantities, A, B, C as unknown quantities, F is a function of A, B, C,

bringing formula 1 into formula 2, which is obtained from formula 2:

the equation of the spatial plane where the measuring point is located can be determined by the formula (3).

Adjacent 3M points on the rotary table_i(x′_i,y′_i,z′_i)，M_i+1(x′_i+1,y′_i+1,z′_i+1)，M_i+2(x′_i+2，y′_i+2，z′_i+2) Refers to three adjacent measuring points M_i、M_i+1、M_iThe spatial coordinate of +2, for example, the center of a spatial circle determined by K, G, J in fig. 5, that is, the coordinate of the rotation center of the turntable in the measurement coordinate system of the laser tracker instrument is O '(u, v, w), and it is easy to know that the center O' is in the plane determined by equation (3), the following equation set can be established:

the position coordinate O' (u, v, w) of the rotation center of the rotary table under the measurement coordinate system of the laser tracker instrument can be determined by solving the formula (4), and the rotation radius corresponding to a series of measurement points on the rotary table is as follows: r⁰The cat eye is the radius of gyration from the center of rotation of the turntable, and the cat eye is fixed on the turntable.

The initial distance from the center of the cat eye to the rotation center of the rotary table can be determined through the formula (5). The coordinates of the midpoint in step S31 are all in the instrument coordinate system.

And S32, determining the initial value of the position of the base station where the laser tracker is located.

The spatial coordinates of the measurement point given by the laser tracker are the coordinates in the instrument coordinate system and this coordinate system is noted as the old coordinate system X ' Y ' Z '. In order to determine the initial value of the base station position, the measuring point needs to be located at the instrument coordinateConverting the coordinates under the system into a turntable coordinate system, and recording the turntable coordinate system as a new coordinate system O₁X₁Y₁Z₁As shown in fig. 3. The change in the coordinate system is shown in fig. 4.

The instrument coordinate system and the turntable coordinate system can be coincided through the rotation and translation transformation of the coordinate system. Assuming that the instrument coordinate system moves by Δ X along the X axis, moves by Δ Y along the Y axis, moves by Δ Z along the Z axis, rotates by an angle α around the X axis, rotates by an angle β around the Y axis, and rotates by γ around the Z axis, the instrument coordinate system can be overlapped with the turntable coordinate system, and the corresponding homogeneous transformation matrix T at this time is the same as equation 27. Assuming that the coordinate of each measurement point in the instrument coordinate system is A'_i(x′_i,y′_i，z′_i) The theoretical reference coordinate of each measurement point in the turntable coordinate is A as in step S31_i(x_i，y_i,z_i) This is a switch from the instrument coordinate system to the spatial coordinates in the turret coordinate system, where x_i＝Rcosθ_i，y_i＝Rsinθ_i，z_i0. Let U ═ x'_i y′_i z′_i 1]^T，U＝[x_i y_i z_i 1]^TAnd constructing a least square model based on the Eugenia distance of the measuring points:

in equation (6), 6 unknown parameters (Δ x, Δ y, Δ z, α, β, γ) are involved, and these six parameters are the unknown quantities of the spatial left transformation, which refer to the displacement quantities in the x, y, z directions and the rotation quantities around the x, y, z axes, respectively. The unknown parameters can be determined by solving through a nonlinear least square algorithm, so that a homogeneous transformation matrix corresponding to the instrument coordinate system transformed to the turntable coordinate system is determined, and the initial position of the laser tracker in the turntable coordinate system can be easily determined.

And S33, calibrating the base station.

During measurement, the cat eye is arranged on the precise rotary table, and the distance between the center of the cat eye and the rotary center of the rotary table is R.

And controlling the turntable to stop moving when the turntable rotates by a certain angle theta, and recording the distance measurement data of the current laser tracker. The distance between the center of the cat eye and the rotation center of the rotary table is R, and the rotary table can measure every theta angle so as to obtain a series of measuring points such as G, H and I. With a precision turntable rotation center O₁As origin of coordinates, with O₁The straight line on which G is located is Y₁Axial, upward direction being Y₁In the positive direction of the axis to pass through O₁Point perpendicular to O₁Straight line of G is X₁Axis, backward direction being X₁The positive direction of the axis. Determining Z according to the right-handed screw rule₁Axes, thereby establishing a base station calibration coordinate system on the turntable, as shown in fig. 5, where xyz is the machine coordinate system.

The position coordinate of the laser tracker under the coordinate system of the turntable is assumed to be P₁(x_p1,y_p1,z_p1) The space coordinate of the laser tracker under the coordinate system of the rotary table, the subscript 1 of P1 indicates the position of the first base station, and the distance measurement reading of the laser tracker is l when the rotary table rotates by theta angles in turn_1iFor the measurement points, the following system of equations can be established: l_1iSubscript 1 indicates the base station position 1, i corresponds to the different measurement points, and indicates the distance reading of the laser tracker when the ith point was measured at base station position 1

Equation (7) is a nonlinear system of redundant equations that can be solved according to the least squares principle.

The residual is recorded as: f. of_iFor residual, i refers to the ith measurement point, θ_iIs the angle value of the turn corresponding to the ith measuring point

Get x⁰，y⁰，z⁰,R⁰Is x_p1，y_p1,z_p1Approximation of RThe values, namely:

x_p1＝x⁰+Δx，y_p1＝y⁰+Δy，z_p1＝z⁰+Δz，R＝R⁰+ΔR (9)

the formula (8) is represented by (x) according to Taylor series⁰，y⁰，z⁰,R⁰) The expansion is obtained by omitting the terms after the first partial derivative in order to eliminate the nonlinear term

Order to

The formula 8, 9 and 11 are changed into 10 after finishing, so that the product can be obtained:

order:

then:

equation (14) is a system of linear redundancy equations with an objective function of:

by the extreme principle, to minimize H, we must have:

at the same time, the user can select the desired position,

it can be seen that the extreme value obtained by each equation in equation (16) is the minimum value, and satisfies the least square condition, and Δ x, Δ y, and Δ z can be obtained by solving equation (16) as shown in equation (18).

When Δ x, Δ y, Δ z, and Δ R are obtained, the position of the first base station where the laser tracker is located and the distance from the center of the cat eye to the rotation center of the turntable can be determined according to equations (19) and (20).

R＝R⁰+ΔR (20)

Through the process, the first base station position P of the laser tracker under the turntable coordinate system can be calibrated₁Repeating the above process to calibrate other base stations P where the laser tracker is located₂、P₃、P₄The position coordinates of (a).

S4, controlling the machine tool translation component to feed along the axis of the machine tool according to a preset path, measuring once when the machine tool moves for a certain distance, sequentially rotating the precision rotary table arranged on the main shaft of the machine tool for 120 degrees, wherein the measuring points are respectively named as A, B, C, and the laser tracker sequentially measures at the positions of four base stations.

In the step, at least four different cat eye positions are measured at each base station position to obtain data redundancy, so that the reliability and the measurement accuracy of the measurement system are improved. When the translational component moves to each measuring point position, controlling the translational component to stop moving; when the rotating part of the precise numerical control rotary table moves to a measuring position, the rotary table is controlled to stop moving. The certain distance is 50 mm.

The laser tracker is located at a first base station position P₁To move the target scope cat eye to A, B, C and measure A, B, C in turn to base station P₁The precise numerical control rotary table is controlled to stop rotating and the machine tool stops axial movement during measurement, then the supporting plate is controlled to move to the next position along the x axis, and the measurement process is repeated until the position P of the first base station is reached₁The measurement of the machine tool movement is completed. Then the laser tracker is respectively moved to P₂、P₃、P₄At the position and repeated at the base station P₁And (4) measuring at the positions until the laser tracker finishes measuring the movement of the machine tool at all the positions of the base station.

And S5, establishing a nonlinear redundancy equation set by using the distance data between the measuring point and each base station measured by the laser tracker, and solving according to the least squares principle to determine the space coordinate of the measuring point.

The base station positions obtained by the assumed calibration are respectively P₁(x_p1,y_p1,z_p1)、P₂(x_p2,y_p2,z_p2)、P₃(x_p3,y_p3,z_p3)、P₄(x_p4,y_p4,z_p4) Measuring point M during movement of machine tool_iDistances to the base stations are l'_1iSubscript 1 indicates the base station position 1, i corresponds to the different measurement points, and indicates the distance reading, l ', of the laser tracker at the point of measurement of the ith at base station position 1'_2i,l′_3i,l′_4i. According to the GPS measurement principle, M in the measurement process is measured_i(x_mi,y_mi,z_mi) The following equation can be established:

by solving the equations (7) in a similar manner, the operation of the machine tool can be determinedActual coordinate M of measuring point in moving process_i(x_mi,y_mi,z_mi)。

S6, using the space coordinates of a series of measuring points when the translational component moves along the axis obtained in the step S5, and 3 vectors formed by two A, B, C three pointsThe direction of the vector will change in the movement process of the machine tool, the change of the direction of the vector is only related to the angular displacement error in the movement of the machine tool, a series of vector establishing equation sets are formed according to the homogeneous transformation matrix and the measuring points caused by the angular displacement error, and the angular displacement errors can be identified by solving the equation sets. Where the letter A, B, C is the coordinate, angle, or straightness error of the point.

At the initial position, assume measurement point A (x)_a0,y_a0,z_a0)，B(x_b0,y_b0,z_b0)，C(x_c0,y_c0,z_c0) At this time, the corresponding direction vectors are:

when the pallet moves a distance L along the x-axis from the initial position, points A, B, C move to the A ', B ', C ' positions in sequence.

In the motion process, the theoretical homogeneous transformation matrix is as follows:

during the movement of the supporting plate, the angle error epsilon_x(x)、ε_y(x) And ε_z(x) The resulting homogeneous transformation matrices are:

positioning error delta_x(x) Straightness error delta_y(x)、δ_z(x) The resulting homogeneous transformation matrix is:

the overall error transformation matrix is:

P＝T₁T₂T₃T₄T₅T₆ (26)

after the P finishing, the method can be simplified as follows:

when the pallet moves a distance L along the x-axis, the total homogeneous transformation matrix is:

suppose A' (x)_a1,y_a1,z_a1)，B′(x_b1,y_b1,z_b1)，C′(x_c1,y_c1,z_c1) Then the following relationship should exist between the point a' and the point a:

by equation (29), we can calculate the spatial coordinates of the measurement point a'. Similarly, the spatial coordinates of the measurement points B ', C' can be calculated sequentially.

After finishing, the product can be obtained

The spatial coordinates of B 'and C' can be calculated sequentially in the same way.

At this time

Then vector during pallet motionDeviation of direction of

The vector can be found by the formula (32) during the movement of the supporting plateThe change of direction being only related to the angular displacement error epsilon_x(x)，ε_y(x)，ε_z(x) Related to the linear displacement error delta_x(x)，δ_y(x)，δ_z(x) Is irrelevant.

In fig. 2, it is assumed that at the initial position When the moving part of the machine tool moves along the x axis L, new vectors are obtained in sequence For vectorsAnd vectorThe following relation should exist

After finishing, the method can be obtained:

in the same way, vectors can be establishedAndand a vectorAnd vectorThe mutual relationship between the two can be obtained by arranging the following components into a matrix form:

the direction value of each vector in equation (35) can be calculated from the coordinates of the measurement point measured by the laser tracker. And solving the equation set by adopting a least square method, so that each angular motion error can be identified. The method increases the redundant data quantity, thereby improving the identification precision of each error.

S7, determining a series of new points by taking the space coordinates of a series of measuring points when the translational component moves along the axis obtained in the step S5 and the angular displacement errors of each item obtained in the step S6, taking the midpoint of A, B as T, the midpoint of B, C as R and the midpoint of C, A as S, respectively establishing a motion error equation for T, R, S to form a redundancy equation set, and substituting the redundancy equation set into the solution of each angular displacement error, so that each linear displacement error can be identified, and the complete identification of the geometric errors of the numerical control machine tool is completed.

At the initial position, the point coordinate in the measurement points A, B is assumed to be T (x)_t,y_t,z_t). When the moving part of the machine moves a distance L along the x-axis, the point T moves to T' (x)_t′,y_t′,z_t') position, the corresponding motion error is then:

Δx_t＝x′_t-x_t，Δy_t＝y′_t-y_t，Δz_t＝z′_t-z_t (36)

Δ x here_t，Δy_tAnd Δ z_tThe method consists of two parts, wherein one part is displacement errors T on an x axis, a y axis and a z axis caused by angular motion errors of a point T in the motion process_x1，t_y1，t_z1Another part is the linear displacement error during motion, e.g. delta_x(x)，δ_y(x)，δ_z(x) In that respect From the previous identification of the angular movement error epsilon during movement_x(x)，ε_y(x)，ε_z(x) The displacement errors T in the x, y and z axes caused by the angular motion errors when the point T moves to T' can be calculated_x1，t_y1，t_z1. The following set of equations may then be established:

in the same way, motion error equations of points R and C in B and C, and points S in A and motion to R 'and S' can be established, a redundant equation set is formed by the motion error equations and the equation (37), and displacement errors of all the lines can be identified by solving.

In the error identification process, the laser tracker can identify corresponding geometric errors only by measuring the movement of the machine tool along a single axis, and the movement track of the machine tool does not need to be in a 3D space, so that the measurement track is simplified.

During measurement, the precise numerical control rotary table is arranged at a proper position near a main shaft of the turning and milling machining center and translates along the axial direction of the numerical control machine along with the machine tool workbench. Controlling the translational part of the numerical control machine tool to feed along the axis direction of the machine tool along a preset path, wherein the motion area is 450mm multiplied by 550mm multiplied by 350mm, and the laser tracker tracks and measures the motion of the translational part in real time.

When the numerical control machine tool translation component feeds along the axis direction of the machine tool, three measuring points are measured respectively when the numerical control machine tool translation component moves for a certain distance.

When the translation component is fed along the x-axis direction, one measuring point is arranged every 50mm of movement, and the total number of the measuring points is 30. When the translation component moves to the position of each measuring point, the machine tool stops moving, the precision rotary table arranged on the main shaft of the machine tool rotates for 120 degrees in sequence, each position controls the precision rotary table to stop for 5 seconds, and distance measuring data of the laser tracker of the cat eye at different positions on the rotary table are recorded.

As shown in fig. 2, after the laser tracker completes the measurement of each measurement point at the first base station, the laser tracker moves to the next base station and repeats the above measurement process until the measurement of the measurement point is completed at all the base station positions. In order to reduce the influence of random errors on the measurement result in the measurement process, three times of measurement are carried out on the measurement point at each base station position, the total measurement time is about 3 hours, and the measurement efficiency is high. By utilizing the measurement algorithm and the machine tool single-term error separation algorithm deduced above, each single-term geometric error of the X axis of the numerical control machine tool can be identified.

The invention realizes the rapid and accurate calibration of the geometric error of the translational shaft of the machine tool, different measuring points are arranged on the moving part of the machine tool and are sequentially connected to form a space vector, the characteristic that the change of the direction of the space vector formed in the moving process of the machine tool is only related to the angular displacement error of the machine tool shaft and is not related to the linear displacement error is utilized, and the measuring point in the moving process of the translational shaft of the machine tool forms the measurement of the change of the direction of the vector formed by the measuring points by using a laser tracker and adopting a multi-station time-sharing measuring principle, thereby sequentially identifying the angular displacement error and the linear displacement error of the translational shaft of the machine tool. A mathematical model for measuring the geometric errors of the translational axis of the numerical control machine tool based on the space vector direction measurement principle is established, and a laser tracker base station position self-calibration algorithm, a measuring point determination algorithm and a machine tool translational axis single error separation algorithm based on precise numerical control turntable measurement are deduced in sequence. The method simplifies the measuring track of the translational axis of the machine tool, reduces the complexity of error identification, and overcomes the defect that the base station calibration precision is influenced by the error of the machine tool. The method has the advantages of high precision, high efficiency, low cost and simple algorithm.

It will be appreciated by those of ordinary skill in the art that the embodiments described herein are intended to assist the reader in understanding the principles of the invention and are to be construed as being without limitation to such specifically recited embodiments and examples. Those skilled in the art can make various other specific changes and combinations based on the teachings of the present invention without departing from the spirit of the invention, and these changes and combinations are within the scope of the invention.