WO2006051651A1 - モーション制御装置とシステム同定方法 - Google Patents

モーション制御装置とシステム同定方法 Download PDF

Info

Publication number
WO2006051651A1
WO2006051651A1 PCT/JP2005/017592 JP2005017592W WO2006051651A1 WO 2006051651 A1 WO2006051651 A1 WO 2006051651A1 JP 2005017592 W JP2005017592 W JP 2005017592W WO 2006051651 A1 WO2006051651 A1 WO 2006051651A1
Authority
WO
WIPO (PCT)
Prior art keywords
hat
value
identification
inertia
viscous friction
Prior art date
Application number
PCT/JP2005/017592
Other languages
English (en)
French (fr)
Inventor
Jun Hagihara
Original Assignee
Kabushiki Kaisha Yaskawa Denki
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 Kabushiki Kaisha Yaskawa Denki filed Critical Kabushiki Kaisha Yaskawa Denki
Priority to JP2006544794A priority Critical patent/JP4919157B2/ja
Publication of WO2006051651A1 publication Critical patent/WO2006051651A1/ja
Priority to US11/797,945 priority patent/US7626351B2/en

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P23/00Arrangements or methods for the control of AC motors characterised by a control method other than vector control
    • H02P23/14Estimation or adaptation of motor parameters, e.g. rotor time constant, flux, speed, current or voltage
    • 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/41Servomotor, servo controller till figures
    • G05B2219/41381Torque disturbance observer to estimate inertia

