WO2017158932A1 - 機械運動軌跡測定装置 - Google Patents

機械運動軌跡測定装置 Download PDF

Info

Publication number
WO2017158932A1
WO2017158932A1 PCT/JP2016/084562 JP2016084562W WO2017158932A1 WO 2017158932 A1 WO2017158932 A1 WO 2017158932A1 JP 2016084562 W JP2016084562 W JP 2016084562W WO 2017158932 A1 WO2017158932 A1 WO 2017158932A1
Authority
WO
WIPO (PCT)
Prior art keywords
motion trajectory
sensor signal
trajectory
unit
acceleration sensor
Prior art date
Application number
PCT/JP2016/084562
Other languages
English (en)
French (fr)
Inventor
智哉 藤田
正行 植松
Original Assignee
三菱電機株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 三菱電機株式会社 filed Critical 三菱電機株式会社
Priority to DE112016006602.3T priority Critical patent/DE112016006602T5/de
Priority to JP2017526604A priority patent/JP6199003B1/ja
Priority to CN201680024212.0A priority patent/CN107533325B/zh
Priority to US15/564,264 priority patent/US10543574B2/en
Publication of WO2017158932A1 publication Critical patent/WO2017158932A1/ja

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical 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/404Numerical 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
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B23MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23QDETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
    • B23Q17/00Arrangements for observing, indicating or measuring on machine tools
    • B23Q17/12Arrangements for observing, indicating or measuring on machine tools for indicating or measuring vibration
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B23MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23QDETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
    • B23Q17/00Arrangements for observing, indicating or measuring on machine tools
    • B23Q17/22Arrangements for observing, indicating or measuring on machine tools for indicating or measuring existing or desired position of tool or work
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical 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/19Numerical 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 positioning or contouring control systems, e.g. to control position from one programmed point to another or to control movement along a programmed continuous path
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B23MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23QDETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
    • B23Q2717/00Arrangements for indicating or measuring
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/37Measurements
    • G05B2219/37388Acceleration or deceleration, inertial measurement
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/37Measurements
    • G05B2219/37392Motion
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/37Measurements
    • G05B2219/37493Use of different frequency band pass filters to separate different signals

