本発明の実施形態による倒立振子型移動体は、全方向移動車両としての移動制御に関する技術を前提とし、他の移動体との連携動作を実現するための構成を特徴部として更に備える。換言すれば、本実施形態による倒立振子型移動体は、全方向移動車両としての基本的な移動制御に関する構成に加え、他の移動体と連携して動作するための構成を備える。
そこで、本実施形態では、「1.前提技術(全方向移動車両)」の項において、前提技術となる着座型の全方向移動車両について一通り説明した後に、「2.連携技術」の項において、松葉杖の様な歩行補助具に適用した場合を例として、本実施形態による倒立振子型移動体の特徴部である連携動作に関わる技術内容を説明する。
1.前提技術(全方向移動車両)
図1〜図6を参照して、本実施形態による倒立振子型移動体の前提技術となる全方向移動車両の構造を説明する。
図1及び図2に示すように、本実施形態における全方向移動車両1は、乗員(運転者)の搭乗部3と、床面に接地しながら該床面上を全方向(前後方向及び左右方向を含む2次元的な全方向)に移動可能な移動動作部5と、この移動動作部5を駆動する動力を該移動動作部5に付与するアクチュエータ装置7と、これらの搭乗部3、移動動作部5及びアクチュエータ装置7が組付けられた基体9とを備える。
ここで、本実施形態の説明では、「前後方向」、「左右方向」は、それぞれ、搭乗部3に標準的な姿勢で搭乗した乗員の上体の前後方向、左右方向に一致もしくはほぼ一致する方向を意味する。なお、「標準的な姿勢」は、搭乗部3に関して設計的に想定されている姿勢であり、乗員の上体の体幹軸を概ね上下方向に向け、且つ、上体を捻ったりしていない姿勢である。
この場合、図1においては、「前後方向」、「左右方向」はそれぞれ、紙面に垂直な方向、紙面の左右方向であり、図2においては、「前後方向」、「左右方向」はそれぞれ、紙面の左右方向、紙面に垂直な方向である。また、本実施形態の説明では、参照符号に付する添え字「R」,「L」は、それぞれ車両1の右側、左側に対応するものという意味で使用する。
基体9は、移動動作部5及びアクチュエータ装置7とが組付けられた下部フレーム11と、この下部フレーム11の上端から上方に延設された支柱フレーム13とを備える。
支柱フレーム13の上部には、該支柱フレーム13から前方側に張り出したシートフレーム15が固定されている。そして、このシートフレーム15上に、乗員が着座するシート3が装着されている。本実施形態では、このシート3が乗員の搭乗部となっている。従って、本実施形態における全方向移動車両1(以降、単に車両1という)は、乗員がシート3に着座した状態で、床面上を移動するものである。
また、シート3の左右には、シート3に着座した乗員が必要に応じて把持するためのグリップ17R,17Lが配置され、これらのグリップ17R,17Lがそれぞれ、支柱フレーム13(又はシートフレーム15)から延設されたブラケット19R,19Lの先端部に固定されている。
下部フレーム11は、左右方向に間隔を存して二股状に対向するように配置された一対のカバー部材21R,21Lを備える。これらのカバー部材21R,21Lの上端部(二股の分岐部分)は、前後方向の軸心を有するヒンジ軸23を介して連結され、カバー部材21R,21Lの一方が他方に対して相対的にヒンジ軸23の周りに揺動可能となっている。この場合、カバー部材21R,21Lは、図示を省略するバネによって、カバー部材21R,21Lの下端部側(二股の先端側)が狭まる方向に付勢されている。
また、カバー部材21R,21Lのそれぞれの外面部には、前記シート3に着座した乗員の右足を載せるステップ25Rと左足を載せるステップ25Lとが各々、右向き、左向きに張り出すように突設されている。
移動動作部5及びアクチュエータ装置7は、下部フレーム11のカバー部材21R,21Lの間に配置されている。これらの移動動作部5及びアクチュエータ装置7の構造を図3〜図6を参照して説明する。
なお、本実施形態で例示する移動動作部5及びアクチュエータ装置7は、例えば前記特許文献2の図1に開示されているものと同じ構造のものである。従って、本実施形態の説明においては、移動動作部5及びアクチュエータ装置7の構成に関して、前記特許文献2に記載された事項については、簡略的な説明に留める。
本実施形態では、移動動作部5は、ゴム状弾性材により円環状に形成された車輪体であり、ほぼ円形の横断面形状を有する。この移動動作部5(以降、車輪体5という)は、その弾性変形によって、図5及び図6の矢印Y1で示す如く、円形の横断面の中心C1(より詳しくは、円形の横断面中心C1を通って、車輪体5の軸心と同心となる円周線)の周りに回転可能となっている。
この車輪体5は、その軸心C2(車輪体5全体の直径方向に直交する軸心C2)を左右方向に向けた状態で、カバー部材21R,21Lの間に配置され、該車輪体5の外周面の下端部にて床面に接地する。
そして、車輪体5は、アクチュエータ装置7による駆動(詳細は後述する)によって、図5の矢印Y2で示す如く車輪体5の軸心C2の周りに回転する動作(床面上を輪転する動作)と、車輪体5の横断面中心C1の周りに回転する動作とを行なうことが可能である。その結果、車輪体5は、それらの回転動作の複合動作によって、床面上を全方向に移動することが可能となっている。
アクチュエータ装置7は、車輪体5と右側のカバー部材21Rとの間に介装される回転部材27R及びフリーローラ29Rと、車輪体5と左側のカバー部材21Lとの間に介装される回転部材27L及びフリーローラ29Lと、回転部材27R及びフリーローラ29Rの上方に配置されたアクチュエータとしての電動モータ31Rと、回転部材27L及びフリーローラ29Lの上方に配置されたアクチュエータとしての電動モータ31Lとを備える。
電動モータ31R,31Lは、それぞれのハウジングがカバー部材21R,21Lに各々取付けられている。なお、図示は省略するが、電動モータ31R,31Lの電源(蓄電器)は、支柱フレーム13等、基体9の適所に搭載されている。
回転部材27Rは、左右方向の軸心を有する支軸33Rを介してカバー部材21Rに回転可能に支持されている。同様に、回転部材27Lは、左右方向の軸心を有する支軸33Lを介してカバー部材21Lに回転可能に支持されている。この場合、回転部材27Rの回転軸心(支軸33Rの軸心)と、回転部材27Lの回転軸心(支軸33Lの軸心)とは同軸心である。
回転部材27R,27Lは、それぞれ電動モータ31R,31Lの出力軸に、減速機としての機能を含む動力伝達機構を介して接続されており、電動モータ31R,31Lからそれぞれ伝達される動力(トルク)によって回転駆動される。各動力伝達機構は、例えばプーリ・ベルト式のものである。すなわち、図3に示す如く、回転部材27Rは、プーリ35Rとベルト37Rとを介して電動モータ31Rの出力軸に接続されている。同様に、回転部材27Lは、プーリ35Lとベルト37Lとを介して電動モータ31Lの出力軸に接続されている。
なお、上記動力伝達機構は、例えば、スプロケットとリンクチェーンとにより構成されるもの、あるいは、複数のギヤにより構成されるものであってもよい。また、例えば、電動モータ31R,31Lを、それぞれの出力軸が各回転部材27R,27Lと同軸心になるように各回転部材27R,27Lに対向させて配置し、電動モータ31R,31Lのそれぞれの出力軸を回転部材27R,27Lに各々、減速機(遊星歯車装置等)を介して連結するようにしてもよい。
各回転部材27R,27Lは、車輪体5側に向かって縮径する円錐台と同様の形状に形成されており、その外周面がテーパ外周面39R,39Lとなっている。
回転部材27Rのテーパ外周面39Rの周囲には、回転部材27Rと同心の円周上に等間隔で並ぶようにして、複数のフリーローラ29Rが配列されている。そして、これらのフリーローラ29Rは、それぞれ、ブラケット41Rを介してテーパ外周面39Rに取付けられ、該ブラケット41Rに回転自在に支承されている。
同様に、回転部材27Lのテーパ外周面39Lの周囲には、回転部材27Lと同心の円周上に等間隔で並ぶようにして、複数(フリーローラ29Rと同数)のフリーローラ29Lが配列されている。そして、これらのフリーローラ29Lは、それぞれ、ブラケット41Lを介してにテーパ外周面39Lに取付けられ、該ブラケット41Lに回転自在に支承されている。
前記車輪体5は、回転部材27R側のフリーローラ29Rと、回転部材27L側のフリーローラ29Lとの間に挟まれるようにして、回転部材27R,27Lと同軸心に配置されている。
この場合、図1及び図6に示すように、各フリーローラ29R,29Lは、その軸心C3が車輪体5の軸心C2に対して傾斜すると共に、車輪体5の直径方向(車輪体5をその軸心C2の方向で見たときに、該軸心C2と各フリーローラ29R,29Lとを結ぶ径方向)に対して傾斜する姿勢で配置されている。そして、このような姿勢で、各フリーローラ29R,29Lのそれぞれの外周面が車輪体5の内周面に斜め方向に圧接されている。
より一般的に言えば、右側のフリーローラ29Rは、回転部材27Rが軸心C2の周りに回転駆動されたときに、車輪体5との接触面で、軸心C2周りの方向の摩擦力成分(車輪体5の内周の接線方向の摩擦力成分)と、車輪体5の前記横断面中心C1の周り方向の摩擦力成分(円形の横断面の接線方向の摩擦力成分)とを車輪体5に作用させ得るような姿勢で、車輪体5の内周面に圧接されている。左側のフリーローラ29Lについても同様である。
この場合、前記したように、カバー部材21R,21Lは、図示しないバネによって、カバー部材21R,21Lの下端部側(二股の先端側)が狭まる方向に付勢されている。このため、この付勢力によって、右側のフリーローラ29Rと左側のフリーローラ29Lとの間に車輪体5が挟持されると共に、車輪体5に対する各フリーローラ29R,29Lの圧接状態(より詳しくはフリーローラ29R,29Lと車輪体5との間で摩擦力が作用し得る圧接状態)が維持される。
以上説明した構造を有する車両1においては、電動モータ31R,31Lによりそれぞれ、回転部材27R,27Lを同方向に等速度で回転駆動した場合には、車輪体5が回転部材27R,27Lと同方向に軸心C2の周りに回転することとなる。これにより、車輪体5が床面上を前後方向に輪転して、車両1の全体が前後方向に移動することとなる。なお、この場合は、車輪体5は、その横断面中心C1の周りには回転しない。
また、例えば、回転部材27R,27Lを互いに逆方向に同じ大きさの速度で回転駆動した場合には、車輪体5は、その横断面中心C1の周りに回転することとなる。これにより、車輪体5がその軸心C2の方向(すなわち左右方向)に移動し、ひいては、車両1の全体が左右方向に移動することとなる。なお、この場合は、車輪体5は、その軸心C2の周りには回転しない。
さらに、回転部材27R,27Lを、互いに異なる速度(方向を含めた速度)で、同方向又は逆方向に回転駆動した場合には、車輪体5は、その軸心C2の周りに回転すると同時に、その横断面中心C1の周りに回転することとなる。
この時、これらの回転動作の複合動作(合成動作)によって、前後方向及び左右方向に対して傾斜した方向に車輪体5が移動し、ひいては、車両1の全体が車輪体5と同方向に移動することとなる。この場合の車輪体5の移動方向は、回転部材27R,27Lの回転方向を含めた回転速度(回転方向に応じて極性が定義された回転速度ベクトル)の差に依存して変化するものとなる。
以上のように車輪体5の移動動作が行なわれるので、電動モータ31R,31Lのそれぞれの回転速度(回転方向を含む)を制御し、ひいては回転部材27R,27Lの回転速度を制御することによって、車両1の移動速度及び移動方向を制御できることとなる。
次に、本実施形態の車両1の動作制御のための構成を説明する。なお、以降の説明では、図1及び図2に示すように、前後方向の水平軸をX軸、左右方向の水平軸をY軸、鉛直方向をZ軸とするXYZ座標系を想定し、前後方向、左右方向をそれぞれX軸方向、Y軸方向と言うことがある。
まず、車両1の概略的な動作制御を説明すると、本実施形態では、基本的には、シート3に着座した乗員がその上体を傾けた場合(詳しくは、乗員と車両1とを合わせた全体の重心点の位置(水平面に投影した位置)を動かすように上体を傾けた場合)に、該上体を傾けた側に基体9がシート3と共に傾動する。そして、この時、基体9が傾いた側に車両1が移動するように、車輪体5の移動動作が制御される。例えば、乗員が上体を前傾させ、ひいては、基体9をシート3と共に前傾させると、車両1が前方に移動するように、車輪体5の移動動作が制御される。
すなわち、本実施形態では、乗員が上体を動かし、ひいては、シート3と共に基体9を傾動させるという動作が、車両1に対する1つの基本的な操縦操作(車両1の動作要求)とされ、その操縦操作に応じて車輪体5の移動動作がアクチュエータ装置7を介して制御される。
ここで、本実施形態の車両1は、その全体の接地面としての車輪体5の接地面が、車両1とこれに搭乗する乗員との全体を床面に投影した領域に比して面積が小さい単一の局所領域となり、その単一の局所領域だけに床反力が作用する。このため、基体9が傾倒しないようにするためには、乗員及び車両1の全体の重心点が車輪体5の接地面のほぼ真上に位置するように、車輪体5を動かす必要がある。
そこで、本実施形態では、乗員及び車両1の全体の重心点が、車輪体5の中心点(軸心C2上の中心点)のほぼ真上に位置する状態(より正確には当該重心点が車輪体5の接地面のほぼ真上に位置する状態)での基体9の姿勢を目標姿勢とし、基本的には、基体9の実際の姿勢を目標姿勢に収束させるように、車輪体5の移動動作が制御される。
また、車両1を発進させる場合等において、アクチュエータ装置7による推進力とは別に、例えば乗員が必要に応じて自身の足により床を蹴り、それにより車両1の移動速度を増速させる推進力(乗員の足平と床との摩擦力による推進力)を、付加的な外力として車両1に作用させた場合には、それに応じて車両1の移動速度(より正確には、乗員及び車両の全体の重心点の移動速度)が増速するように、車輪体5の移動動作が制御される。なお、当該推進力の付加が停止された状態では、車両1の移動速度が一旦、一定速度に保持された後、減衰して、該車両1が停止するように、車輪体5の移動動作が制御される(車輪体5の制動制御が行なわれる)。
さらに、車両1に乗員が搭乗していない状態では、車両1の単体の重心点が、車輪体5の中心点(軸心C2上の中心点)のほぼ真上に位置する状態(より正確には当該重心点が車輪体5の接地面のほぼ真上に位置する状態)での基体9の姿勢を目標姿勢とし、該基体9の実際の姿勢を目標姿勢に収束させ、ひいては、基体9が傾倒することなく車両1が自立するように、車輪体5の移動動作が制御される。
本実施形態では、以上の如き車両1の動作制御を行なうために、図1及び図2に示すように、マイクロコンピュータや電動モータ31R,31Lのドライブ回路ユニットなどを含む電子回路ユニットにより構成された制御ユニット50と、基体9の所定の部位の鉛直方向(重力方向)に対する傾斜角θb及びその変化速度(=dθb/dt)を計測するための傾斜センサ52と、車両1に乗員が搭乗しているか否かを検知するための荷重センサ54と、電動モータ31R,31Lのそれぞれの出力軸の回転角度及び回転角速度を検出するための角度センサとしてのロータリーエンコーダ56R,56Lがそれぞれ、車両1の適所に搭載されている。
この場合、制御ユニット50及び傾斜センサ52は、例えば、基体9の支柱フレーム13の内部に収容された状態で該支柱フレーム13に取付けられている。また、荷重センサ54は、シート3に内蔵されている。また、ロータリーエンコーダ56R,56Lは、それぞれ、電動モータ31R,31Lと一体に設けられている。なお、ロータリーエンコーダ56R,56Lは、それぞれ、回転部材27R,27Lに装着してもよい。
上記傾斜センサ52は、より詳しくは、加速度センサとジャイロセンサ等のレートセンサ(角速度センサ)とから構成され、これらのセンサの検出信号を制御ユニット50に出力する。そして、制御ユニット50が、傾斜センサ52の加速度センサ及びレートセンサの出力を基に、所定の計測演算処理(これは公知の演算処理でよい)を実行することによって、傾斜センサ52を搭載した部位(本実施形態では支柱フレーム13)の、鉛直方向に対する傾斜角度θbの計測値とその変化速度(微分値)である傾斜角速度θbdotの計測値とを算出する。
この場合、計測する傾斜角度θb(以降、基体傾斜角度θbということがある)は、より詳しくは、それぞれ、Y軸周り方向(ピッチ方向)の成分θb_xと、X軸周り方向(ロール方向)の成分θb_yとから成る。同様に、計測する傾斜角速度θbdot(以降、基体傾斜角速度θbdotということがある)も、Y軸周り方向(ピッチ方向)の成分θbdot_x(=dθb_x/dt)と、X軸周り方向(ロール方向)の成分θbdot_y(=dθb_y/dt)とから成る。
なお、本実施形態の説明では、上記基体傾斜角度θbなど、X軸及びY軸の各方向(又は各軸周り方向)の成分を有する運動状態量等の変数、あるいは、該運動状態量に関連する係数等の変数に関しては、その各成分を区別して表記する場合に、該変数の参照符号に、添え字“_x”又は“_y”を付加する。
この場合において、並進速度等の並進運動に係わる変数については、そのX軸方向の成分に添え字“_x”を付加し、Y軸方向の成分に添え字“_y”を付加する。
一方、角度、回転速度(角速度)、角加速度など、回転運動に係わる変数については、並進運動に係わる変数と添え字を揃えるために、便宜上、Y軸周り方向の成分に添え字“_x”を付加し、X軸周り方向の成分に添え字“_y”を付加する。
さらに、X軸方向の成分(又はY軸周り方向の成分)と、Y軸方向の成分(又はX軸周り方向の成分)との組として変数を表記する場合には、該変数の参照符号に添え字“_xy”を付加する。例えば、上記基体傾斜角度θbを、Y軸周り方向の成分θb_xとX軸周り方向の成分θb_yの組として表現する場合には、「基体傾斜角度θb_xy」というように表記する。
前記荷重センサ54は、乗員がシート3に着座した場合に該乗員の重量による荷重を受けるようにシート3に内蔵され、その荷重に応じた検出信号を制御ユニット50に出力する。そして、制御ユニット50が、この荷重センサ54の出力により示される荷重の計測値に基づいて、車両1に乗員が搭乗しているか否かを判断する。
なお、荷重センサ54の代わりに、例えば、乗員がシート3に着座したときにONとなるようなスイッチ式のセンサを用いてもよい。
ロータリーエンコーダ56Rは、電動モータ31Rの出力軸が所定角度回転する毎にパルス信号を発生し、このパルス信号を制御ユニット50に出力する。そして、制御ユニット50が、そのパルス信号を基に、電動モータ53Rの出力軸の回転角度を計測し、さらにその回転角度の計測値の時間的変化率(微分値)を電動モータ53Rの回転角速度として計測する。電動モータ31L側のロータリーエンコーダ56Lについても同様である。
制御ユニット50は、上記の各計測値を用いて所定の演算処理を実行することによって、電動モータ31R,31Lのそれぞれの回転角速度の目標値である速度指令を決定し、その速度指令に従って、電動モータ31R,31Lのそれぞれの回転角速度をフィードバック制御する。
なお、電動モータ31Rの出力軸の回転角速度と、回転部材27Rの回転角速度との間の関係は、該出力軸と回転部材27Rとの間の一定値の減速比に応じた比例関係になるので、本実施形態の説明では、便宜上、電動モータ31Rの回転角速度は、回転部材27Rの回転角速度を意味するものとする。同様に、電動モータ31Lの回転角速度は、回転部材27Lの回転角速度を意味するものとする。
以下に、制御ユニット50の制御処理をさらに詳細に説明する。
制御ユニット50は、所定の制御処理周期で図7のフローチャートに示す処理(メインルーチン処理)を実行する。
まず、ステップS1において、制御ユニット50は、傾斜センサ52の出力を取得する。
次いで、ステップS2に進んで、制御ユニット50は、取得した傾斜センサ52の出力を基に、基体傾斜角度θbの計測値θb_xy_sと、基体傾斜角速度θbdotの計測値θbdot_xy_sとを算出する。
なお、以降の説明では、上記計測値θb_xy_sなど、変数(状態量)の実際の値の観測値(計測値又は推定値)を参照符号により表記する場合に、該変数の参照符号に、添え字“_s”を付加する。
次いで、制御ユニット50は、ステップS3において、荷重センサ54の出力を取得した後、ステップS4の判断処理を実行する。この判断処理においては、制御ユニット50は、取得した荷重センサ54の出力が示す荷重計測値があらかじめ設定された所定値よりも大きいか否かによって、車両1に乗員が搭乗しているか否か(シート3に乗員が着座しているか否か)を判断する。
そして、制御ユニット50は、ステップS4の判断結果が肯定的である場合には、基体傾斜角度θbの目標値θb_xy_objを設定する処理と、車両1の動作制御用の定数パラメータ(各種ゲインの基本値など)の値を設定する処理とを、それぞれステップS5、6で実行する。
ステップS5においては、制御ユニット50は、基体傾斜角度θbの目標値θb_xy_objとして、あらかじめ定められた搭乗モード用の目標値を設定する。
ここで、「搭乗モード」は、車両1に乗員が搭乗している場合での車両1の動作モードを意味する。この搭乗モード用の目標値θb_xy_objは、車両1とシート3に着座した乗員との全体の重心点(以降、車両・乗員全体重心点という)が車輪体5の接地面のほぼ真上に位置する状態となる基体9の姿勢において、傾斜センサ52の出力に基づき計測される基体傾斜角度θbの計測値θb_xy_sに一致又はほぼ一致するようにあらかじめ設定されている。
また、ステップS6においては、制御ユニット50は、車両1の動作制御用の定数パラメータの値として、あらかじめ定められた搭乗モード用の値を設定する。なお、定数パラメータは、後述するhx,hy,Ki_a_x,Ki_b_x,Ki_a_y,Ki_b_y(i=1,2,3)等である。
一方、ステップS4の判断結果が否定的である場合には、制御ユニット50は、基体傾斜角度θb_xyの目標値θb_xy_objを設定する処理と、車両1の動作制御用の定数パラメータの値を設定する処理とを、ステップS7、8で実行する。
ステップS7においては、制御ユニット50は、傾斜角度θbの目標値θb_xy_objとして、あらかじめ定められた自立モード用の目標値を設定する。
ここで、「自立モード」は、車両1に乗員が搭乗していない場合での車両1の動作モードを意味する。この自立モード用の目標値θb_xy_objは、車両1単体の重心点(以降、車両単体重心点という)が車輪体5の接地面のほぼ真上に位置する状態となる基体9の姿勢において、傾斜センサ52の出力に基づき計測される基体傾斜角度θbの計測値θb_xy_sに一致又はほぼ一致するようにあらかじめ設定されている。この自立モード用の目標値θb_xy_objは、搭乗モード用の目標値θb_xy_objと一般的には異なる。
また、ステップS8においては、制御ユニット50は、車両1の動作制御用の定数パラメータの値として、あらかじめ定められた自立モード用の値を設定する。この自立モード用の定数パラメータの値は、搭乗モード用の定数パラメータの値と異なる。
搭乗モードと自立モードとで、上記定数パラメータの値を異ならせるのは、それぞれのモードで上記重心点の高さや、全体質量等が異なることに起因して、制御入力に対する車両1の動作の応答特性が互いに異なるからである。
以上のステップS4〜8の処理によって、搭乗モード及び自立モードの動作モード毎に個別に、基体傾斜角度θb_xyの目標値θb_xy_objと定数パラメータの値とが設定される。
なお、ステップS5,6の処理、又はステップS7,8の処理は、制御処理周期毎に実行することは必須ではなく、ステップS4の判断結果が変化した場合にだけ実行するようにしてもよい。
補足すると、搭乗モード及び自立モードのいずれにおいても、基体傾斜角速度θbdotのY軸周り方向の成分θbdot_xの目標値とX軸周り方向の成分θbdot_yの目標値とは、いずれも“0”である。このため、基体傾斜角速度θbdot_xyの目標値を設定する処理は不要である。
以上の如くステップS5,6の処理、又はステップS7,8の処理を実行した後、制御ユニット50は、次にステップS9において、車両制御演算処理を実行することによって、電動モータ31R,31Lのそれぞれの速度指令を決定する。この車両制御演算処理の詳細は後述する。
次いで、ステップS10に進んで、制御ユニット50は、ステップS9で決定した速度指令に応じて電動モータ31R,31Lの動作制御処理を実行する。この動作制御処理では、制御ユニット50は、ステップS9で決定した電動モータ31Rの速度指令と、ロータリーエンコーダ56Rの出力に基づき計測した電動モータ31Rの回転速度の計測値との偏差に応じて、該偏差を“0”に収束させるように電動モータ31Rの出力トルクの目標値(目標トルク)を決定する。そして、制御ユニット50は、その目標トルクの出力トルクを電動モータ31Rに出力させるように該電動モータ31Rの通電電流を制御する。左側の電動モータ31Lの動作制御についても同様である。
以上が、制御ユニット50が実行する全体的な制御処理である。
次に、上記ステップS9の車両制御演算処理の詳細を説明する。
なお、以降の説明においては、前記搭乗モードにおける車両・乗員全体重心点と、前記自立モードにおける車両単体重心点とを総称的に、車両系重心点という。該車両系重心点は、車両1の動作モードが搭乗モードである場合には、車両・乗員全体重心点を意味し、自立モードである場合には、車両単体重心点を意味する。
また、以降の説明では、制御ユニット50が各制御処理周期で決定する値(更新する値)に関し、現在の(最新の)制御処理周期で決定する値を今回値、その1つ前の制御処理周期で決定した値を前回値ということがある。そして、今回値、前回値を特にことわらない値は、今回値を意味する。
また、X軸方向の速度及び加速度に関しては、前方向きを正の向きとし、Y軸方向の速度及び加速度に関しては、左向きを正の向きとする。
本実施形態では、前記車両系重心点の動力学的な挙動(詳しくは、Y軸方向からこれに直交する面(XZ平面)に投影して見た挙動と、X軸方向からこれに直交する面(YZ平面)に投影して見た挙動)が、近似的に、図8に示すような、倒立振子モデルの挙動(倒立振子の動力学的挙動)によって表現されるものとして、ステップS9の車両制御演算処理が行なわれる。
なお、図8において、括弧を付していない参照符号は、Y軸方向から見た倒立振子モデルに対応する参照符号であり、括弧付きの参照符号は、X軸方向から見た倒立振子モデルに対応する参照符号である。
この場合、Y軸方向から見た挙動を表現する倒立振子モデルは、車両系重心点に位置する質点60_xと、Y軸方向に平行な回転軸62a_xを有して床面上を輪転自在な仮想的な車輪62_x(以降、仮想車輪62_xという)とを備える。そして、質点60_xが、仮想車輪62_xの回転軸62a_xに直線状のロッド64_xを介して支持され、該回転軸62a_xを支点として該回転軸62a_xの周りに揺動自在とされている。
この倒立振子モデルでは、質点60_xの運動が、Y軸方向から見た車両系重心点の運動に相当する。また、鉛直方向に対するロッド64_xの傾斜角度θbe_xがY軸周り方向での基体傾斜角度計測値θb_x_sと基体傾斜角度目標値θb_x_objとの偏差θbe_x_s(=θb_x_s−θb_x_obj)に一致するものとされる。また、ロッド64_xの傾斜角度θbe_xの変化速度(=dθbe_x/dt)がY軸周り方向の基体傾斜角速度計測値θbdot_x_sに一致するものとされる。また、仮想車輪62_xの移動速度Vw_x(X軸方向の並進移動速度)は、車両1の車輪体5のX軸方向の移動速度に一致するものとされる。
同様に、X軸方向から見た挙動を表現する倒立振子モデル(図8の括弧付きの符号を参照)は、車両系重心点に位置する質点60_yと、X軸方向に平行な回転軸62a_yを有して床面上を輪転自在な仮想的な車輪62_y(以降、仮想車輪62_yという)とを備える。そして、質点60_yが、仮想車輪62_yの回転軸62a_yに直線状のロッド64_yを介して支持され、該回転軸62a_yを支点として該回転軸62a_yの周りに揺動自在とされている。
この倒立振子モデルでは、質点60_yの運動が、X軸方向から見た車両系重心点の運動に相当する。また、鉛直方向に対するロッド64_yの傾斜角度θbe_yがX軸周り方向での基体傾斜角度計測値θb_y_sと基体傾斜角度目標値θb_y_objとの偏差θbe_y_s(=θb_y_s−θb_y_obj)に一致するものとされる。また、ロッド64_yの傾斜角度θbe_yの変化速度(=dθbe_y/dt)がX軸周り方向の基体傾斜角速度計測値θbdot_y_sに一致するものとされる。また、仮想車輪62_yの移動速度Vw_y(Y軸方向の並進移動速度)は、車両1の車輪体5のY軸方向の移動速度に一致するものとされる。
なお、仮想車輪62_x,62_yは、それぞれ、あらかじめ定められた所定値Rw_x,Rw_yの半径を有するものとされる。
また、仮想車輪62_x,62_yのそれぞれの回転角速度ωw_x,ωw_yと、電動モータ31R,31Lのそれぞれの回転角速度ω_R,ω_L(より正確には、回転部材27R,27Lのそれぞれの回転角速度ω_R,ω_L)との間には、次式01a,01bの関係が成立するものとされる。
ωw_x=(ω_R+ω_L)/2 ……式01a
ωw_y=C・(ω_R−ω_L)/2 ……式01b
なお、式01bにおける“C”は、前記フリーローラ29R,29Lと車輪体5との間の機構的な関係や滑りに依存する所定値の係数である。
ここで、図8に示す倒立振子モデルの動力学は、次式03x,03yにより表現される。なお、式03xは、Y軸方向から見た倒立振子モデルの動力学を表現する式、式03yは、X軸方向から見た倒立振子モデルの動力学を表現する式である。
d2θbe_x/dt2=α_x・θbe_x+β_x・ωwdot_x ……式03x
d2θbe_y/dt2=α_y・θbe_y+β_y・ωwdot_y ……式03y
式03xにおけるωwdot_xは仮想車輪62_xの回転角加速度(回転角速度ωw_xの1階微分値)、α_xは、質点60_xの質量や高さh_xに依存する係数、β_xは、仮想車輪62_xのイナーシャ(慣性モーメント)や半径Rw_xに依存する係数である。式03yにおけるωwdot_y、α_y、β_yについても上記と同様である。
これらの式03x,03yから判るように、倒立振子の質点60_x,60_yの運動(ひいては車両系重心点の運動)は、それぞれ、仮想車輪62_xの回転角加速度ωwdot_x、仮想車輪62_yの回転角加速度ωwdot_yに依存して規定される。
そこで、本実施形態では、Y軸方向から見た車両系重心点の運動を制御するための操作量(制御入力)として、仮想車輪62_xの回転角加速度ωwdot_xを用いると共に、X軸方向から見た車両系重心点の運動を制御するための操作量(制御入力)として、仮想車輪62_yの回転角加速度ωwdot_yを用いる。
そして、ステップS9の車両制御演算処理を概略的に説明すると、制御ユニット50は、X軸方向で見た質点60_xの運動と、Y軸方向で見た質点60_yの運動とが、車両系重心点の所望の運動に対応する運動となるように、操作量としての上記回転角加速度ωwdot_x,ωwdot_yの指令値(目標値)である仮想車輪回転角加速度指令ωwdot_x_cmd,ωwdot_y_cmdを決定する。さらに、制御ユニット50は、仮想車輪回転角加速度指令ωwdot_x_cmd,ωwdot_y_cmdをそれぞれ積分してなる値を、仮想車輪62_x,62_yのそれぞれの回転角速度ωw_x,ωw_yの指令値(目標値)である仮想車輪回転角速度指令ωw_x_cmd,ωw_y_cmdとして決定する。
そして、制御ユニット50は、仮想車輪回転角速度指令ωw_x_cmdに対応する仮想車輪62_xの移動速度(=Rw_x・ωw_x_cmd)と、仮想車輪回転角速度指令ωw_y_cmdに対応する仮想車輪62_yの移動速度(=Rw_y・ωw_y_cmd)とを、それぞれ、車両1の車輪体5のX軸方向の目標移動速度、Y軸方向の目標移動速度とし、それらの目標移動速度を実現するように、電動モータ31R,31Lのそれぞれの速度指令ω_R_cmd,ω_L_cmdを決定する。
なお、後述する他の移動体との連携動作を考慮に入れなければ、本実施形態では、操作量(制御入力)としての上記仮想車輪回転角加速度指令ωwdot_x_cmd,ωwdot_y_cmdは、それぞれ、後述する式07x,07yに示す如く、3個の操作量成分を加え合わせることによって決定される。
制御ユニット50は、上記の如き、ステップS9の車両制御演算処理を実行するための機能として、図9のブロック図で示す機能を備えている。
なお、図9において、符号70,72,74,76,78,80,82,84が付された要素は、本実施形態による倒立振子型移動体の全方向移動車両としての基本的な移動制御に関わる要素であり、その他の符号701〜708が付された要素は、本実施形態による倒立振子の連携動作に関わる要素である。ここでは、図9の上記要素のうち、全方向移動車両としての基本的な移動制御に関わる要素のみについて説明するものとし、連携動作に関わる要素の詳細については、後述の「2.連携技術」の項で説明する。
制御ユニット50は、基体傾斜角度計測値θb_xy_sと基体傾斜角度目標値θb_xy_objとの偏差である基体傾斜角度偏差計測値θbe_xy_sを算出する偏差演算部70と、前記車両系重心点の移動速度である重心速度Vb_xyの観測値としての重心速度推定値Vb_xy_sを算出する重心速度算出部72と、乗員等による車両1の操縦操作(車両1に推進力を付加する操作)によって要求されていると推定される上記重心速度Vb_xyの要求値としての要求重心速度V_xy_aimを生成する要求重心速度生成部74と、これらの重心速度推定値Vb_xy_s及び要求重心速度V_xy_aimから、電動モータ31R,31Lの回転角速度の許容範囲に応じた制限を加味して、重心速度Vb_xyの目標値としての制御用目標重心速度Vb_xy_mdfdを決定する重心速度制限部76と、後述する式07x,07yのゲイン係数の値を調整するためのゲイン調整パラメータKr_xyを決定するゲイン調整部78とを備える。
制御ユニット50は、さらに、前記仮想車輪回転角速度指令ωw_xy_cmdを算出する姿勢制御演算部80と、この仮想車輪回転角速度指令ωw_xy_cmdを、右側の電動モータ31Rの速度指令ω_R_cmd(回転角速度の指令値)と左側の電動モータ31Lの速度指令ω_L_cmd(回転角速度の指令値)との組に変換するモータ指令演算部82とを備える。
なお、図9中の参照符号84を付したものは、姿勢制御演算部80が制御処理周期毎に算出する仮想車輪回転角速度指令ωw_xy_cmdを入力する遅延要素を示している。該遅延要素84は、各制御処理周期において、仮想車輪回転角速度指令ωw_xy_cmdの前回値ωw_xy_cmd_pを出力する。
前記ステップS9の車両制御演算処理では、これらの上記の各処理部の処理が以下に説明するように実行される。
すなわち、制御ユニット50は、まず、偏差演算部70の処理と重心速度算出部72の処理とを実行する。
偏差演算部70には、前記ステップS2で算出された基体傾斜角度計測値θb_xy_s(θb_x_s及びθb_y_s)と、前記ステップS5又はステップS7で設定された目標値θb_xy_obj(θb_x_obj及びθb_y_obj)とが入力される。そして、偏差演算部70は、θb_x_sからθb_x_objを減算することによって、Y軸周り方向の基体傾斜角度偏差計測値θbe_x_s(=θb_x_s−θb_x_obj)を算出すると共に、θb_y_sからθb_y_objを減算することによって、X軸周り方向の基体傾斜角度偏差計測値θbe_y_s(=θb_y_s−θb_y_obj)を算出する。
なお、偏差演算部70の処理は、ステップS9の車両制御演算処理の前に行なうようにしてもよい。例えば、前記ステップS5又は7の処理の中で、偏差演算部70の処理を実行してもよい。
前記重心速度算出部72には、前記ステップS2で算出された基体傾斜角速度計測値θbdot_xy_s(θbdot_x_s及びθbdot_y_s)の今回値が入力されると共に、仮想車輪速度指令ωw_xy_cmdの前回値ωw_xy_cmd_p(ωw_x_cmd_p及びωw_y_cmd_p)が遅延要素84から入力される。そして、重心速度算出部72は、これらの入力値から、前記倒立振子モデルに基づく所定の演算式によって、重心速度推定値Vb_xy_s(Vb_x_s及びVb_y_s)を算出する。
具体的には、重心速度算出部72は、次式05x,05yにより、Vb_x_s及びVb_y_sをそれぞれ算出する。
Vb_x_s=Rw_x・ωw_x_cmd_p+h_x・θbdot_x_s ……05x
Vb_y_s=Rw_y・ωw_y_cmd_p+h_y・θbdot_y_s ……05y
これらの式05x,05yにおいて、Rw_x,Rw_yは、前記したように、仮想車輪62_x,62_yのそれぞれの半径であり、これらの値は、あらかじめ設定された所定値である。また、h_x,h_yは、それぞれ倒立振子モデルの質点60_x,60_yの高さである。この場合、本実施形態では、車両系重心点の高さは、ほぼ一定に維持されるものとされる。そこで、h_x,h_yの値としては、それぞれ、あらかじめ設定された所定値が用いられる。補足すると、高さh_x,h_yは、前記ステップS6又は8において値を設定する定数パラメータに含まれるものである。
上記式05xの右辺の第1項は、仮想車輪62_xの速度指令の前回値ωw_x_cmd_pに対応する該仮想車輪62_xのX軸方向の移動速度であり、この移動速度は、車輪体5のX軸方向の実際の移動速度の現在値に相当するものである。また、式05xの右辺の第2項は、基体9がY軸周り方向にθbdot_x_sの傾斜角速度で傾動することに起因して生じる車両系重心点のX軸方向の移動速度(車輪体5に対する相対的な移動速度)の現在値に相当するものである。これらのことは、式05yについても同様である。
なお、前記ロータリーエンコーダ56R,56Lの出力を基に計測される電動モータ31R,31Lのそれぞれの回転角速度の計測値(今回値)の組を、仮想車輪62_x,62_yのそれぞれの回転角速度の組に変換し、それらの回転角速度を、式05x、05yのωw_x_cmd_p、ωw_y_cmd_pの代わりに用いてもよい。ただし、回転角速度の計測値に含まれるノイズの影響を排除する上では、目標値であるωw_x_cmd_p、ωw_y_cmd_pを使用することが有利である。
次に、制御ユニット50は、要求重心速度生成部74の処理とゲイン調整部78の処理とを実行する。この場合、要求重心速度生成部74及びゲイン調整部78には、それぞれ、重心速度算出部72で上記の如く算出された重心速度推定値Vb_xy_s(Vb_x_s及びVb_y_s)が入力される。
そして、要求重心速度生成部74は、詳細は後述するが、車両1の動作モードが搭乗モードである場合に、入力された重心速度推定値Vb_xy_s(Vb_x_s及びVb_y_s)を基に、要求重心速度V_xy_aim(V_x_aim,V_y_aim)を決定する。なお、本実施形態では、車両1の動作モードが自立モードである場合には、要求重心速度生成部74は、要求重心速度V_x_aim及びV_y_aimをいずれも“0”とする。
また、ゲイン調整部78は、入力された重心速度推定値Vb_xy_s(Vb_x_s及びVb_y_s)を基に、前記ゲイン調整パラメータKr_xy(Kr_x及びKr_y)を決定する。
このゲイン調整部78の処理を図10及び図11を参照して以下に説明する。
図10に示すように、ゲイン調整部78は、入力された重心速度推定値Vb_x_s,Vb_y_sをリミット処理部86に入力する。このリミット処理部86では、重心速度推定値Vb_x_s,Vb_y_sに、電動モータ31R,31Lのそれぞれの回転角速度の許容範囲に応じた制限を適宜、加えることによって、出力値Vw_x_lim1,Vw_y_lim1を生成する。出力値Vw_x_lim1は、前記仮想車輪62_xのX軸方向の移動速度Vw_xの制限後の値、出力値Vw_y_lim1は、前記仮想車輪62_yのY軸方向の移動速度Vw_yの制限後の値としての意味を持つ。
このリミット処理部86の処理を、図11を参照してさらに詳細に説明する。なお、図11中の括弧付きの参照符号は、後述する重心速度制限部76のリミット処理部104の処理を示すものであり、リミット処理部86の処理に関する説明では無視してよい。
リミット処理部86は、まず、重心速度推定値Vb_x_s,Vb_y_sをそれぞれ処理部86a_x,86a_yに入力する。処理部86a_xは、Vb_x_sを仮想車輪62_xの半径Rw_xで除算することによって、仮想車輪62_xのX軸方向の移動速度をVb_x_sに一致させたと仮定した場合の該仮想車輪62_xの回転角速度ωw_x_sを算出する。同様に、処理部86a_yは、仮想車輪62_yのY軸方向の移動速度をVb_y_sに一致させたと仮定した場合の該仮想車輪62_yの回転角速度ωw_y_s(=Vb_y_s/Rw_y)を算出する。
次いで、リミット処理部86は、ωw_x_s,ωw_y_sの組を、XY−RL変換部86bにより、電動モータ31Rの回転角速度ω_R_sと電動モータ31Lの回転角速度ω_L_sとの組に変換する。
この変換は、本実施形態では、前記式01a,01bのωw_x,ωw_y,ω_R,ω_Lをそれぞれ、ωw_x_s,ωw_y_s,ω_R_s,ω_L_sに置き換えて得られる連立方程式を、ω_R_s,ω_L_sを未知数として解くことにより行なわれる。
次いで、リミット処理部86は、XY−RL変換部86bの出力値ω_R_s,ω_L_sをそれぞれ、リミッタ86c_R,86c_Lに入力する。このとき、リミッタ86c_Rは、ω_R_sが、あらかじめ設定された所定値の上限値(>0)と下限値(<0)とを有する右モータ用許容範囲内に収まっている場合には、ω_R_sをそのまま出力値ω_R_lim1として出力する。また、リミッタ86c_Rは、ω_R_sが、右モータ用許容範囲から逸脱している場合には、該右モータ用許容範囲の上限値と下限値とのうちのω_R_sに近い方の境界値を出力値ω_R_lim1として出力する。これにより、リミッタ86c_Rの出力値ω_R_lim1は、右モータ用許容範囲内の値に制限される。
同様に、リミッタ86c_Lは、ω_L_sが、あらかじめ設定された所定値の上限値(>0)と下限値(<0)とを有する左モータ用許容範囲内に収まっている場合には、ω_L_sをそのまま出力値ω_L_lim1として出力する。また、リミッタ86c_Lは、ω_L_sが、左モータ用許容範囲から逸脱している場合には、該左モータ用許容範囲の上限値と下限値とのうちのω_L_sに近い方の境界値を出力値ω_L_lim1として出力する。これにより、リミッタ86c_Lの出力値ω_L_lim1は、左モータ用許容範囲内の値に制限される。
上記右モータ用許容範囲は右側の電動モータ31Rの回転角速度(絶対値)が高くなり過ぎないようにし、ひいては、電動モータ31Rが出力可能なトルクの最大値が低下するのを防止するために設定された許容範囲である。このことは、左モータ用許容範囲についても同様である。
次いで、リミット処理部86は、リミッタ86c_R,86c_Lのそれぞれの出力値ω_R_lim1,ω_L_lim1の組を、RL−XY変換部86dにより、仮想車輪62_x,62_yのそれぞれの回転角速度ωw_x_lim1,ωw_y_lim1の組に変換する。
この変換は、前記XY−RL変換部86bの変換処理の逆変換の処理である。この処理は、前記式01a,01bのωw_x,ωw_y,ω_R,ω_Lをそれぞれ、ωw_x_lim1,ωw_y_lim1,ω_R_lim1,ω_L_lim1に置き換えて得られる連立方程式を、ωw_x_lim1,ωw_y_lim1を未知数として解くことにより行なわれる。
次いで、リミット処理部86は、RL−XY変換部86dの出力値ωw_x_lim1,ωw_y_lim1をそれぞれ処理部86e_x,86e_yに入力する。処理部86e_xは、ωw_x_lim1に仮想車輪62_xの半径Rw_xを乗じることによって、ωw_x_lim1を仮想車輪62_xの移動速度Vw_x_lim1に変換する。同様に、処理部86e_yは、ωw_y_lim1を仮想車輪62_yの移動速度Vw_y_lim1(=ωw_y_lim1・Rw_y)に変換する。
以上のリミット処理部86の処理によって、仮想車輪62_xのX軸方向の移動速度Vw_xと、仮想車輪62_yのY軸方向の移動速度Vw_yとをそれぞれ重心速度推定値Vb_x_s,Vb_y_sに一致させたと仮定した場合(換言すれば、車輪体5のX軸方向の移動速度とY軸方向の移動速度とをそれぞれ、Vb_x_s,Vb_y_sに一致させたと仮定した場合)に、それらの移動速度を実現するために必要な電動モータ31R,31Lのそれぞれの回転角速度ω_R_s,ω_L_sが、両方とも、許容範囲内に収まっている場合には、Vb_x_s,Vb_y_sにそれぞれ一致する出力値Vw_x_lim1,Vw_y_lim1の組がリミット処理部86から出力される。
一方、電動モータ31R,31Lのそれぞれの回転角速度ω_R_s,ω_L_sの両方又は一方が許容範囲から逸脱している場合には、その両方又は一方の回転角速度が強制的に許容範囲内に制限された上で、その制限後の電動モータ31R,31Lのそれぞれの回転角速度ω_R_lim1,ω_L_lim1の組に対応する、X軸方向及びY軸方向の移動速度Vw_x_lim1,Vw_y_lim1の組がリミット処理部86から出力される。
従って、リミット処理部86は、その出力値Vw_x_lim1,Vw_y_lim1の組に対応する電動モータ31R,31Lのそれぞれの回転角速度が許容範囲を逸脱しないことを必須の必要条件として、その必要条件下で可能な限り、出力値Vw_x_lim1,Vw_y_lim1をそれぞれVb_x_s,Vb_y_sに一致させるように、出力値Vw_x_lim1,Vw_y_lim1の組を生成する。
図10の説明に戻って、ゲイン調整部78は、次に、演算部88_x,88_yの処理を実行する。演算部88_xには、X軸方向の重心速度推定値Vb_x_sと、リミット処理部86の出力値Vw_x_lim1とが入力される。そして、演算部88_xは、Vw_x_lim1からVb_x_sを減算してなる値Vover_xを算出して出力する。また、演算部88_yには、Y軸方向の重心速度推定値Vb_y_sと、リミット処理部86の出力値Vw_y_lim1とが入力される。そして、演算部88_yは、Vw_y_lim1からVb_y_sを減算してなる値Vover_yを算出して出力する。
この場合、リミット処理部86での出力値Vw_x_lim1,Vw_y_lim1の強制的な制限が行なわれなかった場合には、Vw_x_lim1=Vb_x_s、Vw_y_lim1=Vb_y_sとなるので、演算部88_x,88_yのそれぞれの出力値Vover_x,Vover_yはいずれも“0”となる。
一方、リミット処理部86の出力値Vw_x_lim1,Vw_y_lim1が、入力値Vb_x_s,Vb_y_sに対して強制的な制限を施して生成された場合には、Vw_x_lim1のVb_x_sからの修正量(=Vw_x_lim1−Vb_x_s)と、Vw_y_lim1のVb_y_sからの修正量(=Vw_y_lim1−Vb_y_s)とがそれぞれ、演算部88_x,88_yから出力される。
次いで、ゲイン調整部78は、演算部88_xの出力値Vover_xを処理部90_x,92_xに順番に通すことによって、ゲイン調整パラメータKr_xを決定する。また、ゲイン調整部78は、演算部88_yの出力値Vover_yを処理部90_y,92_yに順番に通すことによって、ゲイン調整パラメータKr_yを決定する。なお、ゲイン調整パラメータKr_x,Kr_yは、いずれも“0”から“1”までの範囲内の値である。
上記処理部90_xは、入力されるVover_xの絶対値を算出して出力する。また、処理部92_xは、その出力値Kr_xが入力値|Vover_x|に対して単調に増加し、且つ、飽和特性を有するようにKr_xを生成する。該飽和特性は、入力値がある程度大きくなると、入力値の増加に対する出力値の変化量が“0”になるか、もしくは、“0”に近づく特性である。
この場合、本実施形態では、処理部92_xは、入力値|Vover_x|があらかじめ設定された所定値以下である場合には、該入力値|Vover_x|に所定値の比例係数を乗じてなる値をKr_xとして出力する。また、処理部92_xは、入力値|Vover_x|が所定値よりも大きい場合には、“1”をKr_xとして出力する。なお、上記比例係数は、|Vover_x|が所定値に一致するときに、|Vover_x|と比例係数との積が“1”になるように設定されている。
また、処理部90_y,92_yの処理は、それぞれ上記した処理部90_x,92_xの処理と同様である。
以上説明したゲイン調整部78の処理によって、リミット処理部86での出力値Vw_x_lim1,Vw_y_lim1の強制的な制限が行なわれなかった場合、すなわち、車輪体5のX軸方向及びY軸方向のそれぞれの移動速度Vw_x,Vw_yを、それぞれ、重心速度推定値Vb_x_s,Vb_y_sに一致させるように電動モータ31R,31Lを動作させても、電動モータ31R,31Lのそれぞれの回転角速度が許容範囲内に収まるような場合には、ゲイン調整パラメータKr_x,Kr_yはいずれも“0”に決定される。
一方、リミット処理部86の出力値Vw_x_lim1,Vw_y_lim1が、入力値Vb_x_s,Vb_y_sに対して強制的な制限を施して生成された場合、すなわち、車輪体5のX軸方向及びY軸方向のそれぞれの移動速度Vw_x,Vw_yを、それぞれ、重心速度推定値Vb_x_s,Vb_y_sに一致させるように電動モータ31R,31Lを動作させると、電動モータ31R,31Lのいずれかの回転角速度が許容範囲を逸脱してしまう場合(いずれかの回転角速度の絶対値が高くなり過ぎる場合)には、前記修正量Vover_x,Vover_yのそれぞれの絶対値に応じて、ゲイン調整パラメータKr_x,Kr_yの値がそれぞれ決定される。この場合、Kr_xは、“1”を上限値して、修正量Vx_overの絶対値が大きいほど、大きな値になるように決定される。このことは、Kr_yについても同様である。
図9の説明に戻って、制御ユニット50は、重心速度算出部72及び要求重心速度生成部74の処理を前記した如く実行した後、次に、重心速度制限部76の処理を実行する。
この重心速度制限部76には、重心速度算出部72で算出された重心速度推定値Vb_xy_s(Vb_x_s及びVb_y_s)と、要求重心速度生成部74で決定された要求重心速度V_xy_aim(V_x_aim及びV_y_aim)とが入力される。そして、重心速度制限部76は、これらの入力値を使用して、図12のブロック図で示す処理を実行することによって、制御用目標重心速度V_xy_mdfd(V_x_mdfd及びV_y_mdfd)を決定する。
具体的には、重心速度制限部76は、まず、定常偏差算出部94_x,94_yの処理を実行する。
この場合、定常偏差算出部94_xには、X軸方向の重心速度推定値Vb_x_sが入力されると共に、X軸方向の制御用目標重心速度Vb_x_mdfdの前回値Vb_x_mdfd_pが遅延要素96_xを介して入力される。そして、定常偏差算出部94_xは、まず、入力されるVb_x_sが比例・微分補償要素(PD補償要素)94a_xに入力する。この比例・微分補償要素94_xは、その伝達関数が1+Kd・Sにより表される補償要素であり、入力されるVb_x_sと、その微分値(時間的変化率)に所定値の係数Kdを乗じてなる値とを加算し、その加算結果の値を出力する。
次いで、定常偏差算出部94_xは、入力されるVb_x_mdfd_pを、比例・微分補償要素94_xの出力値から減算してなる値を演算部94b_xにより算出した後、この演算部94b_xの出力値を、位相補償機能を有するローパスフィルタ94c_xに入力する。このローパスフィルタ94c_xは、伝達関数が(1+T2・S)/(1+T1・S)により表されるフィルタである。そして、定常偏差算出部94_xは、このローパスフィルタ94c_xの出力値Vb_x_prdを出力する。
また、定常偏差算出部94_yには、Y軸方向の重心速度推定値Vb_y_sが入力されると共に、Y軸方向の制御用目標重心速度Vb_y_mdfdの前回値Vb_y_mdfd_pが遅延要素96_yを介して入力される。
そして、定常偏差算出部94_yは、上記した定常偏差算出部94_xと同様に、比例・微分補償要素94a_y、演算部94b_y及びローパスフィルタ94c_yの処理を順次実行し、ローパスフィルタ94c_yの出力値Vb_y_prdを出力する。
ここで、定常偏差算出部94_xの出力値Vb_x_prdは、Y軸方向から見た車両系重心点の現在の運動状態(換言すればY軸方向から見た倒立振子モデルの質点60_xの運動状態)から推測される、将来のX軸方向の重心速度推定値の収束予測値の制御用目標重心速度Vb_x_mdfdに対する定常偏差としての意味を持つものである。同様に、定常偏差算出部94_y出力値Vb_y_prdは、X軸方向から見た車両系重心点の現在の運動状態(換言すればX軸方向から見た倒立振子モデルの質点60_yの運動状態)から推測される、将来のY軸方向の重心速度推定値の収束予測値の制御用目標重心速度Vb_y_mdfdに対する定常偏差としての意味を持つものである。以降、定常偏差算出部94_x,94_yのそれぞれの出力値Vb_x_prd,Vb_y_prdを重心速度定常偏差予測値という。
重心速度制限部76は、上記の如く定常偏差算出部94_x,94_yの処理を実行した後、定常偏差算出部94_xの出力値Vb_x_prdに要求重心速度Vb_x_aimを加算する処理と、定常偏差算出部94_yの出力値Vb_y_prdに要求重心速度Vb_y_aimを加算する処理とをそれぞれ、演算部98_x,98_yにより実行する。
従って、演算部98_xの出力値Vb_x_tは、X軸方向の重心速度定常偏差予測値Vb_x_prdに、X軸方向の要求重心速度Vb_x_aimを付加した速度となる。同様に、演算部98_yの出力値Vb_y_tは、Y軸方向の重心速度定常偏差予測値Vb_y_prdに、Y軸方向の要求重心速度Vb_y_aimを付加した速度となる。
なお、車両1の動作モードが自立モードである場合等、X軸方向の要求重心速度Vb_x_aimが“0”である場合には、X軸方向の重心速度定常偏差予測値Vb_x_prdがそのまま、演算部98_xの出力値Vb_x_tとなる。同様に、Y軸方向の要求重心速度Vb_y_aimが“0”である場合には、Y軸方向の重心速度定常偏差予測値Vb_y_prdがそのまま、演算部98_yの出力値Vb_y_tとなる。
次いで、重心速度制限部76は、演算部98_x,98_yのそれぞれの出力値Vb_x_t,Vb_y_tを、リミット処理部100に入力する。このリミット処理部100の処理は、前記したゲイン調整部78のリミット処理部86の処理と同じである。この場合、図11に括弧付きに参照符号で示す如く、リミット処理部100の各処理部の入力値及び出力値だけがリミット処理部86と相違する。
具体的には、リミット処理部100では、前記仮想車輪62_x,62_yのそれぞれの移動速度Vw_x,Vw_yを、Vb_x_t,Vb_y_tにそれぞれ一致させたと仮定した場合の各仮想車輪62_x,62_yの回転角速度ωw_x_t,ωw_y_tがそれぞれ処理部86a_x,86a_yにより算出される。そして、この回転角速度ωw_x_t,ωw_y_tの組が、XY−RL変換部86bにより、電動モータ31R,31Lの回転角速度ω_R_t,ω_L_tの組に変換される。
さらに、これらの回転角速度ω_R_t,ω_L_tが、リミッタ86c_R,86c_Lによって、それぞれ、右モータ用許容範囲内の値と左モータ用許容範囲内の値とに制限される。そして、この制限処理後の値ω_R_lim2,ω_L_lim2が、RL−XY変換部86dによって、仮想車輪62_x,62_yの回転角速度ωw_x_lim2,ωw_y_lim2に変換される。
次いで、この各回転角速度ωw_x_lim2,ωw_y_lim2に対応する各仮想車輪62_x,62_yの移動速度Vw_x_lim2,Vw_y_lim2がそれぞれ処理部86e_x,86e_yによって算出され、これらの移動速度Vw_x_lim2,Vw_y_lim2がリミット処理部100から出力される。
以上のリミット処理部100の処理によって、リミット処理部100は、リミット処理部86と同様に、その出力値Vw_x_lim2,Vw_y_lim2の組に対応する電動モータ31R,31Lのそれぞれの回転角速度が許容範囲を逸脱しないことを必須の必要条件として、その必要条件下で可能な限り、出力値Vw_x_lim2,Vw_y_lim2をそれぞれVb_x_t,Vb_y_tに一致させるように、出力値Vw_x_lim2,Vw_y_lim2の組を生成する。
なお、リミット処理部100における右モータ用及び左モータ用の各許容範囲は、リミット処理部86における各許容範囲と同一である必要はなく、互いに異なる許容範囲に設定されていてもよい。
図12の説明に戻って、重心速度制限部76は、次に、演算部102_x,102_yの処理を実行することによって、それぞれ制御用目標重心速度Vb_x_mdfd,Vb_y_mdfdを算出する。この場合、演算部102_xは、リミット処理部100の出力値Vw_x_lim2から、X軸方向の重心速度定常偏差予測値Vb_x_prdを減算してなる値をX軸方向の制御用目標重心速度Vb_x_mdfdとして算出する。同様に、演算部102_yは、リミット処理部100の出力値Vw_y_lim2から、Y軸方向の重心速度定常偏差予測値Vb_y_prdを減算してなる値をY軸方向の制御用目標重心速度Vb_y_mdfdとして算出する。
以上のようにして決定される制御用目標重心速度Vb_x_mdfd,Vb_y_mdfdは、リミット処理部100での出力値V_x_lim2,V_y_lim2の強制的な制限が行なわれなかった場合、すなわち、車輪体5のX軸方向及びY軸方向のそれぞれの移動速度を、それぞれ、演算部98_xの出力値Vb_x_tと演算部98_yの出力値Vb_y_tとに一致させるように電動モータ31R,31Lを動作させても、電動モータ31R,31Lのそれぞれの回転角速度が許容範囲内に収まるような場合には、要求重心速度Vb_x_aim,Vb_y_aimがそれぞれ、そのまま、制御用目標重心速度Vb_x_mdfd,Vb_y_mdfdとして決定される。
なお、この場合、X軸方向の要求重心速度Vb_x_aimが“0”であれば、X軸方向の制御用目標重心速度Vb_x_mdfdも“0”となり、Y軸方向の要求重心速度Vb_y_aimが“0”であれば、Y軸方向の制御用目標重心速度Vb_y_mdfdも“0”となる。
一方、リミット処理部100の出力値Vw_x_lim2,Vw_y_lim2が、入力値Vb_x_t,Vb_y_tに対して強制的な制限を施して生成された場合、すなわち、車輪体5のX軸方向及びY軸方向のそれぞれの移動速度を、それぞれ、演算部98_xの出力値Vb_x_tと演算部98_yの出力値Vb_y_tとに一致させるように電動モータ31R,31Lを動作させると、電動モータ31R,31Lのいずれかの回転角速度が許容範囲を逸脱してしまう場合(いずれかの回転角速度の絶対値が高くなり過ぎる場合)には、X軸方向については、リミット処理部100の出力値Vw_x_lim2の入力値Vb_x_tからの修正量(=Vw_x_lim2−Vb_x_t)だけ、要求重心速度Vb_x_aimを補正してなる値(当該修正量をVb_x_aimに加算した値)が、X軸方向の制御用目標重心速度Vb_x_mdfdとして決定される。
また、Y軸方向については、リミット処理部100の出力値Vw_y_lim2の入力値Vb_y_tからの修正量(=Vw_y_lim2−Vb_y_t)だけ、要求重心速度Vb_y_aimを補正してなる値(当該修正量をVb_y_aimに加算した値)が、Y軸方向の制御用目標重心速度Vb_y_mdfdとして決定される。
この場合において、例えばX軸方向の速度に関し、要求重心速度Vb_x_aimが“0”でない場合には、制御用目標重心速度Vb_x_mdfdは、要求重心速度Vb_x_aimよりも“0”に近づくか、もしくは、要求重心速度Vb_x_aimと逆向きの速度となる。また、要求重心速度Vb_x_aimが“0”である場合には、制御用目標重心速度Vb_x_mdfdは、定常偏差算出部94_xが出力するX軸方向の重心速度定常偏差予測値Vb_x_prdと逆向きの速度となる。これらのことは、Y軸方向の速度に関しても同様である。
以上が、重心速度制限部76の処理である。
図9の説明に戻って、制御ユニット50は、以上の如く重心速度算出部72、重心速度制限部76、ゲイン調整部78、及び偏差演算部70の処理を実行した後、次に、姿勢制御演算部80の処理を実行する。
この姿勢制御演算部80の処理を、以下に図13を参照して説明する。
なお、図13において、符号80a〜80fが付された要素は、本実施形態による倒立振子型移動体の全方向移動車両としての基本的な移動制御に関わる要素であり、その他の符号80g〜80jが付された要素は、本実施形態による倒立振子型移動体の連携動作に関わる要素である。ここでは、図13の上記要素のうち、全方向移動車両としての基本的な移動制御に関わる要素のみについて説明するものとし、連携動作に関わる要素の詳細については、後述の「2.連携技術」の項で説明する。
また、図13において、括弧を付していない参照符号は、X軸方向に輪転する仮想車輪62_xの回転角速度の目標値である前記仮想車輪回転角速度指令ωw_x_comを決定する処理に係わる参照符号であり、括弧付きの参照符合は、Y軸方向に輪転する仮想車輪62_yの回転角速度の目標値である前記仮想車輪回転角速度指令ωw_y_comを決定する処理に係わる参照符号である。
姿勢制御演算部80には、偏差演算部70で算出された基体傾斜角度偏差計測値θbe_xy_sと、前記ステップS2で算出された基体傾斜角速度計測値θbdot_xy_sと、重心速度算出部72で算出された重心速度推定値Vb_xy_sと、重心速度制限部76で算出された目標重心速度Vb_xy_cmdと、ゲイン調整部78で算出されたゲイン調整パラメータKr_xyとが入力される。
なお、姿勢制御演算部80には、上記の入力値のほか、連携動作を制御するための入力値として、図9に示す後述の偏差演算部702により算出された基体傾斜角度計測値の偏差(θb_xy_s(R)−θb_xy_s(L))と、偏差演算部703により算出された基体傾斜角速度計測値の偏差(θbdot_xy_s(R)−θbdot_xy_s(L))と、偏差演算部705により算出された距離の偏差(Ls−Ltarget)と、外力センサ706により検知された外力Fhandの値と、ゲイン調整部708により生成されたゲインKw_xyが姿勢制御演算部80に与えられるが、その詳細については「2.連携技術」の項で説明するものとし、ここでは、説明の便宜上、連携動作を制御するための入力値を考慮しないで説明することとする。
姿勢制御演算部80は、まず、これらの入力値(連携動作を制御するための入力値を除く)を用いて、次式07x,07yにより、仮想車輪回転角加速度指令ωdotw_xy_comを算出する。
ωwdot_x_cmd=K1_x・θbe_x_s+K2_x・θbdot_x_s
+K3_x・(Vb_x_s−Vb_x_mdfd) ……式07x
ωwdot_y_cmd=K1_y・θbe_y_s+K2_y・θbdot_y_s
+K3_y・(Vb_y_s−Vb_y_mdfd) ……式07y
従って、本実施形態では、Y軸方向から見た倒立振子モデルの質点60_xの運動(ひいては、Y軸方向から見た車両系重心点の運動)を制御するための操作量(制御入力)としての仮想車輪回転角加速度指令ωdotw_x_comと、X軸方向から見た倒立振子モデルの質点60_yの運動(ひいては、X軸方向から見た車両系重心点の運動)を制御するための操作量(制御入力)としての仮想車輪回転角加速度指令ωdotw_y_comとは、それぞれ、3つの操作量成分(式07x,07yの右辺の3つの項)を加え合わせることによって決定される。
この場合、式07xにおける各操作量成分に係わるゲイン係数K1_x,K2_x,K3_xは、ゲイン調整パラメータKr_xに応じて可変的に設定され、式07yにおける各操作量成分に係わるゲイン係数K1_y,K2_y,K3_yは、ゲイン調整パラメータKr_yに応じて可変的に設定される。以降、式07xにおけるゲイン係数K1_x,K2_x,K3_xのそれぞれを第1ゲイン係数K1_x、第2ゲイン係数K2_x、第3ゲイン係数K3_xということがある。このことは、式07yにおけるゲイン係数K1_y,K2_y,K3_yについても同様とする。
式07xにおける第iゲイン係数Ki_x(i=1,2,3)と、式07yにおける第iゲイン係数Ki_y(i=1,2,3)とは、図13中にただし書きで示した如く、次式09x、09yにより、ゲイン調整パラメータKr_x,Kr_yに応じて決定される。
Ki_x=(1−Kr_x)・Ki_a_x+Kr_x・Ki_b_x ……式09x
Ki_y=(1−Kr_y)・Ki_a_y+Kr_y・Ki_b_y ……式09y
(i=1,2,3)
ここで、式09xにおけるKi_a_x、Ki_b_xは、それぞれ、第iゲイン係数Ki_xの最小側(“0”に近い側)のゲイン係数値、最大側(“0”から離れる側)のゲイン係数値としてあらかじめ設定された定数値である。このことは、式09yにおけるKi_a_y、Ki_b_yについても同様である。
従って、式07xの演算に用いる各第iゲイン係数Ki_x(i=1,2,3)は、それぞれに対応する定数値Ki_a_x、Ki_b_xの重み付き平均値として決定される。そして、この場合、Ki_a_x、Ki_b_xにそれぞれ掛かる重みが、ゲイン調整パラメータKr_xに応じて変化させられる。このため、Kr_x=0である場合には、Ki_x=Ki_a_xとなり、Kr_x=1である場合には、Ki_x=Ki_b_xとなる。そして、Kr_xが“0”から“1”に近づくに伴い、第iゲイン係数Ki_xはKi_a_xからKi_b_x近づいていく。
同様に、式07yの演算に用いる各第iゲイン係数Ki_y(i=1,2,3)は、それぞれに対応する定数値Ki_a_y、Ki_b_yの重み付き平均値として決定される。そして、この場合、Ki_a_y、Ki_b_yにそれぞれ掛かる重みが、ゲイン調整パラメータKr_yに応じて変化させられる。このため、Ki_xの場合と同様に、Kr_yの値が“0”から“1”の間で変化するに伴い、第iゲイン係数Ki_yの値が、Ki_a_yとKi_b_yとの間で変化する。
補足すると、上記定数値Ki_a_x、Ki_b_x及びKi_a_y,Ki_b_y(i=1,2,3)は、前記ステップS6又は8において値が設定される定数パラメータに含まれるものである。
姿勢制御演算部80は、上記の如く決定した第1〜第3ゲイン係数K1_x,K2_x,K3_xを用いて前記式07xの演算を行なうことで、X軸方向に輪転する仮想車輪62_xに係わる仮想車輪回転角加速度指令ωwdot_x_cmdを算出する。
さらに詳細には、図13を参照して、姿勢制御演算部80は、基体傾斜角度偏差計測値θbe_x_sに第1ゲイン係数K1_xを乗じてなる操作量成分u1_xと、基体傾斜角速度計測値θbdot_x_sに第2ゲイン係数K2_xを乗じてなる操作量成分u2_xとをそれぞれ、処理部80a,80bで算出する。さらに、姿勢制御演算部80は、重心速度推定値Vb_x_sと制御用目標重心速度Vb_x_mdfdとの偏差(=Vb_x_s−Vb_x_mdfd)を演算部80dで算出し、この偏差に第3ゲイン係数K3_xを乗じてなる操作量成分u3_xを処理部80cで算出する。そして、姿勢制御演算部80は、これらの操作量成分u1_x,u2_x,u3_xを演算部80eにて加え合わせることにより、仮想車輪回転角加速度指令ωwdot_x_comを算出する。
同様に、姿勢制御演算部80は、上記の如く決定した第1〜第3ゲイン係数K1_y,K2_y,K3_yを用いて前記式07yの演算を行なうことで、Y軸方向に輪転する仮想車輪62_yに係わる仮想車輪回転角加速度指令ωwdot_y_cmdを算出する。
この場合には、姿勢制御演算部80は、基体傾斜角度偏差計測値θbe_y_sに第1ゲイン係数K1_yを乗じてなる操作量成分u1_yと、基体傾斜角速度計測値θbdot_y_sに第2ゲイン係数K2_yを乗じてなる操作量成分u2_yとをそれぞれ、処理部80a,80bで算出する。さらに、姿勢制御演算部80は、重心速度推定値Vb_y_sと制御用目標重心速度Vb_y_mdfdとの偏差(=Vb_y_s−Vb_y_mdfd)を演算部80dで算出し、この偏差に第3ゲイン係数K3_yを乗じてなる操作量成分u3_yを処理部80cで算出する。そして、姿勢制御演算部80は、これらの操作量成分u1_y,u2_y,u3_yを演算部80eにて加え合わせることにより、仮想車輪回転角加速度指令ωwdot_x_comを算出する。
ここで、式07xの右辺の第1項(=第1操作量成分u1_x)及び第2項(=第2操作量成分u2_x)は、Y軸周り方向での基体傾斜角度偏差計測値θbe_x_sを、フィードバック制御則としてのPD則(比例・微分則)により“0”に収束させる(基体傾斜角度計測値θb_x_sを目標値θb_x_objに収束させる)ためのフィードバック操作量成分としての意味を持つ。
また、式07xの右辺の第3項(=第3操作量成分u3_x)は、重心速度推定値Vb_x_sと目標重心速度Vb_x_mdfdとの偏差をフィードバック制御則としての比例則により“0”に収束させる(Vb_x_sをVb_x_mdfdに収束させる)ためのフィードバック操作量成分としての意味を持つ。
これらのことは、式07yの右辺の第1〜第3項(第1〜第3操作量成分u1_y,u2_y,u3_y)についても同様である。
姿勢制御演算部80は、上記の如く、仮想車輪回転角加速度指令ωwdot_x_com,ωwdot_y_comを算出した後、次に、これらのωwdot_x_com,ωwdot_y_comをそれぞれ積分器80fにより積分することによって、前記仮想車輪回転速度指令ωw_x_com,ωw_y_comを決定する。
以上が、連携動作に関連する処理を除いた姿勢制御演算部80の処理の詳細である。
補足すると、式07xの右辺の第3項を、Vb_x_sに応じた操作量成分(=K3_x・Vb_x_s)と、Vb_x_mdfdに応じた操作量成分(=−K3_x・Vb_x_mdfd)とに分離した式によって、仮想車輪回転角加速度指令ωdotw_x_comを算出するようにしてよい。同様に、式07yの右辺の第3項を、Vb_y_sに応じた操作量成分(=K3_y・Vb_y_s)と、Vb_y_mdfdに応じた操作量成分(=−K3_y・Vb_y_mdfd)とに分離した式によって、仮想車輪回転角加速度指令ωdotw_y_comを算出するようにしてよい。
また、本実施形態では、車両系重心点の挙動を制御するための操作量(制御入力)として、仮想車輪62_x,62_yの回転角加速度指令ωw_x_cmd,ωw_y_cmdを用いるようにしたが、例えば、仮想車輪62_x,62_yの駆動トルク、あるいは、この駆動トルクに各仮想車輪62_x,62_yの半径Rw_x,Rw_yを乗じてなる並進力(すなわち仮想車輪62_x,62_yと床面との間の摩擦力)を操作量として用いるようにしてもよい。
図9の説明に戻って、制御ユニット50は、次に、姿勢制御演算部80で上記の如く決定した仮想車輪回転速度指令ωw_x_com,ωw_y_comをモータ指令演算部82に入力し、該モータ指令演算部82の処理を実行することによって、電動モータ31Rの速度指令ω_R_comと電動モータ31Lの速度指令ω_L_comとを決定する。このモータ指令演算部82の処理は、前記リミット処理部86(図11参照)のXY−RL変換部86bの処理と同じである。
具体的には、モータ指令演算部82は、前記式01a,01bのωw_x,ωw_y,ω_R,ω_Lをそれぞれ、ωw_x_com,ωw_y_com,ω_R_cmd,ω_L_cmdに置き換えて得られる連立方程式を、ω_R_cmd,ω_L_cmdを未知数として解くことによって、電動モータ31R,31Lのそれぞれの速度指令ω_R_com,ω_L_comを決定する。
以上により前記ステップS9の車両制御演算処理が完了する。
以上説明した如く制御ユニット50が制御演算処理を実行することによって、前記搭乗モード及び自立モードのいずれの動作モードにおいても、基本的には、基体9の姿勢が、前記基体傾斜角度偏差計測値θbe_x_s,θbe_y_sの両方が“0”となる姿勢(以下、この姿勢を基本姿勢という)に保たれるように、換言すれば、車両系重心点(車両・乗員全体重心点又は車両単体重心点)の位置が、車輪体5の接地面のほぼ真上に位置する状態に保たれるように、操作量(制御入力)としての仮想車輪回転角加速度指令ωdotw_xy_comが決定される。より詳しく言えば、基体9の姿勢を前記基本姿勢に保ちつつ、車両系重心点の移動速度の推定値としての重心速度推定値Vb_xy_sを制御用目標重心速度Vb_xy_mdfdに収束させるように、仮想車輪回転角加速度指令ωdotw_xy_comが決定される。なお、制御用目標重心速度Vb_xy_mdfdは、通常は(詳しくは搭乗モードで乗員等が車両1の付加的な推進力を付与しない限り)、“0”である。この場合には、基体9の姿勢を前記基本姿勢に保ちつつ、車両系重心点がほぼ静止するように、仮想車輪回転角加速度指令ωdotw_xy_comが決定されることとなる。
そして、ωdotw_xy_comの各成分を積分してなる仮想車輪回転角速度指令ωw_xy_comを変換してなる電動モータ31R,31Lのそれぞれの回転角速度が、電動モータ31R,31Lの速度指令ω_R_cmd,ω_L_cmdとして決定される。さらに、その速度指令ω_R_cmd,ω_L_cmdに従って、各電動モータ31R,31Lの回転速度が制御される。ひいては車輪体5のX軸方向及びY軸方向のそれぞれの移動速度が、ωw_x_comに対応する仮想車輪62_xの移動速度と、ωw_y_comに対応する仮想車輪62_yの移動速度とに各々一致するように制御される。
このため、例えば、Y軸周り方向で、実際の基体傾斜角度θb_xが目標値θb_x_objから前傾側にずれると、そのずれを解消すべく(θbe_x_sを“0”に収束させるべく)、車輪体5が前方に向かって移動する。同様に、実際のθb_xが目標値θb_x_objから後傾側にずれると、そのずれを解消すべく(θbe_x_sを“0”に収束させるべく)、車輪体5が後方に向かって移動する。
また、例えば、X軸周り方向で、実際の基体傾斜角度θb_yが目標値θb_y_objから右傾側にずれると、そのずれを解消すべく(θbe_y_sを“0”に収束させるべく)、車輪体5が右向きに移動する。同様に、実際のθb_yが目標値θb_y_objから左傾側にずれると、そのずれを解消すべく(θbe_y_sを“0”に収束させるべく)、車輪体5が左向きに移動する。
さらに、実際の基体傾斜角度θb_x,θb_yの両方が、それぞれ目標値θb_x_obj,θb_y_objからずれると、θb_xのずれを解消するための車輪体5の前後方向の移動動作と、θb_yのずれを解消するための車輪体5の左右方向の移動動作とが合成され、車輪体5がX軸方向及びY軸方向の合成方向(X軸方向及びY軸方向の両方向に対して傾斜した方向)に移動することとなる。
このようにして、基体9が前記基本姿勢から傾くと、その傾いた側に向かって、車輪体5が移動することとなる。従って、例えば前記搭乗モードにおいて、乗員が意図的にその上体を傾けると、その傾けた側に、車輪体5が移動することとなる。
なお、制御用目標重心速度Vb_x_mdfd,Vb_y_mdfdが“0”である場合には、基体9の姿勢が基本姿勢に収束すると、車輪体5の移動もほぼ停止する。また、例えば、基体9のY軸周り方向の傾斜角度θb_xを基本姿勢から傾いた一定の角度に維持すると、車輪体5のX軸方向の移動速度は、その角度に対応する一定の移動速度(制御用目標重心速度Vb_x_mdfdと一定の定常偏差を有する移動速度)に収束する。このことは、基体9のX軸周り方向の傾斜角度θb_yを基本姿勢から傾いた一定の角度に維持した場合も同様である。
また、例えば、前記要求重心速度生成部74で生成される要求重心速度Vb_x_aim,Vb_y_aimの両方が“0”となっている状況において、基体9の前記基本姿勢からの傾き量(基体傾斜角度偏差計測値θbe_x_s,θbe_y_s)が比較的大きくなり、それを解消し、もしくはその傾き量を維持するために必要な車輪体5のX軸方向及びY軸方向の一方又は両方の移動速度(これらの移動速度は、それぞれ、図12に示した前記重心速度定常偏差予測値Vb_x_prd、Vb_y_prdに相当する)が、電動モータ31R,31Lの一方又は両方の回転角速度を許容範囲から逸脱させてしまうような、過大な移動速度になるような状況では、該車輪体5の移動速度に対して逆向きとなる速度(詳しくは、Vw_x_lim2−Vb_x_prd及びVw_y_lim2−Vb_y_prd)が制御用目標重心速度Vb_x_mdfd,Vb_y_mdfdとして決定される。そして、制御入力を構成する操作量成分のうちの操作量成分u3_x,u3_yが、この制御用目標重心速度Vb_x_mdfd,Vb_y_mdfdに重心速度推定値Vb_x_s,Vb_y_sをそれぞれ収束させるように決定される。このため、基体9の前記基本姿勢からの傾き量が過大になるのを予防し、ひいては、電動モータ31R,31Lの一方又は両方の回転角速度が高速になり過ぎるのが防止される。
さらに、前記ゲイン調整部78では、重心速度推定値Vb_x_s,Vb_y_sの一方又は両方が大きくなり、ひいては、基体9の前記基本姿勢からの傾きを解消し、もしくはその傾き量を維持するために必要な車輪体5のX軸方向及びY軸方向の一方又は両方の移動速度が、電動モータ31R,31Lの一方又は両方の回転角速度を許容範囲から逸脱させてしまうような、過大な移動速度になる恐れがある状況では、その逸脱が顕著になるほど(詳しくは、図10に示すVover_x,Vover_yの絶対値が大きくなるほど)、前記ゲイン調整パラメータKr_x,Kr_yの一方又は両方が“0”から“1”に近づけられる。
この場合、前記式09xにより算出される各第iゲイン係数Ki_x(i=1,2,3)は、Kr_xが“1”に近づくほど、最小側の定数値Ki_a_xから最大側の定数値Ki_b_xに近づく。このことは、前記式09yにより算出される各第iゲイン係数Ki_y(i=1,2,3)についても同様である。
そして、上記ゲイン係数の絶対値が大きくなることによって、基体9の傾きの変化に対する操作量(仮想車輪回転角加速度指令ωdotw_x_cmd,ωdotw_y_cmd)の感度が高まる。従って、基体9の基本姿勢からの傾き量が大きくなろうとすると、それを素早く解消するように、車輪体5の移動速度が制御されることとなる。従って、基体9が基本姿勢から大きく傾くことが強めに抑制され、ひいては、車輪体5のX軸方向及びY軸方向の一方又は両方の移動速度が、電動モータ31R,31Lの一方又は両方の回転角速度を許容範囲から逸脱させてしまうような、過大な移動速度になるのを防止することができる。
また、搭乗モードにおいて、要求重心速度生成部74が、乗員等の操縦操作による要求に応じて要求重心速度Vb_x_aim,Vb_y_aim(Vb_x_aim,Vb_y_aimの一方又は両方が“0”でない要求重心速度)を生成した場合には、電動モータ31R,31Lの一方又は両方の回転角速度が許容範囲を逸脱するような高速の回転角速度にならない限り(詳しくは図12に示すVw_x_lim2,Vw_y_lim2がVb_x_t,Vb_y_tにそれぞれ一致する限り)、要求重心速度Vb_x_aim,Vb_y_aimがそれぞれ前記制御用目標重心速度Vb_x_mdfd,Vb_y_mdfdとして決定される。このため、要求重心速度Vb_x_aim,Vb_y_aimを実現するように(実際の重心速度が要求重心速度Vb_x_aim,Vb_y_aimに近づくように)、車輪体5の移動速度が制御される。
次に、説明を後回しにした、前記要求重心速度生成部74の処理の詳細を説明する。
要求重心速度生成部74は、本実施形態では、車両1の動作モードが自立モードである場合には、前記したように、要求重心速度Vb_x_aim,Vb_y_aimを“0”とする。
一方、要求重心速度生成部74は、車両1の動作モードが搭乗モードである場合には、乗員等による車両1の操縦操作(車両1に推進力を付加する操作)に応じて、該操作により要求されると推定される要求重心速度Vb_x_aim,Vb_y_aimを決定する。
ここで、例えば、車両1の乗員が、車両1の発進時等において、車両1の移動速度(車両系重心点の移動速度)を積極的に増速させようとする場合には、自身の足により床を蹴り、それにより車両1に移動速度を増速させる推進力(乗員の足平と床との摩擦力による推進力)を車両1に付加する。あるいは、例えば、車両1の乗員の要求に応じて、外部の補助者等が、車両1にその移動速度を増速させる推進力を付加する場合もある。
このような場合に、要求重心速度生成部74は、車両系重心点の実際の速度ベクトル(以降、重心速度ベクトル↑Vbという)の大きさ(絶対値)の時間的変化率に基づいて、車両1の移動速度を増加させる要求としての加速要求が発生したか否かを判断しつつ、それに応じて、↑Vbの目標値としての要求重心速度ベクトル↑Vb_aim(要求重心速度Vb_x_aim,Vb_y_aimを2つの成分とする速度ベクトル)を逐次決定する。
その処理を概略的に説明すると、上記加速要求が発生した場合には、該加速要求が解消するまでの間、要求重心速度ベクトル↑Vb_aimの大きさを増加させるように、要求重心速度ベクトル↑Vb_aimが決定される。そして、上記加速要求が解消すると、要求重心速度ベクトル↑Vb_aimの大きさを段階的に減衰させていくように、要求重心速度ベクトル↑Vb_aimが決定される。この場合、本実施形態では、基本的には、加速要求が解消してから、所定時間の期間、要求重心速度ベクトル↑Vb_aimの大きさが一定に保たれる。そして、その後に、要求重心速度ベクトル↑Vb_aimの大きさが“0”まで連続的に減衰させられる。なお、この減衰時には、適宜、要求重心速度ベクトル↑Vb_aimの向きがX軸方向に近づけられる。
このような処理を実行する要求重心速度生成部74を、図14〜図20のフローチャートを参照して以下に詳細に説明する。
図14を参照して、要求重心速度生成部74は、まず、ステップS21の処理を実行する。この処理では、要求重心速度生成部74は、入力される重心速度推定値Vb_x_s,Vb_y_sを2つの成分とする速度ベクトル(実際の重心速度ベクトル↑Vbの観測値)である推定重心速度ベクトル↑Vb_sの大きさ|↑Vb_s|(=sqrt(Vb_x_s2+Vb_y_s2))の時間的変化率(微分値)DVb_sを算出する。このDVb_sは、実際の重心速度ベクトル↑Vbの大きさの時間的変化率の観測値(推定値)としての意味を持つ。以降、DVb_sを推定重心速度絶対値変化率DVb_sという。なお、上記sqrt( )は、平方根関数である。
さらに、ステップS21では、要求重心速度生成部74は、入力される重心速度推定値Vb_x_s,Vb_y_sのそれぞれの時間的変化率(微分値)である重心加速度推定値Vbdot_x_s,Vvdot_y_sを算出する。なお、Vbdot_x_s,Vvdot_y_sを2つの成分とするベクトルは、車両系重心点の実際の加速度ベクトルの観測値を意味する。
次いで、ステップS22に進んで、要求重心速度生成部74は、要求重心速度Vb_x_aimを算出するための現在の演算処理モードが、どのモードであるかを判断する。
ここで、本実施形態では、要求重心速度生成部74は、要求重心速度ベクトル↑Vb_aimの基本値(以降、基本要求重心速度ベクトル↑Vb_aim1という)を決定した上で、この基本要求重心速度ベクトル↑Vb_aim1に要求重心速度ベクトル↑Vb_aimを追従させるように(定常的には一致させるように)、要求重心速度ベクトル↑Vb_aimを決定する。
上記演算処理モードは、上記基本要求重心速度ベクトル↑Vb_aim1の決定の仕方の種別を表すものである。そして、本実施形態では、該演算処理モードとしては、制動モード、速度追従モード、速度ホールドモードの3種類のモードがある。
制動モードは、基本要求重心速度ベクトル↑Vb_aim1の大きさを“0”に減衰させていくか、又は“0”に保持するように↑Vb_aim1を決定するモードである。また、速度追従モードは、基本要求重心速度ベクトル↑Vb_aim1を、推定重心速度ベクトル↑Vb_sに追従させるように(一致又はほぼ一致させるように)決定するモードである。また、速度ホールドモードは、基本要求重心速度ベクトル↑Vb_aim1の大きさを一定に保つように↑Vb_aim1を決定するモードである。
なお、制御ユニット50の起動時等に該制御ユニット50が初期化された状態での演算処理モード(初期演算処理モード)は、制動モードである。
要求重心速度生成部74は、上記ステップS22において、現在の演算処理モードが制動モードである場合と、速度追従モードである場合と、速度ホールドモードである場合とで、それぞれ、次に、ステップS23の演算処理、ステップS24の演算処理、ステップS25の演算処理を実行し、基本要求重心速度ベクトル↑Vb_aim1を決定する。
これらの各モードに対応する演算処理は、次のように実行される。
ステップS23における制動モードの演算処理は、図15のフローチャートに示す如く実行される。具体的には、要求重心速度生成部74は、まず、前記ステップS21で算出した推定重心速度絶対値変化率DVb_sと、重心加速度推定値Vbdot_x_s,Vbdot_y_sとに関して、DVb_s>DV1、且つ、|Vbdot_x_s|>a1*|Vbdot_y_s|という条件が成立するか否かをステップS23−1で判断する。この判断処理は、車両1の移動速度を車両1の概略前後方向で増速させようとする加速要求が有るか否かを判断する処理である。
ここで、上記DV1は、あらかじめ設定された正の値の第1閾値DV1(>0)である。そして、DVb_s>DV1であるということは、実際の重心速度ベクトル↑Vbの大きさ|↑Vb|が第1閾値DV1よりも大きな時間的変化率で増加している状況を意味する。
また、上記a1は、あらかじめ設定された正の値の係数値である。そして、|Vbdot_x_s|>a1*|Vbdot_y_s|であるということは、車両系重心点の実際の加速度ベクトルが“0”でないX軸方向の成分を有し、且つ、該加速度ベクトルのX軸方向に対する鋭角側の角度(=tan-1(|Vbdot_y_s|/|Vbdot_x_s|)が、所定角度(=tan-1(1/a1))よりも“0”に近い状況であることを意味する。本実施形態では、a1は、例えば“1”又はこれに近い値に設定されている。
従って、ステップS23−1の判断結果が肯定的になる状況は、乗員あるいは外部の補助者等により概略前後方向での重心速度ベクトル↑Vbの大きさを増加させようとする操縦操作(車両1に概略前後方向の推進力を付加する操縦操作)が行なわれている状況である。
ステップS23−1の判断結果が否定的となる場合、すなわち、車両1の加速要求(概略前後方向での車両1の加速要求)が無い場合には、要求重心速度生成部74は、次に、ステップS23−4の判断処理を実行する。
ステップS23−4の判断処理では、要求重心速度生成部74は、ステップS21で算出した推定重心速度絶対値変化率DVb_sが、あらかじめ設定された負の値の第3閾値DV3(<0)よりも小さいか否かを判断する。この判断処理は、車両1の乗員が重心速度ベクトル↑Vbの大きさを積極的に減少させようとする減速要求が発生したか否かを判断するものである。この場合、車両1の乗員が意図的に自身の足を接地させ、自身の足と床との間に車両1の制動方向の摩擦力を発生させたような場合にステップS23−4の判断結果が肯定的になる。
そして、要求重心速度生成部74は、ステップS23−4の判断結果が否定的である場合(減速要求が発生していない場合)には、ステップS23−5で第1制動演算処理を実行することによって、基本要求重心速度ベクトル↑Vb_aim1の大きさ|↑Vb_aim1|(以降、基本要求重心速度ベクトル絶対値|↑Vb_aim1|という)と方位角θvb_aim1(以降、基本要求重心速度ベクトル方位角θvb_aim1という)とを決定し、図15の処理を終了する。また、ステップS23−4の判断結果が肯定的である場合(減速要求が発生した場合には)、要求重心速度生成部74は、ステップS23−6で第2制動演算処理を実行することによって、基本要求重心速度ベクトル絶対値|↑Vb_aim1|と基本要求重心速度ベクトル方位角θvb_aim1とを決定し、図15の処理を終了する。
なお、本実施形態では、基本要求重心速度ベクトル方位角θvb_aim1は、sin(θvb_aim1)=Vb_x_aim1/|↑Vb_aim1|,cos(θvb_aim1)=Vb_y_aim1/|↑Vb_aim1|となる角度(180°≦θVb_aim≦180°)として定義される。また、|↑Vb_aim|=0である場合には、θVb_aim=0°であるとする。
上記ステップS23−5の第1制動演算処理は、図16及び図17のフローチャートに示す如く実行される。
この第1制動演算処理では、要求重心速度生成部74は、まず、ステップS23−5−1において基本要求重心速度ベクトル絶対値|↑Vb_aim1|の前回値|↑Vb_aim_p|からあらかじめ設定された正の所定値ΔVb1だけ減少させてなる値を、|↑Vb_aim1|の候補値ABS_Vbとして算出する。ΔVb1は、制御処理周期毎の|↑Vb_aim1|の減少量(ひいては|↑Vb_aim1|の時間的変化率)を規定する設定値である。
次いで、要求重心速度生成部74は、ステップS23−5−2において、この候補値ABS_Vbと、“0”とのうちの大きい方の値max(0,ABS_Vb)を、|↑Vb_aim1|の今回値として決定する。従って、ABS_Vb≧0である場合には、ABS_Vbがそのまま|↑Vb_aim1|の今回値として決定され、ABS_Vb<0である場合には、|↑Vb_aim1|の今回値は“0”とされる。
次いで、要求重心速度生成部74は、上記の如く決定した|↑Vb_aim1|が“0”であるか否かをステップS23−5−3で判断する。この判断結果が肯定的である場合には、要求重心速度生成部74は、次に、ステップS23−5−4において、θvb_aim1の今回値を0°とし、図16の処理を終了する。
ステップS23−5−3の判断結果が否定的である場合には、要求重心速度生成部74は、次に、ステップS23−5−5からの処理によって、θvb_aim1の前回値θvb_aim1_pの値が、0°≦θvb_aim1_p≦θth1+、θth1-≦θvb_aim1_p<0°、θth2+≦θvb_aim1_p≦180°、−180°≦θvb_aim1_p≦θth2-、θth1+<θvb_aim1_p<θth2+、θth2-<θvb_aim1_p<θth1-のいずれの範囲内の値であるかに応じて、θvb_aim1の今回値を決定する。
ここで、θth1+は、0°と90°との間の値であらかじめ設定された正の方位角閾値、θth1-は、0°と−90°との間の値であらかじめ設定された負の方位角閾値、θth2+は、90°と180°との間の値であらかじめ設定された正の方位角閾値、θth2-は、−90°と−180°との間の値であらかじめ設定された負の方位角閾値である。本実施形態では、θth1+とθth1-とは、それらの絶対値が互いに同一の値(例えば45°もしくはその近辺の角度値)に設定されている。また、θth2+とθth2-とは、それらの絶対値が互いに同一の値(例えば135°もしくはその近辺の角度値)に設定されている。なお、θth1+とθth1-との差(=(θth1+)−(θth1-))と、θth2とθth2-との差(=(θth2+)−(θth2-))とは、互いに同一である必要なない。
ステップS23−5−5からの処理は以下に説明する如く実行される。すなわち、要求重心速度生成部74は、ステップS23−5−5において、0°≦θvb_aim1_p≦θth1+であるか否かを判断する。この判断結果が肯定的である場合には、要求重心速度生成部74は、ステップS23−5−6において、θvb_aim1の前回値θvb_aim1_pからあらかじめ設定された正の所定値Δθvb1だけ減少させてなる値を、θvb_aim1の候補値ANG_Vbとして算出する。Δθvb1は、制御処理周期毎のθvb_aim1の変化量(ひいてはθvb_aim1の時間的変化率)を規定する設定値である。
そして、要求重心速度生成部74は、ステップS23−5−7において、この候補値ANG_Vbと、0°とのうちの大きい方の角度値max(0,ANG_Vb)をθvb_aim1の今回値として決定し、図16の処理を終了する。従って、ANG_Vb≧0°である場合には、ANG_Vbがそのままθvb_aim1の今回値として決定され、ANG_Vb<0°である場合には、θvb_aim1の今回値は0°とされる。
ステップS23−5−5の判断結果が否定的である場合には、要求重心速度生成部74は、次に、ステップS23−5−8において、θth1-≦θvb_aim1_p<0°であるか否かを判断する。この判断結果が肯定的である場合には、要求重心速度生成部74は、ステップS23−5−9において、θvb_aim1の前回値θvb_aim1_pから前記所定値Δθvb1だけ増加させてなる値を、θvb_aim1の候補値ANG_Vbとして算出する。
そして、要求重心速度生成部74は、ステップS23−5−10において、この候補値ANG_Vbと、0°とのうちの小さい方の角度値min(0,ANG_Vb)をθvb_aim1の今回値として決定し、図16の処理を終了する。従って、ANG_Vb≦0°である場合には、ANG_Vbがそのままθvb_aim1の今回値として決定され、ANG_Vb>0°である場合には、θvb_aim1の今回値は0°とされる。
ステップS23−5−8の判断結果が否定的である場合には、要求重心速度生成部74は、次に、図17のステップS23−5−11において、θth2+≦θvb_aim1_p≦180°であるか否かを判断する。この判断結果が肯定的である場合には、要求重心速度生成部74は、ステップS23−5−12において、θvb_aim1の前回値θvb_aim1_pから前記所定値Δθvb1だけ増加させてなる値を、θvb_aim1の候補値ANG_Vbとして算出する。
そして、要求重心速度生成部74は、ステップS23−5−13において、この候補値ANG_Vbと、0°とのうちの小さい方の角度値min(180,ANG_Vb)をθvb_aim1の今回値として決定し、図17の処理を終了する。従って、ANG_Vb≦180°である場合には、ANG_Vbがそのままθvb_aim1の今回値として決定され、ANG_Vb>180°である場合には、θvb_aim1の今回値は180°とされる。
ステップS23−5−11の判断結果が否定的である場合には、要求重心速度生成部74は、次に、ステップS23−5−14において、−180°≦θvb_aim1_p≦θth2-であるか否かを判断する。この判断結果が肯定的である場合には、要求重心速度生成部74は、ステップS23−5−15において、θvb_aim1の前回値θvb_aim1_pから前記所定値Δθvb1だけ減少させてなる値を、θvb_aim1の候補値ANG_Vbとして算出する。
そして、要求重心速度生成部74は、ステップS23−5−16において、この候補値ANG_Vbと、−180°とのうちの大きい方の角度値max(180,ANG_Vb)をθvb_aim1の今回値として決定し、図17の処理を終了する。従って、ANG_Vb≧−180°である場合には、ANG_Vbがそのままθvb_aim1の今回値として決定され、ANG_Vb<−180°である場合には、θvb_aim1の今回値は−180°とされる。
ステップS23−5−14の判断結果が否定的である場合、すなわち、θth1+<θvb_aim1_p<θth2+又はθth2-<θvb_aim1_p<θth1-である場合には、要求重心速度生成部74は、ステップS23−5−17において、θvb_aim1の今回値を、前回値θvb_aim1_pと同じ値に決定し、図17の処理を終了する。
以上が、ステップS23−5の第1制動演算処理の詳細である。
一方、前記ステップS23−6の第2制動演算処理は、図18のフローチャートに示す如く実行される。
この第2制動演算処理では、要求重心速度生成部74は、まず、ステップS23−6−1において基本要求重心速度ベクトル絶対値|↑Vb_aim1|の前回値|↑Vb_aim_p|からあらかじめ設定された正の所定値ΔVb2だけ減少させてなる値を、|↑Vb_aim1|の候補値ABS_Vbとして算出する。ΔVb2は、第2制動演算処理での、制御処理周期毎の|↑Vb_aim1|の減少量(ひいては|↑Vb_aim1|の時間的変化率)を規定する設定値である。この場合、ΔVb2は、前記第1制動演算処理で使用する所定値ΔVb1よりも大きな値に設定されている。
次いで、要求重心速度生成部74は、ステップS23−6−2において、前記ステップS23−5−2と同様の処理を実行し、ステップS23−6−1で算出した候補値ABS_Vbと、“0”とのうちの大きい方の値max(0,ABS_Vb)を、|↑Vb_aim1|の今回値として決定する。
次いで、要求重心速度生成部74は、上記の如く決定した|↑Vb_aim1|が“0”であるか否かをステップS23−6−3で判断する。この判断結果が肯定的である場合には、要求重心速度生成部74は、次に、ステップS23−6−4において、θvb_aim1の今回値を“0”とし、図18の処理を終了する。
また、ステップS23−6−3の判断結果が否定的である場合には、要求重心速度生成部74は、次に、ステップS23−6−5において、θvb_aim1の今回値を、前回値θvb_aim1_pと同じ値に決定し、図18の処理を終了する。
以上が、ステップS23−6の第2制動演算処理の詳細である。
図15の説明に戻って、前記ステップS23−1の判断結果が肯定的である場合、すなわち、概略前後方向での車両1の加速要求が有る場合には、要求重心速度生成部74は、ステップS23−2で基本要求重心速度ベクトル絶対値|↑Vb_aim1|と基本要求重心速度ベクトル方位角θvb_aim1とを決定する。要求重心速度生成部74は、ステップS23−3で演算処理モードを制動モードから速度追従モードに変更し、図15の処理を終了する。
上記ステップS23−2では、具体的には、前記推定重心速度ベクトル↑Vb_s(今回値)の大きさ|↑Vb_s|(=sqrt(Vb_x_s2+Vb_y_s2))に、あらかじめ設定された所定値の比率γを乗じてなる値が、基本要求重心速度ベクトル絶対値|↑Vb_aim1|として決定される。上記比率γは、本実施形態では、“1”よりも若干小さい正の値(例えば0.8)に設定されている。
さらに、ステップS23−2では、推定重心速度ベクトル↑Vb_sの方位角θvb_s(=sin-1(Vb_x_s/|↑Vb_s|))がそのまま、基本要求重心速度ベクトル方位角θvb_aim1として決定される。従って、ステップS23−2では、結果的には、推定重心速度ベクトル↑Vb_sに上記比率γを乗じてなるベクトルが、基本要求重心速度ベクトル↑Vb_aim1として決定されることとなる。
このようなステップS23−2の処理は、|↑Vb_x_aim1|とθvb_aim1との決定の仕方を、次回の制御処理サイクルから開始する速度追従モードに合わせるものである。
なお、上記比率γの値が“1”よりも若干小さいことは必須ではなく、例えば、該比率γの値を“1”あるいはそれよりも若干大きい値に設定してもよい。本実施形態では、乗員が体感的(感覚的)に認識する車両1の移動速度が、実際の移動速度に比べて高速であるかのように認識されるのを防止するために、比率γの値を“1”よりも若干小さい値に設定している。
以上が、ステップS23における制動モードの演算処理である。
なお、ステップS23−1の判断結果が否定的となる場合には、演算処理モードは変更されないので、次回の制御処理周期においても、演算処理モードは、制動モードに維持されることとなる。
次に、ステップS24における速度追従モードの演算処理は、図19のフローチャートに示す如く実行される。具体的には、要求重心速度生成部74は、まず、ステップS24−1において、前記ステップS23−4と同じ判断処理、すなわち、車両1の減速要求が発生したか否かの判断処理を実行する。
この判断結果が肯定的である場合には、要求重心速度生成部74は、次に、ステップS24−6において、前記ステップS23−6と同じ処理(図18のフローチャートに示す処理)を実行することによって、基本要求重心速度ベクトル絶対値|↑Vb_aim1|と基本要求重心速度ベクトル方位角θvb_aim1とを決定する。さらに、要求重心速度生成部74は、ステップS24−7において、演算処理モードを速度追従モードから制動モードに変更し、図19の処理を終了する。
一方、前記ステップS24−1の判断結果が否定的である場合、すなわち、車両1の減速要求が発生していない場合には、要求重心速度生成部74は、次にステップS24−2の処理を実行する。このステップS24−2では、要求重心速度生成部74は、前記ステップS23−2と同じ処理を実行し、基本要求重心速度ベクトル絶対値|↑Vb_aim1|と基本要求重心速度ベクトル方位角θvb_aim1とを決定する。すなわち、|↑Vb_x_s|*γが|↑Vb_aim1|として決定される共に、θvb_sがθvb_aim1として決定される。
次いで、要求重心速度生成部74は、ステップS24−3において、前記推定重心速度絶対値変化率DVb_s(ステップS21で算出した値)が、あらかじめ設定された第2閾値DV2よりも小さいか否かを判断する。この第2閾値DV2は、本実施形態では、前記第3閾値DV3よりも大きい(DV3よりも“0”に近い)負の所定値に設定されている。なお、第2閾値DV2は、“0”もしくは“0”よりも若干大きい正の値(ただし、前記第1閾値DV1よりも小さい値)に設定されていてもよい。
このステップS24−3の判断処理は、速度追従モードから、前記速度ホールドモードへの移行タイミングを判断するものである。そして、要求重心速度生成部74は、ステップS24−3の判断結果が否定的である場合には、そのまま図19の処理を終了する。この場合には、演算処理モードは変更されないので、次回の制御処理周期においても、演算処理モードは、速度追従モードに維持されることとなる。
また、ステップS24−3の判断結果が肯定的である場合には、要求重心速度生成部74は、車両1の加速要求が完了したものとみなして、ステップS24−4において、カウントダウンタイマを初期化する。そして、要求重心速度生成部74は、ステップS24−5にて演算処理モードを速度追従モードから速度ホールドモードに変更し、図19の処理を終了する。
上記カウントダウンタイマは、次回の制御処理周期から開始する速度ホールドモードの開始後の経過時間を計時するタイマである。そして、ステップS24−4では、該タイマの計時値CNTに、あらかじめ設定された初期値Tmがセットされる。初期値Tm_xは、速度ホールドモードを継続させようとする時間の設定値を意味する。
以上が、ステップS24における速度追従モードの演算処理である。
次に、ステップS25における速度ホールドモードの演算処理は、図20のフローチャートに示す如く実行される。具体的には、要求重心速度生成部74は、まず、ステップS25−1において、前記ステップS23−4と同じ判断処理、すなわち、車両1の減速要求が発生したか否かの判断処理を実行する。
そして、このステップS25−1の判断結果が肯定的である場合(車両1の減速要求が発生した場合)には、要求重心速度生成部74は、次に、ステップS25−2において、前記ステップS23−6と同じ処理(図18のフローチャートに示す処理)を実行することによって、基本要求重心速度ベクトル絶対値|↑Vb_aim1|と基本要求重心速度ベクトル方位角θvb_aim1とを決定する。さらに、要求重心速度生成部74は、ステップS25−3において、演算処理モードを速度ホールドモードから制動モードに変更し、図20の処理を終了する。
一方、前記ステップS25−1の判断結果が否定的である場合(車両1の減速要求が発生していない場合)には、要求重心速度生成部74は、前記ステップS23−1と同じ判断処理、すなわち概略前後方向での車両1の加速要求が有るか否かの判断処理をステップS25−4で実行する。
そして、ステップS25−4の判断結果が肯定的である場合(概略前後方向での車両1の加速要求が再び発生した場合)には、要求重心速度生成部74は、ステップS25−5において、前記ステップS23−2と同じ処理を実行することによって、基本要求重心速度ベクトル絶対値|↑Vb_aim1|と基本要求重心速度ベクトル方位角θvb_aim1とを決定する。すなわち、|↑Vb_x_s|*γが|↑Vb_aim1|として決定される共に、θvb_sがθvb_aim1として決定される。
さらに、要求重心速度生成部74は、ステップS25−6において、演算処理モードを速度ホールドモードから速度追従モードに変更し、図20の処理を終了する。
前記ステップS25−4の判断結果が否定的である場合(概略前後方向での加速要求が無い状態が継続している場合)には、要求重心速度生成部74は、ステップS25−7において、前記カウントダウンタイマの計時値CNTをデクリメントする。すなわち、計時値CNTの現在値から所定値ΔT(制御処理周期の時間)を差し引くことによって、計時値CNTを更新する。
次いで、要求重心速度生成部74は、カウントダウンタイマの計時値CNTが“0”よりも大きいか否か、すなわち、カウントダウンタイマの計時が終了したか否かをステップS25−8にて判断する。
このステップS25−8の判断結果が肯定的である場合は、速度ホールドモードが開始してから、カウントダウンタイマの前記初期値Tmにより表される時間が未だ経過していない場合である。この場合には、要求重心速度生成部74は、演算処理モードを速度ホールドモードに維持するものとして、基本要求重心速度ベクトル絶対値|↑Vb_aim1|と基本要求重心速度方位角θvb_aim1とをステップS25−9にて決定し、図20の処理を終了する。
この場合、ステップS25−9では、|↑Vb_aim1|の今回値は、前回値|↑Vb_aim1_p|と同じ値に決定される。また、θvb_aim1の今回値は、前回値θvb_aim1_pと同じ値に決定される。従って、基本要求重心速度ベクトル↑Vb_aim1_pの前回値がそのまま、↑Vb_aim1の今回値の速度ベクトルとして決定される。
なお、ステップS25−8の判断結果が肯定的である場合には、演算処理モードは変更されないので、次回の制御処理周期においても、演算処理モードは、速度ホールドモードに維持されることとなる。
前記ステップS25−8の判断結果が否定的となった場合、すなわち、速度ホールドモードが開始してから、カウントダウンタイマの前記初期値Tmにより表される所定時間が経過した場合には、要求重心速度生成部74は、ステップS25−10において、前記ステップS23−5と同じ処理(図16及び図17のフローチャートの処理)を実行することによって、基本要求重心速度ベクトル絶対値|↑Vb_aim1|と基本要求重心速度方位角θvb_aim1とを決定する。
さらに、要求重心速度生成部74は、ステップS25−11にて、演算処理モードを速度ホールドモードから制動モードに変更し、図20の処理を終了する。
以上が、ステップS25における速度ホールドモードの演算処理である。
図14の説明に戻って、要求重心速度生成部74は、以上の如くステップS23〜25のいずれかの演算処理を実行した後、次に、その演算処理により決定した|↑Vb_aim1|とθvb_aim1とをそれぞれフィルタに入力する処理(フィルタリング処理)をステップS26にて実行する。
ここで、|↑Vb_aim1|とθvb_aim1とをそれぞれ入力するフィルタは、特に、演算処理モードが制動モードから速度追従モードに変更された直後に、要求重心速度ベクトル↑Vb_aimの大きさ|↑Vb_aim|と方位角θvb_aimとがステップ状に急変するのを防止するための一次遅れ特性のローパスフィルタである。この場合、|↑Vb_aim|を入力するフィルタの時定数は比較的短い時定数に設定されており、|↑Vb_aim1|が急変した直後以外の状況では該フィルタの出力値が|↑Vb_aim1|に一致又はほぼ一致するようになっている。θvb_aim1を入力するフィルタについても同様である。
そして、ステップS26では、θvb_aim1を入力したフィルタの出力値が、そのまま、要求重心速度ベクトル↑Vb_aimの方位角θvb_aim(以降、要求重心速度ベクトル方位角θvb_aimという)として決定される。
次いで、ステップS27に進んで、要求重心速度生成部74は、|↑Vb_aim1|を入力したフィルタの出力値を、リミッタに通してなる値を、最終的に要求重心速度ベクトル↑Vb_aimの大きさ|↑Vb_aim|(以降、要求重心速度ベクトル絶対値|↑Vb_aim|という)として決定する。この場合、該リミッタは、|↑Vb_aim|が過大になるのを防止するためのものであり、|↑Vb_aim1|を入力したフィルタの出力値があらかじめ設定された所定の上限値以下である場合には、該フィルタの出力値をそのまま|↑Vb_aim|として出力する。また、該リミッタは、フィルタの出力値が上記上限値を超えている場合には、該上限値を|↑Vb_aim|として出力する。換言すれば、該リミッタは、フィルタの出力値と上記上限値とのうちの小さい方の値を|↑Vb_aim|として出力する。
次いで、ステップS28に進んで、要求重心速度生成部74は、上記の如く決定した|↑Vb_aim|とθvb_aimとから、要求重心速度ベクトル↑Vb_aimのX軸方向成分(X軸方向の要求重心速度)Vb_x_aimと、Y軸方向の成分(Y軸方向の要求重心速度)Vb_y_aimとを算出する。より詳しくは、|↑Vb_aim|*sin(θvb_aim)がVb_x_aimとして算出され、|↑Vb_aim|*cos(θvb_aim)がVb_y_aimとして算出される。
以上が要求重心速度生成部74の処理の詳細である。
以上説明した要求重心速度生成部74の処理によって、要求重心速度ベクトル↑Vb_aimは(ひいては要求重心速度Vb_x_aim,Vb_y_aim)は、以下のような態様で決定されることとなる。
すなわち、例えば、車両1の移動速度を増速するために、乗員が自身の足平で床を蹴ることにって、あるいは、補助者等が車両1を押すことによって、車両1に対して、概略X軸方向の推進力(詳しくは、前記ステップS23−1の判断結果が肯定的となるような推進力)を付加した場合を想定する。
なお、推進力を付加する前の演算処理モードは、前記制動モードであるとする。また、ここでは、理解の便宜上、図14のステップS26で|↑Vb_aim1|を入力するフィルタの出力値は、ステップS27でのリミッタによる強制的な制限がかからない程度の範囲内に収まる値(該リミッタの上限値以下の値)であるとする。同様に、重心速度推定値Vb_x_s,Vb_y_sが、前記リミット処理部104での出力値V_x_lim2,V_y_lim2の強制的な制限が行なわれない程度の範囲内に収まるものとする。
この場合、車両1に推進力を付加することによって、前記ステップS23−1の判断結果が肯定的となると、図15のステップS23−3の処理によって、演算処理モードが制動モードから速度追従モードに変更されることとなる。
この速度追従モードでは、減速要求が発生しない状況(ステップS24−1の判断結果が否定的となる状況)で、推定重心速度ベクトル↑Vb_sの今回値(現在値)に、所定値の比率γを乗じてなるベクトル、すなわち、↑Vb_sよりも大きさが若干小さく、且つ、↑Vb_sと同じ向きの速度ベクトルが、基本要求重心速度ベクトル↑Vb_aim1として逐次決定される。
このため、要求重心速度生成部74が逐次決定する要求重心速度ベクトル↑Vb_aimは、車両1に付加された推進力によって増速する(大きさが増加する)実際の重心速度ベクトル↑Vbにほぼ一致する速度ベクトル↑Vb_aim1(=γ*↑Vb_s)に追従するように、決定されることとなる。
そして、このように決定される要求重心速度ベクトル↑Vb_aimのX軸方向成分とY軸方向成分とが前記制御用目標重心速度Vb_x_mdfd,Vb_y_mdfdとして決定される。さらに、各仮想車輪回転角加速度指令ωdotw_x_cmd,ωdotw_y_cmdにそれぞれ含まれる操作量成分u3_x,u3_yは、制御用目標重心速度Vb_x_mdfd,Vb_y_mdfdに重心速度推定値Vb_x_s,Vb_y_sを各々収束させるように決定される。
この結果、乗員が車両1に付加した推進力による車両系重心点の実際の移動速度の増速(概略前後方向での増速)が該推進力による要求に即して速やかに行なわれるように、車輪体5の移動速度が制御されることとなる。従って、車両1が付加された推進力によって円滑に加速することとなる。
なお、速度追従モードにおいて、車両1に制動力を付加することによって、図19のステップS24−1の判断結果が肯定的になると(減速要求が発生すると)、演算処理モードが前記制動モードに変化する。このため、車両1の移動速度が減衰していくこととなる。この場合、減速要求の発生中は、ステップS23−6の第2制動演算処理(図18の処理)によって、|↑Vb_aim1|とθvb_aim1とが決定されるので、基本要求重心速度ベクトル↑Vb_aim1、あるいはこれに追従する要求重心速度ベクトル↑Vb_aimは、その向きを一定に保ったまま、大きさが一定の時間的変化率(前記所定値ΔVb2により規定される時間的変化率)で減衰していくように決定されることとなる。
次に、速度追従モードにおいて、車両1への推進力の付加が終了し、推定重心速度絶対値変化率DVb_sが前記第2閾値DV2よりも小さくなると(図19のステップS24−3の判断結果が肯定的になると)、図19のステップS24−5の処理によって、演算処理モードが速度追従モードから速度ホールドモードに変更されることとなる。
この速度ホールドモードでは、加速要求及び減速要求が発生しない状況(図19のステップS25−1、25−4の判断結果がいずれも否定的となる状況)で、カウントダウンタイマの計時が終了するまで、基本要求重心速度ベクトル↑Vb_aim1は、前回値の速度ベクトル↑Vb_aim1_pと同じ速度ベクトルに設定される。
従って、速度ホールドモードの開始後、カウントダウンタイマの計時が終了するまでの所定時間(カウントダウンタイマの初期値Tmの時間)の期間において、基本要求重心速度ベクトル↑Vb_aim1は、速度ホールドモードが開始する直前に決定された速度ベクトルと同じ速度ベクトルで一定に保持されることとなる。
このため、↑Vb_aim1に追従するように決定される要求重心速度ベクトル↑Vb_aimも、一定の速度ベクトル(速度ホールドモードが開始する直前に決定された↑Vb_aimと一致又はほぼ一致する速度ベクトル)に保たれるように決定されることとなる。
そして、上記の如く決定される要求重心速度ベクトル↑Vb_aimのX軸方向成分とY軸方向成分とが前記制御用目標重心速度Vb_x_mdfd,Vb_y_mdfdとして決定される。さらに、各仮想車輪回転角加速度指令ωdotw_x_cmd,ωdotw_y_cmdにそれぞれ含まれる操作量成分u3_x,u3_yは、制御用目標重心速度Vb_x_mdfd,Vb_y_mdfdに重心速度推定値Vb_x_s,Vb_y_sを各々収束させるように決定される。
この結果、車両1の増速後に、カウントダウンタイマの計時が終了するまでの期間(前記初期値Tmにより表される時間の期間)において、乗員の上体の姿勢の頻繁な調整を必要とすることなく、車両系重心点の実際の速度ベクトル↑Vbの大きさ及び向きが一定に保たれるように該車輪体5の移動速度が制御されることとなる。従って、この状況での車両1の実際の走行状態は、乗員がその上体を積極的に動かすような操縦操作をせずとも、ほぼ一定の速度ベクトルで滑走するような状態となる。
なお、速度ホールドモードにおいて、再び車両1に概略前後方向の推進力を付加することによって、図20のステップS25−4の判断結果が肯定的になると(加速要求が発生すると)、演算処理モードが前記速度追従モードに戻る。このため、車両1が再び概略前後方向で加速することとなる。
また、速度ホールドモードにおいて、車両1に制動力を付加することによって、図20のステップS25−1の判断結果が肯定的になると(減速要求が発生すると)、演算処理モードが前記制動モードに変化する。このため、車両1の移動速度が減衰していくこととなる。この場合、速度追従モードで減速要求が発生した場合と同様に、減速要求の発生中は、ステップS23−6の第2制動演算処理(図18の処理)によって、|↑Vb_aim1|とθvb_aim1とが決定される
次に、速度ホールドモードにおいて、加速要求及び減速要求が発生しない状況(図20のステップS25−1、25−4の判断結果がいずれも否定的となる状況)が保たれたまま、前記カウントダウンタイマの計時が終了すると、図20のステップS25−11の処理によって、演算処理モードが速度ホールドモードから制動モードに変更されることとなる。
この制動モードでは、加速要求及び減速要求が発生しない状況(図15のステップS23−1,23−4の判断結果がいずれも否定的となる状況)で、図16のステップS23−5−1,23−5−2の処理が制御処理周期毎に実行されることによって、基本要求重心速度ベクトル絶対値|↑Vb_aim1|が、“0”まで一定の時間的変化率(前記ΔVb1により規定される時間的変化率)で連続的に減衰していくこととなる。そして、|↑Vb_aim1|が“0”まで減衰した後には、|↑Vb_aim1|が“0”に保持される。
さらに、制動モードでは、加速要求及び減速要求が発生しない状況で、図16のステップS23−5−3以降の処理が制御処理周期毎に実行される。この場合、速度ホールドモードから制動モードへの移行直前に決定された基本要求重心速度ベクトル↑Vb_aim1の向き(図20のステップS25−8の判断結果が否定的となる制御処理周期の1つ前の制御処理周期で決定された↑Vb_aim1の向き)が、X軸方向と異なり、且つ、比較的X軸方向に近い向きである場合(より正確には、移行直前に決定された↑Vb_aim1の方位角θvb_aim1が、0°<θvb_aim1_p≦θth1+、θth1-≦θvb_aim1_p<0°、θth2+≦θvb_aim1_p<180°、−180°<θvb_aim1_p≦θth2-のいずれかの範囲内の角度値である場合)には、|↑Vb_aim1|が“0”に減衰するまでの期間内において、θvb_aim1が一定の時間的変化率で収束目標角度値としての0°又は180°又は−180°に近づいていき、最終的に、該収束目標角度値に保持される。従って、制動モードの開始後、|↑Vb_aim1|が“0”に減衰するまでの期間内において、基本要求重心速度ベクトル↑Vb_aimの向きが、連続的にX軸方向に近づいていくこととなる。換言すれば、当該期間内において、基本要求重心速度ベクトル↑Vb_aimのX軸方向成分Vb_x_aim1の絶対値に対するY軸方向成分Vb_y_aim1の絶対値の比率が、“0”に近づいていくこととなる。そして、|↑Vb_aim1|が“0”に減衰するまでに、↑Vb_aim1の向きがX軸方向と同じ向きに到達すると(Vb_y_aim1=0になると)、↑Vb_aim1の向きがX軸方向と同じ向きに保持される。
従って、↑Vb_aim1は、その大きさが減衰しつつ、その向きがX軸方向に近づく(収束する)ように決定される。このように↑Vb_aim1が決定される状況では、↑Vb_aim1に追従するように決定される要求重心速度ベクトル↑Vb_aimも、その大きさが減衰しつつ、その向きがX軸方向に近づくように決定されることとなる。
また、速度ホールドモードから制動モードへの移行直前に決定された基本要求重心速度ベクトル↑Vb_aim1の向きが、X軸方向と異なり、且つ、X軸方向から比較的乖離した向きである場合(より正確には、上記移行直前に決定された↑Vb_aim1の方位角θvb_aim1が、θth1+<θvb_aim1_p<θth2+、θth2-<θvb_aim1_p<θth1-のいずれかの範囲内の角度値である場合)には、|↑Vb_aim1|が“0”に減衰するまでの期間内において、θvb_aim1が、上記移行直前に決定された↑Vb_aim1の方位角θvb_aim1と同じ角度値で一定に保持される。
従って、↑Vb_aim1は、その大きさが減衰しつつ、その向きが一定に保たれるように決定される。このように↑Vb_aim1が決定される状況では、↑Vb_aim1に追従するように決定される要求重心速度ベクトル↑Vb_aimも、その大きさが減衰しつつ、その向きが一定に保たれるように決定されることとなる。
なお、本実施形態では、速度ホールドモードでは、↑Vb_aim1は、その大きさ及び向きが一定に保持されるので、速度ホールドモードから制動モードへの移行直前に決定された基本要求重心速度ベクトル↑Vb_aim1は、結果的には、速度追従モードから速度ホールドモードへの移行直前に決定された↑Vb_aim1(本実施形態では、図19のステップS24−3の判断結果が肯定的となる制御処理周期で決定された↑Vb_aim1)に一致する。
そして、制動モードで上記の如く決定される要求重心速度ベクトル↑Vb_aimのX軸方向成分とY軸方向成分とが前記制御用目標重心速度Vb_x_mdfd,Vb_y_mdfdとして決定される。さらに、各仮想車輪回転角加速度指令ωdotw_x_cmd,ωdotw_y_cmdにそれぞれ含まれる操作量成分u3_x,u3_yは、制御用目標重心速度Vb_x_mdfd,Vb_y_mdfdに重心速度推定値Vb_x_s,Vb_y_sを各々収束させるように決定される。
この結果、前記制動モードの前の演算処理モードが速度ホールドモードである場合に、乗員がその上体の動きによる積極的な操縦操作を行なわずとも、車両系重心点の実際の速度ベクルの大きさが速度ホールドモードでの大きさから連続的に減衰していくように車輪体5の移動速度が制御されることとなる。
この場合、速度ホールドモードから制動モードへの移行直前に決定された↑Vb_aim1(=速度追従モードから速度ホールドモードへの移行直前に決定された↑Vb_aim1)の向きが、X軸方向と異なり、且つ、比較的X軸方向に近い向きである場合には、乗員がその上体の動きによる積極的な操縦操作を行なわずとも、車両系重心点の速度ベクトルの大きさが減衰しながら、該速度ベクトルの向きが自動的にX軸方向(乗員の前後方向)に近づいていくこととなる。従って、乗員の前後方向に対する車両1の直進性が高まる。
ここで、車両1を加速しようとする場合、多くの場合、特に乗員の前後方向で車両1を加速することが要求される。この場合、本実施形態の車両1は、上記の如く前後方向に対する直進性が高いので、車両1に付加する推進力の向きが、前後方向から多少ずれていても、その後の速度ホールドモードに続く制動モードにおいて、自動的に車両系重心点の速度ベクトルが、前後方向に向くように車輪体5の移動速度が制御される。
このため、車両1の移動方向のばらつきが生じにくく、乗員の前後方向に対して高い直進性を有する車両1(乗員の前後方向に進み易い車両1)が実現されている。ひいては、車両1を前後方向に移動させる場合に、車両1に付加する推進力を正確に前後方向に向けずとも、該車両1を前後方向に移動させることができる。その結果、車両1を前後方向に移動させる操縦操作が容易になる。
また、速度ホールドモードから制動モードへの移行直前に決定された基本要求重心速度ベクトル↑Vb_aim1(=速度追従モードから速度ホールドモードへの移行直前に決定された↑Vb_aim1)の向きが、X軸方向と異なり、且つ、X軸方向から比較的乖離した向きである場合には、乗員がその上体の動きによる積極的な操縦操作を行なわずとも、車両系重心点の速度ベクトルの大きさが減衰しながら、該速度ベクトルの向きがほぼ一定に保たれる。すなわち、速度ホールドモードから制動モードへの移行直前に決定された↑Vb_aim1の向きが、X軸方向から比較的乖離した向きである場合には、速度追従モードで最終的に乗員が意図した車両系重心の速度ベクトルの向きが、X軸方向と異なる向きである可能性が高い。従って、速度追従モード以後に、車両系重心点が、乗員が意図する方向と乖離する向きに移動しようとするのを防止できる。
以上が、本実施形態による倒立振子型移動体の前提となる全方向移動車両1の詳細である。
2.連携技術
以下では、上述した全方向移動車両の移動制御に関する基本的機能を前提として、本発明の実施形態による倒立振子型移動体が他の移動体と連携動作するための技術内容について、一対の松葉杖のように機能する歩行補助装置に適用した場合を例として説明する。
なお、複数台の移動体が連携動作する形態としては、例えば、各移動体間の位置(距離)を一定に保持するように連携動作する形態、各移動体の姿勢(傾斜角度)を合わせるように連携動作する形態、これらを組み合わせた連携動作の形態等が挙げられる。ここでは、説明の便宜上、各移動体の姿勢(傾斜角度)を一致させるように連携動作する形態を中心に説明し、その他の連携動作の形態については変形例として説明することとする。
図21に、本実施形態による倒立振子型移動体が適用された歩行補助装置1000の使用形態を示す。同図において、歩行補助装置1000は、2台の倒立振子型移動体1000R,1000Lから構成される。このうち、倒立振子型移動体1000Rは、利用者に対して右側(紙面左側)に位置する松葉杖に対応し、倒立振子型移動体1000Lは、この利用者に対して左側(紙面右側)に位置する松葉杖に対応している。これら倒立振子型移動体1000R,1000Lのそれぞれは、前提技術として前述した全方向移動車両1と同様に、倒立振子モデルの動力学に基づく重心運動の制御を実施することにより、自立した状態で床面上を移動するように構成されている。
倒立振子型移動体1000R,1000Lのそれぞれは、連携動作に関する機能を除けば、基本的には前述の全方向移動車両1と同様に構成されるが、歩行補助装置としての機能の必要上、図1に示すシート3、グリップ17R,17L、ステップ25R,25L等に代えて、図21に示す脇当部1000aおよび把持部1000bが設けられている。後述するように、これら2台の倒立振子型移動体1000R,1000Lが互いに連携動作することにより、歩行補助装置1000が松葉杖のように機能し、利用者の歩行を補助するものとなっている。
なお、以下では、2台の倒立振子型移動体1000R,1000Lのうち、利用者の右側に位置する倒立振子型移動体1000Rに着目して各構成要素を説明するものとし、特に明示しない限り、各要素は、倒立振子型移動体1000Rの要素を表すものとする。ただし、左右の概念を入れ替えれば、利用者の左側に位置する倒立振子型移動体1000Lは倒立振子型移動体1000Rと同じである。
また、適宜、2台の倒立振子型移動体1000R,1000Lのうち、本実施形態の説明で着目する構成要素が設けられた倒立振子型移動体1000Rを「自移動体1000R」と称し、この倒立振子型移動体1000Rと連携動作する他方の倒立振子型移動体1000Lを「他移動体1000L」と称す。ここで、「自移動体」なる用語は、連携動作を制御するための構成要素から見て、その構成要素が備えられた制御対象の移動体を意味する。例えば、倒立振子型移動体1000Rに備えられた姿勢制御演算部80から見れば、自移動体は倒立振子型移動体1000Rを意味し、また、倒立振子型移動体1000Lに備えられた姿勢制御演算部80から見れば、自移動体は倒立振子型移動体1000Lを意味する。また、「他移動体」なる用語は、「自移動体」と連携動作する関係にある他の移動体を意味する。このような「自移動体」および「他移動体」の用語の意味は、本発明が規定される請求項においても同様である。
次に、前述の図9および図13を参照して、連携動作に関連する構成を説明する。
図9において、制御ユニット50は、前述の図7に示すステップS9の車両制御演算処理において連携動作に関する処理に必要とされる要素として、通信部701(情報取得手段)、偏差演算部702,703,705、距離センサ704(情報取得手段)、外力センサ706(外力検知手段)、荷重センサ707(荷重検知手段)、ゲイン調整部708を備えている。
ここで、通信部701は、他移動体1000Lから、他移動体1000Lの現在の状態を表す状態情報を取得するためのものであり(情報取得手段)、ここでは、上記状態情報として、他移動体1000Lの現在の姿勢を表す姿勢情報を取得する。具体的には、通信部701は、上記姿勢情報として、他移動体1000Lの取得鉛直方向に対する基体傾斜角度計測値θb_xy_s(L)と基体傾斜角速度計測値θbdot_xy_s(L)とを取得する。なお、通信部701は、無線と有線を問わず、他移動体1000Lから情報を取得することができるものであれば、どのような手段により実現されてもよい。
特に図示していないが、通信部701により、自移動体1000Rの状態(姿勢)を表す状態情報として、自移動体1000Rの基体傾斜角度計測値θb_xy_s(以下、基体傾斜角度計測値θb_xy_s(R)と表記する)と、基体傾斜角速度計測値θbdot_xy_s(以下、基体傾斜角速度計測値θbdot_xy_s(R)と表記する)とが他移動体1000Lに送信され、他移動体1000Lにおける同様の処理に使用される。
偏差演算部702は、自移動体1000Rの基体傾斜角度計測値θb_xy_s(R)と他移動体1000Lの基体傾斜角度計測値θb_xy_s(L)との偏差(θb_xy_s(R)−θb_xy_s(L))を算出するためのものである。偏差演算部703は、自移動体1000Rの基体傾斜角速度計測値θbdot_xy_s(R)と他移動体1000Lの基体傾斜角速度計測値θbdot_xy_s(L)との偏差(θbdot_xy_s(R)−θbdot_xy_s(L))を算出するためのものである。
距離センサ704は、他移動体1000Lの現在の状態を表す状態情報として、自移動体1000Rに対する他移動体1000Lの位置を表す位置情報(または、他移動体1000Lに対する自移動体1000Rの位置を表す位置情報)を取得するためのものである(情報取得手段)。この距離センサ704は、上記位置情報として、図22(a)に示すように、自移動体1000Rと他移動体1000Lとの間の距離Ls_xyを計測する。偏差演算部705は、距離センサ704により測定された距離Ls_xyと、予め設定された目標距離Ltarget_xyとの間の偏差(Ls_xy−Ltarget_xy)を算出するためのものである。
目標距離Ltarget_xyは、利用者によって適宜設定されるものとしてもよい。また、距離センサ704は、赤外線センサ等の公知技術を用いて実現されてもよく、自移動体1000Rに対する他移動体1000Lの位置を把握することができれば、どのような手段により実現されてもよい。
なお、距離センサ704および偏差演算部705は、移動体間の距離を一定に維持する後述の第1の変形例における連携動作において必要とされるものであり、距離を一定に維持する必要のない連携動作の場合には省略してもよい。
外力センサ706は、自移動体1000Rの基体9に作用する水平方向の外力Fhandを検知するためのものである(外力検知手段)。すなわち、この外力センサ706は、利用者が倒立振子型移動体1000Rの把持部1000bに意図的に加える水平方向の外力Fhandを検知する。利用者が把持部1000bに水平方向に外力Fhandを加えることにより、後述するように、利用者により倒立振子型移動体1000Rに移動方向を指示することが可能とされている。
荷重センサ707は、自移動体1000Rの基体9に作用する鉛直方向の荷重Fweightを検知するためのものである(荷重検知手段)。この荷重センサ707は、倒立振子型移動体1000Rの上端に配置された脇当部1000aに装着されており、この脇当部1000aに上方から加えられる利用者の体重による鉛直方向の荷重Fweightを検知する。
ゲイン調整部708は、倒立振子型移動体1000Rの移動制御におけるフィードバックのゲイン調整パラメータKw_xyを生成するためのものである。このゲイン調整パラメータKw_xyは、荷重センサ707により検知された荷重Fweightの増大に対し、後述の図13に示すゲイン係数K6_xyを減少させるためのものであるが、このゲイン係数K6_xyの意味の詳細については後述する。
上述の偏差演算部702により算出された基体傾斜角度計測値の偏差(θb_xy_s(R)−θb_xy_s(L))と、偏差演算部703により算出された基体傾斜角速度計測値の偏差(θbdot_xy_s(R)−θbdot_xy_s(L))と、偏差演算部705により算出された距離の偏差(Ls−Ltarget)と、外力センサ706により検知された外力Fhandの値と、ゲイン調整部708により生成されたゲインKw_xyとが、他移動体1000Lとの連携動作を制御するための入力値として姿勢制御演算部80(移動制御手段)に与えられる。
姿勢制御演算手段80は、上記の状態情報に基づき自移動体1000Rの移動を制御する移動制御手段を構成し、他移動体1000Lの現在の状態に対する自移動体1000Rの状態が所定の条件を満たすように移動を制御する。ここで、所定の条件は、自移動体と他移動体とを連携動作させるために規定された条件であり、具体的には、後述する式071x,071yの右辺の第4項ないし第7項によって与えられる条件である。
続いて、図13を参照して、連係動作に関する姿勢制御演算部80の構成を説明する。
本実施形態による倒立振子型移動体の前提となる全方向移動車両の説明で述べたように、姿勢制御演算部80は、各入力値に対して、仮想車輪回転角加速度指令ωwdot_x_com,ωwdot_y_comを算出するが、ここでは、連携動作に関与する入力値に対する姿勢制御演算部80の処理を説明する。ただし、全方向移動車両の基本的な移動制御に関する前述の式07xおよび式07yに含まれる項に関わる説明は省略する。
姿勢制御演算部80は、入力値を用いて、次式071x,071yにより、連携動作に関与する入力値をフィードバック制御の操作量成分として考慮に入れて、仮想車輪回転角加速度指令ωdotw_xy_comを算出する。
ωwdot_x_cmd=K1_x・θbe_x_s+K2_x・θbdot_x_s
+K3_x・(Vb_x_s−Vb_x_mdfd)
+K4_x・(θb_x_s(R)−θb_x_s(L))
+K5_x・(θbdot_x_s(R)−θbdot_x_s(L))
+K6_x・Fhand_x
+K7_x・(Ls_x−Ltarget_x) ……式071x
ωwdot_y_cmd=K1_y・θbe_y_s+K2_y・θbdot_y_s
+K3_y・(Vb_y_s−Vb_y_mdfd)
+K4_y・(θb_y_s(R)−θb_y_s(L))
+K5_y・(θbdot_y_s(R)−θbdot_x_s(L))
+K6_y・Fhand_y
+K7_y・(Ls_y−Ltarget_y) ……式071y
式071xおよび式071yにおいて、右辺第1項ないし第3項は、前述の式071xおよび式071yと同一であり、右辺第4項ないし第7項が新たに付加されている。
ここで、右辺第4項は、基体傾斜角度計測値の偏差(θb_xy_s(R)−θb_xy_s(L))に対応する項であり、右辺第5項は、基体傾斜角速度計測値の偏差(θbdot_xy_s(R)−θbdot_xy_s(L))に対応する項である。これら右辺第4項および第5項を0に収束させるように倒立振子型移動体1000Rの移動制御が実行される結果、自移動体1000Rと他移動体1000Lとの姿勢(傾斜角度)が一致するようになる。
また、右辺第6項は、把持部1000bを介して基体9に作用する平行方向の外力Fhand_xyに対応する項である。この第6項を0に収束させるようにフィードバック制御が実行される結果、自移動体1000Rは、基体9に加えられる外力Fhand_xyが減少する方向、即ち外力Fhandの方向と略一致する方向に移動する。従って、利用者は、把持部1000bに水平方向の外力Fhand_xyを意図的に加えることにより、倒立振子型移動体1000Rの移動方向を指示することができる。
右辺第7項は、距離の偏差(Ls_xy−Ltarget_xy)に対応する項であり、自移動体1000Rと他移動体1000Lとの間の距離を一定に維持させるための項である。この第7項を0に収束させるように倒立振子型移動体1000Rの移動制御が実行される結果、自移動体1000Rと他移動体1000Lとの間の距離Ls_xyが目標距離Ltarget_xyに維持されることになる。
前述のフィードバック制御則に従って、上述の第1項ないし第7項の全てが0に収束するように仮想車輪回転角加速度指令ωdotw_xy_comが決定される。その結果、自移動体1000Rは、自身の移動制御において、他移動体1000Lとの間の距離を一定に維持しながら、自移動体1000Rの傾斜角度を他移動体1000Lの傾斜角度と一致させるように、自移動体1000Lと連携動作することとなる。他移動体1000Lについても同様のことが言える。
この場合、式071xにおける各操作量成分に係わるゲイン係数K1_x,K2_x,K3_x,K4_x,K5_xは、ゲイン調整パラメータKr_xに応じて可変的に設定され、式071yにおける各操作量成分に係わるゲイン係数K1_y,K2_y,K3_y,K4_y,K5_yは、ゲイン調整パラメータKr_yに応じて可変的に設定される。具体的には、前述した全方向移動車両におけるゲイン係数K1_x,K2_x,K3_xと同様に、ゲイン係数K4_x,K5_xは、Kr_x=0である場合には、K4_x=K4_a_x,K5_x=K5_a_xとなり、Kr_x=1である場合には、K4_x=K4_b_x,K5_x=K5_b_xとなる。従って、Kr_xが“0”から“1”に近づくに伴い、第i(i=1,2,3,4,5)ゲイン係数Ki_xはKi_a_xからKi_b_x近づいていく。
式071yの演算に用いる第iゲイン係数Ki_y(i=1,2,3,4,5)についても同様である。
また、ゲイン係数K6_xyは、ゲイン調整パラメータKw_xyに応じて可変的に設定され、荷重Fweightの増大に対して減少する傾向を示すように制御される。例えば、荷重Fweightが存在しない場合(すなわち、利用者が脇当部1000aに体重を加えていない場合)には、K6_xy=K6_a_xyとなり、利用者が脇当部1000aに体重を徐々に加えると、荷重Fweightの増大に伴ってゲイン係数K6_xyがK6_b_xyに向かって減少する。このように、フィードバック制御におけるゲイン係数K6_xyを調整することにより、利用者が自身の体重を脇当部1000aに加えた状態では、外力Fhandの変化(すなわち、利用者による操作量の変化)に対して倒立振子型移動体1000Rの移動量の変化が抑制される。従って、この場合、倒立振子型移動体1000Rの挙動が穏やかになり、歩行補助装置としての操作性を向上させることができる。
ゲイン係数K7_xyは定数であり、自移動体1000Rと他移動体1000Lとの間の距離を一定に維持するためのフィードバック制御の応答性を定める。すなわち、ゲイン係数K7_xyを増大させると、他移動体1000Lとの間の距離を一定に維持するためのフィードバック制御における応答性が高まり、逆に、ゲイン係数K7_xyを減少させると、上記の応答性が低下する。従って、ゲイン係数K7_xyは、上記の応答性を考慮して適切に設定されるが、利用者が調整可能なものとしてもよい。
本実施形態において、自移動体1000Rの移動のフィードバック制御は、式071x,071yの右辺の各項を“0”に収束させるように実施される。このことは、式071x,071yの右辺の各項によって規定される条件が満たされるようにフィードバック制御が実施されることを意味している。本実施形態では、連携動作の条件(所定の条件)は、式071x,071yの右辺の一部の項によって表現されている。
すなわち、式071xの右辺には、連携動作の条件を規定する項として、第4項(K4_x・(θb_x_s(R)−θb_x_s(L)))、第5項(K5_x・(θbdot_x_s(R)−θbdot_x_s(L)))、第7項(K7_x・(Ls_x−Ltarget_x))が含まれている。このうち、第4項および第5項は、自移動体1000Rの姿勢を他の移動体1000Lの姿勢に合わせるための条件に対応し、第7項は、自移動体1000Rと他移動体1000Lとの間の距離を一定に維持するための条件に対応する。式071yについても同様である。
なお、上述の式071xおよび式071yの右辺第7項は、移動体間の距離を一定に維持する後述の第1の変形例において意味を有するものであるが、ここでは、移動体間の距離については連携動作の対象外としているため、上述の式071xおよび式071yの右辺第7項は不要である。そこで、強制的にK7_xy=0に設定し、式071xおよび式071yにおける右辺第7項を無視するものとする。
姿勢制御演算部80は、前述したたゲイン係数K1_x,K2_x,K3_xと共に、上記の如く決定し連携動作に関与するゲイン係数K4_x,K5_x,K6_x,K7_xを用いて前記式071xの演算を行なうことで、連携動作を考慮したX軸方向に輪転する仮想車輪62_xに係わる仮想車輪回転角加速度指令ωwdot_x_cmdを算出する。
ここで、連携動作に関与する構成について、図13を参照して更に詳細に説明する。
同図に示すように、姿勢制御演算部80は、連携動作に関与する要素として、処理部80g〜80jを備えている。そして、姿勢制御演算部80は、基体傾斜角度偏差計測値の偏差(θb_x_s(R)−θbe_x_s(L))に第1ゲイン係数K4_xを乗じてなる操作量成分u4_xと、基体傾斜角速度計測値の偏差(θbdot_x_s(R)−θbdot_x_s(L))に第5ゲイン係数K5_xを乗じてなる操作量成分u5_xと、外力Fhandの値に第6ゲイン係数K6_xを乗じてなる操作量成分u6_xを、それぞれ、処理部80g,80h,80iで算出する。また、姿勢制御演算部80は、距離の偏差(Ls−Ltarget)にゲイン係数K7_xを乗じてなる操作量成分u7_xを処理部80jで算出する。そして、姿勢制御演算部80は、前述の移動制御に関与する基本的な操作量成分u1_x,u2_x,u3_xと共に、連携動作に関与する操作量成分u4_x,u5_x,u6_x,u7_xを演算部80eにて加え合わせることにより、仮想車輪回転角加速度指令ωwdot_x_comを算出する。
同様に、姿勢制御演算部80は、前述したゲイン係数K1_y,K2_y,K3_yと共に、上記の如く決定した連携動作に関与するゲイン係数K4_y,K5_y,K6_y,K7_yを用いて前記式071yの演算を行なうことで、連携動作を考慮したY軸方向に輪転する仮想車輪62_yに係わる仮想車輪回転角加速度指令ωwdot_y_cmdを算出する。
姿勢制御演算部80は、上記の如く、連携動作を考慮した仮想車輪回転角加速度指令ωwdot_x_com,ωwdot_y_comを算出した後、次に、これらのωwdot_x_com,ωwdot_y_comをそれぞれ積分器80fにより積分することによって、前記仮想車輪回転速度指令ωw_x_com,ωw_y_comを決定する。
以上が連携動作に関与する姿勢制御演算部80の詳細である。その他については、前述の全方向移動車両と同様である。
次に、本実施形態による倒立振子型移動体1000Rおよび倒立振子型移動体1000Lの動作(歩行補助装置としての使用形態)を説明する。
図21に示すように、右足を負傷している利用者は、静止した状態で、右側に倒立振子型移動体1000Rを位置させ、左側に倒立振子型移動体1000Lを位置させ、脇当部1000aを自己の脇に当てるようにして把持部1000bを把持する。そして、自己の体重を脇当部1000aに徐々に加えて、自己の体重を倒立振子型移動体1000R,1000Lで支えるようにする。
この静止状態から、利用者が健常な左足で床を蹴ると、この左足が床から受ける反力により一対の倒立振子型移動体1000R,1000Lが前方に向かって加速され(速度追従モード)、移動を開始する。そして、一定時間にわたって移動速度が一定に保たれ(速度ホールドモード)、その後、徐々に減速される(制動モード)。必要に応じて、利用者が床面を再び蹴れば、同様にして移動が再開される。
ここで、倒立振子型移動体1000R,1000Lは、移動制御の過程で、お互いの姿勢に関する情報を交換し、上述した連携動作に関する処理を実行することにより、お互いの姿勢(傾斜角度)を一致させるように連携動作する。この場合、2台の倒立振子型移動体1000R,1000Lは、互いに並行状態に維持される。したがって、利用者にとって2台の倒立振子型移動体1000R,1000Lが一対の松葉杖と同様に歩行補助装置として機能し、倒立振子型移動体1000R,1000Lが利用者の歩行を補助する。
なお、全方向移動車両の説明で述べたように、利用者が基体9を基本姿勢から傾けるように操作によって、倒立振子型移動体1000R,1000Lを発進させ、停止させることも可能である。
[第1の変形例]
上述の実施形態では、連携動作の条件として、自移動体1000Rと他移動体1000Lの姿勢(傾斜角度)を一致させることを設定したが、この第1の変形例では、連携動作の条件として、図22(a)に示すように、自移動体1000Rと他移動体1000Rと間の距離Ls_xyを一定の目標距離Ltarget_xyに維持することを設定する。この連携動作の条件は、上述の式071x、071yの第7項によって与えられている。
上述の実施形態では、説明の便宜上、右辺第7項を“0”に設定し、この第7項を無視するものとしたが、この第1の変形例では、式071x,071yにおける右辺第7項のゲイン係数K7_xyを、フィードバック制御の応答性を考慮した“0”以外の適切な定数に設定する。
なお、歩行補助装置としての機能を考慮すれば、Y軸方向の目標距離Ltarget_yを利用者の体形に応じた適切な値に設定し、X軸方向の目標距離Ltarget_xを“0”に設定することが望ましい。
このように式071x,071yが設定されれば、自移動体1000Rは、上述した姿勢(傾斜角度)を合わせる連携動作に加えて、他移動体1000Lとの間の距離L_xyを一定の目標値Ltarget_xyに維持するように連携動作するようになる。すなわち、姿勢制御演算部80は、距離センサ704によって取得された距離Ls_xyが、連携動作の条件として与えられた目標値Ltarget_xyに収束するように、すなわち、式071x,071yの右辺第7項が“0”に収束するように、フィードバック制御を実施する。その結果、連携動作の条件が満たされ、自移動体1000Rと他移動体1000Lとの間の距離Ls_xyが一定の目標距離Ltarget_xyに維持される。
この第1の変形例によれば、自移動体1000Rと他移動体1000Lとの間の距離が一定に維持されるので、利用者は、自移動体1000Rと体移動体1000Lとの間隔を修正するための操縦操作に煩わされることがなくなり、従って歩行補助装置としての操作性と利便性を向上させることができる。
なお、第1の変形例では、歩行補助装置としての使用形態を考慮して、X軸およびY軸の各方向の距離を目標値に維持するものとしたが、自移動体1000Rを中心として一定の範囲内に他移動体1000Lが位置するようにしてもよい。
[第2の変形例]
上述の実施形態では、自移動体1000Rの姿勢(傾斜角度)と他移動体1000Lの姿勢(傾斜角度)とを一致させる(姿勢を相互に合わせる)ものとしたが、この第2の変形例では、連携動作の条件として、各移動体の姿勢を設定する。具体的には、連携動作の条件として、図21(b)に示すように、移動体間のX軸周りの傾斜角度の差分(θb_y_s(R)−θb_y_s(L))が所定値となることを設定する。この場合、前述の式071yは次式072yのように表現し直される。なお、ここでは、X軸周りの傾斜角度の差分のみを所定値に維持するものとし、Y軸回りについては上述の実施形態と同様とする。
ωwdot_y_cmd=K1_y・θbe_y_s+K2_y・θbdot_y_s
+K3_y・(Vb_y_s−Vb_y_mdfd)
+K4_y・(θb_y_s(R)−θb_y_s(L)−θc_y)
+K5_y・(θbdot_y_s(R)−θbdot_y_s(L))
+K6_y・Fhand_y
+K7_y・(Ls_y−Ltarget_y) ……式072y
ここで、式071yと式072yの違いは、傾斜角度の偏差に関連する右辺第4項の括弧内に定数「−θc_y」が付加されている点である。この定数「−θc_y」は、2台の倒立振子型移動体1000R,1000LのX軸周りの傾斜角度の差分に対応する所定値である。
このように定数「−θc_y」を含む式072yが設定されれば、第4項を0に収束させるように倒立振子型移動体1000Rの移動制御が実行される結果、自移動体1000Rと他移動体1000Lとの傾斜角度の差分(θb_y_s(R)−θb_y_s(L))が所定値「θc_y」に一致するようになる。従って、2台の倒立振子型移動体1000R,1000LのX軸周りの傾斜角度の差分を所定値「θc_y」に維持するようして、2台の倒立振子型移動体1000R,1000Lの連携動作が行われることになる。
なお、Y軸周りの傾斜角度の差分を所定値に維持するものとする場合には、前述の式071yの第4項に、所定値に対応する定数「−θc_x」を付加すればよい。
第2の変形例によれば、自移動体1000Rの傾斜角度と他移動体1000Lの傾斜角度の差分を所定値に維持するようにしたので、自移動体1000Rと他移動体1000Lとを相互に異なる方向に傾倒させた状態に維持することが可能になる。従って、例えば、一対の倒立振子型移動体の下端側の間隔を上端側の間隔よりも広くした状態に維持させることができ、歩行補助装置としての安定性を向上させることができる。
[第3の変形例]
第3の変形例は、連携動作の条件として、各移動体の姿勢を設定する点では第2の変形例と共通するが、第3の変形例では、各移動体の姿勢を個々に設定する。具体的には、連携動作の条件として、図21(c)に示すように、自移動体1000Rの傾斜角度θb_y_s(R)を所定の目標傾斜角度「+θtarget」に維持し、他移動体1000Lの傾斜角度θb_y_s(L)を所定の目標傾斜角度「−θtarget」に維持することを設定する。この場合、自移動体1000Rにおける前述の式071yは次の式073y(R)のように表現し直され、他移動体1000Lにおける前述の式071yは次の式073y(L)のように表現し直される。なお、ここでは、説明の便宜上、X軸周りの傾斜角度のみを目標傾斜角度に維持するものとし、Y軸回りについては上述の実施形態と同様とする。
ωwdot_y_cmd=K1_y・θbe_y_s+K2_y・θbdot_y_s
+K3_y・(Vb_y_s−Vb_y_mdfd)
+K4_y・(θb_y_s(R)−θtarget)
+K5_y・(θbdot_y_s(R)−θbdot_y_s(L))
+K6_y・Fhand_y
+K7_y・(Ls_y−Ltarget_y) ……式073y(R)
ωwdot_y_cmd=K1_y・θbe_y_s+K2_y・θbdot_y_s
+K3_y・(Vb_y_s−Vb_y_mdfd)
+K4_y・(θb_y_s(L)+θtarget)
+K5_y・(θbdot_y_s(R)−θbdot_y_s(L))
+K6_y・Fhand_y
+K7_y・(Ls_y−Ltarget_y) ……式073y(L)
ここで、式071yと式073y(R)の違いは右辺第4項であり、式071yと式073y(L)の違いについても同様である。式073y(R)によれば、第4項を0に収束させるようにフィードバック制御が実行される結果、自移動体1000Rの傾斜角度(θb_y_s(R))が目標傾斜角度の「+θtarget」に一致するようになる。一方、式074y(L)によれば、第4項を0に収束させるようにフィードバック制御が実行される結果、他移動体1000Lの傾斜角度(θb_y_s(L))が目標傾斜角度の「−θtarget」に一致するようになる。従って、2台の倒立振子型移動体1000R,1000LのX軸周りの各傾斜角度を、それぞれの目標傾斜角度に維持するように、2台の倒立振子型移動体1000R,1000Lの連携動作が行われる。
第3の変形例によれば、2台の倒立振子型移動体1000R,1000LのX軸周りの各傾斜角度が、それぞれの目標傾斜角度に維持されるので、利用者は、2台の倒立振子型移動体1000R,1000Lの傾斜角度の修正に煩わされることがなくなり、歩行補助装置としての操作性と安定性を向上させることができる。
なお、上述の第1の変形例と第3の変形例とを組み合わせてもよい。この場合、自移動体1000Rと他移動体1000Lとの間の距離が一定に維持されると共に、各移動体の傾斜角度が一定に維持されるので、歩行補助装置としての安定性を更に向上させることができる。
以上、本発明の実施形態および変形例を説明したが、本発明は上述の実施形態および変形例に限定されるものではなく、本発明の要旨を逸脱しない範囲で種々の変形が可能である。例えば、上述の実施形態において、荷重センサ704として6軸センサを採用するものとすれば、この6軸センサの各成分を分析することにより、脇当部1000aに加えられる体重による荷重Fweightと、把持部1000bに加えられる外力Fhandとの両方を推定することができる。従って、外力サンサ706を省略し、構成を簡略化することが可能になる。
また、上述の実施形態では、自移動体1000Rと他移動体1000Lとを相互に連携動作させるものとしたが、或る移動体に他の移動体を従属的に連携させるものとしてもよい。この例について更に説明すると、前述の図21において、移動する2台の倒立振子型移動体1000R,1000Lに対して一定の距離を維持しながら、他の倒立振子型移動体(図示なし)が追従するものとしてもよく、2台の倒立振子型移動体1000R,1000Lは、これらに追従する後続の倒立振子型移動体(図示なし)に対しては自身の動作を何ら連携させない。この例では、倒立振子型移動体1000R,1000Lに追従する倒立振子型移動体において、追従すべき対象の方向を特定する必要が生じるが、例えば、先行する倒立振子型移動体1000Rに特定周波数の信号を発信する装置を装着しておけば、その信号の到来方向から、後続の倒立振子型移動体が追従すべき対象の方向を把握することができ、従って倒立振子型移動体1000Rに対する追従動作が可能になる。