Definitions

  • the present invention relates to a motion control device and a system identification method for accurately estimating inertia, viscous friction coefficient, and constant disturbance to be controlled.
  • a conventional device that estimates the inertia of a motor integrates the torque command value and the torque command value of the model within a certain period of time, and multiplies the ratio by the inertia nominal value to perform the estimation.
  • a method is used to remove certain disturbances such as viscous friction, Coulomb friction, and gravity to increase accuracy by using a specific motion. (For example, see Patent Document 1).
  • 3 is an electric motor, and 4 is a machine coupled to the electric motor.
  • a detector 5 is added to the motor 3.
  • Reference numeral 71 denotes a command generation unit that outputs a motor operation speed command vref.
  • Reference numeral 72 denotes a speed control unit that performs proportional-integral control so that the command and motor speed match, and outputs a torque command value Tref.
  • 75 is a current control unit that outputs a current value so that the motor operates according to the torque command Tref.
  • Reference numeral 73 denotes an estimation unit having a motor model, which performs proportional-integral control so that the command and the model speed match, and outputs a model torque command Tref '.
  • the estimation unit 74 is an identification unit, which integrates the actual torque command Tref and the torque command Tref 'input to the model by the estimation unit 73 in the interval of time [a, b], and adds the ratio to the nominal street of the inertia, Multiply to find the estimated inertia.
  • the estimated inertia value theoretically matches the actual inertia J.
  • the integrated value of the torque command value Tref in the time [a, b] section does not include the effect of viscous friction or Coulomb friction. Such constraints are necessary. In addition, special measures are required to remove constant disturbances such as gravity.
  • Patent Document 1 as a method for realizing these, the operation shown in FIG. 10 is used as a constraint condition. From (A) to (F) in Fig. 10, the force is a graph with time on the horizontal axis and speed on the vertical axis.
  • Patent Literature l WO96Z37039 (Page 5-8, Fig. 4 and Fig. 5)
  • the present invention has been made in view of such problems. Even when there is an influence of a constant disturbance such as viscous friction, Coulomb friction, or gravity, regardless of the operation, it is possible to control with high accuracy by simple calculation.
  • Target inertia identification town—hat and viscous friction coefficient identification value D—hat and constant disturbance identification value C—hat can be calculated, and the effect of noise by time differentiation used in subsequent calculations by filtering
  • Inertia can be identified with a simple calculation with only 3 divisions and 5 time differentiation operations, and the friction constant can be identified with a simple calculation with only 1 subtraction and 1 multiplication.
  • the present invention there is a case where there is an influence of a constant disturbance such as viscous friction, Coulomb friction, or gravity.
  • a constant disturbance such as viscous friction, Coulomb friction, or gravity.
  • the noise of the time differential operation used in the subsequent calculations can be calculated.
  • Provides a motor controller and its system identification method that can identify both inertia and friction constant with a simple calculation of 6 times, multiplication 6 times, subtraction 3 times, and division 2 times. For the purpose.
  • the present invention controls the motor based on the command generator, the command output from the command generator, and the detection value obtained from the detector attached to the motor.
  • System with a controller that outputs the current value for detection, and an identifier that identifies the inertia identification city-hat and viscous friction coefficient identification value D-hat and constant disturbance identification value C-hat to be controlled
  • the acceleration detection value afb is calculated by differentiating the speed detection value Vfb with time
  • the signal Fafb obtained by filtering the acceleration detection value afb is calculated
  • the signal Fvfb obtained by filtering the speed detection value Vfb is calculated.
  • Ftref is obtained by filtering the torque command value Tref, and the time derivative and the four arithmetic operations are performed based on Ftref, Fvfb, and Fafb, and the inertia identification ⁇ —hat and viscous friction coefficient identification value D— hat and set Is the was convex steps of calculating a turbulent identified value C-hat.
  • the third aspect of the present invention in the system identification method of the motion control apparatus according to the second aspect, in the process of calculating the viscous friction coefficient identification value D-hat, after calculating the inertia identification value J-hat, W2 The inertia identification ffi-hat is multiplied to calculate W3, and the result of subtracting W3 from W1 is used as the identification value D-hat for the viscous friction coefficient.
  • the present invention described in claim 4 provides a system identification method for the motion control apparatus described in claim 3.
  • the inertia identification ffi-hat and the viscous friction coefficient identification value D-hat are calculated, and then Fvfb is multiplied by the viscous friction identification value D-hat.
  • TF is calculated
  • TA is calculated by multiplying Fafb and inertia identification ffi-hat
  • subtracting TF and TA from Ftref is the constant disturbance identification value C-hat.
  • the filter used in the filtering process is a filter having a denominator order that is at least three orders larger than the numerator order. is there.
  • the present invention according to claim 6 is the system identification method of the motion control apparatus according to claims 1 to 5, wherein the internal state quantity is monitored, and only when the state quantity exceeds a set value.
  • Inertia identification city to be controlled hatch and viscous friction coefficient identification value D—hat and constant disturbance identification value C—hat are calculated. If the state quantity is not equal to or greater than the set value, J-hat, D-hat , C-hat identification values are retained as previous values.
  • the state quantity to be monitored is the speed Vfb, acceleration afb, torque command Tref, or Fvfb obtained by filtering these signals. , Fafb, Ftref, Dvfb, Dw2, or a combination of them.
  • the present invention described in claim 8 is a command generator, a position controller that performs position control based on a command output from the command generator and a detection value obtained from a detector cover added to the motor.
  • the equalizer uses a time differentiator that calculates the acceleration detection value afb by differentiating the speed detection value Vfb with respect to time, and filters the acceleration detection value afb with Fafb
  • a first filter that calculates Ftf by filtering the
  • the present invention according to claim 9 is the motion control apparatus according to claim 8, wherein the JDC estimator includes an inertia estimator that estimates inertia, a viscous friction estimator that estimates viscous friction,
  • a constant disturbance estimator for estimating the constant disturbance is provided.
  • the present invention outputs a current value for driving the motor based on a command generator, a command output from the command generator, and a detection value obtained from a detector force added to the motor.
  • the speed detection value Vfb is time-differentiated to calculate the acceleration detection value afb, the signal Fafb obtained by filtering the acceleration detection value afb is calculated, the signal Fvfb obtained by filtering the speed detection value Vfb is calculated, Calculates the signal Ftref obtained by filtering the torque command value Tre f, calculates Dtref, Dvfb, Dafb, and D2tref, D2vfb, D2afb, which are time-differentiated Ftref, Fvfb, Fafb.
  • Ri the speed detection value
  • J_hat (D2trefDvfb DtrefD2vfb) /
  • the present invention according to claim 11 is the system identification method of the motion control apparatus according to claim 10, wherein the filter used in the filter processing is a filter having a denominator order of at least a third order larger than a numerator order. It is characterized by.
  • the present invention according to claim 12 is the system identification method of the motion control apparatus according to claim 10, wherein the filter used in the filter processing is a series of a plurality of moving average filters connected in series. Is.
  • the present invention according to claim 13 monitors an internal signal of the motion control device, and only when the signal becomes a set value or more, the inertia identification town hat and the viscous friction coefficient identification value of the control object. D-hat and constant disturbance identification value C-hat are calculated. If the state quantity is not equal to or greater than the set value, the previous values are set as the same values for J-hat, D-hat, and C-hat. It is characterized by holding.
  • the signal to be monitored may be velocity Vfb, acceleration afb, torque command Tref, or Fvfb, Fafb, Ftref obtained by filtering these signals, or time differentiation of those signals. It is one of Dtref, Dvfb, Dafb, D2tref, D2vfb, D2afb, or a combination thereof.
  • the present invention includes a command generator for generating a command, and a current controller for generating a current for driving the motor based on the command and a detection value obtained from a detector added to the motor. And a motion control device equipped with an identifier for identifying the inertia identification city—hat and viscous friction coefficient identification value D—hat and constant disturbance identification value C—hat to be controlled.
  • An acceleration calculator for calculating the acceleration detection value afb, a first filter for calculating the signal Fafb obtained by filtering the acceleration detection value afb, and a second filter for calculating the signal Fvfb obtained by filtering the velocity detection value Vfb.
  • a third filter for calculating a signal Ftref obtained by filtering the torque command value Tref, and Dtref, Dvfb, Dafb, and another time derivative obtained by differentiating the Ftref, Fvfb, and Fafb with respect to time.
  • J_hat (D2trefDvfb DtrefD2vfb) /
  • the speed, acceleration, and torque command values can be calculated using an expression that is deformed (using differentiation or division) without approximating the relationship between the torque command, viscous friction, and constant disturbance.
  • the identification calculation is performed simply by performing time differentiation and four arithmetic operations, so even if there is an influence of constant disturbances such as viscous friction, Coulomb friction, and gravity, It is possible to accurately identify the inertia identification city-hat and viscous friction coefficient identification value D-hat and constant disturbance identification value C-hat to be controlled.
  • the effect of noise due to the time differentiation used in subsequent calculations can be reduced by performing filtering.
  • the inertia can be identified by only three divisions and five time differential operations, it is possible to identify with a simple calculation and a small amount of calculation.
  • the identification calculation is not performed when the state quantity being monitored is equal to or smaller than the set value, the accuracy is deteriorated by dividing by a minute value. In any case, the identification accuracy can be maintained.
  • an equation obtained by differentiating both sides of the equation with respect to time there are two formulas: an equation obtained by differentiating both sides of the equation with respect to time, and an equation obtained by differentiating both sides of the equation with respect to time.
  • the identification calculation is performed simply by performing time differentiation and four arithmetic operations using the signals after filtering the speed, acceleration, and torque command values. Even when there is an influence of constant disturbances such as viscous friction, Coulomb friction, and gravity, the inertia identification city-hat and viscous friction coefficient identification value to be controlled accurately.
  • D hatch and constant disturbance identification value C—hat can be accurately identified.
  • the filtering process reduces the effects of noise due to the time differentiation used in subsequent calculations.
  • the denominator order is 3 or more than the numerator order of the filter, it is possible to perform three-time time differentiation performed by the inertia identification calculation without being affected by noise. Can be.
  • the difference in the phase delay of each variable can be eliminated, so that the identification can be performed with higher accuracy.
  • the identification calculation is not performed when the monitored signal is equal to or lower than the set value, the accuracy is prevented from being deteriorated by dividing by a minute value. In any case, the identification accuracy can be maintained.
  • FIG. 1 is a block diagram showing the configuration of a motion control apparatus to which the method of the present invention is applied.
  • FIG. 2 is a block diagram showing a configuration of a second motion control apparatus to which the method of the present invention is applied.
  • FIG. 3 is a block diagram showing a configuration of an identifier of the motion control apparatus to which the first method of the present invention is applied.
  • FIG. 7 is a block diagram showing the configuration of an identifier of a motion control apparatus to which the second method of the present invention is applied.
  • FIG. 9 is a block diagram showing a configuration of a motion control apparatus to which a conventional method is applied.
  • FIG. 1 is a block diagram showing the configuration of a motion control apparatus that implements the method of the present invention.
  • 1 is a command generator that generates a position command Xref.
  • 2 represents the controller, which performs control calculation based on the position command, position detection value Xfb, and speed detection value Vfb, and outputs current I. Any calculation in controller 2 can be performed, but in this embodiment, position controller 7 that outputs speed command Vref based on position command Xref and position detection value Xfb, speed command Vref and speed detection value It consists of a speed controller 8 that outputs a torque command value Tref based on Vfb and a current controller 9 that controls the current to flow according to the torque command value.
  • the speed detection value may be a time derivative of the position detection value.
  • an approximate derivative divided by the post-difference control cycle may be used as the time differentiation calculation.
  • 3 is an electric motor, and 4 machines are connected to the electric motor 3.
  • 5 represents a detector that detects the position and speed of the motor 3.
  • 6 is an identifier that calculates the inertia identification city-hat, viscous friction coefficient identification value D-hat, and constant disturbance identification value C-hat to be controlled based on the torque command Tref and the velocity detection value Vfb.
  • FIG. 3 is a block diagram showing the processing inside the six identifiers.
  • 10 represents a time differentiator
  • the speed detection value Vfb is time-differentiated to calculate the acceleration detection value afb.
  • approximate differentiation that divides by the control period after the difference between the current and previous signals may be used as the time differentiation operation.
  • the filter used here is preferably a filter whose denominator is 3rd order or higher than the numerator order. If you use something like the function Gfil,
  • s represents the Laplace operator.
  • any filter can be used. It is perfectly possible to use three or more primary low-pass filters in series.
  • 14 ⁇ represents a DC estimator.
  • Input Ftref, Fvfb, and Fafb which are obtained by filtering torque command Tref, speed detection value Vfb, and acceleration detection value afb, and identify the inertia to be controlled ffi-hat and viscous friction Coefficient identification value D-hat and constant disturbance identification value C-hat are calculated.
  • the interior of 14 consists of three forces: 15 inertia identifiers, 16 viscous friction identifiers, and 17 constant disturbance identifiers.
  • FIG. 4 is a block diagram showing processing of the inertia identifier 15.
  • 41, 42, 43, 46 and 47 represent time differentiators
  • 44, 45 and 48 represent dividers.
  • Inertia identification town —hat is calculated as follows.
  • Dtref, Dvfb, and Dafb are calculated by differentiating Ftref, Fvfb, and Fafb, respectively.
  • Dtref is divided by the Dvfb to calculate W1
  • Dafb is divided by the Dvfb.
  • W1 is time-differentiated to calculate DW1
  • W2 is time-differentiated to calculate DW2
  • the value obtained by dividing DW1 by DW2 is the inertia identification town-hat and To do.
  • FIG. 5 is a block diagram showing processing of the viscous friction identifier 16.
  • 51 represents a multiplier and 52 represents a subtractor.
  • the viscous friction identification value D-hat is calculated according to the following procedure. First, W2 is multiplied by the inertia identification city hat to calculate W3, and then W1 is subtracted from W3 to obtain the viscous friction coefficient identification value D-hat. This process is expressed by the following equation (10) and force equation (11).
  • FIG. 6 is a block diagram showing processing of the constant disturbance identifier 17.
  • 61 and 62 indicate multipliers, and 63 indicates a subtractor that subtracts two variables from one variable.
  • the constant disturbance identification value D —hat is calculated as follows.
  • TF is calculated by multiplying Fvfb and the viscous friction identification value D-hat, then TA is calculated by multiplying Fafb and inertia identification ⁇ -hat, and then TF and previous TF from Ftref.
  • the value obtained by subtracting TA is the constant disturbance identification value C_hat.
  • the above processing is not always calculated, and the speed Vfb or acceleration afb or one of the torque command values Tref is monitored. Only when the variable is equal to or greater than a preset threshold value, the identification accuracy can be improved by performing the above processing.
  • variable to be monitored is one of Fvfb, Fafb, Ftref, which is a filter of Vfb, afb, Tref, or a combination of them.
  • the last value that was used for the calculation is stored only as the value of the identification value.
  • Each identification value can be updated.
  • the monitored variable is smaller than the threshold value, and sometimes the controlled inertia identification town—hat, viscous friction coefficient identification value D—hat, and constant disturbance identification value C—hat are not calculated. However, it is possible not to do so, but only to perform some operations.
  • a current command may be used instead of a torque command.
  • a torque constant (thrust constant) is used to convert the unit of torque and current.
  • the configuration for performing the position control has been described.
  • the inertia identification town-hat and the control target are identified in exactly the same manner.
  • the viscous friction coefficient identification value D-hat and the constant disturbance identification value C-hat can be identified.
  • the command generator 1 generates the speed command Vref, and only the calculations of the speed controller and current controller are performed inside the controller 2.
  • Expression (16) is obtained by performing both-side filter processing.
  • Equation (22) D can be calculated as in Equation (22).
  • FIG. 1 is a block diagram showing a configuration of a motion control apparatus for carrying out the method of the present invention. Since this diagram has been described in the first embodiment, the description thereof will be omitted.
  • FIG. 7 is a block diagram showing processing inside the identifier 6, which is different from that in FIG. 3 and will be described below.
  • 10 represents a time differentiator
  • the speed detection value Vfb is time differentiated to calculate the acceleration detection value afb.
  • approximate differentiation by dividing by the control period after the difference between the current and previous signals may be used as the time differentiation operation.
  • Reference numerals 11, 12, and 13 denote a first filter, a second filter, and a third filter, respectively.
  • the filter used here is better if the denominator order is 3 or higher than the numerator order. For example, use a filter whose transfer characteristic is expressed by the transfer function Gfil in Eq. (24).
  • s represents the Laplace operator.
  • more than 3 primary low-pass filters can be used. There is no problem even if they are used in series. There is no problem if three or more moving average filters are used in series. In this case, since there is no phase delay due to filtering, the estimation accuracy can be further increased.
  • 14 ⁇ DC estimator, Ftref, Fvfb, Fafb filtered by torque command Tref, speed detection value Vfb, acceleration detection value afb, and input of inertia identification ffi-hat and viscosity Calculate friction coefficient identification value D-hat and constant disturbance identification value C-hat.
  • the 14 JDC estimators consist of two forces: 18 inertias, viscous friction identifiers, and 16 constant disturbance identifiers.
  • FIG. 8 is a block diagram showing the processing of the inertia viscous friction identifier 15.
  • 41, 42, 43, 49, 4A and 4B represent time differentiators
  • 4C, 4D, 4E, 4F, 4G and 4H represent multipliers
  • 41, J and 4K represent subtractors
  • 4L and 4M represent dividers.
  • Inertia identification constant-hat and viscous friction coefficient identification value D-hat are calculated according to the following procedure.
  • Ftref, Fvfb, and Fafb are time-differentiated to calculate Dtref, Dvfb, and Dafb, respectively, and then time-differentiated again to calculate D2tref, D2vfb, and D2afb.
  • approximate differentiation that divides by the control period after the difference between the current and previous signals may be used as the time differentiation operation.
  • equations (25) and (26) are calculated to estimate the inertia identification town hat and the viscous friction coefficient identification value D_hat.
  • J_hat (D2trefDvfb DtrefD2vfb) /
  • FIG. 6 is a block diagram showing processing of the constant disturbance identifier 17. Since the processing inside the constant disturbance identifier has been described in the first embodiment, a description thereof will be omitted.
  • the speed Vfb or acceleration afb is not sufficient for the above processing to always be calculated. Is one of the torque command values Tref, and the combination thereof is monitored, and the identification accuracy can be improved by performing the above processing only when those variables are equal to or greater than a preset threshold value.
  • Variables to be monitored are velocity Vfb, acceleration afb, torque command Tref or Fvfb, Fafb, Ftref obtained by filtering these signals, or D tref, Dvfb, Dafb, D2tref, D2vfb obtained by time-differentiating those signals. Either one of D2afb or its combination.
  • the last value that was used for the calculation is stored only as the value of the identification value. Even if each identification value is updated.
  • the monitored variable is smaller than the threshold value, and sometimes the controlled inertia identification town—hat, viscous friction coefficient identification value D—hat, and constant disturbance identification value C—hat are not calculated. However, it is possible not to do so, but only to perform some operations.
  • a current command may be used instead of a torque command.
  • torque and current units are converted using a torque constant (thrust constant).
  • Equation (30) is obtained.
  • D2tref J-D2af b + D-D2vf b (30)

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Electric Motors In General (AREA)

