Numerical control device and numerical control meihod for machine tool

Abstract

A numerical control device for a machine tool controls a machine tool having a main spindle for attaching a tool, a table holding a workpiece and a jig, three translational axes, and one or more rotation axis. The numerical control device includes an axis-dependent deformation error estimation unit, an input unit, a gravitational deformation estimation unit, a correction value calculation unit, and an addition unit. The correction value calculation unit calculates a correction value of the translational axes and/or the rotation axis with respect to an error of a position and/or a posture of the tool with respect to the workpiece, based on an estimated value of an axis-dependent deformation error, an estimated value of a gravitational deformation error, and command values. The addition unit adds the correction values to the command values.

Claims

1. A numerical control device for a machine tool, the numerical control device controlling a machine tool having a main spindle for attaching a tool, a table holding a workpiece and a jig, three translational axes, and one or more rotation axes, the numerical control device comprising: an axis-dependent deformation error estimation unit that calculates an estimated value of an axis-dependent deformation error caused by a deformation of the machine tool due to operation of at least one of the translational axes and the rotation axes, based on at least an axis-dependent deformation error parameter and command values of the at least one of the translational axes and the rotation axes; an input unit that receives a mass and a barycentric position of at least one of the workpiece and the jig; a gravitational deformation estimation unit that calculates an estimated value of a gravitational deformation error caused by the mass of the at least one of the workpiece and the jig based on at least a gravitational deformation error parameter, the mass and the barycentric position received at the input unit, and the command values; a correction value calculation unit that calculates correction values of the at least one of the translational axes and the rotation axes with respect to an error of a position and a posture of the tool with respect to the workpiece, based on at least one of the estimated value of the axis-dependent deformation error, the estimated value of the gravitational deformation error, and the command values; and an addition unit that adds the correction values to the command values, wherein the axis-dependent deformation error estimation unit calculates the estimated value of the axis-dependent deformation error using a function of a geometric error of six degrees of freedom of the rotation axes taking the command value of the rotation axes as a variable.

2. The numerical control device for the machine tool according to claim 1, wherein the gravitational deformation estimation unit calculates the estimated value of the gravitational deformation error using a function depending on the mass of the at least one of the workpiece and the jig, the barycentric position of the at least one of the workpiece and the jig, and the command value of the rotation axes.

3. The numerical control device for the machine tool according to claim 1, the numerical control device further comprising an error identification unit that measures a position of a target mounted on any one of the main spindle and the table using a measurement unit mounted on another of the main spindle and the table and identifies the axis-dependent deformation error parameter based on three-dimensional coordinate values of a plurality of the measured positions of the target, wherein the axis-dependent deformation error estimation unit calculates the estimated value of the axis-dependent deformation error based on the axis-dependent deformation error parameter identified in the error identification unit and the command value.

4. The numerical control device for the machine tool according to claim 3, wherein the error identification unit takes at least one of the rotation axes as a motion axis and another axis of the rotation axes as a fixed axis, positions the fixed axis at two or more fixed positions, positions the motion axis at each of the respective fixed positions, measures a position of the target with the measurement unit to obtain respective coordinate values of the position of the target, and identifies the axis-dependent deformation error parameter based on the plurality of coordinate values at the respective fixed positions.

5. The numerical control device for the machine tool according to claim 4, wherein the error identification unit identifies an error of six degrees of freedom regarding the motion axis at the respective fixed positions from the plurality of coordinate values at the respective fixed positions, and identifies the axis-dependent deformation error parameter from the identified error of six degrees of freedom.

6. The numerical control device for the machine tool according to claim 3, wherein the error identification unit calculates an estimated value of a position error of the tool with respect to the workpiece caused by the gravitational deformation error from a difference between a position of the tool with respect to the workpiece estimated when the gravitational deformation error exists and an ideal position of the tool with respect to the workpiece without the gravitational deformation error, and identifies the axis-dependent deformation error parameter using a coordinate value obtained by subtracting the estimated value of the position error from a coordinate value of the target measured by the measurement unit.

