CN110109418B - A rapid identification method for geometric errors of a large-scale gantry five-sided machining center - Google Patents
A rapid identification method for geometric errors of a large-scale gantry five-sided machining center Download PDFInfo
- Publication number
- CN110109418B CN110109418B CN201910416007.4A CN201910416007A CN110109418B CN 110109418 B CN110109418 B CN 110109418B CN 201910416007 A CN201910416007 A CN 201910416007A CN 110109418 B CN110109418 B CN 110109418B
- Authority
- CN
- China
- Prior art keywords
- error
- axis
- machine tool
- geometric
- errors
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 60
- 239000002245 particle Substances 0.000 claims abstract description 41
- 238000005259 measurement Methods 0.000 claims abstract description 11
- 238000002922 simulated annealing Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 44
- 230000009466 transformation Effects 0.000 claims description 43
- 238000005096 rolling process Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 2
- 238000000691 measurement method Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000005457 optimization Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000011478 gradient descent method Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 206010063385 Intellectualisation Diseases 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 238000000137 annealing Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000005283 ground state Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
- G05B19/404—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by control arrangements for compensation, e.g. for backlash, overshoot, tool offset, tool wear, temperature, machine construction errors, load, inertia
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/35—Nc in input of data, input till input file format
- G05B2219/35222—From cad derive data points for endball mill, grinder, then radius compensation
Landscapes
- Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Numerical Control (AREA)
Abstract
本发明公开了一种大型龙门五面加工中心的几何误差快速辨识方法，先建立机床的几何误差模型，采用球杆仪对机床进行测量分析，将球杆仪的测量数据代入几何误差模型中，得到多项几何误差参数的超定方程组；对超定方程组进行求解，得到各项几何误差值；其特征在于，先采用模拟退火粒子群算法得到群体最优解，再将群体最优解作为LM算法的初值对上述超定方程组进行仿真求解。本发明具有鲁棒性好，群体搜索能力强，求解速度快，稳定性较好，能够避免陷入局部最小值等优点。
The invention discloses a method for quickly identifying geometric errors of a large-scale gantry five-sided machining center. First, a geometric error model of the machine tool is established, a ballbar is used to measure and analyze the machine tool, and the measurement data of the ballbar is substituted into the geometric error model. Obtain over-determined equations of multiple geometric error parameters; solve the over-determined equations to obtain various geometric error values; it is characterized in that the simulated annealing particle swarm algorithm is used to obtain the optimal solution of the group, and then the optimal solution of the group is obtained. As the initial value of the LM algorithm, the above-mentioned overdetermined equations are simulated and solved. The invention has the advantages of good robustness, strong group search ability, fast solution speed, good stability, and can avoid falling into a local minimum value.
Description
Technical Field
The invention relates to the technical field of machine tool geometric error identification, in particular to a method for quickly identifying geometric errors of a large gantry five-surface machining center.
Background
The numerical control machine tool is a working master machine in the manufacturing industry, the development level of the numerical control machine tool represents the national manufacturing industry level, and the rapid speed of the numerical control machine tool develops towards high precision and intellectualization. A number of studies have shown that: the influence of the geometric error and the thermal error of the machine tool on the precision of the machine tool is up to 70 percent, wherein the geometric error accounts for 35 to 70 percent. Therefore, experts and scholars at home and abroad carry out a great deal of research on geometric error modeling, error identification and error compensation of the numerical control machine tool, and remarkable results are obtained. The identification of the geometric errors of the numerical control machine tool is divided into a single geometric error direct measurement method and a comprehensive geometric error identification method. The single geometric error direct measurement method is to directly utilize a measuring instrument to detect various geometric errors of a machine tool, but has low efficiency and poor precision and is difficult to realize automatic measurement; the comprehensive geometric error identification method is to measure a certain working range of a machine tool and identify the comprehensive errors of measuring points by using a machine tool mathematical model, and generally comprises a grating array method, a one-dimensional spherical array method, a circle measurement method, a 9-line method, a 12-line method and the like of a laser interferometer, and the circle measurement method can comprehensively evaluate the machine tool errors. The most typical method in the circular measurement method is a ball rod instrument, and the ball rod instrument has the advantages of simple measurement method, high precision, high measurement speed and the like.
The key point of geometric error identification based on ball bar instrument measurement lies in accurate and rapid solution of an over-determined equation set, the existing ball bar instrument over-determined equation set solving method is a pure particle swarm algorithm which can effectively solve the solution of the equation set, but the algorithm easily falls into a local optimal solution and is slow in solving speed.
Disclosure of Invention
Aiming at the defects of the prior art, the technical problems to be solved by the invention are as follows: how to provide a method for quickly identifying geometric errors of a large gantry five-surface machining center, which has the advantages of good robustness, strong group search capability, high solving speed and good stability and can avoid falling into a local minimum value.
In order to solve the technical problems, the invention adopts the following technical scheme:
a method for quickly identifying geometric errors of a large gantry five-surface machining center comprises the steps of firstly establishing a geometric error model of a machine tool, adopting a ball bar instrument to carry out measurement analysis on the machine tool, substituting measurement data of the ball bar instrument into the geometric error model, and obtaining an overdetermined equation set of multiple geometric error parameters; solving the overdetermined equation set to obtain each geometric error value; the method is characterized in that a group optimal solution is obtained by adopting a simulated annealing particle swarm algorithm, and then the group optimal solution is used as an initial value of an LM algorithm to carry out simulation solution on the overdetermined equation set.
Further, when solving the over-determined equation set, the following steps are specifically adopted:
s1, randomly initializing the position and the speed of each particle in the population;
s2, evaluating the fitness value of each particle, storing the position and the fitness value of each current particle in Pi of each particle, and storing the position and the fitness value of the individual with the optimal fitness value in all Pbest in Pg;
s3, determining the initial temperature;
s4, determining the fitting value of each Pi at the current temperature according to the following formula:
s5, determining a global optimal substitute value Pg from all Pi by using the roulette strategy, and updating the speed and position of each particle by the following formula:
xij(t+1)＝xij(t)+Vij(t+1)
s6, calculating the target value of each particle, and updating the Pi value of each particle and the Pg value of the population;
s7, calculating a temperature reduction operation;
s8, judging whether the iteration reaches a termination condition according to the preset operation precision or program iteration times, if so, stopping searching and outputting the result, and performing subsequent steps, otherwise, repeating the steps S4-S8;
s9, taking the optimal solution Pg of the population as an initial value of an LM algorithm;
s10, let k equal to 0, calculate
If g (x)k) Stopping calculation and outputting xkIs an approximate minimum point;
s11, solving an equation set
Solved to dk；
S12, let mkTo satisfy
M, let α be a minimum non-negative integer ofk＝βmk，xk+1＝xk+αkdk；
S13, update mukK is made k +1, go to step S10;
wherein
c1And c2Is an acceleration factor; r is1And r2Is a random number between 0 and 1; t is the number of iterations; n is the number of particles; beta and sigma are parameters between 0 and 1; mu.s0Greater than 0; ε is the tolerance error.
Further, the geometric error model of the machine tool is established by adopting the following steps:
the geometric errors of the machine tool are first determined, wherein,
the errors in the movement along the X axis are: positioning error deltaxStraightness error delta in the (X) and Y directionsyStraightness error delta in the (X) and Z directionsz(X), roll error εx(X), yaw error εy(X) and Pitch error εz(X)；
The errors of the movement along the Y axis are: positioning error deltay(Y), straightness error delta in X directionx(Y), straightness error delta in Z directionz(Y), rolling error εy(Y), yaw error εx(Y) and pitch error εz(Y)；
The errors of the movement along the Z axis are: positioning error deltaz(Z), straightness error delta in X directionx(Z), straightness error delta in Y directiony(Z), roll error εz(Z), yaw error εx(Z) and Pitch error εy(Z)；
Moving along the W-axisThe errors are respectively: positioning error deltaz(W), straightness error in X direction δx(W), straightness error in Y direction δy(W), rolling error εz(W), yaw error εx(W) and pitch error εy(W)；
Determining X, Y, Z, W perpendicularity errors due to orthogonality between coordinate axes, including perpendicularity error S between the X-axis and the Y-axisxyError of perpendicularity S between Y-axis and Z-axisyzError of perpendicularity S between Z axis and X axiszxError of perpendicularity S between W axis and Y axisywAnd the perpendicularity error S between the W axis and the X axiswx；
Dividing the machine tool into a workpiece branch and a cutter branch, respectively obtaining a position matrix of the workpiece relative to the machine tool and a position matrix of the cutter equivalent to the machine tool, wherein the difference value of the two matrixes is a relative displacement error matrix of the cutter and the workpiece, namely a geometric error matrix of the machine tool.
Further, the position matrix of the workpiece with respect to the machine tool is obtained as follows:
the homogeneous transformation matrix containing geometric errors of the X-axis with respect to the bed R is:
then the homogeneous coordinate transformation matrix containing geometric errors of the bed R to X axis coordinate system is:
if the position coordinate vector of the theoretical processing point of the workpiece in the workpiece coordinate system is as follows:
TB＝[XB YB ZB 1]T
the homogeneous coordinate transformation matrix of the workpiece B relative to the bed R, which contains geometric errors, is:
the position matrix of the tool corresponding to the machine tool is obtained as follows:
fixing the Z axis, and obtaining homogeneous coordinate transformation from the R axis of the lathe bed to the W axis, homogeneous coordinate transformation from the W axis to the Y axis, and homogeneous coordinate transformation from the Y axis to the cutter respectively as follows:
if the position coordinate vector of the machine tool tip in the tool coordinate system is:
Tt＝[0 0 0 1]T
the homogeneous transformation matrix containing the geometric errors of the tool nose T relative to the lathe bed R is as follows:
fixing the W shaft, and obtaining homogeneous coordinate transformation from the bed body R to the Y shaft, homogeneous coordinate transformation from the Y shaft to the Z shaft and homogeneous coordinate transformation from the Y shaft to the cutter respectively as follows:
if the position coordinate vector of the machine tool tip in the tool coordinate system is:
Tt＝[0 0 0 1]T
the homogeneous transformation matrix containing the geometric errors of the tool nose T relative to the lathe bed R is as follows:
homogeneous coordinate transformation matrix for making tool equivalent to machine tool
Subtracting a homogeneous coordinate transformation matrix of the workpiece relative to the machine tool
Obtaining a relative displacement error matrix of the tool nose and the workpiece, namely a geometric error matrix of the machine tool:
to pair
And (5) asking for help, and removing the error cubic term and the high-order term to obtain a simplified geometric error model moving along the X, Y, Z coordinate axis.
In conclusion, the method has the advantages of good robustness, strong group search capability, high solving speed, good stability, capability of avoiding falling into a local minimum value and the like.
Drawings
Fig. 1 is a schematic structural diagram of a five-surface machining center of a large gantry.
FIG. 2 is a graph comparing the verticality error identified by the simplex particle swarm algorithm with the error measured by the ball arm instrument.
FIG. 3 is a graph comparing the verticality error identified by the SAPSO-LM algorithm with the error measured by the ball bar instrument.
FIG. 4 is a graph of compensated front cue stick counter-clockwise running error data.
FIG. 5 is a graph of compensated front cue stick instrument clockwise running error data.
FIG. 6 is a diagram of error data of counter-clockwise operation of the cue stick instrument compensated by a simplex particle swarm optimization.
FIG. 7 is a diagram of the error data of the club instrument running clockwise after compensation by the simplex particle swarm algorithm.
FIG. 8 is a graph of data for counter-clockwise error in the club instrument compensated using the method of the present invention.
FIG. 9 is a graph of data for error in the clockwise operation of a cue stick apparatus compensated using the method of the present invention.
Detailed Description
The present invention will be described in further detail with reference to examples.
In the specific implementation: as shown in figure 1, a large gantry five-surface machining center is composed of a lathe bed 1, an X axis, a Y axis, a Z axis, a W axis, a cutter 3, a workbench, a workpiece 2 and the like. The machine tool structure has two branches, namely a branch from the machine tool body to a workpiece, namely a machine tool body-X-axis-workpiece, and a branch from the machine tool body to a cutter, namely a machine tool body-W-Y-axis-Z-axis-cutter.
1. Geometric error analysis of large gantry five-surface machining center
In the machine tool structure, a moving pair consists of a guide rail and a slide carriage for a translation shaft. Ideally, the slide carriage only reciprocates in a specified direction, but in actual conditions, the slide carriage and the guide rail have errors such as manufacturing, size, assembly and the like. In addition, errors in the geometric dimensions of the guide surfaces of the guide rails, errors in the parallelism of the guide surfaces, differences in the clearance between the guide rails and the slide carriage, and errors in the shape of the rolling bodies all affect the accuracy of the sliding pairs. Thus, moveThe kinematic pair has errors in directions of various degrees of freedom during the motion process. A kinematic pair has six errors, wherein the six errors comprise 3 translation errors and 3 rotation errors. When moving along the X axis, three translational errors are respectively: positioning error deltaxStraightness error delta in the (X) and Y directionsyStraightness error delta in the (X) and Z directionsz(X); the three rotation errors are respectively: roll error epsilonx(X), yaw error εy(X) and Pitch error εz(X). As can be seen from fig. 1, the W axis specific to the large gantry five-surface machining center is a vertical motion, and the same can be seen: the six errors for the Y axis are: positioning error deltay(Y), straightness error delta in X directionx(Y), straightness error delta in Z directionz(Y), rolling error εy(Y), yaw error εx(Y) and pitch error εz(Y); the six errors for the Z axis are: positioning error deltaz(Z), straightness error delta in X directionx(Z), straightness error delta in Y directiony(Z), roll error εz(Z), yaw error εx(Z) and Pitch error εy(Z); the six errors for the W axis are: positioning error deltaz(W), straightness error in X direction δx(W), straightness error in Y direction δy(W), rolling error εz(W), yaw error εx(W) and pitch error εy(W)。
There are also 5 perpendicularity errors, which refer to errors between the X, Y, Z, W coordinate axes due to orthogonality. In an ideal state, coordinate axes are mutually perpendicular in pairs (except for a Z axis and a W axis), and the two perpendicular axes are not completely perpendicular due to factors such as assembly, so that small-angle deviation is caused. SxyIndicating the error in perpendicularity between the X-axis and the Y-axis, SyzIndicating the error in perpendicularity between the Y axis and the Z axis, SzxIndicating the error in perpendicularity between the Z axis and the X axis, SywIndicating the perpendicularity error between the W and Y axes, SwxIndicating the perpendicularity error between the W axis and the X axis.
2. Large-scale gantry five-surface machining center geometric error rapid identification method combining Simulated Annealing Particle Swarm Optimization (SAPSO) and Levenberg-Marquardt (L-M) algorithm
2.1 modeling of geometric error of large gantry five-surface machining center based on multi-body system
Establishing a universal geometric error model of the machine tool, dividing the machine tool into a workpiece branch and a cutter branch on the basis of a multi-body system theory and a secondary coordinate transformation theory to obtain a position matrix of the cutter relative to the machine tool body and a position matrix of the workpiece relative to the machine tool body, and respectively establishing the position matrix of the cutter relative to the workpiece under an ideal condition and an error state to obtain the universal error model of the machine tool.
2.1.1 homogeneous coordinate transformation of workpiece kinematic chain
(1) Homogeneous coordinate transformation from bed R to X axis
When the X axis moves on the machine tool, 6 geometric errors exist: deltax(X)、δy(X)、δz(X)、εx(X)、εy(X)、εz(X), then the homogeneous transformation matrix containing geometric errors of the X-axis with respect to the bed R is:
then the homogeneous coordinate transformation matrix containing geometric errors of the bed R to X axis coordinate system is:
setting the position coordinate vector of the theoretical processing point of the workpiece in the workpiece coordinate system as follows: t isB＝[XB YB ZB 1]T. The homogeneous coordinate transformation matrix of the workpiece B with respect to the machine bed R, containing geometrical errors, is then:
2.1.2 homogeneous coordinate transformation of tool kinematic chain
The Z axis and the W axis of the large gantry five-surface machining center move in the vertical direction, so that the Z axis and the W axis are respectively fixed for modeling.
(1) Fixed Z axis
The homogeneous coordinate transformation from the lathe bed R to the W axis, the homogeneous coordinate transformation from the W axis to the Y axis and the homogeneous coordinate transformation from the Y axis to the cutter can be respectively as follows:
assuming that the position coordinate vector of the machine tool tip in the tool coordinate system is: t ist＝[0 0 0 1]TThen the homogeneous transformation matrix containing geometric errors of the tool nose T relative to the bed R is:
(2) fixed W shaft
The homogeneous coordinate transformation from the bed body R to the Y axis, the homogeneous coordinate transformation from the Y axis to the Z axis and the homogeneous coordinate transformation from the Y axis to the cutter can be respectively as follows:
assuming that the position coordinate vector of the machine tool tip in the tool coordinate system is: t ist＝[0 0 0 1]TThen the homogeneous transformation matrix containing geometric errors of the tool nose T relative to the bed R is:
2.1.3 geometric error model of large gantry five-surface machining center
Obtaining a homogeneous coordinate transformation matrix of a workpiece relative to a machine tool
Homogeneous coordinate transformation matrix equivalent to machine tool with cutter
And then, the difference value of the two matrixes is a relative displacement error matrix of the tool nose and the workpiece, namely a geometric error matrix of the machine tool:
solving for
When the error cubic term and the higher terms are omitted, a geometric error model which is simplified and moves along the X, Y, Z coordinate axis can be obtained.
2.2 error rapid identification of five-face machining center of large gantry
2.2.1 error identification principle based on ball arm instrument
The research of the geometric error parameter identification method is an important content of the geometric error detection of the numerical control machine tool, and a large amount of research is carried out by scholars at home and abroad. At present, the main methods for identifying relatively mature geometric error parameters include: 9-line method, 12-line method and club instrument identification method.
The ball bar instrument is standard equipment for machine tool precision analysis, and utilizes two shafts of a machine tool to perform circular interpolation by means of a high-precision displacement sensor, analyzes the radius change of a circular arc and the track characteristics of an arc line, and analyzes data by analysis software. The machine tool is measured and analyzed by using a QC20 ball bar instrument of Renishaw, and the QC20 ball bar instrument is arranged on a workbench, and care should be taken to avoid collision during installation, otherwise, the measurement precision is seriously influenced. Adjusting the room temperature to the working temperature of the machine tool, and respectively making circular arc motion in an X-Y, X-Z, Y-Z, X-W, Y-W plane in an unloaded state.
The XY, XZ, YZ, XW, YW plane radial error data are obtained as shown in Table 1.
TABLE 1 XY, XZ, YZ, XW, YW plane radial error values
And substituting the measured data of the ball arm instrument into a geometric error model of the machine tool to obtain an overdetermined equation set about 21 geometric error parameters. The overdetermined equation set refers to an equation set in which the number of the equation sets is more than the number of unknowns. The overdetermined system of equations generally has no exact solution, and in most cases, an approximate solution in some sense is solved. Each geometric error value can be calculated by solving the over-determined equation set. The key point of geometric error identification based on ball bar instrument measurement lies in accurate and rapid solution of an over-determined equation set, the existing ball bar instrument over-determined equation set solving method is a pure particle swarm algorithm which can effectively solve the solution of the equation set, but the algorithm easily falls into a local optimal solution and is slow in solving speed.
2.2.2 geometric error identification algorithm of numerical control machine tool
The particle swarm optimization algorithm is an evolutionary computing technology and is derived from the research on the behavior of predation of birds and animals. The particle updates its velocity and new position according to the following formula.
xij(t+1)＝xij(t)+Vij(t+1)
The particle swarm algorithm has the advantages of memorability, high searching speed, few parameters, simple structure and the like, but is easy to fall into local optimum, so that the convergence precision is low and the convergence is difficult.
The simulated annealing algorithm is a general probability algorithm and is used for finding the optimal solution of the problem in a large search space. The idea is derived from the annealing process of the solid, i.e. the solid is heated to a sufficiently high temperature and then slowly cooled, and finally reaches the ground state at normal temperature, with the internal energy minimized. The simulated annealing algorithm has the advantages of simple calculation process, strong universality and robustness, and can be used for solving the complex nonlinear optimization problem but has the defects of low convergence speed, long execution time and the like.
The L-M algorithm is a combination of the gradient descent method and the Gauss-Newton method. When mu is increased, the algorithm is similar to a gradient descent method, and the global characteristic of the algorithm is exerted; when mu is reduced, the algorithm approaches the Gauss-Newton method and exerts the local convergence characteristic. The L-M algorithm adopts approximate second derivative information, so that the required iteration time is short, the convergence is very quick, the stability of the algorithm is good, and the local minimum value is avoided. The iterative formula is as follows:
μk+1＝μk-(Ak TAk+μl)-1Ak Tek
the L-M algorithm has local convergence of the Gaussian Newton method and the global characteristic of the gradient descent method, but has extremely high requirements on initial values.
Aiming at the advantages and disadvantages of the algorithm, the simulated annealing particle swarm algorithm is combined with the L-M algorithm, and the advantages of the simulated annealing particle swarm algorithm and the L-M algorithm are combined, so that the hybrid algorithm for solving the nonlinear equation set is provided. The mixed algorithm component exerts the group searching capability of the particle swarm algorithm, the simulated annealing algorithm has the probability jump capability in the searching process, and the local minimum solution and the local fine searching performance of the L-M algorithm can be effectively avoided in the searching process.
The algorithm steps combined with the above requirements are as follows:
s1, randomly initializing the position and the speed of each particle in the population;
s2, evaluating the fitness of each particle, storing the position and the fitness value of each current particle in Pi of each particle, and storing the position and the fitness value of the individual with the optimal fitness value in all Pbest in Pg;
s3, determining the initial temperature, specifically:
s4, determining the fitting value of each Pi at the current temperature according to the following formula:
s5, determining a global optimal substitute value Pg from all Pi by using the roulette strategy, and updating the speed and position of each particle by the following formula:
xij(t+l)＝xij(t)+Vj(t+1)
s6, calculating the target value of each particle, and updating the Pi value of each particle and the Pg value of the population;
s7, calculating the temperature reduction operation, specifically Tk+1＝CTkC belongs to (0.5,0.99), and the value determines the cooling process;
s8, if the stop condition is met, stopping searching and outputting the result, and performing subsequent steps, otherwise, repeating the steps S4-S8;
s9, taking the optimal solution Pg of the population as an initial value of an LM algorithm;
s10, let k equal to 0, calculate
If g (x)k) Stopping calculation and outputting xkIs an approximate minimum point;
s11, solving an equation set
Solved to dk；
S12, let mkTo satisfy
M, let α be a minimum non-negative integer ofk＝βmk，xk+1＝xk+αkdk；
S13, update mukK is made k +1, go to step S10;
wherein
c1And c2Is an acceleration factor; r is1And r2Is a random number between 0 and 1; t is the number of iterations; n is the number of particles; beta and sigma are parameters between 0 and 1; mu.s0Greater than 0; ε is the tolerance error.
3. Experiments and analyses
3.1 geometric error of five-face machining center of large gantry
Various geometric error coefficients of the large gantry five-surface machining center can be obtained by writing an MATLAB program by combining simulated annealing particle swarm with an L-M algorithm to solve an overdetermined equation set as shown in a table 2.
TABLE 2 geometric error of five-surface machining center of large gantry solved by SAPSO-LM algorithm
And (3) comparison and analysis: the verticality error in the geometric error of the large-scale gantry pentahedral machining center can be obtained in the measurement of the ball bar instrument QC20, and the results of the simplex particle swarm algorithm and the results of the SAPSO-LM algorithm are respectively compared and shown in the table 3.
TABLE 3 comparison of perpendicularity error obtained by two algorithms with the error measured by a ball arm instrument
From table 3, fig. 2 and fig. 3 can be obtained, and it is obvious that the verticality error identified by the SAPSO-LM algorithm is smaller than the difference of the simplex particle swarm algorithm.
3.2 post-compensation experimental verification
And respectively substituting the identification result of the simplex particle swarm algorithm and the identification result of the SAPSO-LM into the error compensation model to respectively obtain the compensated numerical control program. And then respectively substituting three programs (a pre-compensation program, a compensation program after simplex particle swarm algorithm identification and an SAPSO-LM algorithm identification compensation program) into the machine tool and operating to obtain a ball bar instrument result.
The measured results of the geometrical error compensation front ball arm instrument of the large gantry five-surface machining center are shown in figure 4, figure 5 and table 4,
TABLE 4 Compensation of Forward cue Meter operational error data
The geometrical errors of the large gantry five-surface machining center are compensated by a simplex particle swarm algorithm, and the measured results of the ball arm instrument are shown in figure 6, figure 7 and table 5,
TABLE 5 simplex particle swarm algorithm compensated post-cue instrument operation error data
The measured results of the ball arm instrument after the geometric errors of the large gantry five-surface machining center are compensated by adopting the SAPSO-LM algorithm are shown in figure 8, figure 8 and table 6,
TABLE 6 SAPSO-LM Algorithm compensated post-ball rod instrument operational error data
Comparing fig. 4, fig. 5, fig. 6, fig. 7, fig. 8 and fig. 9, it can be known that the simplex particle swarm algorithm can effectively compensate the geometric error of the numerical control machine, but the compensation effect under the identification of the SAPSO-LM algorithm is obviously better than that of the simplex particle swarm algorithm.
4. Conclusion
In order to improve the geometric precision of the large-scale gantry five-surface machining center, a geometric error model of a machine tool is established, a geometric error coefficient of the machine tool is identified by an SAPSO-LM algorithm, and an error compensation model of the machine tool is established. The effectiveness of the algorithm is determined and is superior to the simplex particle swarm algorithm through the dynamic performance test of the ball arm instrument.
The above description is only exemplary of the present invention and should not be taken as limiting, and any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (1)
1. A method for quickly identifying geometric errors of a large gantry five-surface machining center comprises the steps of firstly establishing a geometric error model of a machine tool, adopting a ball bar instrument to carry out measurement analysis on the machine tool, substituting measurement data of the ball bar instrument into the geometric error model, and obtaining an overdetermined equation set of multiple geometric error parameters; solving the overdetermined equation set to obtain each geometric error value; the method is characterized in that a group optimal solution is obtained by adopting a simulated annealing particle swarm algorithm, and then the group optimal solution is used as an initial value of an LM algorithm to carry out simulation solution on the overdetermined equation set;
the geometric error model of the machine tool is established by adopting the following steps:
the geometric errors of the machine tool are first determined, wherein,
the errors in the movement along the X axis are: positioning error deltaxStraightness error delta in the (X) and Y directionsyStraightness error delta in the (X) and Z directionsz(X), roll error εx(X), yaw error εy(X) and Pitch error εz(X)；
The errors of the movement along the Y axis are: positioning error deltay(Y), straightness error delta in X directionx(Y), straightness error delta in Z directionz(Y), rolling error εy(Y), yaw error εx(Y) and pitch error εz(Y)；
The errors of the movement along the Z axis are: positioning error deltaz(Z), straightness error delta in X directionx(Z), straightness error delta in Y directiony(Z), roll error εz(Z), yaw error εx(Z) and Pitch error εy(Z)；
The errors of the movement along the W axis are: positioning error deltaz(W), straightness error in X direction δx(W), straightness error in Y direction δy(W), rolling error εz(W), yaw error εx(W) and pitch error εy(W)；
Determining X, Y, Z, W perpendicularity errors due to orthogonality between coordinate axes, including perpendicularity error S between the X-axis and the Y-axisxyError of perpendicularity S between Y-axis and Z-axisyzError of perpendicularity S between Z axis and X axiszxError of perpendicularity S between W axis and Y axisywAnd the perpendicularity error S between the W axis and the X axiswx；
Dividing the machine tool into a workpiece branch and a cutter branch, respectively obtaining a position matrix of the workpiece relative to the machine tool and a position matrix of the cutter equivalent to the machine tool, wherein the difference value of the position matrix and the position matrix is a relative displacement error matrix of the cutter and the workpiece, namely a geometric error matrix of the machine tool;
the position matrix of the workpiece relative to the machine tool is obtained as follows:
the homogeneous transformation matrix containing geometric errors of the X-axis with respect to the bed R is:
then the homogeneous coordinate transformation matrix containing geometric errors of the bed R to X axis coordinate system is:
if the position coordinate vector of the theoretical processing point of the workpiece in the workpiece coordinate system is as follows:
TB＝[XB YB ZB 1]T
the homogeneous coordinate transformation matrix of the workpiece B relative to the bed R, which contains geometric errors, is:
the position matrix of the tool corresponding to the machine tool is obtained as follows:
fixing the Z axis, and obtaining homogeneous coordinate transformation from the R axis of the lathe bed to the W axis, homogeneous coordinate transformation from the W axis to the Y axis, and homogeneous coordinate transformation from the Y axis to the cutter respectively as follows:
if the position coordinate vector of the machine tool tip in the tool coordinate system is:
Tt＝[0 0 0 1]T
the homogeneous transformation matrix containing the geometric errors of the tool nose T relative to the lathe bed R is as follows:
fixing the W shaft, and obtaining homogeneous coordinate transformation from the bed body R to the Y shaft, homogeneous coordinate transformation from the Y shaft to the Z shaft and homogeneous coordinate transformation from the Y shaft to the cutter respectively as follows:
if the position coordinate vector of the machine tool tip in the tool coordinate system is:
Tt＝[0 0 0 1]T
the homogeneous transformation matrix containing the geometric errors of the tool nose T relative to the lathe bed R is as follows:
make the cutting tool looksHomogeneous coordinate transformation matrix of machine tool
Subtracting a homogeneous coordinate transformation matrix of the workpiece relative to the machine tool
Obtaining a relative displacement error matrix of the tool nose and the workpiece, namely a geometric error matrix of the machine tool:
to pair
And (5) asking for help, and removing the error cubic term and the high-order term to obtain a simplified geometric error model moving along the X, Y, Z coordinate axis.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910416007.4A CN110109418B (en) | 2019-05-19 | 2019-05-19 | A rapid identification method for geometric errors of a large-scale gantry five-sided machining center |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910416007.4A CN110109418B (en) | 2019-05-19 | 2019-05-19 | A rapid identification method for geometric errors of a large-scale gantry five-sided machining center |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110109418A CN110109418A (en) | 2019-08-09 |
| CN110109418B true CN110109418B (en) | 2021-04-13 |
Family
ID=67490913
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910416007.4A Active CN110109418B (en) | 2019-05-19 | 2019-05-19 | A rapid identification method for geometric errors of a large-scale gantry five-sided machining center |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110109418B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111459095B (en) * | 2020-04-08 | 2021-02-26 | 新代科技(苏州)有限公司 | Method for eliminating machining reversing lines of numerical control machine tool by using ball rod instrument |
| CN111580459B (en) * | 2020-06-01 | 2021-04-13 | 山东大学 | Five-axis machining center assembly error analysis and control method |
| CN112008492B (en) * | 2020-07-29 | 2022-04-08 | 成都飞机工业（集团）有限责任公司 | Method for identifying perpendicularity error of translational shaft of gantry numerical control machine tool |
| CN112346408B (en) * | 2020-09-30 | 2023-11-24 | 北京工业大学 | A geometric error compensation method for CNC grinding machines by constructing a system of linear equations of two variables |
| CN112621385B (en) * | 2020-12-11 | 2022-03-08 | 西安交通大学 | Rapid identification method for geometric errors of linear axis based on 8 displacement line measurements |
| CN112558547B (en) * | 2021-02-19 | 2021-06-08 | 成都飞机工业（集团）有限责任公司 | Quick optimization method for geometric error compensation data of translational shaft of five-axis numerical control machine tool |
| CN115179445B (en) * | 2022-08-09 | 2023-06-27 | 永得利科技（无锡）有限公司 | CNC (computer numerical control) processing method for vacuum cavity sealing piece based on space positioning error model |
| CN115562151A (en) * | 2022-10-13 | 2023-01-03 | 四川普什宁江机床有限公司 | Curve control method based on numerical control machine tool machining guide rail surface |
| CN117348518B (en) * | 2023-12-05 | 2024-03-01 | 广东海洋大学 | A method and device for geometric error identification of CNC machine tools based on the improved nine-line method |
| CN118616924B (en) * | 2024-08-09 | 2024-12-10 | 比亚迪股份有限公司 | Workpiece positioning correction method, correction device, correction system, and storage medium |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006350524A (en) * | 2005-06-14 | 2006-12-28 | Nissan Motor Co Ltd | Process necessity decision apparatus |
| CN103034166A (en) * | 2012-11-26 | 2013-04-10 | 北京工业大学 | Recognition method of critical geometrical error source of machine tool |
| CN103878641B (en) * | 2014-03-14 | 2016-07-06 | 浙江大学 | The rotating shaft geometric error discrimination method that a kind of five-axle number control machine tool is general |
| CN105479268B (en) * | 2016-01-22 | 2019-01-15 | 清华大学 | Five-axle number control machine tool swinging axle geometric error discrimination method based on RTCP |
| CN106141814B (en) * | 2016-08-04 | 2017-12-29 | 天津大学 | The detection of Digit Control Machine Tool translation shaft geometric error and discrimination method based on LaserTRACER |
| CN107066721B (en) * | 2017-04-09 | 2023-09-12 | 北京工业大学 | A multi-axis CNC machine tool C-axis geometric error measurement system and parameter identification method |
| CN107186548B (en) * | 2017-06-08 | 2018-12-18 | 大连理工大学 | A kind of five-axle number control machine tool rotating shaft geometric error detection method |
| CN107081638A (en) * | 2017-06-15 | 2017-08-22 | 天津大学 | A kind of lathe work Space Thermal error identification method for establishing model |
| CN108340210B (en) * | 2018-05-09 | 2019-10-29 | 天津工业大学 | A kind of gang tool geometric error discrimination method based on ball bar measurement |
| CN108655827B (en) * | 2018-07-30 | 2020-09-25 | 天津中德应用技术大学 | Spatial Error Identification Method of Five-axis CNC Machine Tool |
-
2019
- 2019-05-19 CN CN201910416007.4A patent/CN110109418B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110109418A (en) | 2019-08-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110109418B (en) | A rapid identification method for geometric errors of a large-scale gantry five-sided machining center | |
| CN106141814B (en) | The detection of Digit Control Machine Tool translation shaft geometric error and discrimination method based on LaserTRACER | |
| CN106052556B (en) | A kind of three coordinate measuring machine spatial domain coordinates compensation method | |
| CN102699761B (en) | Error identification method of five-axis numerically controlled machine tool based on S-shaped test specimen | |
| CN112008492B (en) | Method for identifying perpendicularity error of translational shaft of gantry numerical control machine tool | |
| CN102785129B (en) | The online test method of the surface machining accuracy of complex parts | |
| CN102062575B (en) | Method for detecting geometric accuracy of numerically-controlled machine tool based on multi-channel laser time-sharing measurement | |
| CN101982726A (en) | Method for detecting geometrical motion error of triaxial numerical control equipment | |
| CN113917888B (en) | Machining precision improving method based on fixed angular calibration and compensation | |
| Hansel et al. | Improving CNC machine tool geometric precision using manufacturing process analysis techniques | |
| Han et al. | A review of geometric error modeling and error detection for CNC machine tool | |
| CN113960959B (en) | A Traceability Method for Key Geometric Errors of Multi-axis CNC Machine Tool | |
| Wu et al. | Sensitivity analysis of geometric errors of two-turntable five-axis machine tool based on S-shaped specimens | |
| CN113985812A (en) | A Method for Predicting Machining Errors of Multi-axis CNC Machine Tools | |
| Merghache et al. | Numerical evaluation of geometrical errors of three-axes CNC machine tool due to cutting forces—case: milling | |
| Chen et al. | An accurate characterization method to tracing the geometric defect of the machined surface for complex five-axis machine tools | |
| CN113467371B (en) | R-test-based five-axis machine tool RTCP parameter calibration method | |
| CN116256186A (en) | A method and system for testing the dynamic performance of a five-axis machine tool | |
| CN115781416A (en) | Method for quickly measuring and compensating geometric errors of machine tool | |
| Erkan et al. | A cluster analysis applied to volumetric errors of five-axis machine tools obtained by probing an uncalibrated artefact | |
| Hu et al. | A parameter dimension reduction-based estimation approach to enhance the kinematic accuracy of a parallel hardware-in-the-loop docking simulator | |
| Kochetkov et al. | Kinematic error modeling of delta 3D printer | |
| Xu et al. | Dynamic synchronous motion accuracy measurement and estimation for a five-axis mirror milling system | |
| CN107038269B (en) | Numerical control machining machine tool optimization method based on X-shaped aviation thin-wall standard test piece | |
| CN113932707A (en) | Method for compensating geometric error of three-coordinate measuring machine based on elastic network algorithm |
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 | ||
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: No.69 Hongguang Avenue, Banan District, Chongqing Patentee after: Chongqing University of Technology Patentee after: Chongqing GAOJIN Industrial Co.,Ltd. Address before: No.69 Hongguang Avenue, Banan District, Chongqing Patentee before: Chongqing University of Technology Patentee before: CHONGQING GAOKIN INDUSTRY Co.,Ltd. |