Abstract

 粘性摩擦や一定外乱が作用する場合も、動作によらず、精度良くシステム定数(イナーシャ、粘性摩擦係数、一定外乱)を同定する。    同定器(6)は速度検出値Vfbを時間微分して加速度検出値afbを算出する時間微分器(10)と、加速度検出値afbをフィルタ処理してFafbを算出する第1のフィルタ(11)と、前記速度検出値Vfbをフィルタ処理してFvfbを算出する第2のフィルタ(12)と、前記トルク指令値Trefをフィルタ処理してFtrefを算出する第3のフィルタ(13)と、前記FtrefとFvfbとFafbをもとに時間微分および四則演算し、前記制御対象のイナーシャ同定値J_hatおよび粘性摩擦係数同定値D_hatおよび一定外乱同定値C_hatを算出するJDC推定器(14)とで構成される。

Description

明 細 書
モーション制御装置とシステム同定方法
技術分野
[0001] 本発明は、制御対象のイナーシャおよび粘性摩擦係数および一定外乱を精度良く 推定するモーション制御装置とシステム同定方法に関する。
背景技術
[0002] 従来のモータのイナ一シャを推定する装置は、トルク指令値とモデルのトルク指令 値をある時間内積分し、それらの比にイナーシャノミナル値を乗じて推定を行ってお り、その際、特定の動作を用いて、粘性摩擦や、クーロン摩擦、重力などの一定外乱 を除去し精度を上げるという方法を用いている。(例えば、特許文献 1参照)。
図 9において、 3は電動機であり、 4は電動機に結合された機械である。電動機 3に は検出器 5が付加されている。 71は指令発生部であり、モータの動作速度指令 vref を出力する。 72は速度制御部であり、指令とモータ速度が一致するように比例積分 制御を行い、トルク指令値 Trefを出力する。 75はトルク指令 Trefどおりにモータが動 作するように電流値を出力する電流制御部である。 73は推定部であり、モータのモ デルを有しており、指令とモデルの速度が一致するように比例積分制御を行い、モデ ルのトルク指令 Tref'を出力する。 74は同定部であり、実際のトルク指令 Trefおよび 、推定部 73でモデルに入力されるトルク指令 Tref'を時間 [a, b]の区間で積分し、 その比にイナ一シャのノミナル街,を乗じてイナーシャ推定値を求める。この方法で は外乱が全く存在しない時には理論的にイナーシャ推定値は実際のイナーシャ Jと一 致する。
しカゝしながら、粘性摩擦やクーロン摩擦、重力などの一定外乱が存在する場合は、 時間 [a, b]区間でのトルク指令値 Trefの積分値が粘性摩擦やクーロン摩擦の影響 を含まないような制約条件が必要になる。また、重力などの一定外乱を除去するため にも特別な工夫が必要になる。
特許文献 1では、これらを実現する方法として、図 10に示すような動作を制約条件 として用いている。 図 10の (A)〜(F)までは、横軸に時間、縦軸に速度をとつたグラフである力 どの 場合も、以下に示す動作の制約条件を満たす動きになっている。
(粘性摩擦やクーロン摩擦を除去するための動作の制約条件) 区間 [a, b]における 速度 Vfbの積分値が零である。これは、例えば、往復動作のような動きである。
(一定外乱を除去するための動作の制約条件) ある速度指令により動作させた際に 区間 [al, bl]より求めたイナーシャ J1と、正負を反転させた逆向きの指令により動作 させた際に区間 [a2, b2]より求めたイナーシャ J2を求め、 J1^J2の平均値を求める。 このように、従来の制御定数同定装置では、上で説明した動作の制約条件を満た す動作をさせたうえでイナーシャ Jを同定して ヽた。
特許文献 l :WO96Z37039号公報(第 5— 8頁、図 4および図 5)
発明の開示
発明が解決しょうとする課題
従来のイナーシャ推定装置は、上に挙げた二つの制約条件があるために、例えば 、 1回の位置決め動作等では、粘性摩擦やクーロン摩擦、重力などの一定外乱の影 響を除去できず、イナ一シャを推定できな!、と 、う問題があった。
また、制約条件を無視してイナ一シャを同定すると、イナーシャ推定誤差が大きくな るという問題を抱えていた。
本発明はこのような問題点に鑑みてなされたものであり、動作によらずに、粘性摩 擦やクーロン摩擦、重力などの一定外乱の影響がある場合も、簡単な計算で精度良 く制御対象のイナーシャ同定街— hatおよび粘性摩擦係数同定値 D— hatおよび一 定外乱同定値 C— hatを算出することができ、フィルタ処理を行うことでその後の計算 で用いる時間微分演算によるノイズの影響を少なくでき、除算 3回と時間微分演算 5 回のみの簡単な計算でイナ一シャを同定することができ、減算 1回と乗算 1回のみの 簡単な計算で計算量も少なく摩擦定数を同定でき、乗算 2回と 3変数の減算 1回の簡 単な計算で計算量も少なく一定外乱を同定でき、フィルタの分子の次数より分母の次 数を 3次以上とするためイナ一シャの同定演算で行う 3回の時間微分が可能となるモ ータ制御装置とシステム同定方法を提供することを目的とする。
また、本発明では、粘性摩擦やクーロン摩擦、重力などの一定外乱の影響がある場 合も、簡単な計算で精度良く制御対象のイナーシャ同定 »— hatおよび粘性摩擦 係数同定値 D— hatを算出することができ、フィルタ処理を行うことでその後の計算で 用いる時間微分演算によるノイズの影響を少なくでき、時間微分演算 6回と乗算 6回 と減算 3回と除算 2回のみの簡単な計算でイナ一シャと摩擦定数の両方を同定できる モータ制御装置とそのシステム同定方法を提供することを目的とする。
課題を解決するための手段
請求項 1記載の本発明は、指令発生器と、前記指令発生器により出力される指令と 電動機に付加された検出器カゝら得られる検出値をもとに制御を行い電動機を駆動す るための電流値を出力する制御器と、制御対象のイナーシャ同定街—hatおよび粘 性摩擦係数同定値 D— hatおよび一定外乱同定値 C— hatを同定する同定器を備え たモーション制御装置のシステム同定方法にぉ 、て、速度検出値 Vfbを時間微分し て加速度検出値 afbを算出し、加速度検出値 afbをフィルタ処理した信号 Fafbを算 出し、速度検出値 Vfbをフィルタ処理した信号 Fvfbを算出し、トルク指令値 Trefをフ ィルタ処理した信号 Ftrefを算出し、 Ftrefと Fvfbと Fafbをもとに時間微分および四 則演算し、制御対象のイナーシャ同定銜—hatおよび粘性摩擦係数同定値 D— hat および一定外乱同定値 C— hatを算出するという手順をとつたのである。
請求項 2記載の本発明は、請求項 1記載のモーション制御装置のシステム同定方 法において、イナ一シャの同定 ffi— hatを算出する処理では、 Ftref、 Fvfb, Fafb をそれぞれ時間微分して Dtref、 Dvfb、 Dafbを算出し、 Dtrefを Dvrefで除算して W 1を算出し、 Dafbを前記 Dvrefで除算して W2を算出し、 W1を時間微分して DW1を 算出し、 W2を時間微分して DW2を算出し、 DW1を前記 DW2で除算した値をイナ ーシャ同定銜 _hatとすると 、う手順をとつたのである。
請求項 3記載の本発明は、請求項 2記載のモーション制御装置のシステム同定方 法において、粘性摩擦係数同定値 D— hatを算出する処理では、イナーシャ同定値 J —hatを算出した後に、 W2とイナーシャ同定 ffi— hatを乗算して W3を算出し、 W1 から W3を減算したものを粘性摩擦係数の同定値 D— hatとするという手順をとつたの である。
請求項 4記載の本発明は、請求項 3記載のモーション制御装置のシステム同定方 法において、一定外乱同定値 C— hatを算出する処理では、イナーシャ同定 ffi— ha tと粘性摩擦係数同定値 D— hatを算出した後に、 Fvfbと前記粘性摩擦同定値 D— hatを乗算して TFを算出し、 Fafbとイナーシャ同定 ffi— hatを乗算して TAを算出し 、 Ftrefから TFと TAを減算したものを一定外乱の同定値 C— hatとすると ヽぅ手順を とったのである。
請求項 5記載の本発明は、請求項 1記載のモーション制御装置のシステム同定方 法において、フィルタ処理で用いるフィルタは、分子の次数に比べて分母の次数が 少なくとも 3次以上大きいフィルタにしたのである。
請求項 6記載の本発明は、請求項 1乃至請求項 5記載のモーション制御装置のシ ステム同定方法において、内部の状態量を監視しておき、状態量が設定値以上にな つた時にのみ、制御対象のイナーシャ同定街— hatおよび粘性摩擦係数同定値 D —hatおよび一定外乱同定値 C— hatを算出する処理を行い、状態量が設定値以上 でない場合には、 J— hat、 D— hat、 C— hatの同定値として、それぞれ前回の値を 保持しておくようにしたのである。
請求項 7記載の本発明は、請求項 6記載のモーション制御装置のシステム同定方 法において、監視しておく状態量は、速度 Vfb、加速度 afb、トルク指令 Trefあるいは それらの信号をフィルタ処理した Fvfb、 Fafb, Ftref、あるいは Dvfb、 Dw2のいずれ か一つ、あるいは、それらの組み合わせであるようにしたのである。
請求項 8記載の本発明は、指令発生器と、指令発生器により出力される指令と電動 機に付加された検出器カゝら得られる検出値をもとに位置制御を行う位置制御器と、速 度制御を行う速度制御器と、電動機を駆動するための電流を出力する電流制御器と 、制御対象のイナーシャ同定街— hatおよび粘性摩擦係数同定値 D— hatおよび一 定外乱同定値 C— hatを同定する同定器を備えたモーション制御装置において、同 定器は速度検出値 Vfbを時間微分して加速度検出値 afbを算出する時間微分器と、 加速度検出値 afbをフィルタ処理して Fafbを算出する第 1のフィルタと、速度検出値 Vfbをフィルタ処理して Fvfbを算出する第 2のフィルタと、トルク指令値 Trefをフィル タ処理して Ftrefを算出する第 3のフィルタと、 Ftrefと Fvfbと Fafbをもとに時間微分 および四則演算し、制御対象のイナーシャ同定街—hatおよび粘性摩擦係数同定 値 D— hatおよび一定外乱同定値 C— hatを算出する JDC推定器を備えるようにした のである。
請求項 9記載の本発明は、請求項 8記載のモーション制御装置において、 JDC推 定器は、イナ一シャを推定するイナーシャ推定器と、粘性摩擦を推定する粘性摩擦 推定器と、
一定外乱を推定する一定外乱推定器とを備えるようにしたのである。
請求項 10記載の本発明は、指令発生器と、前記指令発生器により出力される指令 と電動機に付加された検出器力も得られる検出値をもとに電動機を駆動するための 電流値を出力する制御器と、制御対象のイナーシャ同定街—hatおよび粘性摩擦 係数同定値 D— hatおよび一定外乱同定値 C— hatを同定する同定器を備えたモー シヨン制御装置のシステム同定方法にぉ 、て、速度検出値 Vfbを時間微分して加速 度検出値 afbを算出し、前記加速度検出値 afbをフィルタ処理した信号 Fafbを算出し 、前記速度検出値 Vfbをフィルタ処理した信号 Fvfbを算出し、前記トルク指令値 Tre fをフィルタ処理した信号 Ftrefを算出し、前記 Ftref、 Fvfb, Fafbをそれぞれ時間微 分した Dtref、 Dvfb、 Dafbおよび、もう一度時間微分した D2tref、 D2vfb、 D2afbを 算出し、以下の式により、
J_hat = (D2tref · Dvfb Dtref · D2vfb) /
(Dvfb - D2afb - D2vfb - Dafb) (a)
D hat = (D2afb · Dtref Dafb · D2tref) /
(Dvfb - D2afb - D2vfb - Dafb) (b)
イナーシャ同定銜—hatおよび粘性摩擦係数同定値 D— hatを算出するという手 順をとつたのである。
請求項 11記載の本発明は、請求項 10記載のモーション制御装置のシステム同定 方法において、前記フィルタ処理で用いるフィルタは、分子の次数に比べて分母の 次数が少なくとも 3次以上大きいフィルタであることを特徴とするものである。
請求項 12記載の本発明は、請求項 10記載のモーション制御装置のシステム同定 方法において、前記フィルタ処理で用いるフィルタは、複数の移動平均フィルタを直 列につなげたものであることを特徴とするものである。 請求項 13記載の本発明は、モーション制御装置の内部の信号を監視しておき、前 記信号が設定値以上になった時にのみ、前記制御対象のイナーシャ同定街— hat および粘性摩擦係数同定値 D— hatおよび一定外乱同定値 C— hatを算出する処理 を行い、前記状態量が設定値以上でない場合には、 J—hat、 D— hat、 C— hatの同 定値として、それぞれ前回の値を保持しておくことを特徴とするものである。
請求項 14記載の本発明は、前記監視しておく信号は、速度 Vfb、加速度 afb、トル ク指令 Trefあるいはそれらの信号をフィルタ処理した Fvfb、 Fafb、 Ftref、あるいは、 それらの信号を時間微分した Dtref、 Dvfb、 Dafb、 D2tref、 D2vfb、 D2afbのいず れか一つ、あるいは、それらの組み合わせであることを特徴とするものである。
請求項 15記載の本発明は、指令を生成する指令発生器と、前記指令と電動機に 付加された検出器から得られる検出値をもとに電動機を駆動するための電流を生成 する電流制御器と、制御対象のイナーシャ同定街—hatおよび粘性摩擦係数同定 値 D— hatおよび一定外乱同定値 C— hatを同定する同定器を備えたモーション制 御装置において、速度検出値 Vfbを時間微分して加速度検出値 afbを算出する加速 度算出器と、前記加速度検出値 afbをフィルタ処理した信号 Fafbを算出する第 1のフ ィルタと、前記速度検出値 Vfbをフィルタ処理した信号 Fvfbを算出する第 2のフィル タと、前記トルク指令値 Trefをフィルタ処理した信号 Ftrefを算出する第 3のフィルタ と、前記 Ftref、 Fvfb, Fafbをそれぞれ時間微分した Dtref、 Dvfb、 Dafbおよび、も う一度時間微分した D2tref、 D2vfb、 D2afbを算出する微分器と、以下の式により、
J_hat = (D2tref · Dvfb Dtref · D2vfb) /
(Dvfb - D2afb - D2vfb - Dafb) (a)
D hat = (D2afb · Dtref Dafb · D2tref) /
(Dvfb - D2afb - D2vfb - Dafb) (b)
イナーシャ同定街—hatおよび粘性摩擦係数同定値 D— hatを算出することを特 徴とするちのである。
発明の効果
請求項 1記載の本発明によると、トルク指令と粘性摩擦と一定外乱の関係を近似す ることなく変形 (微分や除算を使用)した式を用いて、速度、加速度、トルク指令値を フィルタ処理した後の信号を用いて、それらを時間微分と四則演算するだけで同定 計算を行うため、動作によらずに、粘性摩擦やクーロン摩擦、重力などの一定外乱の 影響がある場合も、精度良く制御対象のイナーシャ同定街—hatおよび粘性摩擦係 数同定値 D— hatおよび一定外乱同定値 C— hatを精度良く同定することができる。 また、フィルタ処理を行うことでその後の計算で用いる時間微分演算によるノイズの影 響を少なくできる。
請求項 2記載の本発明によると、除算 3回と時間微分演算 5回のみでイナ一シャを 同定できるため、簡単な計算で計算量も少なく同定することができる。
請求項 3記載の本発明によると、減算 1回と乗算 1回のみで摩擦定数を同定できる ため、簡単な計算で計算量も少なく同定することができる。
請求項 4記載の本発明によると、乗算 2回と 3変数の減算 1回のみで一定外乱を同 定できるため、簡単な計算で計算量も少なく同定することができる。
請求項 5記載の本発明によると、フィルタの分子の次数より分母の次数を 3次以上と するため、イナ一シャの同定演算で行う 3回の時間微分が可能になるようにすること ができる。
請求項 6および請求項 7記載の本発明によると、監視している状態量が設定値以下 の場合は、同定計算を行わないため、微小な値で除算することによって精度が悪くな ることを防ぐことができ、どのような場合も同定精度を保つことができる。
請求項 8および請求項 9記載の本発明によると、動作によらずに、粘性摩擦やクー ロン摩擦、重力などの一定外乱の影響がある場合も、精度良く制御対象のイナーシ ャ同定街— hatおよび粘性摩擦係数同定値 D— hatおよび一定外乱同定値 C— ha tを精度良く同定することができるモーション制御装置を提供できる。
請求項 10記載の本発明によると、トルク指令と粘性摩擦と一定外乱の関係を近似 することなく変形した式を両辺時間微分した式と、さらにその式の両辺を時間微分し た式の 2つを連立方程式とみなして解を求め、速度、加速度、トルク指令値をフィルタ 処理した後の信号を用いて、それらを時間微分と四則演算するだけで同定計算を行 うため、動作によらずに、粘性摩擦やクーロン摩擦、重力などの一定外乱の影響があ る場合も、精度良く制御対象のイナーシャ同定街—hatおよび粘性摩擦係数同定値 D— hatおよび一定外乱同定値 C— hatを精度良く同定することができる。また、フィ ルタ処理を行うことでその後の計算で用いる時間微分演算によるノイズの影響を少な くでさる。
請求項 11記載の本発明によると、フィルタの分子の次数より分母の次数を 3次以上 とするため、ノイズの影響を受けずにイナ一シャの同定演算で行う 3回の時間微分が 可能になるようにすることができる。
請求項 12記載の本発明によると、フィルタとして移動平均を用いることで、各変数 の位相遅れの差を無くすことができるため、さらに精度良く同定することができるよう にすることができる。
請求項 13および請求項 14記載の本発明によると、監視している信号が設定値以 下の場合は、同定計算を行わないため、微小な値で除算することによって精度が悪く なることを防ぐことができ、どのような場合も同定精度を保つことができる。
図面の簡単な説明
[0006] [図 1]本発明の方法を適用するモーション制御装置の構成を示すブロック図
[図 2]本発明の方法を適用する第 2のモーション制御装置の構成を示すブロック図 [図 3]本発明の第 1の方法を適用するモーション制御装置の同定器の構成を示すブ ロック図
[図 4]イナーシャ同定部の構成を示すブロック図
[図 5]粘性摩擦同定部の構成を示すブロック図
[図 6]—定外乱同定部の構成を示すブロック図
[図 7]本発明の第 2の方法を適用するモーション制御装置の同定器の構成を示すブ ロック図
[図 8]イナーシャ '粘性摩擦同定部の構成を示すブロック図
[図 9]従来の方法を適用するモーション制御装置の構成を示すブロック図
[図 10]従来の方法で限定されている動作を示す図
符号の説明
[0007] 1 指令発生器
2 制御器 電動機
機械
検出器
同定器
位置制御器
速度制御器
電流制御器
0 時間微分器
1 第 1のフィルタ
2 第 2のフィルタ
3 第 3のフィルタ
4 JDC推定器
5 イナーシャ同定器
6 粘性摩擦同定器
7 一定外乱同定器
8 イナ一シャ'粘性摩擦同定器1, 42, 43 時間微分器
4, 45 除算器
6, 47 時間微分器
8 除算器
9, 4A、4B 時間微分器
C, 4D, 4E, 4F, 4G, 4H 乗算器1, J, 4K 減算器
L, 4M 除算器
1 乗算器
2 減算器
1, 62 乗算器
3 減算器 71 指令発生器
72 速度制御器
73 推定器
74 同定器
75 電流制御器
100 モーション制御装置
発明を実施するための最良の形態
[0008] 以下、本発明の方法の具体的実施例について、図に基づいて説明する。
実施例 1
[0009] 図 1は、本発明の方法を実施するモーション制御装置の構成を示すブロック図であ る。図において 1は位置指令 Xrefを発生する指令発生器となっている。また 2は、制 御器を表し、位置指令と位置検出値 Xfbおよび速度検出値 Vfbをもとに制御演算を 行い電流 Iを出力する。制御器 2内部の演算はどのようなものでも良いが、本実施例 では、位置指令 Xrefと位置検出値 Xfbをもとに速度指令 Vrefを出力する位置制御 器 7と速度指令 Vrefと速度検出値 Vfbをもとにトルク指令値 Trefを出力する速度制 御器 8とトルク指令値通りに電流が流れるように制御を行う電流制御器 9から構成され ている。ここで、速度検出値は位置検出値を時間微分したものでも良ぐデジタル制 御時には、時間微分演算として、差分後制御周期で除算した近似微分を用いてもよ い。 3は電動機であり、 4の機械は電動機 3に結合されている。 5は電動機 3の位置お よび速度を検出する検出器を表す。 6は同定器となっており、トルク指令 Trefおよび 速度検出値 Vfbをもとに制御対象のイナーシャ同定街—hatおよび粘性摩擦係数 同定値 D— hatおよび一定外乱同定値 C— hatを算出する。
[0010] 図 3は 6の同定器内部の処理を示すブロック図である。図において、 10は時間微分 器を表し、速度検出値 Vfbを時間微分し、加速度検出値 afbを計算する。ここで、デ ジタル制御時には時間微分演算として、今回と前回の信号の差分後に制御周期で 除算する近似微分を用いてもょ ヽ。
11、 12、 13はそれぞれ同一のフィルタを示す。ここで用いるフィルタは、分子の次数 より分母の次数が 3次以上のものが良ぐ例えばフィルタの伝達特性が式( 1 )の伝達 関数 Gfilで表されるようなものを用いればょ 、。
Figure imgf000013_0001
ここで、 sはラプラス演算子を表すものとする。
また、用いるフィルタはどのようなものでも良ぐ一次のローパスフィルタを 3つ以上 直列にして使用しても全く問題な 、。
[0011] 14ίお DC推定器を表し、トルク指令 Tref、速度検出値 Vfb、加速度検出値 afbをフ ィルタ処理した Ftref、 Fvfb、 Fafbを入力し、制御対象のイナーシャ同定 ffi— hatお よび粘性摩擦係数同定値 D— hatおよび一定外乱同定値 C— hatを算出する。さら に 14の内部は、 15のイナーシャ同定器、 16の粘性摩擦同定器、 17の一定外乱同 定器の 3つ力 構成されて 、る。
以下に、 15のイナーシャ同定器、 16の粘性摩擦同定器、 17の一定外乱同定器の内 部の処理を詳しく説明する。
[0012] (I)イナーシャ同定器内部の処理
図 4はイナーシャ同定器 15の処理を示すブロック図である。図において、 41, 42, 43, 46、 47は時間微分器を表し、 44, 45, 48は除算器を表す。イナーシャ同定街 —hatは以下の手順で演算する。
まず、 Ftref、 Fvfb、 Fafbをそれぞれ時間微分して Dtref、 Dvfb、 Dafbを算出し、 次に、 Dtrefを前記 Dvfbで除算して W1を算出し、次に、 Dafbを前記 Dvfbで除算し て W2を算出し、次に、 W1を時間微分して DW1を算出し、次に、 W2を時間微分し て DW2を算出し、次に、 DW1を前記 DW2で除算した値をイナーシャ同定街— hat とする。
この処理を数式で表すと式(2)力 式(9)のようになる。
Dtref =dFtref/dt (2)
Dvfb = dFvfb/dt (3)
Dafb = dFafb/dt (4)
Wl = Dtref /Dvfb (5)
W2 = Dafb/Dvfb (6)
DWl = dWl/dt (7) DW2 = dW2/dt (8)
J_hat = DWl/DW2 (9)
ここで、デジタル制御時には時間微分演算として、今回と前回の信号の差分後に制 御周期で除算する近似微分を用いてもょ ヽ。
[0014] (Π)粘性摩擦同定器内部の処理
図 5は粘性摩擦同定器 16の処理を示すブロック図である。図において 51は乗算器 を表し、 52は減算器を表す。粘性摩擦同定値 D— hatは以下の手順で演算する。 まず、 W2と前記イナーシャ同定街— hatを乗算して W3を算出し、次に、 W1から 前記 W3を減算したものを粘性摩擦係数の同定値 D— hatとする。この処理を数式で 表すと式(10)力 式(11)のようになる。
[0015] W3= W2-J_hat (10)
D_hat=Wl -W3 (11)
[0016] (III)一定外乱同定器内部の処理
図 6は一定外乱同定器 17の処理を示すブロック図である。図において、 61、 62は 乗算器を示し、 63は 1変数カゝら 2変数を減算する減算器を示す。一定外乱同定値 D —hatは以下の手順で演算する。
まず、 Fvfbと前記粘性摩擦同定値 D— hatを乗算して TFを算出し、次に、 Fafbと 前記イナーシャ同定銜—hatを乗算して TAを算出し、次に、 Ftrefから前記 TFと前 記 TAを減算したものを一定外乱の同定値 C_hatとする。
この処理を数式で表すと式( 12)力 式( 14)のようになる。
[0017] TF=Fvfb -D_hat (12)
TA= Fafb 'J一 hat (13)
C_hat = Ftref - TF - TA (14)
[0018] このように、以上の処理を同定器 6の中で行うことで、制御対象のイナーシャ同定値 J— hatおよび粘性摩擦係数同定値 D— hatおよび一定外乱同定値 C— hatを同定 することができる。
[0019] また、以上の処理を常時演算するのではなぐ速度 Vfbあるいは加速度 afbあるい はトルク指令値 Trefのうちの一つある 、はその組み合わせを監視しておき、それらの 変数が予め設定した閾値以上の時のみ、上記処理を実施することにより、同定精度 を上げることができる。
また、監視する変数は、 Vfb、 afb、 Trefをフィルタ処理した、 Fvfb、 Fafb、 Ftrefの うちの一つある 、はその組み合わせでもよ 、。
また、監視している変数が閾値よりも小さい時は演算を行わないだけでなぐ演算を 行っていた最後の値を、その同定値の値として保持しておくことにより、演算している 時だけ各同定値が更新されるようにすることができる。
また、監視して!/ヽる変数が閾値よりも小さ!、時に制御対象のイナーシャ同定街— h atおよび粘性摩擦係数同定値 D— hatおよび一定外乱同定値 C— hat全ての演算を 行わな 、ようにするのではなぐそのうちの一部だけ演算を行わな ヽと 、うようにしても よい。
[0020] 本願発明は、トルク指令の変わりに電流指令を用いても良ぐその場合はトルク定 数 (推力定数)を用いて、トルクと電流の単位を変換すればょ 、。
また、ここでは、回転機を例にとって説明したが、並進運動をする機械の場合は、同 様の計算で、イナーシャ Jではなく質量 Mを推定することができる。
[0021] また、本実施例 1では、位置制御を行う構成を述べたが、図 2に示すように、速度制 御を行う場合でも全く同様の方法で、制御対象のイナーシャ同定街—hatおよび粘 性摩擦係数同定値 D— hatおよび一定外乱同定値 C— hatを同定することができる。 この場合は、 1の指令発生器は速度指令 Vrefを発生し、 2の制御器内部では、速度 制御器と電流制御器の演算のみが行われることになる。
[0022] 実施例 1の方式の導出方法を示す。
粘性摩擦 (摩擦係数 D)と一定外乱 Cがイナーシャ Jの機械に作用する場合、そこに 与えられるトルク指令 Trefと、速度 Vfbの関係は式(15)のようになる。
Tref=J-afb + D-Vfb + C (15)
両辺フィルタ処理すると式(16)が得られる。
Ftref = J · Faf b + D · Fvf b + C (16)
両辺時間微分処理すると一定値 Cが消去でき、式(17)が得られる。
Dtref = J - Daf b + D - Dvf b (17) 両辺 Dvfbで除算すると式(18)が得られる。
Dtref /Dvfb =J - Dafb/Dvfb + D (18)
Wl ( = Dtref/Dvfb)、 W2 ( = Dafb/Dvfb)を使って書きなおすと式(19)が得 られる。
W1 =J-W2 + D (19)
両辺時間微分処理すると一定値 Dが消去でき、式(20)が得られる。
DW1 =J-DW2 (20)
よって Jは式(21)で計算できる。
J = DW1/DW2 (21)
式(19)および式(21)より Dは式(22)のように計算できる。
D=W1 -J-W2 (22)
式(16)より Cは式 (23)のように計算できる。
C = Ftref - J - Faf b - D - Fvf b (23)
以上のように、近似を一切使用せず、各値が同定できる式が導出できるのである。 実施例 2
[0023] 次に本発明の第 2の方法について説明する。
図 1は、本発明の方法を実施するモーション制御装置の構成を示すブロック図であ る力 この図については実施例 1で説明したので説明を省略する。
[0024] 図 7は同定器 6内部の処理を示すブロック図であり、図 3のものとは異なっているの で、以下に説明する。図において、 10は時間微分器を表し、速度検出値 Vfbを時間 微分し、加速度検出値 afbを計算する。ここで、デジタル制御時には時間微分演算と して、今回と前回の信号の差分後に制御周期で除算する近似微分を用いてもょ 、。 11、 12、 13はそれぞれ第 1のフィルタ、第 2のフィルタ、第 3のフィルタを示す。ここで 用いるフィルタは、分子の次数より分母の次数が 3次以上のものが良ぐ例えばフィル タの伝達特性が式(24)の伝達関数 Gfilで表されるようなものを用いればょ 、。
Figure imgf000016_0001
ここで、 sはラプラス演算子を表すものとする。
また、用いるフィルタはどのようなものでも良ぐ一次のローパスフィルタを 3つ以上 直列にして使用しても全く問題ない。また、移動平均フィルタを 3つ以上直列にして 使用しても全く問題ない。この場合、フィルタ処理による位相の遅れがないため、さら に推定精度を上げることができる。
[0025] 14〖お DC推定器を表し、トルク指令 Tref、速度検出値 Vfb、加速度検出値 afbをフ ィルタ処理した Ftref、 Fvfb、 Fafbを入力し、制御対象のイナーシャ同定 ffi— hatお よび粘性摩擦係数同定値 D— hatおよび一定外乱同定値 C— hatを算出する。さら に 14の JDC推定器の内部は、 18のイナーシャ '粘性摩擦同定器、 16の一定外乱同 定器の 2つ力 構成されて 、る。
以下に、 18のイナーシャ '粘性摩擦同定器、 16の一定外乱同定器の内部の処理を 詳しく説明する。
[0026] (IV)イナーシャ '粘性摩擦同定器内部の処理
図 8はイナーシ 粘性摩擦同定器 15の処理を示すブロック図である。図において 、 41, 42, 43, 49, 4A, 4Bは時間微分器を表し、 4C, 4D, 4E, 4F, 4G, 4Hは乗 算器を表す。 41, J, 4Kは減算器を表し、 4L, 4Mは除算器を表す。イナ一シャ同 定銜—hatおよび粘性摩擦係数同定値 D— hatは以下の手順で演算する。
[0027] まず、 Ftref、 Fvfb、 Fafbをそれぞれ時間微分して Dtref、 Dvfb、 Dafbを算出し、 次に、もう一度時間微分し, D2tref、 D2vfb、 D2afbを算出する。ここで、デジタル制 御時には時間微分演算として、今回と前回の信号の差分後に制御周期で除算する 近似微分を用いてもよい。
そして、これらの変数を用いて式(25)と(26)の演算を行い、イナーシャ同定街— hatおよび粘性摩擦係数同定値 D_hatを推定する。
J_hat = (D2tref · Dvfb Dtref · D2vfb) /
(Dvfb - D2af b - D2vf b - Dafb) (25)
D— hat = (D2afb · Dtref -Dafb · D2tref) /
(Dvfb - D2af b - D2vf b - Dafb) (26)
[0028] 図 6は一定外乱同定器 17の処理を示すブロック図である。一定外乱同定器内部の 処理につ ヽては、実施例 1で説明したので説明を省略する。
[0029] また、以上の処理を常時演算するのではなぐ速度 Vfbあるいは加速度 afbあるい はトルク指令値 Trefのうちの一つある 、はその組み合わせを監視しておき、それらの 変数が予め設定した閾値以上の時のみ、上記処理を実施することにより、同定精度 を上げることができる。
また、監視する変数は、速度 Vfb、加速度 afb、トルク指令 Trefあるいはそれらの信 号をフィルタ処理した Fvfb、 Fafb、 Ftref、あるいは、それらの信号を時間微分した D tref, Dvfb、 Dafb、 D2tref、 D2vfb、 D2afbのいずれか一つか、その糸且み合わせで ちょい。
また、監視している変数が閾値よりも小さい時は演算を行わないだけでなぐ演算を 行っていた最後の値を、その同定値の値として保持しておくことにより、演算している 時だけ各同定値が更新されるようにしてもょ 、。
また、監視して!/ヽる変数が閾値よりも小さ!、時に制御対象のイナーシャ同定街— h atおよび粘性摩擦係数同定値 D— hatおよび一定外乱同定値 C— hat全ての演算を 行わな 、ようにするのではなぐそのうちの一部だけ演算を行わな ヽと 、うようにしても よい。
[0030] 本発明は、トルク指令の変わりに電流指令を用いても良ぐその場合はトルク定数( 推力定数)を用いて、トルクと電流の単位を変換すればょ 、。
また、ここでは、回転機を例にとって説明したが、並進運動をする機械の場合は、同 様の計算で、イナーシャ Jではなく質量 Mを推定することができる。
[0031] また、本実施例 2では、位置制御を行う構成を述べたが、図 2に示すように、速度制 御を行う場合でも全く同様の方法で、制御対象のイナーシャ同定街—hatおよび粘 性摩擦係数同定値 D— hatおよび一定外乱同定値 C— hatを同定することができる。 この場合は、 1の指令発生器は速度指令 Vrefを発生し、 2の制御器内部では、速度 制御器と電流制御器の演算のみが行われることになる。
[0032] 実施例 2の方式の導出方法を示す。
粘性摩擦 (摩擦係数 D)と一定外乱 Cがイナーシャ Jの機械に作用する場合、そこに 与えられるトルク指令 Trefと、速度 Vfbの関係は式(27)のようになる。
Tref=J - afb + D-Vfb + C (27)
両辺フィルタ処理すると式(28)が得られる。 Ftref = J · Faf b + D · Fvf b + C (28)
両辺時間微分処理すると一定値 Cが消去でき、式(29)が得られる。
Dtref = J - Daf b + D - Dvf b (29)
さらに両辺時間微分処理すると式(30)が得られる。
D2tref =J - D2af b + D - D2vf b (30)
ここで、式 (29)と式 (30)を連立方程式とみなし Jと Dの推定値を計算すると式 (31)と 式(32)が得られる。
J = (D2tref · Dvfb - Dtref · D2vfb) /
(Dvf b - D2afb - D2vfb - Dafb) (31)
D= (D2afb - Dtref - Dafb - D2tref) /
(Dvfb - D2afb - D2vfb - Dafb) (32)
式(27)より Cは式 (33)のように計算できる。
C = Ftref - J - Faf b - D - Fvf b (33)
以上のように、近似を一切使用せず、各値が同定できる式が導出できるのである。 産業上の利用可能性
フィルタ処理した変数を用いて、後は、時間微分処理と四則演算のみで、精度良く 制御対象のイナーシャ同定街— hatおよび粘性摩擦係数同定値 D— hatおよび一 定外乱同定値 C— hatを同定することができて、産業用ロボットや工作機のシステム 同定と 、う用途にも適用できる。

