(57)【要約】【課題】モータドライブシステムの軸トルク加振制御において、共振・非周期外乱・周期外乱の影響をすべて除去しつつ、歪み成分を含んだエンジン加振トルク指令値に追従させる。【解決手段】低域トルク制御器９により、トルク指令値τ*とトルク検出値τdetに基づいて、低域トルク制御器出力τdc*を出力し、加振トルク制御器１１によりトルク指令値τ*とトルク検出値τdetと回転位相検出値θとに基づいて加振トルク指令値τpd*を出力する。そして、高域共振抑制制御器において、低域トルク制御器出力τdc*と加振トルク指令値τpd*とを加算した補正トルク指令値τr*と、トルク検出値τdetと、に基づいてインバータトルク指令値τinv*を出力する。【選択図】図５
2018007287.tif 000002
【特許請求の範囲】
【請求項１】
モータの軸トルクを、インバータを用いて制御するモータドライブシステムであって、
トルク指令値とトルク検出値に基づいて、低域トルク制御器出力を出力する低域トルク制御器と、
前記トルク指令値と前記トルク検出値と回転位相検出値とに基づいて加振トルク指令値を出力する加振トルク制御器と、
前記低域トルク制御器出力と前記加振トルク指令値とを加算した補正トルク指令値と、前記トルク検出値と、に基づいてインバータトルク指令値を出力する高域共振抑制制御器と、を備えたことを特徴とするモータドライブシステム。
【請求項２】
前記高域共振抑制制御器は、
前記補正トルク指令値が、μ設計制御器の指令値入力からμ設計制御器出力までの伝達特性を通った出力と、
前記トルク検出値が、μ設計制御器の検出値入力からμ設計制御器出力までの伝達特性を通った出力と、を加算して、インバータトルク指令値を出力するμ設計制御器を有し、
前記低域トルク制御器はＰＩＤ制御を行うことを特徴とする請求項１記載のモータドライブシステム。
【請求項３】
前記加振トルク制御器は、
前記トルク指令値と、前記回転位相検出値に制御する加振周波数成分とトルクリプル成分の次数ｎを乗算したｎ次回転位相と、に基づいて、前記加振トルク指令値のｎ次周波数成分ベクトルを出力する加振周波数成分抽出器と、
前記トルク検出値と、前記ｎ次回転位相と、に基づいて、周期外乱のｎ次周波数成分ベクトルを出力するリプル抑制周波数成分抽出器と、
前記ｎ次回転位相に基づいて、ｎ次回転周波数を出力する速度変換器と、
前記加振トルク指令値のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗算し、加振による周期外乱指令値のｎ次周波数成分ベクトルを出力する第１逆モデル乗算部と、
前記周期外乱のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗算して操作量推定値のｎ次周波数成分ベクトルを出力する第２逆モデル乗算部と、
前記操作量推定値のｎ次周波数成分ベクトルから、加振トルク制御器出力のｎ次周波数成分ベクトルがローパスフィルタを介した値を減算して周期外乱推定値のｎ次周波数成分ベクトルを出力する第１減算器と、
前記加振による周期外乱指令値のｎ次周波数成分ベクトルから、前記周期外乱推定値のｎ次周波数成分ベクトルを減算して加振トルク制御器出力のｎ次周波数成分ベクトルを出力する第２減算器と、
前記加振トルク制御器出力のｎ次周波数成分ベクトルと前記ｎ次回転位相に基づいて前記加振トルク指令値を出力する補償信号合成部と、
を備えたことを特徴とする請求項１又は２記載のモータドライブシステム。
【請求項４】
前記加振トルク制御器は、
前記トルク指令値と前記トルク検出値とのトルク偏差を算出する第３減算部と、
前記トルク偏差と、回転位相検出値と制御する加振周波数成分とトルクリプル周波数成分の次数ｎとを乗算したｎ次回転位相と、に基づいて周期外乱のｎ次周波数成分ベクトルを出力する周波数成分抽出器と、
前記ｎ次回転位相に基づいて、ｎ次回転周波数を出力する速度変換器と、
前記周期外乱のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗算して操作量推定値のｎ次周波数成分ベクトルを出力する逆モデル乗算部と、
前記操作量推定値のｎ次周波数成分ベクトルと、加振トルク制御器出力のｎ次周波数成分ベクトルがローパスフィルタを介した値と、を加算して加振トルク制御器出力のｎ次周波数成分ベクトルを出力する加算器と、
前記加振トルク制御器出力のｎ次周波数成分ベクトルと前記ｎ次回転位相に基づいて前記加振トルク指令値を出力する補償信号合成部と、
を備えたことを特徴とする請求項１又は２記載のモータドライブシステム。
【請求項５】
前記加振トルク制御器は、
前記トルク指令値と前記トルク検出値とのトルク偏差を算出する第３減算部と、
前記トルク偏差と、回転位相検出値と制御する加振周波数成分とトルクリプル周波数成分の次数ｎとを乗算したｎ次回転位相と、に基づいて、周期外乱のｎ次周波数成分ベクトルを出力する周波数成分抽出器と、
前記ｎ次回転位相に基づいて、ｎ次回転周波数を出力する速度変換器と、
前記周期外乱のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗じて、操作量推定値のｎ次周波数成分ベクトルを求める逆モデル乗算部と、
前記操作量推定値のｎ次周波数成分ベクトルを積分し、加振トルク制御器出力のｎ次周波数成分ベクトルを出力する積分器と、
前記加振トルク制御器出力のｎ次周波数成分ベクトルと前記ｎ次回転位相に基づいて前記加振トルク指令値を出力する補償信号合成部と、
を備えたことを特徴とする請求項１又は２記載のモータドライブシステム。
【請求項６】
前記加振トルク制御器は、
前記トルク指令値と、前記回転位相検出値に制御する加振周波数成分とトルクリプル周波数成分の次数ｎを乗算したｎ次回転位相と、に基づいて、前記加振トルク指令値のｎ次周波数成分ベクトルを出力する加振周波数成分抽出器と、
前記トルク検出値と、前記ｎ次回転位相と、に基づいて、周期外乱のｎ次周波数成分ベクトルを出力するリプル抑制周波数成分抽出器と、
前記ｎ次回転位相に基づいて、ｎ次回転周波数を出力する速度変換器と、
前記加振トルク指令値のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗算した加振による周期外乱指令値のｎ次周波数成分ベクトルを出力する第１逆モデル乗算部と、
前記周期外乱のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗算して操作量推定値のｎ次周波数成分ベクトルを出力する第２逆モデル乗算部と、
前記加振による周期外乱指令値のｎ次周波数成分ベクトルにオブザーバゲインを乗算して第２減算器に出力する第１乗算器と、
前記操作量推定値のｎ次周波数成分ベクトルにオブザーバゲインを乗算して第１減算器に出力する第２乗算器と、
前記第２乗算器の出力から、加振トルク制御器出力のｎ次周波数成分ベクトルがローパスフィルタを介した値を減算して周期外乱推定値のｎ次周波数成分ベクトルを出力する第１減算器と、
前記第１乗算器の出力から、前記周期外乱推定値のｎ次周波数成分ベクトルを減算して加振トルク制御器出力のｎ次周波数成分ベクトルを出力する第２減算器と、
前記加振トルク制御器出力のｎ次周波数成分ベクトルと前記ｎ次回転位相に基づいて前記加振トルク指令値を出力する補償信号合成部と、
を備えたことを特徴とする請求項１又は２記載のモータドライブシステム。
【請求項７】
前記加振トルク制御器は、
前記トルク指令値と前記トルク検出値とのトルク偏差を算出する第３減算部と、
前記トルク偏差と、回転位相検出値と制御する加振周波数成分とトルクリプル周波数成分の次数ｎとを乗算したｎ次回転位相と、に基づいて周期外乱のｎ次周波数成分ベクトルを出力する周波数成分抽出器と、
前記ｎ次回転位相に基づいて、ｎ次回転周波数を出力する速度変換器と、
前記周期外乱のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗算して操作量推定値のｎ次周波数成分ベクトルを出力する逆モデル乗算部と、
前記操作量推定値のｎ次周波数成分ベクトルにオブザーバゲインを乗算して加算器に出力する乗算器と、
前記乗算器の出力と、加振トルク制御器出力のｎ次周波数成分ベクトルがローパスフィルタを通った値と、を加算して加振トルク制御器出力のｎ次周波数成分ベクトルを出力する加算器と、
前記加振トルク制御器出力のｎ次周波数成分ベクトルと前記ｎ次回転位相に基づいて前記加振トルク指令値を出力する補償信号合成部と、
を備えたことを特徴とする請求項１又は２記載のモータドライブシステム。
【請求項８】
前記加振トルク制御器は、
前記トルク指令値と前記トルク検出値とのトルク偏差を算出する第３減算部と、
前記トルク偏差と、回転位相検出値と制御する加振周波数成分とトルクリプル周波数成分の次数ｎとを乗算したｎ次回転位相と、に基づいて周期外乱のｎ次周波数成分ベクトルを出力する周波数成分抽出器と、
前記ｎ次回転位相に基づいて、ｎ次回転周波数を出力する速度変換器と、
前記周期外乱のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗じて、操作量推定値のｎ次周波数成分ベクトルを求める逆モデル乗算部と、
前記操作量推定値のｎ次周波数成分ベクトルにオブザーバゲインを乗算して積分器に出力する乗算器と、
前記乗算器の出力を積分し、加振トルク制御器出力のｎ次周波数成分ベクトルを出力する積分器と、
前記加振トルク制御器出力のｎ次周波数成分ベクトルと前記ｎ次回転位相に基づいて前記加振トルク指令値を出力する補償信号合成部と、
を備えたことを特徴とする請求項１又は２記載のモータドライブシステム。
【請求項９】
前記次数ｎの異なる加振トルク制御器を複数台有し、
各々の加振トルク制御器の出力を合算した値を加振トルク指令値とすることを特徴とする請求項３〜８のうち何れかに記載のモータドライブシステム。
【請求項１０】
加振周波数成分抽出器に前記ｎ次回転位相が入力されず、前記ｎ次回転位相と非同期となる位相が入力される場合、並列段の制御構成の非同期加振周波数成分とｎ次回転周波数を各々監視しておき、これらの周波数が一致した場合、一致した段の制御構成の操作量推定値のｎ次周波数成分ベクトルを０とすることを特徴とする請求項９記載のモータドライブシステム。
【請求項１１】
前記加振トルク制御器の次数ｎは小数も含むことを特徴とする請求項３〜１０記載のモータドライブシステム。
【請求項１２】
供試体の入力側に接続されたエンジン爆発トルクを模擬する駆動モータと、
前記供試体の出力側に接続された車輪および路面の負荷を模擬する吸収モータと、
前記駆動モータのトルク検出値と回転位相検出値とに基づいて、第１インバータトルク指令値を出力する加振制御器と、
前記吸収モータのモータ回転数に基づいて、第２インバータトルク指令値を出力する速度制御器と、
前記第１インバータトルク指令値に基づいて、前記駆動モータを駆動する駆動モータインバータと、
前記第２インバータトルク指令値に基づいて、前記吸収モータを駆動する吸収モータインバータと、
を備えたことを特徴とする請求項１〜１１のうち何れかに記載のモータドライブシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、 ダイナモメータシステム、特に、ドライブトレインベンチにおける加振制御に係り、エンジンのように脈動が発生する複雑なトルク波形を、ダイナモメータの軸トルクで再現するための制御装置ならびに制御手法に関する。
【背景技術】
【０００２】
非特許文献１には、ドライブトレインベンチにおいて、共振抑制と低域トルク制御を実施しながら、加振振幅を制御する方式が開示されている。
【０００３】
特許文献１，２には、機械共振を考慮して周期外乱抑制制御を実施しながら、加振制御を行う手法が開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１１−１７６９５０号公報
【特許文献２】特開２０１１−１７６９５１号公報
【非特許文献】
【０００５】
【非特許文献１】秋山，小川，澤田，山本「ドライブトレインベンチの軸トルク加振制御」電気学会論文誌Ｃ，Ｖｏｌ．１３４，Ｎｏ．７，ｐｐ９０９−９１６，２０１４
【発明の概要】
【発明が解決しようとする課題】
【０００６】
ダイナモメータシステムのように、自動車や自動車部品を試験する計測制御システムでは、制御の安定性だけでなく高速，高応答，大容量，高精度などの高い要求を同時に満たさなければならない。これらの要求を満たすハードウェアの一つとして、高速・高応答を実現する低慣性化されたモータを用いる場合があり、供試体との間で多慣性共振系システムが形成され、計測値である軸トルク検出箇所においては、軸ねじれ共振現象が発生しやすい。
【０００７】
また、仕様・構造上の制約から、モータ極数の減少や電磁気的な設計の都合で、トルクリプルと呼ばれる周期的な外乱も発生しやすい。こうした共振・外乱は、試験・計測器として大きな問題となるため、共振抑制制御や外乱抑制制御が必須である。
【０００８】
特に、自動車の各種変速機，トルクコンバータ等の車両駆動系部品の試験を行うドライブトレインベンチでは、エンジンの爆発トルク波形をモータドライブシステムで模擬することが要求され、低域から高域まで幅広い周波数帯域での軸トルク加振制御が必要となる。
【０００９】
その結果、制御帯域に機械系で発生する共振周波数が含まれることになり、その機械共振は抑制しつつも所望する加振振動成分は発生させなければならないという高度な軸トルク制御が必要になる。
【００１０】
非特許文献１には、ドライブトレインベンチの軸トルク加振制御法が開示されている。ドライブトレインベンチでは、低域の共振点と高域の共振点が存在する。低域の共振は供試体であるトルクコンバータが持つ非線形ばねに起因し、高域の共振は駆動モータ，軸トルクメータおよびそれらをカップリングする機械系の剛性に起因する。
【００１１】
エンジンの加振トルクを模擬するため、軸トルク加振制御では、上述の機械共振による振動拡大を抑制した上で、所望の加振振幅が得られるように制御することが要求される。
【００１２】
そこで、非特許文献１では、供試体の非線形特性に大きく依存する低域の制御にはＩ−ＰＤ制御による軸トルク制御手法を適用し、ダイナモメータ設備に起因する高域の共振についてはμ設計法による共振抑制制御が適用されている。
【００１３】
これらを周波数帯域で分離して組み合わせた制御を実施した上で、軸トルクに所望の加振振幅が得られるように加振振幅の自動調整手法が開示されている。
【００１４】
しかしながら、本手法は単一の周波数成分、すなわち、正弦波による加振制御に限定されており、その正弦波の振幅のみが制御されている。したがって、正弦波の位相については制御されていない。
【００１５】
さらには、エンジンの実際の爆発トルク波形は複数の周波数成分が混在する歪み波形であるため、より正確にエンジン爆発トルクを模擬するためには、複数の周波数成分を同時に加振し、振幅だけでなく位相も合わせ込む必要がある。この非特許文献１の手法では、歪み成分を含んだエンジン爆発トルクを模擬する軸トルク加振制御は実現できない。
【００１６】
また、実際にはモータ構造に起因して発生するトルクリプルが周期外乱として存在するが、μ設計等の共振抑制制御のみでは十分に周期外乱を抑制できず、加振制御による周波数成分と混在して、軸トルク制御をより複雑にする。この非特許文献１では、こうした周期外乱による影響を考慮していないため、正弦波の加振トルク指令を与えたとしても、軸トルクにはエンジン波形とは異なる周波数成分の意図しない歪みが発生することになる。
【００１７】
特許文献１，２では、エンジン爆発トルクに含まれる歪み成分も考慮し、一般化周期外乱オブザーバによる複数の周期成分の制御が実現されている。機械共振の抑制については、予めシステム同定などによって得られたトルク伝達特性の逆特性を用いており、共振抑制テーブルとして実装されている。これにより、各周波数成分の加振成分を残した上で、機械共振や周期外乱（トルクリプル）による影響は除去することができる。
【００１８】
しかしながら、この特許文献１，２で開示されている手法は、予め設定した周波数成分のみを制御対象としているため、共振抑制テーブルについても、トルク指令値に加振周波数成分を重畳したことによる共振の影響を除去するのみである。
【００１９】
したがって、トルクセンサの検出ノイズ等のように非周期的な外乱が含まれると、軸トルク検出に機械共振の影響が残留する。また、共振抑制テーブルの他に、周期外乱を抑制するためのシステム同定テーブルも必要であるため、演算量やメモリの量が増加する。
【００２０】
以上示したように、モータドライブシステムの軸トルク加振制御において、共振・非周期外乱・周期外乱の影響をすべて除去しつつ、歪み成分を含んだエンジン加振トルク指令値に追従させることが課題となる。
【課題を解決するための手段】
【００２１】
本発明は、前記従来の問題に鑑み、案出されたもので、その一態様は、モータの軸トルクを、インバータを用いて制御するモータドライブシステムであって、トルク指令値とトルク検出値に基づいて、低域トルク制御器出力を出力する低域トルク制御器と、前記トルク指令値と前記トルク検出値と回転位相検出値とに基づいて加振トルク指令値を出力する加振トルク制御器と、前記低域トルク制御器出力と前記加振トルク指令値とを加算した補正トルク指令値と、前記トルク検出値と、に基づいてインバータトルク指令値を出力する高域共振抑制制御器と、を備えたことを特徴とする。 
【００２２】
また、その一態様として、前記高域共振抑制制御器は、前記補正トルク指令値が、μ設計制御器の指令値入力からμ設計制御器出力までの伝達特性を通った出力と、前記トルク検出値が、μ設計制御器の検出値入力からμ設計制御器出力までの伝達特性を通った出力と、を加算して、インバータトルク指令値を出力するμ設計制御器を有し、前記低域トルク制御器はＰＩＤ制御を行うことを特徴とする。 
【００２３】
また、その一態様として、前記加振トルク制御器は、前記トルク指令値と、前記回転位相検出値に制御する加振周波数成分とトルクリプル成分の次数ｎを乗算したｎ次回転位相と、に基づいて、前記加振トルク指令値のｎ次周波数成分ベクトルを出力する加振周波数成分抽出器と、前記トルク検出値と、前記ｎ次回転位相と、に基づいて、周期外乱のｎ次周波数成分ベクトルを出力するリプル抑制周波数成分抽出器と、前記ｎ次回転位相に基づいて、ｎ次回転周波数を出力する速度変換器と、前記加振トルク指令値のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗算し、加振による周期外乱指令値のｎ次周波数成分ベクトルを出力する第１逆モデル乗算部と、前記周期外乱のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗算して操作量推定値のｎ次周波数成分ベクトルを出力する第２逆モデル乗算部と、前記操作量推定値のｎ次周波数成分ベクトルから、加振トルク制御器出力のｎ次周波数成分ベクトルがローパスフィルタを介した値を減算して周期外乱推定値のｎ次周波数成分ベクトルを出力する第１減算器と、前記加振による周期外乱指令値のｎ次周波数成分ベクトルから、前記周期外乱推定値のｎ次周波数成分ベクトルを減算して加振トルク制御器出力のｎ次周波数成分ベクトルを出力する第２減算器と、前記加振トルク制御器出力のｎ次周波数成分ベクトルと前記ｎ次回転位相に基づいて前記加振トルク指令値を出力する補償信号合成部と、を備えたことを特徴とする。 
【００２４】
また、その一態様として、前記加振トルク制御器は、前記トルク指令値と前記トルク検出値とのトルク偏差を算出する第３減算部と、前記トルク偏差と、回転位相検出値と制御する加振周波数成分とトルクリプル周波数成分の次数ｎとを乗算したｎ次回転位相と、に基づいて周期外乱のｎ次周波数成分ベクトルを出力する周波数成分抽出器と、前記ｎ次回転位相に基づいて、ｎ次回転周波数を出力する速度変換器と、前記周期外乱のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗算して操作量推定値のｎ次周波数成分ベクトルを出力する逆モデル乗算部と、前記操作量推定値のｎ次周波数成分ベクトルと、加振トルク制御器出力のｎ次周波数成分ベクトルがローパスフィルタを介した値と、を加算して加振トルク制御器出力のｎ次周波数成分ベクトルを出力する加算器と、前記加振トルク制御器出力のｎ次周波数成分ベクトルと前記ｎ次回転位相に基づいて前記加振トルク指令値を出力する補償信号合成部と、を備えたことを特徴とする。 
【００２５】
また、その一態様として、前記加振トルク制御器は、前記トルク指令値と前記トルク検出値とのトルク偏差を算出する第３減算部と、前記トルク偏差と、回転位相検出値と制御する加振周波数成分とトルクリプル周波数成分の次数ｎとを乗算したｎ次回転位相と、に基づいて、周期外乱のｎ次周波数成分ベクトルを出力する周波数成分抽出器と、前記ｎ次回転位相に基づいて、ｎ次回転周波数を出力する速度変換器と、前記周期外乱のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗じて、操作量推定値のｎ次周波数成分ベクトルを求める逆モデル乗算部と、前記操作量推定値のｎ次周波数成分ベクトルを積分し、加振トルク制御器出力のｎ次周波数成分ベクトルを出力する積分器と、前記加振トルク制御器出力のｎ次周波数成分ベクトルと前記ｎ次回転位相に基づいて前記加振トルク指令値を出力する補償信号合成部と、を備えたことを特徴とする。  
【００２６】
また、その一態様として、前記加振トルク制御器は、前記トルク指令値と、前記回転位相検出値に制御する加振周波数成分とトルクリプル周波数成分の次数ｎを乗算したｎ次回転位相と、に基づいて、前記加振トルク指令値のｎ次周波数成分ベクトルを出力する加振周波数成分抽出器と、前記トルク検出値と、前記ｎ次回転位相と、に基づいて、周期外乱のｎ次周波数成分ベクトルを出力するリプル抑制周波数成分抽出器と、前記ｎ次回転位相に基づいて、ｎ次回転周波数を出力する速度変換器と、前記加振トルク指令値のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗算した加振による周期外乱指令値のｎ次周波数成分ベクトルを出力する第１逆モデル乗算部と、前記周期外乱のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗算して操作量推定値のｎ次周波数成分ベクトルを出力する第２逆モデル乗算部と、前記加振による周期外乱指令値のｎ次周波数成分ベクトルにオブザーバゲインを乗算して第２減算器に出力する第１乗算器と、前記操作量推定値のｎ次周波数成分ベクトルにオブザーバゲインを乗算して第１減算器に出力する第２乗算器と、前記第２乗算器の出力から、加振トルク制御器出力のｎ次周波数成分ベクトルがローパスフィルタを介した値を減算して周期外乱推定値のｎ次周波数成分ベクトルを出力する第１減算器と、前記第１乗算器の出力から、前記周期外乱推定値のｎ次周波数成分ベクトルを減算して加振トルク制御器出力のｎ次周波数成分ベクトルを出力する第２減算器と、前記加振トルク制御器出力のｎ次周波数成分ベクトルと前記ｎ次回転位相に基づいて前記加振トルク指令値を出力する補償信号合成部と、を備えたことを特徴とする。 
【００２７】
また、その一態様として、前記加振トルク制御器は、前記トルク指令値と前記トルク検出値とのトルク偏差を算出する第３減算部と、前記トルク偏差と、回転位相検出値と制御する加振周波数成分とトルクリプル周波数成分の次数ｎとを乗算したｎ次回転位相と、に基づいて周期外乱のｎ次周波数成分ベクトルを出力する周波数成分抽出器と、前記ｎ次回転位相に基づいて、ｎ次回転周波数を出力する速度変換器と、前記周期外乱のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗算して操作量推定値のｎ次周波数成分ベクトルを出力する逆モデル乗算部と、前記操作量推定値のｎ次周波数成分ベクトルにオブザーバゲインを乗算して加算器に出力する乗算器と、前記乗算器の出力と、加振トルク制御器出力のｎ次周波数成分ベクトルがローパスフィルタを通った値と、を加算して加振トルク制御器出力のｎ次周波数成分ベクトルを出力する加算器と、前記加振トルク制御器出力のｎ次周波数成分ベクトルと前記ｎ次回転位相に基づいて前記加振トルク指令値を出力する補償信号合成部と、を備えたことを特徴とする。 
【００２８】
また、その一態様として、前記加振トルク制御器は、前記トルク指令値と前記トルク検出値とのトルク偏差を算出する第３減算部と、前記トルク偏差と、回転位相検出値と制御する加振周波数成分とトルクリプル周波数成分の次数ｎとを乗算したｎ次回転位相と、に基づいて周期外乱のｎ次周波数成分ベクトルを出力する周波数成分抽出器と、前記ｎ次回転位相に基づいて、ｎ次回転周波数を出力する速度変換器と、前記周期外乱のｎ次周波数成分ベクトルに、前記ｎ次回転周波数に同期した単一の周波数ベクトルが適用された逆モデルを乗じて、操作量推定値のｎ次周波数成分ベクトルを求める逆モデル乗算部と、前記操作量推定値のｎ次周波数成分ベクトルにオブザーバゲインを乗算して積分器に出力する乗算器と、前記乗算器の出力を積分し、加振トルク制御器出力のｎ次周波数成分ベクトルを出力する積分器と、前記加振トルク制御器出力のｎ次周波数成分ベクトルと前記ｎ次回転位相に基づいて前記加振トルク指令値を出力する補償信号合成部と、を備えたことを特徴とする。 
【００２９】
また、その一態様として、前記次数ｎの異なる加振トルク制御器を複数台有し、各々の加振トルク制御器の出力を合算した値を加振トルク指令値とすることを特徴とする。
【００３０】
また、その一態様として、加振周波数成分抽出器に前記ｎ次回転位相が入力されず、前記ｎ次回転位相と非同期となる位相が入力される場合、並列段の制御構成の非同期加振周波数成分とｎ次回転周波数を各々監視しておき、これらの周波数が一致した場合、一致した段の制御構成の操作量推定値のｎ次周波数成分ベクトルを０とすることを特徴とする請求項９記載のモータドライブシステム。 
【００３１】
また、その一態様として、前記加振トルク制御器の次数ｎは小数も含むことを特徴とする。 
【００３２】
また、その一態様として、供試体の入力側に接続されたエンジン爆発トルクを模擬する駆動モータと、前記供試体の出力側に接続された車輪および路面の負荷を模擬する吸収モータと、前記駆動モータのトルク検出値と回転位相検出値とに基づいて、第１インバータトルク指令値を出力する加振制御器と、前記吸収モータのモータ回転数に基づいて、第２インバータトルク指令値を出力する速度制御器と、前記第１インバータトルク指令値に基づいて、前記駆動モータを駆動する駆動モータインバータと、前記第２インバータトルク指令値に基づいて、前記吸収モータを駆動する吸収モータインバータと、を備えたことを特徴とする。
【発明の効果】
【００３３】
本発明によれば、モータドライブシステムの軸トルク加振制御において、共振・非周期外乱・周期外乱の影響をすべて除去しつつ、歪み成分を含んだエンジン加振トルク指令値に追従させることが可能となる。
【図面の簡単な説明】
【００３４】
【図１】ドライブトレインレンチの基本構成図。
【図２】インバータトルク指令値τinv*→トルク検出値τdetの周波数伝達特性の例を示す図。
【図３】供試体の非線形ばねの特性の例を示す図。
【図４】制御対象の３慣性系近似モデルを示す図。
【図５】実施形態１の基本制御構成図。
【図６】一般化プラントの構成例を示す図。
【図７】μ設計制御器の特性図。
【図８】μ設計制御器の構成図。
【図９】高域共振抑制のみを実施した場合の周波数伝達特性図。
【図１０】低域トルク制御器の設計図。
【図１１】低域トルク制御のみ実施した場合の周波数伝達特性図。
【図１２】加振トルク指令値τpd*→トルク検出値τdetの周波数伝達特性図（低域トルク制御ＯＮ，高域共振抑制制御器ＯＮ）。
【図１３】加振トルク制御器の基本構成図。
【図１４】実施形態１の効果例を示す図。
【図１５】実施形態２の加振トルク制御器の構成図。
【図１６】実施形態３の加振トルク制御器の構成図。
【図１７】実施形態４の加振トルク制御器の構成図。
【図１８】一般化周期外乱オブザーバのモデル誤差に対するロバスト安定性を示す図。
【図１９】１気筒失火モードのトルク波形追従制御例を示す図。
【発明を実施するための形態】
【００３５】
本願は、「歪み波形を含む任意波形追従制御」に関する提案であるため、図１の基本構成には限らないが、実施形態１〜７では、非特許文献１に記載されているドライブトレインベンチと同様の基本構成において、軸トルク加振の任意波形追従制御を実施する例で課題の解決法と適用効果を説明する。
【００３６】
図１は、ドライブトレインベンチの基本構成図である。供試体１はトルクコンバータとし、その入力側の駆動モータ２において、エンジン爆発トルクを模擬する。出力側は、車輪および路面の負荷を模擬する吸収モータ３であり、本構成では回転数制御を行う。供試体１と駆動モータ２の間には軸トルクメータ４を設け、検出されたトルク検出値τdetが、エンジン爆発トルクを模擬したトルク指令値τ*に追従するように制御する。
【００３７】
加振制御器５は、加振周波数成分を含んだトルク指令値τ*に基づいて、加振制御を実施する。本願発明では、非特許文献１と異なり、加振周波数成分の位相も検知するため、回転位置センサ６による回転位相検出値θを用いる。
【００３８】
速度制御器７では、吸収モータ３で回転数制御を行うために、回転位置センサ８のモータ回転数ωmが速度指令値ωm*に追従するように制御する。速度制御器７は一般的なＰＩＤ制御等で実施すれば良い。
【００３９】
加振制御器５，速度制御器７が出力するインバータトルク指令値τinv1*，τinv2*に基づき、駆動モータインバータＩＮＶ１，吸収モータインバータＩＮＶ２を介して駆動モータ２，吸収モータ３を駆動する。
【００４０】
なお、駆動モータインバータＩＮＶ１，吸収モータインバータＩＮＶ２は、直流電力を交流電力に変換する電力変換器である。駆動モータインバータＩＮＶ１，吸収モータインバータＩＮＶ２の交流出力端子を駆動モータ２，吸収モータ３の端子に接続する。駆動モータインバータＩＮＶ１，吸収モータインバータＩＮＶ２は、駆動モータインバータＩＮＶ１，吸収モータインバータＩＮＶ２の交流出力電流（すなわちモータ電流）を制御することでモータの軸トルク（トルク検出値）τdetをインバータトルク指令値τinv1*，τinv2*に制御する機能を持つ。
【００４１】
図２は、駆動モータインバータＩＮＶ１ヘのインバータトルク指令値τinv1*からトルク検出値τdetまでの周波数伝達特性の例である。本願発明では、図１の構成を等価的に３慣性系に近似しており、低周波数領域で変動している共振点は、供試体１であるトルクコンバータの非線形ばねに起因したものである。
【００４２】
参考に非線形ばね特性の例を図３に示す。高周波数領域の共振点は、駆動モータ２・軸トルクメータおよびカップリング等の機械設備の軸ねじれ剛性に起因したものである。
【００４３】
このように、本願発明で説明するドライブトレインベンチでは、低域に非線形な共振特性を持ち、高域に高い共振特性を持つ３慣性系となり、トルク指令値τ*に対して、単純に加振トルクを与えた場合は共振が発生し、所望の加振波形とはならないことを意味する。
【００４４】
なお、広く自動車で用いられている４サイクルエンジンでは、その爆発トルクは、気筒数×０．５×回転数の周波数成分で大きな振動トルクが発生する。例えば、エンジンの回転数を６００ｍｉｎ-1〜６０００ｍｉｎ-1とし、３気筒から８気筒エンジンまで模擬する場合、１５〜４００Ｈｚの帯域で加振制御を実現する必要がある。
【００４５】
また、歪み成分を考慮すれば４００Ｈｚを超えてさらに高い帯域まで必要となるため、機械構成上の限界まで可能な限り加振周波数成分を追従させる必要がある。 
【００４６】
図４に、本構成を３慣性系に近似した場合の制御対象のブロック図を示す。 図４の符号は以下の通りである。
【００４７】
Ｊ１：駆動モータの慣性モーメント，Ｊ２：供試体・カップリングの慣性モーメント，Ｊ３：吸収モータの慣性モーメント，Ｋ１２：カップリングおよび軸トルクメータの軸ねじれ剛性，Ｋ２３およびΔＫ２３：供試体の軸ねじれ剛性（非線形特性を有する：図３），Ｃ２３：供試体の損失，Ｇinv（ｓ）：駆動側インバータの応答伝達関数，Ｇtm（ｓ）：軸トルク検出の応答伝達関数，Ｇθ（ｓ）：回転位相検出の応答伝達関数，ｓ：ラプラス演算子。
【００４８】
前述の通り、供試体１であるトルクコンバータは非線形ばね特性を持つため、Ｋ２３およびΔＫ２３を用いて表現する。また、インバータ，トルクメータ検出，位相検出には、応答遅れや無駄時間が存在するため、近似した応答伝達関数を挿入している。
【００４９】
図４に示した近似モデルのインバータトルク指令値τinv*からトルク検出値τdetまでの周波数伝達特性は図２に示したとおりであるが、加振の周波数帯域で常に０ｄＢ，すなわちトルク指令値τ*からトルク検出値τdetへの伝達特性が１になるようなコントローラを設計することが望ましい。そのためには、低域と高域における共振抑制制御が必要となる。また、振幅特性が０ｄＢより低下している帯域では増幅する必要がある。
【００５０】
［実施形態１］
本実施形態１では、非特許文献１と同様に、低域と高域の共振周波数帯を切り分けて、高域の共振抑制制御と低域の定常トルク制御を行う手法を基本とし、非特許文献１では実現できない「歪み成分を含む加振波形追従制御」を実現するための手法として、一般化周期外乱オブザーバを併用する手法を説明する。
【００５１】
図５は、本実施形態１の基本となる制御構成図である。 図５の符号は以下の通りである。
【００５２】
τ*：トルク指令値，τdc*：低域トルク制御器出力，τpd*：加振トルク指令値，τr*：補正トルク指令値，τinv*：インバータトルク指令値，τdet：トルク検出値，θ：回転位相検出値。
【００５３】
図５に示すように、本実施形態１におけるモータの軸トルクをインバータを用いて制御するモータドライブシステムは、低域トルク制御器９により、トルク指令値τ*とトルク検出値τdetに基づいて低域トルク制御器出力τdc*を出力する。また、加振トルク制御器１１により、トルク指令値τ*とトルク検出値τdetと回転位相検出値θとに基づいて加振トルク指令値τpd*を出力する。高域共振抑制制御器１０は、低域トルク制御器出力τdc*と前記加振トルク指令値τpd*とを加算した補正トルク指令値τr*と、前記トルク検出値τdetと、に基づいてインバータトルク指令値τinv*を出力する。
【００５４】
本実施形態１における「制御対象」は、図４に示した３慣性系近似モデルに該当する。
【００５５】
図５に示すとおり、本実施形態１の機能は「低域トルク制御器９」「高域共振抑制制御器１０」および一般化周期外乱オブザーバによる「加振トルク制御器１１」の３つの機能が存在する。以下、各機能について順次説明する。
【００５６】
・第１の機能「高域共振抑制制御器１０」
図４に示した３慣性系近似モデルでは、低域の共振点と高域の共振点が存在する。仮に、全周波数帯域で共振抑制する包括的なコントローラを設計した場合、保守的な設計となりやすく、特に低域側の非線形性も考慮した場合はロバスト制御性能を満たす設計自体が困難となる。
【００５７】
したがって、このような系に対しては、低域側と高域側の制御系を分離することが有効である。本実施形態１では、高域側の共振抑制制御をロバスト制御手法のひとつであるμ設計制御器で実現し、低域側の供試体１に起因した非線形な共振特性に対しては積極的な共振抑制制御は行わずに定常トルクの追従制御を実現する手法を組み合わせる。
【００５８】
まず、第１の機能である高域共振抑制制御器１０について説明する。高域の共振周波数は、図２の特性から、図４における駆動モータの慣性モーメントＪ１と供試体・カップリングの慣性モーメントＪ２による２慣性共振系の共振周波数に等価できる。したがって、高域の共振周波数は（１）式で計算される共振周波数ｆrHとほぼ一致する。
【００５９】
【数１】
2018007287.tif 000003
【００６０】
そこで、高域の共振を抑制するために近似された２慣性系を考慮し、高域共振抑制制御器の設計に用いる一般化プラントを図６のように構成する。本実施形態１では、ロバスト制御手法のひとつであるμ設計を用いた例で説明するが、Ｈ∞制御などのその他の一般的な共振抑制制御手法を用いても良い。図６の符号は以下の通りである。
【００６１】
Δτ：トルク偏差，Ｇtm（ｓ）：トルクメータ等の検出応答伝達関数，Ｇinv（ｓ）：インバータ応答伝達関数，ｄ１：外乱（周期外乱を含む），ｄ２およびｒ：μ設計制御器の指令値入力，ｄ３：トルク検出ノイズ，ｚ：定常トルク誤差の評価出力，ｗ：定常トルク誤差による外乱入力，ｕ：μ設計制御器出力，ｙ：μ設計制御器の検出器入力，ｅ１：トルク検出値の評価出力，ｅ２：インバータトルク指令の評価出力，ｅ３：μ設計制御器ゲインの評価出力，Ｗn（ｓ）：外乱ｄ１に対する重み関数，Ｗu（ｓ）：μ設計制御器出力ｕに対する重み関数，Ｗe（ｓ）：μ設計制御器ゲインに対する重み関数。
【００６２】
μ設計によるコントローラ設計に際し、機械系パラメータ等の摂動を個別に考慮しても良いが、実用上は物理モデル（バネ・マス要素）を明示的に同定する作業を省略し、トルク入出力の周波数伝達特性から簡易的に同定することが多い。
【００６３】
したがって、ここではインバータから入力されるトルク偏差Δτを摂動項として見積もり、ロバスト制御性能を担保する。
【００６４】
低周波数領域の定常トルク追従制御は後述する低域トルク制御器９で実現するため、高周波数領域の高域共振抑制制御器１０では、低周波数領域におけるコントローラゲインを落とした設計を行う。
【００６５】
すなわち、μ設計制御器の指令値入力ｒからμ設計制御器出力ｕのゲインにおいて、μ設計制御器ゲインに対する重み関数Ｗe（ｓ）を介してμ設計制御器ゲインの評価出力ｅ３を設定しているが、低域トルク制御器９との制御干渉を防ぐために、μ設計制御器ゲインに対する重み関数Ｗe（ｓ）は低域に対して重みを置く。
【００６６】
また、μ設計制御器出力ｕに対する重み関数（インバータトルク指令ｕに対する重み関数）Ｗu（ｓ）は、インバータトルク指令の高周波ゲインを低減するために高周波領域で重み付けする。外乱ｄ１に対する重み関数Ｗn（ｓ）は、周期外乱や非周期的外乱の抑圧性能を向上するように、共振周波数付近に重みを置く。
【００６７】
図６のように構成した一般化プラントにおいて、Ｄ−Ｋイテレーションを実施して得られたμ設計制御器の特性例を図７に示す。伝達特性Ｃtm（ｓ）のゲイン特性を見て分かるように、μ設計制御器ゲインに対する重み関数Ｗe（ｓ）の効果により、トルク検出値τdetに対しては低域（低周波数領域）側で反応しないように，低域のゲインが低下していることが分かる。これにより、後述する低域トルク制御器９との制御干渉を防止できる。
【００６８】
また、伝達特性Ｃref（ｓ）のゲイン特性を見て分かるように共振周波数付近のゲインを低減するように設計しているため、加振トルクの周波数成分による共振を抑制できる。
【００６９】
図７の上段はゲイン線図，下段は位相線図を示している。左側はμ設計制御器の指令値入力ｒからμ設計制御器出力ｕの伝達特性Ｃref（ｓ），右側はμ設計制御器の検出値入力ｙからμ設計制御器出力ｕの伝達特性Ｃtm（ｓ）を示している。
【００７０】
上記で設計された伝達特性Ｃref（ｓ）と伝達特性Ｃtm（ｓ）を有するμ設計制御器は、図８に示す構成で、図５の高域共振抑制制御器１０の部分に実装する。
【００７１】
補正トルク指令値τr*が伝達特性Ｃref（ｓ）を通った出力と、トルク検出値τdetが伝達特性Ｃtm（ｓ）を通った出力と、を加算して、インバータトルク指令値τinv*を生成する。
【００７２】
図９に、高域共振抑制制御のみを実装した場合の、トルク指令値τ*→トルク検出値τdetの周波数伝達特性を示す。高域（高周波数領域）のゲインが０ｄＢ以下に抑制されていることが分かる。
【００７３】
・第２の機能「低域トルク制御器９」
次に、第２の機能として図５の低域トルク制御器９について説明する。低域の共振周波数は、供試体１の非線形ばね特性に依存して変動する。この非線形特性を考慮してロバストな共振抑制制御を実施する手法も考えられるが、供試体の非線形特性に依存した設計となる。実際の供試体１は交換されて試験されることが多く、供試体特性に依存したコントローラを設計した場合は、その都度、再設計やパラメータの再調整が必要となってしまう問題がある。
【００７４】
したがって、本実施形態１では、供試体１が変更されてもその特性に依存しない設計を目指し、低域トルク制御器９では積極的な共振抑制を行わずに定常トルクを安定に追従制御することのみを目的とする。
【００７５】
図４に示した３慣性系近似モデルにおいて、低域の共振周波数に着目すると、図２の特性から、Ｊ１＋Ｊ２，Ｋ２３，Ｊ３による２慣性共振系に近似できる。すなわち、（２）式よって低域の共振周波数ｆrLが概ね分かる。ただし、軸ねじれ損失Ｃ２３の影響については省略する。
【００７６】
【数２】
2018007287.tif 000004
【００７７】
低域の制御対象を２慣性系に近似した場合、低域トルク制御器９は図１０のように近似して設計できる。ここでは一例として、ＰＩＤ制御の比例微分先行型Ｉ−ＰＤ制御を用いており、比例項と微分項に対して１次のローパスフィルタＧF1（ｓ）を適用した形式を採用する。もちろん、図１０の構成に限らず、そのほかのＰＩＤ制御型を用いても実現可能である。
【００７８】
また、高域の共振点による制御スピルオーバを回避するために、Ｉ−ＰＤ制御器の出力に２次のローパスフィルタＧF2（ｓ）を用いている。この２次のローパスフィルタＧF2（ｓ）のカットオフ周波数は、低域のトルク制御特性に大きな影響を与えず、かつ、高域の共振特性を遮断できるカットオフ周波数に設定する。これにより、前述した第１の機能である高域の共振抑制制御系との分離を行って制御干渉を防止する。
【００７９】
図１０の近似された構成におけるトルク指令値τ*からトルク検出値τdetまでの閉ループ伝達特性は４次系となるため、Ｉ−ＰＤ制御のパラメータＫp，Ｋi，Ｋdは、４次標準形の極と一致するモデルマッチング手法などを用いて算出すれば良い。例として、４次バターワース型の特性に極配置する場合は、（３）式のように算出される。ＰＩＤパラメータの決定方法はこれに限らず、種々の他の手法を用いて調整しても良い。
【００８０】
【数３】
2018007287.tif 000005
【００８１】
ａ１，ａ２，ａ３は、 標準形モデルの極配置によって決定される特性方程式の各次数の係数であり、バターワース標準形の場合は、ａ１＝２．６１３１，ａ２＝３．４１４２，ａ３＝２．６１３１である。
【００８２】
また、ｋはトルク追従制御の制御応答を決定するパラメータであり、低域の共振周波数ｆrLに対する係数で指定する。このときの応答周波数ωｃは、ωｃ＝ｋ×（２×π×ｆrL）で決定づけられる。なお、ωｆ１は比例項と微分項に係るローパスフィルタのカットオフ周波数であり、（３）式の通り決定できる。以上のように構成した低域トルク制御を実施した場合のトルク指令値τ*→トルク検出値τdetの周波数伝達特性を図１１に示す。
【００８３】
図１１を見て分かるように、低域トルク制御器９では低域の共振特性をカットオフしているが、所望の指令値振幅が得られるのは低域の共振周波数よりもさらに低い領域のみとなる。
【００８４】
すなわち、定常トルクや低域共振周波数よりも低い領域では安定な追従制御が可能であるが、高い周波数の加振信号をトルク指令値に与えても、反応できない制御器となる。したがって、加振制御については別途、図５に示したように低域トルク制御器９を介さない加振トルク制御器１１を用いて、低域トルク制御器の出力部に重畳する必要がある。
【００８５】
・第３の機能「加振トルク制御器１１」
前述したとおり、低域トルク制御器９では定常トルクを含む低周波領域の追従制御のみを行い、高域共振抑制制御器１０は高周波領域にある共振特性を減衰させることのみに寄与するものである。したがって、所望する広範な周波数帯域で加振制御を実現するために図５の加振トルク制御器１１が必要となる。
【００８６】
加振トルク制御器１１は、加振周波数成分を含むトルク指令値τ*とトルク検出値τdetおよび加振周波数成分を生成するための基準位相θ（モータの回転位相検出値）を入力とする。加振トルク制御器１１から出力される加振トルク指令値τpd*は、低域トルク制御器出力τdc*に重畳し、補正トルク指令値τr*とする。
【００８７】
このとき、加振トルク制御器１１の加振トルク指令値τpd*からトルク検出値τdetまでの伝達特性は、低域トルク制御器９と高域共振抑制制御器１０を含む制御対象の閉ループ伝達特性となり、図１２のようになる。加振周波数帯域において、図１２の振幅（ゲイン）特性が０ｄＢ一定，位相特性が０ｄｅｇ一定になることが望ましいが、実際にそのような周波数特性を形成するコントローラは困難である。
【００８８】
したがって、加振トルク制御器１１では、図１２の伝達特性を考慮して、加振トルクの各周波数成分のゲインおよび位相を調整する必要がある。加振する周波数帯域（１５〜４００Ｈｚ）において、図１２の逆特性を用いて指令値を生成することも考えられるが、供試体１には非線形ばね特性があり、動作状態で変化する。したがって、単純に逆特性を介して指令値を生成する手法では、所望の振幅・位相の加振トルク波形を得ることは難しい。
【００８９】
そこで、一般化周期外乱オブザーバを用いた加振トルク波形の自動調整手法を提案する。図１３に、本実施形態１における加振トルク制御器１１の基本構成図を示す。図１３の符号は以下の通りである。
【００９０】
τ*：トルク指令値（加振周波数成分を含む），τdet：トルク検出値，θ：回転位相検出値，ｎ：次数（制御する加振周波数成分とトルクリプル周波数成分の次数を指定する），ωm：モータ回転数，τrpd*：加振トルク指令値 ，Ｔn*：加振トルク指令値のｎ次周波数成分ベクトル，Ｔn：周期外乱（トルクリプル）のｎ次周波数成分ベクトル，Ｕn＾：操作量推定値のｎ次周波数成分ベクトル（周期外乱を含む推定値），Ｄn*：加振による周期外乱指令値のｎ次周波数成分ベクトル，Ｄn＾：周期外乱推定値のｎ次周波数成分ベクトル，Ｔpdn*：加振トルク制御器出力のｎ次周波数成分ベクトル，ＧF（ｓ）：周波数成分抽出用のローパスフィルタ。
【００９１】
一般化周期外乱オブザーバは、特定の周波数成分に着目して周期外乱を抑制する制御系であるが、本実施形態１では所望の周期的振動を発生させる方式として応用する。本方式は、特定の周波数成分のみに寄与する制御系であり、まず、加振周波数成分抽出器１２で加振したい周波数成分を抽出し、リプル抑制周波数成分抽出器１３で抑制したい周期外乱（トルクリプル）の周波数成分を抽出する。
【００９２】
すなわち、加振周波数成分抽出器１２において、トルク指令値τ*と、回転位相検出値θに次数ｎを乗算したｎ次回転位相ｎθと、に基づいて、加振トルク指令値のｎ次周波数成分ベクトルＴn*を出力する。また、リプル抑制周波数成分抽出器１３において、トルク検出値τdetと、ｎ次回転位相ｎθと、に基づいて、周期外乱のｎ次周波数成分ベクトルＴnを出力する。
【００９３】
これらは、モータ回転数に同期して発生するｎ次周波数成分となるため、回転位相検出値θのｎ倍のｎ次回転位相ｎθを用いて以下のように抽出する。ここで、ｎ次回転位相ｎθの周波数成分に同期した座標をｄnｑn回転座標系として定義し、ｄnを複素ベクトルの実部，ｑnを虚部の軸と置く。
【００９４】
【数４】
2018007287.tif 000006
【００９５】
（４）式はトルク指令値τ*に含まれる加振周波数成分を抽出し、（５）式はトルク検出値τdetに含まれる周期外乱すなわちトルクリプルの周波数成分を抽出するためのものである。厳密なフーリエ変換を用いても良いが、本実施形態１では演算器への実装の容易性を考慮し、ローパスフィルタＧF（ｓ）で周波数成分を抽出している。
【００９６】
速度変換器１４では、ｎ次回転位相ｎθを微分して、ｎ次回転周波数ｎ・ωmを算出する。
【００９７】
次に、逆モデルＱnについて説明する。ｄnｑn回転座標系では特定周波数成分のみに寄与する制御系となるため、加振トルク制御器１１における制御対象モデルは、１次元複素ベクトルで表現できる。ここで、ｎ次回転位相ｎθの周波数 成分に同期した制御対象モデルをＰnとする。ただし、Ｐn＝Ｐdn＋ｊＰqnである。
【００９８】
前述したとおり、加振トルク制御器１１にとっての制御対象システムは、図１２に示した加振トルク指令値τpd*→トルク検出値τdetの周波数伝達特性となるが、このグラフにおけるｎ次回転位相ｎθの周波数成分に同期した振幅・位相特性を、複素ベクトルとして抽出したものが制御対象モデルＰnとなる。したがって、モータ回転数ωmおよび次数ｎに応じて制御対象モデルＰnは変化することを意味する。
【００９９】
例えば、図１２の１〜１０００Ｈｚまでの振幅・位相特性を１Ｈｚ毎に区分した場合、１０００個の複素ベクトルが構成され、その中からモータ回転数によって変化するｎ次回転位相ｎθの周波数成分に同期したものを１つ選択し、制御対象モデルＰnに適用すれば良い。言い換えれば、制御対象モデルＰnはモータ回転数に依存した一種のゲインスケジューリング機能に相当する。
【０１００】
こうして定義したｄnｑn回転座標系における制御対象モデルＰnは、加振ないし抑制したい周波数成分に応じて変化するため、その逆モデルＱnについてもモータ回転数に応じて変化させなければならない。そこで、図１３に示したとおり、速度検出器によって算出されたｎ次回転周波数ｎ・ωmに基づいて、ｎ次回転位相ｎθの周波数成分に同期した逆モデルを選択する。逆モデルＱnは、以下の（６）式で示される。
【０１０１】
【数５】
2018007287.tif 000007
【０１０２】
本制御の目的は、制御対象の出力である軸トルクメータで所望の加振トルクとトルクリプル抑制を実現することであり、制御対象の伝達特性（制御対象モデル）Ｐnを考慮して制御対象の入力（操作量）を決定しなければならない。そこで、図１３のように逆モデルＱnの特性を用いることで、加振トルク指令値のｎ次周波数成分ベクトルＴn*から制御対象モデルＰnの入力部の加振による周期外乱のｎ次周波数成分ベクトルＤn*を算出する。加振による周期外乱のｎ次周波数成分ベクトルＤn*の算出式は、（６−２）式となる。
【０１０３】
【数６】
2018007287.tif 000008
【０１０４】
同様に、周期外乱のｎ次周波数成分ベクトルＴnから操作量推定値のｎ次周波数成分ベクトル（周期外乱を含む推定値）Ｕn＾を推定する。操作量推定値のｎ次周波数成分ベクトルＵn＾の算出式は（６−３）式となる。
【０１０５】
【数７】
2018007287.tif 000009
【０１０６】
ここで、操作量推定値のｎ次周波数成分ベクトルＵn＾には周期外乱であるトルクリプルの成分も含んで推定されているため、制御対象モデルＰnに入力する操作量である加振トルク制御器出力のｎ次周波数成分ベクトルＴpdn*を差し引き、周期外乱推定値のｎ次周波数成分ベクトルＤn＾を推定する。この原理は、従来から広く利用されている外乱オブザーバの手法を踏襲したものである。
【０１０７】
なお、操作量推定値のｎ次周波数成分ベクトルＵn＾から加振トルク制御器出力のｎ次周波数成分ベクトルＴpdn*を差し引く際、加振周波数成分抽出器１２，リプル抑制周波数成分抽出器１３で用いられているローパスフィルタＧF（ｓ）の応答遅れと同期させる目的で、ローパスフィルタＧF（ｓ）を介した加振トルク制御器出力のｎ次周波数成分ベクトルＴpdn*を差し引いている。
【０１０８】
そして、加振による周期外乱指令値のｎ次周波数成分ベクトルＤn*から、周期外乱推定値のｎ次周波数成分ベクトルＤn＾を、さらに差し引くことにより、加振トルク制御器出力のｎ次周波数成分ベクトルＴpdn*が算出される。以上により、加振トルク指令に含まれる周波数成分を残しつつ、周期外乱（トルクリプル）に起因した振動成分は除去することが可能となり、所望する振動成分のみを軸トルク検出部で発生させることができる。
【０１０９】
加振トルク制御器出力のｎ次周波数成分ベクトルＴpdn*は、補償信号合成部１９において、（７）式に基づいてｄnｑn回転座標系から元の時間波形に復元する。なお、ｎに複数の次数を設定して並列に構成することも可能であり、これら各次数の周波数成分を合算することにより、加振トルク指令値τpd*が合成できる。
【０１１０】
【数８】
2018007287.tif 000010
【０１１１】
以上、図５の構成に含まれる３つの基本機能について説明した。この３つの機能を同時に動作させることにより、以下の効果が得られる。
【０１１２】
「低域のトルク追従制御」と「高域の共振抑制制御」と「軸トルク加振制御」が干渉することなく同時に機能し、複数の周波数成分を含んだ歪み波形となるエンジン爆発トルクを模擬することが可能となる。
【０１１３】
加振トルク指令の振動成分は残しつつも、障害となるトルクリプルによる周期外乱のみ抑制することができる。
【０１１４】
高域の共振抑制制御を加えたことで、共振特性による振幅・位相の急激な変化を緩やかにできるため、加振トルク制御器１１における一般化周期外乱オブザーバの逆モデル特性も緩やかになる。このことは、可変速運転時における逆モデルの抽出において、その特性変化が緩やかになることを意味しており、特に共振周波数を交差する場合やモデル化誤差がある場合において、ロバスト安定性の向上に大きく寄与する。
【０１１５】
高域の共振抑制制御は非周期的な外乱も抑制する効果があり、一般化周期外乱オブザーバによる周期外乱抑制効果に加えて、非周期外乱の抑制も同時に実現できる。
【０１１６】
低域トルク追従制御と、高域の共振抑制制御と、周期成分の制御系を分離して設計することで、それぞれが単独では実現困難な包括的な「共振抑制」，「非周期外乱抑制」，「周期外乱抑制」，「加振制御」を同時に実現できる。また、それらが制御干渉しない。 
【０１１７】
図１４に、本実施形態１における効果の例を示す。 図１４上段は、従来の加振制御方式の軸トルク波形である。非特許文献１と同様に「低域トルク制御」および「高域共振抑制制御」を適用しているが、本実施形態１で提案した「加振トルク制御器１１」は機能させずに、歪み成分を含む加振トルク指令値をそのまま加振トルク指令値τpd*に適用した場合である。加振トルク指令値と軸トルク検出値は一致せず、所望する波形が軸トルクでは得られない。
【０１１８】
非特許文献１では、加振振幅の自動調整手法も提案されているが、加振振幅の大きさのみを追従させる方式であり、単一周波数成分の正弦波波形による加振制御を前提としている。したがって、非特許文献１の加振振幅制御を実施しても、歪み波形の形状までは追従制御できず、位相も一致しない。また、トルクリプルによる周期外乱の影響も考慮されていないため、図１４の上段の波形のように意図しない歪み成分が残留することになる。
【０１１９】
図１４下段は、本実施形態１の加振トルク制御器１１によるエンジン加振トルク波形追従制御を実施した場合のトルク波形である。歪み波形を含んだ加振トルク指令値に一致しており、振幅・位相の両方を追従させることができる。また、一般化周期外乱オブザーバによって、加振制御と同時に、周期外乱であるトルクリプルを推定・除去しているため、前述したような周期外乱による意図しない歪みも除去されている。
【０１２０】
なお、モータ回転数（エンジン回転数に相当）やトルクの大きさなど、動作状態が変化した場合であっても、加振トルク制御器が機能して、自動的に形状を追従できる。このときの過渡応答性は一般化周期外乱オブザーバの速応性と同−であり、周波数成分抽出に用いるローパスフィルタＧF（ｓ）によって決定づけられる。数値例として、過渡変化に対しても０．３秒程度で追従可能である。
【０１２１】
なお、実施形態１を図１のモータドライブシステムに適用する場合は、図１の加振制御器５は、図５の低減トルク制御器９と高域共振抑制制御器１０と加振トルク制御器１１によって構成される。
【０１２２】
［実施形態２］
実施形態１の図１３の構成では、トルク指令値τ*に含まれる周波数成分と、トルク検出値τdetに含まれる周期外乱周波数成分を別々に抽出していたが、本実施形態２ではこれらを区別せず、より簡素な構成で加振周波数成分と周期外乱成分を同時に制御する手法を提供する。
【０１２３】
図１５は、本実施形態２における加振トルク制御器１１の構成図である。その他の制御機能については、実施形態１と同様である。
【０１２４】
加振周波数成分を含むトルク指令値τ*とトルク検出値τdetのトルク偏差Δτを取り、その偏差に含まれる周波数成分がゼロとなれば、周期的な振動成分に関しては加振トルク指令値と軸トルク検出値が一致することを意味する。
【０１２５】
したがって、本実施形態２では、加振周波数成分とトルクリプルによる周期外乱周波数成分を区別することなく、統一的にトルク偏差Δτの周波数成分を抽出し、その偏差がなくなるように一般化周期外乱オブザーバを動作させる。実施形態１の図１３においては、指令側の周波数成分ベクトルとして加振による周期外乱指令値のｎ次周波数成分ベクトルＤn*を用いていたが、本実施形態２では偏差をなくすことが目的であるため、加振による周期外乱指令値のｎ次周波数成分ベクトルＤn*＝０とすれば良い。
【０１２６】
図１３において、加振による周期外乱指令値のｎ次周波数成分ベクトルＤn*＝０とした場合の制御ブロック図は、図１５のように等価変換可能であり、比較的簡素な構成で実施形態１と同じ制御効果が得られる。実施形態１の図１３と比べると、加振周波数成分抽出器１２を１つと逆モデルを乗算する逆モデル乗算器１５が削減されるので、マイクロプロセッサ等の計算機の演算負荷を軽減できる。
【０１２７】
また、トルク偏差Δτを用いることにより、定常的には低域トルク制御器９の追従動作によって、直流成分が除去される。このことによる副次的な効果を以下に説明する。
【０１２８】
周波数成分抽出器２２において、仮に直流成分Ｔｄｃが含まれていたとすると、周波数成分ＴdnおよびＴqnは（８）式のように展開される。
【０１２９】
【数９】
2018007287.tif 000011
【０１３０】
（８）式のＴdnｃｏｓ２ｎθ＋Ｔqnｓｉｎ２ｎθ＋２Ｔdcｃｏｓｎθ，Ｔqnｃｏｓ２θ＋Ｔdnｓｉｎ２ｎθ＋２Ｔdcｓｉｎｎθを、ローパスフィルタＧF（ｓ）で除去することにより、Ｔdn，Ｔqnを得ることができる。２ｎθ成分とｎθ成分が含まれており、これらの周波数成分がローパスフィルタＧF（ｓ）のカットオフ周波数に比べて十分に高い場合は良好に周波数成分を抽出できる。
【０１３１】
一方、モータ回転数が低く、ｎ次回転位相ｎθの周波数も低くなる場合は、ローパスフィルタＧF（ｓ）のカットオフ周波数に近づき、Ｔdnｃｏｓ２ｎθ＋Ｔqnｓｉｎ２ｎθ＋２Ｔdcｃｏｓｎθ，Ｔqnｃｏｓ２θ＋Ｔdnｓｉｎ２ｎθ＋２Ｔdcｓｉｎｎθの周波数成分の影響がＴdn，Ｔqnに現れることを意味する。これらは、一般化周期外乱オブザーバのｄnｑn回転座標系においては外乱となり、制御系を不安定化する要因となる。
【０１３２】
したがって、ローパスフィルタＧF（ｓ）のカットオフ周波数は十分に低い値に設定する必要がある。しかしながら、前述したとおり、一般化周期外乱オブザーバの過渡応答性はローパスフィルタＧF（ｓ）の応答性によって決定づけられるため、カットオフ周波数を低くしすぎると、トルク波形追従の速応性が悪化する。
【０１３３】
このトレードオフを緩和するために、（８）式のＴdnｃｏｓ２ｎθ＋Ｔqnｓｉｎ２ｎθ＋２Ｔdcｃｏｓｎθ，Ｔqnｃｏｓ２θ＋Ｔdnｓｉｎ２ｎθ＋２Ｔdcｓｉｎｎθにあるｎθ成分を予め除去しておくことが望ましい。ｓｉｎｎθ，ｃｏｓｎθ成分の係数に着目すると、どちらも２・Ｔdcであること分かる。つまり、入力される信号に含まれる直流成分Ｔdcがゼロであれば、ｎθ成分は発生しない。その結果、（８）式のＴdnｃｏｓ２ｎθ＋Ｔqnｓｉｎ２ｎθ＋２Ｔdcｃｏｓｎθ，Ｔqnｃｏｓ２θ＋Ｔdnｓｉｎ２ｎθ＋２Ｔdcｓｉｎｎθは２ｎθ成分のみとなる。
【０１３４】
２ｎθ成分は、ｎθ成分よりもカットオフが容易になるため、上述したトレードオフの問題が緩和され，ローパスフィルタＧF（ｓ）の設計も容易となる。
【０１３５】
本実施形態２では、（５）式のトルク検出値τdetの代わりにトルク偏差Δτを周波数成分抽出器２２に入力して周期外乱のｎ次周波数成分ベクトルＴnを求める。本実施形態２では、低域トルク制御器９によってトルク偏差Δτの直流成分がなくなり、周波数成分抽出器の入力信号をトルク偏差Δτとしているため、上述した副次的な効果も得られる。
【０１３６】
［実施形態３］
本実施形態３では、周波数成分抽出器２２のローパスフィルタの次数を１に限定した上で、実施形態２の構成をさらに簡素化する。図１６は、本実施形態３の構成図である。
【０１３７】
一般化周期外乱オブザーバの周波数成分抽出器２２に用いるローパスフィルタＧF（ｓ）を（９）式に示す１次のローパスフィルタに限定した場合、実施形態２の図１５は、図１６のように制御ブロック図を等価変換できる。
【０１３８】
【数１０】
2018007287.tif 000012
【０１３９】
図１６は、図１５の周波数成分抽出器２２を周波数成分変換器２４に変更し、ローパスフィルタＧF（ｓ）を積分器２３に変更している。すなわち、図１５の周波数成分抽出器２２の内部にあるローパスフィルタＧF（ｓ）を外に出す。そして、積分器２３は、周波数成分抽出器２２の内部にあるローパスフィルタＧF（ｓ）と後段のローパスフィルタＧF（ｓ）と合成した結果で得られ、カットオフ周波数ωｆと逆モデルをゲインとした単純な積分器２３となる。このように図１５は図１６に等価変換できる。 ただし、図１６の周波数成分変換器２４は、（１０）式を用いる。
【０１４０】
【数１１】
2018007287.tif 000013
【０１４１】
（１０） 式にはローパスフィルタＧF（ｓ）がないため、トルク偏差Δτに含まれる周期外乱の周波数成分を明示的に抽出することをせずに、直接、加振トルク制御器出力のｎ次周波数成分ベクトルＴpdn*を生成することになる。
【０１４２】
本実施形態３によれば、ローパスフィルタＧF（ｓ）を１次に限定した場合、トルク偏差に対する一般化周期外乱オブザーバは、ωｆをゲインとした積分器２３のみで構成される非常に簡単な構成で実現できる。したがって、実施形態２の効果に加えて、制御に係る演算量をさらに低減できる。
【０１４３】
［実施形態４］
本実施形態４では、加振トルク制御器１１内部の一般化周期外乱オブザーバにおいて、オブザーバゲインを追加した構成を説明する。図１７は、実施形態３の構成にオブザーバゲインＫｏｂを追加したものである。もちろん、他の実施形態の構成についても、逆モデルＱnに対してオブザーバゲインＫｏｂを追加する構成を用いれば同様の効果が得られる。
【０１４４】
オブザーバゲインＫｏｂは、逆モデルＱnに対してゲインを乗じるため、逆モデルの伝達特性の振幅に対して、敢えて誤差を与える形になる。通常は、制御対象モデルＰnの逆数をＯｎに設定して、モデル誤差がない真値に可能な限り一致した状態で動作させるが、一般化周期外乱オブザーバはモデル誤差に対する一定のロバスト安定性を有するため、そのロバスト安定範囲内でフィードバックループゲインを高めることを提案する。例として、一般化周期外乱オブザーバの安定性解析例を図１８に示す。
【０１４５】
図１８は、制御対象モデル（真値）Ｐnに対して、その逆モデルＱnが振幅誤差（縦軸）と位相誤差（横軸）を持つ場合の安定範囲を示した数値例である。この数値例では、位相誤差が０に近いほど振幅誤差を大きく許容し，位相誤差が±９０度を超えると不安定となることが分かる。
【０１４６】
また、図１８中の点線で示した部分は、離散系の極配置において最も原点に近づく（最も速応性が良い）場所を示している。すなわち、位相誤差は可能な限り０度としておき、オブザーバゲインＫｏｂを適度に与えることで、過渡変化時の速応性が向上することを意味する。
【０１４７】
本実施形態４では、このように一般化周期外乱オブザーバのロバスト安定範囲を用いて、フィードバックループのオブザーバゲインＫｏｂを適切に与えることで、加振波形追従制御の過渡速応性を向上することができる。
【０１４８】
［実施形態５］
実施形態１〜４の加振トルク制御器１１において、一般化周期外乱オブザーバ方式は、特定周波数成分のみに寄与する周期外乱抑制制御器ではあるが、本実施形態５では、抑制したい次数を別個指定し、それらの一般化周期外乱オブザーバを並列化する。補償信号合成部１９において、並列化した周期外乱補償値を合算することで、複数の周波数成分の周期外乱を同時に抑制することができる。
【０１４９】
トルクリプルはモータの電磁気的な不均一性，機械的アンバランス，インバータのデッドタイム，電流センサ誤差，そのほか非線形性などに起因して発生する周期外乱であり、例えば、抑制対象の次数として、電気的周波数で１次，２次，６次，１２次などを指定する。モータ極数や構造上の特性により、トルクリプルの発生しやすい次数を指定すれば良い。
【０１５０】
例として、図１３の制御構成を２並列する構成を説明する。並列１段目の制御構成のトルク指令値τ*とトルク検出値τdetと並列２段目の制御構成のトルク指令値τ*とトルク検出値τdetは、同じ値が加振周波数成分抽出器１２とリプル抑制周波数成分抽出器１３に入力される。
【０１５１】
また、並列１段目の制御構成の抑制対象次数ｎは６次，並列２段目の制御構成の抑制対象次数ｎは１２次とする。したがって、並列１段目の制御構成のｎ次回転位相ｎθと並列２段目の制御構成のｎ次回転位相ｎθは異なる値となる。
【０１５２】
さらに並列１段目の制御構成の加振トルク指令値τpd*と並列２段目の制御構成の加振トルク指令値τpd*を合算した値を、図５の加振トルク制御器１１の出力（加振トルク指令値τpd*）とする。
【０１５３】
また、加振制御ではエンジントルク脈動を模擬するが、市販されている自動車の多くに使用されている４サイクルエンジンでは、各気筒が２回転に１回燃料を爆発させるため、気筒数×０．５×回転数の周波数を持つ大きな振動トルクが発生する。また、その高調波成分も考慮し、４気筒エンジンでは機械的周波数の２次，４次，６次，８次などを指定してエンジン加振波形を模擬すれば良い。
【０１５４】
これら複数の制御対象次数に対する一般化周期外乱オブザーバを並列化することで、エンジン波形を模擬しながらトルクリプル抑制も同時に実現することができる。
【０１５５】
［実施形態６］
通常、エンジントルクの歪み波形はエンジン回転数に同期して発生する高調波成分であるから、これまでの実施形態１〜５で示したように回転位相検出値θ（＝∫ωmｄｔ）に同期した制御系で加振トルク制御器１１を構成していた。
【０１５６】
ダイナモメータのような試験装置では、通常は上述のように回転同期加振でテストを実施するが、特殊な事例として、モータ回転と非同期な加振制御を実施することも可能である。
【０１５７】
その際、トルクリプル抑制制御で実施する回転に同期した次数とは別に、回転に非同期で発生させる加振周波数成分が重畳されることになる。この場合、図１３の加振周波数成分抽出器１２にはｎ次回転位相ｎθが入力させず、ｎ次回転位相ｎθと非同期となる位相Θ’（＝∫ω’ｄｔ）が入力される。
【０１５８】
ここでω’は加振制御のための周波数成分に相当するモータ回転数ωmとは非同期で異なる値である。ω’は、位相Θ’を微分することによって得られる。位相Θまたはω’は、モータ回転数ωmとは関係なく別途設定する。
【０１５９】
このような加振周波数成分が非同期な条件において、実施形態５のような、図１３の制御構成が２並列している構成（ただし、並列１段目のｎの値（以下、ｎ１と記す）と並列２段目のｎの値（以下、ｎ２と記す）は異なる。並列１段目の加振制御の位相Θ’（以下、Θ１’と記す）と並列２段目の加振制御の位相Θ’（以下、Θ２’と記す）は異なる。）の場合に、モータ回転数ωmによっては並列１段目のトルクリプル抑制制御の周波数成分（ｎ１×ωm）と、並列２段目の加振の周波数成分（ω２’＝ｄΘ２’／ｄｔ）が一致する事例が発生する可能性がある。
【０１６０】
この場合、並列化した両方の制御系で同一の周波数成分を制御することになり、場合によっては制御干渉を引き起こすおそれがある。
【０１６１】
そこで、本実施形態６では、非同期加振のテストを実施する場合に限り、各並列段の制御構成での非同期加振周波数成分（ω１’、ω２’）とｎ次回転周波数すなわちトルクリプル周波数成分（ｎ１×ωm、ｎ２×ωm）を監視しておき、これらの周波数が一致した場合は、一致した段の制御構成のトルクリプル抑制制御をＯＦＦにする。
【０１６２】
例として、並列１段目の（ｎ１×ωm）と並列２段目のω２’が一致した場合では、並列１段目のトルクリプル抑制制御をＯＦＦにする。すなわち、図１３において、Ｕn＾＝０とする。
【０１６３】
本実施形態６によれば、非同期な加振制御を行う場合にトルクリプル抑制制御系との制御干渉を防止することができる。また、トルクリプル抑制制御をＯＦＦにしている期間は、加振トルク制御系がトルクリプル抑制の役割を包含するため、エンジン波形追従制御の結果には影響を与えない。
【０１６４】
［実施形態７］
実際のエンジンは、燃料の燃焼に失敗する失火モードが存在する。例えば、４気筒エンジンで１気筒のみ失火した場合、４回に１回、トルクを出力できない状態となる。
【０１６５】
本実施形態７では、エンジン失火モードに対応するため、小数次数を指定した一般化周期外乱オブザーバで加振トルク制御器１１を実装する。
【０１６６】
例えば、４気筒エンジンで１気筒の失火モードを実施する際は、機械回転数に対して０．５次およびその倍数を指定する。
【０１６７】
本実施形態７の効果の一例を図１９に示す。４回に１回、トルクが低下しており、そのような加振トルク指令値が与えられた場合でも追従できる。
【０１６８】
以上、     本発明において、記載された具体例に対してのみ詳細に説明したが、本発明の技術思想の範囲で多彩な変形および修正が可能であることは、当業者にとって明白なことであり、このような変形および修正が特許請求の範囲に属することは当然のことである。
【符号の説明】
【０１６９】
１…供試体
２…駆動モータ
３…吸収モータ
４…軸トルクメータ
５…加振制御器
６…回転位置センサ
７…速度制御器
９…低域トルク制御器
１０…高域共振抑制制御器
１１…加振トルク制御器
【図１】
2018007287.tif 000014
【図２】
2018007287.tif 000015
【図３】
2018007287.tif 000016
【図４】
2018007287.tif 000017
【図５】
2018007287.tif 000018
【図６】
2018007287.tif 000019
【図７】
2018007287.tif 000020
【図８】
2018007287.tif 000021
【図９】
2018007287.tif 000022
【図１０】
2018007287.tif 000023
【図１１】
2018007287.tif 000024
【図１２】
2018007287.tif 000025
【図１３】
2018007287.tif 000026
【図１４】
2018007287.tif 000027
【図１５】
2018007287.tif 000028
【図１６】
2018007287.tif 000029
【図１７】
2018007287.tif 000030
【図１８】
2018007287.tif 000031
【図１９】
2018007287.tif 000032