7. A numerical control device for a machine tool, the numerical control device controlling a machine tool having a main spindle for attaching a tool, a table holding a workpiece and a jig, three translational axes, and one or more rotation axes, the numerical control device comprising: an axis-dependent deformation error estimation unit that calculates an estimated value of an axis-dependent deformation error caused by a deformation of the machine tool due to operation of at least one of the translational axes and the rotation axes, based on at least an axis-dependent deformation error parameter and command values of the at least one of the translational axes and the rotation axes; an input unit that receives a mass and a barycentric position of at least one of the workpiece and the jig; a gravitational deformation estimation unit that calculates an estimated value of a gravitational deformation error caused by the mass of the at least one of the workpiece and the jig based on a gravitational deformation error parameter, the mass and the barycentric position received at the input unit, and the command values; a correction value calculation unit that calculates correction values of the at least one of the translational axes and the rotation axes with respect to an error of a position and a posture of the tool with respect to the workpiece, based on at least one of the estimated value of the axis-dependent deformation error, the estimated value of the gravitational deformation error, and the command values; an addition unit that adds the correction values to the command values; and an error identification unit that measures a position of a target mounted on any one of the main spindle and the table using a measurement unit mounted on another of the main spindle and the table and identifies the axis-dependent deformation error parameter based on three-dimensional coordinate values of a plurality of the measured positions of the target, wherein the axis-dependent deformation error estimation unit calculates the estimated value of the axis-dependent deformation error based on the axis-dependent deformation error parameter identified in the error identification unit and the command value.

8. The numerical control device for the machine tool according to claim 7, wherein the gravitational deformation estimation unit calculates the estimated value of the gravitational deformation error using a function depending on the mass of the at least one of the workpiece and the jig, the barycentric position of the at least one of the workpiece and the jig, and the command value of the rotation axes.

9. The numerical control device for the machine tool according to claim 7, wherein the error identification unit takes at least one of the rotation axes as a motion axis and another axis of the rotation axes as a fixed axis, positions the fixed axis at two or more fixed positions, positions the motion axis at each of the respective fixed positions, measures a position of the target with the measurement unit to obtain respective coordinate values of the position of the target, and identifies the axis-dependent deformation error parameter based on the plurality of coordinate values at the respective fixed positions.

10. The numerical control device for the machine tool according to claim 9, wherein the error identification unit identifies an error of six degrees of freedom regarding the motion axis at the respective fixed positions from the plurality of coordinate values at the respective fixed positions, and identifies the axis-dependent deformation error parameter from the identified error of six degrees of freedom.

11. The numerical control device for the machine tool according to claim 7, wherein the error identification unit calculates an estimated value of a position error of the tool with respect to the workpiece caused by the gravitational deformation error from a difference between a position of the tool with respect to the workpiece estimated when the gravitational deformation error exists and an ideal position of the tool with respect to the workpiece without the gravitational deformation error, and identifies the axis-dependent deformation error parameter using a coordinate value obtained by subtracting the estimated value of the position error from a coordinate value of the target measured by the measurement unit.

12. A numerical control method for controlling a machine tool having a main spindle for attaching a tool, a table holding a workpiece and a jig, three translational axes, and one or more rotation axes, the numerical control method comprising: calculating an estimated value of an axis-dependent deformation error caused by a deformation of the machine tool due to operation of at least one of the translational axes and the rotation axes, based on at least an axis-dependent deformation error parameter and command values of the at least one of the translational axes and the rotation axes; obtaining a mass and a barycentric position of at least one of the workpiece and the jig; calculating an estimated value of a gravitational deformation error caused by the mass of the at least one of the workpiece and the jig, based on at least a gravitational deformation error parameter, the mass and the barycentric position obtained by the obtaining, and the command values; calculating correction values of the at least one of the translational axes and the rotation axes with respect to an error of at least one of a position and a posture of the tool with respect to the workpiece, based on the estimated value of the axis-dependent deformation error, the estimated value of the gravitational deformation error, and the command values; and adding the correction values to the command values, wherein the axis-dependent deformation error estimation unit calculates the estimated value of the axis-dependent deformation error using a function of a geometric error of six degrees of freedom of the rotation axes taking the command value of the rotation axes as a variable.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a schematic diagram of a five-axis-control machining center.

