JP3861810B2 - 機械定数推定装置 - Google Patents
機械定数推定装置 Download PDFInfo
- Publication number
- JP3861810B2 JP3861810B2 JP2002550395A JP2002550395A JP3861810B2 JP 3861810 B2 JP3861810 B2 JP 3861810B2 JP 2002550395 A JP2002550395 A JP 2002550395A JP 2002550395 A JP2002550395 A JP 2002550395A JP 3861810 B2 JP3861810 B2 JP 3861810B2
- Authority
- JP
- Japan
- Prior art keywords
- torque
- machine
- machine constant
- signal
- disturbance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P23/00—Arrangements or methods for the control of AC motors characterised by a control method other than vector control
- H02P23/14—Estimation or adaptation of motor parameters, e.g. rotor time constant, flux, speed, current or voltage
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P5/00—Arrangements specially adapted for regulating or controlling the speed or torque of two or more electric motors
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P6/00—Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
- H02P6/10—Arrangements for controlling torque ripple, e.g. providing reduced torque ripple
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Electric Motors In General (AREA)
- Feedback Control In General (AREA)
Description
この発明は、モータ等の駆動装置を用いた工作機械やロボット等の駆動機械における慣性モーメント等の機械定数を推定する機械定数推定装置に関するものである。
背景技術
モータ等の駆動装置を用いた工作機械やロボット等の駆動機械の制御を高精度かつ安定に行うためには、駆動機械の慣性モーメント等の機械定数を正確に推定することが必要となる。
第6図は従来の機械定数推定装置を含む速度制御系のブロック図で、例えば特開平11−313495号公報に制御ゲインの収束を設定しうるオートチューニング機能を持つ電動機サーボ系の制御装置として記載された構成と同様のものである。
図において、30は速度指令(図示せず)に基づきトルク指令τrを発生するトルク指令発生部、31はトルク指令τrに応じた駆動トルクτmを発生するモータ等の駆動装置、32は駆動装置31によって駆動される駆動機械、33は駆動機械32の機械速度vmを検出する速度検出器、34は機械速度vmとトルク指令τrとを入力として慣性モーメント推定値Jeを出力する機械定数推定装置、τdは駆動機械32に加わる外乱トルクである。
また、トルク指令発生部30は、機械定数推定部34の出力した慣性モーメント推定値Jeに基づき、内部のゲイン等の定数を変更し(オートチューニング)、慣性モーメント推定値Jeに応じたトルク指令τrを発生する。
第7図は従来の機械定数推定装置の構成を示す図で、電気学会論文誌Vol.114−D,No.4、p.424〜p.431に記載の逐次型最小自乗法を用いた慣性モーメント推定演算式に基づいた機械定数推定装置の構成と同様のものである。
図において、41は機械速度vmを入力し、信号処理を行って加速度変化信号daを生成する加速度変化信号生成部、42はトルク指令τrを入力し、信号処理を行ってトルク変化信号dτを生成するトルク変化信号生成部、43は加速度変化信号daと慣性モーメント推定値ホールド部47にホールドされている前回の慣性モーメント推定値Je〈k−1〉との積を演算する掛算回路、44は掛算回路43で演算した加速度変化信号daと前回の慣性モーメント推定値Je〈k−1〉との積をトルク変化信号dτから減算する演算部、45はトルク変化信号dτから加速度変化信号daと前回の慣性モーメント推定値Je〈k−1〉との積を減算した信号を入力し、前回の慣性モーメント推定値Je〈k−1〉との誤差を出力する推定ゲイン部である。また、46は推定ゲイン部45から出力された誤差に前回の慣性モーメント推定値、Je〈k−1〉を加算して慣性モーメント推定値Je〈k〉を出力する演算部、47は慣性モーメント推定値Jeをホールドする慣性モーメント推定値ホールド部である。
次に、従来の機械定数推定装置の動作について、第6図および第7図により説明する。ここで、駆動装置31で実際に発生する駆動トルクτmがトルク指令τrに一致しているものとする。
速度検出器33で検出された駆動機械32の機械速度をvmとすると、駆動機械32の機械加速度amは下記の式(1)で表される。
上式において、sはラプラス演算子である。
また、駆動機械32に加わる外乱トルクをτd、駆動機械32の慣性モーメント(真値)をJとすると、トルク指令τrは下記の式(2)で表される。
説明の簡単のため、トルク指令τrを連続時間系として表記し、トルク変化信号生成部42は、入力した現時点のトルク指令τrと前回時点のトルク指令τrとの差分を求め、この差分を微分としてトルク変化信号dτを生成するものとすると、トルク変化信号dτは下記の式(3)で表される。
また、速度検出器を用いて検出された機械速度vmには、一般に高周波数のノイズ成分が含まれているので、純粋な微分を用いるとノイズ成分を大きくしてしまい、機械定数の推定誤差の原因となる。そのために、加速度変化信号生成部41では、入力された機械速度vmに対して、純粋微分の代わりにローパスフィルタ特性を付加した疑似微分演算により、加速度変化信号daを求める。ローパスフィルタ特性をF(s)とすると、疑似加速度信号afは下記の式(4)、式(5)で表され、加速度変化信号daは下記の式(6)により演算される。
ここで、式(2)、式(3)、式(5)、式(6)をまとめ、さらにローパスフィルタ特性F(s)を理想的なフィルタとして無視すると、トルク変化信号dτは下記の式(7)で表すことができる。
式(7)において、外乱トルクτdがない定常状態では右辺第1項は0となるので、慣性モーメントJは下記の式(8)で表すことができ、駆動機械32の慣性モーメントJをトルク変化信号dτと加速度変化信号daとの比として推定できる。
第7図に示した機械定数推定装置は、推定精度の向上を図るために、式(8)の使用において逐次最小自乗法を適用したものである。
k時点における慣性モーメント推定値をJe〈k〉、前回時点(〈k−1〉時点)の慣性モーメント推定値をJe〈k−1〉として、逐次最小自乗法に基づいた慣性モーメントの推定演算式を、式(9)〜(11)に示す。
上式において、λは忘却係数と呼ばれる定数で、駆動機械の慣性モーメントの変化に対応するために1より少し小さい値に選ばれる。また、P〈k〉は、機械定数推定演算内部に持つ一種の推定ゲインパラメータであり、加速度変化信号daの大きさに応じて変化しながら更新される。
推定ゲイン部45における推定ゲインG〈k〉は、推定ゲインパラメータP〈k〉の演算式(10)を含んだ式(11)で表される。
前回時点(〈k−1〉時点)の慣性モーメント推定値Je〈k−1〉が真値であると、式(9)の右辺第二項は零となるので、慣性モーメント推定値Jeは更新されない。
しかし、前回時点(〈k−1〉時点)の慣性モーメント推定値Je〈k−1〉が真値と異なる場合には、推定誤差に応じてトルク変化の誤差項(dτ〈k〉−da〈k〉・Je〈k−1〉)が発生するので、式(9)に示すように前回時点の慣性モーメント推定値Je〈k−1〉に、トルク変化の誤差項に推定ゲインG〈k〉を掛けた値を加えた値を、k時点における慣性モーメント推定値Je〈k〉とする。
上述は外乱が存在しない(外乱トルクτd=0)場合の例について説明したが、定常状態においてインパクト的な外乱が加わった瞬間について考える。
定常状態においてモータの回転を妨げる方向にインパクト的な外乱が加わった時に発生するトルク変化をdτ’、加速度変化をda’として、式(9)
に用いると、下記の式(12)となる。
モータの回転を妨げる方向にインパクト外乱が加わった時、モータは一瞬減速させられるため、da’<0の大きな加速度変化を発生させる。一方、τd>0なる外乱トルクが加わった瞬間に、外乱トルク変化s・τd>0が発生し、そのトルク分を補うために、dτ’>0なる大きなトルク変化分が発生する。
式(12)において、トルク変化の誤差項(dτ’〈k〉−da’〈k〉・Je〈k−1〉)は、dτ’>0かつda’<0であることから急激に増加し、またda’<0であることから、右辺第二項は急激に減少するので、慣性モーメント推定値を大幅に小さな値に誤推定してしまうことになる。
慣性モーメント推定値Jeは、第6図に示すように、トルク指令発生部30における制御ゲインのオートチューニングに利用されているため、慣性モーメントを誤推定してしまうと、その結果が制御系のゲインに反映され、所望の動作を行うためのトルク指令を発生することが不可能になるので、制御性能の悪化や駆動機械の発振などの不安定現象に結びつくことがある。
第8図は従来の機械定数推定装置における外乱トルクと慣性モーメント推定値との関係を示す図で、(a)は駆動機械に与えられた外乱トルクと時間との特性を示す図、(b)は慣性モーメントと時間との特性を示す図である。図(a)において、τdは外乱トルク、また図(b)において、実線は慣性モーメントの真値、破線は慣性モーメント推定値である。
従来の機械定数推定装置においては、外乱トルクが加わる以前は、慣性モーメント推定値はほぼ真値を推定しているが、外乱トルクが加わると同時に、慣性モーメント推定値が急減し、慣性モーメント真値に対する誤差が増加している。また、推定値の急減後、もとの慣性モーメント推定値まで回復するまでにも時間がかかっている。
機械定数推定を従来の切削機械やロボットなどの負荷が変動する用途のみでなく、搬送装置において被搬送物を停止したい位置で設備に故意にぶつけて停止させる当て止めなどの用途においても利用したいという要求が生じてきたが、従来の機械定数推定装置では、上述のように、大きなインパクト外乱が駆動機械に与えられた場合には、慣性モーメント推定を正しく動作させることが困難であるため利用することができないという問題点があった。
上述は、第6図に示す一般的なオートチューニング構成における現象であるが、機械定数推定装置として第7図に示すような逐次最小自乗法などの再帰的な演算方式を用いた場合において特に顕著となる。
逐次最小自乗法による機械定数推定演算では忘却係数λを1より十分小さな値に選び、慣性モーメント推定の応答を下げれば、慣性モーメントの誤推定を防ぐことは可能であるが、忘却係数を小さくすると、推定ゲインが小さくなり、推定装置の入力に対する推定感度が鈍るため、インパクト外乱のような急激な変化に対して応答しにくくなり、また慣性モーメントの推定値が真値と異なっている場合の応答についても同様に遅くなってしまい慣性モーメント推定の性能を劣化させるという問題点があった。
第9図は、従来の外乱オブザーバ、機械定数推定装置を含む速度制御系のブロック図である。図において、31〜33、τr、τd、vm、τrは、第6図と同様であり、その説明を省略する。また、35はトルク指令発生部、36は駆動機械32の機械速度vmとインパクト外乱による外乱トルク成分を含まないトルク指令τr’とを入力して、慣性モーメント推定値Jeを出力する機械定数推定装置、37は外乱除去する外乱オブザーバである。
第9図は、大きなインパクト外乱が駆動機械に与えられた場合に、慣性モーメント推定を正しく動作させることが困難であるという従来の機械定数推定装置の問題点を解決するようにしたもので、インパクト外乱の影響を受けずに機械定数推定を行うために、外乱オブザーバにより外乱を除去したトルク指令τr’を機械定数推定装置に入力するようにしたものである。
機械定数推定装置36に入力されるトルク指令として、インパクト外乱による外乱トルク成分を含まないトルク指令τr’を使用することにより、インパクト外乱の影響を受けずに機械定数推定を行うことが可能であると考えられるが、機械定数推定装置による制御系ゲインのオートチューニングと外乱オブザーバによる外乱除去とが干渉し、動作が不安定になる恐れがあるため実用的でなく、また制御系が複雑となるため、収束特性の把握は困難であるという問題点があった。
この発明は、上述のような課題を解決するためになされたものであり、大きなインパクト外乱が駆動機械に与えられた場合においても、慣性モーメント推定を正しく動作させることができる機械定数推定装置を得ることを目的とする。
発明の開示
モータ等の駆動装置を用いた駆動機械の機械速度を入力し、低周波の加速度信号を出力する加速度信号生成部と、トルク指令を入力し雑音除去等を行い、低周波の擬似トルク信号を出力する擬似トルク信号生成部と、前記加速度信号と前回の機械定数推定値との積を演算する掛算回路と、前記低周波の擬似トルク信号からこの掛算回路の演算した積を減算し、トルク誤差信号を出力するトルク誤差信号演算部と、このトルク誤差信号を入力し、トルク誤差変化信号を出力するトルク誤差変化信号生成部と、このトルク誤差変化信号を入力し機械定数を推定する機械定数推定ゲイン部と、この機械定数推定ゲイン部から出力された誤差に前回の機械定数推定値を加算して機械定数推定値を出力する機械定数推定値演算部と、を備え、
前記トルク誤差信号演算部が演算した前記トルク誤差信号を外乱トルク推定値として出力するとともに、
前記機械定数推定ゲイン部は、前記外乱トルク推定値に基づき、機械定数推定演算の更新開始あるいは更新停止あるいは機械定数推定ゲインの変更を行うようにしたので、
駆動機械の機械定数と駆動機械に作用する外乱要素を同時に推定することができ、外乱要素による機械定数の誤推定を防止するとともに機械定数変動による外乱要素の誤推定を防止することができる。
また、フィルタ特性を有し、前記トルク誤差信号演算部が演算した前記トルク誤差信号を入力し、外乱要素の周波数成分を分解して抽出する少なくとも2個以上の外乱トルク抽出部を、備え、
前記機械定数推定ゲイン部は、前記外乱トルク抽出部が出力する周波数成分を分解した外乱トルク成分に基づき、機械定数推定演算の更新開始あるいは更新停止あるいは機械定数推定ゲインの変更を行うようにしたので、
機械または用途に特有の周波数成分の外乱要素を抽出でき、外乱要素による機械定数の誤推定を防止するとともに機械定数変動による外乱要素の誤推定を防止することができる。
さらに、前記機械定数推定ゲイン部は、前記駆動機械に外乱が作用してから定常状態に戻るまでの期間は、機械定数推定演算の更新停止あるいは機械定数推定ゲインの低減を行うようにしたので、
外乱要素による機械定数の誤推定の防止を容易に実現できるとともに機械定数変動による外乱要素の誤推定の防止を容易に実現できる。
さらにまた、前記機械定数推定ゲイン部は、前記駆動機械の機械定数の推定演算において最小自乗法等の統計的処理手法を適用し、前記外乱トルク推定値に基づき、機械定数推定に用いるパラメータ演算の更新開始あるいは更新停止あるいは定数の変更を行うようにしたので、
高応答を維持しながら、外乱要素による機械定数の誤推定の防止を容易に実現できるとともに機械定数変動による外乱要素の誤推定の防止を容易に実現できる。
また、モータ等の駆動装置を用いた駆動機械の機械速度を入力し、加速度変化信号を出力する加速度変化信号生成部と、
トルク指令を入力し雑音除去等を行い、トルク変化信号を出力するトルク変化信号生成部と、
前期加速度変化信号を入力し、加速度変化ベクトルを生成する加速度変化ベクトル生成部と、
機械定数推定値と外乱トルク成分の推定値とからなる前回のパラメータベクトル推定値と前記加速度変化ベクトルとの行列積を演算する掛算回路と、
前記トルク変化信号からこの掛算回路の演算した積を減算し、トルク誤差変化信号を出力する減算器と、
このトルク誤差変化信号を入力し、誤差ベクトルを推定する機械定数推定ゲイン部と、
この機械定数推定ゲイン部から出力された誤差ベクトルに前回のパラメータベクトル推定値を加算して機械定数推定値を出力する加算器と、を備え、機械定数推定値と外乱トルク成分の推定値を一組にしてパラメータベクトル推定値として同時に推定するようにしたので、
外乱トルクによる慣性モーメント推定の誤動作を防ぎ、高精度な推定結果を得ることができる。
発明を実施するための最良の形態
実施の形態1.
第1図はこの発明の実施の形態1に係る機械定数推定装置を含む速度制御系のブロック図である。図において、31,32、33、τr、τd、vmは、従来例である第6図と同様であり、その説明を省略する。また、1はトルク指令τrを発生するトルク指令発生部、2は機械速度vmとトルク信号τrを入力として慣性モーメント推定値と外乱推定値とを出力する機械定数推定装置である。トルク指令発生部1は、機械定数推定装置2の出力した慣性モーメント推定値に応じてゲイン等の調整を行う。
従来例である第6図、第7図および第9図に示した機械定数推定装置は、慣性モーメント推定値を出力するものであったが、実施の形態1における機械定数推定装置は慣性モーメント推定値と外乱推定値とを出力するようにしたものである。
第2図はこの発明の実施の形態1に係る機械定数推定装置2の構成を示す図である。図において、11は機械速度vmを入力し、低周波の加速度信号afを出力する加速度信号生成部、12はトルク指令τrを入力し、雑音除去等を行い低周波の擬似トルク信号τfを出力する擬似トルク信号生成部である。また、13は低周波の加速度信号afと慣性モーメント推定値ホールド部18にホールドされている前回の慣性モーメント推定値Je〈k−1〉との積を演算する掛算回路、14は掛算回路13で演算した低周波の加速度信号afと前回の慣性モーメント推定値Je〈k−1〉との積を低周波の擬似トルク信号τfから減算し、トルク誤差信号τerrを演算するトルク誤差信号演算部である。また、15はトルク誤差信号τerrを入力し、トルク誤差変化信号dτerrを出力するトルク誤差変化信号生成部、16はトルク誤差変化信号dτerrを入力し機械定数を推定する機械定数推定ゲイン部である。また、17は機械定数推定ゲイン部16から出力された機械定数推定値に前回の機械定数推定値としての慣性モーメント推定値Je〈k−1〉を加算して機械定数推定値としての慣性モーメント推定値Je〈k〉を出力する慣性モーメント推定値演算部、18は慣性モーメント推定値Je〈k〉をホールドする慣性モーメント推定値ホールド部である。
次に実施の形態1における機械定数推定装置の動作を第1図および第2図により説明する。
ここで、駆動装置31により実際に発生する駆動トルクτmはトルク指令τrに一致するものとする。駆動機械32に加わる外乱トルクをτd、駆動機械32の機械速度をvm、駆動機械32の慣性モーメント(真値)をJとすると、駆動機械32の機械加速度am、トルク指令τrは下記の式(13)および式(14)で表される。
ただし、sはラプラス演算子を表す。
加速度信号生成部11は、機械速度vmを入力し、下記の式(15)で表す演算を行い、高周波数成分を除去した加速度信号afを出力する。
ここで、F(s)は、機械速度vmに含まれる高周波成分を除去するためのローパスフィルタ特性である。
擬似トルク信号生成部12は、トルク指令τrを入力し、加速度信号生成部11と同一のローパスフィルタ特性を付加して、下式(16)で表す演算を行い、高周波数成分を除去した擬似トルク信号τfを出力する。
加速度信号afと前回の慣性モーメント推定値Je〈k−1〉との積と、擬似トルク信号τfとの差をトルク誤差信号τerrとすると、トルク誤差信号τerrは式(17)のように表される。
ここで、前回の慣性モーメント推定値Je〈k−1)が真値を示していると仮定し、またローパスフィルタ特性が理想的なものであるとして無視すると、式(15)および式(16)を用いて、トルク誤差信号は、式(18)のように表される。
式(18)を式(14)と比較すれば、トルク誤差信号τerrが外乱トルクτdに他ならないことがわかる。つまり、上記の構成により外乱トルクτd0を推定できる。
第7図に示した従来の機械定数推定装置においては、慣性モーメントJを推定するために必要な加速度変化信号da、トルク変化信号dτを、加速度変化信号生成部41およびトルク変化信号生成部42で直接生成していたので、外乱推定値を求めることができなかった。
実施の形態1では、加速度信号生成部11で、高周波数成分を除去した加速度信号afを出力し、擬似トルク信号生成部12で、高周波数成分を除去した擬似トルク信号τfを出力し、トルク誤差信号演算部14でトルク誤差信号τerrを演算した後、トルク誤差変化信号生成部15でトルク誤差変化信号dτerrを生成するようにしたので、一連の処理において、外乱トルク推定値と慣性モーメント推定値とを演算することができる。
また、上式(17)の演算において、前回の慣性モーメント推定値Je〈k−1〉が真値からずれている場合であっても、インパルス外乱などの外乱トルク成分は、慣性モーメントのずれにより生じるトルク誤差成分と周波数成分が異なる場合が多いので、外乱トルク成分を取り出すことが可能である。
トルク誤差変化信号生成部15は、トルク誤差信号τerrを入力とし、これに微分演算を行ってトルク誤差変化信号を生成する。トルク誤差変化信号をdτerrとすると、トルク誤差変化信号生成部15の演算は式(17)を用いて式(19)のように表すことができる。
上式(19)に、式(15)および式(16)を代入し、さらに式(3)および式(6)を利用して変形すると、ローパスフィルタ特性F(s)が理想的であるとして無視した場合、式(19)は式(20)のように置き換えられる。
仮に、機械定数推定ゲイン部16が、従来の機械定数推定装置のものと同様であると仮定するならば、上式(20)に基づいて、慣性モーメント推定演算式は、式(21)と表すことができる。
上式(21)は、従来の機械定数推定装置における慣性モーメントの推定演算式(9)と同様のものであり、従来の機械定数推定装置と同等の機械定数推定演算を行うことも可能である。
機械定数推定ゲイン部16は、推定ゲインG〈k〉、推定ゲインパラメータP〈k〉を下式(22)および下式(23)式に基づいて演算し、トルク誤差変化信号dτerrに推定ゲインを乗算する。
ここで、nf〈k〉は外乱トルク推定値τd0に基づくフラグである。フラグnf〈k〉は、外乱トルクの加わる前の定常状態ではnf〈k〉=1とし、外乱トルクが或るしきい値以上の大きさになった瞬間から定常状態に戻るまでの期間はnf〈k〉=0にする。従って、nf〈k〉=0の期間は推定ゲインG〈k〉=0になるので、式(21)に基づく慣性モーメント推定値の更新は行われない。つまり、機械定数推定ゲイン部16の出力は、外乱トルク成分を含まない慣性モーメント成分が得られることになる。
第3図はこの発明の実施の形態1に係る外乱トルクと慣性モーメント推定値との関係を示す図で、(a)は駆動機械に与えられた外乱トルクと時間との特性を示す図、(b)は慣性モーメントと時間との特性を示す図である。図(a)において、τdは外乱トルク、また図(b)において、実線は慣性モーメントの真値、破線は実施の形態1における慣性モーメント推定値である。
実施の形態1による機械定数推定装置においては、慣性モーメントの変化とそれ以外の外乱トルクを精度良く分離できるようにしたので、点線に示すように外乱トルクが引加されても慣性モーメント推定値はほとんど影響を受けず真値に対してほぼ一致しており、機械定数の誤推定が防止されている。
実施の形態1における機械定数推定装置では、外乱トルクを検出し、外乱トルクが加わるまでの定常状態では従来通りの機械定数推定を行い、外乱トルクによる悪影響の生じうる期間は慣性モーメント推定値の更新を停止するようにしたので、インパクト外乱の影響を受けない高精度な機械定数推定値を得ることができる。また、従来使用できなかった当て止めなどの用途にも利用することができる。
また、実施の形態1における機械定数推定装置では、慣性モーメントが変化するような場合でもその時々の慣性モーメント推定値に基づいて外乱トルクτd0の演算がなされるため、高精度な外乱推定値を得ることができる。
実施の形態2.
第4図はこの発明の実施の形態2に係る機械定数推定装置の構成を示す図である。図において、11〜15、17、18は、第2図と同様であり、その説明を省略する。また、21はローパスフィルタ特性を有し、トルク誤差信号τerrを入力し、低周波の外乱トルク成分τd1を出力する低周波外乱トルク抽出部、22はハイパスフィルタ特性を有し、トルク誤差信号τerrを入力し、高周波の外乱トルク成分τd2を出力する高周波外乱トルク抽出部、23は低周波トルク成分τd1と高周波トルク成分τd2を反映させて慣性モーメント推定演算する機械定数推定ゲイン部である。
なお、低周波外乱トルク抽出部21におけるローパスフィルタの特性は、加速度信号生成部11、擬似トルク信号生成部12で利用しているローパスフィルタの特性とは異なるものであり、加速度信号生成部11、擬似トルク信号生成部12においてローパスフィルタで高調波ノイズ成分を除去した信号の中から、さらに低周波側の成分だけを抽出するものである。また、高周波外乱トルク抽出部22におけるハイパスフィルタは、加速度信号生成部11、擬似トルク信号生成部12においてローパスフィルタで高調波ノイズ成分を除去した信号の中から、高周波側の成分だけを抽出するものである。
実施の形態2における機械定数推定装置の動作を説明する。
駆動機械に作用する外乱τdは、インパクト外乱の場合には外乱トルクの高周波成分が多くなり、摩擦等の場合には低周波成分の外乱が多く含まれるといったように、外乱トルクの種類によって必ずしも特徴が一定でない。このため、実施の形態2では、低周波外乱トルク抽出部21および高周波外乱トルク抽出部22の2つのフィルタ要素を用いて外乱成分を周波数分割して抽出するようにしたものである。
高周波の外乱に対しては、実施の形態1と同様に、外乱トルクが或るしきい値以上の大きさになった瞬間から定常状態に戻るまでの一定期間、推定ゲインを0にし慣性モーメント推定値の更新は行わないこととし、一方、定常的に機械定数推定に作用する低周波の外乱は、低周波の外乱推定値の大きさに応じて、機械定数推定ゲインを小さくする、または機械定数推定ゲイン部23への入力から低周波の外乱推定値を除去する等により低周波の外乱要素による機械定数推定の誤動作を防止する。
実施の形態1では、トルク変化をインパクト外乱として一括して扱っているので、例えば摩擦の大きな駆動機械などの用途においては、低周波の外乱成分により機械定数推定値に誤差を生じさせるが、
実施の形態2では、低周波の外乱成分を分離して抽出するため、低周波の外乱成分による誤差を抑えることが可能となるため、低周波成分の外乱が多く含まれる摩擦の非常に大きな駆動機械などの用途にも使用できる。
上述では、2つのフィルタ要素を用いて外乱成分を周波数分割して抽出する例を示したが、フィルタ要素を追加することにより、さらに細かく外乱の周波数成分を分けて抽出することができ、より高精度に外乱要素の影響を考慮することができる。
実施の形態3.
第5図はこの発明の実施の形態3に係る機械定数推定装置の構成を示す図である。図において、41、42、vm、τr、da、dτは、従来例である第7図と同様であり、その説明を省略する。24は加速度変化信号daを入力し加速度変化ベクトルφを出力する加速度変化ベクトル生成部、25は加速度変化ベクトルφとパラメータベクトル推定値ホールド部29にホールドされている前回のパラメータベクトル推定値θe〈k−1〉との積を演算する掛算回路、26は掛算回路25で演算した加速度変化ベクトルφと前回の慣性モーメント推定値Je〈k−1〉との積をトルク変化信号dτから減算する演算部、27は減算器26から出力されたトルク誤差変化信号を入力しパラメータ誤差ベクトル推定値を出力する機械定数推定ゲイン部、28は機械定数推定ゲイン部27から出力された誤差に前回のパラメータベクトル推定値θe〈k−1〉を加算してパラメータベクトル推定値θe〈k〉を出力する演算部、29はパラメータベクトル推定値θeをホールドするパラメータベクトル推定値ホールド部である。また、θeは慣性モーメント推定値と外乱トルク変化推定値とを一組としたパラメータベクトル推定値である。
また、図において太線で示した加速度変化ベクトル生成部24の出力、機械定数推定ゲイン部27の出力、演算部28の出力、パラメータベクトル推定値ホールド部29の出力はベクトルを表している。
次に、実施の形態3における機械定数推定装置の動作を説明する。
トルク変化信号dτに関する式(7)において、外乱トルクの微分であるs・τdを、改めて外乱トルク変化dτdとおくと、式(7)は式(24)に置き換えられる。
実施の形態1および実施の形態2においては、慣性モーメントに対して最小自乗自乗法を利用して推定を行った例を示したが、実施の形態3では慣性モーメントJと外乱トルク変化dτdの両方に対して同時に最小自乗法を利用して、パラメータベクトル推定値θeとして推定するようにしたものである。
まず式(25)〜式(27)に示すような、定数yとベクトルφ、ベクトルθを導入する。
ここで、式(27)の[]Tはベクトルの転置を表す。上式を用いて式(24)を書き換えると、式(28)が得られる。
式(28)は、トルク変化信号dτと加速度変化信号daより得られるy、φを用いてパラメータベクトルθを推定できることを示している。
第5図に示した機械定数推定装置は、式(28)に対して逐次最小自乗法を適用したものである。
k時点における慣性モーメント推定値をJe〈k〉、外乱トルク変化の推定値をdτde〈k〉とすると、パラメータベクトル推定値θeを式(29)に定義できる。
逐次最小自乗法に基づいた推定パラメータベクトルθeの推定演算式は、式(30)〜(32)で表すことができる。
上式において、P、φ、θeがベクトル量であることを除いては、式内部の演算は従来例における演算と同一である。
従来例における機械定数推定装置おいては、慣性モーメントのみしか推定しないため、外乱トルクが混入した場合、その成分は慣性モーメント誤差成分として発生することになる。
一方、実施の形態3においては、予め慣性モーメントと外乱トルクとの関係式を与え、各々が別々に存在するものとして推定を行う。例えば、インパクト外乱成分が入力信号に含まれた場合、その成分は主に外乱トルク変化成分として推定出力されるため、慣性モーメント推定値にはその成分がほとんど出現してこない。結果として外乱トルクによる慣性モーメント推定の誤動作を防ぎ、高精度な推定結果を得ることができる。
従って実施の形態3の機械定数推定装置においては、外乱トルクの大きさにより慣性モーメント推定動作の停止等を行うという方法を用いなくても、実施の形態1の機械定数推定装置と同じ効果を得ることが可能である。
ところで、上述の実施の形態1、実施の形態2および実施の形態3では、機械定数として慣性モーメントを推定する例を説明したが、その他の機械定数として、弾性係数や粘性定数等の定数をも推定するようにしてもよい。
また、上述では、外乱トルクが或るしきい値以上の大きさになった場合に、機械定数推定の演算を停止する例を説明したが、例えば外乱トルクの大きさに応じて、機械定数推定ゲインを小さくするというようにゲイン変更を行っても良い。
また、上述では、外乱トルクが或るしきい値以上の大きさになった場合に、機械定数推定の演算を停止する例を説明したが、推定ゲインパラメータP〈k〉の演算についても、停止もしくはゲイン変更などを行っても良い。
また、上述では、外乱トルクが或るしきい値以上の大きさになった瞬間から定常状態に戻るまでの一定期間、フラグnf〈k〉=0にする例を説明したが、この定常状態に戻る時点の検出は、機械速度vmやトルク指令τr等の制御応答波形を観察し、定常状態に戻る時点を実際に検出しても良いし、外乱トルクが或るしきい値以下に戻ってから一定の遅延時間を持たせることによって得ても良い。
また、上述では、速度検出器33を用いて機械速度vmを検出する例を説明したが、位置検出器を用いて機械位置を検出し、機械位置を差分あるいは微分した信号を機械速度vmとするような、等価な構成にしても良い。
また、上述では、逐次最小自乗法により駆動機械の慣性モーメント推定値Jeを計算した例を説明したが、電気学会論文誌Vol.114−D、No.4、P.424−p.431に記載の固定トレース法や、「MATLABによる制御のためのシステム同定」(東京電気大学出版局)、p.6、に記載のように非正規化勾配法など他の統計処理手法を用いても良い。
また、上述では、トルク誤差変化信号生成部15は、加速度信号afと前回の慣性モーメント推定値との積信号と擬似トルク信号τfとの差信号τerrを入力とし、これに微分演算を行ってトルク誤差変化信号dτerrを生成した例を説明したが、二次のローパスフィルタ特性等を付加することにより高周波成分を除去するようにしてもよい。
また、上述では、機械定数と外乱を同一のアルゴリズムにおいて推定する例を説明したが、外乱オブザーバを別に用意してその結果を機械定数推定ゲインに反映させる形でも、同一の効果が得られる。
産業上の利用可能性
以上のように、本発明に係る機械定数推定装置は、大きなインパクト外乱が駆動機械に与えられた場合でも、慣性モーメント推定を正しく動作させることができるので、搬送装置において被搬送物を停止したい位置で設備に故意にぶつけて停止させる当て止めなどの用途において用いられるのに適している。
【図面の簡単な説明】
第1図はこの発明の実施の形態1に係る機械定数推定装置を含む速度制御系のブロック図である。
第2図はこの発明の実施の形態1に係る機械定数推定装置2の構成を示す図である。
第3図はこの発明の実施の形態1に係る外乱トルクと慣性モーメント推定値との関係を示す図である。
第4図はこの発明の実施の形態2に係る機械定数推定装置の構成を示す図である。
第5図はこの発明の実施の形態3に係る機械定数推定装置の構成を示す図である。
第6図は従来の機械定数推定装置を含む速度制御系のブロック図である。
第7図は従来の機械定数推定装置34の構成を示す図である。
第8図は従来の機械定数推定装置における外乱トルクと慣性モーメント推定値との関係を示す図である。
第9図は従来の外乱オブザーバ、機械定数推定装置を含む速度制御系のブロック図である。
Claims (5)
- モータ等の駆動装置を用いた駆動機械の機械速度を入力し、低周波の加速度信号を出力する加速度信号生成部と、
トルク指令を入力し雑音除去等を行い、低周波の擬似トルク信号を出力する擬似トルク信号生成部と、
前記低周波の加速度信号と前回の機械定数推定値との積を演算する掛算回路と、
前記低周波の擬似トルク信号からこの掛算回路の演算した積を減算し、トルク誤差信号を出力するトルク誤差信号演算部と、
このトルク誤差信号を入力し、トルク誤差変化信号を出力するトルク誤差変化信号生成部と、
このトルク誤差変化信号を入力し機械定数を推定する機械定数推定ゲイン部と、
この機械定数推定ゲイン部から出力された誤差に前回の機械定数推定値を加算して機械定数推定値を出力する機械定数推定値演算部と、を備え、
前記トルク誤差信号演算部が演算した前記トルク誤差信号を外乱トルク推定値として出力するとともに、
前記機械定数推定ゲイン部は、前記外乱トルク推定値に基づき、機械定数推定演算の更新開始あるいは更新停止あるいは機械定数推定ゲインの変更を行うようにしたことを特徴とする機械定数推定装置。 - フィルタ特性を有し、前記トルク誤差信号演算部が演算した前記トルク誤差信号を入力し、外乱トルク成分の周波数成分を分解して抽出する少なくとも2個以上の外乱トルク抽出部を、備え、
前記機械定数推定ゲイン部は、この外乱トルク抽出部が出力する周波数成分を分解した外乱トルク成分に基づき、機械定数推定演算の更新開始あるいは更新停止あるいは機械定数推定ゲインの変更を行うようにしたことを特徴とする請求の範囲第1項に記載の機械定数推定装置。 - 前記機械定数推定ゲイン部は、前記駆動機械に外乱が作用してから定常状態に戻るまでの期間は、機械定数推定演算の更新停止あるいは機械定数推定ゲインの低減を行うようにしたことを特徴とする請求の範囲第1項ないし第2項に記載の機械定数推定装置。
- 前記機械定数推定ゲイン部は、前記駆動機械の機械定数の推定演算において最小自乗法等の統計的処理手法を適用し、前記外乱トルク推定値に基づき、機械定数推定に用いるパラメータ演算の更新開始あるいは更新停止あるいは定数の変更を行うようにしたことを特徴とする請求の範囲第1項ないし第2項に記載の機械定数推定装置。
- モータ等の駆動装置を用いた駆動機械の機械速度を入力し、加速度変化信号を出力する加速度変化信号生成部と、
トルク指令を入力し雑音除去等を行い、トルク変化信号を出力するトルク変化信号生成部と、
前期加速度変化信号を入力し、加速度変化ベクトルを生成する加速度変化ベクトル生成部と、
機械定数推定値と外乱トルク成分の推定値とからなる前回のパラメータベクトル推定値と前記加速度変化ベクトルとの行列積を演算する掛算回路と、
前記トルク変化信号からこの掛算回路の演算した積を減算し、トルク誤差変化信号を出力する減算器と、
このトルク誤差変化信号を入力し、誤差ベクトルを推定する機械定数推定ゲイン部と、
この機械定数推定ゲイン部から出力された誤差ベクトルに前回のパラメータベクトル推定値を加算して機械定数推定値を出力する加算器と、を備え、
機械定数推定値と外乱トルク成分の推定値を一組にしてパラメータベクトル推定値として同時に推定することを特徴とする機械定数推定装置。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2000/008865 WO2002049201A1 (fr) | 2000-12-14 | 2000-12-14 | Dispositif d'evaluation des constantes d'une machine |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2002049201A1 JPWO2002049201A1 (ja) | 2004-04-15 |
JP3861810B2 true JP3861810B2 (ja) | 2006-12-27 |
Family
ID=11736799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002550395A Expired - Fee Related JP3861810B2 (ja) | 2000-12-14 | 2000-12-14 | 機械定数推定装置 |
Country Status (6)
Country | Link |
---|---|
US (1) | US6736018B2 (ja) |
JP (1) | JP3861810B2 (ja) |
DE (1) | DE10085405T1 (ja) |
GB (1) | GB2376303B (ja) |
TW (1) | TW494055B (ja) |
WO (1) | WO2002049201A1 (ja) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004030500A (ja) * | 2002-06-28 | 2004-01-29 | Fanuc Ltd | モータ制御装置 |
US7587442B2 (en) * | 2004-05-14 | 2009-09-08 | Gm Global Technology Operations, Inc. | Method of determining the derivative of an input signal |
JP4591176B2 (ja) * | 2005-04-20 | 2010-12-01 | 株式会社明電舎 | エンジンの慣性モーメント測定装置 |
JP3796261B1 (ja) * | 2005-12-02 | 2006-07-12 | 山洋電気株式会社 | モータの負荷イナーシャ推定方法 |
JP4565034B2 (ja) * | 2008-12-16 | 2010-10-20 | ファナック株式会社 | イナーシャ推定を行う制御装置及び制御システム |
DE102010044644B4 (de) | 2010-09-07 | 2018-12-27 | Robert Bosch Gmbh | Verfahren zur Kollisionserkennung für eine Antriebseinheit |
JP5657425B2 (ja) * | 2011-02-25 | 2015-01-21 | Ntn株式会社 | 電気自動車 |
US8315823B2 (en) | 2011-04-20 | 2012-11-20 | Bertec Corporation | Force and/or motion measurement system having inertial compensation and method thereof |
US8315822B2 (en) | 2011-04-20 | 2012-11-20 | Bertec Corporation | Force measurement system having inertial compensation |
JP6433387B2 (ja) * | 2015-07-13 | 2018-12-05 | 三菱電機株式会社 | 交流回転機の制御装置及び交流回転機の慣性モーメント演算方法 |
JP6698733B2 (ja) * | 2018-04-06 | 2020-05-27 | ファナック株式会社 | モータエンコーダ及びセンサを用いて学習制御を行うロボットシステム |
CN113084828B (zh) * | 2021-04-02 | 2022-11-11 | 上海电气集团股份有限公司 | 一种运动控制方法、装置、设备和存储介质 |
KR20230037722A (ko) * | 2021-09-09 | 2023-03-17 | 현대자동차주식회사 | 차량의 주행 상태 판단 방법 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0678579A (ja) * | 1992-06-26 | 1994-03-18 | Meidensha Corp | モータイナーシャ推定装置およびオートチューニング方法 |
JPH0695703A (ja) | 1992-09-10 | 1994-04-08 | Toyo Electric Mfg Co Ltd | 機械定数の推定方法 |
JP3545487B2 (ja) | 1995-03-09 | 2004-07-21 | ファナック株式会社 | モータで駆動する機械系のイナーシャおよび摩擦トルクの推定方法 |
JPH11313495A (ja) | 1998-04-27 | 1999-11-09 | Mitsubishi Heavy Ind Ltd | 電動機サーボ系の制御装置 |
US6827051B2 (en) * | 1999-12-03 | 2004-12-07 | Nissan Motor Co., Ltd. | Internal EGR quantity estimation, cylinder intake air quantity calculation, valve timing control, and ignition timing control |
US6226585B1 (en) * | 2000-04-18 | 2001-05-01 | Ford Global Technologies, Inc. | Torque estimation method for an internal combustion engine |
US6379283B1 (en) * | 2000-04-18 | 2002-04-30 | Ford Global Technologies, Inc. | Torque estimation method for an internal combustion engine |
JP3467455B2 (ja) * | 2000-07-17 | 2003-11-17 | 本田技研工業株式会社 | 内燃機関の気筒別空燃比推定装置 |
-
2000
- 2000-12-14 JP JP2002550395A patent/JP3861810B2/ja not_active Expired - Fee Related
- 2000-12-14 WO PCT/JP2000/008865 patent/WO2002049201A1/ja active IP Right Grant
- 2000-12-14 DE DE10085405T patent/DE10085405T1/de not_active Withdrawn
- 2000-12-14 US US10/182,794 patent/US6736018B2/en not_active Expired - Fee Related
- 2000-12-14 GB GB0218435A patent/GB2376303B/en not_active Expired - Fee Related
- 2000-12-29 TW TW089128280A patent/TW494055B/zh not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
US20030233887A1 (en) | 2003-12-25 |
GB0218435D0 (en) | 2002-09-18 |
JPWO2002049201A1 (ja) | 2004-04-15 |
DE10085405T1 (de) | 2003-08-21 |
WO2002049201A1 (fr) | 2002-06-20 |
US6736018B2 (en) | 2004-05-18 |
GB2376303B (en) | 2004-07-14 |
GB2376303A (en) | 2002-12-11 |
TW494055B (en) | 2002-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3861810B2 (ja) | 機械定数推定装置 | |
JP3302235B2 (ja) | 電動機の位置制御装置 | |
JP2007156699A (ja) | モータの負荷イナーシャ推定方法 | |
JP5322534B2 (ja) | 制御装置およびモータ制御装置 | |
JP2001352773A (ja) | オンラインオートチューニングサーボ制御装置 | |
Chan et al. | Application of Unscented Kalman filter to sensorless permanent-magnet synchronous motor drive | |
JP2021125898A (ja) | モータ制御装置 | |
JP5515885B2 (ja) | 電気車制御装置 | |
JP4305698B2 (ja) | 同期電動機の位置および速度推定装置 | |
JP5256704B2 (ja) | 慣性モーメント推定装置 | |
JP3683121B2 (ja) | 駆動機械の機械定数推定装置 | |
JP4533677B2 (ja) | 電力変換装置 | |
Jang et al. | Stator Flux Linkage Estimation of Synchronous Machines Based on Integration Error Estimation for Improved Transient Performance | |
JP2001069777A (ja) | 同期電動機の位置および速度推定装置 | |
KR100491693B1 (ko) | 기계정수 추정장치 | |
EP4102716A1 (en) | Motor control device | |
JP2634333B2 (ja) | 誘導電動機の制御装置 | |
JP7459906B2 (ja) | モータ制御装置 | |
Xu et al. | An optimized nonlinear extended state observer for sensorless IPM synchronous motor drives | |
JPH1198891A (ja) | 誘導電動機のトルク制御装置 | |
JPH04304181A (ja) | 軸ねじり振動抑制制御装置 | |
TWI720662B (zh) | 伺服放大器和伺服系統 | |
JPH08168300A (ja) | 誘導電動機のベクトル制御装置 | |
JPH07177778A (ja) | モータ制御装置 | |
Zerbo et al. | Real time flux and torque estimator for induction machines |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7421 Effective date: 20040629 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20060905 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060918 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091006 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101006 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111006 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121006 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131006 Year of fee payment: 7 |
|
LAPS | Cancellation because of no payment of annual fees |