Claims

請求の範囲
[1] 指令発生器と、前記指令発生器により出力される指令と電動機に付加された検出 器力 得られる検出値をもとに電動機を駆動するための電流値を出力する制御器と
、制御対象のイナーシャ同定街— hatおよび粘性摩擦係数同定値 D— hatおよび一 定外乱同定値 C— hatを同定する同定器を備えたモーション制御装置のシステム同 定方法において、
速度検出値 Vfbを時間微分して加速度検出値 afbを算出し、
前記加速度検出値 afbをフィルタ処理した信号 Fafbを算出し、
前記速度検出値 Vfbをフィルタ処理した信号 Fvfbを算出し、
前記トルク指令値 Trefをフィルタ処理した信号 Ftrefを算出し、
前記 Ftrefと Fvfbと Fafbをもとに時間微分および四則演算し、
前記制御対象のイナーシャ同定街—hatおよび粘性摩擦係数同定値 D— hatおよ び一定外乱同定値 C_hatを算出する、
という手順で処理することを特徴とするモーション制御装置のシステム同定方法。
[2] 前記イナ一シャの同定街— hatを算出する処理では、
前記 Ftref、 Fvfb, Fafbをそれぞれ時間微分して Dtref、 Dvfb、 Dafbを算出し、 前記 Dtrefを前記 Dvfbで除算して W1を算出し、
前記 Dafbを前記 Dvfbで除算して W2を算出し、
前記 W1を時間微分して DW1を算出し、
前記 W2を時間微分して DW2を算出し、
前記 DW1を前記 DW2で除算した値をイナーシャ同定街—hatとすること特徴とす る請求項 1記載のモーション制御装置のシステム同定方法。
[3] 前記粘性摩擦係数同定値 D— hatを算出する処理では、
前記イナーシャ同定街— hatを算出した後に前記 W2と前記イナーシャ同定 ffi— hatを乗算して W3を算出し、
前記 W1から前記 W3を減算したものを粘性摩擦係数の同定値 D— hatとすることを 特徴とする請求項 1記載のモーション制御装置のシステム同定方法。
[4] 前記一定外乱同定値 C— hatを算出する処理では、 前記イナーシャ同定銜—hatと前記粘性摩擦係数同定値 D— hatを算出した後に
、前記 Fvfbと前記粘性摩擦同定値 D— hatを乗算して TFを算出し、
前記 Fafbと前記イナーシャ同定 ffi— hatを乗算して TAを算出し、
前記 Ftrefから前記 TFと前記 TAを減算したものを一定外乱の同定値 C— hatとす ることを特徴とする請求項 1記載のモーション制御装置のシステム同定方法。
[5] 前記フィルタ処理で用いるフィルタは、分子の次数に比べて分母の次数が少なくと も 3次以上大きいフィルタであることを特徴とする請求項 1記載のモーション制御装置 のシステム同定方法。
[6] 内部の状態量を監視しておき、前記状態量が設定値以上になった時にのみ、前記 制御対象のイナーシャ同定街— hatおよび粘性摩擦係数同定値 D— hatおよび一 定外乱同定値 C— hatを算出する処理を行い、前記状態量が設定値以上でない場 合には、 J— hat、 D— hat、 C— hatの同定値として、それぞれ前回の値を保持してお くことを特徴とする請求項 1乃至請求項 5記載のモーション制御装置のシステム同定 方法。
[7] 前記監視しておく状態量は、速度 Vfb、加速度 afb、トルク指令 Trefある ヽはそれら の信号をフィルタ処理した Fvfb、 Fafb、 Ftref、あるいは前記 Dvfb、 Dw2のいずれ か一つ、あるいは、それらの組み合わせであることを特徴とする請求項 6記載のモー シヨン制御装置のシステム同定方法。
[8] 指令発生器と、前記指令発生器により出力される指令と電動機に付加された検出 器力 得られる検出値をもとに制御を行い電動機を駆動するための電流値を出力す る電流制御器と、制御対象のイナーシャ同定街—hatおよび粘性摩擦係数同定値 D— hatおよび一定外乱同定値 C— hatを同定する同定器を備えたモーション制御 装置において、
前記同定器は
速度検出値 Vfbを時間微分して加速度検出値 afbを算出する時間微分器と、 前記加速度検出値 afbをフィルタ処理して Fafbを算出する第 1のフィルタと、 前記速度検出値 Vfbをフィルタ処理して Fvfbを算出する第 2のフィルタと、 前記トルク指令値 Trefをフィルタ処理して Ftrefを算出する第 3のフィルタと、 前記 Ftrefと Fvfbと Fafbをもとに時間微分および四則演算し、前記制御対象のィ ナーンャ同定街—hatおよび粘性摩擦係数同定値 D— hatおよび一定外乱同定値 C_hatを算出する JDC推定器とで構成されることを特徴とするモーション制御装置。
[9] 前 ¾JDC推定器は、
イナ一シャを推定するイナーシャ推定器と、
粘性摩擦を推定する粘性摩擦推定器と、
一定外乱を推定する一定外乱推定器とを備えることを特徴とする請求項 8記載のモ ーシヨン制御装置。
[10] 指令を生成する指令発生器と、前記指令と電動機に付加された検出器から得られ る検出値をもとに電動機を駆動するための電流を生成する電流制御器と、制御対象 のイナーシャ同定 ffi— hatおよび粘性摩擦係数同定値 D— hatおよび一定外乱同 定値 C— hatを同定する同定器を備えたモーション制御装置のシステム同定方法に おいて、
速度検出値 Vfbを時間微分して加速度検出値 afbを算出し、
前記加速度検出値 afbをフィルタ処理した信号 Fafbを算出し、
前記速度検出値 Vfbをフィルタ処理した信号 Fvfbを算出し、
前記トルク指令値 Trefをフィルタ処理した信号 Ftrefを算出し、
前記 Ftref、 Fvfb, Fafbをそれぞれ時間微分した Dtref、 Dvfb、 Dafbおよび、もう 一度時間微分した D2tref、 D2vfb、 D2afbを算出し、
以下の式により、
J_hat = (D2tref · Dvfb Dtref · D2vfb) /
(Dvfb - D2afb - D2vfb - Dafb) (a)
D hat = (D2afb · Dtref Dafb · D2tref) /
(Dvfb - D2afb - D2vfb - Dafb) (b)
イナーシャ同定 iBj_hatおよび粘性摩擦係数同定値 D_hatを算出するという手順 で処理することを特徴とするモーション制御装置のシステム同定方法。
[11] 前記フィルタ処理で用いるフィルタは、分子の次数に比べて分母の次数が少なくと も 3次以上大きいフィルタであることを特徴とする請求項 10記載のモーション制御装 置のシステム同定方法。
[12] 前記フィルタ処理で用いるフィルタは、複数の移動平均フィルタを直列につなげた ものであることを特徴とする請求項 10記載のモーション制御装置のシステム同定方 法。
[13] モーション制御装置の内部の信号を監視しておき、前記信号が設定値以上になつ た時にのみ、前記制御対象のイナーシャ同定街—hatおよび粘性摩擦係数同定値 D— hatおよび一定外乱同定値 C— hatを算出する処理を行い、前記状態量が設定 値以上でない場合には、 J— hat、 D— hat、 C— hatの同定値として、それぞれ前回 の値を保持しておくことを特徴とする請求項 10乃至請求項 12記載のモーション制御 装置のシステム同定方法。
[14] 前記監視しておく信号は、速度 Vfb、加速度 afb、トルク指令 Trefある 、はそれらの 信号をフィルタ処理した Fvfb、 Fafb、 Ftref、あるいは、それらの信号を時間微分し た Dtref、 Dvfb、 Dafb、 D2tref、 D2vfb、 D2afbのいずれか一つ、あるいは、それら の組み合わせであることを特徴とする請求項 13記載のモーション制御装置のシステ ム同定方法。
[15] 指令を生成する指令発生器と、前記指令と電動機に付加された検出器から得られ る検出値をもとに電動機を駆動するための電流を生成する電流制御器と、制御対象 のイナーシャ同定 ffi— hatおよび粘性摩擦係数同定値 D— hatおよび一定外乱同 定値 C— hatを同定する同定器を備えたモーション制御装置において、
速度検出値 Vfbを時間微分して加速度検出値 afbを算出する加速度算出器と、 前記加速度検出値 afbをフィルタ処理した信号 Fafbを算出する第 1のフィルタと、 前記速度検出値 Vfbをフィルタ処理した信号 Fvfbを算出する第 2のフィルタと、 前記トルク指令値 Trefをフィルタ処理した信号 Ftrefを算出する第 3のフィルタと、 前記 Ftref、 Fvfb, Fafbをそれぞれ時間微分した Dtref、 Dvfb、 Dafbおよび、もう 一度時間微分した D2tref、 D2vfb、 D2afbを算出する微分器と、
を備え、以下の式により、
J_hat = (D2tref · Dvfb Dtref · D2vfb) /
(Dvfb - D2af b - D2vf b - Dafb) (a) D hat = (D2afb · Dtref Dafb . D2tref) /
(Dvfb - D2af b - D2vf b - Daf b) (b)
イナーシャ同定街 _hatおよび粘性摩擦係数同定値 D_hatを算出することを特徴 とするモーション制御装置。
PCT/JP2005/017592 2004-11-12 2005-09-26 モーション制御装置とシステム同定方法 WO2006051651A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2006544794A JP4919157B2 (ja) 2004-11-12 2005-09-26 モーション制御装置とシステム同定方法
US11/797,945 US7626351B2 (en) 2004-11-12 2007-05-09 Motion controller and system identifying method

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2004328574 2004-11-12
JP2004-328574 2004-11-12
JP2005019283 2005-01-27
JP2005-019283 2005-01-27

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US11/797,945 Continuation-In-Part US7626351B2 (en) 2004-11-12 2007-05-09 Motion controller and system identifying method

