(57)【特許請求の範囲】
【請求項１】
一対の電磁石と、前記一対の電磁石の間に配置され回転する回転体とを有し、消費電力削減を目的としたゼロバイアス線形制御を行う磁気軸受装置であって、
前記ゼロバイアス線形制御は、一方の電磁石に流す電流をｉ１、他方の電磁石に流す電流をｉ３、ノーマルエアギャップをｘ０、電磁石の変位をｘｕ、電磁石が発生する磁気吸引力の合力をＦｘｕ、定数をＫｕとする下記式に従い前記一対の電磁石のうち一方に電流を流し、他方の電磁石への電流をゼロにして前記回転体の回転制御を行うものであり、
ジャイロ効果の影響による不安定を抑制するために、運動方程式に回転周波数に依存する項を付加した線形パラメータ変動システムを含んでおり、凸補間した制御器を予め複数用意しておき、回転周波数によって前記制御器を切り替えるゲインスケジュール制御を行うことを特徴とする磁気軸受装置。
【数１】
0005445887.tif 000147
【請求項２】
一対の電磁石と、前記一対の電磁石の間に配置され回転する回転体とを有し、消費電力削減を目的としたゼロバイアス線形制御を行う磁気軸受装置であって、
前記ゼロバイアス線形制御は、一方の電磁石に流す電流をｉ１、他方の電磁石に流す電流をｉ３、ノーマルエアギャップをｘ０、電磁石の変位をｘｕ、電磁石が発生する磁気吸引力の合力をＦｘｕ、定数をＫｕとする下記式に従い前記一対の電磁石のうち一方に電流を流し、他方の電磁石への電流をゼロにして前記回転体の回転制御を行うものであり、
ジャイロ効果の影響による不安定を抑制するために、運動方程式に回転周波数に依存する項を付加した線形パラメータ変動システムを含んでおり、変動パラメータによって常に変動する超平面に状態量を追従させるジャイロ保障のためのゲインスケジュール制御を行うことを特徴とする磁気軸受装置。
【数２】
0005445887.tif 000148
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、フライホイール式電力貯蔵装置などに使用され、回転体を磁気軸受で非接触支持して電動モータにより回転させる磁気軸受装置に関する。
【背景技術】
【０００２】
この種の磁気軸受装置として、機械本体とコントローラがケーブルで接続されており、機械本体に、回転体と、電磁石の磁気吸引力により回転体をアキシアル方向およびラジアル方向に非接触支持する制御型アキシアル磁気軸受およびラジアル磁気軸受と、回転体のアキシアル方向およびラジアル方向の変位を検出するための変位センサと、回転体を回転させる電動モータと、回転体の回転数を検出するための回転センサと、タッチダウン用の保護軸受とが設けられ、コントローラが、変位センサの出力に基づいて磁気軸受の電磁石を制御するとともに、回転センサの出力に基づいて電動モータを制御するものが知られている。
【発明の開示】
【発明が解決しようとする課題】
【０００３】
磁気軸受は強い非線形性を有しているため、線形理論を適用して制御を行うためには、平衡点近傍においてテイラー展開による一次近似を行うという線形化の手法がとられてきた。しかし、この線形化法は、回転体（ロータ）を常時浮上させるためには常にバイアス電流を供給しなければならない。さらに、バイアス電流を供給する手法では、回転時に回転体が磁束を常に横切ることとなり、相対的な回転体の磁束変動を発生させる。このため、磁性材料内部でのヒステリシス損失、渦電流損失が大きくなる。
【０００４】
また、フライホイール式電力貯蔵装置の場合など、回転体にフライホイールが取り付けられている磁気軸受装置は、回転周波数の変化とともに、ジャイロ効果による前向き振れ回り運動や後向き振れ回り運動が発生し、回転体の軌跡に花びらのような振れ回りが発生し、安定性を損なう。
【０００５】
ゼロバイアス制御磁気軸受の従来技術は、制御器の次数を低次元化するため、ジャイロ効果を考慮していないモデリングであった。そのため、実際に回転させると、速度上昇に伴ってジャイロ効果が大きくなり、システムの安定性が悪くなる。
【０００６】
この発明の目的は、上記の問題を解決し、消費電力削減を目的としたバイアス電流を用いない線形化を行い、ジャイロ効果の影響による安定性の劣化を抑制できる磁気軸受装置を提供することにある。
【０００７】
この発明の目的は、上記の問題を解決し、モデルベースからジャイロ効果を考慮し、回転速度によりコントローラのジャイロ補償のゲインをスケジューリングして、回転体の安定性を向上させた磁気軸受装置を提供することにある。
【課題を解決するための手段】
【０００８】
この発明による磁気軸受装置は、消費電力削減を目的としたゼロバイアス線形制御を行う磁気軸受装置において、ジャイロ効果の影響による不安定を抑制するために、運動方程式に回転周波数に依存する項を付加した線形パラメータ変動システムを含んでおり、凸補間した制御器を予め複数用意しておき、回転周波数によって前記制御器を切り替えるゲインスケジュール制御を行うことを特徴とするものである。
【０００９】
磁気軸受が非線形性を有するのは、吸引力と電磁石への供給電流との間の関係式であり、吸引力を制御入力と考えることで、磁気軸受は線形な制御対象とみなすことができる。つまり、回転体が１つの電磁石に接近するとき、接近した電磁石の制御電流を零にして、反対側の電磁石にのみ制御電流を供給し、回転体が完全に平衡点になれば、一切の電流を流さないということである。
【００１０】
さらに、モデルを回転周波数によりシステムが変動する線形パラメータ依存モデルであると考え、ゲインスケジュール制御器を設計する。実験システムへゲインスケジュール制御器の実装を考えたとき、サンプリングごとの離散化は計算処理が遅れるという問題に対しては、凸補間された制御器に対して、予め回転数ごとに制御器を用意しておき、スケジューリングパラメータによって制御器を切り替えるという方法を提案する。
【００１１】
吸引力を制御入力と考えることで磁気軸受は線形な制御対象とする線形化方法を用いて対向配置する電磁石には一方しか電流を流さない。
【００１２】
運動方程式に回転周波数に依存する項を付加し、線形パラメータ変動システムを構築することで、ジャイロ効果による影響をモデルに加える。
【００１３】
実験にゲインスケジュール制御を適用する際は、予め凸補間した制御器をいくつか用意しておき、回転周波数によって制御器を切り替える方法を適用する。
【００１４】
電流切換により、ゼロバイアス制御を実現した。さらに、フライホイール・磁気軸受系に対してテーブルを用いたゲインスケジュール制御を提案したことで、浮上時だけではなく、回転実験においても、ロータは剛性モードおよび弾性モードの後向き振れ回り運動による花びらのような軌跡を十分に制御によって抑えることができる。
【００１５】
この発明による磁気軸受装置は、消費電力削減を目的としたゼロバイアス線形制御を行う磁気軸受装置において、ジャイロ効果の影響による不安定を抑制するために、運動方程式にジャイロ効果の項を付加した線形パラメータ変動システムを含んでおり、変動パラメータによって常に変動する超平面に状態量を追従させるジャイロ補償のためのゲインスケジュール制御を行うことを特徴とするものである。
【００１６】
運動方程式にジャイロ効果の項を付加し、線形パラメータ変動システムを構築した。制御するため、状態量を変動パラメータによって常に変動する超平面に追従させるジャイロ補償のゲインスケジュール型スライディングモード制御則を考案した。
【００１７】
電流切換により、ゼロバイアス制御を実現した。さらに、根軌跡から見ると、回転速度に係わらずほとんど動かない状態で安定性を保持できた。
【発明の効果】
【００１８】
この発明の磁気軸受装置によれば、ジャイロ効果の影響による安定性の劣化を防止することができる。
【００１９】
この発明の磁気軸受装置によれば、回転体の安定性を向上させることができる。
【発明を実施するための最良の形態】
【００２０】
以下、図面を参照して、この発明をフライホイール式電力貯蔵装置に使用される５軸制御型磁気軸受装置に適用した実施形態について説明する。
【００２１】
図１は磁気軸受装置の全体構成を概略的に示すブロック図、図２は磁気軸受装置の機械的部分の主要部を示す縦断面図、図３はさらに図２の主要部を示す斜視図である。
【００２２】
磁気軸受装置は、ケーブルにより接続された機械本体(1)およびコントローラ（制御手段）(2)を備えている。
【００２３】
磁気軸受装置は、鉛直円筒状のケーシング(3)の内側で鉛直軸状の回転体（ロータ）(4)が回転する縦型のものであり、ケーシング(3)より突出した回転体(4)の上端部にフライホイール(4a)が固定されている。以下の説明において、回転体(4)の鉛直な軸方向（アキシアル方向）の制御軸（アキシアル制御軸）をＺ軸、Ｚ軸と直交するとともに互いに直交する２つの水平な径方向（ラジアル方向）の制御軸（ラジアル制御軸）をＸ軸およびＹ軸とする。
【００２４】
機械本体(1)には、回転体(4)を軸方向に非接触支持する１組の制御型アキシアル磁気軸受(5)、回転体(4)を径方向に非接触支持する上下２組の制御型ラジアル磁気軸受(6)(7)、回転体(4)の軸方向および径方向の変位を検出するための変位検出部(8)、回転体(4)を高速回転させるためのビルトイン型電動モータ(9)、回転体(4)の回転数を検出するための回転センサ(10)、ならびに回転体(4)の軸方向および径方向の可動範囲を規制して回転体(4)を磁気軸受(5)(6)(7)で支持していないときに回転体(4)を機械的に支持する上下２組のタッチダウン用の保護軸受(11)(12)が設けられている。
【００２５】
コントローラ(2)には、センサ回路(13)(14)、電磁石駆動回路(15)、インバータ(16)およびＤＳＰボード(17)が設けられ、ＤＳＰボード(17)には、ソフトウェアプログラムが可能なディジタル処理手段としてのＤＳＰ(18)、ＲＯＭ(19)、不揮発性記憶装置としてのＲＡＭ(20)、ＡＤ変換器(21)(22)およびＤＡ変換器(23)(24)が設けられている。
【００２６】
変位検出部(8)は、回転体(4)の軸方向の変位を検出するための１個のアキシアル変位センサ(25)、および回転体(4)の径方向の変位を検出するための上下２組のラジアル変位センサユニット(26)(27)を備えている。
【００２７】
アキシアル磁気軸受(5)は、回転体(4)の下部に一体に形成されたフランジ部(4b)をＺ軸方向の両側から挟むように配置された１対のアキシアル電磁石(28a)(28b)を備えている。アキシアル電磁石は、符号(28)で総称する。
【００２８】
アキシアル変位センサ(25)は、回転体(4)の下端面にＺ軸方向の下側から対向するように配置され、回転体(4)の下端面との距離（空隙）に比例する距離信号を出力する。
【００２９】
２組のラジアル磁気軸受(6)(7)は、アキシアル磁気軸受(5)の上側において上下方向に所定の距離をおいて配置されており、これらの間にモータ(9)が配置されている。上側のラジアル磁気軸受(6)は、回転体(4)をＸ軸方向の両側から挟むように配置された１対のラジアル電磁石(29a)(29b)、および回転体(4)をＹ軸方向の両側から挟むように配置された１対のラジアル電磁石(29c)(29d)を備えている。これらのラジアル電磁石は、符号(29)で総称する。同様に、下側のラジアル磁気軸受(7)も、２対のラジアル電磁石(30a)(30b)(30c)(30d)を備えている。これらの電磁石も、符号(30)で総称する。
【００３０】
上側のラジアル変位センサユニット(26)は、上側のラジアル磁気軸受(6)の近傍に配置されており、Ｘ軸方向の電磁石(29a)(29b)の近傍においてＸ軸方向の両側から回転体(4)を挟むように配置された１対のラジアル変位センサ(31a)(31b)、およびＹ軸方向の電磁石(29c)(29d)の近傍においてＹ軸方向の両側から回転体(4)を挟むように配置された１対のラジアル変位センサ(31c)(31d)を備えている。これらのラジアル変位センサは、符号(31)で総称する。同様に、下側のラジアル変位センサユニット(27)は、下側のラジアル磁気軸受(7)の近傍に配置されており、２対のラジアル変位センサ(32a)(32b)(32c)(32d)を備えている。これらのラジアル変位センサも、符号(32)で総称する。各ラジアル変位センサ(31)(32)は、回転体(4)の外周面との距離に比例する距離信号を出力する。
【００３１】
電磁石(28)(29)(30)、変位センサ(25)(31)(32)は、ケーシング(3)に固定されている。
【００３２】
保護軸受(11)(12)はアンギュラ玉軸受などの転がり軸受よりなり、各保護軸受(11)(12)の外輪がケーシング(3)に固定され、内輪が回転体(4)の周囲に所定の隙間をあけて配置されている。２組の保護軸受(11)(12)はいずれも径方向の支持が可能なものであり、少なくとも１組は軸方向の支持も可能なものである。この例では、上側の保護軸受(11)は径方向の支持のみを行い、下側の保護軸受(12)は径方向の支持と軸方向の支持を行うようになっている。
【００３３】
コントローラ(2)のＲＯＭ(19)には、ＤＳＰ(18)における処理プログラムなどが格納されている。ＲＡＭ(20)には、磁気軸受の制御パラメータなどが記憶されている。
【００３４】
センサ回路(13)は、変位検出部(8)の各変位センサ(25)(31)(32)を駆動し、各変位センサ(25)(31)(32)の出力に基づいて、回転体(4)の軸方向の変位、ならびに上下のラジアル変位センサユニット(26)(27)の部分におけるＸ軸方向およびＹ軸方向の変位を演算し、その演算結果である変位信号をＡＤ変換器(21)を介してＤＳＰ(18)に出力する。変位検出部(8)およびセンサ回路(13)により変位検出装置が構成され、変位検出部(8)のラジアル変位センサユニット(26)(27)およびセンサ回路(13)のラジアル変位センサ(31)(32)に関する部分によりラジアル変位検出装置が構成されている。
【００３５】
センサ回路(14)は、回転センサ(10)を駆動し、回転センサ(10)の出力を回転体(4)の回転数に対応する回転数信号に変換し、これをＡＤ変換器(22)を介してＤＳＰ(18)に出力する。
【００３６】
ＤＳＰ(18)は、ＡＤ変換器(21)から入力する変位信号に基づいて、各磁気軸受(5)(6)(7)の各電磁石(28)(29)(30)に対する制御電流値を求め、制御電流値に対応する励磁電流信号をＤＡ変換器(21)を介して磁気軸受駆動回路(15)に出力する。そして、駆動回路(15)は、ＤＳＰ(18)からの励磁電流信号に基づく励磁電流を対応する磁気軸受(5)(6)(7)の電磁石(28)(29)(30)に供給し、これにより、回転体(4)が所定の目標位置に非接触される。ＤＳＰ(18)は、また、ＡＤ変換器(22)から入力する回転センサ(10)からの回転数信号に基づいて、モータ(9)に対する電流制御信号をＤＡ変換器(24)を介してインバータ(16)に出力し、インバータ(16)は、この信号に基づいて、モータ(9)の電流を制御することにより、回転数を制御する。そして、その結果、回転体(4)が、磁気軸受(5)(6)(7)により目標位置に非接触支持された状態で、モータ(9)により高速回転させられる。
【００３７】
次に、ＤＳＰ(18)による上下のラジアル電磁石(29)(30)の制御について説明する。
【００３８】
フライホイール・磁気軸受系のゲインスケジュールＨ∞制御
磁気軸受は不安定な系であることから、フィードバック制御による安定化が課題である。さらに、磁気軸受は、強い非線形性を有しているため、線形理論を適用して制御を行うためには、平衡点近傍においてテイラー展開による一次近似を行うという線形化の手法がとられてきた。しかし、この線形化法は、ロータを常時浮上させるためには常にバイアス電流を供給しなければならない。さらに、バイアス電流を供給する手法では、回転時にロータが磁束を常に横切ることとなり、相対的なロータの磁束変動を発生させる。このため、磁性材料内部でのヒステリシス損失、渦電流損失が大きくなる。このため、エネルギ貯蔵用フライホイールなどにおけるエネルギ収支の観点から消費電力、回転損失が大きいという欠点を持つ。
【００３９】
そこで、近年、磁気軸受において消費電力の最小化を目指し、バイアス電流を用いない制御方式が盛んに研究されている。ある研究では、従来の電流を制御入力と考え、ゼロパワー制御を行っているが、制御電流に定常偏差が残る結果となった。従来、磁気軸受において、制御入力は電流として考えてきた。しかし、磁気軸受が非線形性を有するのは、吸引力と電磁石への供給電流との間の関係式であり、吸引力を制御入力と考えることで、磁気軸受は線形な制御対象とみなすことができる。
【００４０】
先行研究では、ロバスト制御理論の適用が行われ、その有効性が確認されている。また、非線形な制御型磁気軸受の特性に対して、吸引力を制御入力と考えるゼロバイアス手法を提案し、消費電力の低減化に成功した。
【００４１】
しかし、実験において、ロータの軌跡に乱れが発生し、１００Ｈｚ付近にてタッチダウンに至る。これは、慣性力を増加させることを目的にフライホイールがロータに取り付けられているために、ジャイロ効果による前向き触れ回り運動や後向き触れ回り運動の発生、およびそれらの固有振動数の変化が高速回転を妨げていることが原因であると考えられる。
【００４２】
そこで、本論文では、まず、ジャイロ効果の特性を考慮したモデリングを行った。プラントのＡ行列に回転周波数ωを定数として考慮し、Ｈ∞固定制御器の性能をシミュレーションおよび実験にて検証する。
【００４３】
次に、回転周波数によりシステムが変動する線形パラメータ依存モデル（ＬＰＶ）であると考え、ゲインスケジュール（ＧＳ）制御器を設計し、Ｈ∞固定制御器との比較をシミュレーションおよび固有値解析にて行う。実験システムへＧＳ制御器の実装を考えたとき、サンプリングごとの離散化は計算処理が遅れるという問題に対しては、凸補間された制御器に対して、予め回転数ごとに制御器を用意しておき、ＧＳパラメータによって制御器を切り替えるという方法を提案する。また、テーブルを用いることでＧＳ制御器の閉ループ安定性が損なわれないかを、初期応答により検証する。
【００４４】
最後に、実験において制御性能を劣化させていると考えられる曲げ一次弾性バックモードを含む低次元化モデルを用いて検証する。シミュレーションを実験の不一致問題に対して、原因がアクチュエータの諸元によるものと考え、ＧＳ制御器にαという未知パラメータ係数を加えたときの、閉ループの固有値解析を行う。
【００４５】
本論文は、０．５ＫＷｈ級磁気軸受フライホイールシステムに対して、ＬＭＩに基づくＧＳ制御の有効性について論じたものである。
【００４６】
制御対象は、図１〜図３に示す磁気軸受装置である。
【００４７】
３．制御理論
３．１  概念
Ｈ∞制御では、種々の制御問題を統一的な枠組みで扱えるよう、図４に示すフィードバック系が用いられる。図４において、Ｇはプラントと呼ばれ、次式で示す入出力信号を持つ伝達関数として定義される。
【数１】
0005445887.tif 000002
【００４８】
また、Ｇの状態空間実現は、次のように定義される。
【数２】
0005445887.tif 000003
【００４９】
ここで、ｗ∈Ｒｍ１は外部入力と呼ばれ、参照信号や外乱、センサノイズなど、制御系に外部から加わる入力を表し、ｚ∈Ｒｐ１は制御量と呼ばれ、制御偏差や制御入力など、制御によって小さくしたい量を表す。また、ｕ∈Ｒｍ２およびｙ∈Ｒｐ２は、制御入力と観測出力で、それぞれ、補償器からの出力および入力となる量である。今、一般化プラントＧに対して補償器
ｕ＝Ｋｙ    (3.3)
を用いてフィードバック制御を行うと、ｗからｚまでの伝達関数は次式となる。
【００５０】
ｚ＝Ｇｚｗｗ
ただし、
Ｇｚｗ：＝Ｇ１１＋Ｇ１２Ｋ(Ａ−Ｇ２２Ｋ)−１Ｇ２１    (3.4)
である。制御目的は外部入力ｗに対して、制御量ｚをなるべく小さく抑えることである。したがって、伝達関数Ｇｚｗの大きさを何らかの意味で小さくする補償器Ｋを設計すればよい。Ｈ∞制御においては、このＧｚｗの大きさの尺度としてＨ∞ノルムというものを用いる。Ｇｚｗが安定であるとき、そのＨ∞ノルムを‖Ｇｚｗ‖∞と書き、次のように定義する。
【数３】
0005445887.tif 000004
【００５１】
ここで、σは行列の最大特異値を表し、次式で定義される。
【数４】
0005445887.tif 000005
【００５２】
ここで、λｍａｘは行列の最大固有値、
【数５】
0005445887.tif 000006
【００５３】
は行列Ｇｚｗ(ｊｗ)の複素共役転置である。このとき、Ｈ∞ノルムは入出力のエネルギ比の最大値を意味する。
【００５４】
以上のことを踏まえて、Ｈ∞制御問題は、図４の一般化プラントＧに対してｕ＝Ｋｙなるフィードバック制御を行うことにより、閉ループを安定化し、かつ、ある与えられた整数γに対して
‖Ｇｚｗ(ｓ)‖∞＜γ    (3.7)
を満たす補償器Ｋを求めることである。このようにＨ∞制御問題は非常にシンプルなものとして定式化されており、Ｇの構成を変えることで、つまり外部入力ｗ、制御量ｚを制御目的に合わせて選ぶことにより、ロバスト制御、外乱抑圧制御、感度最適化、目標値追従、そして閉ループ整形問題など様々な問題に適用することができる。また、式(3.2)の定義、さらに以下の仮定の下ではＨ∞制御問題は、標準Ｈ∞制御問題と呼ばれる。
【００５５】
仮定１  (Ａ,Ｂ２)は可安定、かつ(Ｃ２,Ａ)は可検出
仮定２  Ｄ１２は縦長列フルランク、かつＤ２１は横長行フルランク
仮定３  Ｇ１２は虚軸上に不変零点を持たない、すなわち、すべてのｗに対して、
【数６】
0005445887.tif 000007
【００５６】
はフルランク
仮定４  Ｇ２１は虚軸上に不変零点を持たない、すなわち、すべてのｗに対して、
【数７】
0005445887.tif 000008
【００５７】
はフルランク
３．２  外乱抑圧制御とＨ∞制御
フィードバック制御は制御対象を安定化させることだけが目的ではなく、抑制問題、たとえば外乱抑制特性、目標値追従特性を達成することも目的である。前節のＨ∞制御問題においてｗを実外乱、ｚを外乱の影響を抑えたい物理量にとれば外乱抑制制御そのものになる。今、図５(a)のように、制御対象の入力端に加わる外乱ｗを出力端ｚで抑圧する問題を考える。ｚ、ｕはそれぞれ
ｚ＝ｙ＝Ｐ(ｕ＋ｗ)  ,  ｕ＝Ｋｙ    (3.8)
と表せる。これより、ｗからｚまでの伝達関数は
ｚ＝(Ａ−ＰＫ)−１Ｐｗ    (3.9)
と求まる。つまり、図５(a)における外乱抑圧問題は、なるべく小さな正数γに対して、
‖(Ａ−ＰＫ)−１Ｐ‖∞＜γ    (3.10)
を満たす補償器Ｋを求める問題となる。ただし、式(3.10)の評価では、Ｈ∞ノルムの定義よりｗからｚまでの伝達関数のゲイン（多入出力系ならば最大特異値）をすべての周波数帯域にわたってγ以下に押さえることを意味し、一般には厳しい要求となる。しかし、実際には、すべての周波数帯域で周波数スペクトルが大きいという外乱は見られない。そこで、外乱ｗの周波数スペクトルが大きな周波数帯域でのみ(Ａ−ＰＫ)−１Ｐを小さくできれば十分となる。そこで、小さくしたい周波数帯域でのみ十分大きなゲインを持つ伝達関数ＷＭを導入し、式(3.10)の代わりに、
‖ＷＭ(Ａ−ＰＫ)−１Ｐ‖∞＜γ    (3.11)
とする。今、
Ｍｒ：＝(Ａ−ＰＫ)−１Ｐ    (3.12)
と定義し、Ｍｒが１入出力の場合を考えると、Ｈ∞ノルムの定義より、式(3.11)は
｜ＷＭ(ｊｗ)｜｜Ｍｒ(ｊｗ)｜＜γ,∀ｗ    (3.13)
に等価であるので、｜ＷＭ(ｊｗ)｜が大きな周波数帯域では、｜Ｍｒ(ｊｗ)｜をとくに小さくすることができることがわかる。
【００５８】
３．３  目標値追従制御とＨ∞制御
次に、ｗを目標値、ｚを出力信号ととれば、目標値追従制御そのものになる。今、図６において、ｗからｚまでの伝達関数は
ｚ＝(Ａ−ＰＫ)−１ｗ    (3.14)
と求まる。つまり、図６(a)における目標値追従問題は、なるべく小さな正数γに対して、
‖(Ａ−ＰＫ)−１‖∞＜γ    (3.15)
を満たす補償器Ｋを求める問題となる。ただし、前節同様、式(3.15)の評価では、ｗからｚまでの伝達関数のゲイン（多入出力系ならば最大特異値）をすべての周波数帯域にわたってγ以下に押さえるは、一般には厳しい要求となる。そこで、小さくしたい周波数帯域でのみ十分大きなゲインを持つ伝達関数ＷＭを導入し、式(3.15)の代わりに、
‖ＷＭ(Ａ−ＰＫ)−１‖∞＜γ    (3.16)
とする。今、
Ｍｔ：＝(Ａ−ＰＫ)−１    (3.17)
と定義し、Ｍｔが１入出力の場合を考えると、式(3.16)は
｜ＷＭ(ｊｗ)｜｜Ｍｔ(ｊｗ)｜＜γ,∀ｗ    (3.18)
に等価であるので、｜ＷＭ(ｊｗ)｜が大きな周波数帯域では、｜Ｍｔ(ｊｗ)｜をとくに小さくすることができることがわかる。また、一般に式(3.17)は感度関数と呼ばれる。
【００５９】
３．４  ロバスト安定化問題
Ｈ∞制御ではロバスト安定化問題を陽に取り扱うことができるが、これはロバスト安定性を保証する際に用いるスモールゲイン定理がＨ∞ノルム条件で記述されているからである。
【００６０】
定理３．１（スモールゲイン定理）
図７において、Ａ(ｓ)およびＢ(ｓ)は安定でプロパな伝達関数とする。このとき、
‖Ａ(ｓ)Ｂ(ｓ)‖∞＜１    (3.19)
を満たすと、図７の閉ループは安定となる。
【００６１】
３．５  乗法的誤差に対するロバスト安定化
図８(a)において、フィードバック系のロバスト安定化問題を考える。実制御対象Ｐは公称モデルＰとそれに対する乗法的誤差Δを用いて
Ｐ＝(Ｉ＋Δ)Ｐ    (3.20)
と表される。ここで、Δは簡単のため安定と仮定する。また、図８(a)において、点ａから点ｂまでの伝達関数を求めると
Ｔｍ：＝(Ｉ−ＰＫ)−１ＰＫ    (3.21)
となる。ここで、式(3.21)を相補感度関数と呼ぶ。このとき、図８(a)と図８(b)は等価であるので、スモールゲイン定理を適用すると、安定となるための条件
‖ΔＴｍ‖∞＜１    (3.22)
を得る。しかし、モデル化誤差を正確に見積もることは不可能である。したがって、
σ{Δ(ｊｗ)}＜｜ｗｍ(ｊｗ)｜  ,  ∀ｗ    (3.23)
を満たす既知のスカラ伝達関数ｗｍを用いることで
‖ｗｍＴｍ‖∞＜１    (3.24)
が成り立つことが、ロバスト安定のための条件となる。したがって、図８(c)の一般化プラントを構成し、ｗからｚまでのＨ∞ノルムを１未満とする補償器Ｋを求めればよい。
【００６２】
３．６  混合感度問題
本来、フィードバック制御を行う目的は、ある種の制御性能を達成しながらも、ロバスト安定化を図ることである。つまり、制御性能のみを考慮しても、必ず存在するモデル化誤差を保証しない。また、ロバスト安定化のみを図っても、それは十分条件でしかない。やはり、制御性能、ロバスト安定化の両方を考慮した設計をしなければならない。今、図９のブロック線図を考える。ここで、ｗ１、ｗ２は外部入力、ｚ１、ｚ２は制御量である。ｗ２からｚ１までの伝達関数Ｇｚ１ｗ２は
Ｇｚ１ｗ２＝Ｗｓ(ｓ)Ｍｔ(ｓ)    (3.25)
ここで、
Ｍｔ(ｓ)＝(Ｉ＋ＰＫ)−１    (3.26)
と表され、目標値追従性を達成するものであり、式(3.26)は重み付き感度関数でもある。これを小さくすることが感度低減問題となる。一般に、低周波域に重みを付けて感度を下げることが多い。したがって、感度低減問題は次式の評価関数となり、できるだけ小さいγを探索する問題となる。
【００６３】
‖Ｇｚ１ｗ２‖∞＜γ    (3.27)
また、外乱抑圧特性を達成するためには、ｗ１からｚ１までの伝達関数Ｇｚ１ｗ１
Ｇｚ１ｗ１＝Ｗｓ(ｓ)Ｍｒ(ｓ)    (3.28)
ここで、
Ｍｒ(ｓ)＝(Ｉ＋ＰＫ)−１Ｐ    (3.29)
においても、できるだけ小さいγを探索する問題となる。
【００６４】
‖Ｇｚ１ｗ１‖∞＜γ    (3.30)
一方、ロバスト安定性を確保するためには、ｗ１からｚ２までの伝達関数Ｇｚ２ｗ１は
Ｇｚ２ｗ１＝Ｗｔ(ｓ)Ｔｍ(ｓ)    (3.31)
ここで、
Ｔｍ(ｓ)＝(Ｉ＋ＰＫ)−１Ｋ    (3.32)
であるので、
‖Ｇｚ２ｗ１‖∞＜１    (3.33)
を満たすように設計しなければならない。さらに、ｗ２からｚ２までの伝達関数Ｇｚ２ｗ２
Ｇｚ２ｗ２＝Ｗｔ(ｓ)Ｔａ(ｓ)    (3.34)
ここで、
Ｔａ(ｓ)＝(Ｉ＋ＰＫ)−１Ｋ    (3.35)
においても
‖Ｇｚ２ｗ２‖∞＜１    (3.36)
を満たすように設計しなければならない。なお、式(3.35)を一般的に準相補感度関数と呼ぶ。式(3.27)、式(3.30)、式(3.33)、式(3.36)を同時に満たすのは相反する要求であり、トレードオフを考えなければならない。現実には、Ｓの重みＷｓ(ｊｗ)は低周波域で大きくとればよく、制御対象の不確定性を反映するＷｔ(ｊｗ)は高周波域で大きくとればよい。よって、低周波域でＭｔ(ｓ)、Ｍｒ(ｓ)を、高周波域でＴｍ(ｓ)、Ｔａ(ｓ)を小さくするという設計指針が導かれる。このような問題を混合感度問題という。式(3.27)、式(3.30)、式(3.33)、式(3.36)の組み合わせでは
【数８】
0005445887.tif 000009
【００６５】
となり、式(3.37)を同時に満たす補償器Ｋ(ｓ)を設計する問題となる。今、ｗ＝[ｗ１ ｗ２]Ｔ、[ｚ１ ｚ２]Ｔなる変数を定義すると、ｗからｚまでの伝達関数は
【数９】
0005445887.tif 000010
【００６６】
となり、式(3.37)の代わりに
【数１０】
0005445887.tif 000011
【００６７】
を満たす補償器Ｋ(ｓ)を求める問題となる。
【００６８】
３．７  ＬＭＩに基づくゲインスケジュールＨ∞制御
ＬＭＩベースのＧＳＨ∞制御の定式化は文献に詳しく述べられている。本研究ではこの定式化に基づき制御系を設計することである。制御対象プラントと制御器の状態空間表現を次式で表す。
【数１１】
0005445887.tif 000012
【数１２】
0005445887.tif 000013
【００６９】
ここで、ｘおよびｘｋはそれぞれプラントと制御器の状態ベクトル、ｚとｙは制御量と測定された出力ベクトルを表す。ｕは制御入力、ｗは外乱入力ベクトルである。2つの式を組み合わせると、次式の閉ループ系の式が得られる。
【数１３】
0005445887.tif 000014
【００７０】
閉ループ行列Ａｃｌ、Ｂｃｌ、Ｃｃｌ、Ｄｃｌは
【数１４】
0005445887.tif 000015
【００７１】
ここで
【数１５】
0005445887.tif 000016
【００７２】
制御器行列は単位行列Ωにまとめられることに注目しよう。リアプノフ関数Ｖ(ｘ)＝ｘＴＰｘ、Ｐ＞０が式(3.42)の閉ループ系に対して定義されれば、大域的漸近安定性を保証できる。今、ＬＴＩシステムのｗからｚまでのＬ２誘導ノルムは
‖ｚ‖２＜γ‖ｗ‖２    (3.45)
と与えられる。これから最終的に、次式を満たす正定なリアプノフ関数Ｖ(ｘ)＝ｘＴＰｘが存在する。
【数１６】
0005445887.tif 000017
【００７３】
Ｈ∞準最適制御問題は、Ｘｃｌ＞０に対して次の不等式の解が存在することと等価である。
【数１７】
0005445887.tif 000018
【００７４】
さらに、ＬＭＩ式(3.47)の解は次の不等式を満足する２つの対称行列ＲおよびＳが存在することと等価である。
【数１８】
0005445887.tif 000019
【数１９】
0005445887.tif 000020
【数２０】
0005445887.tif 000021
【００７５】
ここで、ＮＲおよびＮＳはそれぞれ
【数２１】
0005445887.tif 000022
【００７６】
および(Ｃ２,Ｄ１２)の零空間を表す。
【００７７】
上記のＨ∞制御問題はＬＴＩシステムにのみ有効であるが、ＬＰＶシステムにまで拡張することも可能である。ＬＰＶプラントを状態空間モデルで表すと
【数２２】
0005445887.tif 000023
【００７８】
ここで、ｘ、ｙ、ｕはそれぞれ状態ベクトル、観測出力ベクトル、制御入力ベクトルを表す。ｐは時変プラントのパラメータベクトルであり、この結果、制御対象のプラントはｐの関数となる。実際にｐは速度や減衰、剛性といった時間依存の物理パラメータとなっており、次のように変動幅の最小値と最大値で与えられる。
【００７９】
ｐｍｉｎ≦ｐｉ(ｔ)≦ｐｍａｘ    (3.52)
ｐ(ｔ)が制御中に大きな変動を受けるとき、固定されたロバストＬＴＩ制御器で高性能な制御性能を得ることはできない。
【００８０】
もし、ｐ(ｔ)の実測値が制御中に実時間で得られるならば、制御器の動特性を実測値ｐ(ｔ)に依存して実時間で変化させることができる。このとき、制御器は次式で示される。
【数２３】
0005445887.tif 000024
【００８１】
ここで、ｙは観測ベクトル、ｕは制御入力である。時変パラメータを正確に実測できれば、制御器をプラントダイナミクスの変動に追従させて連続的に調整することにより大域的安定性と良好な制御性能を実現できる。
【００８２】
現在のパラメータ値ｐ(ｔ)の凸分解
【数２４】
0005445887.tif 000025
【００８３】
を与えると、Ａｋ(ｐ)、Ｂｋ(ｐ)、Ｃｋ(ｐ)、Ｄｋ(ｐ)の値は、
【数２５】
0005445887.tif 000026
【００８４】
によって、パラメータボックスの端点での、
【数２６】
0005445887.tif 000027
【００８５】
の値から導出される。言い換えれば、動作点ｐ(ｔ)でのコントローラの状態空間行列は、ＬＴＩ端点コントローラ
【数２７】
0005445887.tif 000028
【００８６】
の凸補間で得られる。これにより、時変パラメータの実測値ｐ(ｔ)によるコントローラ行列のスムーズなスケジューリングが行える。
【００８７】
今、ＬＰＶプラントの状態空間表現を一般的な表現で考える。
【数２８】
0005445887.tif 000029
【００８８】
行列Ａ(・)、Ｂ(・)、Ｃ(・)、Ｄ(・)は変動パラメータｐの関数である。また、Ｂ２、Ｃ２、Ｄ１２、Ｄ２１行列はｐの影響を受けない。
【００８９】
コントローラに対して、図１０の相互接続に関するシンセシス問題を考える。ＬＰＶプラントに対するＨ∞制御問題の解は、次式のようにＬＴＩプラントと同じ形をしている。
【数２９】
0005445887.tif 000030
【数３０】
0005445887.tif 000031
【数３１】
0005445887.tif 000032
【００９０】
ここで、Ａｉ、Ｂ１ｉ、Ｃ１ｉ、Ｄ１１ｉはポリトープ形パラメータｐ＝ＰｉでのＡ(ｐ)、Ｂ１(ｐ)、Ｃ１(ｐ)、Ｄ１１(ｐ)のパラメータ値である。不等式(3.57)〜(3.59)は凸最適化アルゴリズムを使ったツールボックスで解くことができる。行列ＲとＳから得られた制御器行列Ωの構造も同じ凸プログラムで得られる。
【００９１】
４．モデリング
４．１  ジャイロモーメント
ジャイロモーメントとは、高速回転円板が傾き運動するときに発生する作用であり、円板の傾きと直角方向に発生するモーメントのことである。ジャイロモーメントは以下の式により求められる。
【数３２】
0005445887.tif 000033
【００９２】
この効果により、静止時には１つである固有振動数が前向き振れ回りモードと後向き振れ回りモードに分かれ、ジャイロ効果の大小によって変化幅が決まるという性質がある。
【００９３】
ジャイロ効果により、ロータは磁気軸受にとって回転数に伴って振動特性が変化するとともにＸ−Ｙ方向が連成する制御対象になるので、磁気軸受の制御特性を決める上で大きな障害となる。したがって、回転周波数に対してロータの振動特性がどのように変化するかを、予め精度良く求めておくことが必要である。
【００９４】
極軌跡の虚部、すなわち固有振動数の変化を縦軸に、回転周波数を横軸にとった固有振動数線図を図１１に示す。式(4.1)から回転数と傾き角速度が増加するほどジャイロモーメントは増加することより、２次モードの方がジャイロモーメントが大きく、前章での実験におけるタッチダウンの原因であると考えられる。したがって、剛性２次モードと運搬線の交点での危険速度を越えるためには、ジャイロ効果を考慮した設計が必要である。
【００９５】
４．２  ジャイロ効果を考慮したモデリング
本研究が対象とする磁気軸受は、図１２のような５軸制御型フライホイールであり、そのパラメータを表４．１に示す。アキシアル方向はＰＩＤ制御がなされているので、ラジアル方向についてのみモデル化を考える。
【表１】
0005445887.tif 000034
【００９６】
図１２、図１３の図中の記号は、
Ｇ：ロータの重心位置
θｘ,θｙ：各軸周りの回転角度
Ｌｕ,Ｌｌ：上部、下部の電磁石までの距離
Ｌ１,Ｌ２：上部、下部のセンサまでの距離
ｉ１〜ｉ８：各電磁石における制御電流
ｘｇ,ｘｕ,ｘｌ：Ｘ−Ｚ平面におけるロータ重心、上部電磁石、下部電磁石位置での変位
を示す。ジャイロ効果
【数３３】
0005445887.tif 000035
【００９７】
を考慮したロータの重心に対する並進方向、回転方向の運動方程式は式(4.2)となる。
【数３４】
0005445887.tif 000036
【００９８】
出力ｙは
【数３５】
0005445887.tif 000037
【００９９】
であるので、状態変数、入力、出力を
【数３６】
0005445887.tif 000038
【０１００】
とすると、その状態空間実現は
【数３７】
0005445887.tif 000039
【数３８】
0005445887.tif 000040
【数３９】
0005445887.tif 000041
【０１０１】
である。
【０１０２】
５．制御系設計と性能検証１
５．１  Ｈ∞制御器設計
設計において、一般化プラントは図１４の修正混合感度問題を使用した。外乱抑圧特性を得るために仮想外乱をＰの直前に導入し、目標値追従特性を得るために仮想外乱をＰの直後に導入している。なお、εは仮想外乱を観測ノイズとみなすためのものでもあり、標準Ｈ∞制御問題の可解条件を満たすためのものである。重み関数は上側、下側の質量の割合が異なるため、別々に設計した。設計に用いた重み関数の周波数応答を図１４(a)に示す。また、ランク条件を満たすために導入したε値を０．０４とした。γ値は０．８９である。
【０１０３】
Ｗｔはロバスト性に対する周波数重み関数であり、高周波数領域での制御入力の影響を低減する。また、Ｗｓは制御性能に対する周波数重み関数であり、ロータを原点に固定することを制御指針として低周波領域でサーボ系の形になるように設計した。Ｙ軸の上部、下部それぞれに対する周波数重み関数はＸ軸と等しいものを使用している。
【０１０４】
前節で求めた補償器は、制御入力が吸引力であるため、そのまま実験機に適用することはできない。つまり、切換条件を基に力制御入力から制御電流への変換が必要である。本研究では、この切換を力制御入力の正負のみにより行う。つまり、力制御入力の正負により２対の電磁石の切換を行い、制御器が求めた吸引力ｆを再現する電流を逆算する。上側、下側はそれぞれ独立に切換を考えるので、上部側について詳細に述べる。
【０１０５】
1）Ｆｘｕ＞０のとき
図１５(c)において、電磁石１が吸引すればよい。そのときの力と電流ｉ１の関係式をｉ３＝０として
【数４０】
0005445887.tif 000042
【０１０６】
から逆算すると、次式となる。
【数４１】
0005445887.tif 000043
【０１０７】
2）Ｆｘｕ＜０のとき
図１５(d)において、電磁石３が吸引すればよい。そのときの力を電流ｉ３の関係式をｉ１＝０として
【数４２】
0005445887.tif 000044
【０１０８】
から逆算すると、次式となる。
【数４３】
0005445887.tif 000045
【０１０９】
下側についても同様に
1）Ｆｘｌ＞０のとき
【数４４】
0005445887.tif 000046
【０１１０】
2）Ｆｘｌ＜０のとき
【数４５】
0005445887.tif 000047
【０１１１】
この切換条件を用いて力制御入力を各電磁石の制御電流に変換することで制御を行う。Ｙ−Ｚ平面は、Ｘ−Ｚ平面と同様であり、以下のとおりである。
【０１１２】
上側
1）Ｆｙｕ＞０のとき
【数４６】
0005445887.tif 000048
【０１１３】
2）Ｆｙｕ＜０のとき
【数４７】
0005445887.tif 000049
【０１１４】
下側
1）Ｆｙｌ＞０のとき
【数４８】
0005445887.tif 000050
【０１１５】
2）Ｆｙｌ＜０のとき
【数４９】
0005445887.tif 000051
【０１１６】
５．２  制御電流の切換法
ゼロパワー切換制御の基本的な方法は、フライホイールシステムのロータ部を剛体とみなし、ロータが１つの電磁石に接近するとき、接近した電磁石の制御電流を零にして、反対側にのみ制御電流を供給して、ロータが完全に平衡点になれば、一切の電流（パワー）を流さないということである。さらに、制御器を設計するとき、制御入力は対向した電磁石間の吸引力として仮定する。求めた吸引力から切換法によって、実際の制御電流を求める。
【０１１７】
上側磁気軸受、下側磁気軸受の２対の電磁石が発生する磁気吸引力の合力をそれぞれＦｘｕ、Ｆｘｌとする。電磁石の磁気吸引力の合力は非線形性を有しているため、次式のように吸引力を制御入力と考えることで線形な制御対象として扱う。
【数５０】
0005445887.tif 000052
【０１１８】
ここで
ｋ：吸引力定数
ｉ：入力電流
ｘｇ：ノミナルエアギャップ
ｘ：ロータ変位
である。式(5.2)をみると、磁気吸引力は入力電流の二乗に比例し、電磁石とロータの距離の二乗に反比例するという形になっており、非常に強い非線形性を有している。従来の方法は、対向している両方の電磁石にバイアス電流を流すことによって、吸引力と電流の関係は小さな変位の領域で線形化できる。まず、入力電流と吸引力の関係は図１６(a)(b)の点線のようになっている。ただ、Ｆ１、Ｆ２は対向している電磁石の吸引力特性である。両方の電磁石にバイアス電流を流しておくと、図１６(c)の実線のように電磁石の吸引力の差が線形な形になる。
【０１１９】
この線形化による線形制御方法を適用し、コントローラを設計することが容易となる。しかし、ロータが安定に浮上しても、常時バイアス電流を流し続けなければならない。バイアス電流による渦電流損失は無視できず、回転抵抗や発熱などの問題を引き起こし、消費電力が多くなるという短所がある。また、この線形化法はモデリングにおいては、一般的な平衡点近傍でテイラー展開の一次近似を行うため、平衡点近傍では所望の制御性能が得られるが、大変位のときに制御性能が劣化するという問題もある。
【０１２０】
５．３  シミュレーション
実験の前段階として、制御器の閉ループ安定性を検証するため、ＭＡＴＬＡＢ  Ｓｉｍｕｌｉｎｋによるシミュレーションを行った。シミュレーションは図１７のようにＳｉｍｕｌｉｎｋにより構成した。より実験機システムに近づけるために、一度切換条件を基に力制御入力から制御電流への変換を行い、力制御入力に変換しなおしている。物理パラメータにより求めたモデルは制御入力を吸引力とみなしているため、制御電流から力制御入力に再度変換しなおす形となる。なお、制御電流は実際の実験と同様、１．５Ａのリミッタを設けている。また、サンプリング周期は実際の実験と同様１０ｋＨｚとし、補償器は双一次変換により離散化したものを用いた。Ｘ軸の上側に１０μｍのステップ目標値を０．０１ｓから加えている。このときの初期値応答を図１８に、各電磁石の制御電流を図１９に示す。
【０１２１】
図１８(a)において、オーバーシュートは少なく、整定時間も早く、良好な結果が得られていることがわかる。また、図１９から、制御入力が正負に変動する際に、電磁石の電流もしっかり切り換わっている。さらに、整定後は電流は流れておらず、ゼロパワー制御が行われていることがわかり、良好な結果となっている。
【０１２２】
５．４  実験
シミュレーションで用いた制御器を実際の磁気軸受フライホイールに適用し、９３Ｈｚ（５５８０ｒｐｍ）までの回転実験を行った。実験システムは図２０のように構成し、位置情報をセンサで取得し、Ａ／Ｄボードを通してＤＳＰに取り込み、力制御入力を計算した後、切換条件を基に各電磁石の制御電流に変換し、Ｄ／Ａボードを通して電磁石に出力される仕組みとなっている。サンプリング時間は０．１ｍｓとする。また、ＭＡＴＬＡＢからＤＳＰへはリアルタイムで情報を取得できるようになっており、図２１のようにｆｉｇｕｒｅにて各センサ情報を画面に表示することが可能である。インバータにはＭｙＷａｙ技研製のスピードコントローラを用いており、５１２ｒｐｍ／ｓの一定加速度を加えている。
【０１２３】
回転周波数０〜９０Ｈｚでの上側、下側におけるロータの軌跡および各電磁石の制御電流を図２３〜図４２に示す。また、上側の１対の電磁石の制御電流を拡大したものを図２２に示す。
【０１２４】
図２２より、制御電流においては、制御電流は切換が行われていることが確認できる。また、回転周波数に同期した半波波が制御電流に見られることがわかった。これは回転することで一種の調和外乱が印加されたような状態となり、補償器がそれを安定化しようと回転周波数に同期した制御電流を供給するためと考えられる。
【０１２５】
下側のロータの軌道は上側より大きい。これは、重心から下側磁気軸受までの距離が上側に比べ約１．８倍長いため、そしてアキシアル方向のＰＩＤ補償器との干渉のためと考えられる。
【０１２６】
浮上時には安定していたシステムが回転周波数９３Ｈｚで下側の振れ回りを大きくし、タッチダウンする結果となった。また、このモデルにおいては、回転周波数６０Ｈｚ、８０Ｈｚ付近で花びらのようなロータの軌跡が確認できる。モデルの剛性モード２次の前向き、後向き振れ回りが存在する。この結果から、モデルに連成を考慮しても固定制御器のロバスト性だけでは、ジャイロ効果による影響を低減する効果は少ないと考えられる。
【０１２７】
６．制御系設計と性能検証２
６．１  線形パラメータ依存モデル
実験では回転周波数８０Ｈｚ近傍でジャイロ効果が著しく発生することから、回転周波数によりシステムが変動する線形パラメータ依存モデル（ＬＰＶ）であると考えられる。ＭＡＴＬＡＢ  ＬＭＩ  ｃｏｎｔｒｏｌ  ｔｏｏｌｂｏｘを用いてＧＳ制御器を設計するにあたり、線形パラメータによって変動するシステム行列Ａを以下のように時変システムと線形時不変システム（ＬＴＩ）に分解して考える。
【数５１】
0005445887.tif 000053
【数５２】
0005445887.tif 000054
【数５３】
0005445887.tif 000055
【０１２８】
ここで、Ａ２２は、ジャイロモーメントの項を含むパラメータ変動行列である。回転周波数をスケジューリングパラメータとし、次の範囲で変化するものとする。
【０１２９】
ｐ＝ｆ  ｐ∈[０．１,１００]    (6.4)
６．２  ゲインスケジュール設計器設計
求めたＬＰＶモデルに対してＬＭＩ制御理論により制御系設計を行う。スケジューリングパラメータに対応する最適な制御器は、回転周波数０Ｈｚと１００Ｈｚに対応するＬＴＩ端点制御器（Ｋｍａｘ,Ｋｍｉｎ）を式(6.5)により凸補間することで求められる。
【数５４】
0005445887.tif 000056
【０１３０】
ＧＳ制御器の周波数応答を図４４に示す。モデルは４入力４出力で、さらに連成を考慮したことにより１６個の周波数応答が存在するが、ここでは例として、Ｘ軸上側からの入力に対する周波数応答が示されている。また、設計の際に用いた重み関数を図４３に示す。ロータの重心が磁気軸受の上側に存在し、バランスをとるために式(6.6)、式(6.7)のように、上側、下側で別の重みを用いている。
【数５５】
0005445887.tif 000057
【数５６】
0005445887.tif 000058
【０１３１】
図４４(a)、図４４(b)より、Ｙ軸の出力に対するコントローラのゲインが大きくなっていることが確認できる。これは、回転周波数が増加するとともに、連成方向へ影響に対する制御入力が大きくなることを意味する。
【０１３２】
６．３  シミュレーション
６．３．１  Ｈ∞制御器との比較
ＧＳ制御器の閉ループ安定性を検証するため、ジャイロ効果を考慮したモデルでシミュレーションを行った。初期値応答にてＨ∞制御器とＧＳ制御器の性能比較を行った。より実験機システムに近づけるために、一度切換条件を基に力制御入力から制御電流へと変換し、再び力制御入力に変換しなおしている。
【０１３３】
なお、制御電流は実験と同様に１．５Ａのリミッタを設けており、サンプリング周波数は実際の実験と同じ１０ｋＨｚで行い、補償器は双一次変換により離散化したものを用いた。回転周波数はパラメータの変動に対する応答をとるため、浮上時の０Ｈｚおよびジャイロ効果が著しく発生する８０Ｈｚ（４８００ｒｐｍ）とした。浮上時０Ｈｚでの初期値応答による出力変位の比較を図４５に、そのときの制御入力を図４６に示す。また、ジャイロ効果が著しく発生する８０Ｈｚ（４８００ｒｐｍ）での初期値応答による出力変位を図４７に、そのときの制御入力を図４８に示す。シミュレーション条件として、Ｘ軸の上側に１０μｍのステップ目標値を０．０１ｓから加えている。
【０１３４】
図４５、図４７より、Ｈ∞制御器、ＧＳ制御器ともに安定化しているが、収束の速さが異なる結果となった。とくにＧＳ制御器では、ジャイロ効果を考慮した影響により、Ｙ軸方向の下側の応答では振動が少なく良好な結果を得ることができた。また、図４５においても、Ｈ∞制御器に比べ素早く目標値に追従する良好な結果が得られた。ロバスト性を必要とするＨ∞制御器に比べ、ＧＳ制御器の端点コントローラ（Ｋｍｉｎ）は０Ｈｚで最適な設計がなされているためだと考えられる。
【０１３５】
６．３．２  テーブルを用いたＧＳ制御
実験システムへの実装を考えたとき、高周波領域を制御するために非常に短いサンプリング時間としていることと、４入力４出力システムであることに加え、剛性２次モードを制御する必要があるためモデルの低次元化を行うことができず、次数の大きい制御器を用いている関係から、サンプリングごとの離散化は計算処理が遅れるという問題がある。そこで、０Ｈｚから１００Ｈｚまで１０Ｈｚ間隔で１１個の制御器を予め凸補間により求めておき、オフラインでそれぞれを離散化した制御器をテーブルとして用いて、スケジューリングパラメータによって切り換えるという方法を提案する（図１３９参照）。実装する制御器の数はＤＳＰメモリ容量の関係より選択した。
【０１３６】
テーブルを用いることで、ＧＳ制御器の閉ループ安定性が損なわれないかを、初期応答により検証した。シミュレーション条件は前節と同じである。浮上時０Ｈｚの入出力を図４９、図５０に示す。また、制御器の切換時の応答を比較するため、７０Ｈｚ（４２００ｒｐｍ）と８０Ｈｚ（４８００ｒｐｍ）でそれぞれ凸補間して求めた制御器を用意しておき、０．１ｓ地点で制御器の切換を行ったときの入出力を図５１、図５２に示す。ジャイロの影響が強く出る８０Ｈｚ付近では若干の違いがあったが、安定性に問題がなく、良好な結果を得ることができた。
【０１３７】
６．４  実験
テーブルを用いたＧＳ制御器を用いて回転実験を行った。サンプリング周波数などの条件は、Ｈ∞制御器と同様とする。回転周波数０〜９０Ｈｚ（５４００ｒｐｍ）での上側、下側におけるロータの軌跡および制御入力を図５４〜図７１に示す。
【０１３８】
制御入力は、調和外乱に同期した半波が崩れたような波形となった。これは、回転周波数に依存する同期信号だけでなく、連成方向に対する入力が含まれるためであると考えられる。
【０１３９】
浮上時だけではなく、回転実験においてもロータは原点付近に固定されており、ジャイロモーメントによる影響を十分に制御できていることを実験にて確認できた。また、剛性モードおよび弾性モードの後向き振れ回り運動による花びらのような軌跡が発生していない。以上のことから、ジャイロ効果による影響に対してＧＳ制御器の有効性が示された。
【０１４０】
しかし、回転周波数９５Ｈｚ（５７００ｒｐｍ）でＨ∞制御器と同様にタッチダウンする結果となったことから、ジャイロ効果以外にも制御性能を劣化させる原因があると考えられる。
【０１４１】
６．５  固有値解析
実験結果を考察するため、実験に用いたＧＳ制御器について、パラメータ変動に対する閉ループ系固有値解析を行った。回転周波数を０Ｈｚから１００Ｈｚまで５Ｈｚずつ増加させたときのＨ∞制御器に対する複素平面を図５３(a)に、ジャイロ効果を考慮したＧＳ制御器に対する複素平面を図５３(b)に示す。
【０１４２】
ＧＳ制御器、Ｈ∞制御器はともに複素平面の左半面に位置しており、システムは安定である。しかし、プラントのパラメータ変動により、ＧＳ制御器の極がほぼ固定されているのに対して、いくつかの固有値が右平面に移動していることが確認される。図５３(a)と図５３(b)を比較すると、ＧＳ制御器がパラメータ変動に対しての安定性を補償していることを確認できる。しかし、回転周波数を２００Ｈｚにまで増加しても相変わらずシステムは安定しているため、実験において制御性能を劣化させている原因をモデル化する必要がある。
【０１４３】
７．弾性モデルによる検証
７．１  閉ループ系安定性の考察
剛性モデルでは、ジャイロ効果を無視した制御器を用いても、考慮した制御器を用いても、閉ループはともに安定であった。しかし、回転実験では、１００Ｈｚ付近で不安定化している。この不一致の問題に対して、弾性曲げモードも考慮して考察する。弾性モデルの導出は有限要素法によって行う。フライホイールを集中質量として、ロータを図７２のように８要素に分割する。要素２と３、要素７と８の間に電磁石が配置されている。ここで、減衰と不釣合い振動は考慮しないものとする。図７２のモデリングについて、ジャイロ効果を考慮したロータの運動方程式は次式となる。
【数５７】
0005445887.tif 000059
【０１４４】
Ｍは質量行列、Ｇはジャイロ反対称行列、Ｋは剛性行列、Ｕは制御入力ベクトル、Ｔは電磁石の設置場所を表す行列（Ｍ,Ｇ,Ｋ∈Ｒ３６×３６,ＴinＲ３６×４）、ｑは状態ベクトルである（ｑinＲ３６×１）。まず、フリーフリーロータに対して、固有振動数と軸回転数の関係を図７３で示す。曲げ一次モードの前向き振れ回りモードと後向き振れ回りモードは回転後大きな分岐をすることがわかった。このことから、曲げ一次弾性モードの後向き振れ回りモードは、ジャイロ効果の影響を強く受けて１００Ｈｚ付近まで下がってきている。この結果、運搬線と後向き振れ回りモードの交差周波数は１１０Ｈｚ付近で、ロータが１００Ｈｚを越えて回転すると、後向き曲げ共振点に接近していることがわかった。
【０１４５】
７．２  モード分離と低次元化
１００Ｈｚを越えた高速回転の実現には、ＦＥＭモデルの剛性モードと曲げ一次モードのジャイロ効果による影響が重要となる。ここで、式(7.1)には反対称のジャイロ行列Ｇが存在するので、標準形のモード解析法を簡単に使うことができない。モード分離を行うため、新しいベクトル
【数５８】
0005445887.tif 000060
【０１４６】
を導入すると次式となる。
【数５９】
0005445887.tif 000061
【数６０】
0005445887.tif 000062
【０１４７】
式(7.2)の固有値問題を次式により考察する。
【数６１】
0005445887.tif 000063
【０１４８】
式(7.4)は標準形式に似ているが、ＫとＧの特性からＧ＊は相変わらず反対称行列になっている。しかし、このような固有値問題の解はｎ組の共役純虚数
【数６２】
0005445887.tif 000064
【０１４９】
とそれに対応しているｎ組の共役複素数の固有ベクトル
【数６３】
0005445887.tif 000065
【０１５０】
を持っていることがわかっている。そして、ｓ＝ｉｗとｚ＝ｕ＋ｉｕを式(7.4)に代入し、実部と虚部を分けることにする。最後に、ｕとｖに関連する式は次のような２つの固有値問題形式が得られる。
【数６４】
0005445887.tif 000066
【０１５１】
ここで、Ｋ＊は正定対称行列である。式(7.5)はｎ組の重要なλｒを持ち、ｕｒとｖｒは同じλｒに対応していることがわかる。次に、Ｍ＊に対してコレスキー分解を行うと（Ｍ＊＝ＲＲ−Ｔ）、次の標準固有値問題になる。
【数６５】
0005445887.tif 000067
【０１５２】
ここで、Φは対称正定行列である。λｒに対応している固有ベクトルはＰ＝[ｕｒ ｖｒ]となる。座標変換式ｚ＝Ｒ−ＴＰξを設定し、式(7.2)に代入すると、モード分離した方程式(7.7)が得られる。ξ＝｜ηｒ ζｒ｜はモード座標である。Θは式(7.7)で表すように、ブロック対角行列となる。
【数６６】
0005445887.tif 000068
【数６７】
0005445887.tif 000069
【０１５３】
式(7.8)に基づいて、剛性モードと曲げ一次弾性バックモードを含むモデルの低次元化は次式となる。
【数６８】
0005445887.tif 000070
【０１５４】
ただし、Ａｒ∈Ｒ１０×１０,Ｂｒ∈Ｒ１０×４,Ｃｒ∈Ｒ４×１０,Ｄｒ∈Ｒ４×４,Ｕは制御入力ベクトル（Ｕ∈Ｒ４×１）である。低次元化した前後のモードのボード線図を図７４で示す。
【０１５５】
７．３  シミュレーション
前節で作成した剛性モードと曲げ一次弾性バックモードを含む低次元化モデルが、実験と一致するモデルであるかを検証する。制御器は剛性モードのみを考慮して設計したＧＳ制御器を用いた。サンプリング周波数は１０ｋＨｚとし、Ｘ軸上側に１０μｍのステップ目標値を与えたときの初期値応答シミュレーションと実験の比較を図７５に示す。
【０１５６】
弾性モデルを用いることで、実験結果と同じ高周波がシミュレーションで励起されるようになった。浮上時における弾性モードの固有振動数と一致することから、低次元化モデルに曲げ一次弾性バックモードが含まれていることが初期値応答からも確認できる。しかし、低周波成分や収束時間が大幅に異なり、実験システムよりも安定化しやすいモデルであるといえる。この原因がインダクタンスや渦電流損失などアクチュエータの諸元によるものと考え、図７６に示すようにＧＳ制御器によって求められた制御電流にαという係数を掛ける。αを加えた閉ループでのシミュレーションを図７７に示す。低周波成分は完全に一致はしていないが、オーバーシュートの振幅や立ち上がり時の振幅が実験結果に近づく有効な結果となった。なお、αの値は試行錯誤的に０．２２とした。
【０１５７】
７．４  固有値解析
係数αを加えたときと、加えないときの安定性の比較を、図７８に示す閉ループの固有値解析にて行った。図７８によりαを加えないときは固有値が実軸からかなり離れた地点に存在するため、安定であった。しかし、実験システムは安定ではあるものの、極は実軸付近に存在する。この現象は、シミュレーションと実験の不一致の原因の１つであったと考えられる。αを加えることで、極は右側に移動し、実験システムに近い極配置となった。また、前節のシミュレーション結果では、剛性２次モードは制御されていて、時系列応答には励起されていなかった。これについても、係数αをっ加えたときの極解析では剛性２次モードの極が安定化され、実験結果と一致する。周波数を実験でタッチダウンする１００Ｈｚとし、Ｘ軸上側に１０μｍのステップ目標値を与えたときの初期値応答を図７９に示す。これまでのシミュレーションで１００Ｈｚにおいても安定していたＧＳ制御器を用いたにもかかわらず、発散する結果となった。以上のことから、係数αを加えることで、低次元化弾性モデルは、実システムに近い応答を示すことが確認できる。
【０１５８】
８．結論
本論文では、磁気軸受系型電力貯蔵フライホイールの制御に関する研究を行った。慣性力を増加させることを目的にフライホイールがロータに取り付けられているために、ジャイロ効果による前向き振れ回り運動や後向き振れ回り運動の発生、およびそれらの固有振動の変化が高速回転を妨げている。
【０１５９】
そこで、本研究では、まず、ジャイロ効果の特性を考慮したモデリングを行った。プラントのＡ行列に回転周波数ωを定数として考慮し、Ｈ∞制御器を設計した。浮上時には安定していたシステムが回転周波数９３Ｈｚで下側の振れ回りを大きくし、タッチダウンする結果となった。下側が先にタッチダウンした原因は、重心から下側磁気軸受までの距離が上側に比べ約１．８倍長いため、下側のロータの軌道は上側より大きくなるためである。また、回転周波数６０Ｈｚ、８０Ｈｚ付近で、ロータは花びらのような軌跡を描いた。モデルの剛性モード２次の前向き、後向き振れ回りが存在することが確認された。この結果から、モデルに連成を考慮した固定制御器のロバスト性だけでは、ジャイロ効果による影響を低減する効果は少ないと考えられる。
【０１６０】
次に、回転周波数によりシステムが変動する線形パラメータ依存モデル（ＬＰＶ）であると考え、ＧＳ制御器を設計し、シミュレーションを行った。浮上時では、Ｈ∞制御器、ＧＳ制御器ともに安定化しているが、収束の速さが異なる結果となった。とくに、ＧＳ制御器では、ジャイロ効果を考慮した影響により、Ｙ軸方向の下側の応答においても振動が少なく、素早く目標値に追従する良好な結果が得られた。ロバスト性を必要とするＨ∞制御器に比べ、ＧＳ制御器の端点コントローラ（Ｋｍｉｎ）は０Ｈｚで最適な設計がなされているためだと考えられる。
【０１６１】
実験システムへの実装を考えたとき、高周波領域を制御するために非常に短いサンプリング時間としていることと、４入力４出力システムであることに加え、剛性２次モードを制御する必要があるため、モデルの低次元化を行うことができず、次数の大きい制御器を用いている関係から、サンプリングごとの離散化は計算処理が遅れるという問題が発生した。そこで、０Ｈｚから１００Ｈｚまで１０Ｈｚ間隔で制御器を予め凸補間により求めておき、オフラインでそれぞれを離散化した制御器をテーブルとして用いて、スケジューリングパラメータによって切り換えるという方法を提案した。実装する制御器の数はＤＳＰのメモリ容量の制約より選択した。テーブルを用いることで、ＧＳ制御器の閉ループ安定性が損なわれないかを、初期応答により検証した。ジャイロの影響が強く出る８０Ｈｚ付近での切換では、収束速度に若干の違いがあったが、安定性に問題がなく、良好な結果を得ることができた。
【０１６２】
浮上時だけではなく、回転実験においても、ロータは原点付近に固定されており、ジャイロモーメントによる影響を十分に制御できていることを実験にて確認できた。また、剛性モードおよび弾性モードの後向き振れ回り運動による花びらのような軌跡が発生していない。以上のことから、ジャイロ効果による影響に対して、ＧＳ制御器の有効性が示された。しかし、回転周波数９５Ｈｚ（５７００ｒｐｍ）で、Ｈ∞制御器と同様にタッチダウンする結果となった。
【０１６３】
実験結果を考察するため、実験に用いたＧＳ制御器について、パラメータ変動に対する閉ループ系固有値解析を行った。ＧＳ制御器、Ｈ∞制御器の極はともに複素平面の左半面に位置しており、システムは安定である。しかし、プラントのパラメータ変動した際、ＧＳ制御器の極は固定されているのに対して、Ｈ∞制御器ではいくつかの極が右平面に移動していることが確認された。しかし、ジャイロ効果を考慮した剛性モデルの回転周波数を２００Ｈｚにまで増加しても、閉ループの安定性は保証される。そこで、実験において制御性能を劣化させていると考えられる曲げ一次弾性バックモードを含む低次元化モデルを用いて検証を行った。弾性モデルを用いることで、実験結果と同じ高周波がシミュレーションで励起されるようになった。また、浮上時における弾性モードの固有振動数と一致することから、低次元モデルに曲げ一次弾性バックモードが含まれていることが初期値応答からも確認できる。しかし、低周波成分や収束時間が大幅に異なり、実験システムよりも安定化しやすいモデルであるといえる。この原因がアクチュエータの諸元によるものと考え、ＧＳ制御器によって求められた制御電流にαという係数を掛け、再びシミュレーションを行った。低周波成分は完全に一致はしていないが、オーバーシュートの振幅や立ち上がり時の振幅が実験結果に近づく有効な結果となった。
【０１６４】
係数αを加えたときと、加えないときの安定性の比較を、閉ループの固有値解析にて行った。αを加えないときは、固有値が実軸からかなり離れた地点に存在するため、安定であった。しかし、実験システムは安定であるものの、極は実軸付近に存在する。この現象は、シミュレーションと実験の不一致の原因の１つであると考えられる。αを加えることで、極は右側に移動し、実験システムに近い極配置となった。また、前節のシミュレーション結果では、剛性２次モードは制御されていて、時系列応答には励起されていなかった。これについても、係数αを加えたときの極解析では、剛性２次モードの極が安定化され、実験結果と一致する。また、初期値応答を行ったところ、これまでのシミュレーションで１００Ｈｚにおいても安定していたＧＳ制御器を用いたにもかかわらず、発散する結果となった。以上のことから、係数αを加えた低次元化弾性モデルは、実システムに近いモデルであるといえる。
【０１６５】
今後の課題をしては、弾性モデルに係数αを加えたモデルを用いて、より良い性能を有する制御器の設計が必要と思われる。
【０１６６】
ジャイロ効果を考慮した磁気軸受系のゼロパワー制御
１．３  本研究の目的
本研究の目的は、バイアス電流を流さないで、非線形を有したまま、非線形制御理論に基づいて、ゼロパワー切換制御法によって安定化を図り、かつ、消費電力を低減して、ロータが高速回転できる制御器を設計し、有効性を検証する。
【０１６７】
１．３．１  モデリング簡単化
磁気軸受系では、運動方程式を求めるとき、ロータ系を考慮して、重心に関する状態方程式を得ることが一般である。しかし、実験装置においては、センサの観測点は重心の位置と異なっており、ロータの上、下部電磁石に近いところにおかれている。それで、制御器を設計するとき、センサからのフィードバック状態量は重心に関するものに変換しなければならないので、演算時間もかかり、精度も下がる。このような欠点に対して、重心位置でのロータ変位と電磁石に対応する点でのロータ変位、回転角と変位間の関係により、すべての電磁石に対してロータ重心に加えている電磁力と回転モーメントをロータ上、下部の電磁石により分離して考えると、最後に求めた状態方程式が簡単化される。すべての状態量はセンサからの変位信号のみに関連するものである。
【０１６８】
１．３．２  ゼロパワー切換制御法
ゼロパワー切換制御の基本的な方法はフライホイールシステムのロータ部を剛体とみなし、ロータが１つの電磁石に接近するとき、接近した電磁石の制御電流を零にして、反対側にのみ制御電流を供給して、ロータが完全に平衡点になれば、一切の電流（パワー）を流さないといういうことである。さらに、制御器を設計するとき、制御入力は対向した電磁石間の吸引力として仮定する。求めた吸引力から切換法によって、実際の制御電流を求める。
【０１６９】
１．３．３  非線形制御器の設計と性能検証
本論文は０．５ＫＷｈ級磁気軸受フライホイールシステムに対して、評価関数の最適化を目指した制御理論に基づいて、ＬＱＲ法により制御器を設計して、実験で検証する。ＬＱＲ法のレギュレータ制御器を設計した上で、定常偏差をなくすため、サーボ系を付加したＬＱＩ制御器を設計する。実験の結果から、サーボ系を付加した場合は、回転速度は１００Ｈｚ付近に達した。しかし、ジャイロ効果を考慮すると、システムが回転速度という変動パラメータに依存するＬＰＶシステムになったため、普通のＬＴＩコントローラでは厳密には補償できない。したがって、ゲインスケジュール型スライディング制御を提案する。実験の結果、ジャイロ効果の影響を補償したが、ＬＱＩコントローラと同じ所でタッチダウンしてしまった。つまり、ある周波数に到達すれば、弾性モードを考慮しないと、システムが不安定になることがわかった。そのため、浮上状態の加振実験により、回転周波数と固有振動数の関係を求める。この結果より、回転周波数が１００Ｈｚ付近から、大きなピークが発生している。すなわち、１００Ｈｚは弾性モードのバックワード共振点と認められる。本論文では、弾性モデルの次数が高いため、非モデルベースのＰＩＤ制御で１３０Ｈｚまで増速できた。しかし、ＰＩＤ制御の弾性モデルのシミュレーションと実験結果が一致しない問題が発生しており、これについては今後考察する。
【０１７０】
制御対象は、図１〜図３に示す磁気軸受装置である。
【０１７１】
３．制御理論
本磁気軸受系では、バイアス電流を完全に回避できるため、電磁力と変位間の非線形関係を有したまま、非線形理論に基づいて、制御器を設計する。
【０１７２】
３．１  線形２次形式評価模範による制御系設計
３．１．１  最適レギュレータ
評価関数を設け、それを最小とするように状態フィードバックによる最適制御入力を決定する設計法について述べる。可制御で線形な時不変システム
【数６９】
0005445887.tif 000071
【０１７３】
のもとで、ある評価関数を最小にする制御則が線形な状態フィードバック制御則となるためには、評価関数を以下のように２次形式とする必要がある。
【数７０】
0005445887.tif 000072
【０１７４】
ここで、制御目的の重み行列Ｑは非負定な対称行列、制御入力の重み行列Ｒは正定対称行列とする。２次形式評価関数式(3.2)で終端時間ｔｆ→∞とした式(3.3)の評価関数
【数７１】
0005445887.tif 000073
【０１７５】
を最小にする最適制御入力は
【数７２】
0005445887.tif 000074
【０１７６】
で与えられる。ここで、Ｐは式(3.5)のリカッチ型代数方程式の正定対称行列であり、一意に決定される。
【数７３】
0005445887.tif 000075
【０１７７】
式(3.3)からわかるように、重み行列の定性的な設計指針としては、Ｑを大きくとれば制御性能が増し、Ｒを大きくとれば制御入力を小さくできることになる。目的関数と制御入力の間には、トレードオフの関係がある。
【０１７８】
制御の目的関数として何を採用するかは、制御仕様によるが、振動系の制御問題では、位置エネルギ：(１／２)ｋｘ２、運動エネルギ：(１／２)ｍｘ２を採用する場合がある。また、たとえば加速度を抑制することを目的とする場合には、加速度は一般に
ｙ＝Ｃｘ＋Ｄｕ    (3.6)
と、制御入力の直達項を含めて表されるから、
【数７４】
0005445887.tif 000076
【０１７９】
を評価関数に取り入れる必要がある。この場合には、クロスターム項を加えた次に２次形式評価関数を採用する。
【数７５】
0005445887.tif 000077
【０１８０】
この場合の最適制御入力は以下のように与えられる。
【数７６】
0005445887.tif 000078
【０１８１】
Ｐは以下のリカッチ方程式の一意正定対称行列である。
【数７７】
0005445887.tif 000079
【０１８２】
評価関数に２次形式を用いて線形な状態フィードバック制御則を導き、状態を原点に戻す閉ループ系を構成する制御系を線形２次形式レギュレータという。
【０１８３】
３．２  リアプノフの安定性理論
３．２．１  内部安定性の定義
非線形時変システムのふるまいは、線形時不変システムのように指数関数で表されるわけではなく、一般に複雑である。したがって、安定性の精密な定義を行えば、多様な概念を考えることができる。連続系の場合と同様、離散値系にもリアプノフの安定の概念が適用される。今、次のようなシステムを考える。ここでは、そのうちのいくつかを説明する。
【０１８４】
ｘ＝ｆ(ｘ,ｔ)    （3.11)
ここで、ｘはｎ次元状態ベクトル、ｆは時間ｔとｘに関して連続な関数である。式(3.11)は任意の初期時刻ｔ０と初期状態ｘ０に対して、すべての時刻にわたって一意解を有する。また、ｔ≧０における値をｘ(ｔ,ｘ０,ｔ０)と表す。ただし、ｘ(ｔ０,ｘ０,ｔ０)＝ｘ０である。このシステムにおいて
ｘ(ｔ,ｘ０,ｔｅ)≡ｘｅ    (3.12)
が成り立つ状態ｘｅを平衡状態と呼ぶ。平衡状態とは、式(3.11)の定数解であるから、すべてのｔで
ｆ(ｔ,ｘｅ)＝０    (3.13)
を満たすものである。この式は非線形時変方程式であるから、定数解をいつも持つとは限らない。また、持つとしても唯一とは限らない。システムの平衡状態ｘｅの周りの振る舞いを考えるとき、ｘｅが状態空間の原点ｘ＝０であるとしても一般性は失われない。もしｘｅ≠０ならば、状態空間の座標系を移動して、新たなｆ関数に対する微分方程式は原点ｘ＝０を平衡状態として持ち、その解の振る舞いは式(3.11)と同じである。
【０１８５】
安定性
安定性とは初期状態ｘ０を平衡状態ｘ＝０に近く取ることによって、解ｘ(ｔ,ｔ０,ｘ０)をｘ＝０の近くに止めることができるという概念である。
【０１８６】
定義１：任意の正数εと初期時刻ｔ０に対して、ある正数δ(ｔ０,ε)が存在し、‖ｘ０‖≦δなる初期状態ｘ０について、すべての時刻ｔ≧ｔ０で、‖ｘ(ｔ,ｔ０,ｘ０)‖≦εが成立するとき、平衡状態ｘ＝０は安定である。図８０(a)を参照する。
【０１８７】
一様安定性
安定性の定義では、解ｘ(ｔ)が平衡状態の近傍に止まるためには、どれくらい初期状態が平衡状態に近くなければならないかの距離が初期状態に依存して、初期時刻にかかわらず、だいたい同じであることを一様と定義する。
【０１８８】
定義２：定義１において、正数δとして、初期時刻に独立なものδ(ε)が存在するとき、平衡状態ｘ＝０は一様安定であるという。微分方程式(3.11)が時不変、つまりその右辺の関数ｆが時間に独立で、非線形関数ｆ(ｘ)ならば、その解は初期時刻に存在しない。ゆえに、このときの安定性と一様安定性は等価である。
【０１８９】
一様漸近安定性
一様安定性に加えて、解が平衡状態に時間とともに収束するとき、一様漸近安定性という。
【０１９０】
定義３：一様安定したうえで、ある正数ｒと任意的な小さい正数μに対して時間Ｔ(μ,ｒ)が存在し、‖ｘ０‖≦ｒなるすべての初期状態ｘ０と任意の初期時刻ｔ０について、すべての時刻ｔ≧ｔ０＋Ｔで
‖ｘ(ｔ,ｔ０,ｘ０)‖≦μ    (3.14)
が成立するとき、一様漸近安定であるという。図８０(c)を参照する。
【０１９１】
大域的一様安定性
解が平衡状態に収束するとき、その収束解の出発する初期状態の集合を吸収領域と呼ぶ。この吸収領域が状態空間の全体であることを大域的という語で表す。図８０(b)を参照する。
【０１９２】
定義４：平衡状態ｘ＝０は、任意の大きな正数ｒに対する初期状態ｘ０の集合‖ｘ０‖≦ｒを考えても、一様漸近安定である。
【０１９３】
任意の正数ηに対してある正数ξ(η)が存在し、‖ｘ０‖≦ηなる任意の初期状態ｘ０とすべての初期時刻ｔ０について、すべての時刻ｔ≧ｔ０で、
‖ｘ(ｔ,ｔ０,ｘ０)‖≦ξ    (3.15)
が成立するとき、大域的一様漸近安定であるという。
【０１９４】
大域的一様漸近安定の定義は、単に一様漸近安定性の吸収領域を大域的にしただけでなく、条件が追加されている。この条件は一様有界性と呼ばれている。一様有界性は一様安定性に似ているが、一様安定性では、解の存在範囲に対する初期状態の範囲を考えて、一様有界性では、初期状態の範囲に対する解の存在範囲を見ている点である。したがって、一様有界であって一様安定でない平衡状態や、その逆のものも存在する。
【０１９５】
３．２．２  リアプノフの安定判別法
微分方程式の平衡状態が、前項で述べた安定性のどの性質を持つかを調べる方法としてリアプノフの安定判別法と呼ばれる方法がある。解の平衡状態からの距離を時間的変化が状態のスカラ関数を用いて評価されるもので、解を求めることなく、解析的に安定性を調べることができる。このスカラ関数は、平衡状態からの距離の概念を拡張したもの、または、状態が物理的変数ならば、その変数に対応する要素が持つポテンシャルエネルギの概念を拡張したものと見なすことができる。
【０１９６】
今、時間ｔと状態ｘを変数とするスカラ関数(ｔ,ｘ)を考える。ｔとｘに関して連続微分可能で、すべてのｔでＶ(ｔ,０)＝０とする。そして、この関数に対して、以下の性質を定義する。
【０１９７】
α(σ)は連続かつσに関して非減少の１変数スカラ関数で、α(０)＝０かつσ＞０に対してα(σ)＞０である。すべてのｔ,ｘについて
α(‖ｘ‖)≦Ｖ(ｔ,ｘ)    (3.16)
が成立するようなα(σ)が存在するとき、Ｖ(ｔ,ｘ)は正定であるという。
【０１９８】
式(3.16)を満たすα(σ)で
【数７８】
0005445887.tif 000080
【０１９９】
となるものが存在するとき、Ｖ(ｔ,ｘ)は半径方向に非有界であるという。
【０２００】
α(σ)と同様の性質を持つ関数β(σ)が存在し、すべてのｔ,ｘについて
Ｖ(ｔ,ｘ)≦β(‖ｘ‖)    (3.18)
が成立するとき、Ｖ(ｔ,ｘ)はデレクセントであるという。これは、‖ｘ‖が減少することについて、Ｖ(ｔ,ｘ)も一様に減少するという意味である。
【０２０１】
次に、時刻ｔで状態がｘであるとして、その点に式(3.11)の解に沿ったＶ(ｔ,ｘ)の時間微分Ｖ(ｔ,ｘ)を考える。Ｖ(ｔ,ｘ)は
【数７９】
0005445887.tif 000081
【０２０２】
と表され、次の性質を定義する。
【０２０３】
すべてのｔ,ｘについて
Ｖ(ｔ,ｘ)≦０    (3.20)
であるとき、Ｖ(ｔ,ｘ)は準負定であるという。
【０２０４】
α(σ)と同様の性質を持つ関数γ(σ)が存在し、すべてのｔ,ｘについて
Ｖ(ｔ,ｘ)≦−γ(‖ｘ‖)    (3.21)
が成立するとき、Ｖ(ｔ,ｘ)は負定であるという。
【０２０５】
以上の準備のもとに、リアプノフの考えに基づく安定定理を述べる。
【０２０６】
以下の性質を持つ関数Ｖ(ｔ,ｘ)が存在すれば、式(3.11)の平衡状態ｘ＝０について、次のことが言える。
【０２０７】
Ｖ(ｔ,ｘ)が正定で、Ｖ(ｔ,ｘ)が準負定ならば、安定である。
【０２０８】
Ｖ(ｔ,ｘ)が正定かつデレクセントで、Ｖ(ｔ,ｘ)が準負定ならば、一様安定である。
【０２０９】
Ｖ(ｔ,ｘ)が正定かつデレクセントで、Ｖ(ｔ,ｘ)が負定ならば、一様漸近安定である。
【０２１０】
Ｖ(ｔ,ｘ)正定かつデレクセント、かつ半径方向に非有界で、Ｖ(ｔ,ｘ)が負定ならば、大域的一様漸近安定である。
【０２１１】
定義から明らかなように、安定性、一様安定性、一様漸近安定性は平衡状態近傍の解の性質である。したがって、この定理の条件は、平衡状態の近傍で成立すれば十分である。この定理が重要なのは、非線形微分方程式を解くことなく、平衡状態を知ることができる点である。これらの定理の条件を満たす関数をリアプノフ関数という。
【０２１２】
３．２．３  コントロールリアプノフ関数（ＣＬＦ）
次に、前項のリアプノフ関数をより拡張した、コントロールリアプノフ関数について述べる。
【０２１３】
次のような時不変システムを考える。
【０２１４】
ｘ＝ｆ(ｘ,ｕ)，ｘ∈Ｒｎ，ｕ∈Ｒ，ｆ(０,０)＝０    (3.22)
ここで、制御入力ｕにフィードバック制御則αを代入して得られる閉ループ系
ｘ＝ｆ(ｘ,α(ｘ))    (3.23)
が平衡点ｘ＝０において大域的漸近安定となるようなα(ｘ)を設計することを考える。Ｖ(ｘ)をリアプノフ関数の候補とすると、式(3.23)が安定であるためには、Ｖ(ｘ)≦−Ｗ(ｘ)を満たす必要がある。ここで、Ｗ(ｘ)は正定関数である。すべてのｘ∈Ｒｎに対して
【数８０】
0005445887.tif 000082
【０２１５】
を満たすようなαを見つけるという問題に帰着する。しかし、この問題は難しい問題である。なぜなら、式(3.22)を安定化する制御則は存在するかもしれないが、Ｖ(ｘ)やＷ(ｘ)の選び方によっては式(3.24)を満たさない可能性があるからである。そこで、望ましいＶ(ｘ)やＷ(ｘ)を選択することができるシステムは「コントロールリアプノフ関数を所有する」と言われる。以下のその概念を明確にする。
【０２１６】
定義２．４  連続正定関数で半径方向に非有界な関数Ｖ：Ｒｎ→Ｒ＋が次式を満たすとき、Ｖ(ｘ)をコントロールリアプノフ関数と呼ぶ。
【数８１】
0005445887.tif 000083
【０２１７】
また、次式で表されるようなシステム
ｘ＝ｆ(ｘ)＋ｇ(ｘ)ｕ，ｆ(０)＝０    (3.26)
の場合、コントロールリアプノフ関数の不等式は以下のようになる。
【数８２】
0005445887.tif 000084
【０２１８】
もしＶ(ｘ)を式(3.26)のコントロールリアプノフ関数とすると、ソンタックの公式（Sontag's formula）により、すべてのｘ≠０に対して連続な安定化制御則α(ｘ)は以下のように一意に与えられる。
【数８３】
0005445887.tif 000085
【０２１９】
ここで、式(3.27)は
【数８４】
0005445887.tif 000086
【０２２０】
のときのみ条件を満たすことに注意しなければならない。この場合、式(3.28)は
【数８５】
0005445887.tif 000087
【０２２１】
となる。
【０２２２】
さらに、式(3.26)に対する安定化制御則の特徴として、α(ｘ)はコントロールリアプノフ関数Ｖ(ｘ)が次のような条件を満たすときのみ、ｘ＝０で連続となる。どのε＞０に対してもｘ≠０のとき、｜ｘ｜＜δを満たすδ(ε)＞０が存在し、｜ｕ｜＜εなるｕが次式を満たすとき
【数８６】
0005445887.tif 000088
【０２２３】
制御系設計を用いる際のコントロールリアプノフ関数の主な欠点は、ほとんどの非線形システムはコントロールリアプノフ関数が未知であるということである。また、それと同じくらい安定化制御則を設計するのは困難である。しかし、バックステッピング法を用いれば、コントロールリアプノフ関数を見つけることと、安定化制御則を設計するという２つの課題を同時に解決することができる。この手法を始めるにあたって、少なくともまず、スカラシステムに対してコントロールリアプノフ関数Ｖ(ｘ)と安定化制御則α(ｘ)を見つけられるようにならなければならない。幸運なことに、スカラシステムに対してはＶ(ｘ)＝ｘ２／２がいつも妥当なコントロールリアプノフ関数であり、不等式(3.27)を容易に満たす。
【０２２４】
３．３  スライディングモード法
制御系の構造を変える理論は可変構造制御理論と呼ばれ、その中で最も理論的に体系化されているのがスライディングモード制御理論である。その特徴として優れたロバスト制御系が構築できるといったことが挙げられ、今日システムの安定化やサーボ系をはじめとして、様々な制御目的に適応できるまでに理論が進んでいる。
【０２２５】
３．３．１  スライディングモードの存在条件
スライディングモード制御では、システムの位相空間上での挙動は到達モードとスライディングモード２つに分けられ、状態が超平面の方向に向かい、そして超平面に辿り着く条件到達条件と呼ばれる。この到達条件を満足していれば、到達モードは有限時間の間に超平面に到達し、マッチング条件を満たす外乱に対しロバストになる。これらのことから、スライディングモード制御を実現するためには以下の手順が必要となる。
【０２２６】
(1) 与えられたプラントより低次数で、一般には（システム次元：ｎ)（入力次元は：ｍ）の望ましいシステムダイナミクスを表すための切換面を設計すること。
【０２２７】
(2) 切換面の外にあるどのような状態量も、有限時間に切換面に到達するような可変構造制御入力ｕ(ｔ)を設計すること。
【０２２８】
これらは、「スライディングモードコントローラを設計する」ことに他ならない。スライディングモード到達条件を満足するために、様々なアプローチが提案されているが、本研究ではリアプノフ関数法を適用する。今、リアプノフ関数の候補を次のように選ぶ。
【０２２９】
Ｖ＝σＴ・σ／２    (3.32)
このとき、包括的なスライディングモード到達条件は次式で与えられる。
【０２３０】
Ｖ＝σＴ・σ≦０    (3.33)
この到達則を満足させるように可変構造制御入力を決定することが重要である。
【０２３１】
３．３．２  等価制御法による解析
スライディングモードが存在しているとき、システムは非線形性の最も強いスイッチング入力のため、解析などが著しく困難となる。このスイッチング入力を連続入力で置き換えることにより、解析および設計の見通しがよくなる。システムは
ｘ＝Ａｘ＋Ｂｕ
σ＝Ｓｘ    (3.34)
と定義し、Ａ∈Ｒｎで、Ｂ∈ＲｍとＳ∈Ｒｍをプルランクとする。式(3.35)において、スライディングモードが存在すると
σ＝０    (3.35)
式(3.36)の関係を得られる。
【０２３２】
σ＝Ｓｘ＝Ｓ(Ａｘ＋Ｂｕ)＝０    (3.36)
よりｄｅｔ(ＳＢ)≠０ならば等価制御入力が以下のように求まる。
【０２３３】
ｕｅｑ＝−(ＳＢ)−１ＳＡｘ    (3.37)
この等価制御入力ｕｅｑはシステムの理想的な連続制御入力であり、ｕｅｑを元のシステムの式(3.34)に代入することによって、入力の数だけ低次元化された以下のようなシステムが得られる。このシステムを等価制御系という。
【０２３４】
ｘ＝{Ｉ−Ｂ(ＳＢ)−１Ｓ}Ａｘ    (3.38)
この閉ループの固有値は伝達関数Ｓ(ｓＩ−Ａ)−１Ｂの（ｎ−ｍ）個の零点とｍ個の原点極からなっている。スライディングモード制御の考え方法を示すために、次式のように表すことができる線形時不変系の制御対象を考える。簡単のため、制御対象は１入力（ｍ＝１）、システムの次数はｎ次とし、正準形に変換されているものとする。
【数８７】
0005445887.tif 000089
【０２３５】
スライディングモード制御系の構造は、次式のような切換関数σによって支配される。
【数８８】
0005445887.tif 000090
【０２３６】
ここで、正則行列ｓ２は任意であるので、ｓ２＝Ｉｍとする。切換関数は、位相平面をσの符号が異なるように分割する。この切換関数の符号によって制御入力を切り換えることにより、システムの状態を切換超平面上に拘束し、入力の数だけ低次元化された望ましいシステムを構成する。切換超平面上に状態を拘束されたシステムを等価制御系といい、この系について考えることにより、スライディングモードにあるシステムの動特性を解析できる。
【０２３７】
超平面上においては式(3.40)に示す切換関数σは０であるため、
ｘ２(ｔ)＝−ｓ１ｘ１(ｔ)    (3.41)
となり、ｘ２がｘ１のサブシステムへの入力となる。すなわち、式(3.39)は
ｘ１(ｔ)＝(Ａ１１−Ａ１２ｓ１)ｘ１(ｔ)σ
＝ｓ１ｘ１(ｔ)＋ｓ２ｘ２(ｔ)＝０    (3.42)
となり、低次元化された遅いｘ１のシステムと早いシステムであるσ＝０とに分離された、マッチング条件を満たすパラメータ変動や外乱の影響を受けないシステムとなる。上式に示されたように、スライディングモードの動特性は行列Ｓにより決められる。
【０２３８】
３．３．３  スライディングモード到達条件とチャタリング低減
ここでは、超平面上にシステムの状態を到達させ、拘束するための条件、すなわちスライディングモード到達条件を満たす制御入力を求める。先にも示したようにリアプノフ関数を次式のように置く。
【０２３９】
Ｖ＝σ２／２    (3.43)
このとき、スライディングモード到達条件は、次式を満足すればよい。
【０２４０】
ｖ＝σσ＜０    (3.44)
切換関数は次式で表すことができる。
【数８９】
0005445887.tif 000091
【０２４１】
ただし、ｚ(ｔ)は目標値をｒとして出力との誤差積分値を算出したものとする。このとき式(3.44)のスライディングモード到達条件を満足するためには、制御入力ｕ(ｔ)が以下の条件を満たさなければならない。
【数９０】
0005445887.tif 000092
【０２４２】
制御入力を次式のような形とする。
【数９１】
0005445887.tif 000093
【０２４３】
ここでｓｑｎ(σ)に関する符号関数である。このとき、超平面上でのシステムの状態、すなわち等価制御系の応答と外乱が既知であれば、式(3.40)を常に満足させることが可能となる。しかし、実際に外乱などの観測が不可能な場合が多く、そのため、切換ゲインｋを適当な定数とするがいっぱいである。切換幅をあまり大きく取りすぎると、制御入力にチャタリングを引き起こす原因となる。また、実際の連続系のプラントに使用するとき、アナログ切換装置には何らかのおそれが存在する。また、ディジタルで使用する場合には、サンプリング周期を無限小にすることができず、一種の遅れが存在する。このような理想的でない装置を使うと、切換周波数を無限に大きくすることは不可能であり、現実のスライディングモードは切換面をすべることにはならず、その近傍でチャタリングすることになる。その制御入力のチャタリングにより制御対象をモデル化する際に無視した高周波領域を励振し、スピルオーバを引き起こす原因となる可能性がある。このようなチャタリングを低減させるため、本研究では次のような平滑化定数を制御入力に導入する。
【数９２】
0005445887.tif 000094
【０２４４】
切換制御入力を平滑化することによりチャタリングの低減が可能となるが、超平面近傍でスライディングモード到達条件を満たすことができず、切換超平面からの偏差が発生する。このとき、式(3.50)は
【数９３】
0005445887.tif 000095
【０２４５】
となる。上式において定常的な外乱などが存在するとき、切換関数σの値がσ→∞となるとき、式(3.50)の右辺→０となり、切換幅ｋを定数とする場合には左辺＝右辺となるσで釣り合ってしまう。これにより、マッチング条件を満たすシステムパラメータ変動および外乱に対し不感であると言うスライディングモード制御の特徴が失われることになるため、平滑化定数δは必要最小限の大きさにするべきである。また、このときのはサンプリングタイムと密接な関係をとっており、離散時間系の設計をする必要があるものと思われる。
【０２４６】
４.モデリング
４．１  ジャイロ効果を考慮しない場合のモデリング
使用する実験装置は図８１の左図で説明したような５軸制御型であるが、ｚ方向の浮上はＰＩＤ制御によって制御されており、残ったラジアル方向の４自由度だけを本研究の制御対象と考えている。右図はロータのＸ、Ｙ軸方向の運動形式が全く同じ対称なので、Ｘ−Ｚ平面だけを示す。式(4.1)はロータの重心に対する運動方程式である。
【数９４】
0005445887.tif 000096
【０２４７】
ここで、Ｆｘｕ、Ｆｘｌ、Ｆｙｕ、ＦｙｌはロータのＸ、Ｙ軸方向の上下部に対向した電磁石間の吸引力差である。それぞれは式(4.2)のようになる。
【数９５】
0005445887.tif 000097
【０２４８】
重心に対する運動方程式から、状態空間モデルを得られるが、状態量は重心位置でのロータ変位と回転角であり、しかし、実験において、センサの観測点は重心の位置と異なって、ロータの上下部の電磁石に近いところに置かれている。その問題を解消するため、モデルの簡単化を行う。まず、図８１の右側の図を見ると、重心位置でのロータ変位と電磁石に対応する点でのロータ変位、さらにそれぞれの回転角と変位間の関係は式(4.3)のようになる。
【数９６】
0005445887.tif 000098
【０２４９】
さらに、すべての電磁石に対してロータ重心に加えている電磁石と回転モーメントをロータ上下部の電磁石により分離して考えると、式(4.3)を式(4.1)に代入して、最後に運動方程式は式(4.4)のように簡単となる。
【数９７】
0005445887.tif 000099
【０２５０】
係数ａｕ、ａｌは式(4.5)となる。
【数９８】
0005445887.tif 000100
【０２５１】
式(4.4)によって状態方程式は以下のように表される。
【数９９】
0005445887.tif 000101
【０２５２】
ここで、ｘ１、ｘ２は変位と速度の状態量、制御入力Ｕは対向した電磁石間の吸引力差である。
【数１００】
0005445887.tif 000102
【０２５３】
それぞれのパラメータを表４．１に示す。
【表２】
0005445887.tif 000103
【０２５４】
５．非線形制御器設計と性能検証１
５．１  ＬＱＩ制御器設計
５．１．１  レギュレータ制御器設計
まず、第３章に述べたバックステッピング制御理論により、状態変位はｘ１、速度はｘ２とおき、システムは式(5.1)のような２次系となる。
【０２５５】
ｘ１＝ｘ２
ｘ２＝Ｂ２Ｕ    (5.1)
ここで、ｘ１＝(ｘｕ ｘｌ ｙｕ ｙｌ)Ｔ、ｘ２＝(ｘｕ ｘｌ ｙｕ ｙｌ)Ｔ
Ｕ＝(Ｆｘｕ  Ｆｘｌ  Ｆｙｕ  Ｆｙｌ)Ｔ
評価関数は式(5.2)となる。
【数１０１】
0005445887.tif 000104
【０２５６】
ここで、
【数１０２】
0005445887.tif 000105
【０２５７】
ように設定する。さらに、目標値入力に定常偏差なく追従するサーボ系の設計を行う。ここで、ｒを一定目標値とするとき、制御問題は閉ループ系を安定し、かつ、出力ｙを目標値に一致させることである。この場合、目標値と出力の誤差ｅ＝ｒ−ｙの積分ｘｒを含めた式(5.3)のような拡大状態方程式を得る。
【数１０３】
0005445887.tif 000106
【０２５８】
ここで、ｘａ＝[ｘＴｘｒｔ]とすると、評価関数は式(5.4)となる。
【数１０４】
0005445887.tif 000107
【０２５９】
評価関数を最小にする状態フィードバック制御を考える。状態フィードバック制御入力は式(5.5)となる。
【数１０５】
0005445887.tif 000108
【０２６０】
ブロック線図を図８２に示す。
【０２６１】
５．１．２  制御電流の切換法
実験するとき、制御入力は対向した両電磁石間の吸引力ではなくて、端的な電磁石に流している電流であり、それで、求められた吸引力の正負により２対の電磁石の切換方法に基づいて電流を逆算する。ここで、ロータのＸ軸方向の電磁石の切換法を詳細に述べる。
【０２６２】
1) Ｆｘｕ≧０のとき
Ｆｘｕ≧０のときは、図８１において、電磁石１が吸引すればよい。そのときの力と電流ｉ１の関係式（ｉ３＝０として）
【数１０６】
0005445887.tif 000109
【０２６３】
から、ｉ１を以下のように示す。
【数１０７】
0005445887.tif 000110
【０２６４】
2) Ｆｘｕ＜０のとき
Ｆｘｕ＜０のときは、図８１において、電磁石３が吸引すればよい。そのときの力と電流ｉ３の関係式（ｉ１＝０として）
【数１０８】
0005445887.tif 000111
【０２６５】
から、ｉ３を以下のように示す。
【数１０９】
0005445887.tif 000112
【０２６６】
下部についても同様に考えることができる。
【０２６７】
1) Ｆｘｌ≧０のとき
Ｆｘｌ≧０のときは、図８１において、電磁石５が吸引すればよい。そのときの力と電流ｉ５の関係式（ｉ７＝０として）
【数１１０】
0005445887.tif 000113
【０２６８】
から、ｉ５を以下のように示す。
【数１１１】
0005445887.tif 000114
【０２６９】
2) Ｆｘｌ＜０のとき
Ｆｘｌ＜０のときは、図８１において、電磁石７が吸引すればよい。そのときの力と電流ｉ７の関係式（ｉ５＝０として）
【数１１２】
0005445887.tif 000115
【０２７０】
から、ｉ７を以下のように示す。
【数１１３】
0005445887.tif 000116
【０２７１】
この切換条件によって吸引力を各電磁石の電流に変換し、電磁石の出力とすることで制御を行う。
【０２７２】
５．１．３  シミュレーション
図８３(a)はタッチダウンから平衡点までの初期値と０．０５秒時に、５０Ｎの外乱を印加したインパルス応答。この結果より、システムが整定時間が早い、ロータが中心に行って、良い安定性がわかった。図８３(b)は制御器設計するとき、仮想的な力制御入力を示す。この図より、ロータが平衡点のとき、対向した電磁石間の電磁石は零になり、その後、一瞬にインパルス外乱を入れると、外乱と等値な電磁力を同時に発生して、影響を打ち消した、良いロバスト性を有することがわかった。図８４は、切換条件に基づいて、力制御入力から実際の制御電流に引き換えることを示す。ここで、ｉ１とｉ３、ｉ５とｉ７はロータの上部に対向している電磁石中の電流。ｉ２とｉ４、ｉ６とｉ８は下部の電磁石電流。０．０５秒にインパルス外乱の影響で電流が急に発生した、再び零になることから、最終にシステムの消費電力が零になることを示す。図８５〜図８８は、回転速度を５０Ｈｚ、１００Ｈｚに設定したときに対応しているそれぞれの制御器のシミュレーション結果。
【０２７３】
５．１．４  実験結果
実験では、用いた制御器はシミュレーションと完全に同じものである。制御器をコンパイルして、サンプリング時間は０．１２５ｍｓとする。ＤＳＰを用いて、４つの変位センサから、磁気軸受の位置情報をＡＤ変換器を通して入手し、切換条件によって８つの制御電流に変換して、アンプに入り増幅された後アクチュエータを駆動して、実験を行う。実験での各装置の関係と行う方法を図８９で説明している。
【０２７４】
実験では、Ｚ軸方向がもうＰＩＤ制御されたと前提して、ロータはタッチダウン状態からの浮上と回転実験を行う。図９０はロータが浮上する状況を示す。上はリサージュ波形であり、真ん中は位相のパワースペクトルである。そして、下はＸ軸方向の制御電流である。この図からみると、上側と下側両方とも中心の一点に収束し、定常偏差がなく、安定に浮上することがわかる。図９１〜図９７は、ロータが回転しているとき、各回転段階での軌跡、変位のパワースペクトルと制御電流である。結果より、上下部のロータ軌跡が中心に行ったり、戻ったり繰り返しているが、８５Ｈｚの回転周波数を超えて、５７００ｒｐｍの回転速度に達した。さらに、図９４と図９５の制御電流から見ると、電流ｉ１の谷はちょうど対向するｉ３の波の頂上である。すなわち、電流がきちんと切り換えしていることがわかる。しかし、ステップ応答が遅く、下部ロータの軌跡が大きくなることについて、ゲインのパラメータを調節し、下部制御入力を増大させるが、５０Ｈｚ付近でジャイロ効果が出始め、安定性が悪くなった。
【０２７５】
５．２  実験結果の分析
図９０〜図９７の結果から、ロータを定常偏差なく、中心の平衡点に収束させることができ、電流もきちんと切り換えしていることがわかる。しかし、回転時には、上下部のロータ軌跡は中心に行ったり戻ったりして、９７Ｈｚでタッチダウンした。図９７の軌跡から見ると、花びらのような軌跡が発生しており、それはジャイロ効果により後向き歳差運動が発生していると考えられる。つまり、回転速度が上がることによってジャイロ効果が強くなって安定性が崩れてしまったことがタッチダウンする原因と思われる。その問題を解決するため、ジャイロ効果を考慮したモデルを作った。
【０２７６】
５．３  ジャイロ効果
磁気軸受を用いるロータは、比較的高速回転のものが多いので、ジャイロ効果が問題となる。ジャイロ効果というのは、高速フライホイールが傾き運動をするときに発生する作用である。この効果により、次の特徴的な性質が発生する。
【０２７７】
・停止時には１つである固有振動数が前回り振れ回りモード（フォワードモード）と後回り振れ回りモード（バックワードモード）に２つに分かれ、ジャイロ効果の大小によって変化幅が決まる。
【０２７８】
・ロータのＸ、Ｙ軸の振動が連成し、振れ回り軌跡が円に近くなる。このジャイロ効果より、ロータ系は磁気軸受にとって回転数に伴って振動特性が変化するとともに、Ｘ、Ｙ方向が連成する制御対象となるので、磁気軸受の制御特性を決める上で大きな障害となる。
【０２７９】
ここで、ジャイロモーメントを詳しく解析する。まず、図９８に示すように、フライホイールがＸ軸（Ｙ軸）に向け撓んで、そして傾いたとき、自らのＺ軸に沿って回っているので、Ｙ軸（Ｘ軸）に関する回転モーメントを生じることである。この効果は、回転速度の増大に従って強くなっていく。結局は、フライホイールは前向き固有振動と後向き固有振動を有することになる。
【０２８０】
新しい座標系によると、角運動量の成分は次式で表される。
【数１１４】
0005445887.tif 000117
【０２８１】
ここに、θｘ、θｙ、θｚはフライホイールの角変位、ＩｒはＸ軸、Ｙ軸に関する慣性モーメント、ＩｚはＺ軸に関する極慣性モーメントである。空間に固定された座標系（ｘ,ｙ,ｚ）に関して運動方程式を立てるので、角運動量の成分を固定座標系へ投影すると、その成分は図９８より
【数１１５】
0005445887.tif 000118
【０２８２】
となる。ただし、角度はすべて微小とした。角運動量の法則による式(5.7)を時間で微分して、さらにθｚ＝−ωを代入すると、Ｘ軸、Ｙ軸に関する回転モーメントは以下のように得られる。
【数１１６】
0005445887.tif 000119
【０２８３】
ここで、−ωＩｚθｙ、ωＩｚθｘの項はいわゆるジャイロ効果の影響を示すものである。
【０２８４】
５．４  ジャイロ効果を考慮した場合のモデリング
ジャイロモーメントを考慮したロータの重心に対する運動方程式は式(5.9)である。
【数１１７】
0005445887.tif 000120
【０２８５】
前節にも述べたようにモデルを簡単化して、最終のモデルは式(5.10)で示す。
【数１１８】
0005445887.tif 000121
【０２８６】
式(5.10)によって状態方程式は以下のように表される。
【数１１９】
0005445887.tif 000122
【０２８７】
ここで、ｘ１、ｘ２は変位と速度の状態量、制御入力Ｕは対向した電磁石間の吸引力差である。
【数１２０】
0005445887.tif 000123
【０２８８】
そして、ジャイロ効果を考慮したモデリングに今のＬＱＩ制御器を適用して、根軌跡は図９９に示す。
【０２８９】
この図から、回転速度が増加することにより閉ループの極はますます虚軸に近づくことがわかる。つまり、閉ループのＡ行列の固有値は回転速度により変わっていくＬＰＶシステム（Linear Parameter Varying system）である。
【０２９０】
６．非線形制御器設計と性能検証２
６．１  ゲインスケジュール型スライディングモード制御系の設計
前章で求めたモデリングはＬＰＶシステムである。このＬＰＶシステムに対してゲインスケジュール型スライディングモード制御を提案する。
【０２９１】
６．１．１  ゲインスケジュール型Ｈ∞制御超平面の設計
次の正準系システムのプラントモデルを考える。
【数１２１】
0005445887.tif 000124
【０２９２】
とくに、制御システムをサーボ系と見なすとすれば、誤差信号は次のように表現できる。
【０２９３】
ｅ＝ｒ−Ｃｘ１    (6.2)
この制御系において、超平面を次式のように定義する。
【０２９４】
φ＝Ｓ(ｘ１)＋ｘ２    (6.3)
ここで、Ｓ(ｘ１)はｘ１の線形オペレータである。定義された切換関数は従来の切換関数と比べると新しい状態変数ｚによるダイナミクスを有する。
【０２９５】
ｚ＝Ｆ(ｗ)ｚ＋Ｇ(ｗ)ｅ    (6.4)
Ｓ(ｘ１,ｚ)＝−Ｈ(ｗ)−Ｌ(ｗ)ｅ    (6.5)
そして、式(6.1)と式(6.5)を組み合わせると、拡大システムは次のように得られる。
【数１２２】
0005445887.tif 000125
【０２９６】
また、式(6.3)から
ｘ２＝Ｈｚ−ＬＣｘ１＋Ｌｒ    (6.8)
を得られる。よって、超平面上ではシステムの状態方程式は次式のように低次元化される。
【数１２３】
0005445887.tif 000126
【０２９７】
時変パラメータを用いた制御系は図１００のとおりである。
【０２９８】
しかし、実システムではジャイロ項はＡ２２に入っているから、そのままでは適用できない。そこで、下のようなフィルタを使うことにする。
【０２９９】
１／(０．００１ｓ＋１）    (6.10)
プラントの状態空間モデルとフィルタを組み合わせると、次の拡大系を得る。
【数１２４】
0005445887.tif 000127
【０３００】
そして、図１００は図１０１のようになる。
【０３０１】
式(6.7)の超平面に対して、リアプノフ関数Ｖを次式のように定義する。
【０３０２】
Ｖ＝φＴφ／２    (6.12)
閉ループ系を安定するために次式を満足しなければならない。
【０３０３】
Ｖ＝φＴφ＜０    (6.13)
超平面上に状態が拘束されているとすると
φ＝φ＝０    (6.14)
上の式より、制御入力ｕｅｑは次のように与えられる。
【数１２５】
0005445887.tif 000128
【０３０４】
本研究では、Ｌ(ω)は０にする。さらに、式(6.3)により
【数１２６】
0005445887.tif 000129
【０３０５】
式(6.16)で与えられる安定条件を満足する制御入力を次のように選ぶ。
【数１２７】
0005445887.tif 000130
【０３０６】
ここで、非線形入力は次のように選ぶ。
【数１２８】
0005445887.tif 000131
【０３０７】
６．１．２  コントローラの計算
式(6.15)の中のＨ(ω)、Ｆ(ω)、Ｇ(ω)、Ｌ(ω)つまり制御器をＭＡＴＬＡＢのＬＭＩコントロールツールボックスを使用して設計した。ＬＰＶシステムのボード線図を図１０２に示す。ここで、一番目の入力対４つの出力のボード線図しか載せていないが、その他はほぼ同じため、省略する。破線は変動パラメータが最小値のプラントのボード線図、実線は変動パラメータが最大値のボード線図。
【０３０８】
変動パラメータが最大値のプラントのボード線図から２００Ｈｚ付近にピークがあるので、それに応じる重み関数とコントローラのボード線図を図１０３と図１０４に表す。
【０３０９】
６．２  シミュレーション結果
６．２．１  台形近似法によるリアルタイム離散系
制御器の離散化を行うため、連続時間での制御器のダイナミクスを離散時間でのダイナミクスで近似する必要がある。前節に設計した制御器について、サンプリング時間をＴとし、時刻ｋＴにおける状態をｘ(ｋＴ)＝ｘｋと表すことにする。また、時刻[ｋＴ,(ｋ＋１)Ｔ]において、スケジューリングパラメータα(ｔ)および観測量ｙ(ｔ)は、時刻ｋＴでの値で近似できると仮定する。すなわち、ｋＴ≦ｔ＜(ｋ＋１)Ｔにおいて
【数１２９】
0005445887.tif 000132
【０３１０】
とできると仮定する。このときに、次の台形近似の定理が成り立つ。
【数１３０】
0005445887.tif 000133
【０３１１】
この定理を証明すると、仮定より、時刻(ｋ＋１)Ｔにおける状態ｘｋ＋１について
【数１３１】
0005445887.tif 000134
【０３１２】
が成り立つ。この式を台形近似すると、
【数１３２】
0005445887.tif 000135
【０３１３】
となる。さらに、時刻ｋＴの項と(ｋ＋１)Ｔの項に分けると、
【数１３３】
0005445887.tif 000136
【０３１４】
となる。ここで、
【数１３４】
0005445887.tif 000137
【０３１５】
とおくと、式(6.24)は
【数１３５】
0005445887.tif 000138
【０３１６】
と変形される。さて、今
【数１３６】
0005445887.tif 000139
【０３１７】
が成立する。したがって、
【数１３７】
0005445887.tif 000140
【０３１８】
が存在すれば、式(6.26)は
【数１３８】
0005445887.tif 000141
【０３１９】
ここで、式(6.25)を適用すると、
【数１３９】
0005445887.tif 000142
【０３２０】
となる。これより、
【数１４０】
0005445887.tif 000143
【０３２１】
が導かれる。
【０３２２】
また、式(6.21)も
【数１４１】
0005445887.tif 000144
【０３２３】
が存在すれば、式(6.25)の関係を用いて、同様に導くことができる。図１０５は、タッチダウンから平衡点までの初期値応答と０．０５秒に５０Ｎ外乱を印加したときのインパルス応答を示す。さらに、０．１秒にバックステッピング制御器と同じようにロータ上部のＸ軸方向に０．０２ｍｍ、Ｙ軸方向に０．０１ｍｍの目標値を入れたステップ応答を示す。結果より、ＬＱＩ制御器とほぼ同じ性能を有し、さらに、ステップ応答は速い。図１０６は、４個の切換関数は零になり、状態量が超平面に拘束されて、平衡点に拘束されたことを示す。外乱があるとき、切換関数が上がることが明らかにした。図１０７は仮想的な力制御入力であり、図１０８および図１０９はその力制御入力から切換法によって逆算した制御電流を示す。０．０５秒と０．１秒の電流ピークはそれぞれのインパルス外乱とステップ信号で引き起こされている。
【０３２４】
６．３  実験結果
さらに計算を簡単にするため、式(6.30)は
【数１４２】
0005445887.tif 000145
【０３２５】
に変化する。しかし、次数が高いとサンプリング時間が短いため、ＤＳＰの演算が間に合わなくなって実験できなくなった。その問題に対して、ＭＡＴＬＡＢで離散化してから凸補間するという方法を提案する。
【０３２６】
前のシミュレーションと同じように図１１０はタッチダウンから平衡点までの初期値応答と０．０５秒に５０Ｎ外乱を印加したときのインパルス応答を示す。さらに０．１秒にバックステッピング制御器と同じようにロータ上部のＸ軸方向に０．０２ｍｍ、Ｙ軸方向に０．０１ｍｍの目標値を入れたステップ応答を示す。それと図１０５を比べて見たら、やはり離散化してから凸補間する方法の方に性能が若干悪くなってとくにステップ応答が遅れた。しかし、まだ許容範囲であると思われる。実験においては、制御器はシミュレーションの場合と完全に同じものを用いている。実験では、Ｚ軸方向がもうＰＩＤ制御されたと前提して、ロータはタッチダウン状態からの浮上と回転実験を行う。図１１１はロータが浮上する状況を示す。上はリサージュ波形であり、真ん中は位相のパワースペクトルである。そして、下はＸ軸方向の制御電流である。この図からみると、上側と下側両方とも中心の一点に収束し、定常偏差がなく、安定に浮上することがわかる。図１１２〜図１２０は、ロータが回転しているとき、各回転段階での軌跡、変位のパワースペクトルと制御電流である。
【０３２７】
６．４  実験結果の分析
確かに三角波のような形で、正負にスイッチングをしている。右側図は力制御入力から逆算した制御電流結果より、どんな回転周波数にしても、ＬＱＩ制御器に比べると、電流もはっきりオン・オフしているし、円軌跡もきれいになっていることがわかる。さらに、図１２１に示す今回の根軌跡と図９９と比べてみると、固定コントローラよりゲインスケジュール制御の方が根の動きが小さいことが見える。
【０３２８】
以上の結果より、ゲインスケジュール型スライディングモード制御の安定性はＬＱＩ制御器より良くなっていることがわかった。けれども、ＬＱＩ制御と同じように９０Ｈｚから安定性が悪くなって、１００Ｈｚ付近ではタッチダウンした。前に残していた１００Ｈｚを超えない問題を解決していなかった。
【０３２９】
７．弾性モデルとＰＩＤ制御
７．１  弾性モデリングの導入
剛性モデルでは、ジャイロ効果を無視した補償器を用いても、考慮した補償器を用いても、閉ループ系はともに安定であった。しかし、回転実験では、１００Ｈｚ付近で不安定化している。この不一致の問題に対して、弾性曲げモードも考慮して考察する。弾性モデルの導出は有限要素法によって行う。フライホイールを集中質量として、ロータを図１２２のように８要素に分割する。要素２と３、要素７と８の間に電磁石が配置されている。ここで、減衰と不釣合い振動は考慮しないものとする。
【０３３０】
図１２２のモデリングについて、ジャイロ効果を考慮したロータの運動方程式は次式となる。
【数１４３】
0005445887.tif 000146
【０３３１】
Ｍは質量行列、Ｇはジャイロ反対称行列、Ｋは剛性行列、Ｕは制御入力ベクトル、Ｔは電磁石の設置場所を表す行列（Ｍ,Ｇ,Ｋ∈Ｒ３６×３６,Ｔ∈Ｒ３６×４）、ｑは状態ベクトルである（ｑ∈Ｒ３６×１）。
【０３３２】
まず、フリーフリーロータに対して、固有振動数と軸回転数の関係を図１２３で示す。曲げ一次モードの前向きと後向きは回転後大きな分岐をすることがわかった。このことから、曲げ一次弾性モードの後向き固有値は、ジャイロ効果の影響を強く受けて１００Ｈｚ付近まで下がってきている。この結果、運搬線と後向きモードの交差周波数は１００Ｈｚ付近で、ロータが１００Ｈｚを越えて回転すると、後向き曲げ一次共振点に接近していることがわかった。
【０３３３】
７．２  ＰＩＤ制御
求めた弾性モデリングは７２次であり、モデルベースの制御を実現できない。そのため、非モデルベースのＰＩＤ制御を提案する。ブロック線図は図１２４に示す。
【０３３４】
ＰＩＤ１はシステムを安定することを図るコントローラである。すなわち、ＰＩＤ１の制御でタッチダウン状態から浮上する状態になる。ジャイロ効果を考慮すれば、Ｘ、Ｙ方向がお互いに干渉する。ＰＩＤ２はＸ、Ｙ方向の信号を交差してフィードバックするというクロスフィードバックコントローラである。つまり、システムが安定する上ジャイロ効果を補償するためにお互いに干渉する信号をクロスフィードバックさせている。用いたＰＩＤ制御器のボード線図を図１２５に示す。
【０３３５】
７．２．１  ＰＩＤ制御の実験結果
図１２６はロータが９０Ｈｚで回転している状況を示す。上はリサージュ波形であり、真ん中は位相のパワースペクトルである。そして、下はＸ軸方向の制御電流である。この図から見ると、上側と下側両方とも中心の一点に収束し、定常偏差がなく、安定に浮上することがわかる。図１２７〜図１３３は、ロータが回転しているとき、各回転段階での軌跡、変位のパワースペクトルと制御電流である。
【０３３６】
７．２．２  実験結果の分析
ＰＩＤ制御の結果を見ると、固定コントローラであるため円軌跡と電流の形がどっちでもきれいではないけれども、回転速度が前より３０Ｈｚくらい上がった。クロスフィードバックコントローラがある程度に弾性後向き一次モードの影響を補償していることがわかった。１３０Ｈｚまで増速できるＰＩＤコントローラにより、低次元化した弾性モデルのシミュレーションと実験の比較を図１３４に示す。実線が回転数が０Ｈｚに設定した弾性モデルのシミュレーションであり、破線がＰＩＤ制御器で浮上時実験データである。両方とも、０．０３７秒に０．１＊１０−４ｍのステップ信号を入れてその応答を見る。
【０３３７】
さらに、右下の図を拡大して図１３５に示す。
【０３３８】
図１３５が示した高周波数の振動が弾性モードと思われ、数えてみると１７０Ｈｚくらいになっており、図１２３の一番上の弾性曲げ一次モードの０Ｈｚのときの初期値とほぼ一致しているが、ステップ応答の結果からシミュレーションと実験結果に一致しないことが発生しており、これについても考察を行う。図１３６は今の入力を０．２倍にしたシミュレーション結果である。
【０３３９】
この図から見ると、制御入力を０．２倍にすると実験結果がほぼ一致することがわかる。さらに、全入力と制御入力０．２倍にする閉ループの０Ｈｚから１５０Ｈｚまでの根軌跡を図１３７と図１３８に示す。
【０３４０】
２つの図を比較して見ると、全入力の方に剛性一次モードは６５Ｈｚであり、剛性二次モードは８０Ｈｚである。弾性一次曲げ一次モードに対応する根は虚軸に近づくが、全然安定の領域である。それは実験結果と一致しない所である。逆に０．２倍入力の場合は、剛性一次モードが２５Ｈｚくらいであり、剛性二次モードが見えない。それと図１３５の実験データを一致する。さらに、１４０Ｈｚになると、弾性一次曲げモードに対応根が虚軸に近すぎで、安定性が悪くなる。それも実験の方に一致する。しかし、なぜ働いている入力が実際の０．２倍しかないのかがまだ考察中である。おそらく、ハードウェアとか電気系とかの問題と思われる。
【０３４１】
８．まとめ
本論文では、磁気軸受支持型電力貯蔵フライホイールの制御に関する研究を行った。大きなフライホイールが付けられているので、きわめて不安定なシステムである。従来の方法では、ロータが平衡点に達しても、バイアス電流を流さなければならないので、消費電流が大きくなる。大変位のとき、制御性能が悪くなるという欠点がある。この問題に対して、本研究はバイアス電流を使わないで、非線形を有したまま、現代制御の非線形理論に基づいて、制御器を設計、有効性を検証した。
【０３４２】
本研究においては、まず、線形２次形式評価規範によるサーボ付き最適レギュレータ制御器を設計して、シミュレーションから、良い結果が得られたが、実験の結果からロータを定常偏差なく、中心の平衡点に収束させることができ、電流もきちんと切り換えしていることがわかる。しかし、回転時には、上下部のロータ軌跡は中心に行ったり戻ったりして、９３Ｈｚでタッチダウンした。図９７の軌跡から見ると、花びらのような軌跡が出てきた。それはジャイロ効果による後向き歳差運動軌跡と考えられる。つまり、回転速度が上がることによってジャイロ効果が強くなって安定性が崩れてしまったことがタッチダウンする原因と思われる。その問題を解決するため、ジャイロ効果を考慮したモデルを作った。そして、ジャイロ効果を考慮したモデリングに今のＬＱＩ制御器を適用して、図９９に示した根軌跡から回転速度が増加することにより閉ル−プの極はますます虚軸に近づくことがわかる。つまり、閉ループのＡ行列の固有値は回転速度により変わっていくＬＰＶシステム（Linear Parameter Varying system）であり、固定コントローラも補償できないことがわかる。このＬＰＶシステムに対してゲインスケジュール型スライディングモード制御を提案する。スライディングモードの超平面がスケジューリングして状態量を常に変動している超平面に収束させて安定性を図る方法である。その実験結果とＬＱＩ制御を比べてみると、きれいな軌跡となっており、非同期成分が消えている。電流信号も規則的な信号となっている。ジャイロ効果を補償すると安定性がアップしたことがわかる。しかし、ＬＱＩと同じように１００Ｈｚ付近では安定性が悪くなって、タッチダウンしてしまった。本システムは１００Ｈｚのところに弾性モードの共振点が出てきた。弾性モードを考慮していないＬＱＩコントローラおよびゲインスケジュール型スライディングモードコントローラがその共振点にぶつかって安定性が崩れて、タッチダウンしてしまう。弾性モデリングを実行したが、次数が高いため、モデルベースの制御を実現できない。そこで、非モデルベースのＰＩＤ制御で１３０Ｈｚまで増速できたが、ＰＩＤコントローラの閉ループ系と、低次元化した弾性モデルのシミュレーションと実験の結果が一致しないという問題がある。制御入力を元の０．２倍にすると実験とよく一致したが、なぜこのようにすると合致するのか検討中である。ハードウェアおよび電気系の問題と思われる。今後の課題としては以下のようになる。
【０３４３】
１．ＰＩＤコントローラをもっとチューニングする。上下のバランスを良くして、さらに周波数によって切り換えして２００Ｈｚまで増速する。
【０３４４】
２．実験結果とシミュレーションの一致しない原因を調べる。低次化した弾性モデルを修正して、ゲインスケジュール型スライディングモード制御系を設計する。
【０３４５】
磁気軸受装置の全体構成あるいは各部の構成、制御の方法などは、上記実施形態のものに限らず、適宜変更可能である。
【図面の簡単な説明】
【０３４６】
【図１】図１は、この発明の実施形態における磁気軸受装置の全体構成を概略的に示すブロック図である。
【図２】図２は、図１の磁気軸受装置の機械的部分の主要部を示す縦断面図である。
【図３】図３は、図２の主要部の斜視図である。
【図４】図４は、一般化プラントを示すブロック図である。
【図５】図５は、Ｈ∞制御のブロック図である。
【図６】図６は、Ｈ∞制御のブロック図である。
【図７】図７は、閉ループ系のブロック図である。
【図８】図８は、ロバスト安定化問題を表すフィードバック系のブロック図である。
【図９】図９は、混合感度問題を表すフィードバック系のブロック図である。
【図１０】図１０は、クローズドループのブロック図である。
【図１１】図１１は、固有振動線図である。
【図１２】図１２は、フライホイール磁気軸受のモデルを示す斜視図である。
【図１３】図１３は、図１２のＸ−Ｚ平面における図である。
【図１４】図１４は、一般化プラントを示すブロック図である。
【図１５】図１５は、上部磁気軸受を示す説明図である。
【図１６】図１６は電磁石の吸引力を示すグラフである。
【図１７】図１７は、閉ループシステムのブロック図である。
【図１８】図１８は、初期応答を示す図である。
【図１９】図１９は、制御電流を示す図である。
【図２０】図２０は、実験装置を示すブロック図である。
【図２１】図２１は、実験装置のモニタを示す図である。
【図２２】図２２は、電磁石の制御電流を示す図である。
【図２３】図２３は、実験結果を示す図である。
【図２４】図２４は、入力電流を示す図である。
【図２５】図２３は、実験結果を示す図である。
【図２６】図２６は、入力電流を示す図である。
【図２７】図２７は、実験結果を示す図である。
【図２８】図２８は、入力電流を示す図である。
【図２９】図２９は、実験結果を示す図である。
【図３０】図３０は、入力電流を示す図である。
【図３１】図３１は、実験結果を示す図である。
【図３２】図３２は、入力電流を示す図である。
【図３３】図３３は、実験結果を示す図である。
【図３４】図３４は、入力電流を示す図である。
【図３５】図３５は、実験結果を示す図である。
【図３６】図３６は、入力電流を示す図である。
【図３７】図３７は、実験結果を示す図である。
【図３８】図３８は、入力電流を示す図である。
【図３９】図３９は、実験結果を示す図である。
【図４０】図４０は、入力電流を示す図である。
【図４１】図４１は、実験結果を示す図である。
【図４２】図４２は、入力電流を示す図である。
【図４３】図４３は、重み関数を示す図である。
【図４４】図４４は、ＬＴＩ制御器の周波数応答を示す図である。
【図４５】図４５は、初期値応答を示す図である。
【図４６】図４６は、入力電流を示す図である。
【図４７】図４７は、初期値応答を示す図である。
【図４８】図４８は、入力電流を示す図である。
【図４９】図４９は、初期値応答を示す図である。
【図５０】図５０は、入力電流を示す図である。
【図５１】図５１は、初期値応答を示す図である。
【図５２】図５２は、入力電流を示す図である。
【図５３】図５３は、閉ループの根軌跡の図である。
【図５４】図５４は、実験結果を示す図である。
【図５５】図５５は、入力電流を示す図である。
【図５６】図５６は、実験結果を示す図である。
【図５７】図５７は、入力電流を示す図である。
【図５８】図５８は、実験結果を示す図である。
【図５９】図５９は、入力電流を示す図である。
【図６０】図６０は、実験結果を示す図である。
【図６１】図６１は、入力電流を示す図である。
【図６２】図６２は、実験結果を示す図である。
【図６３】図６３は、入力電流を示す図である。
【図６４】図６４は、実験結果を示す図である。
【図６５】図６５は、入力電流を示す図である。
【図６６】図６６は、実験結果を示す図である。
【図６７】図６７は、入力電流を示す図である。
【図６８】図６８は、実験結果を示す図である。
【図６９】図６９は、入力電流を示す図である。
【図７０】図７０は、実験結果を示す図である。
【図７１】図７１は、入力電流を示す図である。
【図７２】図７２は、磁気軸受装置のロータを８要素に分割した図である。
【図７３】図７３は、固有振動数と軸回転数の関係を示す図である。
【図７４】図７４は、ボード線図である。
【図７５】図７５は、低次元化モデルの初期値応答を示す図である。
【図７６】図７６は、閉ループシステムのブロック図である。
【図７７】図７７は、低次元化モデルの初期値応答を示す図である。
【図７８】図７８は、閉ループの固有値解析を示す図である。
【図７９】図７９は、低次元化モデルの初期値応答を示す図である。
【図８０】図８０は、安定性の概念を示す図である。
【図８１】図８１は、５軸制御型磁気軸受装置のモデルを示す図である。
【図８２】図８２は、ＬＱＩ制御系のブロック図である。
【図８３】図８３は、初期応答を示す図である。
【図８４】図８４は、電磁石の制御電流を示す図である。
【図８５】図８５は、初期応答を示す図である。
【図８６】図８６は、電磁石の制御電流を示す図である。
【図８７】図８８は、初期応答を示す図である。
【図８８】図８８は、電磁石の制御電流を示す図である。
【図８９】図８９は、実験装置のブロック図である。
【図９０】図９０は、実験結果を示す図である。
【図９１】図９１は、実験結果を示す図である。
【図９２】図９２は、実験結果を示す図である。
【図９３】図９３は、実験結果を示す図である。
【図９４】図９４は、実験結果を示す図である。
【図９５】図９５は、実験結果を示す図である。
【図９６】図９６は、実験結果を示す図である。
【図９７】図９７は、実験結果を示す図である。
【図９８】図９８は、ジャイロモーメントの説明図である。
【図９９】図９９は、閉ループの根軌跡を示す図である。
【図１００】図１００は、ＬＰＶ低次元化システムのブロック図である。
【図１０１】図１０１は、フィルタ付きＬＰＶ低次元化システムのブロック図である。
【図１０２】図１０２は、ＬＰＶシステムのボード線図である。
【図１０３】図１０３は、重み関数の図である。
【図１０４】図１０４は、コントローラのボード線図を示す図である。
【図１０５】図１０５は、初期値応答、インパルス応答、ステップ応答を示す図である。
【図１０６】図１０６は、スイッチング関数の図である。
【図１０７】図１０７は、力制御入力の図である。
【図１０８】図１０８は、Ｘ軸の制御電流を示す図である。
【図１０９】図１０９は、Ｙ軸の制御電流を示す図である。
【図１１０】図１１０は、初期値応答、インパルス応答、ステップ応答を示す図である。
【図１１１】図１１１は、実験結果を示す図である。
【図１１２】図１１２は、実験結果を示す図である。
【図１１３】図１１３は、実験結果を示す図である。
【図１１４】図１１４は、実験結果を示す図である。
【図１１５】図１１５は、実験結果を示す図である。
【図１１６】図１１６は、実験結果を示す図である。
【図１１７】図１１７は、実験結果を示す図である。
【図１１８】図１１８は、実験結果を示す図である。
【図１１９】図１１９は、実験結果を示す図である。
【図１２０】図１２０は、実験結果を示す図である。
【図１２１】図１２１は、閉ループの根軌跡を示す図である。
【図１２２】図１２２は、磁気軸受装置のロータを８要素に分割した図である。
【図１２３】図１２３は、固有振動数と軸回転数の関係を示す図である。
【図１２４】図１２４は、ＰＩＤコントローラの閉ループのブロック図である。
【図１２５】図１２５は、ＰＩＤコントローラのボード線図である。
【図１２６】図１２６は、実験結果を示す図である。
【図１２７】図１２７は、実験結果を示す図である。
【図１２８】図１２８は、実験結果を示す図である。
【図１２９】図１２９は、実験結果を示す図である。
【図１３０】図１３０は、実験結果を示す図である。
【図１３１】図１３１は、実験結果を示す図である。
【図１３２】図１３２は、実験結果を示す図である。
【図１３３】図１３３は、実験結果を示す図である。
【図１３４】図１３４は、シミュレーションと実験の比較結果を示す図である。
【図１３５】図１３５は、図１３４の右下の図の拡大図である。
【図１３６】図１３６は、シミュレーションと実験の比較結果を示す図である。
【図１３７】図１３７は、閉ループの根軌跡を示す図である。
【図１３８】図１３８は、閉ループの根軌跡を示す図である。
【図１３９】図１３９は、テーブルを用いたＧＳ制御を説明する図である。
【符号の説明】
【０３４７】
(1)     機械本体
(2)     コントローラ
(4)     回転体
(4a)    フライホイール
(5)     アキシアル磁気軸受
(6)(7)    ラジアル磁気軸受
(8)     変位検出部
(9)     電動モータ
(10)    回転センサ
(11)(12)    保護軸受
(25)    アキシアル変位センサ
(26)(27)    ラジアル変位センサユニット
(28a)(28b)    アキシアル電磁石
(29a)(29b)(29c)(29d)(30a)(30b)(30c)(30d)    ラジアル電磁石
(31a)(31b)(31c)(31d)(32a)(32b)(32c)(32d)    ラジアル変位センサ
【図１】
0005445887.tif 000149
【図２】
0005445887.tif 000150
【図３】
0005445887.tif 000151
【図４】
0005445887.tif 000152
【図５】
0005445887.tif 000153
【図６】
0005445887.tif 000154
【図７】
0005445887.tif 000155
【図８】
0005445887.tif 000156
【図９】
0005445887.tif 000157
【図１０】
0005445887.tif 000158
【図１１】
0005445887.tif 000159
【図１２】
0005445887.tif 000160
【図１３】
0005445887.tif 000161
【図１４】
0005445887.tif 000162
【図１５】
0005445887.tif 000163
【図１６】
0005445887.tif 000164
【図１７】
0005445887.tif 000165
【図１８】
0005445887.tif 000166
【図１９】
0005445887.tif 000167
【図２０】
0005445887.tif 000168
【図２１】
0005445887.tif 000169
【図２２】
0005445887.tif 000170
【図２３】
0005445887.tif 000171
【図２４】
0005445887.tif 000172
【図２５】
0005445887.tif 000173
【図２６】
0005445887.tif 000174
【図２７】
0005445887.tif 000175
【図２８】
0005445887.tif 000176
【図２９】
0005445887.tif 000177
【図３０】
0005445887.tif 000178
【図３１】
0005445887.tif 000179
【図３２】
0005445887.tif 000180
【図３３】
0005445887.tif 000181
【図３４】
0005445887.tif 000182
【図３５】
0005445887.tif 000183
【図３６】
0005445887.tif 000184
【図３７】
0005445887.tif 000185
【図３８】
0005445887.tif 000186
【図３９】
0005445887.tif 000187
【図４０】
0005445887.tif 000188
【図４１】
0005445887.tif 000189
【図４２】
0005445887.tif 000190
【図４３】
0005445887.tif 000191
【図４４】
0005445887.tif 000192
【図４５】
0005445887.tif 000193
【図４６】
0005445887.tif 000194
【図４７】
0005445887.tif 000195
【図４８】
0005445887.tif 000196
【図４９】
0005445887.tif 000197
【図５０】
0005445887.tif 000198
【図５１】
0005445887.tif 000199
【図５２】
0005445887.tif 000200
【図５３】
0005445887.tif 000201
【図５４】
0005445887.tif 000202
【図５５】
0005445887.tif 000203
【図５６】
0005445887.tif 000204
【図５７】
0005445887.tif 000205
【図５８】
0005445887.tif 000206
【図５９】
0005445887.tif 000207
【図６０】
0005445887.tif 000208
【図６１】
0005445887.tif 000209
【図６２】
0005445887.tif 000210
【図６３】
0005445887.tif 000211
【図６４】
0005445887.tif 000212
【図６５】
0005445887.tif 000213
【図６６】
0005445887.tif 000214
【図６７】
0005445887.tif 000215
【図６８】
0005445887.tif 000216
【図６９】
0005445887.tif 000217
【図７０】
0005445887.tif 000218
【図７１】
0005445887.tif 000219
【図７２】
0005445887.tif 000220
【図７３】
0005445887.tif 000221
【図７４】
0005445887.tif 000222
【図７５】
0005445887.tif 000223
【図７６】
0005445887.tif 000224
【図７７】
0005445887.tif 000225
【図７８】
0005445887.tif 000226
【図７９】
0005445887.tif 000227
【図８０】
0005445887.tif 000228
【図８１】
0005445887.tif 000229
【図８２】
0005445887.tif 000230
【図８３】
0005445887.tif 000231
【図８４】
0005445887.tif 000232
【図８５】
0005445887.tif 000233
【図８６】
0005445887.tif 000234
【図８７】
0005445887.tif 000235
【図８８】
0005445887.tif 000236
【図８９】
0005445887.tif 000237
【図９０】
0005445887.tif 000238
【図９１】
0005445887.tif 000239
【図９２】
0005445887.tif 000240
【図９３】
0005445887.tif 000241
【図９４】
0005445887.tif 000242
【図９５】
0005445887.tif 000243
【図９６】
0005445887.tif 000244
【図９７】
0005445887.tif 000245
【図９８】
0005445887.tif 000246
【図９９】
0005445887.tif 000247
【図１００】
0005445887.tif 000248
【図１０１】
0005445887.tif 000249
【図１０２】
0005445887.tif 000250
【図１０３】
0005445887.tif 000251
【図１０４】
0005445887.tif 000252
【図１０５】
0005445887.tif 000253
【図１０６】
0005445887.tif 000254
【図１０７】
0005445887.tif 000255
【図１０８】
0005445887.tif 000256
【図１０９】
0005445887.tif 000257
【図１１０】
0005445887.tif 000258
【図１１１】
0005445887.tif 000259
【図１１２】
0005445887.tif 000260
【図１１３】
0005445887.tif 000261
【図１１４】
0005445887.tif 000262
【図１１５】
0005445887.tif 000263
【図１１６】
0005445887.tif 000264
【図１１７】
0005445887.tif 000265
【図１１８】
0005445887.tif 000266
【図１１９】
0005445887.tif 000267
【図１２０】
0005445887.tif 000268
【図１２１】
0005445887.tif 000269
【図１２２】
0005445887.tif 000270
【図１２３】
0005445887.tif 000271
【図１２４】
0005445887.tif 000272
【図１２５】
0005445887.tif 000273
【図１２６】
0005445887.tif 000274
【図１２７】
0005445887.tif 000275
【図１２８】
0005445887.tif 000276
【図１２９】
0005445887.tif 000277
【図１３０】
0005445887.tif 000278
【図１３１】
0005445887.tif 000279
【図１３２】
0005445887.tif 000280
【図１３３】
0005445887.tif 000281
【図１３４】
0005445887.tif 000282
【図１３５】
0005445887.tif 000283
【図１３６】
0005445887.tif 000284
【図１３７】
0005445887.tif 000285
【図１３８】
0005445887.tif 000286
【図１３９】
0005445887.tif 000287