(2) FIG. 2 is a block diagram of a conventional numerical control device.

(3) FIG. 3 is a block diagram of a numerical control device that performs a numerical control method of the disclosure.

(4) FIG. 4 is a flowchart to identify a geometric error.

DETAILED DESCRIPTION OF THE INVENTION

(5) The following describes embodiments of the disclosure based on the drawings.

(6) In an exemplary embodiment of the disclosure, FIG. 1 illustrates a five-axis-control machining center (hereinafter referred to as a “five-axis machine”) as an example. This time, an embodiment for a trunnion table type machining center whose axis configuration from a workpiece to a tool is a C-axis, A-axis, Y-axis, X-axis, and Z-axis is taken. However, the disclosure does not limit the machine tool to a machining center and a five-axis machine insofar as a machine tool having three translational axes and at least one rotation axis.

(7) Next, a description will be given of a geometric error.

(8) The geometric error is defined to have total six components (δx, δy, δz, α, β, and γ) of relative translational errors in three directions and relative rotational errors in three directions between respective axes. In the five-axis machine in this example, a connection of the axes from a workpiece (a material to be worked) fixed to the table 3 to a tool attached to a main spindle of the main spindle head 2 has an order of the C-axis, the A-axis, the Y-axis, the X-axis, and the Z-axis. There are a total of 36 pieces of geometric errors, considering between the Z-axis and the tool, and between the workpiece and the C-axis. However, there are a plurality of ones having redundant relationships among the 36 pieces of geometric errors. Thus, final geometric errors with these ones removed are total 13 pieces: δx5, δy5, α5, β5, δy4, δz4, β4, γ4, γ3, α2, β2, α1, and β1 when the order from a side of the tool of the respective geometric errors is represented as subscripts. These respective geometric errors mean the X-direction error of the C-axis center position, the offset error between the C and A axes, the angular offset error of the A-axis, the squareness between the C and A axes, the Y-direction error of the A-axis center position, the Z-direction error of the A-axis center position, the squareness between the A and X axes, the squareness between the A and Y axes, the squareness between the X and Y axes, the squareness between the Y and Z axes, the squareness between the Z and X axes, the squareness between the main spindle and the Y-axis, and the squareness between the main spindle and the X-axis.

(9) FIG. 2 illustrates an exemplary conventional numerical control device.

(10) After a machining program 11 is input, command value generation means 12 generates a command value of each drive axis.

(11) Next, correction value calculation means 16 calculates a correction value of each axis based on the generated command value. Then, a total value of the command value and the correction value is transmitted to servo command value transformation means 13 to calculate a servo command value. The servo command values of the respective axes are transmitted to servo amplifiers 14a to 14e of the respective axes. The servo amplifiers 14a to 14e of the respective axes drive servo motors 15a to 15e respectively to control a relative position and posture of the main spindle head 2 with respect to the table 3.

(12) Here, a description will be given of a calculation method of the correction value with respect to the geometric error. In order to transform a tool center point vector P.sub.T in a main spindle coordinate system on the main spindle head 2 into a vector in a workpiece coordinate system on the table 3, when assuming that a length of a tool to be used is t, and command positions of the X, Y, Z, A, and C axes are x, y, z, a, and c respectively, the tool center vector in the workpiece coordinate system can be obtained by performing a homogeneous coordinate transformation using the following formula 1. That is, a tool center point vector P.sub.1 in the workpiece coordinate system without the geometric error is obtained.