Publications (1)

Publication Number Publication Date
WO2006051651A1 true WO2006051651A1 (ja) 2006-05-18

Family

ID=36336334

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2005/017592 WO2006051651A1 (ja) 2004-11-12 2005-09-26 モーション制御装置とシステム同定方法

Country Status (3)

Country Link
US (1) US7626351B2 (ja)
JP (1) JP4919157B2 (ja)
WO (1) WO2006051651A1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4501117B2 (ja) * 2005-08-10 2010-07-14 株式会社安川電機 システム同定装置
JP4291344B2 (ja) * 2006-08-31 2009-07-08 ファナック株式会社 産業用ロボット
US8232758B2 (en) * 2009-08-28 2012-07-31 Fanuc Ltd Controller of electric motor having function of estimating inertia and friction simultaneously
EP3183818A1 (en) 2014-08-20 2017-06-28 Wright State University Fractional scaling digital signal processing

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1994014234A1 (en) * 1992-12-16 1994-06-23 Kabushiki Kaisha Yaskawa Denki Method and apparatus for determining constants of functions expressing characteristics of controlled system
WO1996037039A1 (fr) * 1995-05-17 1996-11-21 Kabushiki Kaisha Yaskawa Denki Appareil de determination de constantes de commande
JP2001238477A (ja) * 2000-02-25 2001-08-31 Yaskawa Electric Corp 制御定数調整装置

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0683403A (ja) * 1992-07-17 1994-03-25 Fanuc Ltd 適応pi制御方式
JPH07110717A (ja) * 1993-08-19 1995-04-25 Fanuc Ltd モータの制御方式
JPH0793003A (ja) * 1993-09-27 1995-04-07 Mitsubishi Heavy Ind Ltd 電動機の制御装置
JP3166446B2 (ja) * 1993-10-26 2001-05-14 株式会社明電舎 速度推定オブザーバ
GB9714696D0 (en) * 1997-07-11 1997-09-17 Switched Reluctance Drives Ltd Exercise apparatus
FI112734B (fi) * 1998-05-20 2003-12-31 Abb Oy Menetelmä ja sovitelma kuormitusmuutosten adaptiiviseksi kompensoinniksi

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1994014234A1 (en) * 1992-12-16 1994-06-23 Kabushiki Kaisha Yaskawa Denki Method and apparatus for determining constants of functions expressing characteristics of controlled system
WO1996037039A1 (fr) * 1995-05-17 1996-11-21 Kabushiki Kaisha Yaskawa Denki Appareil de determination de constantes de commande
JP2001238477A (ja) * 2000-02-25 2001-08-31 Yaskawa Electric Corp 制御定数調整装置

