WO2013038856A1 - アクチュエータの制御方法及びアクチュエータの制御装置 - Google Patents

アクチュエータの制御方法及びアクチュエータの制御装置 Download PDF

Info

Publication number
WO2013038856A1
WO2013038856A1 PCT/JP2012/070480 JP2012070480W WO2013038856A1 WO 2013038856 A1 WO2013038856 A1 WO 2013038856A1 JP 2012070480 W JP2012070480 W JP 2012070480W WO 2013038856 A1 WO2013038856 A1 WO 2013038856A1
Authority
WO
WIPO (PCT)
Prior art keywords
control
time
output
maximum
actuator
Prior art date
Application number
PCT/JP2012/070480
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 EP12831914.2A priority Critical patent/EP2757427B1/en
Priority to US14/345,167 priority patent/US9519275B2/en
Priority to CN201280044822.9A priority patent/CN103797429B/zh
Publication of WO2013038856A1 publication Critical patent/WO2013038856A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00—Automatic controllers
    • G05B11/01—Automatic controllers electric
    • G05B11/14—Automatic controllers electric in which the output signal represents a discontinuous function of the deviation from the desired value, i.e. discontinuous controllers
    • G05B11/16—Two-step controllers, e.g. with on/off action
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00—Automatic controllers
    • G05B11/01—Automatic controllers electric
    • G05B11/36—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/0205—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system
    • G05B13/021—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system in which a variable is automatically adjusted to optimise the performance
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00—Programme-control systems
    • G05B19/02—Programme-control systems electric
    • G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/19—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by 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
    • G05B19/21—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by 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 using an incremental digital measuring device
    • G05B19/23—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by 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 using an incremental digital measuring device for point-to-point control
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00—Programme-control systems
    • G05B19/02—Programme-control systems electric
    • G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/19—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by 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
    • G05B19/27—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by 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 using an absolute digital measuring device
    • G05B19/29—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by 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 using an absolute digital measuring device for point-to-point control
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00—Programme-control systems
    • G05B19/02—Programme-control systems electric
    • G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/19—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by 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
    • G05B19/33—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by 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 using an analogue measuring device
    • G05B19/35—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by 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 using an analogue measuring device for point-to-point control
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00—Programme-control systems
    • G05B19/02—Programme-control systems electric
    • G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/416—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by control of velocity, acceleration or deceleration
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00—Program-control systems
    • G05B2219/30—Nc systems
    • G05B2219/41—Servomotor, servo controller till figures
    • G05B2219/41086—Bang bang control
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00—Program-control systems
    • G05B2219/30—Nc systems
    • G05B2219/41—Servomotor, servo controller till figures
    • G05B2219/41087—Determine switch point
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00—Program-control systems
    • G05B2219/30—Nc systems
    • G05B2219/41—Servomotor, servo controller till figures
    • G05B2219/41414—Time delay control, estimate non linear dynamics, correct with time delayed input
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00—Program-control systems
    • G05B2219/30—Nc systems
    • G05B2219/43—Speed, acceleration, deceleration control ADC
    • G05B2219/43034—Form of profile, ramp, trapezoid, S-curve, exponential
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00—Program-control systems
    • G05B2219/30—Nc systems
    • G05B2219/43—Speed, acceleration, deceleration control ADC
    • G05B2219/43062—Maximum acceleration, limit

