(57)【要約】【課題】  風洞試験に用いる磁力支持天秤装置のような磁気支持装置における少なくとも位置・姿勢角についてフィードバック制御を行う制御系において、多数に上る制御定数の調整の容易化を図ることができる磁力支持装置における制御定数決定方法を提供する。【解決手段】  供試体の少なくとも位置・姿勢角についてのフィードバック制御をする制御装系において、外乱入力から位置・姿勢角の出力までの伝達関数に関連する感度関数と、また外乱入力からコイル・電磁石に印加される電流の出力までの伝達関数に関連する相補感度関数とが定められる。感度関数を、最大値を抑制し且つ振動周波数域を高周波数域にシフトさせるように整形することで、外乱に対して供試体の振動を抑えることができる。相補感度関数は、それに応じて高周波数域にシフトする整形がされる。伝達関数が有する多数の制御定数は、両関数に基づいて、ロバスト制御で定義されるＨ∞ノルムを最小化する値として決定される。【選択図】  図１
2005157816.tif 000002
【特許請求の範囲】
【請求項１】
金属及び磁石を有する供試体の位置・姿勢角を検出する位置センサ、前記供試体に磁気力を作用させるための１つ又は複数個のコイル・電磁石、及び前記位置センサにより検出された前記供試体の位置・姿勢角に基づいて前記供試体の位置・姿勢角についてのフィードバック制御をするため前記コイル・電磁石に印加する電流を演算・出力する制御装置を有し、前記供試体を前記磁気力により空中の任意の位置・姿勢角に保持する磁力支持装置において、外乱入力から前記位置・姿勢角の出力までの伝達関数に関連して定められる感度関数と、前記外乱入力から前記コイル・電磁石に印加される前記電流の出力までの伝達関数に関連して定められる相補感度関数とに基づいて、前記伝達関数が有する制御定数を、ロバスト制御で定義されるＨ∞ノルムを最小化する値に定めることから成る磁力支持装置における制御定数決定方法。
【請求項２】
前記感度関数の最大値を抑えると共に高周波数側にシフトさせる整形を施す重み関数と、前記相補感度関数を高周波数側にシフトさせる整形を施す重み関数をそれぞれ定め、前記各重み関数によって整形された前記感度関数と前記相補感度関数とに対して前記Ｈ∞ノルムを定めることから成る請求項１に記載の磁力支持装置における制御定数決定方法。
【請求項３】
前記制御装置は、前記コイル・電磁石により形成される磁場を検出する磁気センサと前記磁気センサが検出した前記磁場に基づいて前記コイル・電磁石を制御する磁場制御器とを備え、前記供試体の前記位置・姿勢角についてのフィードバック制御と同時に、前記磁気センサが検出した前記磁場に基づいて前記磁場についてのフィードバック制御を行っており、前記磁場についてのフィードバック制御を行うときの前記感度関数及び前記相補感度関数は、前記位置・姿勢角についてのフィードバック制御のみを行うときの前記磁場の伝達関数を前記磁場についてのフィードバック制御系全体の伝達関数に置換することで定められることから成る請求項１又は２に記載の磁力支持装置における制御定数決定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、金属及び磁石を備える供試体を磁気力により空中の任意の位置と姿勢角に制御可能に支持する磁力支持装置に関する。
【背景技術】
【０００２】
航空機やスペースシャトルのような宇宙往復機の開発・設計に際しては、気流が機体に及ぼす影響が時間を追う毎に変化する動的特性を知ることが不可欠である。動的特性が求められる場合としては、例えば、突風や気流の乱れ等のように気流の流れが変化する場合、気流の密度や温度等の特性自体が変化する場合、或いは操舵によって気流に対して航空機の姿勢が変更される場合等がある。このような状況の変化に応じて模型に作用する揚力、抗力、モーメント等の動的に変化する力を知ることで、実機の位置・姿勢角に応じた動的特性を知ることができ、そうした実機の動的特性は、航空機の翼や胴体等の外形形状、構造強度、エンジン特性、或いは操舵装置等の開発・設計に資することができる。
【０００３】
従来、物体の空力的な特性を模型で試験する風洞試験においては、風洞設備の測定部での模型を支持体で支持することが一般的に行われてきたが、この模型の支持形態は支持体自体が模型の周りにおける空気流れに影響を及ぼす。したがって、風洞試験結果をそのまま模型（ひいては物体）の空力特性として採用することができない。そこで、風洞試験において、模型を磁力で支持することが提案されている。模型を磁力支持することによって支持体が不要となるので、支持体が存在することに起因した模型への空力的な影響が取り除かれ、理想的な風洞試験を行うことができる。
【０００４】
風洞模型のための磁力支持装置は、風洞試験において模型の周りを流れる気流が模型に作用する揚力、抗力、ピッチング（縦揺れ）モーメント等の静的又は動的な空力特性を、模型の内部に設けられる磁石と相互作用する磁気力を生じさせるために設けられているコイルに流す電流の大きさに置き換えて測定する装置である。こうした空気力とコイル電流の大きさとの関係を調べて予めマップ、関数、表等の対応関係を用意しておくことにより、コイル電流を測定することで模型に作用する空力特性を知ることができる。
【０００５】
図４を参照して、磁力支持型風洞及び磁力支持装置の概要を説明する。図４は磁力支持装置を備えた風洞、及び磁力支持装置に用いられる計測系と制御系の概要を示す斜視図である。図４に示す磁力支持装置２０は、アルミなどの非磁性材料で所望の形状に成形されている風洞模型（以下、単に「模型」という）７を磁気の力で気流中に浮揚支持する装置である。模型７にはその内部に、磁化された物質、超伝導コイルのような電流を流し続けているコイル、或いは永久磁石等から成る強力な円柱状の磁石体が搭載される。模型７の磁石体には、風洞の測定部の周りに配置したコイルに通電することにより生じた磁場との磁気作用によって磁気力が生じる。この磁気力によって、模型７を重力や空気力に抗して浮上支持させて、支持体による干渉がない風洞試験を実現することができる。
【０００６】
模型に力を作用させるための磁場は、磁気支持コイルとしてのコイル２２〜２５とコイル２６〜２９とから成る二つの磁気回路、及びその外側に配置された同じく磁気支持コイルとしての空芯コイル２１，３０によって発生される。磁気回路の各コイル２２〜２９に流れる電流を調節することにより、磁気回路内のｙ−ｚ面内での磁場の強さと方向及びそれらのｘ軸方向の変化率を連続的に変化させることができる。また、空芯コイル２１，３０に流れる電流を調節することによりｘ軸方向磁場の強さのｘ軸方向で見た変化率を制御することができ、都合５軸の制御が可能である。即ち、磁気回路のコイル２２〜２９は、模型７に働く揚力と縦揺れモーメントとに対抗する磁気力を与える揚力用コイルとして機能し、空芯コイル２１，３０は模型７に働く抗力に対抗する磁気力を与える抗力対抗用コイルとして機能している。
【０００７】
風洞には、模型７とコイル２１〜３０の他に、各コイルを駆動する電源系、模型７の位置と姿勢とを計測する計測系（カメラのような光学センサ９）、及び模型７の位置と姿勢とを制御する制御系が組み込まれている。図４に示すように、計測系である光学センサ９が検出した模型７の位置姿勢に関する計測データは、パーソナルコンピュータ等の計算機１３に送信され、計算機１３での演算結果が電源アンプ１４によって増幅された後、各コイル２１〜３０に制御された指令電流として送られる。
【０００８】
磁力支持装置２０は、本質的に不安定なシステムであるので、模型７の位置・姿勢角について常にフィードバック制御が行われている。即ち、指令電流によって模型７に磁力が作用し、その結果、出力として模型７の状態変数である位置・姿勢角に変化が表れる。この模型７の位置・姿勢角は、光学センサ９によって検出され、検出値に応じて制御装置である計算機１３が指令電流を出力することで、位置・姿勢角についてフィードバック制御（以下、単に、「位置フィードバック制御」と略す。詳細については後述する。）が行なわれている。
【０００９】
従来の制御において用いられる座標系は、図４に示すように、風洞測定部中心を原点とし、流れの上流方向をｘ軸（逆抗力方向）、鉛直上方（揚力方向）をｚ軸、これらの軸と右手系を作る軸をｙ軸と定め、ｘ軸回りのロール角、ｙ軸回りのピッチ角、ｚ軸回りのヨー角をそれぞれ、φ，θ，ψとしている。これらの軸方向の磁気力及び軸回りのモーメントを作用させるために、表１に示すようにコイルを組み合わせている。表１に記載したコイル番号は、図７に示すコイルに付した番号に対応している。
【表１】
2005157816.tif 000003
【００１０】
模型７が保持する磁気モーメントベクトルをＭ、模型７周り磁場強さベクトルをＨとすると、模型７に作用する磁気力Ｆ、モーメントＮはそれぞれ、式（１）、式（２）のように示される。
【数１】
2005157816.tif 000004
【数２】
2005157816.tif 000005
従来の位置フィードバック制御系では、軸間干渉と、ロール軸（ロール角φ）方向の運動とを考慮せず、また、ピッチ角θとヨー角ψは十分小さいとの仮定の下で、式（１）及び式（２）を各軸方向の電流について線形化し、電流入力に対する模型７の位置・姿勢角の伝達関数を求めて、模型７の位置・姿勢角についてのフィードバック制御が行われている。伝達関数Ｇｍは、模型７に内挿された磁石のｘ軸方向の磁気モーメントをＭｘ、模型７も全体質量又は慣性モーメントをｍ、単位電流当たりの磁場勾配をｈとして、式（３）で表される。
【数３】
2005157816.tif 000006
コイルの巻き数を少なくし、大きな電流を流すことにより、コイルの応答速度を向上させることができる。即ち、コイルのダイナミックスを式（４）で表される一次遅れ系で近似したとき、その時定数γが非常に小さくすることができる（一例として、５軸すべてで最大でおよそ９．６ｍｓｅｃ）。
【数４】
2005157816.tif 000007
【００１１】
上記の各式から、模型系とコイル系とを含めたシステムの状態方程式は、次の式（５）〜式（８）で表すことができる。但し、Ｏ（ベクトル）は適当な次元を持ったゼロ行列である。
【数５】
2005157816.tif 000008
【数６】
2005157816.tif 000009
【数７】
2005157816.tif 000010
【数８】
2005157816.tif 000011
【００１２】
模型７の位置フィードバック制御は、図５に示される制御システムで行われている。磁力支持装置２０の制御システムは、本来、不安定であるため、図５にブロック線図として示すように模型７の位置・姿勢角についてフィードバック制御が行われている。即ち、現在値（状態変数）としての位置・姿勢角Ｙに基づく信号と、目標値設定器１における位置・姿勢角の各目標値Ｒとの偏差Ｅが伝達関数Ｋを持つＰＩ制御器（比例積分コントローラ）２に入力されて制御量が計算される。計算された制御量は、更にバイアス付与器３によるバイアスＢが付加された後、Ｄ／Ａ変換器（デジタル／アナログ変換器）４を経て連続量に変換されて指令信号Ｕが作られる。バイアスＢは、模型７が重力と釣り合うために必要な電流値である。指令信号Ｕは電源アンプ１４で増幅された後に磁力支持装置２０のコイル２１〜３０（図４参照）からなるコイル駆動系５に供給され、コイル駆動系５が生成する磁場（伝達関数Ｇｃ）はコイル磁気力として模型７に作用する。この磁場と模型７の磁石体との相互作用による磁気力と、空気力や重力等の外力源６からの外乱Ｄとが模型７に作用し、その結果、模型７の位置・姿勢角Ｙが出力として得られる。
【００１３】
制御システムには、模型７の現在値としての位置・姿勢角Ｙを検出するため、光学センサ９が備わっている。検出された位置・姿勢角Ｙにはノイズ源８からセンサノイズｎが混入するため、ノイズ除去用に６０Ｈｚ程度の適当なカットオフ周波数を持ち、伝達関数Ｈｎで表される２次バターワース型のノイズカットフィルタ１０にかけられる。光学センサ９は、ＣＣＤセンサのような電荷の蓄積に要する時間として検出から出力までに６ミリ秒程度の遅延時間（むだ時間）Ｌがあるため、検出された模型７の位置・姿勢角の情報は実際にその瞬間の位置・姿勢角とは異なる。また計算機１３はコイル電流やコイル電圧といったコイル２１〜３０に供すべき指令電流Ｕの計算に時間を要するため、力を作用させようとする瞬間には模型７の位置・姿勢角は更に変動している。そのため伝達関数Ｈｐ（ｓ）が式（９）で表される二重位相進み器１１が、これらの時間遅れを補償するのに用いられている。
【数９】
2005157816.tif 000012
光学センサ９による模型７の位置・姿勢角の検出信号は、ノイズカットフィルタ１０を経て、更に二重位相進み要素１１を経るときに位相を進められるが、この位相の二重進みがシステムの安定性に重要な役割を果たしている。この補償された位置・姿勢角信号と目標値Ｒとの偏差ＥがＰＩ制御器２に入力される。ＰＩ制御器２の伝達関数Ｋ（ｓ）が、式（１０）に示されている。このように、磁力支持装置２０は、検出された模型７の位置・姿勢角Ｙに基づいて、模型７の位置・姿勢角についてのフィードバック制御を行っている。
【数１０】
2005157816.tif 000013
【００１４】
上記のように、磁力支持装置２０は、模型７の位置・姿勢角を検出するための光学的な位置センサ（この場合の「位置」は、姿勢角も含む広義の用語）、模型７に磁気力を作用させるための１つ又は複数個のコイル・電磁石２１〜３０、コイル・電磁石２１〜３０に印加する電流値又は電圧値を演算し出力する制御装置（計算機１３等）を基にして構成されており、風洞内において磁気力によって模型７を気流中に空中の任意の位置・姿勢角に保持することができる装置であって、模型７に接触して支持する支持体を設けた場合に生じる支持体と気流との干渉の問題を避けることができ、模型７の動的試験をも可能にする有用な装置である（例えば、非特許文献１参照）。
【非特許文献１】澤田秀夫、国益徹也、「低速風洞用６０ｃｍ磁力支持装置の開発」、日本航空宇宙学会論文集、社団法人日本航空宇宙学会、平成１４年５月、第５０巻、第５８０号、ｐ．１８８−１９５
【００１５】
このように位置・姿勢角をフィードバックする制御系を備えた磁力支持装置は、供試体（模型）の位置・姿勢角を随時検出し、検出した位置・姿勢角に基づいて制御量を計算し、得られた制御量をコイル・電磁石に印加することにより供試体を安定に浮上させることができる。特に非特許文献１に記載のような二重位相進み器とＰＩ制御器とを併用した制御系は、供試体を安定に浮上させることができる制御系の一つであり、供試体を浮上させている状態でも制御定数を変更することが容易にできる実用的な制御系である。この制御系は、供試体が静的な場合又はステップ応答やヒービングといった一次元の簡単な運動をさせる場合には問題は無いが、このときの供試体の運動は制御量に対して位相がずれていることが観測されている。一方、供試体の運動は、コイル・電磁石が形成する磁場に対しては同位相であることも観測されている。この性質を利用して供試体に作用する力を高精度で検出する目的で、磁力支持装置にホール素子を配設する場合もある（例えば、非特許文献２）。
【非特許文献２】澤田秀夫、他３名、“１９９５年の航空宇宙技術研究所における磁力支持天秤装置の研究（STATUS OF MSBS STUDY AT NAL IN 1995 ）”、第３回磁力支持に関する国際シンポジウム（Third International Symposium on Magnetic Suspension Technology ）、（米国）、ｐ．５０５−５１９
【００１６】
一方、磁力支持装置と関連した公知の技術の中には、ホール素子によって磁場を検出し、検出した磁場に基づいて制御量を計算し、それを電磁石に印加することによって物体を一定の位置に保持する装置がある（例えば、特許文献１、特許文献２）。特許文献１には磁気軸受に関する発明についての記載があり、一方、特許文献２には磁気浮上に関する発明についての記載がある。これらの装置の中には光学センサは含まれておらず、また３次元空間の任意の位置に供試体を保持するものではない。
【特許文献１】特開平９−２５９３４号公報（第３頁、第１図）
【特許文献２】米国特許出願公開第６１５４３５３号明細書
【００１７】
このような磁力支持装置において、上記の制御量に対する供試体運動の位相の遅れは、供試体の応答速度の遅れを招くため、制御性能を悪化させる要因の一つとなっている。したがって、この影響が供試体の動安定性を損なうという形で顕在化する可能性が非常に大きく、最悪の場合、供試体の落下、破損といった事故を招きかねない状態にある。供試体をコイル・電磁石からの磁場との相互作用によって浮揚支持させ、供試体の位置・姿勢角についての位置フィードバック制御を行う磁力支持装置において、供試体の位置・姿勢角の制御にはその位置・姿勢角についての情報のみならず、位置・姿勢角を変化させる磁気力を決める磁場が重要な要素になっている。当発明者等は、高速な流れに対しても制御能力を向上させるため、位置フィードバック制御に対して磁場フィードバック制御を追加することにより、速応性の向上に成果を挙げている（非特許文献３参照）。更に、磁場の検出速度を向上させて高速な流れに対する制御能力を向上させる目的で、ホール素子をフィードバック制御に利用する提案もなされている（非特許文献４参照）。
【非特許文献３】須田信一、澤田秀夫、国益徹也、「磁場フィードバック回路の追加による航技研磁力支持天秤装置の制御性能改善」、２００２年１１月２６日、第４５回自動制御連合講演会、計測自動制御学会、ｐ．３７３−３７６
【非特許文献４】Chin E. Lin 他３名、“Improvement on Drag and Control Performance in NCKU Magnetic Suspension Wind Tunnel ”, Proceedings of National Science Council, Republic of China, Part A: Physical Science and Engineering, Vol. 24, No. 5, 2000, （台湾）、pp. 330-340
【００１８】
図６は、磁力支持装置における制御装置の別の例を示すブロック線図である。図６に示す磁力支持装置において、供試体の位置及び姿勢角を制御する位置フィードバック制御は、図５に示すブロック線図で示したフィードバック制御と同様であるので、同じ機能を奏する計測・制御要素には同じ符号を付すことによって再度の説明を省略する。
【００１９】
図６に示す磁気支持装置におけるフィードバック制御のためのブロック線図は、位置・姿勢角についてのフィードバック制御に対して、磁場フィードバック制御を更に組み入れている点で異なっている。磁場フィードバック制御のために新たに追加される計測制御要素は、磁気センサとしてのホール素子１２、プレフィルタ１５、磁場制御器１６及び１７である。磁力支持装置内で発生した磁場はホール素子１２により計測される。プレフィルタ１５は、電流から磁場強さへの変換係数を含み、磁場制御器１６の出力と同次元の量へ変換するとともに、これを増幅する制御要素である。ＰＩ制御器２から出力された制御電流は、プレフィルタ１５により、磁場強さ或いは磁場強さの指標となる量へと変換された後、磁場制御器１６からの出力と差が取られる。磁場制御器１７は、上記差を増幅するとともに、磁場強さ或いは磁場強さの指標となる量からそれに相当する電流値へと変換する。変換後は、バイアス３からのバイアスＢを加えてコイル駆動系５へと出力される。
【００２０】
コイル駆動系５が発生させる磁場を検出する磁気センサとして、ホール素子は、図４に示した磁力支持装置のコイルの配置図を援用すれば、コイル２１〜３０に配設される。ホール素子は、コイル２１〜３０のそれぞれにできるだけ近い場所に配置すること及び気流を乱さないことを理由として、測定部内側に突出しないように測定部壁の内部に埋設されている。ホール素子の埋設位置は、好ましくは、コイル断面の中心位置である。ただし空芯コイル２１，３０については、ホール素子は、コイル２１においては図面に向かって右上、コイル３０においては図面に向かって左下というように、測定部座標原点に対して対称となる位置に配設されている。また、風洞では気流速度を測定するためのピトー静圧管を設置している場合があり、このピトー管にもｘ軸方向の磁場を計測するためのホール素子が配設されている。
【００２１】
上記ホール素子を用いて全５軸の磁場を計測するためには、ホール素子の出力を適当に組み合せて各軸独立となる量を利用する必要がある。この実施例では、風洞の主気流方向であるｘ軸方向の磁気量ｈｘ、ｘ軸方向と直交する横軸であるｙ軸方向の磁気量ｈｙ、ｘ軸方向とｙ軸方向とにそれぞれ直交するｚ軸方向の磁気量ｈｚ、ｙ軸回りのピッチ軸θ方向の磁気量ｈθ、及びｚ軸回りのヨー軸ψ方向の磁気量ｈψを、それぞれ次の式（１１）の（ａ）〜（ｅ）に示すように定義した。
【数１１】
2005157816.tif 000014
ここで右辺のｈはホール素子の出力で、数字の添え字は４に示すコイル２１〜３０を０から９まで順番に番号を振り直したものである。
【００２２】
ホール素子出力から各軸方向の磁場の大きさを求めるに際して、従来は、式（１１）によって求めていたが、より一般的には、式（１２）で表される式の変換行列Ａとオフセット行列Ｂとを求める方法がある。ホール素子の数を６個とするとき、変換行列Ａは５×６の行列であり、オフセット行列Ｂは５×１のベクトルである。
【数１２】
2005157816.tif 000015
ここでの磁場は単位が任意であるため、電流と大きさを同じにすることで制御プログラムを簡単化することができる。即ち、式（１２）において、磁気量ｈを式（１３）に示すように電流値Ｉで置き換え、変換行列Ａとオフセット行列Ｂとを最小二乗法によって求めることができる。
【数１３】
2005157816.tif 000016
【００２３】
例としてｙ軸方向の電流を変化させた際の式（１１）における各磁気量の変化が、図７にグラフとして示されている。図７からは、ｙ軸方向の制御電流の変化に対してｈy のみが大きく変化しており、式（１２）をｙ軸方向の磁場を表す量として利用できることがわかる。その他の軸についても軸同士の対応が取れていることが確認された。また、ｘ軸方向の磁場を表す量として、式（１１）のみならず、式（１４）に示すように空芯コイル２１、３０に配設したホール素子の平均値が利用可能である。
【数１４】
2005157816.tif 000017
また、式（１４）に代えて式（１５）に示すように、ピトー静圧管に配設したホール素子からの出力も利用可能である。
【数１５】
2005157816.tif 000018
【００２４】
各制御要素について、上記のような役割に加え、それぞれがどのような周波数特性を持つべきかを検討する。実施例としてコイル駆動系５が一次遅れ要素と位相遅れ要素の積で表現される場合について述べる。ここでは簡単のため式（１６）のような最も次数の低いものを使って説明する。
【数１６】
2005157816.tif 000019
ここで、コイル駆動系５の極、零点の大きさはａ0 ＞ｂ1 ＞ａ1 となると仮定する。つまり時定数が比較的小さく、それに若干の位相遅れが付加されている形である。コイル駆動系５は式（１６）のような形であってもａ0 が運動周波数に対して充分大きく、位相の遅れも充分小さければ実用上は問題が無い。そのためコイル駆動系５に対する改善点は、コイル駆動系５をフィードバック回路にしたときの時定数を小さくすること、かつｂ1 ，ａ1 の較差を小さくして位相の遅れ度合いを小さくすることの２点に絞られる。
【００２５】
例として、式（１７）に示すように、磁場制御器が積分でありその他の要素は定数である場合を考える。
【数１７】
2005157816.tif 000020
式（１７）に示す制御系の場合、コイル駆動系５の部分の閉ループ伝達関数は、式（１８）に示すように１次／３次の伝達関数になる。
【数１８】
2005157816.tif 000021
この伝達関数は２次遅れのような振る舞いをするので、コイル駆動系５の動特性を変化させ、制御系の調整を複雑にする。そのため、磁場制御器１７の伝達関数Ｋm は分母と分子が同次の有理関数であることが望ましく、その候補として定数、位相進み、近似微分が挙げられる。その中で最も簡単な定数でさえも時定数を小さくし、かつ位相遅れの大きさを小さくすることが可能である。さらに閉ループ化したコイル駆動系のゲインはプレフィルタ１５の伝達関数Ｆm で調整することが効果的であることもわかる。これは式（１８）の中で、磁場制御器１６の伝達関数Ｋm のゲインは分母、分子の両方に現れるのに対して、プレフィルタ１５の伝達関数Ｆm のゲインは分子のみに現れるからである。以上のことから、磁場フィードバック制御の制御要素は全て定数で充分で、ゲインの調整は主にプレフィルタ１５で行うことが効果的であることがわかる。
【００２６】
アンペア単位の電流を入力とし、ボルト単位の磁場を出力とするコイル駆動系５のステップ応答を計測した結果、式（１１）においては単位電流当たり、式（１９）で示す出力をする。
【数１９】
2005157816.tif 000022
ここでｄｉａｇは各要素を対角要素にもつ対角行列を意味する。プレフィルタ１５の伝達関数Ｆm 、磁場制御器１６の伝達関数Ｋm とも、これを基準に調整する。
【００２７】
磁場制御器１７の伝達関数は、これらの逆数を取って式（２０）となる。
【数２０】
2005157816.tif 000023
また、プレフィルタ１５の伝達関数Ｆm は、式（１９）の磁場変換行列１７の行列Ｖm に定数ゲインを掛けた値になる。
【００２８】
数値シミュレーション及び浮揚試験前のプログラム動作確認時に制御電流の大きさからゲインを調整し、実際にはプレフィルタ１５の伝達関数Ｆm 、磁場制御器１６の伝達関数Ｋm は次の式（２１）や式（２２）のように与えた。
【数２１】
2005157816.tif 000024
【数２２】
2005157816.tif 000025
また、磁場フィードバックの加え合わせ点における単位は、磁場強さではなく、磁場強さの指標となる量であるホール素子出力電圧であるため、実際には磁場変換行列(磁場制御器１７に関連する)の行列Ｖm は対角要素が全て１の単位行列を用いた。
【００２９】
このとき、例として、ｙ軸についてはコイル駆動系５は、開ループの場合の式（２３）で表される伝達関数Ｍop（ｓ）から、閉ループの場合の式（２４）で表される伝達関数Ｍcl（ｓ）のようになり、時定数が小さくなりかつ位相遅れ度合いが小さくなっていることが分かる。即ち、コイル駆動系５の位相遅れが改善されていることが判る。
【数２３】
2005157816.tif 000026
【数２４】
2005157816.tif 000027
また、これらの伝達関数の比較を図８に示す。
【００３０】
実際に試験を行い、従来の位置フィードバック制御のみを行う制御系と、磁場フィードバックをも行う制御系とのステップ応答での比較を取ってみたところ、ｙ軸の整定時間は０．７０秒から０．４０秒へ大幅に短縮されたことが確認された。図９は、この試験結果の時間履歴を示した図である。
【００３１】
磁力支持天秤装置において、風洞模型には、常に、空気力が外乱として作用するため、速応性のみならず、外乱に対するロバスト性を考慮する必要がある。また、調整すべきパラメータの数は、磁場フィードバック制御系を含まない位置フィードバック制御だけの場合においても、ＰＩ制御の場合で比例ゲインｋｐと積分時間Ｔｉで二つ（ＰＩＤ制御の場合には更に微分ゲインを含む三つ）、二重位相進み（ｎｐ，Ｔｐ）で二つあり、これらのパラメータがロール軸を除く５軸について存在する。調整すべきパラメータが多いので、ボード線図などの古典的な手法では設計が複雑になり、パラメータ空間法では３つ以上のパラメータについては図的表現が難しく、制御定数の調整が非常に煩雑になると共に各軸の制御性能も統一しづらいという問題がある。また、位置フィードバック制御に磁場フィードバック制御を加える場合には、更に調整すべきパラメータが多くなり過ぎるという問題がある。また、磁場フィードバック制御系は定数ゲインのみで構成されているが、それ以外の制御器を用いた際の制御特性が明確にされていない。即ち、磁場フィードバック制御を含んだ制御では、パラメータ数は、ＰＩ制御と二重位相進み要素の他に更に、電流から磁場強さへの変換を行うプレフィルタの伝達関数及びホール素子が検出した磁場強さの入力を受けて磁場を制御する磁場制御器の伝達関数のパラメータも存在する。更に、調整した制御定数が最適かあるいはそれに近い定数になっているか否かは保証されていない。
【発明の開示】
【発明が解決しようとする課題】
【００３２】
そこで、風洞試験に用いる磁力支持天秤装置のような磁気支持装置における制御系の設計において、ロバスト制御の手法を利用して、多数に上る制御定数の調整の容易化を図る点で解決すべき課題がある。
【００３３】
本発明の目的は、風洞模型を磁力で支持する磁力支持天秤装置のような、少なくとも位置・姿勢角についてのフィードバック制御を行う制御系を備える磁力支持装置において、多数の複雑で且つ最適か否か保証するのが困難な制御定数を持つ制御系の設計を行う際に、ロバスト制御の手法を利用して混合感度問題を設定し、制御定数を合理的に決定することができる磁力支持装置における制御定数決定方法を提供することである。
【課題を解決するための手段】
【００３４】
上記課題を解決するために、本発明による磁力支持装置における制御定数決定方法は、金属及び磁石を有する供試体の位置・姿勢角を検出する位置センサ、前記供試体に磁気力を作用させるための１つ又は複数個のコイル・電磁石、及び前記位置センサにより検出された前記供試体の位置・姿勢角に基づいて前記供試体の位置・姿勢角についてのフィードバック制御をするため前記コイル・電磁石に印加する電流を演算・出力する制御装置を有し、前記供試体を前記磁気力により空中の任意の位置・姿勢角に保持する磁力支持装置において、外乱入力から前記位置・姿勢角の出力までの伝達関数に関連して定められる感度関数と、前記外乱入力から前記コイル・電磁石に印加される前記電流の出力までの伝達関数に関連して定められる相補感度関数とに基づいて、前記伝達関数が有する制御定数を、ロバスト制御で定義されるＨ∞ノルムを最小化する値に定めることから成っている。
【００３５】
この磁力支持装置における制御定数決定方法によれば、供試体の位置・姿勢角についてのフィードバック制御をするための制御装置において、外乱入力から位置・姿勢角の出力までの伝達関数に関連して感度関数が定められ、外乱入力からコイル・電磁石に印加される電流の出力までの伝達関数に関連して相補感度関数が定められる。そのようにして定められた感度関数と相補感度関数とに基づいて、伝達関数が有する制御定数は、ロバスト制御で定義されるＨ∞ノルムを最小化する値として決定される。ロバスト制御は、動特性モデルの同定が不確実である場合やランダムな外乱がある場合に、それらの変動に対して強い特性を示す制御形態である。本発明では、磁力支持装置がパラメータ変動や外乱に曝され易い状態にあることに着目して、その制御装置に用いられるフィードバック制御系にロバスト制御を適用して、フィードバック制御系から定まる感度関数と相補感度関数とに基づいて、フィードバック制御系を構成する制御要素に含まれる多数の制御定数を合理的に定めることができる。
【００３６】
この磁力支持装置における制御定数決定方法において、前記感度関数の最大値を抑えると共に高周波数側にシフトさせる整形を施す重み関数と、前記相補感度関数を高周波数側にシフトさせる整形を施す重み関数をそれぞれ定め、前記各重み関数によって整形された前記感度関数と前記相補感度関数とに対して前記Ｈ∞ノルムを定めることができる。感度関数については、重み関数を用いて、最大値を抑え且つ全体に高周波数側にシフトさせる整形を施すことで、従来よりも高い周波数領域の外乱に対して模型の振動を抑えることが可能になる。また、相補感度関数については、感度関数の整形に応じて、全体に高周波数側にシフトさせる整形を施すことにより、制御電流に高周波成分が必要となることに対応する。
【００３７】
この磁力支持装置における制御定数決定方法において、前記制御装置は、前記コイル・電磁石により形成される磁場を検出する磁気センサと前記磁気センサが検出した前記磁場に基づいて前記コイル・電磁石を制御する磁場制御器とを備え、前記供試体の前記位置・姿勢角についてのフィードバック制御と同時に、前記磁気センサが検出した前記磁場に基づいて前記磁場についてのフィードバック制御を行っており、前記磁場についてのフィードバック制御を行うときの前記感度関数及び前記相補感度関数は、前記位置・姿勢角についてのフィードバック制御のみを行うときの前記磁場の伝達関数を前記磁場についてのフィードバック制御系全体の伝達関数に置換することで定めることができる。制御装置が磁場についてフィードバック制御を行う場合、磁場についてのフィードバック制御は、位置・姿勢角についてのフィードバック制御の内部で行われるので、位置・姿勢角についてのフィードバック制御のみを行うときの磁場の伝達関数は、磁場についてのフィードバック制御系全体の伝達関数に対応している。したがって、対応関係にある伝達関数を置換することで、磁場についてのフィードバック制御を行う磁力支持装置において、感度関数及び相補感度関数を容易に定めることができる。
【発明の効果】
【００３８】
本発明によるにおける制御定数決定方法によれば、風洞模型を磁力で支持する磁力支持天秤装置のような、少なくとも位置・姿勢角についてのフィードバック制御を行う制御系を備える磁力支持装置において、ロバスト制御の手法を利用して混合感度問題を設定し、Ｈ∞ノルムを最小化する値として制御要素の制御定数を定めているので、多数の複雑で且つ最適か否か保証するのが困難な制御定数を持つ制御系の設計を行う際の制御定数を合理的に決定することができる。また、ロバスト制御の手法を利用しているので、決定した制御定数は、最適か、またはそれに十分近いことが保証される。
【発明を実施するための最良の形態】
【００３９】
次に、磁力支持装置において本発明によるロバスト制御の手法を利用した制御定数の決定方法について、実施例を挙げて説明する。図５に示す位置フィードバック制御のみからなる制御系においては、システムの感度関数Ｓ（ｓ）と相補感度関数Ｔ（ｓ）とをそれぞれ式（２５）及び式（２６）に示すように定義する。
【数２５】
2005157816.tif 000028
【数２６】
2005157816.tif 000029
感度関数Ｓ（ｓ）は、外乱入力から位置出力までの伝達関数に関連した関数である。また、相補感度関数は、外乱入力から制御電流出力までの伝達関数に関連する関数である。図６に示すような位置フィードバック制御に加えて磁場フィードバック制御をも行う制御系においては、式（２５）及び式（２６）の中で、磁場の伝達関数Ｇｃ（ｓ）が式（２７）で示す伝達関数Ｇ’ｃ（ｓ）に置き換えられる。
【数２７】
2005157816.tif 000030
【００４０】
図１及び図２は、感度関数と相補感度関数の一例を示す周波数応答線図である。横軸に周波数が採られ、縦軸に大きさ（ゲイン）が採られている。図中の点線で示すグラフが磁場フィードバック制御を含まない図５に示すような制御のシステム場合を示す。また、実線で示すグラフは図５に示すような磁場フィードバック制御を含む制御系の場合を示している。また、太線は、整形用のが重み関数Ｗｓ（ｓ），Ｗｔ（ｓ）を示している。図１に示す感度関数Ｓ（ｓ）によれば、種々の周波数を持つ外乱に対して、模型には約数十Ｈｚの振動が大きく現れることを示している。そこで、感度関数Ｓ（ｓ）に対して、最大値を抑え、且つ全体を高周波数側へシフトさせることによって、従来の制御系よりも、より高い周波数の領域を含む範囲の周波数を持つ外乱に対して、模型の振動を抑えることが可能となる。これに対応して、制御電流もより一層高い周波数成分が必要となるため、相補感度関数Ｔ（ｓ）を全体的に高周波数側へシフトさせることが必要となる。
【００４１】
重み関数Ｗｓ（ｓ），Ｗｔ（ｓ）は、感度関数Ｓ（ｓ）と相補感度関数Ｔ（ｓ）とを最大値を抑え、全体を高周波数側へシフトさせるという整形するための関数であり、それぞれ式（２８）及び式（２９）に示す形に設定する。
【数２８】
2005157816.tif 000031
【数２９】
2005157816.tif 000032
この発明による制御定数の決定方法によれば、ロバスト制御理論から、ベクトルＰ（ｓ）を式（３０）に示すように定義するとき、その絶対値を最小化するように各制御要素に現れる制御定数を決定する。
【数３０】
2005157816.tif 000033
ここで、絶対値‖Ｐ（ｓ）‖∞は、ロバスト制御では、Ｈ∞ノルムと称される量であり、Ｐ（ｓ）が１次元の場合には、そのゲインの最大値となり、多次元の場合には、最大特異値となる。ロバスト制御の理論から、‖Ｐ（ｓ）‖∞が小さければその制御系は良い制御性Ｗを備えているということが言える。
【００４２】
数回の試行の結果、図５に示すような位置フィードバック制御のみを備えた制御系の場合には、ＰＩ制御器２の伝達関数Ｋ（ｓ）、二重位相進み器１１の伝達関数Ｈｐ（ｓ）は、それぞれ式（３１）、式（３２）に示すものに決定された。
【数３１】
2005157816.tif 000034
【数３２】
2005157816.tif 000035
【００４３】
また、図６に示すような、位置フィードバック制御のみならず磁場フィードバック制御をも備えた制御系の場合には、ＰＩ制御器２の伝達関数Ｋ（ｓ）、二重位相進み器１１の伝達関数Ｈｐ（ｓ）、並びに、磁場フィードバック制御系の磁場制御器１６，１７の伝達関数Ｋｍ（ｓ）及びＨｍ（ｓ）は、それぞれ式（３３）、式（３４）、式（３５）及び式（３６）に示すものに決定された。
【数３３】
2005157816.tif 000036
【数３４】
2005157816.tif 000037
【数３５】
2005157816.tif 000038
【数３６】
2005157816.tif 000039
ここで、重み関数Ｗｓ（ｓ），Ｗｔ（ｓ）は、両方の制御系に対して式（３７）、式（３８）に示すような共通した形に設定された。
【数３７】
2005157816.tif 000040
【数３８】
2005157816.tif 000041
図３には、模型位置の時間変化を示す数値シミュレーションの結果が示されている。図３に示すように、従来の位置フィードバック制御のみの制御系と、磁場フィードバック制御も行う制御系とでは、外乱に対しての揺動には大きな差は見受けられなかったが、その大きさは非常に小さく抑えることが期待できる。
【００４４】
このように、種々の周波数を含む外乱を受けやすい磁気支持装置において、ロバスト制御の手法を適用することで、制御素子の制御定数を最適な又はそれに近い状態に容易に且つ効率的に自動調整・決定をすることができる。本決定方法は、制御定数の個数が多数ある場合に、非常に効率的である。また、模型の支持軸は５軸（又は６軸）存在するが、全ての軸について、制御性能を統一化しやすい。更に、与えられた制御性能要求（重み関数）に対して最適な制御定数を求めることができる。
【図面の簡単な説明】
【００４５】
【図１】この発明による磁気支持装置における制御定数決定方法で用いられる感度関数の一例を示す線図である。
【図２】この発明による磁気支持装置における制御定数決定方法で用いられる相補感度関数の一例を示す線図である。
【図３】位置フィードバック制御のみの制御系と磁場フィードバック制御も行う制御系とにおける模型位置の時間変化を示す数値シミュレーションの結果を示すグラフである。
【図４】磁力支持型風洞と磁力支持装置、及び磁力支持装置に用いられる計測系と制御系の概要を示す斜視図である。
【図５】位置フィードバック制御のみを行う磁気支持装置の制御ブロックの一例を示す図である。
【図６】位置フィードバック制御と磁場フィードバック制御とを行う磁気支持装置の制御装置の一例を示すブロック線図である。
【図７】ｙ軸方向の電流を変化させた際の各磁気量の変化を示すグラフである。
【図８】磁気支持装置のコイル駆動系のｙ軸について、開ループ及び閉ループの場合の伝達関数の比較を示す図である。
【図９】位置フィードバック制御のみを行う制御系と、磁場フィードバックをも行う制御系とのステップ応答を示すグラフである。
【符号の説明】
【００４６】
１  目標値設定器                    ２  ＰＩ制御器（伝達関数Ｋ）
３  バイアス付与器                  ４  Ｄ／Ａ変換器
５  コイル駆動系
７  風洞模型                        ８  ノイズ源
９  光学センサ（遅延時間Ｌ）        １０  ノイズカットフィルタ（伝達関数Ｈｎ）  １１  二重位相進み器（伝達関数Ｈｐ）    １２  ホール素子（磁気センサ）
１３  計算機                        １４  電源アンプ
１５  プレフィルタ（伝達関数Ｆm ）  １６  磁場制御器（伝達関数Ｋm ）
１７  磁場制御器（伝達関数Ｈm ）    
２０  磁力支持装置                  ２１，３０  空芯コイル
２２〜２９  コイル
ｘ，ｙ，ｚ  右手系を構成する座標軸
φ  ロール角                  θ  ピッチ角          ψ  ヨー角
Ｙ  位置・姿勢角              Ｒ  目標値            Ｅ  偏差
Ｂ  バイアス                  Ｄ  外乱              Ｕ  指令信号
Ｎ  センサノイズ  
Ｍｇ    磁場の伝達関数
ｈｘ  ｈｙ  ｈｚ  ｈθ  ｈψ  磁気量
Ｍop  コイル駆動系５が開ループの場合の伝達関数
Ｍcl  コイル駆動系５が閉ループの場合の伝達関数
Ｓ（ｓ）  感度関数                  Ｔ（ｓ）  相補感度関数
Ｗｓ（ｓ），Ｗｔ（ｓ）  重み関数
【図１】
2005157816.tif 000042
【図２】
2005157816.tif 000043
【図３】
2005157816.tif 000044
【図４】
2005157816.tif 000045
【図５】
2005157816.tif 000046
【図６】
2005157816.tif 000047
【図７】
2005157816.tif 000048
【図８】
2005157816.tif 000049
【図９】
2005157816.tif 000050
【手続補正書】
【提出日】平成15年12月26日(2003.12.26)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】０００９
【補正方法】変更
【補正の内容】
【０００９】
従来の制御において用いられる座標系は、図４に示すように、風洞測定部中心を原点とし、流れの上流方向をｘ軸（逆抗力方向）、鉛直上方（揚力方向）をｚ軸、これらの軸と右手系を作る軸をｙ軸と定め、ｘ軸回りのロール角、ｙ軸回りのピッチ角、ｚ軸回りのヨー角をそれぞれ、φ，θ，ψとしている。これらの軸方向の磁気力及び軸回りのモーメントを作用させるために、表１に示すようにコイルを組み合わせている。表１に記載したコイル番号は、図７に示すコイルに付した番号に対応している。
【表１】
2005157816.tif 000051