Also Published As

Publication number Publication date
US20070216333A1 (en) 2007-09-20
US7626351B2 (en) 2009-12-01
JPWO2006051651A1 (ja) 2008-05-29
JP4919157B2 (ja) 2012-04-18

Similar Documents

Publication Publication Date Title
CN102969971B (zh) 电动机控制装置
CN102904521B (zh) 具有同时推定惯量、摩擦和弹性的功能的电动机控制装置
CA3021274C (en) Control method for electric vehicle and control device for electric vehicle
JP5591400B2 (ja) 駆動機械の負荷特性推定装置
EP2815914A1 (en) Vibration suppression control device for electric motor-driven vehicle and method for controlling vibration suppression
JP3892823B2 (ja) モータの速度制御装置
KR101220915B1 (ko) 활성화 함수와 토크 보상기를 이용한 속도 리플 억제 방법
WO2006051651A1 (ja) モーション制御装置とシステム同定方法
WO2014115626A1 (ja) 誘導モータ制御装置および誘導モータ制御方法
JP2020010454A (ja) 電動車両の制御方法、及び、制御装置
EP2597062B1 (en) Elevator controlling method, elevator controlling device, and elevator device using the same
JP2006217729A (ja) モータ制御装置とその制御方法
JP4507110B2 (ja) ディジタルサーボ制御装置
JP6586856B2 (ja) 電動車両の制御方法、及び、制御装置
JP2001051721A (ja) 衝突検出装置
JP5256704B2 (ja) 慣性モーメント推定装置
JP5407435B2 (ja) モータ制御装置
KR101417486B1 (ko) 착용로봇의 의도토크 추출방법 및 추출시스템
JP3892824B2 (ja) モータの位置制御装置
EP0413030A1 (en) Servo control method by the disturbance estimation observer
JP2023167670A (ja) 電動車両制御方法及び電動車両制御装置
JP3856215B2 (ja) 速度制御装置
WO2023171122A1 (ja) モータ制御装置、およびその自動調整方法
JP6697678B2 (ja) 誘導電動機の制御装置
JPH11313495A (ja) 電動機サーボ系の制御装置

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KM KP KR KZ LC LK LR LS LT LU LV LY MA MD MG MK MN MW MX MZ NA NG NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SM SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): BW GH GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LT LU LV MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 2006544794

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 200580035784.0

Country of ref document: CN

WWE Wipo information: entry into national phase

Ref document number: 11797945

Country of ref document: US

WWP Wipo information: published in national office

Ref document number: 11797945

Country of ref document: US

122 Ep: pct application non-entry in european phase

Ref document number: 05785280

Country of ref document: EP

Kind code of ref document: A1