(57)【要約】【課題】高速回転する球状物体の角速度ベクトルを非接触で計測する手法およびその機構に関するものである。【解決手段】球体に、目的に応じて選択した適当な物理量の偏向を与えるマーカーをつけ、球体の周囲に配置したセンサによってこれを計測することによって、球体の角速度ベクトルを算出するものである。また、従来の手法においては限定的であった、回転軸方向の計測可能範囲を大幅に拡大するものである。【選択図】図１
2010025856.tif 000002
【特許請求の範囲】
【請求項１】
表面上の基準点から一定の距離以内にある領域により構成される第１領域と、前記表面上の該第１領域を除く領域により構成される第２領域と、を区別するマーカーが付された回転する球体と、
該球体に対向するように非接触で配置され、前記マーカーを認識する少なくとも３つのセンサと、
各センサの認識の結果に基づいて、前記球体における前記基準点と該球体の回転軸とを結ぶ大円がなす基準線が各センサを通過する時刻を求め、該時刻に基づいて、前記球体の回転軸方向を算出する算出手段と、
を具備することを特徴とする回転軸方向算出装置。
【請求項２】
表面上の基準点からの距離に従って同心円状に物理量が単調増加又は単調減少するマーカーが付された回転する球体と、
該球体に対向するように非接触で配置され、前記マーカーを認識する少なくとも３つのセンサと、
各センサの認識の結果に基づいて、前記球体における前記基準点と該球体の回転軸とを結ぶ大円がなす基準線が各センサを通過する時刻を求め、該時刻に基づいて、前記球体の回転軸方向を算出する算出手段と、
を具備することを特徴とする回転軸方向算出装置。
【請求項３】
前記球体に前記マーカーを付することが、該球体の表面の光の反射率を変えることであり、
前記センサが光学センサである、請求項１又は請求項２に記載の回転軸方向算出装置。
【請求項４】
前記球体に前記マーカーを付することが、該球体を磁極化させることであり、
前記センサがホール素子である、請求項１又は請求項２に記載の回転軸方向算出装置
【請求項５】
前記球体に前記マーカーを付することが、該球体の表面の導電率を変えることであり、
前記センサが渦電流センサである、請求項１又は請求項２に記載の回転軸方向算出装置。
【請求項６】
前記球体に前記マーカーを付することが、該球体の表面の音の反射率を変えることであり、
前記センサが音センサである、請求項１又は請求項２に記載の回転軸方向算出装置。
【請求項７】
前記球体に前記マーカーを付することが、該球体の電波の吸収率を変えることであり、
前記センサが電波センサである、請求項１又は請求項２に記載の回転軸方向算出装置。
【請求項８】
前記算出手段が、前記球体の角速度の大きさも求めることにより、該球体の角速度ベクトルも算出する、請求項１から請求項７のいずれかに記載の回転軸方向算出装置。
【請求項９】
前記算出手段が、前記時刻を用いて、各センサの前記球体の回転軸に対する該球体の球面上の位相角を求め、該位相角と前記球面上における各センサの位置とによって、前記球体の回転軸方向を算出する、請求項１から請求項８のいずれかに記載の回転軸方向算出装置。
【請求項１０】
前記センサの数を増やすことにより、前記回転軸方向の算出可能範囲を増加させている、請求項１から請求項９のいずれかに記載の回転軸方向算出装置。
【請求項１１】
前記球体には、前記第１領域の大きさと前記第２領域の大きさとが等しくなるように前記マーカーが付されていることにより、前記回転軸方向の算出可能範囲を増加させている、請求項１に記載の回転軸方向算出装置。
【請求項１２】
前記基準点に対して直交する方向に前記球体の最大慣性主軸を近づけることよって、前記球体の回転軸の収束方向を前記直交する方向に近づけ、前記基準線が各センサを通過する確率を高くしている、請求項１に記載の回転軸方向算出装置。
【請求項１３】
表面上の基準点から一定の距離以内にある領域により構成される第１領域と、前記表面上の該第１領域を除く領域により構成される第２領域と、を区別するマーカーが付された回転する球体と、
該球体に対向するように非接触で配置され、前記マーカーを認識する少なくとも３つのセンサと、を含む計測装置に接続可能な算出装置であって、
各センサの認識の結果を取得する手段と、
該認識の結果に基づいて、前記球体における前記基準点と該球体の回転軸とを結ぶ大円がなす基準線が各センサを通過する時刻を求め、該時刻に基づいて、前記球体の回転軸方向を算出する手段と、
を具備することを特徴とする算出装置。
【請求項１４】
表面上の基準点からの距離に従って同心円状に物理量が単調増加又は単調減少するマーカーが付された回転する球体と、
該球体に対向するように非接触で配置され、前記マーカーを認識する少なくとも３つのセンサと、を含む計測装置に接続可能な算出装置であって、
各センサの認識の結果を取得する手段と、
該認識の結果に基づいて、前記球体における前記基準点と該球体の回転軸とを結ぶ大円がなす基準線が各センサを通過する時刻を求め、該時刻に基づいて、前記球体の回転軸方向を算出する手段と、
を具備することを特徴とする算出装置。
【請求項１５】
表面上の基準点から一定の距離以内にある領域により構成される第１領域と、前記表面上の該第１領域を除く領域により構成される第２領域と、を区別するマーカーが付された回転する球体を準備するステップと、
該球体に対向するように非接触で配置され、前記マーカーを認識する少なくとも３つのセンサを準備するステップと、
各センサの認識の結果に基づいて、前記球体における前記基準点と該球体の回転軸とを結ぶ大円がなす基準線が各センサを通過する時刻を求め、該時刻に基づいて、前記球体の回転軸方向を算出するステップと、
を含むことを特徴とする算出方法。
【請求項１６】
表面上の基準点からの距離に従って同心円状に物理量が単調増加又は単調減少するマーカーが付された回転する球体を準備するステップと、
該球体に対向するように非接触で配置され、前記マーカーを認識する少なくとも３つのセンサを準備するステップと、
各センサの認識の結果に基づいて、前記球体における前記基準点と該球体の回転軸とを結ぶ大円がなす基準線が各センサを通過する時刻を求め、該時刻に基づいて、前記球体の回転軸方向を算出するステップと、
を含むことを特徴とする算出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は球状の物体の角速度ベクトルを非接触で外部から計測する手法および装置に関するものである。
【背景技術】
【０００２】
従来、球形の物体の角速度ベクトルを外部から非接触で計測する手法はいくつか考えられてきた。１つには、対象物を磁極化させ、ホール素子を用いてその磁場の変化を計測することにより角速度ベクトルを求める手法がある。
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかし、磁場を発生することが許されない対象や磁化できない対象においては、この手法は適用できなかった。また、光学センサを用いた手法も存在するが、許容される対象の回転軸の範囲が有限である場合に限定されている。画像処理によって回転軸を求め、角速度ベクトルを算出する手法もあるが、装置の規模が大きくなり、また処理に時間がかかるため高速回転する物体には適さない。このように、従来の角速度ベクトルを計測する手法には特徴があるが、磁場を用いずに、任意の軸周りに高速回転する球状物体の角速度ベクトルを非接触で計測する有効な手法は存在していない。
【０００４】
本発明は、高速回転する球状物体の角速度ベクトルを非接触で計測する手法およびその機構に関するものである。球体に、目的に応じて選択した適当な物理量の偏向を与えるマーカーをつけることによって、１回の計測期間にわたって球体の回転が一定であると仮定して球体の角速度ベクトルを算出するものである。また、従来の手法においては限定的であった、回転軸方向の計測可能範囲を大幅に拡大するものである。
【課題を解決するための手段】
【０００５】
本発明に係る回転軸方向算出装置は、表面上の基準点から一定の距離以内にある領域により構成される第１領域と、前記表面上の該第１領域を除く領域により構成される第２領域と、を区別するマーカーが付された回転する球体と、該球体に対向するように非接触で配置され、前記マーカーを認識する少なくとも３つのセンサと、各センサの認識の結果に基づいて、前記球体における前記基準点と該球体の回転軸とを結ぶ大円がなす基準線が各センサを通過する時刻を求め、該時刻に基づいて、前記球体の回転軸方向を算出する算出手段と、を具備することを特徴とする。
【０００６】
本発明に係る回転軸方向算出装置は、表面上の基準点からの距離に従って同心円状に物理量が単調増加又は単調減少するマーカーが付された回転する球体と、該球体に対向するように非接触で配置され、前記マーカーを認識する少なくとも３つのセンサと、各センサの認識の結果に基づいて、前記球体における前記基準点と該球体の回転軸とを結ぶ大円がなす基準線が各センサを通過する時刻を求め、該時刻に基づいて、前記球体の回転軸方向を算出する算出手段と、を具備することを特徴とする。
【０００７】
本発明に係る算出装置は、表面上の基準点から一定の距離以内にある領域により構成される第１領域と、前記表面上の該第１領域を除く領域により構成される第２領域と、を区別するマーカーが付された回転する球体と、該球体に対向するように非接触で配置され、前記マーカーを認識する少なくとも３つのセンサと、を含む計測装置に接続可能であって、各センサの認識の結果を取得する手段と、該認識の結果に基づいて、前記球体における前記基準点と該球体の回転軸とを結ぶ大円がなす基準線が各センサを通過する時刻を求め、該時刻に基づいて、前記球体の回転軸方向を算出する手段と、を具備することを特徴とする。
【０００８】
本発明に係る算出装置は、表面上の基準点からの距離に従って同心円状に物理量が単調増加又は単調減少するマーカーが付された回転する球体と、該球体に対向するように非接触で配置され、前記マーカーを認識する少なくとも３つのセンサと、を含む計測装置に接続可能であって、各センサの認識の結果を取得する手段と、該認識の結果に基づいて、前記球体における前記基準点と該球体の回転軸とを結ぶ大円がなす基準線が各センサを通過する時刻を求め、該時刻に基づいて、前記球体の回転軸方向を算出する手段と、を具備することを特徴とする。
【０００９】
本発明に係る算出方法は、表面上の基準点から一定の距離以内にある領域により構成される第１領域と、前記表面上の該第１領域を除く領域により構成される第２領域と、を区別するマーカーが付された回転する球体を準備するステップと、該球体に対向するように非接触で配置され、前記マーカーを認識する少なくとも３つのセンサを準備するステップと、各センサの認識の結果に基づいて、前記球体における前記基準点と該球体の回転軸とを結ぶ大円がなす基準線が各センサを通過する時刻を求め、該時刻に基づいて、前記球体の回転軸方向を算出するステップと、を含むことを特徴とする。
【００１０】
本発明に係る算出方法は、表面上の基準点からの距離に従って同心円状に物理量が単調増加又は単調減少するマーカーが付された回転する球体を準備するステップと、該球体に対向するように非接触で配置され、前記マーカーを認識する少なくとも３つのセンサを準備するステップと、各センサの認識の結果に基づいて、前記球体における前記基準点と該球体の回転軸とを結ぶ大円がなす基準線が各センサを通過する時刻を求め、該時刻に基づいて、前記球体の回転軸方向を算出するステップと、を含むことを特徴とする。
【発明を実施するための最良の形態】
【００１１】
以下、本発明の実施の形態について図面を参照して説明する。
【００１２】
球１について、その表面上の適当な基準点２から一定の距離以内にある特定領域によって、球面をマーク面（上記特定領域により構成される「第１領域」）３と非マーク面４（上記特定領域を除く領域により構成される「第２領域」）の２つの部分に区別できるようなマーカーをつけ、それぞれの面を区別して認識できるセンサ５を球１の周辺に慣性系に固定して配置する（図１）。このマーカーとセンサの組み合わせは、要求に応じて変更可能であり、センサに反応しかつ明確に２面を区別できるものであれば何でもよい。球１が回転する（図２におけるステップ２０１）と、センサ５はその面している球面のマーカーを認識してON/OFFのデジタル出力を与える（ステップ２０２、２０３）。回転によってセンサ５がマーク面３、非マーク面４を交互に認識する場合は、球１の回転速度Ωに一致した周波数のON/OFFの矩形波が得られる（図３）。これにより、球１の角速度ベクトルの絶対値Ωを得ることができる（ステップ２０４）。
【００１３】
次に、慣性系における球の回転軸方向Ｓを求める。まず、球の回転軸方向Ｓと基準点方位Ｏを結ぶ球面上の基準線７を考える。さらに、あるセンサの出力がOFFからONへ変わる時刻における計測点をＡ、ONからOFFへと変わる時刻における計測点をＢとし、この間センサが通過する球表面上の点の軌跡ＡＢについて考える（図４）。まず球面三角法において、マーカーは基準点２を中心として描かれているので、辺ＯＡと辺ＯＢの角は等しい。一方、回転軸Ｓとセンサ方向のなす角は一定であるため、辺ＳＡと辺ＳＢの角は等しい。よって、△ＯＳＡ≡△ＯＳＢであり、∠ＯＳＡ＝∠ＯＳＢである。これより、辺ＳＯの延長線（すなわち基準線７）と軌跡ＡＢの交点をＭとおくと、∠ＭＳＡ＝∠ＭＳＢも成立する。従って、球の回転が一定であるとすると、出力がOFFからONへ変わる時刻とONからOFFへと変わる時刻の中間の時刻においては、センサはちょうど点Ｍを計測していることになり、センサの直下を基準線７が通過しているといえる。これより、各センサのこの中間時刻を計算することにより基準線７が各センサの直下を通過する時刻の相対関係がわかる。この時刻の相対関係により、２つのセンサｉ、ｊについて、センサｉの通過時刻からセンサｊの通過時刻までの時刻をｔijとすれば（ステップ２０５）、回転軸に対する２つのセンサの相対位相角θijを次式で求めることができる（ステップ２０６）。
【数１】
2010025856.tif 000003
ただしＴは回転の周期であり、
【数２】
2010025856.tif 000004
である。
【００１４】
球１の中心を原点とした、慣性系に固定された座標系を考える。球１の周辺にｎ個のセンサを慣性系に固定して配置する。これらのセンサの方向を、単位ベクトルＳk(ｋ＝１，２，…ｎ)で表す。ただしｎは３以上の整数である。またいずれの２つのセンサ方向も球の中心と１直線上に乗らないように配置する。
【００１５】
２つのセンサＳｉ、Ｓｊと回転軸Ｓを頂点とする球面三角形について、図５のように各種記号を定義する。ＳとＳi間の辺の角をαi、ＳとＳj間の辺の角をαjとする。また、２つのセンサ間の辺の角をΨijとおく。回転軸方向の単位ベクトルをＳ＝[ｓx,ｓy,ｓz]T、センサｋ方向の単位ベクトルをＳk＝[ｓxk,ｓyk,ｓzk]Tとおくと、次式が成立する（ステップ２０９）。
【数３】
2010025856.tif 000005
ここで、センサ方向は既知であるので、この式の未知数はｓx、ｓy、ｓzの３個である。よって３個のセンサについてαを求めることができれば回転軸Ｓは(3)式を連立させて解くことにより求めることができる。
【００１６】
一方、球面三角法により次式が得られる（ステップ２０７）。
【数４】
2010025856.tif 000006
(4)式により得られる式の数は、ｉ、ｊの組み合わせの数であり、
【数５】
2010025856.tif 000007
である。これらの式が互いに独立であるとして連立方程式をたて解くことにより、ｎ≧３で各センサについてcosαは求めることができる（ステップ２０８）。すなわち(3)式を解くことができ、回転軸を求めることができる（ステップ２１０）。
【００１７】
３つのセンサａ、ｂ、ｃを直交３座標軸上に配置した場合には連立方程式は解析的に解くことができ、反復計算をすることなくcosαが求められ、次式で直接回転軸ベクトルＳを求めることができる。
【数６】
2010025856.tif 000008
ここで各成分の符号は次式で決定する。
【数７】
2010025856.tif 000009
【００１８】
以上で求めた角速度ベクトルのノルムΩと回転軸Ｓにより、球の角速度ベクトルを算出できる（ステップ２１１）。
【００１９】
センサはマーカーを感知してON/OFF信号を出力するだけでよいので、高速でのサンプリングも可能である。このため高速回転する場合についても高精度で計測することが可能となる。
【００２０】
この手法は、最低３つのセンサがあれば構成できるが、センサが感知できない軸周りに回転すると計測はできなくなる。このような計測できない領域の評価の仕方として、単位球表面においてセンサがマーカーのON/OFFを感知できる領域の割合を、基準点と回転軸の位置関係に対して表す手法を用いる。球面において、基準点と回転軸方向がなす辺の角をΦｓ、マーカーの半径をΦｍとおくと、図６より、センサがマーカーのON/OFFを感知できる領域は、回転軸方向からの角度Φによって次式で表される。
【数８】
2010025856.tif 000010
この領域の表面積Ｓｅを求めると、
【数９】
2010025856.tif 000011
となる。Φｓ、Φｍはいずれも0°〜180°の値であり、両者とも90°のときＳｅは最大となる。すなわち、Φｓ、Φｍを90°に近づければ、センサが感知できない領域を減らすことができる。
【００２１】
Φｍを90°に近づけることは、マーカーが球を半分ずつに分割するように付することで実現できる。
Φｓを90°に近づけることは、回転球の最大慣性主軸を基準点方向と直交するように与えるよう調整して、エネルギー消散によって回転軸を収束させることによって可能である。
【００２２】
以上の手法に加えて、センサの数を増やすことによって、いずれか３つのセンサが感知できる確率を増やす手法が存在する。この極限として、全方向を感知できる３次元センサを用いれば、Φｍが0°あるいは180°の場合のみセンサのON/OFFが感知できない点となる。逆に、感知できないことで角速度の大きさは求められないが回転軸をこの点（基準点）に限定することは可能となる。すなわち、どの方向に回転軸があっても、回転軸を特定できることになる。
【００２３】
マーカーとセンサの組み合わせについて、マーカーを基準点からの角度に依存して変化する物理量の連続的な変化とし、これを連続的に感知できるセンサを用いれば、前述のようなセンサが感知できない領域についての問題は回避できる。
球について、その表面上の適当な基準点を中心とし、この基準点からの球面上の距離に依存して適当な物理量が単調増加あるいは単調減少するマーカーをつけ、その物理量を計測できるセンサを球の周辺に慣性系に固定して配置する。このマーカーとセンサの組み合わせは、要求に応じて変更可能である。球が回転すると、センサが回転軸方向にある場合と、回転軸方向と基準点方向が一致している場合を除いて、センサはその面している球面のマーカーの物理量を計測し、球の回転速度Ωに一致した周波数の出力が得られる（図７）。これにより、球の角速度ベクトルの絶対値Ωを得ることができる。
【００２４】
次に、慣性系における球の回転軸方向Ｓを求める。球の回転軸と基準点を結ぶ球面上の基準線を考える。すると、出力が最大となる時刻においては、そのセンサの直下を基準線が通過しているといえる。すなわち、基準線が各センサの直下を通過する時刻の相対関係がわかる。この時刻の相対関係により、２つのセンサｉ、ｊについて、センサｉの通過時刻からセンサｊの通過時刻までの時刻をｔijとすれば、回転軸に対する２つのセンサの相対位相角θijを求めることができ、以下ON/OFF出力の場合と同様の議論で回転軸を求めることができる。
【００２５】
一例として、人工衛星の姿勢制御装置（3次元球体リアクションホイール）の試験機として製作した，直径11/8インチの鋼球を３自由度回転させる装置においてその角速度ベクトルを計測する装置の例を図１に示す。この装置は、ベアリングに用いられる高真球度の鋼球を浮上させ３自由度回転させることで任意の軸周りに反トルクを得ることができる装置である。その姿勢制御装置としての役割から、数千rpmで回転する剛球の角速度ベクトルを計測してフィードバックする必要があるため、本発明における計測手法を導入した。この回転する球はブラスト加工によって半面ずつ反射率を変えてある。すなわちΦｍは90°である。さらに、反射率の境界面上に軸穴をあけ、最大慣性主軸方向が基準点方向と直交するようにしている。これにより、Φｓが90°に収束するように工夫している。ブラスト加工による反射率の変化を検出するため、球の周辺にトランジスタ出力タイプのフォトダイオードを配置し、出力電圧をコンパレータによりデジタル信号に変換する。このパルス信号はFPGAによってサンプリングする。FPGAは、各センサの反応時間履歴から、パルス周期、継続してONとなる時間、を記録し、さらに各センサ間のONになるタイミングのずれの時間も記録する。この情報はCPUへと送られ、この情報を元に角速度ベクトルを算出する。
【００２６】
なお、図１では、本発明を分かり易く説明するために、回転球１を非接触で支持する軸受けの表現は省略されている。
【００２７】
また、フォトダイオードにより得られた信号をサンプリングするＦＰＧＡ、及び、このＦＰＧＡから受信した情報に基づいて回転球１の回転軸方向及び角速度ベクトルを算出するＣＰＵ（算出手段）の表現も、説明の便宜上、省略されている。これらＦＰＧＡ及びＣＰＵは、図１に示した計測装置の内部に搭載することもでき、或いは、この計測装置とは別の装置（算出装置）に搭載することもできる。
【００２８】
この他にも、マーカー、センサの組み合わせは多岐にわたって考えられる。そのマーカー―センサの組み合わせ例は以下のようにあげられる(括弧内は観測する物理量を示す)。
・各種光反射塗料―光学センサ(光)
・球の磁極化―ホール素子(磁場)
・導電性素材―渦電流センサ(導電率)
・音の反射板―超音波センサ(音波)
・電波吸収体―電波センサ(電波)
【産業上の利用可能性】
【００２９】
本発明は、慣性諸元測定装置、３次元自由運動模擬装置（無重力模擬）、人工衛星の姿勢制御装置、全方向移動車輪の回転検出、トラックボールの回転計測、各種球技におけるボールの回転計測、FA(ファクトリーオートメーション)における製造、検査の工程での球体の回転計測等において、利用することができる。
【図面の簡単な説明】
【００３０】
【図１】回転球と回転軸、マーカーの関係を表した計測装置の模式図。
【図２】角速度ベクトルの導出の流れ図。
【図３】球の回転により得られるセンサ出力の時間履歴の模式図。
【図４】センサ出力がOFFからONへ変わる時刻とONからOFFへと変わる時刻との中間の時刻に基準線７がセンサの直下を通過することを示す説明図。
【図５】球面三角法に用いる記号の定義の説明図。
【図６】センサが感知できる領域の説明図。
【図７】マーカーの物理量が連続的に変化する場合の説明図。
【符号の説明】
【００３１】
１  回転球
２  基準点
３  マーク面
４  非マーク面
５  センサ
６  回転軸
７  基準線
【図２】
2010025856.tif 000012
【図１】
2010025856.tif 000013
【図３】
2010025856.tif 000014
【図４】
2010025856.tif 000015
【図５】
2010025856.tif 000016
【図６】
2010025856.tif 000017
【図７】
2010025856.tif 000018