Definitions

  • the present invention relates to an actuator control method and an actuator control device, and more particularly, in a shortest time control, a feedback element can be incorporated by using a control law that corrects an output pattern at every preset time.
  • the present invention relates to an actuator control method and an actuator control apparatus that can converge the control output after the control is completed.
  • the simplest shortest time control is a control in which the control target is accelerated with the maximum thrust in half of the stroke to the target position, and the remaining half of the stroke is decelerated with the maximum deceleration, and the control target is stopped at the target position. Since this output pattern is determined before the start of control, this shortest time control can be said to be feedforward control.
  • the shortest time control is a control method in which the actuator is moved with the maximum driving force of the actuator and stopped with the maximum braking force. Theoretically, the control object can be stopped at the target in the shortest time. In other words, the shortest time control is a control method that completely meets the purpose of the control of “stopping the controlled object at the target position in the shortest time”.
  • a control means for controlling a servo motor and a value at no load are used as a reference.
  • Correspondence relationship storage means for storing the relationship between the controlled amount and the load weight, load estimation calculation means, acceleration / deceleration constant determination means for determining an acceleration / deceleration constant based on the work information calculated by the load estimation calculation means, , Equipped with command creation means to create a command to pay out to the servo control means using the determined acceleration / deceleration constant, when the workpiece is gripped, the acceleration time becomes long, and when the workpiece is not gripped, A robot minimum time control device that shortens the acceleration time has been proposed.
  • the output pattern is an open control in which the initial speed, maximum acceleration, and maximum deceleration are determined in consideration of the feedback element. Therefore, when the target and the control amount do not match, there is no correction method and it is difficult to accurately match the target and the control amount. Therefore, it is rarely used for actual control. Met.
  • the present inventor uses the shortest time control and calculates for control using the maximum acceleration ⁇ p and the maximum deceleration ⁇ m at the maximum output of the measured control force of the actuator in advance.
  • the switching time t1 for switching from the acceleration output to the deceleration output and the end time t2 of the deceleration output are calculated, and the control force of the actuator is maximized from the calculation time t0 to the switching time t1.
  • the acceleration output is set to the maximum deceleration output from the switching time t1 to the end time t2, the output of the control force is ended at the end time t2, and the calculation step is performed every preset time.
  • an actuator control method and an actuator control device that calculate and update the switching time t1 and the end time t2 were considered.
  • This control method and control apparatus incorporates a feedback element that updates the switching time and the end time by adding a deviation between the target amount and the control amount at each time of control at a preset time. Even if there is a change, and even if the control time interval is not increased, a stable control result can always be obtained.
  • the present invention has been made in view of the above situation, and an object of the present invention is to incorporate an element of feedback control into the shortest time control and to converge a control output after the end of control and an actuator control method. It is to provide a control device.
  • the actuator control method of the present invention for achieving the above object uses the shortest time control and uses the maximum acceleration ⁇ p and the maximum deceleration ⁇ m at the maximum output of the measured control force of the actuator in advance.
  • the control force of the actuator is set to the maximum acceleration output until the switching time t1, the control force of the actuator is set to the maximum deceleration output from the switching time t1 to the end time t2, and the output of the control force is ended at the end time t2.
  • a control output step, and the calculation step is repeated at predetermined time intervals to change the switching time t1 and the end time.
  • a method comprising: an updating step for calculating and updating t2, and a control output decreasing step for reducing the control output together with a decrease in residual energy that is the sum of the remaining work and kinetic energy of the control system. It is.
  • the target trajectory of control from the position of the control target object to the target position at the calculation time t0 is recalculated in consideration of the deviation and speed every predetermined period or every irregular time, and the switching time Update t1 and end time t2.
  • the switching time t1 and the end time t2 are elapsed times from the recalculated calculation time t0.
  • the maximum acceleration ⁇ p or the maximum deceleration ⁇ m is the second-order differential value of the control value or the second-order difference value of the control value, and the control value can be applied to a temperature that does not involve movement. it can.
  • the control since the control is the shortest time control, the control speed becomes high, the maximum acceleration ⁇ p and the maximum deceleration ⁇ m can be set by measurement, and there is no adjustment term other than this, so it is necessary to adjust the control gain. No. In addition, since it is ON / OFF control, there is no need to output an intermediate output, and the controller and driver can be simplified.
  • the calculation is performed again with the deviation X between the target amount and the control amount at each control time, and a switching time t1 and an end time t2 that are elapsed times from the recalculation calculation time t0 are displayed. Since a feedback element of updating is incorporated, a stable control result can always be obtained even if the external force changes or the control time interval is not increased. As a result, it is possible to achieve both “control speed” and “stability”, both of which are contradictory issues in the conventional control law.
  • control output can be converged to zero after the control is completed.
  • each of the maximum acceleration output and the maximum deceleration output is multiplied by the residual energy, and this multiplied value is further applied to the residual energy.
  • the value obtained by multiplying the energy gain, which is a coefficient for determining the magnitude of the control output is the output of the control force, the control output can be easily converged to zero after the control is completed.
  • the initial speed is represented by a combination of two quadratic curves that contact each other from the position at the calculation time t0 to the target position at the end time t2.
  • V0 and the deviation between the target position and the control amount are X
  • the switching time t1 and the end time t2 are When calculating with the formulas (1) and (2), the switching time t1 and the end time t2 can be easily calculated.
  • ⁇ 1 is the maximum acceleration ⁇ p or maximum deceleration ⁇ m from calculation time t0 to switching time t1
  • ⁇ 2 is the maximum deceleration ⁇ m or maximum acceleration ⁇ p from switching time t1 to end time t2.
  • the actuator control apparatus of the present invention for achieving the above object uses the shortest time control, and the maximum acceleration ⁇ p and the maximum deceleration ⁇ m at the maximum output of the control force of the actuator measured in advance. And a calculation means for calculating a switching time t1 for switching from the acceleration output to the deceleration output, an end time t2 of the deceleration output, with the elapsed time display from the calculation time t0 when the calculation for the control is performed, and the calculation time t0 From the switching time t1 to the end time t1, the control force of the actuator is set to the maximum acceleration output, from the switching time t1 to the end time t2, the control force of the actuator is set to the maximum deceleration output, and the control force is output at the end time t2.
  • a control output means for ending is provided, and the switching time t1 and the end time t2 are repeated by the calculating means at every preset time.
  • control output reducing means multiplies each of the maximum acceleration output and the maximum deceleration output by the residual energy, and further multiplies the multiplied value with respect to the residual energy.
  • a value obtained by multiplying the energy gain, which is a coefficient for determining the magnitude of the control output, is configured as the output of the control force.
  • the update means recalculates a value obtained by multiplying the maximum acceleration and the maximum deceleration by a control output limit value as a new maximum acceleration and maximum deceleration. Is done.
  • the calculation means may represent the initial speed by expressing the trajectory from the position at the calculation time t0 to the target position at the end time t2 as a combination of two quadratic curves.
  • V0 and the deviation between the target position and the control amount are X
  • the switching time t1 and the end time t2 are (1) and (2).
  • ⁇ 1 is the maximum acceleration ⁇ p or maximum deceleration ⁇ m from calculation time t0 to switching time t1
  • ⁇ 2 is the maximum deceleration ⁇ m or maximum acceleration ⁇ p from switching time t1 to end time t2.
  • the acceleration ⁇ p between the calculation time t0 and the switching time t1 and the acceleration ⁇ m between the switching time t1 and the end time t2 are determined according to Table 1.
  • the above-described actuator control method can be implemented, and the same effects can be achieved.
  • the actuator control method and the actuator control apparatus of the present invention since the control is the shortest time control, the control speed is high, the maximum acceleration and the maximum deceleration can be set by measurement, and there is no adjustment term other than this. Therefore, there is no need to adjust the control gain. In addition, since it is ON / OFF control, there is no need to output an intermediate output, and the controller and driver can be simplified.
  • control output can be converged after the end of the control by reducing the control output along with the decrease of the residual energy that is the sum of the remaining work and the kinetic energy of the control system.
  • FIG. 1 is a diagram showing a shortest time control model for explaining the shortest time control used in the actuator control method according to the embodiment of the present invention.
  • FIG. 2 is a diagram for explaining the calculation conditions of the control trajectory.
  • FIG. 3 is a diagram for explaining a change in the control trajectory in the recalculation.
  • FIG. 4 is a diagram showing target trajectories that f (t) and g (t) can take.
  • FIG. 5 is a diagram illustrating an example of a control flow of the actuator control method according to the embodiment of the present invention.
  • FIG. 6 is a diagram showing a detailed control flow in step S18 of the control flow in FIG.
  • FIG. 7 is a diagram showing residual energy in the shortest time control in the actuator control method according to the embodiment of the present invention.
  • FIG. 1 is a diagram showing a shortest time control model for explaining the shortest time control used in the actuator control method according to the embodiment of the present invention.
  • FIG. 2 is a diagram for explaining the calculation
  • FIG. 8 is a diagram illustrating an example of a relationship between residual energy and control output (control output limit value) in the actuator control method according to the embodiment of the present invention.
  • FIG. 9 is a diagram showing another example of the relationship between the residual energy and the control output (control output limit value) in the actuator control method according to the embodiment of the present invention.
  • FIG. 10 is a diagram showing a control result when there is no damping force in the mechanical system of the feedback shortest time control in which the energy gain of the actuator control method according to the embodiment of the present invention is introduced.
  • FIG. 11 is a diagram illustrating a control result when there is a damping force of the mechanical system of the feedback shortest time control in which the energy gain of the actuator control method according to the embodiment of the present invention is introduced.
  • FIG. 10 is a diagram showing a control result when there is no damping force in the mechanical system of the feedback shortest time control in which the energy gain of the actuator control method according to the embodiment of the present invention is introduced.
  • FIG. 11 is a diagram
  • FIG. 12 is a diagram illustrating the control results of the example and the comparative example when the damping coefficient of the mechanical system is 1.
  • FIG. 13 is a diagram illustrating the control results of the example and the comparative example when the damping coefficient of the mechanical system is 0.1.
  • FIG. 14 is a diagram illustrating the control results of the example and the comparative example when the damping coefficient of the mechanical system is 2.
  • FIG. 15 is a diagram illustrating a simulation result of the example.
  • FIG. 16 is a diagram illustrating a simulation result of the comparative example.
  • FIG. 17 is a diagram illustrating a control result when there is no damping force of the mechanical system of the feedback shortest time control without introducing the energy gain of the actuator control method which is a reference of the present invention.
  • FIG. 18 is a diagram showing a control result when there is a damping force of the mechanical system of the feedback shortest time control without introducing the energy gain of the actuator control method which is a reference of the present invention.
  • FIG. 19 is a diagram illustrating a control result when there is no damping force of the mechanical system of the shortest time control according to the conventional technique.
  • FIG. 20 is a diagram showing a control result when there is a damping force of the mechanical system of the shortest time control according to the conventional technique.
  • FIG. 21 is a diagram illustrating a PID control model for explaining PID control.
  • FIG. 22 is a diagram illustrating the residual energy of PID control.
  • the actuator control apparatus uses the shortest time control, and includes a calculation unit, a control force output unit, an update unit, and a control output decrease unit.
  • This calculation means displays the elapsed time from the calculation time t0 when calculation for control is performed using the maximum acceleration ⁇ p and the maximum deceleration rate ⁇ m at the maximum output of the control force of the actuator measured in advance.
  • the switching time t1 for switching to the deceleration output and the deceleration output end time t2 are calculated.
  • the control output means sets the control force of the actuator to the maximum acceleration output from the calculation time t0 to the switching time t1, sets the control force of the actuator to the maximum deceleration output from the switching time t1 to the end time t2, and at the end time t2. Ends output of control force.
  • the updating means is configured to repeatedly calculate and update the switching time t1 and the end time t2 by the calculating means at a predetermined period or every predetermined time that is irregular.
  • the control output reducing means reduces the control output as the residual energy that is the sum of the remaining work and the kinetic energy of the control system decreases.
  • This calculation means represents the trajectory from the recalculated position at the calculation time t0 to the target position by a combination of two quadratic curves that are adjacent to each other, so that the initial speed is V0, and the deviation between the target amount and the control amount is X ,
  • the switching time t1 and the end time t2 are (1) and (2).
  • ⁇ 1 is the maximum acceleration ⁇ p or maximum deceleration ⁇ m from the calculation time t0 to the switching time t1
  • ⁇ 2 is the maximum deceleration ⁇ m or maximum acceleration ⁇ p from the switching time t1 to the end time t2.
  • the shortest time control used in the actuator control method according to the embodiment of the present invention is equivalent to a machine model in which the weight m shown in FIG.
  • the weight when the weight is lifted and released from the slope, the weight descends the slope with a constant acceleration due to gravity, and then rises on the opposite slope with a constant deceleration. If there is no energy loss such as friction, the weight m reaches zero at the same height as the beginning and stops for a moment. If gravity disappears at this time, the weight m will remain stationary at that position.
  • control object is stationary when In this way, if the position where the control object is stationary is the target position, the control ends.
  • the object to be controlled is accelerated / decelerated with the maximum thrust of the actuator and controlled to the target position, so that it can theoretically be controlled in the shortest time. Further, since the control output pattern at this time is determined before the start of control, the shortest time control is feedforward control.
  • the conventional PID control is based on the basic mass m, spring, and damper damping vibration model as shown in FIG. 21, and the role of the spring is the P term and the dash pod is the D term. Fulfill.
  • the zero point correction term is the I term, but here the I term has little physical meaning.
  • the shortest time control and the PID control are fundamentally different controls, but on the other hand, if the P gain and D gain of the PID control are greatly increased and the upper limit of the output is cut by the maximum output of the actuator, Since the control output waveform of the PID control approaches the control output waveform of the shortest time control, the shortest time control can be regarded as a control in which the PD gain of the PID control is increased to the limit.
  • the control is diverged if the PD gain is increased too much due to delays in calculation and control. This is because the PID control is feedback control, and the control output is always determined after the phenomenon. Therefore, if this delay becomes too large, the control diverges.
  • the shortest time control is feedforward control, and the control output is always determined prior to the phenomenon from the start of movement until it stops. Therefore, the control is stable, and the PD gain is increased by the shortest time control. Even if a control output equivalent to the PID control increased to the limit is output, the control does not diverge.
  • FIG. 19 and 20 show the simplest example of the shortest time control.
  • FIG. 19 shows a case where there is no mechanical damping force
  • FIG. 20 shows a case where there is a mechanical damping force.
  • the actuator is operated toward the target with the maximum thrust from the calculation time t0, the actuator is operated with the maximum deceleration at the switching time t1, and then the actuator thrust is set to zero at the end time t2. Yes.
  • the switching time t1 and the end time t2 are determined.
  • the switching time t1 and the end time t2 can be determined based on the following calculation conditions.
  • control trajectory is composed of two quadratic curves f (t) and g (t), and the control time t1 for switching the control output and the control based on the following conditions:
  • the end time t2 at which the output ends is obtained.
  • accelerations ⁇ 1 and ⁇ 2 are determined from Table 1.
  • V0 is a first-order differential value (or difference value) of the controlled variable obtained during the control.
  • the calculation conditions are the following (1) to (7).
  • the maximum acceleration ⁇ p and the maximum deceleration ⁇ m that can be generated at the maximum output of the actuator are known, that is, obtained from acceleration measured in advance.
  • the velocity V0 at the calculation time t0 is known, that is, obtained by first-order differentiation (or difference) of the measured value.
  • the value of the first quadratic curve f (t0) is zero at the calculation time t0.
  • the first-order differential value of the first quadratic curve f (t0) at the calculation time t0 is the speed (initial speed) at the time t0.
  • the second quadratic curve g (t) is in contact with the first quadratic curve f (t).
  • the value of the second quadratic curve g (t2) is the target value at the end time t2.
  • the first-order differential value of the second quadratic curve g (t) is zero at the end time t2.
  • A2 is a state in which when X> 0 and V0> 0, the target is not exceeded if the actuator is decelerated at the maximum deceleration ⁇ m generated.
  • the condition for not overshooting at time t3 is V0 2 / 2 ⁇ m ⁇ X.
  • A3 is a condition of X> 0, V0 ⁇ 0, and A4 is a condition of X ⁇ 0, V0> 0.
  • A5 is a state in which when X ⁇ 0, V0 ⁇ 0, the target will not be exceeded if the actuator is decelerated at the maximum deceleration ⁇ p generated.
  • the condition for not overshooting at time t3 is V0 2 / 2 ⁇ p> X.
  • A6 is a state in which, when X ⁇ 0, V0 ⁇ 0, the target is exceeded even if the actuator is decelerated at the maximum deceleration ⁇ p generated.
  • the condition for overshooting at time t3 is V0 2 / 2 ⁇ p ⁇ X.
  • the shortest time control used in the actuator control method according to the embodiment of the present invention can be performed by a control flow as shown in FIGS.
  • this control flow starts, data of maximum acceleration ⁇ p and maximum deceleration ⁇ m are read in step S11.
  • step S12 data of the target amount (target value) Tx and the control amount (control value) x is read.
  • counting of elapsed time t and elapsed time tc for recalculation is started.
  • step S13 it is determined whether or not the period of the orbit recalculation is reached, that is, whether or not the recalculated elapsed time tc after the orbit calculation is equal to or longer than the orbit recalculation period tcr.
  • step S13 If the trajectory recalculation period tcr is reached in step S13 (YES), the trajectory calculation is performed again in step S14, and then the process goes to step S15. In step S13, the start recalculation period is not reached. (NO), the trajectory calculation in step S14 is bypassed and the process goes to step S18.
  • the orbit recalculation period tcr is preferably about 1/10 of the control period. However, even if the orbit recalculation period tcr is the same as the control period, no problem occurs.
  • step S15 an acceleration limit value (control output limit value) for limiting the acceleration is calculated.
  • This acceleration limit value is a value obtained by multiplying the sum (residual energy) of the kinetic energy E and the absolute value
  • That is, “acceleration limit value (E +
  • step S17 Since the maximum value of the acceleration limit value is 1, it is determined in step S17 whether or not the acceleration limit value is greater than 1. If it is greater than 1 (YES), the acceleration limit value is set to 1 in step S17. After setting, go to step S18. If the acceleration limit value is not greater than 1 in step S17 (NO), leave it as is, and go to step S18.
  • step S18 when entering step S18, in the first step S18a, it is determined whether or not the elapsed time t is smaller than the switching time t1, and if it is smaller (YES), the process goes to step S18b. Go and set the output acceleration to “ ⁇ 1 ⁇ acceleration limit value”, then go to step S19.
  • step S18a if the elapsed time t is not smaller than the switching time t1 (NO), go to step S18c to determine whether the elapsed time t is smaller than the end time t2, and if smaller (YES), Go to step S18d, set the output acceleration to “ ⁇ 2 ⁇ acceleration limit value”, and then go to step S19.
  • step S18c if the elapsed time t is not smaller than the end time t2 (NO), the process goes to step S18e, the output acceleration is set to “zero”, and the process goes to step S19.
  • step S19 an actuator thrust corresponding to the output acceleration is generated for a preset time (time related to various determination intervals), and the control target is controlled. Also, the elapsed time t and the recalculated elapsed time tc are counted. Thereafter, the process returns to step S12, and steps S12 to S19 are repeated.
  • the thrust of the actuator can be controlled while the recalculated elapsed time tc recalculates the trajectory calculation in step S14 for each trajectory recalculation period tcr.
  • the target trajectory is corrected at predetermined intervals or at predetermined intervals.
  • FIG. 3 shows how the control trajectory changes due to this recalculation.
  • the control trajectory (dotted line) of the first calculation is actually scheduled to be X1
  • a new control trajectory (solid line) is calculated by recalculation, and the switching time t1 and The end time t2 also becomes a new value, and the control force is controlled based on the new switching time t1 and end time t2.
  • FIG. 17 and FIG. 18 show the control results when the switching time t1 and the end time t2 are recalculated at regular intervals. Recalculation is performed at the rising edge of the signal pulse.
  • a feedback (FB) shortest time control method In order to distinguish this shortest time control method from the shortest time control method of the prior art, here, it is referred to as a feedback (FB) shortest time control method.
  • FIG. 17 shows a case where there is no mechanical damping force
  • FIG. 18 shows a case where there is a mechanical damping force.
  • the control target can be matched with the target position with a typical trajectory.
  • an energy gain (E gain) is used as a coefficient for determining the magnitude of the control output with respect to the residual energy.
  • control output becomes zero at the end of the control in order to reduce the stability and energy consumption of the control. Therefore, by multiplying the control output by an appropriate residual energy as an evaluation value for the end of control, the control output is made zero simultaneously with the end of control.
  • the control output limit value is determined according to the map data obtained by converting the relationship of “residual energy ⁇ E gain and control output limit value” as shown in FIGS. 8 and 9, and this value is set to the maximum thrust.
  • the multiplied value is used as the control output.
  • the maximum value of this control output limit value is 1. Further, when recalculation is performed at predetermined intervals or at predetermined intervals, a value obtained by multiplying the maximum acceleration and maximum deceleration by the control output limit value is recalculated as the maximum acceleration and maximum deceleration.
  • the residual energy in PID control can be used as an evaluation index for control used in the evaluation of control used for quality engineering. It can be defined as the sum of potential energy and kinetic energy.
  • FIG. 12 to 14 show the comparison results between the feedback shortest time control example and the PID control comparison example.
  • the control gain was adjusted with reference to the case where the attenuation coefficient of the mechanical system shown in FIG. 13 and 14, only the damping coefficient of the mechanical system is changed to 0.1 and 2 while the control gain is fixed.
  • FIG. 12 shows almost the same result, although the comparative example of PID control is slightly faster.
  • FIG. 13 in which the attenuation coefficient is reduced, in the embodiment of the FB shortest time control, a control result without overshoot is obtained, whereas in the comparative example of PID control, a large overshoot occurs. .
  • FIG. 13 in which the attenuation coefficient is reduced
  • FIG. 13 in which the attenuation coefficient is reduced
  • FB shortest time control a control result without overshoot is obtained
  • a large overshoot occurs.
  • FIG. 13 in which the attenuation coefficient is reduced
  • the example of the feedback shortest time control quickly reaches the target, whereas the comparative example of PID control has a speed greatly reduced immediately before the target, and the response to the target is It is late.
  • the example of the feedback shortest time control it can be seen that a more stable control result can be obtained than the comparative example of the PID control with respect to the resistance change represented by the change of the attenuation coefficient of the controlled object.
  • FIG. 15 shows a simulation result of an example of feedback shortest time control
  • FIG. 16 shows a simulation result of a comparative example of PID control.
  • the control result greatly changes with respect to the load fluctuation of the control system in the comparative example of the PID control shown in FIG. 16, whereas in the embodiment of the feedback shortest time control shown in FIG. There is no disturbance in the result, and a stable control result is always obtained. That is, in the comparative example of FIG. 16, the control result changes greatly with respect to the load fluctuation of the control system, whereas in the embodiment of FIG. 15, the control result is not disturbed and a stable control result is always obtained. ing.
  • the actuator control method and the actuator control device described above since the control is performed for the shortest time, the control speed becomes high, the maximum acceleration and the maximum deceleration can be set by measurement, and there is no adjustment term other than this. Therefore, there is no need to adjust the control gain. In addition, since it is ON / OFF control, there is no need to output an intermediate output, and the controller and driver can be simplified.
  • control output can be converged to zero after the control is completed.
  • the control speed is high, and it is not necessary to adjust the control gain, and it is ON / OFF control, so there is no need to output an intermediate output.
  • the controller and driver can be simplified, and a feedback element is used to update the switching time and end time by adding a deviation between the target amount and control amount at each control time for each preset time. Even if the external force changes and the control cycle is not set to high speed, a stable control result can always be obtained, and the residual energy, which is the sum of the remaining work and kinetic energy of the control system, can be reduced.
  • the control output can be converged to zero after the end of the control.
  • Vessel such as it can be utilized as a control device for controlling the method and the actuator of a number of actuators.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Feedback Control In General (AREA)
  • Control Of Position Or Direction (AREA)

Abstract

アクチュエータの制御方法及び制御装置において、最短時間制御を用いると共に、制御力の最大出力時の最大加速度（αｐ）と最大減速度（αｍ）を用いて、加速出力から減速出力へ切り替える切替時刻（ｔ１）と、減速出力の終了時刻（ｔ２）を算出し、計算時刻（ｔ０）から切替時刻（ｔ１）までは、アクチュエータの制御力を最大加速出力とし、切替時刻（ｔ１）から終了時刻（ｔ２）まではアクチュエータの制御力を最大減速力とし、終了時刻（ｔ２）では制御力の出力を終了すると共に、予め設定した時間毎に切替時刻（ｔ１）と終了時刻（ｔ２）を繰り返して算出して更新し、更に、制御系が持つ残り仕事と運動エネルギーの和である残留エネルギーの減少と共に、制御出力を小さくする。これにより、最短時間制御にフィードバック制御の要素を取り入れると共に、制御終了後に制御出力を収束させる。

Description

アクチュエータの制御方法及びアクチュエータの制御装置
　本発明は、アクチュエータの制御方法及びアクチュエータの制御装置に関し、より詳細には、最短時間制御において、予め設定した時間毎に出力パターンを修正する制御則を用いることでフィードバックの要素を取り入れることができると共に、制御終了後に制御出力を収束させることができるアクチュエータの制御方法及びアクチュエータの制御装置に関する。
　従来の制御ではＰＩＤ制御のフィードバック制御が一般的に用いられてきている。このＰＩＤ制御では、現象に対して必ず遅れて制御出力が決定されるため、制御速度を上げようとして、ＰＩＤの各制御ゲインを大きくすると制御が現象に対し間に合わなくなり、制御が不安定となる。特に制御対象の機械的減衰力が著しく低下すると制御が不安定になり易く、制御が発散する場合もある。この制御の不安定を避けるためのＰＩＤ制御の各制御ゲインを決定する方法として、制御の安定性を保証できるＨ∞等の制御理論を適用している。しかし、ＰＩＤ制御という制約のもとでは、負荷変動によりオーバーシュートや制御遅れが生じる。
　ＰＩＤ制御においても、スライディングモード制御を使用すると、制御ゲインを制御状態によって切り替えることで、負荷変動の影響を理論的には排除できる。しかし、制御周期が遅くなると、この制御は振動したまま収束しなくなる。そのため、負荷変動の影響を完全に排除するためには、無限に速く制御ゲインを切り替える必要があり、現象に対し無限に速いと言える程度まで、高速なコントロールが必要になる。更に、ＰＩＤ等の各制御ゲインの調整が必要であり、この制御ゲインの調整の善し悪しが制御の善し悪しを決定してしまうため、この制御ゲインの調整が非常に重要な要素となってしまっている。
　また、これらの制御理論は、ＰＩＤ制御の欠点を補うもので、「最短の時間で目標位置に制御対象を静止させる」ことを制御の目的とするために構築された手法ではない。そのため、この単純な目的に対してはＰＩＤ制御よりも、むしろ、最短時間制御がより目的に合致した制御方法であるといえる。
　最も単純な最短時間制御は、目標位置までの行程の半分で制御対象を最大推力で加速し、残り半分の行程を最大減速度で減速して目標位置に制御対象を静止させる制御である。この出力パターンは制御開始前に決定されるので、この最短時間制御は、フィードフォワード制御と言える。
　言い換えれば、最短時間制御とはアクチュエータの最大駆動力で動かし、最大制動力で停止させる制御方法であり、理論的には最短時間で制御対象を目標に静止させることが可能な制御である。つまり、最短時間制御は「最短の時間で目標位置に制御対象を静止させる」という制御の目的に完全に合致する制御方法である。
　例えば、日本出願の特開２０００－９４３７１号公報に記載されているように、この最短時間制御を用いた制御装置として、サーボモータを制御するための制御手段と、無負荷時の値を基準とした制御量と負荷重量の関係を記憶しておく対応関係記憶手段と、負荷推定計算手段と、負荷推定計算手段で計算されたワーク情報を基に加減速定数を決定する加減速定数決定手段と、決定された加減速定数を使用してサーボ制御手段に払い出す指令を作成する指令作成手段を備えて、ワークを把持した時は、加速時間が長くなり、ワークを把持していない時は、加速時間を短くするロボットの最短時間制御装置が提案されている。
　しかしながら、最短時間制御は、理論的に最短の時間で制御可能な理想的な制御である一方、出力パターンが初速や最大加速度、最大減速度を考慮して決定されるオープン制御であり、フィードバック要素が無いため、目標と制御量が合致しないときに、修正の方法がなく、正確に目標と制御量を一致させることが困難であるという問題があるため、実際の制御に採用されることは稀であった。
日本出願の特開２０００－９４３７１号公報
　これに対して、本発明者は、最短時間制御を用いると共に、予め、計測されたアクチュエータの制御力の最大出力時の最大加速度αｐと最大減速度αｍを用いて、制御のための計算を行う計算時刻ｔ０からの経過時間表示で、加速出力から減速出力へ切り替える切替時刻ｔ１と、減速出力の終了時刻ｔ２を算出し、前記計算時刻ｔ０から前記切替時刻ｔ１までは前記アクチュエータの制御力を最大加速出力とし、前記切替時刻ｔ１から前記終了時刻ｔ２までは前記アクチュエータの制御力を最大減速出力とし、前記終了時刻ｔ２で制御力の出力を終了すると共に、前記算出ステップを予め設定した時間毎に繰り返し、前記切替時刻ｔ１と前記終了時刻ｔ２を算出して更新するアクチュエータの制御方法とアクチュエータの制御装置を考えた。
　この制御方法と制御装置は、予め設定した時間毎に、制御の各時刻における目標量と制御量の偏差を入れて切替時刻と終了時刻を更新するというフィードバックの要素を取り入れているため、外力が変化しても、また、制御の時間間隔を高速にしなくても、常に安定した制御結果を得ることができるという効果を奏することができた。
　しかしながら、この制御方法と制御装置においては、制御量が制御目標と合致してもパルス状の制御出力が出たままとなる問題が発生した。
　本発明は、上記の状況を鑑みてなされたものであり、その目的は、最短時間制御にフィードバック制御の要素を取り入れると共に、制御終了後に制御出力を収束させることができるアクチュエータの制御方法及びアクチュエータの制御装置を提供することにある。
　上記のような目的を達成するための本発明のアクチュエータの制御方法は、最短時間制御を用いると共に、予め、計測されたアクチュエータの制御力の最大出力時の最大加速度αｐと最大減速度αｍを用いて、制御のための計算を行う計算時刻ｔ０からの経過時間表示で、加速出力から減速出力へ切り替える切替時刻ｔ１と、減速出力の終了時刻ｔ２を算出する算出ステップと、前記計算時刻ｔ０から前記切替時刻ｔ１までは前記アクチュエータの制御力を最大加速出力とし、前記切替時刻ｔ１から前記終了時刻ｔ２までは前記アクチュエータの制御力を最大減速出力とし、前記終了時刻ｔ２で制御力の出力を終了する制御出力ステップを備えると共に、前記算出ステップを予め設定した時間毎に繰り返し、前記切替時刻ｔ１と前記終了時刻ｔ２を算出して更新する更新ステップと、更に、制御系が持つ残り仕事と運動エネルギーの和である残留エネルギーの減少と共に、制御出力を小さくする制御出力減少ステップを備えたことを特徴とする方法である。
　つまり、計算時刻ｔ０の制御対象物の位置から目標位置までの制御の目標軌跡を予め設定した時間毎である一定周期毎又は不定期毎に偏差、速度を考慮して再計算して、切替時刻ｔ１と終了時刻ｔ２を更新する。この切替時刻ｔ１と終了時刻ｔ２は再計算した計算時刻ｔ０のからの経過時間となる。
　なお、この最大加速度αｐ、又は、最大減速度αｍは、制御値の２階の時間微分値若しくは制御値の２階の差分値であり、制御値が運動を伴わない温度のようなものでも適用できる。
　この方法によれば、最短時間制御であるために、制御速度が高速となり、最大加速度αｐと最大減速度αｍは計測によって設定でき、これ以外に調整項が無いため、制御ゲインを調整する必要が無い。また、ＯＮ／ＯＦＦ制御であるため、中間的な出力を出す必要がなく、コントローラ、ドライバーを簡略化できる。
　また、予め設定した時間毎に、制御の各時刻における目標量と制御量の偏差Ｘを入れて計算し直して再計算の計算時刻ｔ０からの経過時間表示である切替時刻ｔ１と終了時刻ｔ２を更新するというフィードバックの要素を取り入れているため、外力が変化しても、また、制御の時間間隔を高速にしなくても、常に安定した制御結果を得ることができる。その結果、従来の制御則では相反する大きな課題であった「制御の速度」と「安定性」を両立させることができる。
　その上、制御系が持つ残り仕事と運動エネルギーの和である残留エネルギーの減少と共に、制御出力を小さくすることにより、制御終了後に制御出力を収束させて、ゼロにすることができる。
　また、上記のアクチュエータの制御方法において、前記制御出力減少ステップで、前記最大加速出力及び前記最大減速出力のそれぞれに、前記残留エネルギーを乗じ、この乗じた値に、更に、前記残留エネルギーに対して制御出力の大きさを決定するための係数であるエネルギーゲインを乗じた値を制御力の出力とすると、容易に制御終了後に制御出力を収束させて、ゼロにすることができる。
　また、上記のアクチュエータの制御方法において、前記更新ステップで、前記最大加速度、前記最大減速度に、制御出力制限値を掛け合わせた値を新たな最大加速度、最大減速度として再計算すると、容易に制御終了後に制御出力を収束させて、ゼロにすることができる。
　また、上記のアクチュエータの制御方法において、前記算出ステップにおいて、前記計算時刻ｔ０での位置から前記終了時刻ｔ２の目標位置までの軌道を相接する２つの２次曲線の組み合わせで表して、初速をＶ０、目標位置と制御量の偏差をＸとした場合に、前記切替時刻ｔ１と、前記終了時刻ｔ２を、
Figure JPOXMLDOC01-appb-M000006
Figure JPOXMLDOC01-appb-M000007
 　の（１）式と（２）式で算出すると、容易に、切替時刻ｔ１と記終了時刻ｔ２を算出できる。なお、α１は計算時刻ｔ０～切替時刻ｔ１の最大加速度αｐ又は最大減速度αｍであり、α２は切替時刻ｔ１～終了時刻ｔ２の最大減速度αｍ又は最大加速度αｐである。　　　
　そして、上記のような目的を達成するための本発明のアクチュエータの制御装置は、最短時間制御を用いると共に、予め、計測されたアクチュエータの制御力の最大出力時の最大加速度αｐと最大減速度αｍを用いて、制御のための計算を行う計算時刻ｔ０からの経過時間表示で、加速出力から減速出力へ切り替える切替時刻ｔ１と、減速出力の終了時刻ｔ２を算出する算出手段と、前記計算時刻ｔ０から前記切替時刻ｔ１までは前記アクチュエータの制御力を最大加速出力とし、前記切替時刻ｔ１から前記終了時刻ｔ２までは前記アクチュエータの制御力を最大減速出力とし、前記終了時刻ｔ２で制御力の出力を終了する制御出力手段を備えると共に、予め設定した時間毎に前記算出手段により前記切替時刻ｔ１と前記終了時刻ｔ２を繰り返して算出し、更新する更新手段と、更に、制御系が持つ残り仕事と運動エネルギーの和である残留エネルギーの減少と共に、制御出力を小さくする制御出力減少手段を備えて構成される。
　また、上記のアクチュエータの制御装置において、前記制御出力減少手段が、前記最大加速出力及び前記最大減速出力のそれぞれに、前記残留エネルギーを乗じ、この乗じた値に、更に、前記残留エネルギーに対して制御出力の大きさを決定するための係数であるエネルギーゲインを乗じた値を制御力の出力とするように構成される。
　また、上記のアクチュエータの制御装置において、前記更新手段が、前記最大加速度、前記最大減速度に、制御出力制限値を掛け合わせた値を新たな最大加速度、最大減速度として再計算するように構成される。
　また、上記のアクチュエータの制御装置において、前記算出手段が、前記計算時刻ｔ０での位置から前記終了時刻ｔ２の目標位置までの軌道を相接する２つの２次曲線の組み合わせで表して、初速をＶ０、目標位置と制御量の偏差をＸとした場合に、前記切替時刻ｔ１と、前記終了時刻ｔ２を、
Figure JPOXMLDOC01-appb-M000008
Figure JPOXMLDOC01-appb-M000009
 　の（１）式と（２）式で算出するように構成する。なお、α１は計算時刻ｔ０～切替時刻ｔ１の最大加速度αｐ又は最大減速度αｍであり、α２は切替時刻ｔ１～終了時刻ｔ２の最大減速度αｍ又は最大加速度αｐである。
　また、上記のアクチュエータの制御装置において、計算時刻ｔ０から切替時刻ｔ１の間の加速度αｐと切替時刻ｔ１から終了時刻ｔ２の間の加速度αｍを表１によって決定する。
Figure JPOXMLDOC01-appb-T000010
　これらの構成のアクチュエータの制御装置によれば、上記のアクチュエータの制御方法を実施でき、同様の効果を奏することができる。
　本発明に係るアクチュエータの制御方法及びアクチュエータの制御装置によれば、最短時間制御であるために、制御速度が高速となり、最大加速度と最大減速度は計測によって設定でき、これ以外に調整項はないため、制御ゲインを調整する必要が無い。また、ＯＮ／ＯＦＦ制御であるため、中間的な出力を出す必要がなく、コントローラ、ドライバーを簡略化できる。
　また、予め設定した時間毎に、制御の各時刻における目標量と制御量の偏差を入れて切替時刻と終了時刻を更新するというフィードバックの要素を取り入れているため、外力が変化しても、また、制御の時間間隔を高速にしなくても、常に安定した制御結果を得ることができる。その結果、従来の制御則では相反する大きな課題であった「制御の速度」と「安定性」を両立させることができる。
　その上、制御系が持つ残り仕事と運動エネルギーの和である残留エネルギーの減少と共に、制御出力を小さくすることにより、制御終了後に制御出力を収束させることができる。
図１は、本発明の実施の形態のアクチュエータの制御方法で用いられる最短時間制御を説明するための最短時間制御モデルを示す図である。 図２は、制御軌跡の計算条件を説明するための図である。 図３は、再計算における制御軌跡の変化を説明するための図である。 図４は、ｆ（ｔ）、ｇ（ｔ）の取り得る目標軌跡を示す図である。 図５は、本発明の実施の形態のアクチュエータの制御方法の制御フローの一例を示す図である。 図６は、図５の制御フローのステップＳ１８の詳細な制御フローを示す図である。 図７は、本発明の実施の形態のアクチュエータの制御方法における、最短時間制御における残留エネルギーを示す図である。 図８は、本発明の実施の形態のアクチュエータの制御方法における残留エネルギーと制御出力（制御出力制限値）の関係の一例を示す図である。 図９は、本発明の実施の形態のアクチュエータの制御方法における残留エネルギーと制御出力（制御出力制限値）の関係の他の例を示す図である。 図１０は、本発明の実施の形態のアクチュエータの制御方法のエネルギーゲインを導入したフィードバック最短時間制御の機械系の減衰力が無い場合の制御結果を示す図である。 図１１は、本発明の実施の形態のアクチュエータの制御方法のエネルギーゲインを導入したフィードバック最短時間制御の機械系の減衰力が有る場合の制御結果を示す図である。 図１２は、機械系の減衰係数が１の場合の実施例と比較例の制御結果を示す図である。 図１３は、機械系の減衰係数が０．１の場合の実施例と比較例の制御結果を示す図である。 図１４は、機械系の減衰係数が２の場合の実施例と比較例の制御結果を示す図である。 図１５は、実施例のシミュレーション結果を示す図である。 図１６は、比較例のシミュレーション結果を示す図である。 図１７は、本発明の参考となるアクチュエータの制御方法のエネルギーゲインを導入しないフィードバック最短時間制御の機械系の減衰力が無い場合の制御結果を示す図である。 図１８は、本発明の参考となるアクチュエータの制御方法のエネルギーゲインを導入しないフィードバック最短時間制御の機械系の減衰力が有る場合の制御結果を示す図である。 図１９は、従来技術の最短時間制御の機械系の減衰力が無い場合の制御結果を示す図である。 図２０は、従来技術の最短時間制御の機械系の減衰力が有る場合の制御結果を示す図である。 図２１は、ＰＩＤ制御を説明するためのＰＩＤ制御モデルを示す図である。 図２２は、ＰＩＤ制御の残留エネルギーを示す図である。
　以下、本発明に係る実施の形態のアクチュエータの制御方法及びアクチュエータの制御装置について、図面を参照しながら説明する。なお、ここでは、本発明の最短時間制御を明確にするために、ＰＩＤ制御との比較を行いながら説明する。
　本発明に係る実施の形態のアクチュエータの制御装置は、最短時間制御を用いると共に、算出手段と、制御力出力手段と、更新手段と、制御出力減少手段を備えて構成される。
　この算出手段は、予め計測されたアクチュエータの制御力の最大出力時の最大加速度αｐと最大減速度αｍを用いて、制御のための計算を行う計算時刻ｔ０からの経過時間表示で、加速出力から減速出力へ切り替える切替時刻ｔ１と、減速出力の終了時刻ｔ２を算出する。
　また、制御出力手段は、計算時刻ｔ０から切替時刻ｔ１までは、アクチュエータの制御力を最大加速出力とし、切替時刻ｔ１から終了時刻ｔ２まではアクチュエータの制御力を最大減速出力とし、終了時刻ｔ２で制御力の出力を終了する。
　更に、更新手段は、一定周期又は不定期な予め設定した時間毎に算出手段で切替時刻ｔ１と終了時刻ｔ２を繰り返し算出して更新するように構成される。また、制御出力減少手段は、制御系が持つ残り仕事と運動エネルギーの和である残留エネルギーの減少と共に、制御出力を小さくする。
　この算出手段は、再計算される計算時刻ｔ０の位置から目標位置までの軌道を相接する２つの２次曲線の組み合わせで表すことにして、初速をＶ０、目標量と制御量の偏差をＸとした場合に、切替時刻ｔ１と、終了時刻ｔ２を、
Figure JPOXMLDOC01-appb-M000011
Figure JPOXMLDOC01-appb-M000012
 　の（１）式と（２）式で算出する。ここで、α１は計算時刻ｔ０～切替時刻ｔ１の最大加速度αｐ又は最大減速度αｍであり、α２は切替時刻ｔ１～終了時刻ｔ２の最大減速度αｍ又は最大加速度αｐである。
　また、本発明の実施の形態のアクチュエータの制御方法で用いる最短時間制御は、図１に示す錘ｍが斜面を下って上がるような機械モデルと等価となる。図１のように、斜面に対し錘を持ち上げてから放すと、重力のために一定の加速度で錘は斜面を下り、続いて反対側の斜面を一定の減速度で上がる。摩擦等のエネルギー損失がなければ、錘ｍは始めと同じ高さまで登った所で速度ゼロとなり、一瞬停止する。この時に重力が無くなれば錘ｍはその位置に静止し続ける。
　これをアクチュエータで制御対象を移動させる制御システムに置き換えると、アクチュエータの最大推力で錘を押し、続いてアクチュエータの最大推力で錘を引き戻すと、推す時に入力した仕事量と、引き戻す時に入力した仕事量が釣り合った時に制御対象が静止する。この様にして制御対象が静止した位置が目標位置となっていれば制御終了となる。
　つまり、最短時間制御ではアクチュエータの最大推力で制御対象を加減速させ、目標位置に制御するため、理論的に最短の時間で制御できる。また、この時の制御出力パターンは、制御開始前に決定されるため、最短時間制御はフィードフォワード制御である。
　これに対し、従来技術のＰＩＤ制御は、図２１に示すような基本的な質量ｍ、ばね、ダンパ系の減衰振動モデルを基本としたものとなり、ばねがＰ項、ダッシュポッドがＤ項の役割を果たす。なお、ゼロ点の補正項がＩ項となるが、ここではＩ項には物理的な意味はあまりない。
　この制御をエネルギー変換で考えると、図２１の中段に示すように、目標位置を入力することで、目標変位としてばねを撓ませて、ばねの歪エネルギーを制御系に入力すると、ばね力に引かれて錘（質量ｍ）が動き始める。この時、歪エネルギーから運動エネルギーへの変換が行われる。錘ｍが運動を始めるとダッシュポッドにより運動エネルギーが熱エネルギーに変換され、系の外へエネルギーが放出される。初めに入力されたばねの歪エネルギーが全て熱エネルギーに変換されると、図２１の下段に示すように、制御対象である錘ｍは目標位置で静止することになる。なお、ＰＩＤ制御では制御対象の運動状態から制御出力が決定されるため、フィードバック制御である。
　このように最短時間制御とＰＩＤ制御は根本的に異なる制御であると言えるが、一方で、ＰＩＤ制御のＰゲイン、Ｄゲインを非常に大きくし、アクチュエータの最大出力で出力の上限をカットすると、ＰＩＤ制御の制御出力波形は最短時間制御の制御出力波形に近づいていくため、最短時間制御はＰＩＤ制御のＰＤゲインを極限まで大きくした制御であると見ることもできる。しかし、ＰＩＤ制御では、通常は計算や制御の遅れにより、ＰＤゲインをあまり大きくすると制御が発散してしまうことになる。これはＰＩＤ制御がフィードバック制御であり常に現象に遅れて制御出力が決定されるため、この遅れが大きくなりすぎると制御が発散するからである。
　これに対し、最短時間制御はフィードフォワード制御であり、動かし始めてから止めるまでを考慮した制御出力が常に現象よりも先に決定されるため、制御が安定であり、最短時間制御で、ＰＤゲインを極限まで大きくしたＰＩＤ制御と同等の制御出力を出しても制御は発散しない。
　この最短時間制御の最も単純な例を図１９及び図２０に示す。図１９は、機械系の減衰力が無い場合を、図２０は機械系の減衰力がある場合を示す。図１９及び図２０では、計算時刻ｔ０から最大推力で目標に向かってアクチュエータを作動させ、切替時刻ｔ１で今度は最大減速度でアクチュエータを作動させ、続いて終了時刻ｔ２でアクチュエータ推力をゼロにしている。このように最短時間制御では切替時刻ｔ１と終了時刻ｔ２を決定すれば制御が可能である。切替時刻ｔ１と終了時刻ｔ２は以下の計算条件を元に決定することができる。
　この計算条件では、制御対象が相接する２つの２次曲線上を通って目標に到達すると仮定する。そして、アクチュエータが発生可能な最大加速度αｐと、アクチュエータが発生可能な最大減速度αｍと、計算時刻ｔ０での目標量と制御量の偏差Ｘ（＝目標量－制御量：目標位置と制御時の位置の差）と、計算時刻ｔ０での制御対象の速度Ｖ０を既知の値とし、切替時刻ｔ１で２つの２次曲線が接することと、終了時刻ｔ２で制御対象の速度をゼロとすること（Ｖ＝０）、終了時刻ｔ２で偏差Ｘをゼロとすること（Ｘ＝０）を拘束条件として、切替時刻ｔ１と終了時刻ｔ２を求める。
　制御軌跡は、図２に示すように、２つの２次曲線ｆ（ｔ）、ｇ（ｔ）で構成されるものと仮定し、以下の条件を元に、制御出力を切り替える切替時刻ｔ１と制御出力を終了する終了時刻ｔ２求める。また加速度α１、α２を表１より決定する。また、Ｖ０は制御中に得られた制御量の１階の微分値（または差分値）である。
　計算条件は、次の（１）から（７）である。
（１）アクチュエータの最大出力時に発生可能な最大加速度αｐ、最大減速度αｍは既知、すなわち、予め測定した加速度より得られる。
（２）計算時刻ｔ０における速度Ｖ０は既知、すなわち、測定値の１階の微分（または差分）により得られる。
（３）計算時刻ｔ０で第１の２次曲線ｆ（ｔ０）の値がゼロである。
（４）計算時刻ｔ０で第１の２次曲線ｆ（ｔ０）の１階の微分値が時刻ｔ０での速度（初速）である。
（５）切替時刻ｔ１で第２の２次曲線ｇ（ｔ）が、１つ目の２次曲線ｆ（ｔ）に接する。（６）終了時刻ｔ２で第２の２次曲線ｇ（ｔ２）の値が目標値である。
（７）終了時刻ｔ２で第２の２次曲線ｇ（ｔ）の１階の微分値がゼロである。
　以上の条件から、下記の（３）式～（１３）式が得られる。これらを連立させて、出力の切替時刻ｔ１と終了時刻ｔ２を求め、加速度α１、α２を求める。ここでα１は計算時刻ｔ０から切替時刻ｔ１の加速度であり、α２は切替時刻ｔ１から終了時刻ｔ２の加速度で、α１＝αｐのとき、α２＝αｍ、α１＝αｍのときα２＝αｐである。
Figure JPOXMLDOC01-appb-M000013
　（３）式と（５）式より、
Figure JPOXMLDOC01-appb-M000014
　（３）式と（６）式と（１３）式より、
Figure JPOXMLDOC01-appb-M000015
　（３）式と（７）式より、
Figure JPOXMLDOC01-appb-M000016
　（４）式と（１２）式より、
Figure JPOXMLDOC01-appb-M000017
　（４）式と（１１）式と（１６）式より、
Figure JPOXMLDOC01-appb-M000018
　（４）式と（１０）式と（１６）式と（１７）式より、
Figure JPOXMLDOC01-appb-M000019
　（１８）式を変形して、
Figure JPOXMLDOC01-appb-M000020
　（３）式と（５）式と（９）式と（１５）式と（１６）式と（１７）式と（１９）式より、
Figure JPOXMLDOC01-appb-M000021
　（３）式と（５）式と（８）式と（１５）式と（１６）式と（１７）式と（１９）式と（２０）式より、
Figure JPOXMLDOC01-appb-M000022
　（２１）式を変形して、
Figure JPOXMLDOC01-appb-M000023
　（２２）に２次方程式の解の公式を適用して、
Figure JPOXMLDOC01-appb-M000024
　（２３）式を（２０）式に代入してｔ２を得る。
　ここで、ｆ（ｔ）、ｇ（ｔ）の取り得る軌跡は図４に示すＡ１～Ａ６の６通りとなる。この６通りを場合分けする。
　Ａ１は、Ｘ＞０，Ｖ０＞０のとき、アクチュエータ発生最大減速度αｍで減速しても目標をオーバーする状態である。最大減速度αｍで減速した時に速度がゼロになるまでの時間をｔ３とすると、ｔ３＝Ｖ０／αｍとなり、時刻ｔ３においてオーバーシュートする条件は、Ｖ０×ｔ３／２＝Ｖ０2／２αｍ＞Ｘである。
　Ａ２は、Ｘ＞０，Ｖ０＞０のとき、アクチュエータ発生最大減速度αｍで減速すれば、目標をオーバーしない状態である。Ａ１と同様に計算すると、時刻ｔ３においてオーバーシュートしない条件は、Ｖ０2／２αｍ＜Ｘである。
　Ａ３は、Ｘ＞０，Ｖ０＜０の条件であり、Ａ４は、Ｘ＜０，Ｖ０＞０の条件である。　また、Ａ５は、Ｘ＜０，Ｖ０＜０のとき、アクチュエータ発生最大減速度αｐで減速すれば目標をオーバーしない状態である。Ａ１と同様に計算すると、時刻ｔ３においてオーバーシュートしない条件は、Ｖ０2／２αｐ＞Ｘである。
　Ａ６は、Ｘ＜０，Ｖ０＜０のとき、アクチュエータ発生最大減速度αｐで減速しても目標をオーバーする状態である。Ａ１と同様に計算すると、時刻ｔ３においてオーバーシュートする条件は、Ｖ０2／２αｐ＜Ｘである。
　Ａ１、Ａ２、Ａ３が上凸から下凸への変化であり、α１＝αｍ＜０、α２＝αｐ＞０であり、Ａ４、Ａ５、Ａ６が下凸から上凸への変化であり、α１＝αｐ＞０、α２＝αｍ＜０である。これにより、α１とα２を決定する。
この場合分けを表１に示す。
Figure JPOXMLDOC01-appb-T000025
　この様にして求めたα１、α２に仮想質量ｍを乗じ、アクチュエータ推力を求める。つまり、計算時刻ｔ０～切替時刻ｔ１の間はアクチュエータ推力１（＝α１×仮想質量）で、切替時刻ｔ１～終了時刻ｔ２の間はアクチュエータ推力２（＝α２×仮想質量）となる。
　この本発明の実施の形態のアクチュエータの制御方法で用いる最短時間制御は、図５及び図６に示すような制御フローによって行うことができる。この制御フローがスタートすると、ステップＳ１１で、最大加速度αｐと最大減速度αｍのデータを読み込む。次のステップＳ１２で、目標量（目標値）Ｔｘと制御量（制御値）ｘのデータを読み込む。それとともに、経過時間ｔと再計算用経過時間ｔｃのカウントを始める。ステップＳ１３で，軌道再計算の周期になっているか否か、即ち、軌道計算してからの再計算用経過時間ｔｃが軌道再計算の周期ｔｃｒ以上になったか否かを判定する。ステップＳ１３で，軌道再計算の周期ｔｃｒになった場合には（ＹＥＳ）、ステップＳ１４で、軌道計算し直してから、ステップＳ１５に行き、ステップＳ１３で，起動再計算の周期になっていない場合には（ＮＯ）、ステップＳ１４の軌道計算を迂回してステップＳ１８に行く。なお、軌道再計算の周期ｔｃｒは、制御周期の１／１０程度の周期とするのが好ましいが、軌道再計算の周期ｔｃｒを制御周期と同じにしても問題は生じない。
　ステップＳ１４の軌道計算では、初速Ｖ０を「Ｖ０＝（ｘ－ｘ-1）／ｔｃｒ」（ｘ-1は、計算周期前の制御量）で、運動エネルギーＥを「Ｅ＝Ｖ０2／２」で算出し、表１からα１とα２を決定し、（２３）式と（２０）式より、切替時刻ｔ１と終了時刻ｔ２を計算する。また、再計算用経過時間ｔｃをリセットしてゼロにする（Ｔｃ＝０）。
　次にステップＳ１５では、加速度の制限を行うための加速度制限値（制御出力制限値）を計算する。この加速度制限値は、運動エネルギーＥと残り仕事Ｗの絶対値｜Ｗ｜の和（残留エネルギー）にＥゲインを乗じた値である。つまり、「加速度制限値＝（Ｅ＋｜Ｗ｜）×Ｅゲイン」となる。
　この加速度制限値の最大値は１であるので、ステップＳ１７で、加速度制限値が１より大きいか否かを判定し、１より大きい場合は（ＹＥＳ）、ステップＳ１７で、加速度制限値を１に設定してから、ステップＳ１８に行く。ステップＳ１７で、加速度制限値が１より大きくない場合は（ＮＯ）、そのままとし、ステップＳ１８に行く。
　ステップＳ１８では、図６に示すように、ステップＳ１８に入ると、最初のステップＳ１８ａでは、経過時間ｔが切替時刻ｔ１より小さいか否かを判定し、小さい場合には（ＹＥＳ）、ステップＳ１８ｂに行き、出力加速度を「α１×加速度制限値」にしてから、ステップＳ１９に行く。ステップＳ１８ａで、経過時間ｔが切替時刻ｔ１より小さくない場合には（ＮＯ）、ステップＳ１８ｃに行き、経過時間ｔが終了時刻ｔ２より小さいか否かを判定し、小さい場合には（ＹＥＳ）、ステップＳ１８ｄに行き、出力加速度を「α２×加速度制限値」にしてから、ステップＳ１９に行く。また、ステップＳ１８ｃで、経過時間ｔが終了時刻ｔ２より小さくない場合には（ＮＯ）、ステップＳ１８ｅに行き、出力加速度を「ゼロ」にしてから、ステップＳ１９に行く。
　ステップＳ１９では、出力加速度に相当するアクチュエータ推力を予め設定された時間（各種の判定のインターバルに関係する時間）の間発生し、制御対象を制御する。また、経過時間ｔと再計算用経過時間ｔｃをカウントする。その後、ステップＳ１２に戻り、ステップＳ１２～ステップＳ１９を繰り返す。これにより、再計算用経過時間ｔｃが軌道再計算の周期ｔｃｒ毎にステップＳ１４の軌道計算を計算し直しながら、アクチュエータの推力を制御することができる。
　なお、経過時間ｔが予め設定してある時間を超えたり、制御を終了するスイッチ信号が入力されたりする等の、この図５の制御を終了すべき事情が発生すると、この制御フローのどのステップを実施していても、割り込みが発生し、リターンに行き、上位の制御フローに戻り、上位の制御フローの終了と共に、図５の制御フローも終了する。
　上記の様にして求めた切替時刻ｔ１と終了時刻ｔ２でアクチュエータの動作を切り替えることで、図１９に示すような理想的な制御結果が得られる。しかし、この結果は、摩擦や減衰、誤差などが無い理想的な状況下での結果である。制御対象に機械的な減衰力が存在すると制御結果は図２０に示すように目標と一致しなくなってしまう。
　この問題を解決するために、本発明のアクチュエータの制御方法では、一定周期又は不定期な予め設定した時間毎に目標軌跡を修正する。この再計算による制御軌跡の変化の様子を図３に示す。この図３では、最初の計算の制御軌道（点線）ではＸ１になる予定が実際にはＸ２であった場合に、再計算により、新たな制御軌道（実線）が計算されて、切替時刻ｔ１と終了時刻ｔ２も新たな値となり、この新たな切替時刻ｔ１と終了時刻ｔ２に基づいて制御力が制御されることになる。
　図１７及び図１８に一定周期毎に切替時刻ｔ１と終了時刻ｔ２を計算し直したときの制御結果を示す。信号パルスの立ち上がりのタイミングで再計算を実施している。この最短時間制御方法を、従来技術の最短時間制御方法と区別するために、ここでは、フィードバック（ＦＢ）最短時間制御方法ということにする。図１７に機械系の減衰力が無い場合を、図１８に機械系の減衰力がある場合を示すが、結果として図１８に示すような制御対象に機械的な減衰力がある場合でもほぼ理想的な軌跡で制御対象を目標位置に合致させることができることが分かる。
　しかしながら、制御量が制御目標に合致してもパルス状の制御出力が出たままとなることが分かった。
　本発明では、これを解決して、制御終了後に制御出力を収束させる方法として、制御系が持つ残り仕事と運動エネルギーの和として定義される残留エネルギーの減少と共に、制御出力を小さくしていく方法を採用する。また、この残留エネルギーに対する制御出力の大きさを決定するための係数としてのエネルギーゲイン（Ｅゲイン）を用いる。
　ここでは、目標位置近傍での加速度α２と偏差Ｘ及び質量ｍを掛け合わせた値を残り仕事として定義し、残留エネルギーを残り仕事と運動エネルギーの和として定義する。つまり、「残り仕事＝α２×Ｘ×ｍ」であり、「残留エネルギー＝｜残り仕事｜＋運動エネルギー」となる。ここで、「ｍ」は制御系の質量（仮想的な質量でもよい）であり、「ｍ＝アクチュエータ最大推力／α２」で定義される。
　このように定義した残留エネルギーの時間変化を図７に示す。この図７に示すように、偏差＝０且つ速度＝０のときのみ残留エネルギーがゼロとなるため、制御の評価値として利用することが可能となる。
　そして、制御では、制御の安定性やエネルギー消費量を低減するために、制御終了と共に制御出力がゼロになることが望ましい。そこで、制御終了の評価値として適切な残留エネルギーを制御出力に乗ずることで制御終了と同時に制御出力をゼロにする。具体的には、図８及び図９に示すような「残留エネルギー×Ｅゲインと制御出力制限値」の関係をデータ化したマップデータに従って、制御出力制限値を決定し、この値を最大推力に掛け合わせた値を制御出力とする。この制御出力制限値の最大値は１である。また、一定周期や不定期の予め設定された時間毎における再計算の時に、最大加速度、最大減速度に、この制御出力制限値を掛け合わせた値を最大加速度、最大減速度として再計算する。
　これにより、アクチュエータ推力は、「アクチュエータ推力＝アクチュエータ最大推力×残留エネルギー×Ｅゲイン」となる。図１０及び図１１にこの制御の結果を示す。制御量が制御目標になり、制御対象が目標に合致すると制御出力もゼロになっていることが分かる。
　なお、参考までに、最短時間制御ではベースとするモデルが異なるが、ＰＩＤ制御における残留エネルギーは、品質工学に用いた制御の評価を実施する際に用いる制御の評価指標として使用でき、制御系の持つポテンシャルエネルギーと運動エネルギーの和として定義できる。
　このＰＩＤ制御における制御系のポテンシャルエネルギーは、ばね系の歪エネルギーのことで、目標と制御値の差を歪量Ｘとし、制御系のばね定数ｋ（Ｐ項に相当）とすると「（１／２）・ｋｘ2＝（１／２）・Ｐｘ2」として定義できる。また制御系の持つ運動エネルギーとは制御系の質量をｍ（仮想的な質量でも構わない）と、速度をｄｘ／ｄｔ＝ｖとすると「（１／２）・ｍｖ2」として定義できる。
　ここで制御が終了するとはどういうことか考えてみると、偏差がゼロで且つ速度もゼロの時であることは明らかである。言い換えると、ポテンシャルエネルギー＝０、運動エネルギー＝０の時である。つまり、これらの和である残留エネルギー＝０のときに、制御が終了したと判断できる。残留エネルギーで一般的な制御波形をみると図２２のようになり、オーバーシュートが有ったとしても制御開始から残留エネルギーは単調に減少し、制御終了と同時にゼロとなるので、制御の評価値として適切であることが分かる。
　図１２～図１４に、フィードバック最短時間制御の実施例とＰＩＤ制御の比較例との比較結果を示す。両制御とも、図１２に示す機械系の減衰係数を１とした時を基準に制御ゲインを調整した。図１３と図１４では制御ゲイン固定のまま機械系の減衰係数のみ０．１と２に変更した。まず、図１２をみるとＰＩＤ制御の比較例の方が若干速いものの、ほぼ同等の結果を示している。これに対し減衰係数を小さくした図１３の場合、ＦＢ最短時間制御の実施例では、オーバーシュートの無い制御結果が得られたのに対し、ＰＩＤ制御の比較例では大きなオーバーシュートが発生している。一方、減衰係数を大きくした図１４の場合、フィードバック最短時間制御の実施例は素早く目標に達しているのに対し、ＰＩＤ制御の比較例は目標直前で速度が大きく低下し、目標への応答が遅くなっている。この様にフィードバック最短時間制御の実施例では、制御対象の減衰係数の変化で表わされる抵抗変化に対し、ＰＩＤ制御の比較例よりも安定した制御結果が得られることが分かる。　　　
　更に、図１５にフィードバック最短時間制御の実施例のシミュレーション結果を示し、図１６にＰＩＤ制御の比較例のシミュレーション結果を示す。この結果を比較すると、図１６に示すＰＩＤ制御の比較例では制御系の負荷変動に対し、制御結果が大きく変化してしまうのに対し、図１５に示すフィードバック最短時間制御の実施例では、制御結果に乱れがなく常に安定した制御結果が得られている。つまり、図１６の比較例では、制御系の負荷変動に対し、制御結果が大きく変化してしまうのに対し、図１５の実施例では制御結果に乱れが無く、常に安定した制御結果が得られている。
　従って、上記のアクチュエータの制御方法及びアクチュエータの制御装置によれば、最短時間制御であるために、制御速度が高速となり、最大加速度と最大減速度は計測によって設定でき、これ以外に調整項はないため、制御ゲインを調整する必要が無い。また、ＯＮ／ＯＦＦ制御であるため、中間的な出力を出す必要がなく、コントローラ、ドライバーを簡略化できる。
　また、予め設定した時間毎に、制御の各時刻における目標量と制御量の偏差Ｘを入れて切替時刻ｔ１と終了時刻ｔ２を更新するというフィードバックの要素を取り入れているため、外力が変化しても、また、制御周期を高速にしなくても、常に安定した制御結果を得ることができる。その結果、従来の制御則では相反する大きな課題であった「制御の速度」と「安定性」を両立させることができる。
　その上、制御系が持つ残り仕事と運動エネルギーの和である残留エネルギーの減少と共に、制御出力を小さくすることにより、制御終了後に制御出力を収束させて、ゼロにすることができる。
　本発明のアクチュエータの制御方法及びアクチュエータの制御装置によれば、制御速度が高速となり、しかも、制御ゲインを調整する必要が無く、ＯＮ／ＯＦＦ制御であるため、中間的な出力を出す必要がなく、コントローラ、ドライバーを簡略化でき、また、予め設定した時間毎に、制御の各時刻における目標量と制御量の偏差を入れて切替時刻と終了時刻を更新するというフィードバックの要素を取り入れているため、外力が変化しても、また、制御周期を高速にしなくても、常に安定した制御結果を得ることができ、その上、制御系が持つ残り仕事と運動エネルギーの和である残留エネルギーの減少と共に、制御出力を小さくすることにより、制御終了後に制御出力を収束させて、ゼロにすることができるので、自動車等に搭載した機器など、数多くのアクチュエータの制御方法及びアクチュエータの制御装置として利用できる。
ｔ０　計算時刻
ｔ１　切替時刻
ｔ２　終了時刻
Ｖ０　初速
Ｘ　目標量と制御量の偏差
αｐ　最大加速度
αｍ　最大減速度

Claims (9)

  1. 　最短時間制御を用いると共に、
    　予め、計測されたアクチュエータの制御力の最大出力時の最大加速度（αｐ）と最大減速度（αｍ）を用いて、制御のための計算を行う計算時刻（ｔ０）からの経過時間表示で、加速出力から減速出力へ切り替える切替時刻（ｔ１）と、減速出力の終了時刻（ｔ２）を算出する算出ステップと、
    　前記計算時刻（ｔ０）から前記切替時刻（ｔ１）までは前記アクチュエータの制御力を最大加速出力とし、前記切替時刻（ｔ１）から前記終了時刻（ｔ２）までは前記アクチュエータの制御力を最大減速出力とし、前記終了時刻（ｔ２）で制御力の出力を終了する制御出力ステップを備えると共に、
    　前記算出ステップを予め設定した時間毎に繰り返し、前記切替時刻（ｔ１）と前記終了時刻（ｔ２）を算出して更新する更新ステップと、
    　更に、制御系が持つ残り仕事と運動エネルギーの和である残留エネルギーの減少と共に、制御出力を小さくする制御出力減少ステップを備えたことを特徴とするアクチュエータの制御方法。
  2. 　前記制御出力減少ステップで、前記最大加速出力及び前記最大減速出力のそれぞれに、前記残留エネルギーを乗じ、この乗じた値に、更に、前記残留エネルギーに対して制御出力の大きさを決定するための係数であるエネルギーゲインを乗じた値を制御力の出力とすることを特徴とする請求項１に記載のアクチュエータの制御方法。
  3. 　前記更新ステップで、前記最大加速度、前記最大減速度に、制御出力制限値を掛け合わせた値を新たな最大加速度、最大減速度として再計算することを特徴とする請求項１又は２に記載のアクチュエータの制御方法。
  4. 　前記算出ステップにおいて、前記計算時刻（ｔ０）での位置から前記終了時刻（ｔ２）の目標位置までの軌道を相接する２つの２次曲線の組み合わせで表して、初速をＶ０、目標位置と制御量の偏差をＸとした場合に、前記切替時刻（ｔ１）と、前記終了時刻（ｔ２）を、
    Figure JPOXMLDOC01-appb-M000001
    Figure JPOXMLDOC01-appb-M000002
     　ここで、α１は計算時刻（ｔ０）～切替時刻（ｔ１）の最大加速度（αｐ）又は最大減速度（αｍ）であり、（α２）は切替時刻（ｔ１）～終了時刻（ｔ２）の最大減速度（αｍ）又は最大加速度（αｐ）であり、この（１）式と（２）式で算出することを特徴とする請求項１～３のいずれか１項に記載のアクチュエータの制御方法。
  5. 　最短時間制御を用いると共に、
    　予め、計測されたアクチュエータの制御力の最大出力時の最大加速度（αｐ）と最大減速度（αｍ）を用いて、制御のための計算を行う計算時刻（ｔ０）からの経過時間表示で、加速出力から減速出力へ切り替える切替時刻（ｔ１）と、減速出力の終了時刻（ｔ２）を算出する算出手段と、
    　前記計算時刻（ｔ０）から前記切替時刻（ｔ１）までは前記アクチュエータの制御力を最大加速出力とし、前記切替時刻（ｔ１）から前記終了時刻（ｔ２）までは前記アクチュエータの制御力を最大減速出力とし、前記終了時刻（ｔ２）で制御力の出力を終了する制御出力手段を備えると共に、
    　予め設定した時間毎に前記算出手段により前記切替時刻（ｔ１）と前記終了時刻（ｔ２）を繰り返して算出し、更新する更新手段と、
    　更に、制御系が持つ残り仕事と運動エネルギーの和である残留エネルギーの減少と共に、制御出力を小さくする制御出力減少手段を備えたことを特徴とするアクチュエータの制御装置。
  6. 　前記制御出力減少手段が、前記最大加速出力及び前記最大減速出力のそれぞれに、前記残留エネルギーを乗じ、この乗じた値に、更に、前記残留エネルギーに対して制御出力の大きさを決定するための係数であるエネルギーゲインを乗じた値を制御力の出力とすることを特徴とする請求項５に記載のアクチュエータの制御装置。
  7. 　前記更新手段が、前記最大加速度、前記最大減速度に、制御出力制限値を掛け合わせた値を新たな最大加速度、最大減速度として再計算することを特徴とする請求項５又は６に記載のアクチュエータの制御装置。
  8. 　前記算出手段が、前記計算時刻（ｔ０）での位置から前記終了時刻（ｔ２）の目標位置までの軌道を相接する２つの２次曲線の組み合わせで表して、初速をＶ０、目標位置と制御量の偏差をＸとした場合に、前記切替時刻（ｔ１）と、前記終了時刻（ｔ２）を、
    Figure JPOXMLDOC01-appb-M000003
    Figure JPOXMLDOC01-appb-M000004
     　ここで、α１は計算時刻（ｔ０）～切替時刻（ｔ１）の最大加速度（αｐ）又は最大減速度（αｍ）であり、α２は切替時刻（ｔ１）～終了時刻（ｔ２）の最大減速度（αｍ）又は最大加速度（αｐ）であり、この（１）式と（２）式で算出することを特徴とする請求項５～７のいずれか１項に記載のアクチュエータの制御装置。
  9. 　計算時刻（ｔ０）から切替時刻（ｔ１）の間の加速度（αｐ）と切替時刻（ｔ１）から終了時刻（ｔ２）の間の加速度（αｍ）を表１によって決定する
    Figure JPOXMLDOC01-appb-T000005
     　ことを特徴とする請求項１～８のいずれか１項に記載のアクチュエータの制御装置。
PCT/JP2012/070480 2011-09-16 2012-08-10 アクチュエータの制御方法及びアクチュエータの制御装置 WO2013038856A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
EP12831914.2A EP2757427B1 (en) 2011-09-16 2012-08-10 Actuator control method and actuator control device
US14/345,167 US9519275B2 (en) 2011-09-16 2012-08-10 Actuator control method and actuator control device
CN201280044822.9A CN103797429B (zh) 2011-09-16 2012-08-10 致动器的控制方法以及致动器的控制装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2011-202769 2011-09-16
JP2011202769A JP5838681B2 (ja) 2011-09-16 2011-09-16 アクチュエータの制御方法及びアクチュエータの制御装置

Publications (1)

Publication Number Publication Date
WO2013038856A1 true WO2013038856A1 (ja) 2013-03-21

Family

ID=47883092

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2012/070480 WO2013038856A1 (ja) 2011-09-16 2012-08-10 アクチュエータの制御方法及びアクチュエータの制御装置

Country Status (5)

Country Link
US (1) US9519275B2 (ja)
EP (1) EP2757427B1 (ja)
JP (1) JP5838681B2 (ja)
CN (1) CN103797429B (ja)
WO (1) WO2013038856A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180304048A1 (en) * 2015-04-09 2018-10-25 Vigmed Ab Needle tip shielding device and catheter hub therefore

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI578123B (zh) * 2016-08-17 2017-04-11 東台精機股份有限公司 粉末積層製造之檢測修補裝置及其方法
US10792773B2 (en) * 2016-09-26 2020-10-06 Esko-Graphics Kongsberg As Actuators and methods for controlling tools
CN108060417A (zh) 2016-11-07 2018-05-22 东台精机股份有限公司 粉末积层制造的检测修补装置及其方法
JP7119827B2 (ja) * 2018-09-21 2022-08-17 トヨタ自動車株式会社 異常診断方法、異常診断装置及びプログラム
JP7271834B2 (ja) * 2019-03-29 2023-05-12 ニデックパワートレインシステムズ株式会社 電動アクチュエータ、およびアクチュエータ装置
US11346377B2 (en) 2020-08-24 2022-05-31 Epiroc Drilling Solutions, Llc System and method for automatic calibration of actuators

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS58119004A (ja) * 1982-01-11 1983-07-15 Hitachi Ltd バングバング制御方法
JPS6177905A (ja) * 1984-09-26 1986-04-21 Mitsubishi Electric Corp 位置決め制御回路
JPS62120515A (ja) * 1985-11-20 1987-06-01 Fuji Electric Co Ltd 位置決め制御方式
JPH05224702A (ja) * 1992-02-14 1993-09-03 Omron Corp バンバン制御装置
JPH11296234A (ja) * 1998-04-10 1999-10-29 Matsushita Electric Ind Co Ltd モータ類の位置制御時の加重制御方法およびその装置
JP2000094371A (ja) 1998-09-18 2000-04-04 Yaskawa Electric Corp ロボットの最短時間制御装置

Family Cites Families (39)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3835360A (en) * 1970-08-26 1974-09-10 Tektronix Inc Numerical control system
US4200827A (en) * 1977-06-29 1980-04-29 International Business Machines Corporation Positioning system employing feedforward and feedback control
JPS57120112A (en) * 1980-12-22 1982-07-27 Fujitsu Ltd Locus control method of arm
EP0103493B1 (en) * 1982-09-15 1987-12-16 Ampex Corporation Improved bang-bang controlled velocity command generator
JPH01164280A (ja) * 1987-12-21 1989-06-28 Fanuc Ltd 加減速制御方式
CN1051300A (zh) * 1989-10-28 1991-05-15 李庆东 一种将鬼针草制造成口服液的方法
US5119250A (en) 1990-06-15 1992-06-02 International Business Machines Corporation Method and apparatus for performing a seek
EP0507907B1 (en) 1990-09-18 1999-01-13 Rodime PLC Digital servo control system for use in disk drives
US5331264A (en) * 1993-04-15 1994-07-19 Fanuc Robotics North America, Inc. Method and device for generating an input command for a motion control system
US5508596A (en) * 1993-10-07 1996-04-16 Omax Corporation Motion control with precomputation
JP3159287B2 (ja) 1994-05-02 2001-04-23 キヤノン株式会社 位置決め制御装置
US5677609A (en) 1994-07-28 1997-10-14 National Semiconductor Corporation Intelligent servomechanism controller
JPH0916241A (ja) * 1995-06-29 1997-01-17 Fanuc Ltd ロボットの加減速動作の設定方法
JP2812920B2 (ja) * 1996-05-30 1998-10-22 川崎重工業株式会社 ロボットの動作指令作成方法
US6002184A (en) 1997-09-17 1999-12-14 Coactive Drive Corporation Actuator with opposing repulsive magnetic forces
JP2000190262A (ja) * 1998-12-22 2000-07-11 Denso Corp ロボットの制御装置
JP2001100802A (ja) * 1999-09-30 2001-04-13 Matsushita Electric Works Ltd 加減速制御装置
US7135003B2 (en) * 2001-06-29 2006-11-14 Honda Giken Kogyo Kabushiki Kaisha Feedback estimation of joint forces and joint moments
US7623944B2 (en) * 2001-06-29 2009-11-24 Honda Motor Co., Ltd. System and method of estimating joint loads in a three-dimensional system
US6686716B1 (en) 2001-07-18 2004-02-03 Itt Manufacturing Enterprises, Inc. Tuned open-loop switched to closed-loop method for rapid point-to-point movement of a periodic motion control system
JP3881593B2 (ja) * 2002-05-17 2007-02-14 株式会社山武 制御装置
JP4112900B2 (ja) * 2002-05-21 2008-07-02 株式会社山武 制御方法及び制御装置
JP4059213B2 (ja) * 2004-02-27 2008-03-12 ブラザー工業株式会社 モータ制御装置およびモータ制御方法
DE102004054776B3 (de) * 2004-11-12 2006-03-16 Bayerische Motoren Werke Ag Verfahren zur Kalibrierung eines Wegsensors einer Drehaktuatorvorrichtung zur Ansteuerung eines Gaswechselventils einer Brennkraftmaschine
US20060190136A1 (en) * 2005-02-24 2006-08-24 Paragon Technical Services, Inc. Method and system for controlling a manipulator
SE528539C2 (sv) * 2005-05-20 2006-12-12 Aros Electronics Ab Metod och anordning för återvinningsbart upptagande av elektrisk retardationsenergi från ett industrirobotsystem
US8082062B2 (en) * 2005-06-10 2011-12-20 Honda Motor Co., Ltd. Regenerative actuation in motion control
EP1815972B1 (en) * 2006-02-06 2013-12-18 ABB Research Ltd. Press line system and method
JP2008090513A (ja) * 2006-09-29 2008-04-17 Fanuc Ltd 数値制御装置
US8024052B1 (en) 2007-03-30 2011-09-20 Tim Hakala Adaptive mapping of device output change to amounts of control effort
US8335600B2 (en) * 2007-12-14 2012-12-18 The Boeing Company Regenerative integrated actuation system and associated method
JP4920054B2 (ja) * 2009-03-30 2012-04-18 株式会社日立製作所 車両運動制御装置
US8209100B2 (en) * 2009-05-06 2012-06-26 GM Global Technology Operations LLC Methods and systems for controlling braking in vehicles
JP4768086B2 (ja) * 2009-07-02 2011-09-07 パナソニック株式会社 ロボット、ロボットアームの制御装置、及びロボットアームの制御プログラム
JP5445932B2 (ja) 2009-10-13 2014-03-19 横河電機株式会社 バルブ制御装置およびバルブ制御方法
US20110276216A1 (en) * 2010-05-07 2011-11-10 Texas Instruments Incorporated Automotive cruise controls, circuits, systems and processes
EP2422935B1 (en) * 2010-08-31 2015-02-25 Kabushiki Kaisha Yaskawa Denki Robot, robot system, robot control device, and state determining method
WO2013086002A1 (en) 2011-12-05 2013-06-13 Cellworks Research India Private Limited Compositions, process of preparation of said compositions and method of treating cancer
US9228645B2 (en) * 2013-06-11 2016-01-05 Raytheon Company Vacuum stable mechanism drive arm

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS58119004A (ja) * 1982-01-11 1983-07-15 Hitachi Ltd バングバング制御方法
JPS6177905A (ja) * 1984-09-26 1986-04-21 Mitsubishi Electric Corp 位置決め制御回路
JPS62120515A (ja) * 1985-11-20 1987-06-01 Fuji Electric Co Ltd 位置決め制御方式
JPH05224702A (ja) * 1992-02-14 1993-09-03 Omron Corp バンバン制御装置
JPH11296234A (ja) * 1998-04-10 1999-10-29 Matsushita Electric Ind Co Ltd モータ類の位置制御時の加重制御方法およびその装置
JP2000094371A (ja) 1998-09-18 2000-04-04 Yaskawa Electric Corp ロボットの最短時間制御装置

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180304048A1 (en) * 2015-04-09 2018-10-25 Vigmed Ab Needle tip shielding device and catheter hub therefore
US11202887B2 (en) * 2015-04-09 2021-12-21 Vigmed Ab Needle tip shielding device and catheter hub therefore

Also Published As

Publication number Publication date
US20140350730A1 (en) 2014-11-27
CN103797429B (zh) 2016-03-16
US9519275B2 (en) 2016-12-13
EP2757427B1 (en) 2019-11-20
EP2757427A4 (en) 2016-05-18
JP5838681B2 (ja) 2016-01-06
EP2757427A1 (en) 2014-07-23
JP2013065152A (ja) 2013-04-11
CN103797429A (zh) 2014-05-14

Similar Documents

Publication Publication Date Title
JP5838681B2 (ja) アクチュエータの制御方法及びアクチュエータの制御装置
CN108319144B (zh) 一种机器人轨迹跟踪控制方法及***
US7902785B2 (en) Method and device for guiding the movement of a moveable machine element of a machine
JP6209176B2 (ja) 位置制御装置
JP2004213472A (ja) 制御装置
JP6020537B2 (ja) モータ制御装置及びモータ制御方法
JP2737669B2 (ja) 自動演奏ピアノの鍵盤駆動装置
CN108983703A (zh) 超精密运动***前馈控制器参数整定方法
JP2016134975A (ja) 振動を抑制するモータ制御装置
JP5803361B2 (ja) アクチュエータの制御方法及びアクチュエータの制御装置
JPWO2007105527A1 (ja) 位置決め機構の制御方法および制御装置
JP5834628B2 (ja) アクチュエータの制御方法及びアクチュエータの制御装置
Bucci et al. Nonlinear control algorithm for improving settling time in systems with friction
JP2002186269A (ja) 電動機制御方法
JP4183057B2 (ja) 数値制御システム
WO2002025390A1 (fr) Procede de servocommande
JP6340160B2 (ja) 送り駆動系の設計方法
JP6886148B2 (ja) アクティブマスダンパおよびアクティブマスダンパの制御装置
JP5870668B2 (ja) エネルギ評価制御方法及びエネルギ評価制御装置
JP6631745B1 (ja) リニアモータステージの制御方法及び制御装置
Sievert et al. Numerical and Experimental Investigation of Time-Varying Vibration Control for Beam Subjected to Moving Masses
JP5862265B2 (ja) エネルギ評価制御方法及びエネルギ評価制御装置
JP5862264B2 (ja) エネルギ評価制御方法及びエネルギ評価制御装置

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 201280044822.9

Country of ref document: CN

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

Ref document number: 12831914

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 14345167

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2012831914

Country of ref document: EP