(13) $[Formula 1]$ $P_{1} = M_{5} .Math. M_{4} .Math. M_{3} .Math. M_{2} .Math. M_{1} .Math. P_{T}$ $where$ $M_{5} = [\begin{matrix} \cos c & - \sin c & 0 & 0 \\ \sin c & \cos c & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{matrix}], M_{4} = [\begin{matrix} 1 & 0 & 0 & 0 \\ 0 & \cos a & - \sin a & 0 \\ 0 & \sin a & \cos a & 0 \\ 0 & 0 & 0 & 1 \end{matrix}], M_{3} = [\begin{matrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & y \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{matrix}], M_{2} [\begin{matrix} 1 & 0 & 0 & x \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{matrix}], M_{1} = [\begin{matrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & z \\ 0 & 0 & 0 & 1 \end{matrix}], P_{T} = [\begin{matrix} 0 \\ 0 \\ - t \\ 1 \end{matrix}]$

(14) On the other hand, as in the following formula 2, a tool center point vector P.sub.R in the workpiece coordinate system with the geometric error is obtained by arranging the respective geometric error as transformation matrices between transformation matrices of the respective axes in the formula 1. The formula 2 is an approximation formula that assumes a product of the geometric errors is 0 taking the geometric errors minute.

(15) $\begin{matrix} P_{R} = {.Math.}_{6} .Math. M_{5} .Math. {.Math.}_{5} .Math. M_{4} .Math. {.Math.}_{4} .Math. M_{3} .Math. {.Math.}_{3} .Math. M_{2} .Math. {.Math.}_{2} .Math. M_{1} .Math. {.Math.}_{1} .Math. P_{T} where {.Math.}_{j} = [\begin{matrix} 1 & - γ_{j} & β_{j} & δ x_{j} \\ γ_{j} & 1 & - α_{j} & δ y_{j} \\ - β_{j} & α_{j} & 1 & δ z_{j} \\ 0 & 0 & 0 & 1 \end{matrix}] (j = 1, 2, .Math., 6) & [Formula 2] \end{matrix}$

(16) Therefore, a position error ΔP.sub.W=(δx, δy, δz) at a tool center point in the workpiece coordinate system is expressed by the following formula 3.
ΔP.sub.W=P.sub.R−P.sub.I [Formula 3]

(17) Furthermore, an error ΔP.sub.O on the command value can be obtained by performing a coordinate transformation as in the following formula 4 on the position error ΔP.sub.W at the tool center point in the workpiece coordinate system.
ΔP.sub.O=M.sub.4.sup.−1.Math.M.sub.5.sup.−1.Math.ΔP.sub.W [Formula 4]

(18) Accordingly, a correction value ΔP=(Δx, Δy, Δz) of the X, Y, and Z axes is obtained in the following formula 5 that takes the command values of the respective axes and the geometric error in the above-described formula as preliminarily measured and identified parameters (geometrical parameters).
ΔP=−ΔP.sub.O=M.sub.3.Math.M.sub.2.Math.M.sub.1.Math.P.sub.T−M.sub.4.sup.−1.Math.M.sub.5.sup.−1.Math.ε.sub.6.Math.M.sub.5.Math.ε.sub.5.Math.M.sub.4.Math.ε.sub.4.Math.M.sub.3.Math.ε.sub.3.Math.M.sub.2.Math.ε.sub.2.Math.M.sub.1.Math.ε.sub.1.Math.P.sub.T [Formula 5]

(19) When a posture error is corrected, as in the following formula 6, the correction values Δa and Δc of the rotation axes, A-axis and C-axis are obtained. The correction value ΔP=(Δx, Δy, Δz) of the X, Y, and Z axes at this time can be obtained as in the following formula 7.

(20) $\begin{matrix} {\begin{matrix} Δ c = - {.Math.}_{j = 5}^{6} γ_{j} \\ Δ a = - {.Math.}_{j = 1}^{5} α_{j} \end{matrix} & [Formula 6] \\ Δ P = - Δ P_{0} = M_{3} .Math. M_{2} .Math. M_{1} .Math. P_{T} - M_{4}^{- 1} .Math. M_{5}^{- 1} .Math. {.Math.}_{6} .Math. ξ_{C} .Math. M_{5} .Math. {.Math.}_{5} .Math. ξ_{A} .Math. M_{4} .Math. {.Math.}_{4} .Math. M_{3} .Math. {.Math.}_{3} .Math. M_{2} .Math. {.Math.}_{2} .Math. M_{1} .Math. {.Math.}_{1} .Math. P_{T} where ξ_{C} = [\begin{matrix} 1 & - Δ c & 0 & 0 \\ Δ c & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{matrix}], ξ_{A} = [\begin{matrix} 1 & 0 & 0 & 0 \\ 0 & 1 & - Δ a & 0 \\ 0 & Δ a & 1 & 0 \\ 0 & 0 & 0 & 1 \end{matrix}] & [Formula 7] \end{matrix}$

(21) Thus, the correction values of the respective translational axes are added to the command values of the respective corresponding translational axes to command, thus ensuring the correction of the position error and the posture error at the tool center point due to the geometric error.

(22) On the other hand, when the deformation error due to the external force exists, the deformation error is assumed to be a function f.sub.G that indicate this error with a curved line or a line segment group. The function f.sub.G of the deformation error can be obtained as in the following formula 8 by adding an estimated value f.sub.G_sim and an identified value f.sub.G_idt that are described later. When the function f.sub.G is a discrete value, an interpolation between points is performed with a straight line and a curved line.
f.sub.G=f.sub.G_sim+f.sub.G_idt [Formula 8]

(23) This deformation error can be handled as a partial geometric error. For example, when the table 3 deforms by the load of the workpiece, its deformation error can be handled as a geometric error between the C-axis and the A-axis. A tool center point vector .sup.WP.sub.D in the workpiece coordinate system when the deformation error exists can be obtained from the following formula 9.

(24) $\begin{matrix} {}^{W}P_{D} = M_{5} .Math. ξ_{G} .Math. M_{4} .Math. M_{3} .Math. M_{2} .Math. M_{1} .Math.^{T} P = M_{D} .Math.^{T} P where ξ_{G} = [\begin{matrix} 1 & - f_{GC} & f_{GB} & f_{GX} \\ f_{GC} & 1 & - f_{GA} & f_{GY} \\ - f_{GB} & f_{GA} & 1 & f_{GZ} \\ 0 & 0 & 0 & 1 \end{matrix}] & [Formula 9] \end{matrix}$

(25) Next, FIG. 3 illustrates an exemplary numerical control device for performing a numerical control method in the disclosure. The command value generation means 12 and the servo command value transformation means 13 are identical to those in FIG. 2. The machining program 11, the servo amplifiers 14a to 14e, and the servo motors 15a to 15e are omitted for simplification of the drawing.

(26) In this numerical control device, based on the command value calculated by the command value generation means 12, gravitational deformation estimation means 18 and the axis-dependent deformation error estimation means 20 each calculate an estimated value of the deformation error (a gravitational deformation estimating step and an axis-dependent deformation error estimating step). The correction value for the deformation error is obtained by adding the estimated values of the respective deformation errors calculated in the gravitational deformation estimation means 18 and the axis-dependent deformation error estimation means 20 to input its total value to the correction value calculation means 16 (a correction value calculating step). The calculated correction value is added to the command value in an adder 21 as addition means to input it into the servo command value transformation means 13 (an adding step).

(27) A description will be given of a sequence of flows that input the mass and the barycentric position of the workpiece into the gravitational deformation estimation means 18, and calculate the estimated value of the deformation error (hereinafter referred to as a “gravitational deformation error”) based on the input mass and barycentric position of the workpiece.

(28) First, input means 17 arranged in the numerical control device inputs a mass m.sub.w and a barycentric position v.sub.w of the workpiece into the gravitational deformation estimation means 18 (a workpiece information obtaining step). The input means 17 includes a method to input the mass and the barycentric position of the workpiece that have been manually input or automatically identified.

(29) The gravitational deformation estimation means 18 calculates the gravitational deformation error in accordance with a function f.sub.G_sim(a, c, m.sub.w, v.sub.w), based on the input mass and barycentric position, and the command value generated in the command value generation means 12. The function f.sub.G_sim(a, c, m.sub.w, v.sub.w) of the deformation error is a function that takes the command values of the A-axis and the C-axis, the mass m.sub.w, and the barycentric position v.sub.w as gravitational deformation error parameters. This function is obtained in advance by, for example, an actually measured value and an analytical value.

(30) Next, a description will be given of a sequence of flows that identify the function f.sub.G_idt of the deformation error (hereinafter referred to as an “axis-dependent deformation error function”) in the axis-dependent deformation error estimation means 20 based on the geometric error identified in error identification means 19, and calculate the estimated value of the deformation error (hereinafter referred to as an “axis-dependent deformation error”).

(31) FIG. 4 is a sequence of processes that perform a measurement using the measuring device in the error identification means 19 to identify the geometric error.

(32) First, before the measurement of the geometric error, an initial and reference position of a target ball installed on a table is identified by a touch trigger probe as measurement means attached to the main spindle (S1-1).

(33) Next, an index measurement for identifying the geometric error is executed the number of times set under a measurement condition (S1-2).

(34) Next, in accordance with index conditions of a swivel axis (the A-axis) and a rotation axis (the C-axis) preliminarily set by the input means 17 and the like (S1-3), index operations of the A-axis as an inclined axis and the C-axis as the rotation axis are executed (S1-4).

(35) Then, a process (S1-5) to measure a coordinate of the target ball under the respective index conditions is repeated (S1-6). From its result, the geometric error is identified (S1-7).

(36) Next, the following describes details of Step S1-7.

(37) Under one measurement condition, it is assumed that one of the rotation axes is fixed and the other is indexed to a plurality of angles, and the center position of the target ball is measured. Difference vectors of the measurement values at the ball center position with respect to the command values under this measurement condition can be distributed to radial direction, axial direction, and tangent direction components of an index axis. These respective components can be approximated as Fourier series of a zero-order component (radius error), a first-order component (a center deviation), and a second-order component (elliptical shape), namely, an arc with error by, for example, least-square method.

(38) A radial direction component dRr.sub.i, an axial direction component dRa.sub.i, and a tangential direction component dRt.sub.i of the measurement values at a k-th index angle θ.sub.ijk of a j-th rotation axis under a measurement condition i can be expressed as the following formula 10.
dRr.sub.i=ra0.sub.i+ra1.sub.i*cos(θ.sub.ijk)+rb1.sub.i*sin(θ.sub.ijk)+ra2.sub.i*cos(2θ.sub.ijk)+rb2.sub.i*sin(2θ.sub.ijk)
dRa.sub.i=aa0.sub.i+aa1.sub.i*cos(θ.sub.ijk)+ab1.sub.i*sin(θ.sub.ijk)+aa2.sub.i*cos(2θ.sub.ijk)+ab2.sub.i*sin(2θ.sub.ijk)
dRt.sub.i=ta0.sub.i+ta1.sub.i*cos(θ.sub.ijk)+tb1.sub.i*sin(θ.sub.ijk)+ta2.sub.i*cos(2θ.sub.ijk)+tb2.sub.i*sin(2θ.sub.ijk) [Formula 10]

(39) The deformation error is handled as the geometric error between the C-axis and the A-axis. An axis-dependent deformation error when an A-axis angle is a is assumed to be f.sub.G_idt(a). When it is assumed that a fixed axis is the A-axis, the index axis is the C-axis, and the A-axis angle under the measurement condition i is an axis-dependent deformation error f.sub.G_idt(a.sub.i) when the geometric error between the C-axis and the A-axis under the measurement condition i, that is, the A-axis angle is a.sub.i is identified as in the following formula 11. Here, H is a coordinate in the Z direction of the target ball center in a table coordinate system (a coordinate system on the table that takes an intersection point of the A and C axes as an origin in an ideal state without the geometric error, and whose X-axis is parallel to the X-axis of the machine).
f.sub.G_idt X(a.sub.i)=−ra1.sub.i−f.sub.G_idt B(a.sub.i)*H
f.sub.G_idt Y(a.sub.i)=rb1.sub.i+f.sub.G_idt A(a.sub.i)*H
f.sub.G_idt Z(a.sub.i)=aa0.sub.i
f.sub.G_idt A(a.sub.i)=aa1.sub.i
f.sub.G_idt B(a.sub.i)=aa1.sub.i
f.sub.G_idt C(a.sub.i)=ta0.sub.i

(40) Under a plurality of measurement conditions, a geometric error of six degrees of freedom is identified using a function of the geometric error of six degrees of freedom taking the command value of the A-axis as a variable as in the formula 11, thus obtaining the axis-dependent deformation error f.sub.G_idt(a.sub.i) with respect to each A-axis angle a.sub.i under the plurality of measurement conditions. An axis-dependent deformation error function f.sub.G_idt(a) is obtained by interpolating the obtained axis-dependent deformation error f.sub.G_idt(a.sub.i) with respect to each A-axis angle a, with the straight line and the curved line, and by approximating the discrete value with, for example, a trigonometric function.

(41) An example of the approximation with the trigonometric function will be described. When it is assumed that an X-direction component f.sub.G_idtx(a) of the axis-dependent deformation error function can be expressed in a first-order trigonometric function, the X-direction component f.sub.G_idtx(a) of the axis-dependent deformation error function can be expressed as the following formula 12 using X-direction components δx.sub.a and δx.sub.b of an axis-dependent deformation error parameter.
f.sub.G_idt X(a)=δx.sub.a(1−cos(a))+δx.sub.b sin(a) [Formula 12]

(42) The formula 12 can be solved for the X-direction axis-dependent deformation error parameters δx.sub.a and δx.sub.b using, for example, the least-square method, by identifying an X-direction axis-dependent deformation error f.sub.G_idtx(a.sub.i) with respect to a plurality of A-axis angles a.sub.i. Five components other than the X-direction component are similarly calculated to identify the axis-dependent deformation error parameters of six degrees of freedom.

(43) Thus, the axis-dependent deformation error estimation means 20 calculates the estimated value of the axis-dependent deformation error based on the axis-dependent deformation error parameter identified in the error identification means 19 and the command value. This estimated value is combined with the estimated value of the gravitational deformation error calculated in the gravitational deformation estimation means 18 to input its total value into the correction value calculation means 16, thus obtaining the correction value to be added to the command value.

(44) Thus, according to the numerical control device and the numerical control method of the above-described five-axis machine, the estimated value of the axis-dependent deformation error caused by the deformation of the five-axis machine such that the translational axis and/or the rotation axis operates is calculated based on the axis-dependent deformation error parameter and the command value of the translational axis and/or the rotation axis by the axis-dependent deformation error estimation means 20. The mass and the barycentric position of the workpiece are input by the input means 17. The estimated value of the gravitational deformation error caused by the mass of the workpiece is calculated based on the gravitational deformation error parameter, the mass and the barycentric position received at the input means 17, and this command value by the gravitational deformation estimation means 18. The correction value of the translational axis and/or the rotation axis with respect to the error of the position and/or the posture of the tool with respect to the workpiece is calculated based on the estimated value of the axis-dependent deformation error, the estimated value of the gravitational deformation error, and this command value by the correction value calculation means 16. The correction value is added to the command value by the adder 21.

(45) Thus, the deformation error can be identified with the measuring device even when the deformation behavior changes depending on a variation for each machine table and a usage condition of the machine. For example, when a large workpiece has been loaded to make difficult to perform the measurement with an interference, the deformation error is identified in a measurable state, such as a state where only the jig has been loaded. After the deformation error is identified, information on the mass and the barycentric position of a jig and a workpiece that have been additionally loaded from that point of time is obtained. Based on the information, the deformation is estimated more accurately. Further, when the barycentric position of the workpiece is eccentric from the center of the rotation axis, C-axis, the deformation error that changes depending on the position of the C-axis can be estimated based on the mass and the barycentric position of the workpiece. The change in the deformation error depending on the position of the C-axis is also corrected by a correction based on the estimated deformation error. Further, the use of the measuring device can identify the deformation error that was difficult to be estimated.

(46) As described above, the correction of the identified and estimated error ensures the machining with higher accuracy with respect to the various workpieces.

(47) As another embodiment, a description will be given of an example that identifies the axis-dependent deformation error when the jig and the workpiece are not loaded on the table, by identifying the axis-dependent deformation error excepting the gravitational deformation error f.sub.G_sim in a state where the jig and the workpiece have been loaded on the table.

(48) The position errors of the tool center points in the X, Y, and Z-axes directions caused by the gravitational deformation error f.sub.G_sim are obtained in the following formula 13.

(49) $\begin{matrix} Δ P_{0 G_sim} = M_{3} .Math. M_{2} .Math. M_{1} .Math. P_{T} - M_{4}^{- 1} .Math. ξ_{G_sim} .Math. M_{4} .Math. M_{3} .Math. M_{2} .Math. M_{1} .Math. P_{T} where ξ_{G_sim} = [\begin{matrix} 1 & - f_{G_simC} & f_{G_simB} & f_{G_simX} \\ f_{G_simC} & 1 & - f_{G_simA} & f_{G_simY} \\ - f_{G_simB} & f_{G_simA} & 1 & f_{G_simZ} \\ 0 & 0 & 0 & 1 \end{matrix}] & [Formula 13] \end{matrix}$

(50) Accordingly, in the error identification means 19, after the position errors of the tool center points in the X, Y, and Z axes directions obtained in the formula 13 are subtracted from the measurement value on the ball center position, the calculation is performed in accordance with the formula 10, the formula 11, and the formula 12. Thus, the axis-dependent deformation error parameter excepting the gravitational deformation error f.sub.G_sim can be identified.

(51) The gravitational deformation error calculated in the gravitational deformation estimation means 18 and the axis-dependent deformation error calculated in the axis-dependent deformation error estimation means 20 are summed as the formula 8 to input its total value into the correction value calculation means 16. Thus, the correction values of the position error and the posture error of the tool center point are calculated. The calculated correction value is added to the command value to correct the position error and the posture error of the tool center point.

(52) In the above-described configuration and the other embodiment, the mass and the barycentric position of the workpiece are input. However, it may be a mass and a barycentric position including the workpiece and the jig, or may be a mass and a barycentric position of only the jig.

(53) It is explicitly stated that all features disclosed in the description and/or the claims are intended to be disclosed separately and independently from each other for the purpose of original disclosure as well as for the purpose of restricting the claimed invention independent of the composition of the features in the embodiments and/or the claims. It is explicitly stated that all value ranges or indications of groups of entities disclose every possible intermediate value or intermediate entity for the purpose of original disclosure as well as for the purpose of restricting the claimed invention, in particular as limits of value ranges.

Numerical control device and numerical control meihod for machine tool

Assignee

Inventors

Cpc classification

Classification Explorer

B23Q17/2275

PERFORMING OPERATIONS; TRANSPORTING

Classification Explorer

G05B19/404

PHYSICS

Classification Explorer

B23Q15/24

PERFORMING OPERATIONS; TRANSPORTING

Classification Explorer

B25J9/0096

PERFORMING OPERATIONS; TRANSPORTING

Classification Explorer

B23Q15/26

PERFORMING OPERATIONS; TRANSPORTING

Classification Explorer

G05B2219/41114

PHYSICS

Classification Explorer

G05B2219/49186

PHYSICS

International classification

Classification Explorer

B23Q17/22

PERFORMING OPERATIONS; TRANSPORTING

Classification Explorer

B23Q15/24

PERFORMING OPERATIONS; TRANSPORTING

Classification Explorer

B23Q15/26

PERFORMING OPERATIONS; TRANSPORTING

Classification Explorer

B25J9/00

PERFORMING OPERATIONS; TRANSPORTING

Abstract

Claims

Description