Definitions

  • the present invention relates to a machine motion trajectory measuring apparatus for measuring a motion trajectory of a mechanical device such as a numerically controlled machine tool, an industrial machine, a robot, or a conveyor.
  • the servo control device is a device that performs control using an actuator so that the position of a driven body that is a drive target detected using a position detector matches a command position.
  • a mechanical device having a multi-degree of freedom such as a two-dimensional plane or a three-dimensional space such as a numerically controlled machine tool, an industrial machine, a robot, or a conveyor has a servo control device called an axis of one degree of freedom. Have multiple.
  • an actuator attached to each axis drives and controls a driven body of each axis.
  • the mechanical device combines these axial motions to achieve multi-degree-of-freedom mechanical motion.
  • trajectory control Servo control performed so that the motion trajectory accurately follows the command trajectory that is the commanded path. If a response error occurs in the servo controller of each axis due to disturbance factors such as friction generated during trajectory control or mechanical structure vibration generated during trajectory control, the motion trajectory of the driven body deviates from the command trajectory. A trajectory error occurs.
  • Patent Documents 1, 2, and 3 disclose a method for measuring and detecting a machine motion trajectory.
  • this displacement is read when an arc motion is performed to keep the relative distance between two steel balls coupled via a displacement detector constant.
  • This measuring method is called a ball bar method and is widely used.
  • the general ball bar method two steel balls are mechanically connected, and the displacement of the moving bar is measured with a displacement sensor, thereby measuring the trajectory error during arc movement with an accuracy of several micrometers. it can.
  • Patent Document 2 detects either the speed or acceleration of a control target, and compares any of the detected speed or acceleration with a planned speed or acceleration, thereby detecting a trajectory error. Detect outbreaks.
  • Patent Document 2 discloses a method for suppressing a trajectory error by feeding back a speed or acceleration component of the detected trajectory error.
  • Patent Document 3 discloses a method for estimating the motion trajectory of a driven body by second-order integration of the acceleration of the driven body when performing sinusoidal motion using one or more axes. Further, Patent Document 3 discloses a method for estimating a motion trajectory with high accuracy by changing a sensitivity coefficient so that an error that occurs during integration does not exceed a threshold value.
  • the motion accuracy test method disclosed in Patent Document 1 has a problem that a measurable locus is limited to an arc.
  • the position of the rotation center is calculated in advance to set up the apparatus, and the installation position of one sphere is matched with the rotation center at the time of trajectory measurement. There is a problem that it takes time and effort to set up.
  • the method for suppressing the trajectory error disclosed in Patent Document 2 there is no limit on the measurable trajectory.
  • the locus error suppression method disclosed in Patent Document 2 since the occurrence of an error is detected using acceleration or velocity, there is a problem that the amount of locus error from the command locus of the driven body cannot be obtained. Further, in the locus error suppression method disclosed in Patent Document 2, the locus error can be suppressed by feeding back acceleration or velocity.
  • the trajectory error suppression method disclosed in Patent Document 2 can suppress trajectory errors by feeding back acceleration or velocity, but has a problem that it cannot be used for parameter adjustment of other trajectory error suppression methods.
  • Another trajectory error suppression method is a correction method that suppresses trajectory errors derived from friction, which is called lost motion using a friction model or stick motion.
  • the motion trajectory estimation method disclosed in Patent Document 3 has a problem that measurable trajectories are limited to arc trajectories, elliptical trajectories, and spherical trajectories that can be realized by combining sine waves or sine waves.
  • the accuracy of the relative motion trajectory between two points on the movable side and the fixed side may be important. More specifically, in a machine in which a tool is mounted on the fixed side of a casting, which is a mechanical structure, and the servo controller for each axis and the workpiece are mounted on the casting, the workpiece is contour controlled. Then, the process of removing the material from the workpiece is performed by causing the tool and the workpiece to interfere with each other. In such a case, the driving reaction force generated by the reaction of the driving force of the shaft that is the movable part is transmitted to the fixed side via the mechanical structure, and deformation or vibration may occur on the fixed side. Therefore, it is necessary to measure the relative motion trajectory between the movable side and the fixed side.
  • the methods disclosed in Patent Documents 1 to 3 have a problem that the measurement accuracy is low because the relative locus error on the fixed side cannot be taken into consideration.
  • the present invention has been made in view of the above, and an object thereof is to obtain a mechanical motion trajectory measuring apparatus capable of measuring a trajectory error in an arbitrary motion trajectory with high accuracy with a simple setup.
  • the mechanical motion trajectory measuring apparatus of the present invention has an actuator, and the position of the actuator or the driven state is such that the motion trajectory of the driven body follows the command trajectory.
  • a machine motion trajectory measuring device that feeds back a detected position signal output from a position detector that detects the position of the body and measures the motion trajectory of the mechanical device that drives the actuator, and measures the acceleration of the motion trajectory measurement object
  • an acceleration sensor that outputs the acceleration sensor signal
  • a sensor signal separation unit that separates the acceleration sensor signal into two or more frequency bands.
  • the machine movement trajectory measuring device separates the detected position signal into the same frequency band as the sensor signal separating unit, the acceleration sensor signal separated by the sensor signal separating unit, and the detected position signal separating unit.
  • a data calibration unit that calibrates an acceleration sensor signal in each of two or more frequency bands using the detected detection position signal to obtain a motion trajectory component of each of the two or more frequency bands, and two or more frequencies
  • a motion trajectory calculation unit that combines the motion trajectory components of the bands and outputs the motion trajectory components is provided.
  • FIG. 1 is a perspective view of a machine movement trajectory measuring apparatus and a numerically controlled machine tool according to a first embodiment.
  • 1 is a diagram showing a machine motion trajectory measuring apparatus according to a first embodiment, a numerically controlled machine tool to which the machine motion trajectory measuring apparatus according to the first embodiment is applied, and a servo control device. Configuration diagram of the motor drive section shown in FIG.
  • a plot of the operation trajectory at the detector position, the motion trajectory to be controlled, and the command trajectory during arc motion on the xy plane The figure which shows the acceleration waveform of the x-axis direction when noise does not exist, and the displacement amount of the x-axis direction calculated from the acceleration.
  • trajectory error in xy plane on the conditions similar to FIG. The figure which shows the acceleration waveform of the x-axis direction in case noise exists, and the displacement amount of the x-axis direction calculated from the acceleration
  • Configuration diagram of the data calibration unit shown in FIG. Configuration diagram of the motion trajectory calculation unit shown in FIG.
  • FIG. 4 The figure which shows the machine motion locus measuring apparatus which concerns on Embodiment 4, the numerical control machine tool to which the machine motion locus measuring apparatus which concerns on Embodiment 4 is applied, and a servo control apparatus.
  • the block diagram of the data calibration part with which the mechanical-motion-trajectory measuring apparatus which concerns on Embodiment 5 is provided
  • FIG. 1 is a perspective view of a machine motion trajectory measuring apparatus and a numerically controlled machine tool according to the first embodiment.
  • FIG. 2 is a diagram showing a detailed configuration of the x-axis drive mechanism shown in FIG. 1 and a configuration of a servo control device that controls the operation of the x-axis drive mechanism.
  • a numerically controlled machine tool 99 shown in FIG. 1 is an example of an object whose machine motion trajectory measuring apparatus measures a mechanical motion trajectory.
  • the numerically controlled machine tool 99 is a vertical three-axis vertical machine tool, and has a machine structure called a C column structure.
  • the numerically controlled machine tool 99 includes a gantry 21, a saddle 24 installed on the gantry 21 and driven in the y-axis direction, a work table 4 installed on the saddle 24, and fixed to the gantry 21 and above the gantry 21. And a column 5 extending.
  • a ram 6 is attached to the column 5, and a workpiece 17, which is an object for measuring a movement locus, is installed on the work table 4.
  • the numerically controlled machine tool 99 is an x-axis drive mechanism 15x that is an actuator that is attached to the saddle 24 and drives the work table 4 in the x-axis direction, and an actuator that is attached to the mount 21 and drives the saddle 24 in the y-axis direction.
  • a y-axis drive mechanism 15y and a z-axis drive mechanism 15z, which is an actuator attached to the column 5 and driving the ram 6 in the z-axis direction, are provided.
  • the x-axis drive mechanism 15x includes an x-axis motor 1x, a feed screw 2x that is a feed shaft driven by the x-axis motor 1x, and a rotation angle detector 3x that detects a rotation angle of the feed screw 2x.
  • the y-axis drive mechanism 15y includes a y-axis motor 1y, a feed screw 2y that is a feed shaft driven by the y-axis motor 1y, and a rotation angle detector 3y that detects a rotation angle of the feed screw 2y.
  • the z-axis drive mechanism 15z includes a z-axis motor 1z, a feed screw 2z driven by the z-axis motor 1z, and a rotation angle detector 3z that detects a rotation angle of the feed screw 2z.
  • the work table 4 is driven by the x-axis drive mechanism 15x, and the saddle 24 and the x-axis drive mechanism 15x installed thereon are driven by the y-axis drive mechanism 15y.
  • the ram 6, the main shaft 7 and the tool 16 are driven by the z-axis drive mechanism 15 z attached to the column 5, and the workpiece 17 is machined by the tool 16.
  • FIG. 2 shows a detailed configuration of the x-axis drive mechanism 15x shown in FIG. 1, and shows a servo control device 101 for controlling the position of the driven body in the x-axis direction.
  • FIG. 2 shows only the configuration of the x-axis drive mechanism 15x among the three drive mechanisms shown in FIG.
  • the configurations of the y-axis drive mechanism 15y and the z-axis drive mechanism 15z are the same as the configurations of the x-axis drive mechanism 15x.
  • the x-axis drive mechanism 15x, the y-axis drive mechanism 15y, and the z-axis drive mechanism 15z differ in the following points.
  • the driven body that is the control target of the x-axis drive mechanism 15 x is the work table 4, whereas the driven body that is the control target of the y-axis drive mechanism 15 y is located above the saddle 24 and the saddle 24. It differs from the attached x-axis drive mechanism 15x in that the driven body that is the control target of the z-axis drive mechanism 15z is the column 5 and the main shaft 7 attached to the column 5.
  • the x-axis drive mechanism 15x shown in FIG. 2 includes a saddle 24, and two support bearings 10x each outer ring fixed to the saddle 24 and each inner ring rotatably supporting the feed screw 2x.
  • the x-axis drive mechanism 15x fixes the nut 9x meshed with the feed screw 2x, the work table 4 that moves in the axial direction of the feed screw 2x by the nut 9x, and one end of the feed screw 2x when the x-axis motor 1x rotates. Coupling 8x.
  • Rotational motion of the x-axis motor 1x is transmitted to the feed screw 2x via the coupling 8x, and the rotational motion of the feed screw 2x is converted into linear motion by the nut 9x.
  • the work table 4 fixed to the nut 9x moves along the axial direction of the feed screw 2x.
  • the feed screw 2x is supported by the support bearing 10x, the x-axis motor 1x and the support bearing 10x are supported by the saddle 24, and the rotation angle detector 3x is supported by the x-axis motor 1x.
  • the servo control device 101 for controlling the position of the driven body in the x-axis direction includes a position command generation section 11x that generates a position command for controlling the position of the driven body in the x-axis direction, and a position command generation section 11x. And a motor drive unit 12x that controls the rotation angle of the x-axis motor 1x according to the generated position command.
  • FIG. 2 shows only the servo control device 101 for controlling the position of the driven body in the x-axis direction, but the servo control device for controlling the position of the driven body in each of the y-axis direction and the z-axis direction. Is similarly configured.
  • the position command generated by the position command generation unit 11x is transmitted to the motor drive unit 12x, and the motor drive unit 12x that has received the position command sets the screw of the feed screw 2x to the motor rotation angle detected by the rotation angle detector 3x.
  • the rotation angle of the x-axis motor 1x is controlled so that the error between the detection position obtained by multiplying the pitch and the position command becomes small.
  • a linear motor may be used as the movable shaft of the numerically controlled machine tool 99 instead of the x-axis motor 1x and the feed screw 2x.
  • a linear encoder or a laser displacement meter that can directly detect the position of the work table 4 may be used for the movable shaft of the numerically controlled machine tool 99 instead of the rotation angle detector 3x.
  • the relative displacement between the tool 16 and the workpiece 17 is important. If a relative displacement occurs between the tool 16 and the workpiece 17 during machining while performing contour control, a machining error occurs if material remains uncut or overcut in the workpiece 17. It is. In the servo control device 101, feedback control is performed in the motor drive unit 12x in order to prevent such processing errors from occurring.
  • Dynamic disturbances include elastic deformation caused by the coupling 8x, the feed screw 2x or the support bearing 10x, vibration caused by the feed screw 2x or the support bearing 10x, and a change in the posture of the column 5 or the ram 6 shown in FIG.
  • vibrations of the column 5 or ram 6 and errors due to the frictional force of the movable shaft include the mass of the workpiece 17, the aging of the machine, the wear of the feed screw 2x or the support bearing 10x, the amount of lubricating oil on each movable shaft, the temperature change, and other factors in the factory. It is known to change depending on the operating status of the machine.
  • the position command generation unit 11x and the motor drive unit 12x use a friction and vibration model to maintain an accurate amount of error in order to maintain highly accurate contour control performance.
  • a model-based error correction method in which a correction amount for canceling the error is generated and feedforward control is performed is used.
  • a triaxial acceleration sensor 13 which is an example of an object acceleration sensor for measuring the acceleration of the workpiece 17, is attached near the workpiece 17 that is an object for measuring the movement locus.
  • the triaxial acceleration sensor 13 measures accelerations in three orthogonal directions.
  • the triaxial acceleration sensor 13 and the machine motion trajectory measuring apparatus 100 are connected to each other via a cable 40.
  • Acceleration sensor signals indicating the respective accelerations in the three-axis directions detected by the three-axis acceleration sensor 13 are input to the machine motion trajectory measuring apparatus 100 via the cable 40.
  • the respective acceleration components in the three-axis direction detected by the three-axis acceleration sensor 13 are referred to as an x-axis direction acceleration component, a y-axis direction acceleration component, and a z-axis direction acceleration component.
  • the three-axis acceleration sensor 13 can measure the acceleration in the three-axis directions orthogonal to each other by using one sensor, the acceleration in three dimensions can be measured by using one three-axis acceleration sensor 13.
  • three uniaxial acceleration sensors that measure the respective accelerations in the three axial directions may be used instead of the triaxial acceleration sensor 13.
  • the triaxial acceleration sensor 13 there are methods for attaching the triaxial acceleration sensor 13 such as attachment using a magnetic force by a magnet, fastening using a jig and a screw, fixation using wax, or fixation using an adhesive. Since measurement is possible simply by fixing to the work table 4 which is a measurement object, it can be attached by a simple operation immediately before measurement. In addition, since no special adjustment work is required when installing the acceleration sensor, even a beginner can easily install it.
  • the number of cableless acceleration sensors using wireless has increased, and when a wireless acceleration sensor is used, it is not necessary to consider the handling of the cable 40, so that it can be installed more easily.
  • the triaxial acceleration sensor 13 is always attached to the back side or inside the work table 4. Also good.
  • the machine motion trajectory measuring apparatus 100 may be mounted outside the numerically controlled machine tool 99, or may be mounted inside the numerically controlled machine tool 99.
  • FIG. 3 is a diagram showing a machine motion trajectory measuring apparatus according to the first embodiment, a numerically controlled machine tool to which the machine motion trajectory measuring apparatus according to the first embodiment is applied, and a servo control apparatus.
  • FIG. 4 is a block diagram of the motor drive unit shown in FIG.
  • FIG. 3 differs from the servo control device 101 shown in FIG. 2 in that the servo control device 101 shown in FIG. 3 generates a position command in addition to the position command generation unit 11x and the motor drive unit 12x. A portion 11y and a motor drive portion 12y.
  • the motor drive unit 12x is inputted position command Pix generated by the position command generating unit 11x, a motor rotation angle detected by the rotation angle detector 3x feedback position Pd x is inputted.
  • the motor drive unit 12x based on the position command Pix and feedback position Pd x, and outputs a torque command Icx a drive command of the x-axis motor 1x.
  • the motor drive unit 12y the position command Piy generated by the position command generating unit 11y is input, a motor rotation angle detected by the rotation angle detector 3y feedback position Pd y is input.
  • the motor drive unit 12y based on the position command Piy and feedback position Pd y, and outputs a torque command Icy a drive command y-axis motor 1y.
  • FIG. 4 shows the configuration of the motor drive unit 12x shown in FIG.
  • the motor drive unit 12x includes a subtraction unit 124a for determining the position deviation which is a response error between a response position as command Pix feedback position Pd x, velocity command executes the control with respect to the position deviation subtraction unit 124a is determined And a position control unit 123 for generating.
  • the motor drive unit 12x performs differential operation of the feedback position Pd x, a differentiating unit 122 that calculates the velocity feedback value, obtained by the speed command and the differentiating unit 122 obtained by the position control unit 123
  • An adder / subtractor 124b that calculates a speed deviation from the feedback value and a speed controller 121 that outputs a torque command Icx that is a drive command are provided.
  • the addition / subtraction unit 124 a obtains a position deviation that is a deviation between the position command Pix and the feedback position Pd x and outputs the position deviation to the position control unit 123.
  • the position control unit 123 performs position control processing such as proportional control, PI (Proportional Integral Differential) control, or PID (Proportional Integral Differential) control so as to reduce the positional deviation output from the adder / subtractor 124a, thereby reducing the positional deviation.
  • the speed command to be output is output.
  • an actual speed obtained by differentiating the feedback position Pd x is obtained.
  • the adder / subtractor 124 b obtains a speed deviation that is a deviation between the speed command obtained by the position controller 123 and the actual speed of the feedback position Pd x obtained by the differential calculator 122, and outputs it to the speed controller 121.
  • the speed control unit 121 performs a PI control speed control process so as to reduce the speed deviation output from the addition / subtraction unit 124b, and outputs a torque command Icx. 3 is configured similarly to the motor drive unit 12x shown in FIG.
  • Control that uses the rotation angle detectors 3x and 3y to obtain the feedback positions Pd x and Pd y is called semi-closed loop control.
  • the rotation angle detector 3x is positioned between the mounting position of the rotation angle detector 3x and the work table 4 that is the true control object, or the rotation angle detector 3x.
  • the rotation angle detector 3x Between the mounting position and the workpiece 17 shown in FIG. 1, there are mechanical structures such as a feed screw 2x and a nut 9x. Therefore, a mechanical transmission delay occurs between the rotation angle detector 3x and the work table 4, or between the rotation angle detector 3x and the workpiece 17. Therefore, the feedback position Pd x detected by the rotation angle detector 3x does not coincide with the movement locus of the work table 4 or the workpiece 17.
  • a sudden disturbance force acts on the movable shaft by reversing the direction of the applied friction. Therefore, it is known that in the feedback control, the generated position deviation cannot be completely suppressed and an error occurs.
  • a typical example is a trajectory error called lost motion or stick motion during arc motion.
  • a sine wave command having a phase difference of 90 degrees is commanded between the x-axis motor 1x and the y-axis motor 1y.
  • FIG. 5 is a diagram in which an operation trajectory at the detector position, a motion trajectory to be controlled, and a command trajectory during arc motion are plotted on an xy plane.
  • the detector position represents the feedback position Pd x and the feedback position Pd y described above.
  • the operation trajectory at the detector position is indicated by a dotted line
  • the motion trajectory to be controlled is indicated by a solid line
  • the command trajectory during the arc motion is indicated by a one-dot chain line.
  • the trajectory error is enlarged and displayed 400 times in the radial direction.
  • the command trajectory is a circle, whereas the motion trajectory at the detector position and the motion trajectory of the workpiece 17 to be controlled indicate a protruding error pattern at the quadrant switching position.
  • stick motion occurs because feedback control is delayed due to reversal of the direction in which friction is applied.
  • the motion trajectory of the control target shows a trajectory that jumps outward at the quadrant switching position and then bites inward. This is a phenomenon called lost motion that occurs due to a delay in the response of the controlled object because mechanical elements such as the feed screw 2x and the nut 9x are present in addition to the stick motion caused by friction.
  • the triaxial acceleration sensor 13 is located near the workpiece 17 on the work table 4. is set up.
  • the machine motion trajectory measuring apparatus 100 includes an acceleration component a x in the x-axis direction measured by the three-axis acceleration sensor 13, an acceleration component a y in the y-axis direction measured by the three-axis acceleration sensor 13, and an x-axis motor.
  • a feedback position Pd x that is a 1x detector position and a feedback position Pd y that is a detector position of the y-axis motor 1y are input.
  • the acceleration component in the z-axis direction can be ignored when adjusting the parameter for friction correction.
  • the acceleration component a x and the acceleration component a y may be simply referred to as an acceleration sensor signal.
  • the feedback positions Pd x and Pd y may be referred to as detector positions.
  • a machine motion trajectory measuring apparatus 100 shown in FIG. 3 includes a sensor signal separation unit 30 that separates each of an acceleration component a x in the x- axis direction and an acceleration component a y in the y-axis direction into two or more frequency bands, and a feedback.
  • a motor signal separation unit 31 that is a detection position signal separation unit that separates each of the positions Pd x and Pd y into two or more frequency bands that are the same as the sensor signal separation unit 30 is provided.
  • the machine motion trajectory measuring apparatus 100 uses the feedback positions Pd x and Pd y as reference signals for the separated frequency bands of the sensor signal separation unit 30 and the motor signal separation unit 31 to obtain acceleration components a x and a y .
  • the data calibration unit 32 that calibrates and calculates and outputs a position response in each frequency band, and adds the position response for each frequency band calculated by the data calibration unit 32, and operates on the x-axis and the y-axis.
  • a motion trajectory calculation unit 33 that calculates a motion trajectory of the workpiece 17 and outputs the motion trajectory as information indicating the motion trajectory.
  • Information indicating the driving trajectory calculated by the motion trajectory calculating unit 33 is output to a trajectory display device 35 connected to the mechanical motion trajectory measuring device 100.
  • the trajectory display device 35 is a display device represented by a video monitor, and the driving trajectory calculated by the trajectory display device 35 is displayed on the screen of the display device. As a result, the driving trajectory is presented to the experimenter.
  • the trajectory display device 35 inputs the motion trajectory and the command trajectory at the feedback positions Pd x and Pd y , and displays the measurement result of the command trajectory and the motion trajectory at a plurality of positions as shown in FIG. Also good.
  • the velocity Vx can be calculated from the numerical sequence shown in the following equation (1) using the acceleration component a x and the sampling time dt.
  • the displacement amount Px is calculated from the numerical sequence shown in the following equation (2) using the velocity Vx and the sampling time dt.
  • the initial values of speed Vx and Px are set to zero.
  • the displacement amount Py can be calculated from the acceleration with respect to the y-axis, the motion locus can be calculated from the displacement amounts Px and Py.
  • FIG. 6 is a diagram showing an acceleration waveform in the x-axis direction when there is no noise, and a displacement amount in the x-axis direction calculated from the acceleration.
  • the horizontal axis represents time
  • the vertical axis represents acceleration.
  • the horizontal axis represents time
  • the vertical axis represents the amount of displacement.
  • the respective waveforms of acceleration and displacement shown in FIG. 6 are obtained when the acceleration sampling time dt is 1 ms.
  • FIG. 7 is a diagram showing the result of calculating the trajectory error in the xy plane under the same conditions as in FIG.
  • the actual trajectory that is an actual driving trajectory is indicated by a dotted line
  • the calculation result of the trajectory error calculated using the above formulas (1) and (2) is indicated by a solid line.
  • the dotted real trajectory overlaps with the solid line calculation result.
  • the motion trajectory can be calculated from the acceleration with a simple calculation.
  • Noise refers to all signal components other than the acceleration component of the motion trajectory, and is not strictly distinguished. Further, since it is difficult to actually observe the acceleration component of the true motion trajectory from which noise has been completely removed, an ideal acceleration sensor signal is generated by simulation in FIG.
  • FIG. 8 is a diagram showing an acceleration waveform in the x-axis direction when noise is present, and a displacement amount in the x-axis direction calculated from the acceleration.
  • the horizontal axis represents time
  • the vertical axis represents acceleration.
  • the horizontal axis represents time
  • the vertical axis represents the amount of displacement.
  • FIG. 9 is a diagram showing the result of calculating the trajectory error in the xy plane under the same conditions as in FIG.
  • the actual trajectory that is the actual driving trajectory is indicated by a dotted line
  • the calculation result of the trajectory error is indicated by a solid line.
  • FIG. 9 shows that the trajectory error calculated by the calculation does not match the actual trajectory.
  • a calculation error resulting from such an integration operation is called an integration error.
  • the sensor signal separation unit 30 and the motor signal separation unit 31 separate frequency bands, and the data calibration unit 32 generates an integration error by using an acceleration sensor signal for each frequency band. Reduce.
  • the configuration of the sensor signal separation unit 30, the motor signal separation unit 31, and the data calibration unit 32 will be specifically described below.
  • FIG. 10 is a block diagram of the sensor signal separation unit shown in FIG.
  • the sensor signal separation unit 30 separates the acceleration sensor signal into a noise band in order to remove noise components included in the acceleration components a x and a y that are acceleration sensor signals.
  • the noise band is a noise component of the acceleration sensor signal. By removing the noise component, the error during acceleration compensation can be reduced.
  • the acceleration components a x and a y input to the sensor signal separation unit 30 are separated from the noise components a xn and a yn by the first noise removal unit 301.
  • signal components a xs and a ys other than the noise components a xn and a yn are extracted.
  • First noise removing unit 301 as a low-pass filter described by the transfer function G filt represented by the following formula (3), is mounted on the sensor signal separating unit 30.
  • Tfilt in the following formula (3) is a cutoff time constant of the low-pass filter.
  • Each of the signal components a xs and a ys output from the first noise removal unit 301 is converted into a low frequency band component a xl and a yl and a high frequency by the first signal extraction unit 302 and the second signal extraction unit 303.
  • the band components are separated into a xh and a yh .
  • the first signal extraction unit 302 is a low-pass filter described by a transfer function Gl expressed by the following equation (4).
  • T cutoff in the following equation (4) is a cutoff time constant. However, in order to effectively remove the noise component, it is preferable that Tfilt ⁇ Tcutoff .
  • the second signal extraction unit 303 is a high-pass filter described by a transfer function Gh expressed by the following equation (5).
  • each filter is designed so that the sum of the transfer function of the first signal extraction unit 302 and the transfer function of the second signal extraction unit 303 is 1, the signals of the signal components a xs and a ys Can be extracted without excess or deficiency. That is, in order to prevent the signal from becoming excessive or insufficient before and after the frequency separation, it is necessary to design each filter so as to satisfy the following relational expression (6). That is, the sensor signal separation unit 30 includes a high frequency band filter designed so that the sum of transfer functions with the low frequency band filter is 1.
  • the high frequency band is a vibration frequency component of the mechanical device that drives the actuator.
  • FIG. 11 is a block diagram of the motor signal separation unit shown in FIG. Motor signal separation section 31, feedback position Pd x, each of the noise component Pd xn of Pd y, the Pd yn removed, the noise component Pd xn, signal components Pd xs other than Pd yn, and outputs the separated Pd ys a second noise removing unit 311, the signal component Pd xs, low frequency band component Pd xl of Pd ys, the first signal extraction unit 312 for extracting a Pd yl, signal components Pd xs, Pd ys of the high frequency band components Pd and a second signal extraction unit 313 for extracting xh and Pd yh .
  • the transfer function Gl is equal to the transfer function Gh expressed by the above equation (5).
  • FIG. 12 is a block diagram of the data calibration unit shown in FIG.
  • the data calibration unit 32 includes the low frequency band components a xl and a yl separated by the sensor signal separation unit 30 shown in FIG. 10 and the low frequency band components Pd xl and Pd separated by the motor signal separation unit 31 shown in FIG.
  • the first data calibration unit 321 for calculating the low frequency band components PT xl and PT yl of the motion trajectory using yl and the high frequency band components a xh , a separated by the sensor signal separation unit 30 shown in FIG. yh and high frequency band components Pd xh separated in the motor signal separation unit 31 shown in FIG. 11, by using the Pd yh, frequency band components PT xh movement trajectory, and the second data calibration unit 322 which calculates the PT yh Is provided.
  • the displacement amounts Pxl and Pyl calculated by numerically integrating the low frequency band components a xl and a yl of acceleration, and the low frequency band Since the difference between the components Pd xl and Pd yl is an integration error, the integration error can be reduced by compensating for the difference.
  • a compensation method of the integral error there is a method of calculating the low frequency band components PT xl and PT yl by approximating the integral error with a polynomial and subtracting the approximated integral error from the displacement amounts Pxl and Pyl.
  • the unknown parameter of the approximate expression may be determined by the least square method, or a numerical solution method such as the downhill simplex method may be used.
  • an integration error compensation method in the second data calibration unit 322 an integration error compensation method in a low frequency band may be used.
  • the low frequency band is a control band of the mechanical device that drives the actuator.
  • the second data calibration unit 322 uses the high frequency band components Pd xh and Pd yh calculated by numerical integration as low frequency band components. May be passed through a high-pass filter that excludes and outputs the result as high-frequency band components PT xh and PT yh .
  • the low frequency band components PT xl and PT yl and the high frequency band components PT xh and PT yh calculated by the data calibration unit 32 are input to the motion trajectory calculation unit 33 shown in FIG.
  • FIG. 13 is a block diagram of the motion trajectory calculation unit shown in FIG.
  • the motion trajectory calculation unit 33 includes a trajectory coupling unit 331 and a trajectory error calculation unit 332.
  • the trajectory coupling unit 331 couples the low frequency band components PT xl and PT yl and the high frequency band components PT xh and PT yh of the motion trajectory calculated by the data calibration unit 32 shown in FIG. That is, the trajectory coupling unit 331 calculates the x-axis direction component PT x of the motion trajectory at the mounting position of the triaxial acceleration sensor 13 based on the following equation (7), and the triaxial acceleration sensor 13 of the triaxial acceleration sensor 13 based on the following equation (8). calculating the y-axis direction component PT y motion trajectory in the mounting position.
  • Trajectory error calculation unit 332 based on the calculated x-axis direction component PT x and y-axis direction component PT y a locus coupling portion 331, x-axis direction of the trajectory error PT 'x and y-axis direction of the trajectory error PT' Calculate y and output.
  • a method for displaying the trajectory error of the arc motion a method of displaying the trajectory error in the radial direction as shown in the following formulas (9) and (10) is known.
  • Rcom shown in the following formulas (9) and (10) represents a command radius
  • MAG represents an enlarged display magnification of a trajectory error.
  • FIG. 14 is a diagram comparing the actual motion trajectory, the calculation result of the trajectory error calculated by the mechanical motion trajectory measuring apparatus according to Embodiment 1, and the operation trajectory at the detector position.
  • the actual motion trajectory is indicated by a one-dot chain line
  • the calculation result of the trajectory error is indicated by a solid line
  • the operation trajectory at the detector position is indicated by a dotted line.
  • FIG. 14 by using the machine motion trajectory measuring apparatus 100 according to the first embodiment, it is possible to calculate the motion trajectory of the workpiece 17 with an accuracy of 1 micrometer.
  • the acceleration sensor signal and the detected position signal are separated into two or more bands and calibrated, so that an arbitrary setup is possible with a simple setup.
  • the trajectory error in the motion trajectory can be measured with high accuracy.
  • the data calibration unit corrects the calculation error of the motion trajectory component for each frequency band from the acceleration sensor signal separated by the sensor signal separation unit, and the frequency band separated by the detection position signal separation unit. Each detection position signal is used.
  • the calculation error can be effectively corrected by using the detected position signal as a reference signal.
  • the data calibration unit is configured to perform second-order integration of the acceleration input signal at the time of data calibration and correct the integration error so that the average value of the integration results becomes zero.
  • the reference position does not move, it is only necessary to correct the integration error so that the average value of the integration results becomes 0, and the integration error can be easily compensated.
  • the data calibration unit according to Embodiment 1 has a low frequency that cannot be passed by the high-pass filter used in the sensor signal separation unit with respect to the signal output as an integration result when performing data calibration in the high frequency band. It is configured to cut the components and output them. It is possible to remove the integration error component by cutting and outputting the low frequency component caused by the integration error.
  • the difference between the second-order integration result of the acceleration sensor signal separated by the sensor signal separation unit and the low frequency band component of the detection position signal when the data calibration of the low frequency band is performed.
  • An error is corrected so as not to exceed the allowable value, and the motion trajectory of the object in the low frequency band is calculated.
  • the machine motion trajectory measuring apparatus 100 includes a motion trajectory display unit that displays a motion trajectory of the mechanical device, and the motion trajectory display unit obtains a target trajectory from a target position input to the position command generation unit.
  • the command locus is synthesized from the command position output from the position command generation unit
  • the detection portion trajectory is synthesized from the detection position signal
  • one of the target locus, the command locus, and the detection portion locus is The motion trajectory is displayed on the motion trajectory display unit in an overlapping manner.
  • FIG. FIG. 15 is a diagram showing a machine motion trajectory measuring apparatus according to the second embodiment, a numerically controlled machine tool to which the machine motion trajectory measuring apparatus according to the second embodiment is applied, and a servo control apparatus. Differences between the first embodiment and the second embodiment are as follows. In the following description, the same parts as those in the first embodiment are denoted by the same reference numerals, and the description thereof is omitted. Only different parts will be described here. (1) The machine motion trajectory measuring apparatus 100 according to the second embodiment includes a machine response simulation unit 34 in addition to the sensor signal separation unit 30, the motor signal separation unit 31, the data calibration unit 32, and the motion trajectory calculation unit 33. thing.
  • the mechanical motion trajectory measuring apparatus 100 instead of the feedback position Pd y of the feedback position Pd x and y-axis motor 1y of the x-axis motor 1x, generated by the position command generating unit 11x.
  • the position command Pix and the position command Py generated by the position command generation unit 11y are input.
  • the data calibration unit 32 of the machine motion trajectory measuring apparatus 100 according to the second embodiment calibrates the acceleration components a x and a y using the feedback position Pm calculated by the machine response simulation unit 34 as a reference signal, Calculating the position response in the frequency band of.
  • FIG. 16 is a block diagram of the machine response simulation unit shown in FIG.
  • the machine response simulation unit 34 includes a position control unit 123A, a speed control unit 121, a differential operation unit 122, and a machine model 341 configured in the same manner as the motor drive unit 12x illustrated in FIG.
  • Examples of the machine model 341 include a 2-inertia model and a 3-inertia model.
  • the 2-inertia model is a model obtained by approximating the reciprocal 1 / J of the load inertia J of the inertia model movable shaft or the inertia of the motor and the driven body by a 2-inertia vibration system.
  • the three-inertia model is a model in which the motor inertia, the driven body, and the inertia of the feed screw are approximated by a three-inertia vibration system.
  • the machine response simulation unit 34 that generates the feedback position Pm using the position command Pix, but the machine response simulation unit 34 that generates the feedback position Pm using the position command Piy has the same configuration.
  • the position control unit 123 is configured similarly to the motor drive unit 12y shown in FIG.
  • the rotation angle detectors 3x and 3y have a low resolution, and thus the rotation angle detector. Even when the feedback positions output from 3x and 3y cannot be used as reference signals, the feedback positions calculated by the machine model 341 can be used as reference signals.
  • the mechanical motion trajectory measuring apparatus 100 according to the second embodiment can virtually calculate the driven body position in the semi-closed loop control motor driving units 12x and 12y, and can use the calculated driven body position as a reference signal. .
  • Embodiment 3 The difference between the first embodiment and the third embodiment is that, in the mechanical motion trajectory measuring apparatus 100 according to the third embodiment, the filter of the first signal extraction unit 302 shown in FIG. 10 is a motor drive unit 12x, 12y. This is a point simulating the response of the position control system.
  • the transfer function of the position control unit 123 shown in FIG. 4 is Gp (s)
  • the response of the position control system can be approximated by the following equation (11).
  • An example of the simple position control unit 123 is a proportional controller represented by the following equation (12).
  • the relationship between the first signal extraction unit 302 and the second signal extraction unit 303 is as described above ( Since it is necessary to satisfy the equation (6), the transfer function of the second signal extraction unit 303 may be expressed by the following equation (13).
  • the acceleration component extracted by the second signal extraction unit 303 does not include a drive component.
  • the high frequency band components Pd xh and Pd yh extracted by the second signal extraction unit 313 shown in FIG. 11 become 0, the high frequency band components are used in the data calibration in the second data calibration unit 322 shown in FIG. There is no need to use Pd xh and Pd yh .
  • the average value of the displacement amounts Pxh and Pyh after the integration operation becomes 0, so that the error correction process is simplified.
  • the drive component included in the trajectory component that is, the component whose final value of the signal when integrated for one period is not necessarily zero is extracted as the low-frequency component. it can.
  • FIG. 17 is a diagram showing a machine motion trajectory measuring apparatus according to the fourth embodiment, a numerically controlled machine tool to which the machine motion trajectory measuring apparatus according to the fourth embodiment is applied, and a servo control apparatus.
  • FIG. 18 is a block diagram of the sensor signal separation unit shown in FIG. Differences between the first embodiment and the fourth embodiment are as follows. In the following description, the same parts as those in the first embodiment are denoted by the same reference numerals, and the description thereof is omitted. Only different parts will be described here. (1) The motion trajectory to be measured by the mechanical motion trajectory measuring apparatus 100 according to the fourth embodiment is an assumed motion trajectory between the tool 16 and the workpiece 17.
  • a numerically controlled machine tool 99 to which the machine motion trajectory measuring apparatus 100 according to the fourth embodiment is applied includes an object acceleration sensor 13a and a reference acceleration sensor 13b.
  • the object acceleration sensor 13a is an acceleration sensor installed on the work table 4 near the workpiece 17, and corresponds to the three-axis acceleration sensor 13 of the first embodiment.
  • the reference acceleration sensor 13 b is installed in the column 5 near the tool 16.
  • the column 5 is vibrated by the reaction force of the work table 4 and the tool 16 may vibrate due to the vibration of the column 5. Since vibrations generated during machining of the workpiece 17 cause machining errors, it is necessary to measure the relative displacement between the tool 16 and the workpiece 17 when vibration occurs in the tool 16.
  • the mechanical motion trajectory measuring apparatus 100 according to the fourth embodiment is close to the object acceleration sensor 13a installed on the work table 4 near the workpiece 17 and the tool 16.
  • a reference acceleration sensor 13b installed in the column 5 is used.
  • FIG. 18 is a block diagram of the sensor signal separation unit shown in FIG.
  • the difference between the sensor signal separation unit 30 of the first embodiment and the sensor signal separation unit 30 of the fourth embodiment is that the sensor signal separation unit 30 of the fourth embodiment includes the first noise removal unit 301, the first noise removal unit 301, In addition to the signal extraction unit 302 and the second signal extraction unit 303, a relative acceleration calculation unit 304 is provided.
  • the relative acceleration calculation unit 304 calculates the relative acceleration between the acceleration measured by the reference acceleration sensor 13b and the acceleration measured by the object acceleration sensor 13a in order to calculate the relative acceleration. That relative acceleration computing unit 304, from the difference between the measured acceleration component a x in the acceleration component measured by the reference acceleration sensor 13b A x and the object acceleration sensor 13a, obtains the relative acceleration component in the x-axis direction, the reference The relative acceleration component in the y-axis direction is obtained from the difference between the acceleration component Ay measured by the acceleration sensor 13b and the acceleration component ay measured by the object acceleration sensor 13a.
  • the relative acceleration calculated by the relative acceleration calculation unit 304 is input to the first noise removal unit 301.
  • the first noise removing unit 301 removes the noise components a xn and a yn from the relative acceleration components of the x-axis and the y-axis calculated by the relative acceleration calculating unit 304, and other than the noise components a xn and a yn
  • the signal components a xs and a ys are separated and output.
  • the signal components a xs, low frequency band component a xl of a ys, a yl is extracted
  • the second signal extraction section 303 the signal components a xs, frequency band components of a ys a xh and a yh are extracted.
  • the relative displacement amount can be calculated from the relative acceleration, the relative motion trajectory can be measured.
  • Embodiment 5 When the frequency of the generated mechanical vibration is less than the cut-off frequency of the first signal extraction unit 302, the motion trajectory is calculated by the method shown in the fourth embodiment, and when the integration error is compensated in the data calibration unit 32, It is not possible to distinguish between vibration components and integration errors. In such a case, the amplitude of the vibration component may be underestimated in the motion trajectory of the calculation result.
  • the reference acceleration sensor 13b and the object acceleration sensor 13a are used as in the fourth embodiment.
  • a sensor signal separation unit 30 configured similarly to the sensor signal separation unit 30 of the first embodiment is used.
  • the acceleration components A x and A y measured by the reference acceleration sensor 13b are converted into the high frequency band components A xh and A yh and the low frequency band components A xl and A yl . To be separated.
  • the acceleration components a x and a y measured by the object acceleration sensor 13a are divided into the high frequency band components a xh and a yh and the low frequency band components a xl and a yl. And separated.
  • FIG. 19 is a configuration diagram of a data calibration unit provided in the mechanical motion trajectory measuring apparatus according to the fifth embodiment.
  • the data calibration unit 32 illustrated in FIG. 19 includes a first reference data calibration unit 323 and a second reference data calibration unit 324 in addition to the first data calibration unit 321 and the second data calibration unit 322.
  • the first reference data calibration unit 323 receives the low frequency band components Ad xl and Ad yl of the reference acceleration separated by the sensor signal separation unit 30.
  • the second reference data calibration unit 324 receives the high frequency band components Ad xh and Ad yh of the reference acceleration separated by the sensor signal separation unit 30. Since the reference point is not a movable part but a fixed point, there is no need to use a motor detector signal as a reference signal.
  • the first reference data calibration unit 323 calculates reference point displacements PS xl and PS yl from the low frequency band components Ad xl and Ad yl of the reference acceleration.
  • the second reference data calibration unit 324 calculates reference point displacements PS xh and PS yh from the high frequency band components Ad xh and Ad yh of the reference acceleration.
  • a method for calculating the reference point displacements PS xl and PS yl will be described.
  • Axl (t) which is time series data, is converted into a frequency region of AXL (j ⁇ ) by Fourier transform.
  • j is an imaginary unit
  • is an angular frequency.
  • the integration operation in the frequency domain corresponds to calculating PSXL (j ⁇ ) shown in the following equation (14) at the angular frequency ⁇ .
  • PSXL (j ⁇ ) shown in the following equation (14) at the angular frequency ⁇ .
  • the frequency component is regarded as a noise component, and if the integration operation is not performed, the noise Accumulation of components can be prevented.
  • the low frequency band component PRxl of the relative displacement between the tool 16 and the workpiece 17 can be calculated from the difference between the reference point displacement PS xl and the low frequency band component PT xl .
  • the high frequency band component PRxh of the relative displacement can be calculated from the reference point displacement PS xh and the high frequency band component PT xh .
  • the relative displacement PRx between the tool 16 and the workpiece 17 can be calculated from the low frequency band component PRxl and the high frequency band component PRxh.
  • FIG. 20 is a diagram comparing the calculation result when the relative motion trajectory is measured by the mechanical motion trajectory measuring apparatus according to the fifth embodiment, the actual relative motion trajectory, and the motion trajectory at the detector position.
  • the actual driving locus is indicated by a one-dot chain line
  • the calculation result of the locus error calculated by the mechanical movement locus measuring apparatus 100 is indicated by a solid line
  • the driving locus at the detector position is indicated by a dotted line.
  • the relative displacement between the actual tool 16 and the workpiece 17 can be measured with an error of 1 micrometer, and the relative displacement can be measured with higher accuracy. is there.
  • the data calibration unit 32 of the mechanical motion trajectory measuring apparatus 100 performs the object motion trajectory for each frequency band from the acceleration sensor signal and the reference acceleration sensor signal separated by the sensor signal separation unit.
  • the calculation error of the object motion trajectory component and the reference motion trajectory component is corrected using the detection position signal for each frequency band separated by the detection position signal separation unit, and the target object is corrected.
  • the difference between the motion trajectory component and the reference motion trajectory component is output as a relative motion trajectory in that band. It is possible to calculate a relative motion trajectory with respect to the reference position by calculating the reference motion trajectory.
  • Embodiment 6 The configuration of the mechanical motion trajectory measuring apparatus 100 according to the sixth embodiment is the same as that of the first embodiment.
  • the difference from the first embodiment is that the velocity Vxl, Vyl calculated by integrating the low frequency band components a xl , a yl of the acceleration into the first rank and the low frequency band component of the motor signal when correcting the integration error.
  • the difference between Vdxl and Vdyl calculated by first-order differentiation of Pd xl and Pd yl is compensated as an integral error of speed, and Pxl ′, Pyl ′ and Pd xl calculated by further first-order integration after compensation. , Pd yl is compensated for.
  • error compensation can be performed before the integration error is accumulated. Therefore, even if the integration error cannot be compensated after the second-order integration with a large noise level, the integration error can be corrected with high accuracy.
  • the data calibration unit converts an input signal into a low frequency domain signal when performing high frequency band data calibration, performs an integration operation in the low frequency domain, and then converts the signal into the time domain by inverse transformation.
  • the output signal is output as a high frequency band component of the motion trajectory. It is possible to remove the integration error component by cutting and outputting the low frequency component caused by the integration error.
  • the data calibration unit of the sixth embodiment performs the first-order integration result of the acceleration sensor signal separated by the sensor signal separation unit and the first-order differentiation result of the low-frequency band component of the detection position signal at the time of data calibration in the low frequency band.
  • the error is corrected so that the difference between and does not exceed the allowable value.
  • Embodiment 7 FIG.
  • the configuration of the mechanical motion trajectory measuring apparatus 100 according to the seventh embodiment is the same as the configuration of the first embodiment.
  • the difference from the first embodiment is that the motion trajectory to be measured in the mechanical motion trajectory measuring apparatus 100 according to the seventh embodiment is not a circle but a quadrangle motion trajectory.
  • FIG. 21 shows a comparison between the actual trajectory when a square motion trajectory is commanded, the calculation result of the trajectory error calculated by the mechanical motion trajectory measuring apparatus according to the seventh embodiment, and the operation trajectory at the detector position.
  • FIG. FIG. 21 is an enlarged view of a corner portion when a square motion trajectory is commanded.
  • the actual trajectory is indicated by a solid line
  • the calculation result of the trajectory error calculated by the mechanical motion trajectory measuring apparatus 100 is indicated by a dotted line
  • the operation trajectory at the detector position is indicated by a one-dot chain line.
  • FIG. 21 shows that the actual trajectory and the motion trajectory of the calculation result are in good agreement.
  • the mechanical motion trajectory measuring apparatus 100 according to Embodiment 7 can be used not only for motion trajectory measurement in a two-dimensional plane but also for motion trajectory measurement in a three-dimensional space, and for any motion waveform.
  • the motion trajectory can be calculated with high accuracy.
  • FIG. FIG. 22 is a diagram showing a machine motion trajectory measuring apparatus according to the eighth embodiment, a numerically controlled machine tool to which the machine motion trajectory measuring apparatus according to the eighth embodiment is applied, and a servo control apparatus.
  • the difference from the second embodiment is that the machine motion trajectory measuring apparatus 100 according to the eighth embodiment has a parameter setting unit 500.
  • the parameter setting unit 500 notifies the sensor signal separation unit 30 and the motor signal separation unit 31 of a filter parameter setting command that is a filter design parameter.
  • the sensor signal separation unit 30 and the motor signal separation unit 31 generate a filter to be used based on the received filter parameter setting command. For example, when using a second-order low-pass filter to separate a low-frequency band signal in the first signal extraction unit 312, the parameters notified here are filter coefficients a, b.
  • the notified parameter is a filter cutoff which is a physical parameter expressing a filter characteristic formulated as shown in the following equation (16) instead of the filter coefficients a and b in the above equation (15). It may be the frequency ⁇ and the attenuation ⁇ .
  • the parameter setting unit 500 designs filters of the nth order shown in the following equation (17) in the sensor signal separation unit 30 and the motor signal separation unit 31, and transmits the filter parameters.
  • a filter having an arbitrary characteristic expressed by a transfer function up to the nth order is set by setting the high-order coefficients of the denominator and the numerator in equation (17) to 0 according to the required filter order. realizable.
  • the optimum filter characteristics may differ depending on the conditions. In such a case, the optimum filter can always be used by changing the parameter for each measurement condition from the parameter setting unit 500 every time the measurement condition is changed.
  • the parameter setting unit 500 sets the filter parameters, and notifies the sensor signal separation unit 30 and the motor signal separation unit 31 of the parameters, thereby The trouble of setting the filters of the signal separation unit 30 and the motor signal separation unit 31 individually can be saved.
  • the filter characteristics can be changed efficiently when the filter settings are changed frequently.
  • FIG. 23 is a diagram showing a machine motion trajectory measuring apparatus according to the ninth embodiment, a numerically controlled machine tool to which the machine motion trajectory measuring apparatus according to the ninth embodiment is applied, and a servo control device.
  • FIG. 24 is a schematic view showing an example of a numerically controlled machine tool to which the machine motion trajectory measuring apparatus according to the ninth embodiment is applied. The difference from the first embodiment is that the three-axis acceleration sensor 13 is mounted on the rotary table 501 of the numerically controlled machine tool 99b.
  • the triaxial acceleration sensor 13 is installed such that one axial direction is a normal direction of the rotary table 501, that is, a direction in which a straight line (not shown) that passes through the center of rotation of the circular table extends, and the other axis is a tangential direction of rotation, that is, It is installed so as to coincide with the radial direction of the rotary table 501.
  • the orientation of the sensor does not change during one round of circular motion.
  • the orientation of the three-axis acceleration sensor 13 attached to the rotary table 501 changes with rotation.
  • FIG. 25 is a diagram showing an example of the measurement result of the normal direction acceleration at measured when the rotary table is driven.
  • FIG. 26 is a diagram illustrating an example of a measurement result of the tangential acceleration ar measured when the rotary table is driven.
  • the vertical axis in FIG. 25 represents the normal direction acceleration at, and the horizontal axis in FIG. 25 represents time.
  • the vertical axis in FIG. 26 represents the tangential acceleration ar, and the horizontal axis in FIG. 26 represents time.
  • the normal direction acceleration at is influenced by the centripetal acceleration, and a vibration error component is superimposed around the nominal centripetal acceleration during circular motion.
  • a vibration error component is superimposed around zero.
  • the motor drive unit 12y outputs a tangential direction component Pdr from the feedback position Pd y as shown in the following equation (18).
  • the feedback position Pdy represents the rotation angle of the table, and the unit is radians.
  • the motor drive unit 12y outputs the normal direction component Pdt of the motion locus from the distance Rcom as shown in the following equation (19).
  • the distance Rcom is the distance from the rotation center of the rotary table 501 to the attachment of the triaxial acceleration sensor 13.
  • the sensor signal separation unit 30 and the motor signal separation unit 31 separate frequency bands using the tangential acceleration ar, the normal acceleration at, the tangential component Pdr, and the normal component Pdt. I do.
  • the trajectory coupling unit 331 of the motion trajectory calculation unit 33 calculates a tangential direction trajectory PTr and a normal direction trajectory PTt.
  • the X-direction component PT x is calculated by the calculation formula shown in the following formula (20)
  • the Y-direction component PT y is calculated by the calculation formula shown in the following formula (21).
  • the motion trajectory calculation unit 33 calculates the trajectory error in the same manner as the above equations (9) and (10) shown in the first embodiment.
  • the machine motion trajectory measuring apparatus 100 it is possible to calculate the machine motion trajectory even in the machine configuration driven by the rotary table 501.
  • the machine motion trajectory measuring apparatus 100 can be applied to a mechanical apparatus such as an industrial machine, a robot, or a transporter that has one or more movable axes and drives an object to be controlled using an actuator. It is.
  • the configuration described in the above embodiment shows an example of the contents of the present invention, and can be combined with another known technique, or a part of the configuration can be used without departing from the gist of the present invention. It can be omitted or changed.

Landscapes

  • Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Mechanical Engineering (AREA)
  • Numerical Control (AREA)
  • Machine Tool Sensing Apparatuses (AREA)
  • Control Of Position Or Direction (AREA)

Abstract

機械運動軌跡測定装置(100)は、運動軌跡測定対象物の加速度を測定し、加速度センサ信号として出力する3軸加速度センサ(13)と、加速度センサ信号を、2つ以上の周波数帯域に分離するセンサ信号分離部(30)と、検出位置信号を、センサ信号分離部(30)と同じ周波数帯域に分離するモータ信号分離部(31)と、センサ信号分離部(30)で分離された加速度センサ信号と、モータ信号分離部(31)で分離された検出位置信号とを用いて、2つ以上の周波数帯域の各々において加速度センサ信号を較正し、2つ以上の周波数帯域の各々の運動軌跡成分を得るデータ較正部(32)と、2つ以上の周波数帯域の各々の運動軌跡成分を結合し、運動軌跡として出力する運動軌跡算出部(33)とを備える。

Description

機械運動軌跡測定装置
 本発明は、数値制御工作機械、産業用機械、ロボットまたは搬送機といった機械装置の運動軌跡を測定する機械運動軌跡測定装置に関する。
 サーボ制御装置は、位置検出器を用いて検出された駆動対象である被駆動体の位置が指令位置に一致するようにアクチュエータを用いて制御を行う装置である。数値制御工作機械、産業用機械、ロボットまたは搬送機のように、可動範囲が2次元平面または3次元空間のような多自由度を有する機械装置は、1自由度の軸と呼ばれるサーボ制御装置を複数有する。このような機械装置では、各軸の各々に取付けられたアクチュエータが各軸の被駆動体を駆動制御する。機械装置がそれらの軸運動を組み合わせることにより、多自由度の機械の運動が実現される。
 運動軌跡が指令された経路である指令軌跡に正確に追従するように行うサーボ制御は、軌跡制御または輪郭運動制御と呼ばれる。軌跡制御中に発生する摩擦、または軌跡制御中に発生する機械構造の振動といった外乱要因によって、各軸のサーボ制御装置に応答誤差が発生すると、被駆動体の運動軌跡が指令軌跡を外れるため、軌跡誤差が生じる。
 切削工具を用いて各軸の運動を加工対象の工作物に転写することによって形状を創生する数値制御工作機械においては、数十マイクロメートルの軌跡誤差であっても加工不良と判断されることがある。そのため、位置検出器を用いてフィードバック制御を行う場合においても、位置検出器の取付け位置が被駆動体の位置と一致していない場合、または被駆動体の位置が真の制御対象の位置、すなわち工具またはワークの位置と一致していない場合、位置検出器で正確に制御対象の位置を検出できない。そのため、制御対象の運動軌跡に軌跡誤差が生じる。
 運動軌跡に軌跡誤差が発生した場合には、機械の運動軌跡を測定して問題の原因を探り、機械の設計変更、またはサーボ制御装置の各種のパラメータ調整が行われる。各種のパラメータ調整を行うためには、指令軌跡からの運動軌跡のずれ量を数マイクロメートルから数百マイクロメートルの精度で測定する必要がある。特許文献1,2,3には機械運動軌跡の測定および検出の方法が開示されている。
 特許文献1に開示される運動精度試験方法では、変位検出器を介して結合された2つの鋼球の間の相対距離を一定に保つような円弧運動を行わせたとき、この変位を読み取ることにより、2つの鋼球の間の相対運動軌跡が測定される。この測定方法は、ボールバー法と呼ばれ、広く普及している。一般的なボールバー法では、2つの鋼球の間を機械的に接続し、運動中のバーの変位を変位センサで測定することにより、円弧運動時の軌跡誤差を数マイクロメートルの精度で測定できる。
 特許文献2に開示される数値制御装置は、制御対象の速度および加速度の何れかを検出し、検出した速度および加速度の何れかを、予定された速度または加速度と比較することにより、軌跡誤差の発生を検出する。また特許文献2には、検出した軌跡誤差の速度または加速度成分をフィードバックすることにより、軌跡誤差を抑制する方法が開示されている。
 特許文献3には、1つ以上の軸を用いて正弦波運動を行ったときの被駆動体の加速度を2階積分することにより、被駆動体の運動軌跡を推定する方式が開示されている。さらに特許文献3には、積分の際に生じる誤差がしきい値を超えないように感度係数を変更することにより、高精度に運動軌跡を推定する方式が開示されている。
特開昭61-209857号公報 特開昭61-007907号公報 特開2015-182141号公報
 しかしながら、特許文献1に開示された運動精度試験方法では、測定可能な軌跡が円弧に限られるという課題がある。また、特許文献1に開示された運動精度試験方法では、装置のセットアップのために、回転中心の位置を予め計算しておき、片方の球の設置位置と軌跡測定時の回転中心とを一致させてセットアップする必要があり、セットアップに手間がかかるという課題がある。
 特許文献2に開示された軌跡誤差の抑制方法を用いれば、測定可能な軌跡の制限はない。しかしながら特許文献2に開示された軌跡誤差の抑制方法では、加速度または速度を用いて誤差の発生を検知するため、被駆動体の指令軌跡からの軌跡誤差量を求めることができないという課題がある。また特許文献2に開示された軌跡誤差の抑制方法では、加速度または速度をフィードバックすることにより、軌跡誤差を抑制することができる。しかしながら、特許文献2に開示された軌跡誤差の抑制方法は、加速度または速度をフィードバックすることで軌跡誤差を抑制できるが、その他の軌跡誤差抑制方式のパラメータ調整に用いることができないという課題がある。その他の軌跡誤差抑制方式としては、摩擦モデルを用いたロストモーショ、またはスティックモーションとよばれる、摩擦に由来する軌跡誤差を抑制する補正方式である。
 特許文献3に開示された運動軌跡の推定方式では、測定可能な軌跡が、正弦波または正弦波を組み合わせて実現可能な円弧軌跡、楕円軌跡、および球軌跡に限られるという課題があった。
 さらに数値制御工作機械およびロボットにおいては、可動側と固定側の2点間の相対運動軌跡の精度が重要となる場合がある。具体的に説明すると、機械構造物である鋳物の固定側に工具が取付けられており、鋳物の上に各軸のサーボ制御装置と工作物とが取付けられるタイプの機械では、工作物を輪郭制御し、工具と工作物とを干渉させることにより、工作物から材料を除去するという加工が行われる。このような場合、可動部である軸の駆動力の反作用で生じる駆動反力が、機械構造物を経由して固定側に伝達され、固定側に変形または振動が発生することがある。そのため可動側と固定側との間の相対運動軌跡を測定する必要がある。ところが特許文献1から3に開示された方法では、上記固定側の相対軌跡誤差を考慮することができないため、測定精度が低いという課題があった。
 本発明は、上記に鑑みてなされたものであって、簡易なセットアップで、任意の運動軌跡における軌跡誤差を高精度に測定可能な機械運動軌跡測定装置を得ることを目的とする。
 上述した課題を解決し、目的を達成するために、本発明の機械運動軌跡測定装置は、アクチュエータを有し、被駆動体の運動軌跡が指令軌跡に追従するように、アクチュエータの位置または被駆動体の位置を検出する位置検出器から出力される検出位置信号をフィードバックし、アクチュエータを駆動する機械装置の運動軌跡を測定する機械運動軌跡測定装置であって、運動軌跡測定対象物の加速度を測定し、加速度センサ信号として出力する加速度センサと、加速度センサ信号を、2つ以上の周波数帯域に分離するセンサ信号分離部とを備える。機械運動軌跡測定装置は、検出位置信号を、センサ信号分離部と同じ周波数帯域に分離する検出位置信号分離部と、センサ信号分離部で分離された加速度センサ信号と、検出位置信号分離部で分離された検出位置信号とを用いて、2つ以上の周波数帯域の各々において加速度センサ信号を較正し、2つ以上の周波数帯域の各々の運動軌跡成分を得るデータ較正部と、2つ以上の周波数帯域の各々の運動軌跡成分を結合し、運動軌跡として出力する運動軌跡算出部とを備えることを特徴とする。
 本発明によれば、簡易なセットアップで、任意の運動軌跡における軌跡誤差を高精度に測定できる、という効果を奏する。
実施の形態1に係る機械運動軌跡測定装置と数値制御工作機械との斜視図 図1に示すx軸駆動機構の詳細構成と、x軸駆動機構の動作を制御するサーボ制御装置の構成とを示す図 実施の形態1に係る機械運動軌跡測定装置と、実施の形態1に係る機械運動軌跡測定装置が適用される数値制御工作機械と、サーボ制御装置とを示す図 図3に示すモータ駆動部の構成図 検出器位置における運転軌跡と制御対象の運動軌跡と円弧運動時の指令軌跡とをxy平面でプロットした図 ノイズが存在しない場合のx軸方向の加速度波形と、その加速度から算出したx軸方向の変位量とを示す図 図6と同様の条件でxy平面における軌跡誤差を計算した結果を示す図 ノイズが存在する場合のx軸方向の加速度波形と、その加速度から算出したx軸方向の変位量とを示す図 図8と同様の条件でxy平面における軌跡誤差を計算した結果を示す図 図3に示すセンサ信号分離部の構成図 図3に示すモータ信号分離部の構成図 図3に示すデータ較正部の構成図 図3に示す運動軌跡算出部の構成図 実際の運動軌跡と、実施の形態1に係る機械運動軌跡測定装置で計算された軌跡誤差の計算結果と、検出器位置における運転軌跡とを比較して示す図 実施の形態2に係る機械運動軌跡測定装置と、実施の形態2に係る機械運動軌跡測定装置が適用される数値制御工作機械と、サーボ制御装置とを示す図 図15に示す機械応答模擬部の構成図 実施の形態4に係る機械運動軌跡測定装置と、実施の形態4に係る機械運動軌跡測定装置が適用される数値制御工作機械と、サーボ制御装置とを示す図 図17に示すセンサ信号分離部の構成図 実施の形態5に係る機械運動軌跡測定装置が備えるデータ較正部の構成図 実施の形態5に係る機械運動軌跡測定装置で相対運動軌跡を測定した場合の計算結果と実際の相対運動軌跡と検出器位置における運動軌跡を比較した図 四角形の運動軌跡を指令したときの実軌跡と、実施の形態7に係る機械運動軌跡測定装置で計算された軌跡誤差の計算結果と、検出器位置における運転軌跡とを比較して示す図 実施の形態8に係る機械運動軌跡測定装置と、実施の形態8に係る機械運動軌跡測定装置が適用される数値制御工作機械と、サーボ制御装置とを示す図 実施の形態9に係る機械運動軌跡測定装置と、実施の形態9に係る機械運動軌跡測定装置が適用される数値制御工作機械と、サーボ制御装置とを示す図 実施の形態9に係る機械運動軌跡測定装置が適用される数値制御工作機械の一例を示す概略図 ロータリテーブルを駆動したときに測定される法線方向加速度の測定結果の一例を示す図 ロータリテーブルを駆動したときに測定される接線方向加速度の測定結果の一例を示す図
 以下に、本発明の実施の形態に係る機械運動軌跡測定装置を図面に基づいて詳細に説明する。なお、この実施の形態によりこの発明が限定されるものではない。
実施の形態1.
 図1は実施の形態1に係る機械運動軌跡測定装置と数値制御工作機械との斜視図である。図2は図1に示すx軸駆動機構の詳細構成と、x軸駆動機構の動作を制御するサーボ制御装置の構成とを示す図である。
 図1に示す数値制御工作機械99は、機械運動軌跡測定装置が機械運動軌跡を測定する対象の一例である。数値制御工作機械99は、直交3軸の立型工作機械であり、Cコラム構造と呼ばれる機械構造を有する工作機械である。数値制御工作機械99は、架台21と、架台21上に設置されy軸方向に駆動されるサドル24と、サドル24上に設置されるワークテーブル4と、架台21に固定され架台21の上方へ延びるコラム5とを備える。コラム5にはラム6が取付けられ、ワークテーブル4上には、運動軌跡測定対象物である工作物17が設置される。
 また数値制御工作機械99は、サドル24に取付けられワークテーブル4をx軸方向に駆動するアクチュエータであるx軸駆動機構15xと、架台21に取付けられサドル24をy軸方向に駆動するアクチュエータであるy軸駆動機構15yと、コラム5に取付けられラム6をz軸方向に駆動するアクチュエータであるz軸駆動機構15zとを備える。
 x軸駆動機構15xは、x軸モータ1xと、x軸モータ1xにより駆動される送り軸である送りねじ2xと、送りねじ2xの回転角度を検出する回転角検出器3xとを備える。y軸駆動機構15yは、y軸モータ1yと、y軸モータ1yにより駆動される送り軸である送りねじ2yと、送りねじ2yの回転角度を検出する回転角検出器3yとを備える。z軸駆動機構15zは、z軸モータ1zと、z軸モータ1zにより駆動される送りねじ2zと、送りねじ2zの回転角度を検出する回転角検出器3zとを備える。
 x軸駆動機構15xによりワークテーブル4が駆動され、y軸駆動機構15yによりサドル24とその上部に設置されたx軸駆動機構15xとが駆動される。コラム5に取付けられたz軸駆動機構15zによりラム6、主軸7および工具16が駆動され、工具16により工作物17が加工される。結果として工作物17のxy平面内の2自由度運動と、主軸7の先端に取付けられた工具16のz軸方向の1自由度の運動とを組み合わせて、xyzの3次元空間内、すなわち3自由度において、工具16と工作物17とが干渉した部分である工作物17の表面の材料が除去される。これにより3次元形状が創成される。3つの回転角検出器3x,3y,3zの各々により検出されたモータ回転角度は、図2に示すサーボ制御装置101にフィードバックされる。
 図2には、図1に示すx軸駆動機構15xの詳細構成が示され、被駆動体のx軸方向における位置を制御するためのサーボ制御装置101が示される。図2には、図1に示す3つの駆動機構の内、x軸駆動機構15xの構成のみが示される。y軸駆動機構15yおよびz軸駆動機構15zの構成は、x軸駆動機構15xの構成と同様である。ただし、x軸駆動機構15x、y軸駆動機構15yおよびz軸駆動機構15zは以下の点が異なる。すなわち、x軸駆動機構15xの制御対象である被駆動体は、ワークテーブル4であるのに対して、y軸駆動機構15yの制御対象である被駆動体は、サドル24とサドル24の上部に取付けられたx軸駆動機構15xとであり、z軸駆動機構15zの制御対象である被駆動体は、コラム5とコラム5に取付けられた主軸7とであることが異なる。
 図2に示すx軸駆動機構15xは、サドル24と、各々の外輪がサドル24に固定され、各々の内輪が送りねじ2xを回転自在に支持する2つのサポートベアリング10xとを備える。またx軸駆動機構15xは、送りねじ2xに噛み合わされたナット9xと、ナット9xによって送りねじ2xの軸線方向へ移動するワークテーブル4と、x軸モータ1xの回転時に送りねじ2xの一端を固定するカップリング8xとを備える。
 x軸モータ1xの回転運動はカップリング8xを介して送りねじ2xに伝達され、送りねじ2xの回転運動はナット9xにより直線運動に変換される。これによりナット9xに固定されたワークテーブル4が送りねじ2xの軸方向に沿って移動する。
 送りねじ2xはサポートベアリング10xに支持され、x軸モータ1xおよびサポートベアリング10xはサドル24に支持され、回転角検出器3xはx軸モータ1xに支持されている。
 被駆動体のx軸方向における位置を制御するためのサーボ制御装置101は、被駆動体のx軸方向における位置を制御する位置指令を生成する位置指令生成部11xと、位置指令生成部11xで生成された位置指令によりx軸モータ1xの回転角を制御するモータ駆動部12xとを備える。図2には、被駆動体のx軸方向における位置を制御するためのサーボ制御装置101のみ示すが、被駆動体のy軸方向およびz軸方向の各々における位置を制御するためのサーボ制御装置も同様に構成されている。
 位置指令生成部11xで生成された位置指令はモータ駆動部12xに送信され、位置指令を受信したモータ駆動部12xは、回転角検出器3xにより検出されたモータ回転角度に、送りねじ2xのねじピッチを乗じて得られる検出位置と位置指令との誤差が小さくなるように、x軸モータ1xの回転角を制御する。
 なお数値制御工作機械99の可動軸には、x軸モータ1xおよび送りねじ2xの代わりにリニアモータを用いる場合もある。また数値制御工作機械99の可動軸には、回転角検出器3xの代わりにワークテーブル4の位置を直接検出することが可能なリニアエンコーダまたはレーザ変位計を用いる場合もある。
 数値制御工作機械99においては、工具16と工作物17との相対変位が重要である。輪郭制御をしながら加工している最中に、工具16と工作物17との間に相対変位が生じた場合、工作物17において材料の削り残しおよび削りすぎが発生すると、加工誤差が生じるからである。サーボ制御装置101では、このような加工誤差が発生しないようにするため、モータ駆動部12xにおいて、フィードバック制御を行っている。
 しかしながら検出器位置を工具16または工作物17と厳密に一致させることはできないため、検出できない軌跡誤差が生じる場合がある。また軌跡制御中に発生する摩擦、または軌跡制御中に発生する機械構造の振動といった動的な外乱によって、モータ駆動部12xのフィードバック制御が追いつかない場合、すなわちサーボ制御装置101に応答誤差が発生した場合、被駆動体の運動軌跡が指令軌跡を外れるため、工具16と工作物17との間に相対変位が生じる。
 動的な外乱としては、カップリング8x、送りねじ2xまたはサポートベアリング10xで生じる弾性変形と、送りねじ2xまたはサポートベアリング10xで生じる振動と、図1に示すコラム5またはラム6の姿勢変化と、コラム5またはラム6の振動と、可動軸の摩擦力による誤差といったものが知られている。これらの動的な外乱は、工作物17の質量と、機械の経年変化と、送りねじ2xまたはサポートベアリング10xの摩耗と、各可動軸の潤滑油量と、気温の変化と、工場内の他の機械の稼働状況とによって変化することが知られている。このような要因により動的な外乱が変化するため、高精度な輪郭制御性能を維持するために位置指令生成部11xおよびモータ駆動部12xには、摩擦および振動モデルを用いて、発生する誤差量を予見し、誤差を打ち消すための補正量を生成し、フィードフォワード制御するというモデルベースの誤差補正方式を用いることが多い。
 しかしながら、動的な誤差に影響する因子は多く、それらすべての因子を考慮したモデルを構築するのは現実的ではない。そこで作業開始前または仕上げ加工直前に、定期的または継続的に機械の運動軌跡を測定し、各種補正パラメータをその都度調整できることが望ましい。
 図1において、運動軌跡を測定する対象物である工作物17の近くには、工作物17の加速度を測定する対象物加速度センサの一例である3軸加速度センサ13が取付けられている。3軸加速度センサ13は、直交する3軸方向の加速度を各々測定する。3軸加速度センサ13と機械運動軌跡測定装置100はケーブル40を介して相互に接続される。3軸加速度センサ13で検出された3軸方向の各々の加速度を示す加速度センサ信号は、ケーブル40を介して、機械運動軌跡測定装置100に入力される。以下では、3軸加速度センサ13で検出された3軸方向の各々の加速度の成分を、x軸方向の加速度成分、y軸方向の加速度成分およびz軸方向の加速度成分と称する。
 3軸加速度センサ13は、1つのセンサで直交する3軸方向の加速度を測定できるため、1つの3軸加速度センサ13を用いることにより、3次元内における加速度を測定できる。ただし、本実施の形態では、3軸加速度センサ13の代わりに、3軸方向の各々の加速度を測定する1軸加速度センサを3つ用いてもよい。
 3軸加速度センサ13の取付け方法には、マグネットによる磁力を用いた取付け、ジグおよびねじを用いた締結、ワックスによる固定、または接着剤を用いた固定といった方法がある。測定対象物であるワークテーブル4に固定するだけで測定が可能なため、測定直前に簡単な作業で取付けることができる。また加速度センサの取付けに際しては特別な調整作業が不要なため、初習の作業者でも簡単に取付けることができる。
 近年は、無線を用いたケーブルレスの加速度センサも増加しており、無線の加速度センサを使用した場合、ケーブル40の取り回しを考慮する必要がないため、より簡易に設置できる。また、加工中の切りくずまたは切削油によって、3軸加速度センサ13およびケーブル40が破損することを防止するために、ワークテーブル4の裏側または内部に、3軸加速度センサ13を常時取付けておいてもよい。なお機械運動軌跡測定装置100は、数値制御工作機械99の外部に実装してもよいし、数値制御工作機械99の内部に実装してもよい。
 次にサーボ制御装置101によってx軸およびy軸を駆動したときの摩擦による運動軌跡への影響を調査する際、数値制御工作機械99をセットアップする構成例を説明する。
 図3は実施の形態1に係る機械運動軌跡測定装置と、実施の形態1に係る機械運動軌跡測定装置が適用される数値制御工作機械と、サーボ制御装置とを示す図である。図4は図3に示すモータ駆動部の構成図である。
 図3に示すサーボ制御装置101と図2に示すサーボ制御装置101との異なる点は、図3に示すサーボ制御装置101は、位置指令生成部11xおよびモータ駆動部12xに加えて、位置指令生成部11yおよびモータ駆動部12yを備えることである。
 モータ駆動部12xには、位置指令生成部11xで生成された位置指令Pixが入力され、回転角検出器3xで検出されたモータ回転角度であるフィードバック位置Pdが入力される。モータ駆動部12xは、位置指令Pixおよびフィードバック位置Pdに基づき、x軸モータ1xの駆動指令であるトルク指令Icxを出力する。
 モータ駆動部12yには、位置指令生成部11yで生成された位置指令Piyが入力され、回転角検出器3yで検出されたモータ回転角度であるフィードバック位置Pdが入力される。モータ駆動部12yは、位置指令Piyおよびフィードバック位置Pdに基づき、y軸モータ1yの駆動指令であるトルク指令Icyを出力する。
 図4には図3に示すモータ駆動部12xの構成が示される。モータ駆動部12xは、位置指令Pixと応答位置であるフィードバック位置Pdとの応答誤差である位置偏差を求める加減算部124aと、加減算部124aが求めた位置偏差に対して制御を実行し速度指令を生成する位置制御部123とを備える。またモータ駆動部12xは、フィードバック位置Pdの微分演算を実行し、速度フィードバック値を計算する微分演算部122と、位置制御部123で求められた速度指令と微分演算部122で求められた速度フィードバック値との速度偏差を求める加減算部124bと、駆動指令であるトルク指令Icxを出力する速度制御部121とを備える。
 加減算部124aは、位置指令Pixとフィードバック位置Pdとの偏差である位置偏差を求め、位置制御部123へ出力する。位置制御部123は、加減算部124aから出力される位置偏差を小さくするように、比例制御、PI(Proportional Integral)制御またはPID(Proportional Integral Differential)制御といった位置制御処理を実行し、位置偏差を小さくする速度指令を出力する。
 微分演算部122では、フィードバック位置Pdを微分した実速度が求められる。加減算部124bは、位置制御部123で求められた速度指令と微分演算部122で求められたフィードバック位置Pdの実速度との偏差である速度偏差を求め、速度制御部121へ出力する。速度制御部121では、加減算部124bから出力される速度偏差を小さくするように、PI制御の速度制御処理が行われ、トルク指令Icxが出力される。なお図3に示すモータ駆動部12yは、図4に示すモータ駆動部12xと同様に構成されているものとする。
 フィードバック位置Pd,Pdを求めるために回転角検出器3x,3yを使用する制御をセミクローズドループ制御と呼ぶ。図2に示すサーボ制御装置101およびx軸駆動機構15xによるセミクローズドループ制御では、回転角検出器3xの取付け位置と真の制御対象であるワークテーブル4との間、または回転角検出器3xの取付け位置と図1に示す工作物17との間に、送りねじ2xおよびナット9xといった機械構造物が存在する。そのため、回転角検出器3xとワークテーブル4との間、または回転角検出器3xと工作物17との間には、機械的な伝達遅れが生じる。従って、回転角検出器3xで検出されたフィードバック位置Pdは、ワークテーブル4または工作物17の運動軌跡と一致しない。
 また運動方向反転位置では、加わる摩擦の向きが反転することにより、可動軸には急激な外乱力が作用する。そのためフィードバック制御では、発生する位置偏差を完全に抑制することはできず、誤差が発生することが知られている。典型的な例は、円弧運動時のロストモーションまたはスティックモーションと呼ばれる軌跡誤差である。x軸とy軸を用いた円弧運動では、x軸モータ1xとy軸モータ1yとの間に、位相が90度ずれた正弦波指令を指令する。
 図5は検出器位置における運転軌跡と制御対象の運動軌跡と円弧運動時の指令軌跡とをxy平面でプロットした図である。検出器位置とは、前述したフィードバック位置Pdおよびフィードバック位置Pdを表す。図5では、検出器位置における運転軌跡が点線で示され、制御対象の運動軌跡が実線で示され、円弧運動時の指令軌跡が一点鎖線で示される。なお図5では軌跡誤差を半径方向に400倍に拡大して表示している。指令軌跡は円であるのに対して、検出器位置における運動軌跡と制御対象である工作物17の運動軌跡とは、象限切り替え位置で突起状の誤差パターンを示している。
 検出器位置では、象限切り替え位置において、運動軌跡が円の外側に飛び出した形状をしており、スティックモーションと呼ばれる。スティックモーションは、摩擦の加わる方向の反転によってフィードバック制御に遅れが生じるために発生することが明らかになっている。
 制御対象の運動軌跡は、象限切り替え位置で外側に飛び出した後、遅れて内側に食い込むような軌跡を示している。これは、前述の摩擦によるスティックモーションに加えて、送りねじ2xおよびナット9xといった機械要素が介在するために、制御対象の応答が遅れることによって発生するロストモーションと呼ばれる現象である。
 ロストモーションによる軌跡誤差を抑制するためには、位置指令生成部11x,11yまたはモータ駆動部12x,12yにおいて、発生するロストモーションを打ち消すように補正指令を加算する必要がある。そのためには、工作物17の運動軌跡を正確に測定する必要がある。
 図1および図3に示すように、本実施の形態に係る機械運動軌跡測定装置100が適用される数値制御工作機械99では、ワークテーブル4上の工作物17の近くに3軸加速度センサ13が設置されている。機械運動軌跡測定装置100には、3軸加速度センサ13で測定されたx軸方向の加速度成分aと、3軸加速度センサ13で測定されたy軸方向の加速度成分aと、x軸モータ1xの検出器位置であるフィードバック位置Pdと、y軸モータ1yの検出器位置であるフィードバック位置Pdとが入力される。ただし、摩擦による運動誤差は、可動平面内で完結することが知られているため、摩擦補正に関するパラメータを調整する場合、z軸方向の加速度成分については無視できる。なお、以下では、加速度成分aおよび加速度成分aを単に、加速度センサ信号と称する場合がある。またフィードバック位置Pd,Pdを検出器位置と称する場合がある。
 図3に示す機械運動軌跡測定装置100は、x軸方向の加速度成分aおよびy軸方向の加速度成分aの各々を、2つ以上の周波数帯域に分離するセンサ信号分離部30と、フィードバック位置Pd,Pdの各々を、センサ信号分離部30と同一の2つ以上の周波数帯域に分離する検出位置信号分離部であるモータ信号分離部31とを備える。
 また機械運動軌跡測定装置100は、センサ信号分離部30およびモータ信号分離部31の各々の分離された周波数帯域毎に、フィードバック位置Pd,Pdを参照信号として加速度成分a,aを較正し、各々の周波数帯域における位置応答を計算して出力するデータ較正部32と、データ較正部32で計算された各々の周波数帯域毎の位置応答を加算し、x軸およびy軸に稼働時の工作物17の運動軌跡を算出して運動軌跡を示す情報として出力する運動軌跡算出部33とを備える。
 運動軌跡算出部33で計算された運転軌跡を示す情報は、機械運動軌跡測定装置100に接続された軌跡表示装置35に出力される。軌跡表示装置35はビデオモニタを代表とする表示器であり、表示器の画面には軌跡表示装置35で計算された運転軌跡が表示される。これにより運転軌跡が実験者に提示される。このとき、軌跡表示装置35は、フィードバック位置Pd,Pdにおける運動軌跡および指令軌跡を入力し、図5に示すように指令軌跡と複数位置における運動軌跡との測定結果を重ねて表示してもよい。
 以下、加速度センサ信号から運動軌跡を測定する方式について述べる。xy平面で、半径を2mmとし、送り速度を500mm/minとして円弧軌跡運動を繰り返したときの運動軌跡測定を考える。ワークテーブル4上に設置された3軸加速度センサ13により工作物17の運動軌跡が測定される。このとき、3軸加速度センサ13で測定された加速度成分aを2階積分すれば、3軸加速度センサ13の取付け位置における運動軌跡である変位量Pxが算出できる。
 加速度成分aのデータが離散的なサンプリングによって取得された場合、速度Vxは、加速度成分aとサンプリング時間dtとを用いて、下記(1)式に示す数列から算出できる。さらに、変位量Pxは、速度Vxとサンプリング時間dtとを用いて、下記(2)式に示す数列より計算される。ただし速度VxおよびPx初期値は0にする。同様にy軸についても加速度から変位量Pyを算出できるので、変位量Px,Pyから運動軌跡を算出できる。
Figure JPOXMLDOC01-appb-M000001
Figure JPOXMLDOC01-appb-M000002
 図6はノイズが存在しない場合のx軸方向の加速度波形と、その加速度から算出したx軸方向の変位量とを示す図である。紙面上側の図では、横軸は時間を表し、縦軸は加速度を表す。紙面下側の図では、横軸は時間を表し、縦軸は変位量を表す。図6に示す加速度および変位量のそれぞれの波形は、加速度のサンプリング時間dtを1msとしたときのものである。
 図7は図6と同様の条件でxy平面における軌跡誤差を計算した結果を示す図である。図7では、実際の運転軌跡である実軌跡は点線で示され、上記(1)式および上記(2)式を用いて計算した軌跡誤差の計算結果は実線で示される。ノイズが存在しない理想的な状態では実軌跡と計算結果が一致するため、図7では点線の実軌跡が実線の計算結果と重複している。このようにノイズが存在しない理想的な条件では、簡単な計算で加速度から運動軌跡を算出できる。
 しかしながら、実際の測定環境では、電磁ノイズの影響、量子化誤差、またはエイリアシングの影響により、3軸加速度センサ13で測定された加速度センサ信号に様々なノイズ成分が重畳する。ノイズは、運動軌跡の加速度成分以外の信号成分の全てを指し、厳密には区別しない。また現実において完全にノイズを除去した真の運動軌跡の加速度成分を観測することは難しいため、図6では、理想的な加速度センサ信号をシミュレーションによって生成している。
 図8はノイズが存在する場合のx軸方向の加速度波形と、その加速度から算出したx軸方向の変位量とを示す図である。紙面上側の図では、横軸は時間を表し、縦軸は加速度を表す。紙面下側の図では、横軸は時間を表し、縦軸は変位量を表す。
 図9は図8と同様の条件でxy平面における軌跡誤差を計算した結果を示す図である。図9では、実際の運転軌跡である実軌跡は点線で示され、軌跡誤差の計算結果は実線で示される。図9より、計算によって算出した軌跡誤差は実軌跡と一致しないことがわかる。このような積分操作に起因する計算誤差を積分誤差と呼ぶ。加速度を積分して変位を算出する場合、加速度の積分誤差の影響を補償しないときには、数値制御工作機械99を代表とする機械の運動軌跡測定において要求される精度を実現することが難しい。
 本実施の形態に係る機械運動軌跡測定装置100では、センサ信号分離部30およびモータ信号分離部31が周波数帯域の分離を行い、データ較正部32が周波数帯域毎の加速度センサ信号により、積分誤差を低減する。以下にセンサ信号分離部30、モータ信号分離部31およびデータ較正部32の構成を具体的に説明する。
 図10は図3に示すセンサ信号分離部の構成図である。センサ信号分離部30は、加速度センサ信号である加速度成分a,aに含まれるノイズ成分を除去するために、加速度センサ信号をノイズ帯域に分離する。ノイズ帯域は加速度センサ信号のノイズ成分である。ノイズ成分を除去することで加速度補償時の誤差を小さくすることができる。具体的には、センサ信号分離部30は、加速度成分a,aの各々のノイズ成分axn,aynを除去し、ノイズ成分axn,ayn以外の信号成分axs,aysを分離して出力する第1のノイズ除去部301と、信号成分axs,aysの低周波帯域成分axl,aylを抽出する第1の信号抽出部302と、信号成分axs,aysの高周波帯域成分axh,ayhを抽出する第2の信号抽出部303とを備える。
 センサ信号分離部30に入力された加速度成分a,aは、第1のノイズ除去部301によって各々のノイズ成分axn,aynが分離される。これによりノイズ成分axn,ayn以外の信号成分axs,aysが抽出される。第1のノイズ除去部301は、下記(3)式で示される伝達関数Gfiltで記述されるローパスフィルタとして、センサ信号分離部30に実装される。下記(3)式のTfiltは、ローパスフィルタのカットオフ時定数である。ローパスフィルタを用いることで低域側の信号を効率的に分離できる。
Figure JPOXMLDOC01-appb-M000003
 第1のノイズ除去部301から出力された信号成分axs,aysの各々は、第1の信号抽出部302と第2の信号抽出部303によって、低周波帯域成分axl,aylと高周波帯域成分axh,ayhに分離される。
 第1の信号抽出部302は、下記(4)式で示される伝達関数Glで記述されるローパスフィルタである。下記(4)式のTcutoffは、カットオフ時定数である。ただし効果的にノイズ成分を除去するためには、Tfilt<Tcutoffとなっていることが好ましい。
Figure JPOXMLDOC01-appb-M000004
 第2の信号抽出部303は、下記(5)式で示される伝達関数Ghで記述されるハイパスフィルタである。
Figure JPOXMLDOC01-appb-M000005
 このとき、第1の信号抽出部302の伝達関数と第2の信号抽出部303の伝達関数との和が1となるように、各々のフィルタを設計すると、信号成分axs,aysの信号を過不足なく抽出できる。すなわち、周波数分離の前後で信号の過不足が生じないようにするためには、下記(6)式の関係式を満たすように各々のフィルタを設計する必要がある。すなわちセンサ信号分離部30は、低周波帯域用のフィルタとの伝達関数の和が1となるように設計された高周波帯域用のフィルタを有する。高周波帯域はアクチュエータを駆動する機械装置の振動周波数成分である。伝達関数の和を1とすることで、周波数分離を行う前後で、信号が失われまたは信号が増幅される帯域がなくなる。
Figure JPOXMLDOC01-appb-M000006
 図11は図3に示すモータ信号分離部の構成図である。モータ信号分離部31は、フィードバック位置Pd,Pdの各々のノイズ成分Pdxn,Pdynを除去し、ノイズ成分Pdxn,Pdyn以外の信号成分Pdxs,Pdysを分離して出力する第2のノイズ除去部311と、信号成分Pdxs,Pdysの低周波帯域成分Pdxl,Pdylを抽出する第1の信号抽出部312と、信号成分Pdxs,Pdysの高周波帯域成分Pdxh,Pdyhを抽出する第2の信号抽出部313とを備える。
 第2のノイズ除去部311、第1の信号抽出部312および第2の信号抽出部313の各々の伝達関数は、上記(3)式で示される伝達関数Gfilt、上記(4)式で示される伝達関数Gl、上記(5)式で示される伝達関数Ghに等しい。
 図12は図3に示すデータ較正部の構成図である。データ較正部32は、図10に示すセンサ信号分離部30で分離された低周波帯域成分axl,aylと図11に示すモータ信号分離部31で分離された低周波帯域成分Pdxl,Pdylとを用いて、運動軌跡の低周波帯域成分PTxl,PTylを算出する第1のデータ較正部321と、図10に示すセンサ信号分離部30で分離された高周波帯域成分axh,ayhと図11に示すモータ信号分離部31で分離された高周波帯域成分Pdxh,Pdyhとを用いて、運動軌跡の高周波帯域成分PTxh,PTyhを算出する第2のデータ較正部322とを備える。
 第1のデータ較正部321が低周波帯域成分PTxl,PTylを算出する際、加速度の低周波帯域成分axl,aylを数値積分して算出した変位量Pxl,Pylと、低周波帯域成分Pdxl,Pdylとの差分は、積分誤差であるので、この差分を補償すれば積分誤差を低減できる。積分誤差の補償方法としては、積分誤差を多項式で近似し、変位量Pxl,Pylから近似された積分誤差分を減算することにより、低周波帯域成分PTxl,PTylを算出する方法がある。積分誤差波形を近似する際には、最小二乗法により近似式の未知パラメータを決定してもよいし、滑降シンプレックス法のような数値解法を使ってもよい。
 第2のデータ較正部322における積分誤差の補償方法としては、低周波帯域における積分誤差の補償方法を用いてもよい。低周波帯域は、アクチュエータを駆動する機械装置の制御帯域である。また、積分誤差は、低次の次数で発生することが知られているので、第2のデータ較正部322は、数値積分により算出された高周波帯域成分Pdxh,Pdyhを、低周波帯域成分を除外するハイパスフィルタに通過させ、その結果を高周波帯域成分PTxh,PTyhとして出力する構成でもよい。
 データ較正部32で計算された低周波帯域成分PTxl,PTylおよび高周波帯域成分PTxh,PTyhは、図3に示す運動軌跡算出部33へ入力される。
 図13は図3に示す運動軌跡算出部の構成図である。運動軌跡算出部33は、軌跡結合部331と軌跡誤差計算部332とを備える。軌跡結合部331は、図12に示すデータ較正部32で計算された運動軌跡の低周波帯域成分PTxl,PTylと高周波帯域成分PTxh,PTyhとを結合する。すなわち、軌跡結合部331は、下記(7)式に基づき3軸加速度センサ13の取付け位置における運動軌跡のx軸方向成分PTを計算し、下記(8)式に基づき3軸加速度センサ13の取付け位置における運動軌跡のy軸方向成分PTを計算する。
Figure JPOXMLDOC01-appb-M000007
Figure JPOXMLDOC01-appb-M000008
 軌跡誤差計算部332は、軌跡結合部331で計算されたx軸方向成分PTおよびy軸方向成分PTに基づいて、x軸方向の軌跡誤差PT’とy軸方向の軌跡誤差PT’とを計算して出力する。円弧運動の軌跡誤差の表示方式としては、下記(9)式および(10)式に示されるように、半径方向に軌跡誤差を拡大して表示する方法が知られている。下記(9)式および(10)式に示されるRcomは、指令半径を表し、MAGは軌跡誤差の拡大表示倍率を表す。
Figure JPOXMLDOC01-appb-M000009
Figure JPOXMLDOC01-appb-M000010
 図14は実際の運動軌跡と、実施の形態1に係る機械運動軌跡測定装置で計算された軌跡誤差の計算結果と、検出器位置における運転軌跡とを比較して示す図である。図14では、実際の運動軌跡が一点鎖線で示され、軌跡誤差の計算結果が実線で示され、検出器位置における運転軌跡が点線で示される。図14に示すように、実施の形態1に係る機械運動軌跡測定装置100を用いることにより、1マイクロメートルの精度で工作物17の運動軌跡を算出することが可能である。
 以上に説明したように実施の形態1に係る機械運動軌跡測定装置100によれば、加速度センサ信号と検出位置信号とを2つ以上の帯域に分離して較正することにより、簡易なセットアップで任意の運動軌跡における軌跡誤差を高精度に測定可能である。
 また実施の形態1に係るデータ較正部は、センサ信号分離部で分離された加速度センサ信号から周波数帯域毎に運動軌跡成分の計算誤差を補正するとき、検出位置信号分離部で分離された周波数帯域毎の検出位置信号を用いるように構成されている。検出位置信号を参照信号として使用することで計算誤差を効果的に補正することが可能となる。
 また実施の形態1に係るデータ較正部は、データ較正の際に加速度の入力信号を2階積分し、積分結果の平均値が0となるように積分誤差を補正するように構成されている。参照位置が動かない場合、積分結果の平均値が0になるように積分誤差を補正すればよく、簡単に積分誤差の補償ができる。
 また実施の形態1に係るデータ較正部は、高周波帯域のデータ較正の際、積分結果として出力される信号に対して、センサ信号分離部で使用されたハイパスフィルタで通過しえない低域の周波数成分をカットして出力するように構成されている。積分誤差によって生じる低域の周波数成分をカットして出力することで積分誤差成分を除去することが可能となる。
 また実施の形態1に係るデータ較正部は、低周波帯域のデータ較正の際、センサ信号分離部で分離された加速度センサ信号の2階積分結果と検出位置信号の低周波帯域成分との差が許容値を超えないように誤差を補正し、低周波帯域における対象物の運動軌跡を計算するように構成されている。2階積分後に積分誤差を補正することで、最終結果に誤差成分が残らないようにできる。
 また実施の形態1に係る機械運動軌跡測定装置100は、機械装置の運動軌跡を表示する運動軌跡表示部を備え、運動軌跡表示部は、位置指令生成部に入力される目標位置から目標軌跡を合成し、位置指令生成部から出力される指令位置から指令軌跡を合成し、検出位置信号から検出部軌跡を合成し、目標軌跡、指令軌跡および検出部軌跡の何れか一の軌跡と機械装置の運動軌跡とを重ねて運動軌跡表示部に表示させる。表示装置を有することで、操作者に直感的に測定結果を表示できる。
実施の形態2.
 図15は実施の形態2に係る機械運動軌跡測定装置と、実施の形態2に係る機械運動軌跡測定装置が適用される数値制御工作機械と、サーボ制御装置とを示す図である。実施の形態1と実施の形態2との相違点は以下の通りである。なお以下では実施の形態1と同一部分には同一符号を付してその説明を省略し、ここでは異なる部分についてのみ述べる。
(1)実施の形態2に係る機械運動軌跡測定装置100は、センサ信号分離部30、モータ信号分離部31、データ較正部32および運動軌跡算出部33に加えて、機械応答模擬部34を備えること。
(2)実施の形態2に係る機械運動軌跡測定装置100には、x軸モータ1xのフィードバック位置Pdとy軸モータ1yのフィードバック位置Pdとの代わりに、位置指令生成部11xで生成された位置指令Pixと位置指令生成部11yで生成された位置指令Piyとが入力されること。
(3)実施の形態2に係る機械運動軌跡測定装置100のデータ較正部32は、機械応答模擬部34で計算されたフィードバック位置Pmを参照信号として加速度成分a,aを較正し、各々の周波数帯域における位置応答を計算すること。
 図16は図15に示す機械応答模擬部の構成図である。機械応答模擬部34は、図15に示すモータ駆動部12xと同様に構成された位置制御部123Aと、速度制御部121と、微分演算部122と、機械モデル341とを備える。
 機械モデル341の一例としては2慣性モデルおよび3慣性モデルがある。2慣性モデルは、慣性モデル可動軸の負荷イナーシャJの逆数1/Jまたはモータのイナーシャと被駆動体とを、2慣性の振動系で近似したモデルである。3慣性モデルは、モータイナーシャと被駆動体と送りねじのイナーシャとを3慣性の振動系で近似したモデルである。なお図16では、位置指令Pixを用いてフィードバック位置Pmを生成する機械応答模擬部34の構成例を示すが、位置指令Piyを用いてフィードバック位置Pmを生成する機械応答模擬部34も同様に構成され、この場合、位置制御部123は図15に示すモータ駆動部12yと同様に構成されている。
 実施の形態2に係る機械運動軌跡測定装置100によれば、機械モデル341で計算したフィードバック位置Pmを参照信号として使用することにより、回転角検出器3x,3yの分解能が低いため回転角検出器3x,3yから出力されるフィードバック位置を参照信号として使用できない場合でも、機械モデル341で計算したフィードバック位置を参照信号として代用できる。また実施の形態2に係る機械運動軌跡測定装置100は、セミクローズドループ制御のモータ駆動部12x,12yにおいて仮想的に被駆動体位置を計算し、計算した被駆動体位置を参照信号として利用できる。
実施の形態3.
 実施の形態1と実施の形態3との相違点は、実施の形態3に係る機械運動軌跡測定装置100では、図10に示す第1の信号抽出部302のフィルタが、モータ駆動部12x,12yの位置制御系の応答を模擬している点である。図4に示す位置制御部123の伝達関数をGp(s)とした場合、位置制御系の応答は、下記(11)式のように近似できる。簡単な位置制御部123の例は下記(12)式に示す比例制御器である。
Figure JPOXMLDOC01-appb-M000011
Figure JPOXMLDOC01-appb-M000012
 また図10に示すセンサ信号分離部30において、周波数分離の前後で信号の過不足が生じなくするためには、第1の信号抽出部302と第2の信号抽出部303の関係は、上記(6)式を満たす必要があるので、第2の信号抽出部303の伝達関数は下記(13)式のようにすればよい。
Figure JPOXMLDOC01-appb-M000013
 第1の信号抽出部302のローパスフィルタに位置制御系の応答を模擬したモデルを使用することにより、第2の信号抽出部303で抽出される加速度成分に駆動成分が含まれなくなる。また図11に示す第2の信号抽出部313において抽出される高周波帯域成分Pdxh,Pdyhが0になるため、図12に示す第2のデータ較正部322におけるデータ較正の際、高周波帯域成分Pdxh,Pdyhを使用する必要がなくなる。これにより積分操作後の変位量Pxh,Pyhの平均値が0になるため、誤差補正処理が簡単になる。またローパスフィルタを位置制御系の伝達関数と一致させることで、軌跡成分に含まれる駆動成分、すなわち1周期分積分した時の信号の最終値が必ずしも0とはならない成分を、低周波成分として抽出できる。
実施の形態4.
 図17は実施の形態4に係る機械運動軌跡測定装置と、実施の形態4に係る機械運動軌跡測定装置が適用される数値制御工作機械と、サーボ制御装置とを示す図である。図18は図17に示すセンサ信号分離部の構成図である。実施の形態1と実施の形態4との相違点は以下の通りである。なお以下では実施の形態1と同一部分には同一符号を付してその説明を省略し、ここでは異なる部分についてのみ述べる。
(1)実施の形態4に係る機械運動軌跡測定装置100の測定対象となる運動軌跡は、工具16と工作物17との間の想定運動軌跡であること。
(2)実施の形態4に係る機械運動軌跡測定装置100が適用される数値制御工作機械99は、対象物加速度センサ13aおよび基準加速度センサ13bを備えること。対象物加速度センサ13aは、工作物17の近くにおいてワークテーブル4に設置された加速度センサであり、実施の形態1の3軸加速度センサ13に相当する。基準加速度センサ13bは、工具16の近くにおいてコラム5に設置されている。
(3)実施の形態4に係る機械運動軌跡測定装置100では、対象物加速度センサ13aで測定された加速度成分aおよび加速度成分aに加えて、基準加速度センサ13bで測定された基準加速度センサ信号である加速度成分Aおよび加速度成分Aを用いていること。
 実際の数値制御工作機械99においては、ワークテーブル4の駆動反力でコラム5に振動が発生し、コラム5の振動により工具16が振動する時がある。工作物17の加工中に発生した振動は加工誤差の原因となるため、工具16に振動が発生した場合、工具16と工作物17との間の相対変位を測定する必要がある。このような相対変位の測定を行うため、実施の形態4に係る機械運動軌跡測定装置100は、工作物17の近くにおいてワークテーブル4に設置された対象物加速度センサ13aと、工具16の近くにおいてコラム5に設置された基準加速度センサ13bとを用いる。
 図18は図17に示すセンサ信号分離部の構成図である。実施の形態1のセンサ信号分離部30と実施の形態4のセンサ信号分離部30との相違点は、実施の形態4のセンサ信号分離部30は、第1のノイズ除去部301、第1の信号抽出部302および第2の信号抽出部303に加えて、相対加速度計算部304を備えることである。
 相対加速度計算部304では、相対加速度を算出するために、基準加速度センサ13bで測定された加速度と、対象物加速度センサ13aで計測された加速度との相対加速度が計算される。すなわち相対加速度計算部304は、基準加速度センサ13bで測定された加速度成分Aと対象物加速度センサ13aで測定された加速度成分aとの差分から、x軸方向における相対加速度成分を求め、基準加速度センサ13bで測定された加速度成分Aと対象物加速度センサ13aで測定された加速度成分aとの差分から、y軸方向における相対加速度成分を求める。相対加速度計算部304で計算された相対加速度は第1のノイズ除去部301に入力される。
 第1のノイズ除去部301は、相対加速度計算部304で計算されたx軸およびy軸の各々の相対加速度成分からノイズ成分axn,aynを除去し、ノイズ成分axn,ayn以外の信号成分axs,aysを分離して出力する。第1の信号抽出部302では、信号成分axs,aysの低周波帯域成分axl,aylが抽出され、第2の信号抽出部303では、信号成分axs,aysの高周波帯域成分axh,ayhが抽出される。実施の形態4によれば、相対加速度から相対変位量を計算できるため、相対運動軌跡の測定が可能になる。
実施の形態5.
 発生する機械振動の周波数が第1の信号抽出部302のカットオフ周波数未満の場合、実施の形態4に示した方法で運動軌跡を計算すると、データ較正部32における積分誤差の補償の際に、振動成分と積分誤差とを区別することができない。このような場合、計算結果の運動軌跡において振動成分の振幅が過少に評価される場合がある。
 なお実施の形態5では、実施の形態4と同様に基準加速度センサ13bおよび対象物加速度センサ13aが用いられる。一方、実施の形態5では、実施の形態1のセンサ信号分離部30と同様に構成されるセンサ信号分離部30が用いられる。ただし実施の形態5のセンサ信号分離部30では、基準加速度センサ13bで測定された加速度成分A,Aが、高周波帯域成分Axh,Ayhと低周波帯域成分Axl,Aylとに分離される。さらに、実施の形態5のセンサ信号分離部30では、対象物加速度センサ13aで測定された加速度成分a,aが、高周波帯域成分axh,ayhと低周波帯域成分axl,aylとに分離される。
 図19は実施の形態5に係る機械運動軌跡測定装置が備えるデータ較正部の構成図である。図19に示すデータ較正部32は、第1のデータ較正部321および第2のデータ較正部322に加えて、第1の基準データ較正部323および第2の基準データ較正部324を備える。
 第1の基準データ較正部323には、センサ信号分離部30で分離された基準加速度の低周波帯域成分Adxl,Adylが入力される。第2の基準データ較正部324には、センサ信号分離部30で分離された基準加速度の高周波帯域成分Adxh,Adyhが入力される。基準点は、可動部ではなく不動点であるので、モータ検出器信号を参照信号として用いる必要がない。第1の基準データ較正部323は、基準加速度の低周波帯域成分Adxl,Adylから基準点変位PSxl,PSylを計算する。第2の基準データ較正部324は、基準加速度の高周波帯域成分Adxh,Adyhから基準点変位PSxh,PSyhを計算する。
 基準点変位PSxl,PSylの計算方法について述べる。低周波帯域成分Axlについて、時系列データであるAxl(t)は、フーリエ変換によってAXL(jω)の周波数領域に変換される。jは虚数単位、ωは、角周波数である。このとき、周波数領域での積分操作は、角周波数ωにおいて下記(14)式に示したPSXL(jω)を計算することに相当する。PSXL(jω)を逆フーリエ変換して基準点変位PSxlを算出すれば、周波数領域での積分操作が行える。このとき、しきい値を設定しておき、フーリエ変換後のAXLの絶対値がしきい値未満となる場合、その周波数成分は、ノイズ成分であるとみなして、積分操作を行わなければ、ノイズ成分の蓄積を防止することができる。
Figure JPOXMLDOC01-appb-M000014
 基準点変位PSxlと低周波帯域成分PTxlとの差から、工具16と工作物17との相対変位の低周波帯域成分PRxlが算出できる。同様に、高周波帯域においても基準点変位PSxhと高周波帯域成分PTxhとから、相対変位の高周波帯域成分PRxhが算出できる。さらに低周波帯域成分PRxlと高周波帯域成分PRxhとから、工具16と工作物17との間の相対変位PRxが算出できる。
 図20は実施の形態5に係る機械運動軌跡測定装置で相対運動軌跡を測定した場合の計算結果と実際の相対運動軌跡と検出器位置における運動軌跡とを比較した図である。図20では、実際の運転軌跡が一点鎖線で示され、機械運動軌跡測定装置100で計算した軌跡誤差の計算結果が実線で示され、検出器位置における運転軌跡が点線で示される。実施の形態5に係る機械運動軌跡測定装置100によれば、実際の工具16と工作物17との間の相対変位を1マイクロメートルの誤差で測定でき、相対変位をより高精度に測定可能である。
 以上に説明したように実施の形態5に係る機械運動軌跡測定装置100のデータ較正部32は、センサ信号分離部で分離された加速度センサ信号および基準加速度センサ信号から周波数帯域毎に対象物運動軌跡成分および基準運動軌跡成分の計算を行うとき、検出位置信号分離部で分離された周波数帯域毎の検出位置信号を用いて対象物運動軌跡成分および基準運動軌跡成分の計算誤差を補正し、対象物運動軌跡成分と基準運動軌跡成分との差をその帯域における相対運動軌跡として出力する。基準運動軌跡を計算することで基準位置に対する相対運動軌跡を算出することが可能となる。
実施の形態6.
 実施の形態6に係る機械運動軌跡測定装置100の構成は、実施の形態1と同じである。実施の形態1との相違点は、積分誤差を補正する際に、加速度の低周波帯域成分axl,aylを1階数積分して算出した速度Vxl,Vylと、モータ信号の低周波帯域成分であるPdxl,Pdylを1階微分して算出したVdxl,Vdylとの差分を、速度の積分誤差として補償し、補償後にさらに1階積分して計算したPxl‘、Pyl‘と、Pdxl,Pdylとの差分を補償することである。これにより積分誤差が蓄積する前に誤差補償ができるため、ノイズレベルが大きく2階積分した後では、積分誤差を補償できない場合においても、高精度に積分誤差を補正できる。
 実施の形態6のデータ較正部は、高周波帯域のデータ較正の際、入力された信号を低周波数領域の信号に変換した後、低周波数領域で積分操作を行い、その後逆変換によって時間領域に変換した信号を運動軌跡の高周波帯域成分として出力するように構成されている。積分誤差によって生じる低域の周波数成分をカットして出力することで積分誤差成分を除去することが可能となる。
 また実施の形態6のデータ較正部は、低周波帯域のデータ較正の際、センサ信号分離部で分離された加速度センサ信号の1階積分結果と検出位置信号の低周波帯域成分の1階微分結果との差が許容値を超えないように誤差を補正するように構成されている。1階積分したところで積分誤差を補正することで次の積分操作で誤差の累積を最小限にできる。
実施の形態7.
 実施の形態7に係る機械運動軌跡測定装置100の構成は実施の形態1の構成と同じである。実施の形態1との相違点は、実施の形態7の機械運動軌跡測定装置100において測定対象となる運動軌跡は、円ではなく四角形の運動軌跡であることである。
 図21は四角形の運動軌跡を指令したときの実軌跡と、実施の形態7に係る機械運動軌跡測定装置で計算された軌跡誤差の計算結果と、検出器位置における運転軌跡とを比較して示す図である。図21には、四角形の運動軌跡を指令したときのコーナ部分が拡大して示される。図21では、実軌跡が実線で示され、機械運動軌跡測定装置100で計算した軌跡誤差の計算結果が点線で示され、検出器位置における運転軌跡が一点鎖線で示される。
 図21によれば、実軌跡と計算結果の運動軌跡とが良く一致することが分かる。また実施の形態7に係る機械運動軌跡測定装置100は、2次元平面内の運動軌跡測定だけに限らず、3次元空間内の運動軌跡測定にも使用することができ、任意の運動波形に対して高精度に運動軌跡を算出できる。
実施の形態8.
 図22は実施の形態8に係る機械運動軌跡測定装置と、実施の形態8に係る機械運動軌跡測定装置が適用される数値制御工作機械と、サーボ制御装置とを示す図である。実施の形態2との相違点は、実施の形態8に係る機械運動軌跡測定装置100がパラメータ設定部500を有する点である。
 パラメータ設定部500は、フィルタの設計パラメータであるフィルタパラメータ設定指令をセンサ信号分離部30とモータ信号分離部31に通知する。センサ信号分離部30とモータ信号分離部31は、受け取ったフィルタパラメータ設定指令に基づき、使用するフィルタを生成する。ここで通知されるパラメータは、例えば、第1の信号抽出部312において低周波帯域の信号を分離するために、2次のローパスフィルタを使用する場合、下記(15)式に示すフィルタ係数a,bである。
Figure JPOXMLDOC01-appb-M000015
 また前記の通知されるパラメータは、上記(15)式のフィルタ係数a,bに代えて、下記(16)式に示すように定式化されるフィルタ特性を表現する物理パラメータであるフィルタのカットオフ周波数ωと減衰ζであってもよい。
Figure JPOXMLDOC01-appb-M000016
 同様にパラメータ設定部500は、センサ信号分離部30とモータ信号分離部31に下記(17)式に示すn次の次数のフィルタを設計しておき、フィルタのパラメータを送信する。
Figure JPOXMLDOC01-appb-M000017
 この際に必要なフィルタ次数に応じて、上記(17)式における分母と分子それぞれの高次の係数を0と設定することで、n次まで伝達関数で表現される任意の特性を有するフィルタを実現できる。複数の測定条件で機械運動軌跡を測定する場合、条件によって最適なフィルタ特性が異なる場合がある。このような場合は、測定条件が変更される度に、パラメータ設定部500から測定条件毎にパラメータを変更することで常に最適なフィルタを使用できる。
 実施の形態8に係る機械運動軌跡測定装置100によれば、パラメータ設定部500でフィルタのパラメータを設定すると共に、当該パラメータをセンサ信号分離部30とモータ信号分離部31に通知することにより、センサ信号分離部30とモータ信号分離部31のフィルタを個別に設定する手間を省くことができる。また、フィルタ設定を頻繁に変更する場合にフィルタ特性の効率的な変更が可能となる。
実施の形態9.
 図23は実施の形態9に係る機械運動軌跡測定装置と、実施の形態9に係る機械運動軌跡測定装置が適用される数値制御工作機械と、サーボ制御装置とを示す図である。また、図24は実施の形態9に係る機械運動軌跡測定装置が適用される数値制御工作機械の一例を示す概略図である。実施の形態1との相違点は、数値制御工作機械99bのロータリテーブル501上に3軸加速度センサ13が取付けられている点である。
 3軸加速度センサ13は、1つの軸方向がロータリテーブル501の法線方向、すなわち円テーブルの回転中心を通る不図示の直線が延びる方向に設置され、もう1つの軸が回転の接線方向、すなわちロータリテーブル501の半径方向と一致するように設置される。
 実施の形態8までで使用された数値制御工作機械99のような直交2軸のテーブルに3軸加速度センサ13を取付ける場合は、円運動で1周周回する間にセンサの向きは変わらない。一方で、ロータリテーブル501に取付けられた3軸加速度センサ13は回転に伴ってセンサの向きが変わる。
 図25はロータリテーブルを駆動したときに測定される法線方向加速度atの測定結果の一例を示す図である。図26はロータリテーブルを駆動したときに測定される接線方向加速度arの測定結果の一例を示す図である。図25の縦軸は法線方向加速度atを表し、図25の横軸は時間を表す。図26の縦軸は接線方向加速度arを表し、図26の横軸は時間を表す。
 法線方向加速度atは、向心加速度の影響で、円運動中はノミナルな向心加速度を中心に振動の誤差成分が重畳する。一方で、接線方向加速度arは、ゼロを中心に振動の誤差成分が重畳する。モータ駆動部12yは、フィードバック位置Pdから下記(18)式のように接線方向成分Pdrを出力する。ただし、ロータリテーブルでは、フィードバック位置Pdyは、テーブルの回転角度を表し、単位はラジアンである。またモータ駆動部12yは、距離Rcomから下記(19)式のように運動軌跡の法線方向成分Pdtを出力する。ただし距離Rcomは、ロータリテーブル501の回転中心から3軸加速度センサ13の取付けまでの距離である。
Figure JPOXMLDOC01-appb-M000018
Figure JPOXMLDOC01-appb-M000019
 センサ信号分離部30とモータ信号分離部31は、実施の形態1と同様に、接線方向加速度ar、法線方向加速度at、接線方向成分Pdrおよび法線方向成分Pdtを用いて、周波数帯の分離を行う。運動軌跡算出部33の軌跡結合部331では、接線方向軌跡PTrと法線方向軌跡PTtが計算される。運動軌跡算出部33の軌跡誤差計算部332では、下記(20)式に示す計算式でX方向成分PTが算出され、下記(21)式に示す計算式でY方向成分PTが算出される。これにより運動軌跡算出部33では、実施の形態1に示した上記(9)式および上記(10)式と同様に軌跡誤差が計算される。
Figure JPOXMLDOC01-appb-M000020
Figure JPOXMLDOC01-appb-M000021
 実施の形態9に係る機械運動軌跡測定装置100によれば、ロータリテーブル501で駆動される機械構成においても機械運動軌跡を計算することが可能となる。
 なお実施の形態1から9に係る機械運動軌跡測定装置100は、1つ以上の可動軸を有しアクチュエータを用いて制御対象を駆動する産業用機械、ロボットまたは搬送機といった機械装置にも適用可能である。
 以上の実施の形態に示した構成は、本発明の内容の一例を示すものであり、別の公知の技術と組み合わせることもできるし、本発明の要旨を逸脱しない範囲で、構成の一部を省略、変更することもできる。
 1x x軸モータ、1y y軸モータ、1z z軸モータ、2x,2y,2z 送りねじ、3x,3y,3z 回転角検出器、4 ワークテーブル、5 コラム、6 ラム、7 主軸、8x カップリング、9x ナット、10x サポートベアリング、11x,11y 位置指令生成部、12x,12y モータ駆動部、13 3軸加速度センサ、13a 対象物加速度センサ、13b 基準加速度センサ、15x x軸駆動機構、15y y軸駆動機構、15z z軸駆動機構、16 工具、17 工作物、21 架台、24 サドル、30 センサ信号分離部、31 モータ信号分離部、32 データ較正部、33 運動軌跡算出部、34 機械応答模擬部、35 軌跡表示装置、40 ケーブル、99,99b 数値制御工作機械、100 機械運動軌跡測定装置、101 サーボ制御装置、121 速度制御部、122 微分演算部、123,123A 位置制御部、124a,124b 加減算部、301 第1のノイズ除去部、302,312 第1の信号抽出部、303,313 第2の信号抽出部、304 相対加速度計算部、311 第2のノイズ除去部、321 第1のデータ較正部、322 第2のデータ較正部、323 第1の基準データ較正部、324 第2の基準データ較正部、331 軌跡結合部、332 軌跡誤差計算部、341 機械モデル、500 パラメータ設定部、501 ロータリテーブル。

Claims (16)

  1.  アクチュエータを有し、被駆動体の運動軌跡が指令軌跡に追従するように、前記アクチュエータの位置または前記被駆動体の位置を検出する位置検出器から出力される検出位置信号をフィードバックし、前記アクチュエータを駆動する機械装置の運動軌跡を測定する機械運動軌跡測定装置であって、
     運動軌跡測定対象物の加速度を測定し、加速度センサ信号として出力する加速度センサと、
     前記加速度センサ信号を、2つ以上の周波数帯域に分離するセンサ信号分離部と、
     前記検出位置信号を、前記センサ信号分離部と同じ周波数帯域に分離する検出位置信号分離部と、
     前記センサ信号分離部で分離された前記加速度センサ信号と、前記検出位置信号分離部で分離された前記検出位置信号とを用いて、前記2つ以上の周波数帯域の各々において前記加速度センサ信号を較正し、前記2つ以上の周波数帯域の各々の運動軌跡成分を得るデータ較正部と、
     前記2つ以上の周波数帯域の各々の運動軌跡成分を結合し、運動軌跡として出力する運動軌跡算出部と
     を備えることを特徴とする機械運動軌跡測定装置。
  2.  アクチュエータを有し、被駆動体の運動軌跡が指令軌跡に追従するように、前記アクチュエータの位置または前記被駆動体の位置を検出する位置検出器から出力される検出位置信号をフィードバックし、前記アクチュエータを駆動する機械装置の運動軌跡を測定する機械運動軌跡測定装置であって、
     運動軌跡測定対象物の加速度を測定し、加速度センサ信号として出力する加速度センサと、
     前記運動軌跡測定対象物の軌跡測定の基準となる点の加速度を測定し、基準加速度センサ信号として出力する基準加速度センサと、
     前記加速度センサ信号と基準加速度センサ信号とを2つ以上の周波数帯域に分離するセンサ信号分離部と、
     前記検出位置信号を、前記センサ信号分離部と同じ周波数帯域に分離する検出位置信号分離部と、
     前記センサ信号分離部で分離された前記加速度センサ信号と、前記センサ信号分離部で分離された前記基準加速度センサ信号と、前記検出位置信号分離部で分離された前記検出位置信号とを用いて、前記2つ以上の周波数帯域の各々において前記加速度センサ信号を較正し、前記2つ以上の周波数帯域の各々の運動軌跡成分を得るデータ較正部と、
     前記2つ以上の周波数帯域の各々の運動軌跡成分を結合し、運動軌跡を示す情報として出力する運動軌跡算出部と
     を備えることを特徴とする機械運動軌跡測定装置。
  3.  前記センサ信号分離部は、前記加速度センサ信号を高周波帯域と低周波帯域の2つの帯域に分離し、周波数分離後の前記加速度センサ信号を出力することを特徴とする請求項1に記載の機械運動軌跡測定装置。
  4.  前記センサ信号分離部は、前記加速度センサ信号および前記基準加速度センサ信号の何れかを高周波帯域と低周波帯域の2つの帯域に分離し、周波数分離後の前記加速度センサ信号および前記基準加速度センサ信号の何れかを出力することを特徴とする請求項2に記載の機械運動軌跡測定装置。
  5.  前記センサ信号分離部は、前記高周波帯域、前記低周波帯域、およびノイズ帯域の3つの帯域に分離することを特徴とする請求項3に記載の機械運動軌跡測定装置。
  6.  前記低周波帯域は、前記アクチュエータを駆動する機械装置の制御帯域であり、前記高周波帯域は、前記アクチュエータを駆動する機械装置の振動周波数成分であり、前記ノイズ帯域は、前記加速度センサ信号のノイズ成分であることを特徴とする請求項5に記載の機械運動軌跡測定装置。
  7.  前記センサ信号分離部は、前記アクチュエータを駆動する機械装置が備える送り軸の位置制御系の応答を模擬した伝達関数を有することを特徴とする請求項1から請求項5の何れか一項に記載の機械運動軌跡測定装置。
  8.  前記センサ信号分離部は、前記低周波帯域用のフィルタとの伝達関数の和が1となるように設計された高周波帯域用のフィルタを有することを特徴とする請求項6に記載の機械運動軌跡測定装置。
  9.  前記データ較正部は、前記センサ信号分離部で分離された前記加速度センサ信号から周波数帯域毎に運動軌跡成分の計算誤差を補正するとき、前記検出位置信号分離部で分離された周波数帯域毎の前記検出位置信号を用いることを特徴とする請求項1に記載の機械運動軌跡測定装置。
  10.  前記データ較正部は、前記センサ信号分離部で分離された前記加速度センサ信号および前記基準加速度センサ信号から周波数帯域毎に対象物運動軌跡成分および基準運動軌跡成分の計算を行うとき、前記検出位置信号分離部で分離された周波数帯域毎の前記検出位置信号を用いて前記対象物運動軌跡成分および前記基準運動軌跡成分の計算誤差を補正し、前記対象物運動軌跡成分と前記基準運動軌跡成分との差をその帯域における相対運動軌跡として出力することを特徴とする請求項2に記載の機械運動軌跡測定装置。
  11.  前記データ較正部は、データ較正の際に加速度の入力信号を2階積分し、積分結果の平均値が0となるように積分誤差を補正することを特徴とする請求項9または請求項10に記載の機械運動軌跡測定装置。
  12.  前記データ較正部は、高周波帯域のデータ較正の際、前記積分結果として出力される信号に対して、ハイパスフィルタの通過帯域以下の周波成分を除去して出力することを特徴とする請求項11に記載の機械運動軌跡測定装置。
  13.  前記データ較正部は、高周波帯域のデータ較正の際、入力された信号を低周波数領域の信号に変換した後、低周波数領域で積分操作を行い、その後逆変換によって時間領域に変換した信号を運動軌跡の高周波帯域成分として出力することを特徴とする請求項9または請求項10に記載の機械運動軌跡測定装置。
  14.  前記データ較正部は、低周波帯域のデータ較正の際、前記センサ信号分離部で分離された加速度センサ信号の1階積分結果と前記検出位置信号の低周波帯域成分の1階微分結果との差が許容値を超えないように誤差を補正することを特徴とする請求項9または請求項10に記載の機械運動軌跡測定装置。
  15.  前記データ較正部は、低周波帯域のデータ較正の際、前記センサ信号分離部で分離された加速度センサ信号の2階積分結果と前記検出位置信号の低周波帯域成分との差が許容値を超えないように誤差を補正し、低周波帯域における対象物の運動軌跡を計算することを特徴とする請求項9、請求項10または請求項14に記載の機械運動軌跡測定装置。
  16.  前記機械装置の運動軌跡を表示する運動軌跡表示部を備え、
     前記運動軌跡表示部は、位置指令生成部に入力される目標位置から目標軌跡を合成し、前記位置指令生成部から出力される指令位置から指令軌跡を合成し、前記検出位置信号から検出部軌跡を合成し、前記目標軌跡、前記指令軌跡および前記検出部軌跡の何れか一の軌跡と前記機械装置の運動軌跡とを重ねて前記運動軌跡表示部に表示させることを特徴とする請求項1または請求項2に記載の機械運動軌跡測定装置。
PCT/JP2016/084562 2016-03-16 2016-11-22 機械運動軌跡測定装置 WO2017158932A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
DE112016006602.3T DE112016006602T5 (de) 2016-03-16 2016-11-22 Maschinenbewegungsbahnverlaufsmessvorrichtung
JP2017526604A JP6199003B1 (ja) 2016-03-16 2016-11-22 機械運動軌跡測定装置
CN201680024212.0A CN107533325B (zh) 2016-03-16 2016-11-22 机械运动轨迹测定装置
US15/564,264 US10543574B2 (en) 2016-03-16 2016-11-22 Machine motion trajectory measuring apparatus

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2016-053139 2016-03-16
JP2016053139 2016-03-16

Publications (1)

Publication Number Publication Date
WO2017158932A1 true WO2017158932A1 (ja) 2017-09-21

Family

ID=59850199

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2016/084562 WO2017158932A1 (ja) 2016-03-16 2016-11-22 機械運動軌跡測定装置

Country Status (5)

Country Link
US (1) US10543574B2 (ja)
JP (1) JP6199003B1 (ja)
CN (1) CN107533325B (ja)
DE (1) DE112016006602T5 (ja)
WO (1) WO2017158932A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020217597A1 (ja) * 2019-04-22 2020-10-29 三菱電機株式会社 サーボ制御装置

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6606054B2 (ja) * 2016-12-22 2019-11-13 Dmg森精機株式会社 工作機械の運動誤差同定方法
JP6538772B2 (ja) * 2017-07-26 2019-07-03 ファナック株式会社 数値制御装置
WO2019090487A1 (zh) * 2017-11-07 2019-05-16 大连理工大学 数控机床高动态大范围任意轮廓误差单目六维测量方法
US20190196422A1 (en) * 2017-12-22 2019-06-27 Industrial Technology Research Institute Tuning system, simulation unit and tuning and simulation method thereof
CN109000702A (zh) * 2018-05-16 2018-12-14 苏州汇川技术有限公司 编码器校正系统及方法
DE102018122759A1 (de) * 2018-09-17 2020-03-19 Chiron-Werke Gmbh & Co. Kg Verfahren und Vorrichtung für eine Verschleißanalyse an einer Werkzeugmaschine
CN108983699B (zh) * 2018-09-25 2023-08-15 辽宁工业大学 基于plc的三轴数控钻床控制系统及其控制方法
WO2020080181A1 (ja) * 2018-10-19 2020-04-23 Thk株式会社 アクチュエータ
JP7170510B2 (ja) * 2018-11-12 2022-11-14 オークマ株式会社 モニタリング装置及びモニタリング方法
US11198227B2 (en) 2018-12-04 2021-12-14 Raytheon Company Adjustable ballast system and method for same
US11027435B2 (en) * 2018-12-04 2021-06-08 Raytheon Company Automated work piece testing system and method for same
PL3864474T3 (pl) 2019-01-18 2024-03-11 Károly SZIPKA System pomiarowy i sposób powiązany z systemem pomiarowym
IT201900016778A1 (it) * 2019-09-19 2021-03-19 Springa S R L Sistema di localizzazione per macchina utensile e metodo di localizzazione realizzato da detto sistema
CN113205119B (zh) * 2021-04-20 2023-10-31 北京百度网讯科技有限公司 数据标注方法、装置、电子设备及可读存储介质
CN113129340B (zh) * 2021-06-15 2021-09-28 萱闱(北京)生物科技有限公司 操作设备的运动轨迹分析方法、装置、介质和计算设备
CN114193456B (zh) * 2022-01-04 2023-12-26 上海博远升控制技术有限公司 一种使用多峰输入整形法的抑振机器人及其控制方法
CN115047215B (zh) * 2022-02-18 2023-03-21 太原理工大学 三轴振动加速度信号耦合修正检测系统及方法
CN116871709B (zh) * 2023-09-08 2023-12-08 常州市鹏瑞自动化设备有限公司 激光切割异形圆管件的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006011002A (ja) * 2004-06-25 2006-01-12 Advanced Telecommunication Research Institute International 音声応答装置、音声応答方法および音声応答プログラム
WO2010067651A1 (ja) * 2008-12-09 2010-06-17 三菱電機株式会社 機械運動軌跡測定装置、数値制御工作機械および機械運動軌跡測定方法
JP2011115885A (ja) * 2009-12-02 2011-06-16 Mitsubishi Electric Corp 軌跡測定装置
JP2012196749A (ja) * 2011-03-23 2012-10-18 Seiko Epson Corp ロボット装置

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59225315A (ja) 1983-06-06 1984-12-18 Tokyo Keiki Co Ltd 座標位置検出装置
JPS60244859A (ja) * 1984-05-18 1985-12-04 Omron Tateisi Electronics Co 工具折損検出装置
KR900007293B1 (ko) * 1984-05-18 1990-10-08 다데이시 덴기 가부시기가이샤 공구 감시시스템
JPS617907A (ja) * 1984-06-22 1986-01-14 Sanesu Shoko:Kk 数値制御装置
JPS61209857A (ja) * 1985-03-08 1986-09-18 Yoshiaki Kakino Nc工作機械の運動精度試験方法および装置
JPS63315909A (ja) 1987-06-19 1988-12-23 Fujitsu Ltd 慣性航法装置
US5847960A (en) 1995-03-20 1998-12-08 Electro Scientific Industries, Inc. Multi-tool positioning system
DE10246093C1 (de) * 2002-10-02 2003-11-27 Siemens Ag Verfahren zur Dämpfung mechanischer Schwingungen von Achsen von Werkzeugmaschinen, Produktionsmaschinen oder Robotern
JP2006110702A (ja) * 2004-10-18 2006-04-27 Fanuc Ltd 学習制御機能を備えたロボット及びロボットの制御方法
CN101328535A (zh) * 2007-06-18 2008-12-24 北京翔博科技有限责任公司 振动消除应力多任务控制系统
JP4944806B2 (ja) * 2007-10-09 2012-06-06 オークマ株式会社 位置制御装置
JP4609735B2 (ja) * 2008-03-31 2011-01-12 ソニー株式会社 情報処理装置および方法、プログラム、並びに記録再生装置
JP5479855B2 (ja) 2009-11-10 2014-04-23 アピックヤマダ株式会社 切削装置及び切削方法
CN101967553A (zh) * 2010-10-12 2011-02-09 北京翔博科技有限责任公司 三维振动消除应力多任务采集和控制系统及其方法
CN101984379B (zh) * 2010-10-25 2012-03-28 武汉华中数控股份有限公司 一种基于指令序列的数控车削加工过程振动抑制方法
JP5374613B2 (ja) * 2012-05-30 2013-12-25 株式会社神戸製鋼所 多関節ロボットの弾性変形補償制御装置および制御方法
WO2014141515A1 (ja) * 2013-03-14 2014-09-18 三菱電機株式会社 サーボ制御装置
US9557799B2 (en) * 2013-07-02 2017-01-31 Infineon Technologies Ag Sensor interface systems and methods
CN103792888B (zh) * 2014-01-26 2016-12-07 西安交通大学苏州研究院 基于进给系统位置反馈信号的振动抑制控制器
JP6298672B2 (ja) * 2014-03-20 2018-03-20 Dmg森精機株式会社 運動誤差測定装置及びこれを備えた工作機械
CN104181863A (zh) * 2014-07-28 2014-12-03 上海理工大学 基于位置反馈信号的数控机床误差补偿系统
CN105511399B (zh) * 2015-12-02 2016-09-21 山东科技大学 一种结构优化的伺服电机速度闭环控制方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006011002A (ja) * 2004-06-25 2006-01-12 Advanced Telecommunication Research Institute International 音声応答装置、音声応答方法および音声応答プログラム
WO2010067651A1 (ja) * 2008-12-09 2010-06-17 三菱電機株式会社 機械運動軌跡測定装置、数値制御工作機械および機械運動軌跡測定方法
JP2011115885A (ja) * 2009-12-02 2011-06-16 Mitsubishi Electric Corp 軌跡測定装置
JP2012196749A (ja) * 2011-03-23 2012-10-18 Seiko Epson Corp ロボット装置

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020217597A1 (ja) * 2019-04-22 2020-10-29 三菱電機株式会社 サーボ制御装置
WO2020217282A1 (ja) * 2019-04-22 2020-10-29 三菱電機株式会社 サーボ制御装置
JPWO2020217597A1 (ja) * 2019-04-22 2021-10-28 三菱電機株式会社 サーボ制御装置
JP7034383B2 (ja) 2019-04-22 2022-03-11 三菱電機株式会社 サーボ制御装置

Also Published As

Publication number Publication date
US20180133860A1 (en) 2018-05-17
CN107533325A (zh) 2018-01-02
DE112016006602T5 (de) 2018-12-13
JP6199003B1 (ja) 2017-09-20
CN107533325B (zh) 2020-01-21
JPWO2017158932A1 (ja) 2018-03-22
US10543574B2 (en) 2020-01-28

Similar Documents

Publication Publication Date Title
JP6199003B1 (ja) 機械運動軌跡測定装置
US9144869B2 (en) Machine motion trajectory measuring device, numerically controlled machine tool, and machine motion trajectory measuring method
JP6316323B2 (ja) モータ制御装置
JP6275245B2 (ja) 送り軸の制御方法および数値制御工作機械
US9304504B2 (en) Servo controller for reducing interference between axes in machining
JP5388823B2 (ja) 軌跡測定装置
JPWO2012057235A1 (ja) 数値制御方法
WO2010150663A1 (ja) サーボ制御装置
JP2006216051A (ja) 未加工品と工作機械の工具との間の相対運動の運動分割方法
JP5249452B1 (ja) 補正データを考慮した軌跡表示装置
JP2019141967A (ja) 振動解析装置および振動解析方法
JP5972553B2 (ja) 位置決め制御装置、これを備えた工作機械
JP5404507B2 (ja) 補正パラメータ調整装置
US9921568B2 (en) Trajectory measuring device, numerical control device, and trajectory measuring method
CN113942008A (zh) 过冲量检测方法及机器人系统
EP3864474B1 (en) Measurement system, and a method in relation to the measurement system
WO2020217282A1 (ja) サーボ制御装置
WO2018066072A1 (ja) 制御装置
JP5225060B2 (ja) 機械運動測定装置
Hatwesh et al. Expedient identification of the dynamics of ball-screw drives using online-validation
TW201617753A (zh) 適應性凹陷濾波器減振控制裝置

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref document number: 2017526604

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 15564264

Country of ref document: US

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 16894549

Country of ref document: EP

Kind code of ref document: A1

122 Ep: pct application non-entry in european phase

Ref document number: 16894549

Country of ref document: EP

Kind code of ref document: A1