(57)【要約】【課題】筋電位のフィードバックにより動作した場合も、安定にアシスト力を生成することが可能なパワーアシストロボットを提供する。【解決手段】  外骨格型ロボット１を構成する各能動関節は、ユーザの筋電位を計測する筋電位センサと、能動関節を駆動するアクチュエータを備えている。制御部１３４は、外部制御装置２０からのコマンドにしたがって、筋電位のフィードバックに基づくアシスト力のフィードバック制御において、能動関節が発生すべき目標トルクを、筋電位から推定される関節駆動トルクに第１のフィードバック係数を乗じた値とアシスト力を発生するためのアシストトルクに第２のフィードバック係数を乗じた値の和として制御する。【選択図】図３
2014155998.tif 000002
【特許請求の範囲】
【請求項１】
ユーザの筋骨格系運動の支援をするためのパワーアシストロボットであって、
外骨格に対応するフレーム構造と、
前記筋骨格系運動において対象となる前記ユーザの関節にサポート力を与えるように配置される能動関節
と、
前記能動関節による前記サポート力を受ける前記ユーザの筋肉の筋電位を計測するための筋電位計測手段と、
前記能動関節に対するトルク制御により前記能動関節を動作させる制御部とを備え、
前記制御部は、前記筋電位のフィードバックに基づくアシスト力のフィードバック制御において、前記能動関節が発生すべき目標トルクを、前記筋電位から推定される前記ユーザの関節駆動トルクに第１のフィードバック係数を乗じた値と前記アシスト力を発生するためのアシストトルクに第２のフィードバック係数を乗じた値の和として制御する、パワーアシストロボット。
【請求項２】
前記制御部は、
前記計測された筋電位を非線形モデルにより筋張力に変換し、前記筋張力を線形モデルにより前記関節駆動トルクに変換する関節トルク変換手段と、
前記関節駆動トルクと、前記能動関節の発生するアシストトルクとに基づき、前記目標トルクを算出する目標トルク算出手段と、
前記目標トルクから前記能動関節の駆動を制御する制御信号を生成する制御信号生成手段とを含む、請求項１記載のパワーアシストロボット。
【請求項３】
前記目標トルク算出手段は、前記アシストトルクを、前記目標トルクに、前記目標トルクから前記アシスト力へ変換するヤコビアンと、ユーザの関節角変化に対する前記ユーザと前記パワーアシストロボットとのインタラクション部分の位置変化を表すヤコビアンの転置とを乗じた値として算出する、請求項２記載のパワーアシストロボット。
【請求項４】
各前記能動関節は、
エアマッスルと、
電動モータと、
前記エアマッスルからの第１の駆動力と前記電動モータからの第２の駆動力とを合成して前記能動関節の屈伸を駆動するための駆動力合成手段と、
前記エアマッスルからの前記第１の駆動力を前記駆動力合成手段に伝達するためのワイヤとを含み、
前記制御部は、
前記ワイヤを所定のバネ定数で表されるバネとした力学モデルと前記エアマッスルの収縮割合に対する前記第１の駆動力の所定の関係を表す駆動力モデルとの組合せの逆モデルにより、前記第１の駆動力を生成するために前記エアマッスルに加える圧力を制御する、請求項１〜３のいずれか１項に記載のパワーアシストロボット。
【請求項５】
各前記能動関節の関節角を検出するための第１のセンサをさらに備え、
前記制御部は、
前記第１のセンサの検出結果により、前記エアマッスルの前記収縮割合を検知する、請求項４記載のパワーアシストロボット。
【請求項６】
各前記エアマッスルにおいて、前記エアマッスルと前記ワイヤとの接合部に設けられ、前記エアマッスルからの前記第１の駆動力の大きさを検出するための第２のセンサをさらに備え、
前記制御部は、
前記第２のセンサの検出結果に応じて、前記エアマッスルの前記第１の駆動力および前記電動モータからの前記第２の駆動力とを制御する、請求項４または５に載のパワーアシストロボット。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、使用者（ユーザ）の行う運動を支援するパワーアシストロボットに関するものである。
【背景技術】
【０００２】
日本をはじめ多くの国で少子高齢化社会が問題になる中、ロボティクス技術を応用したアシスト機器への要望が高まっている。一方、バランスや歩行が可能なロボットが開発されてきている。例えば、運動に必要な作用力を空間上の任意の複数接触点に最適に配分し，ヒトと同じように各関節のトルクを発生できるロボットが存在する（特許文献１参照）。
【０００３】
また、近年では、下肢・体幹運動の支援をめざした外骨格型ロボットのようなリハビリテーションを支援するロボットの開発をますます要求することとなってきている（非特許文献１）。たとえば、外骨格型ロボットは、患者の自立生活を促進するリハビリテーションにおいて、脊髄損傷患者のために使用される（たとえば、非特許文献２、３、４、５、６を参照）。
【０００４】
これらの外骨格型ロボットのうち、体重を支えるタイプでは、大型の電動モータまたは油圧シリンダが用いられており、ロボット本体の自重も相当大きい。一方、体重を支えないリハビリ用途としては、空圧式の人工筋肉（以下、「エアマッスル」という。）を用いるものがあった。
【０００５】
エアマッスルは、その他のアクチュエータと比べて軽量で安価に製造可能で、空圧独特の柔らかさを持つ。その反面、遅れ、非線形性、個体のばらつきなどから、正確な制御には向いていないとされてきた。
【０００６】
このようなエアマッスルの問題に対処するために、空圧式のエアマッスルと電動モータとを組み合わせた空電ハイブリッド式アクチュエータも提案されている（特許文献２を参照）。
【０００７】
一方で、たとえば、脳卒中片麻痺患者などに対しての下肢のリハビリテーション支援を目的として、外骨格型ロボットを制御するには、筋電位（Electromyogram：ＥＭＧ）をもちいることが、有効な手段である。特に上肢を対象とした外骨格型ロボットでは運動意図の推定情報源として広く用いられている。
【０００８】
ここで、従来、筋電位（ＥＭＧ）はヒトの運動意図を推定するための情報源として利用され、その応用にあたっては、ＥＭＧより関節の駆動トルクτemgを推定することが行われてきた。すなわち、以下のようにして、ＥＭＧ信号の値にフィードバックゲインをかけるフィードバック方法が一般的である（非特許文献７）。
【０００９】
【数１】
2014155998.tif 000003
ここでα1は筋電位推定トルクのフィードバック係数である。（変数の上付きのアスタリスク（＊）はタスクの目標値を意味する。）
【先行技術文献】
【特許文献】
【００１０】
【特許文献１】ＷＯ２００７／１３９１３５号公報
【特許文献２】特開２０１２−０４５１９４号公報
【非特許文献】
【００１１】
【非特許文献１】K. Suzuki, M. G. Kawamoto, H. Hasegarwa, and Y. Sankai, “Intension-based walking support for paraplegia patients with Robot Suit HAL,” Advanced Robotics, vol. 21, no. 12, pp. 1441-1469, 2007.
【非特許文献２】K. Suzuki, M. G. Kawamoto, H. Hasegarwa, and Y. Sankai. Intension-based walking support for paraplegia patients with Robot Suit HAL. Advanced Robotics, Vol. 21, No. 12, pp. 1441-1469, 2007.
【非特許文献３】S. K. Au, P. Dilworth, and H. Herr. An ankle-foot emulation system for the study of human walking biomechanics. In IEEE International Conference on Robotics and Automation, pp. 2939{2945, 2006.
【非特許文献４】H. Kobayashi, A. Takamitsu, and T. Hashimoto. Muscle Suit Development and Factory Application. International Journal of Automation Technology, Vol. 3, No. 6, pp. 709-715, 2009.
【非特許文献５】G. Yamamoto and S. Toyama. Development of Wearable Agri-Robot- Mechanism for Agricultural Work. In IEEE/RSJ International Conference on Intelligent Robots and System, pp. 5801{5806, 2009.
【非特許文献６】T. Kagawa and Y. Uno. Gait pattern generation for a power-assist device of paraplegic gait. In The 18th IEEE International Symposium on Robot and Human Interactive Communication, pp. 633{638, 2009.
【非特許文献７】C. Fleischer, and G. Hommel, “A human-exoskeleton interface utilizing electromyography,” IEEE TRANSACTION ON ROBOTICS, vol. 24, no. 4, pp. 872-882, August 2008.
【発明の概要】
【発明が解決しようとする課題】
【００１２】
しかしながら、このような従来の手法の問題は、人の筋電位からロボットの駆動トルクの一方向のみを考えており、人-ロボット間インタラクションを考慮していないことである。
【００１３】
つまり、従来技術においては、ＥＭＧより推定される関節駆動力τemgが人が必要とするτｈ＊であると解釈されている背景がある。
【００１４】
そして、一般的に、ＥＭＧのフィードバックゲインを高くすると、ヒト-ロボットインタラクションを含む系全体が不安定になるため、実用的にはゲインの調整が不可欠である。しかしながら、このようなフィードバックゲインとして、どのようなゲインが適切であるかについては、その試行の都度、調整を行う必要があるという問題があった。
【００１５】
また、空電ハイブリッド式アクチュエータは、相補的なアクチュエータとして、２つの異なるエネルギー源で駆動される駆動機構を利用し、弾性を有する柔軟なトルク制御が可能である。１つの駆動機構は、空圧式エアマッスルを用いた空圧式アクチュエータであって、軽量化に向いている。しかしながら、上述のとおり、空圧式アクチュエータの圧縮性は、コントロールする際には、レスポンスが遅く、非線形力学を引き起こすために制御が難しい。
【００１６】
一方で、他の１つの駆動機構である電動式アクチュエータは、速く制御しやすい特性を有している。しかしながら、電動モータを用いた電動式アクチュエータは、大きなトルクを制御しようとすると、アクチュエータのヒートアップや、あるいは、重量の増加を招いてしまう。このため、空電ハイブリッド式アクチュエータにおいても、電動式アクチュエータを大型化することは軽量化に反する結果となる。
【００１７】
ここで、上述したようなリハビリテーションを目的とした外骨格型ロボットの場合は特に、軽量化は重要な課題であるので、駆動力源である空圧式エアマッスルや電動モータの軽量化は、もちろんで必要であるが、これらの駆動力源からの力を外骨格型ロボットの関節部に伝達するための伝達機構の軽量化も重要な課題である。
【００１８】
しかしながら、伝達機構の軽量化は、制御対象の力学をより複雑化して、制御の困難性を増大させてしまう恐れがある。
【００１９】
本発明は、上記のような問題点を解決するためになされたものであって、その目的は、筋電位のフィードバックにより動作した場合も、安定にアシスト力を生成することが可能なパワーアシストロボットを提供することである。
【００２０】
この発明の他の目的は、空電ハイブリッド式アクチュエータを備えるパワーアシストロボットを軽量化しつつ、正確な制御性を達成することが可能なパワーアシストロボットを提供することである。
【課題を解決するための手段】
【００２１】
この発明の１つの局面に従うと、ユーザの筋骨格系運動の支援をするためのパワーアシストロボットであって、外骨格に対応するフレーム構造と、筋骨格系運動において対象となるユーザの関節にサポート力を与えるように配置される能動関節と、能動関節によるサポート力を受けるユーザの筋肉の筋電位を計測するための筋電位計測手段と、能動関節に対するトルク制御により能動関節を動作させる制御部とを備え、制御部は、筋電位のフィードバックに基づくアシスト力のフィードバック制御において、能動関節が発生すべき目標トルクを、筋電位から推定されるユーザの関節駆動トルクに第１のフィードバック係数を乗じた値とアシスト力を発生するためのアシストトルクに第２のフィードバック係数を乗じた値の和として制御する。
【００２２】
好ましくは、制御部は、計測された筋電位を非線形モデルにより筋張力に変換し、筋張力を線形モデルにより関節駆動トルクに変換する関節トルク変換手段と、関節駆動トルクと、能動関節の発生するアシストトルクとに基づき、目標トルクを算出する目標トルク算出手段と、目標トルクから能動関節の駆動を制御する制御信号を生成する制御信号生成手段とを含む。
【００２３】
好ましくは、目標トルク算出手段は、アシストトルクを、目標トルクに、目標トルクからアシスト力へ変換するヤコビアンと、ユーザの関節角変化に対するユーザとパワーアシストロボットとのインタラクション部分の位置変化を表すヤコビアンの転置とを乗じた値として算出する。
【００２４】
好ましくは、各能動関節は、エアマッスルと、電動モータと、エアマッスルからの第１の駆動力と電動モータからの第２の駆動力とを合成して能動関節の屈伸を駆動するための駆動力合成手段と、エアマッスルからの第１の駆動力を駆動力合成手段に伝達するためのワイヤとを含み、制御部は、ワイヤを所定のバネ定数で表されるバネとした力学モデルとエアマッスルの収縮割合に対する第１の駆動力の所定の関係を表す駆動力モデルとの組合せの逆モデルにより、第１の駆動力を生成するためにエアマッスルに加える圧力を制御する。
【００２５】
好ましくは、各能動関節の関節角を検出するための第１のセンサをさらに備え、制御部は、第１のセンサの検出結果により、エアマッスルの収縮割合を検知する。
【００２６】
好ましくは、各エアマッスルにおいて、エアマッスルとワイヤとの接合部に設けられ、エアマッスルからの第１の駆動力の大きさを検出するための第２のセンサをさらに備え、制御部は、第２のセンサの検出結果に応じて、エアマッスルの第１の駆動力および電動モータからの第２の駆動力とを制御する。
【発明の効果】
【００２７】
本発明のパワーアシストロボットによれば、筋電位のフィードバックにより動作した場合も、安定にアシスト力を生成することが可能である。
【００２８】
また、この発明によれば、エアマッスルと電動モータとを組み合わせ、駆動力の伝達機構を軽量化して、人間の運動のアシストに必要な高負荷トルクを正確に制御できるパワーアシストロボットを提供できる。
【図面の簡単な説明】
【００２９】
【図１】本実施の形態における外骨格型ロボット１の構成事例を示す図である。
【図２】外骨格型ロボット１の自由度の構成を示す図である。
【図３】外骨格型ロボット１のブロック図の例である。
【図４】本実施の形態の外骨格型ロボットを使用者が装着した場合の構成を示す図である。
【図５】本実施の形態において、筋電位から関節駆動トルクを推定して、能動関節のアクチュエータ（空電ハイブリッド式アクチュエータ）を制御するための制御ブロックを示す図である。
【図６】ＥＭＧ信号フィルタ部２１０２による筋電位デジタルフィルタの処理を示すフローチャートである。
【図７】シミュレーションでの関節角度、人の筋活性、１軸アシスト装置によるアシスト力、およびロボットと人のトルクを比較して示す図である。
【図８】外骨格型ロボットの構成において、１自由度分の動作を行う空電ハイブリッドアクチュエータの部分を抜き出した外観を示す図である。
【図９】図８に示した１自由度分の動作を行う空電ハイブリッドアクチュエータのシステムを制御する構成を説明するための機能ブロック図である。
【図１０】各アクチュエーターが、関節を動かす態様を説明するための概念図である。
【図１１】１軸アシスト装置の実験セットアップの外観を示す図である。
【図１２】筋電位計測のための体表電極の位置の例を示す図である。
【図１３】筋電位計測により推定された被験者の関節駆動トルクτhと１軸アシスト装置の駆動トルクτrを示す図である。
【発明を実施するための形態】
【００３０】
以下、本発明の実施の形態の外骨格型ロボットの構成について、図に従って説明する。なお、以下の実施の形態において、同じ符号を付した構成要素および処理工程は、同一または相当するものであり、必要でない場合は、その説明は繰り返さない。
【００３１】
また、外骨格型ロボットの関節を駆動するためのアクチュエータとしては、一例として、以下に説明する「空電ハイブリッド式のアクチュエータ」を例として説明する。
【００３２】
そこで、以下、本実施の形態において、歩行・姿勢リハビリテーションのための空電ハイブリッド式のアクチュエータによる外骨格型ロボットについて説明する。
【００３３】
ただし、本発明の空電ハイブリッド式外骨格型ロボットは、下肢の運動をアシストするための外骨格型ロボットに対してだけでなく、上肢の運動をアシストする外骨格型ロボットとしても使用することが可能である。
【００３４】
また、以下の説明では、下肢の対としての運動をアシストする外骨格型ロボットについて説明するが、下肢のうちのいずれか一方、または、上肢のうちのいずれか一方の運動をアシストする外骨格型ロボットとして使用することも可能である。
【００３５】
さらに、本発明の空電ハイブリッド式外骨格型ロボットは、対象となる人間の筋骨格系の運動をアシストするのであれば、上述したような「下肢のうちの少なくともいずれか一方、または、上肢のうちの少なくともいずれか一方の運動」に限定されるものではなく、たとえば、対象となる人間の腰の運動のみをアシストするものであってもよいし、歩行または走行時において下肢の運動と連動して腰の運動をアシストするものであってもよい。本明細書では、このような対象となる人間の運動のアシストを総称して、「対象となる人間の筋骨格系運動の支援」と呼ぶことにする。
【００３６】
本実施の形態の外骨格型ロボットは、外骨格を有する。「外骨格」とは、人間の骨格構造に対応してロボットが有する骨格構造のことである。より特定的には、「外骨格」とは、外骨格型ロボットを装着する人間の体の一部を、外部から支えるフレーム（枠組み）構造のことをいう。
【００３７】
このフレーム構造には、さらに、フレーム構造の各部を人間の骨格構造に基づく運動に応じて動かすための関節が設けられる。
【００３８】
特に、下肢の運動をアシストする外骨格型ロボットは、ベースと下半身とを有し、足首、膝、腰の左右の位置に、能動６自由度の関節を有するロボットである。また、当該６つの関節は、空電ハイブリッド駆動の関節である。以下、このように、外骨格型ロボットにおいて、ユーザの関節に対してサポート力を与えるためにアクチュエータにより駆動される関節のことを「能動関節」と呼ぶ。
【００３９】
図１は、本実施の形態における外骨格型ロボット１の構成事例を示す図である。本外骨格型ロボット１は、１０自由度である。
【００４０】
図１において、図１（ａ）は、外骨格型ロボットの外観を示す図であり、図１（ｂ）は、図１（ａ）の外観において、外骨格型ロボット１の主要部を抽出して示す斜視図である。
【００４１】
図１（ｂ）において、外骨格型ロボット１は、両脚に対応したフレーム構造、バックパック１０１、柔軟シート１０２、ＨＡＡ拮抗筋１０３、ＨＦＥ伸筋１０４、ＨＦＥモータ１１１、ＫＦＥ伸筋１０５、 ＫＦＥモータ１０６、ＡＦＥ伸筋・ＡＡＡ拮抗筋１０７、ＡＦＥ屈筋１０８、ユニバーサルジョイント１０９、フレーム構造に設けられたプーリー付回転関節１１０を備える。
【００４２】
なお、図１（ｂ）では、バックパック１０１が運動を支援する構造に直接とりつけられているが、図１（ａ）に示すように、バックパック１０１は、この構造から取り外されていてもよい。
【００４３】
また、ユニバーサルジョイント１０９には、たとえば、光学式エンコーダを回転軸に取り付け、関節角度を計測する。プーリー付回転関節１１０も同様に光学式エンコーダを取り付ける。光学式エンコーダは、軸に取り付けるのではなく、軸に巻かれたベルトの移動方向と移動量を読み取る構成としてもよい。なお、ハイブリッド関節であるＨＦＥおよびＫＦＥ関節においては、モータ付属のエンコーダを用いて関節角度を計測してもよい。
【００４４】
また、図２は、外骨格型ロボット１の自由度の構成を示す図である。
【００４５】
図２において、各関節において、「Ｒ＿」との表示は、右側の関節であることを示し、「Ｌ＿」との表示は、左側の関節であることを示す。
【００４６】
図１および図２を参照して、全１０自由度のうち、ＨＦＥ関節とＫＦＥ関節はハイブリッド駆動としている。また、図２において、全１０自由度のうち、左右のＡＦＥ関節は伸筋と屈筋による拮抗駆動を採用している。ハイブリッド駆動および拮抗駆動以外の関節は、パッシブな駆動である。ただし、より多くの関節、たとえば、全ての関節をハイブリッド駆動としてもよい。
【００４７】
図１において、両脚が接続する胴体部には姿勢センサを搭載してベース部の姿勢を検出している。また、全ての関節にワイヤ式エンコーダ（またはモータ付属のエンコーダ）を取り付け、関節角度を計測できるようにしている。関節角度ならびに後に説明するように下肢の筋電位を検出することで、各関節に発生させる目標トルクが算出できる。
【００４８】
また、足底部には、床反力センサを搭載し、接触を想定する足底部が実際に接触しているかどうかを判定したり、ヤコビ行列に含まれるモデル誤差を修正するために補助的に使用する構成としてもよい。
【００４９】
また、バックパック１０１内には制御器の他、エアマッスルのバルブおよび電動モータのドライバを内蔵している。
【００５０】
また、バックパック１０１内に、バッテリーと圧搾した空気ボンベ（またはＣＯ２ガスボンベでもよい）、レギュレータを搭載し、電源ラインとエア供給が断絶した場合に備え、短時間の自律駆動を可能にする構成であってもよい。
【００５１】
また、図３は、外骨格型ロボット１のブロック図の例である。
【００５２】
外骨格型ロボット１を制御するためのコマンドが、外部制御装置２０から、通信経路を介して外骨格型ロボットに与えられる。特に限定されないが、外部制御装置２０は、汎用のパーソナルコンピュータを用いることが可能であり、通信経路としては、イーサネット（登録商標）ケーブルを用いることができる。もちろん、通信経路としては、その他の規格の有線通信の経路の他、無線による通信経路、たとえば、無線ＬＡＮ（Local Area Network）や他の通信規格の無線などを使用してもよい。
【００５３】
外部制御装置２０は、ユーザからの指示入力を受ける入力部２０８と、コマンドを生成するためのプログラムや、様々な制御パラメータなど制御のために必要とされるデータが記録された不揮発性の記憶装置２０６と、外部制御装置２０を起動するためのファームウェアが記憶されたＲＯＭ（Read Only Memory）や、ワーキングメモリとして動作するＲＡＭ（Random Access Memory）などを含むメモリ２０４と、プログラムに応じて、コマンドを生成する処理を実行する演算装置２１０と、コマンドを通信経路を介して、外骨格型ロボットに送信するためのインタフェース（Ｉ／Ｆ）部２０２と、演算装置２１０の制御の下で、外骨格型ロボット１への制御の状態に関する情報などを表示するための表示装置２１２とを備える。
【００５４】
上述のとおり、外部制御装置２０が、汎用のパーソナルコンピュータである場合は、演算装置２１０は、ＣＰＵ（Central Processing Unit）で構成され、不揮発性の記憶装置２０６としては、ハードディスクドライブやソリッドステートドライブなどを用いることができる。ただし、外部制御装置２０の機能ブロックの一部または全部は、専用のハードウェアにより構成されてもよい。
【００５５】
また、外部制御装置２０は、たとえば、外骨格型ロボットが装着される使用者について、筋電位から関節トルクを推定するためのモデルをキャリブレーション時に構成する処理を行う。
【００５６】
外骨格型ロボット１は、外骨格部１２、内部制御装置１０を備える。
【００５７】
外骨格部１２は、ベース１２１、下半身１２２、能動関節１２３、検出機構１２４を備える。さらに、能動関節１２３は、エアマッスル１２３１（図示せず）、電動モータ１２３２（図示せず）を備える。
【００５８】
また、内部制御装置１０は、Ｉ／Ｆ部１１、記録装置１３１、記憶装置 １３２、計測装置１３３、制御部１３４、出力装置１３５を備える。
【００５９】
Ｉ／Ｆ部１１は、外部制御蔵置２０から指令されたトルクまたは位置指令等を受け付けることができる。
【００６０】
なお、ベース１２１は、腰の位置の骨格、腰の位置の能動関節１２３を含むと考えても良いし、腰の位置の骨格のみであると考えても良い。
【００６１】
下半身１２２は、腿や足の位置の骨格、腿や足の位置の能動関節１２３を含むと考えても良いし、腿や足の位置の骨格のみであると考えても良い。
【００６２】
能動関節１２３は、左右の足首、左右の膝、および腰の左右の各位置に配置されている能動の関節である。ここで、能動関節１２３とは、アクチュエータで能動的に動作することのできる関節である。つまり、能動関節１２３は、アクチュエータを備える。
【００６３】
また、ここでの１以上の能動関節１２３は、ハイブリッド型である。つまり、能動関節１２３の少なくとも一部のものは、エアマッスル１２３１、電動モータ１２３２を備えるハイブリッド型である。なお、アクチュエータは、制御目標値となるトルク値を駆動信号として受け付け、受け付けたトルク値に基づいて制御する機能を有している。
【００６４】
アクチュエータとして、サーボモータを使用する場合、アクチュエータは、例えば、電流制御が可能な駆動回路を有し、電流に比例したトルクを発生させるサーボモータは、制御目標値として入力されたトルク値に、ギヤ比により決定されるトルク定数を乗じて駆動回路に指令することで入力されたトルクを発生させるトルク制御を実現する。特に、能動関節１２３にトルクセンサを配設し、当該トルクセンサにより検出した値を駆動回路にフィードバックすることにより、高精度のトルク制御が可能となる。
【００６５】
検出機構１２４は、ロボットの状態を検出する。検出機構１２４は、例えば、各関節に配置されたエンコーダ、足平に配置された床反力センサ、骨盤部に配置された姿勢検出のためのジャイロセンサ、各エアマッスルの駆動力を検知するロードセルなどである。検出機構１２４は、関節の角度を検出する角度センサや、ロボットの姿勢を取得する姿勢センサ、外力センサなどでも良い。
【００６６】
また、検出機構１２４は、使用者の脚部に装着された筋電センサにより使用者の筋電位を検出する。
【００６７】
内部制御装置１０は、能動関節１２３を動作させる。内部制御装置１０は、Ｉ／Ｆ部１１が受け付けた目標トルクまたは位置指令等に対応して、能動関節１２３を動作させる。
【００６８】
計測装置１３３は、センサ等の検出機構１２４から検出結果を示す様々な信号（データ）を受け付ける。
【００６９】
制御部１３４は、制御目標値の算出等の様々な演算を行う。制御部１３４が行う演算は後に説明する。
【００７０】
出力装置１３５は、能動関節１２３に制御信号を出力する。出力装置１３５は、例えば、目標とするエアマッスルの圧力値やモータの制御値を能動関節１２３に出力する。より詳しくは、このような制御値は、モータドライバや圧力制御バルブなどに出力され、物理的なエネルギー（圧縮空気や電流）がモータやエアマッスルに供給される。
【００７１】
図４は、本実施の形態の外骨格型ロボットを使用者が装着した場合の構成を示す図である。
【００７２】
外骨格型ロボットは、自律して立位姿勢を保持することが可能である。
【００７３】
（筋電位フィードバックによる能動関節の制御）
上述のとおり、従来のように、式（１）に従うフィードバックを行った場合、ＥＭＧ信号により推定される関節駆動力τemgが人が必要とするτｈ＊であると解釈されている背景がある。そして、一般的に、ＥＭＧのフィードバックゲインを高くすると、ヒト-ロボットインタラクションを含む系全体が不安定になるという問題がある。
【００７４】
これは、例えば、ある時点でシステムが完全に筋力をサポートすると、次の瞬間には必要力が０になるため、式（１）よりサポート力は０となる。その結果、ロボットからのアシスト力は０である。つまり、人の必要力は増大する。これは人の筋電位調整とシステムの応答が理想的なケースであるが、実際は拘束部分の緩みやシステムの遅れなどからα1を増加すると振動的な挙動となる。そのため、α1は小さく設定しなければならなかった、という事情による。
【００７５】
本実施の形態では、人-ロボット間インタラクションによる発生力を考慮することで、筋電位から予測するフィードバックトルクを調整することをコンセプトとする。なぜならば、ヒトはアシストされた分だけ筋の発生力を減らすよう調整するからである。
【００７６】
【数２】
2014155998.tif 000004
ここで、τｈ＊はタスクに必要なトルク、τemg＊はヒトが発生すべき関節トルクの目標値であり、τrhはロボットとのインタラクションによるアシストトルクである。
【００７７】
この式（２）より、ロボットの発生力の目標値を以下のように調整する。
【００７８】
【数３】
2014155998.tif 000005
ここで、式（３）の第１項は式（１）と共通した項であり、第２項はインタラクションの調整項である。つまり、人-ロボット間のインタラクションの結果アシストされる駆動トルクτrhの分だけ必要トルクが増分されていることがわかる。この増分のフィードバック係数がα2である。ただし、τrhは計測や評価が、一般には、非常に難しいため、式（３）は概念上の式である。そこで、第２項については、実測可能な量への置き換えを行う。このような第２項のτrhの本実施の形態での実装に関しては後述することとする。
【００７９】
図５は、本実施の形態において、筋電位から関節駆動トルクを推定して、能動関節のアクチュエータ（ここでは、一例として、空電ハイブリッド式アクチュエータ）を制御するための制御ブロックを示す図である。
【００８０】
図５を参照して、関節トルク変換部２１００は、アシスト対象である使用者の下肢２１４０から測定された筋電位（ＥＭＧ）信号は、以下に説明するような手続きにより、非線形なモデルにより筋張力を算出し、線形モデルにより、関節トルクを計算する。
【００８１】
より具体的には、関節トルク変換部２１００は、ＥＭＧ信号に対して所定のフィルタリング処理を行うＥＭＧ信号フィルタ部２１０２と、フィルタリングにより得られる筋ｉから出力ｕi（ｔ）から、筋の活性度を算出して、筋の活性度ｑiから非線形モデルにより筋張力ζiを算出する筋張力算出部２１０４と、筋張力ζiから関節駆動トルクτemgを推定する関節トルク算出部２１０６とを含む。
【００８２】
さらに、目標トルク算出部２１１０は、推定されたτemgからアクチュエータの生成するべき目標トルクτｒ＊を算出する。
【００８３】
目標トルク算出部２１１０は、逆動力学作用モデル算出部２１１２を含んでおり、逆動力学作用モデル算出部２１１２は、後に説明するように、目標トルクτｒ＊からのフィードバック２１１４により、目標トルクτｒ＊を算出する。
【００８４】
腱スプリングモデル算出部２１２０は、エアマッスルに供給する圧力ｐ＊およびモータへの制御値を算出する。制御対象２１３０は、後に説明するような能動関節に設けられた空電ハイブリッド式アクチュエータである。
［筋電位からの関節駆動力τemgの推定］
次に、以下では、筋電位から関節駆動力τemgを推定することを考える。
（筋電位に基づく関節駆動トルクの推定）
特に限定されないが、たとえば、使用者の脚部には、臀部の関節角を検出するための角度計と、膝の関節角を検出するための角度計と、くるぶし部分の角度を検出するための角度系とが装着されており、また、使用者の脚部には、後に説明するように、使用者の脚部の筋電位を計測するための筋電センサが装着されているものとする。
【００８５】
さて、筋電位に基づくアシストロボットの制御は、長い間、取り組まれている課題である。ある筋の活性度をq、筋張力をζとし、qからζおよびτemgの変換について考える。
【００８６】
この変換には非線形があることが知られており、以下の文献に加持されたHill-stroeveモデルが広く普及している。
【００８７】
公知文献１：A. Hill, “The heat of shortening and the dynamic constants of muscle,” Proceedings of the Royal Society of London. Series B, Biological Sciences, vol. 126, no. 843, pp. 136-195, 1938.
公知文献２：S. Stroeve, “Learning combined feedback and feedforward control of a musculoskeletal system,” Biological Cybernetics, vol. 75, no. 1, pp. 73-83, 1996.
この非線形変換を機械学習（ニューラルネットワークなど）で学習するアプローチでは、筋電位から直接関節トルクτemgへの変換を実験データより推定し、フィードバック制御が試みられている。
【００８８】
しかし、これらの手法では筋電位を貼り直すたびに再学習が必要であり、キャリブレーションに多く労力が必要となる。一方で、安定性や実用性の観点からは線形モデルによる筋電位フォードバック制御も用いられているが、さまざまな動き（筋の収縮率・速度）が含まれる場合は推定誤差が増大する。
【００８９】
Fleischerらが、以下の文献で提案した手法は筋張力ζから関節駆動トルクτemgの変換は線形変換で、活性度qとζは非線形変換を行っている。非線形変換は実験ごとにパラメータを共通化できるため、キャリブレーションに必要なデータ量を少なくできる。
【００９０】
公知文献３：C. Fleischer, C. Reinicke, and G. Hommel, “Predicting the intended motion with emg signals for an exoskeleton orthosis controller,” Intelligent Robots and Systems, 2005.(IROS 2005). 2005 IEEE/RSJ International Conference on, pp. 2029-2034, 2005.
［ＥＭＧ信号から膝駆動トルクを算出するモデルの構築（キャリブレーション）］
本実施の形態では、生体力学的な知見に基づき筋活性度から筋張力の非線形変換を考え、線形モデルを仮定したキャリブレーションにより、筋の太さや種類およびインピーダンスや電極位置の変化といった実験設定によって変化する対象を考慮する。本実施の形態では、筋張力算出部２１０４において実行される、筋ｉ
の活動度qiから筋張力ζiへの変換は、以下の近似関数を用いる。ここで、「筋ｉ」とは、筋電位測定のためにユーザの体表につけられた電極のチャネルを示す。そして、電極のチャネル番号と、その電極がいずれの筋肉につけられているかについての対応関係については、記憶装置内に予め格納されているものとする。
【００９１】
【数４】
2014155998.tif 000006
ここで、記号の意味は以下の通りである。
【００９２】
【数５】
2014155998.tif 000007
非線形変換では、関節トルク推定処理
の線形化を目的とするため、解剖学的な筋の最大力は１とした。ｋ（…）とｈ（…，…）は、以下のように、筋繊維に依存する２つのパラメータ（ａ1とａ2）を持つ関数である。
【００９３】
【数６】
2014155998.tif 000008
ここで、パラメータａ1とａ2は、筋繊維に依存するパラメータであり、実験的に、たとえば、ａ1＝２．７およびａ2＝０．１を用いた。
【００９４】
また、それぞれの変換写像ｋ（），ｈ（）は、たとえば、以下の文献に開示されている。
【００９５】
公知文献４：H. Hatze：“A myocybernetic control model of skeletal muscle” Biological Cybernetics,25,103/119,1977
（筋電位フィルタリング）
図６は、ＥＭＧ信号フィルタ部２１０２による筋電位デジタルフィルタの処理を示すフローチャートである。
【００９６】
筋電位を測定する電極からアンプを通した筋電位は所定の時間間隔でンプリングされる（Ｓ１００）。この時間サンプリングをｓtとする。
【００９７】
続いて、時間サンプリングをｓtは、ハイパスフィルタにより直流成分が除去されるとともに、正規化される（Ｓ１０２）。
【００９８】
続いて、ＥＭＧ信号フィルタ部２１０２は、信号の整流（つまり絶対値とる処理）を実行する（Ｓ１０４）。
【００９９】
ここで、このようなステップＳ１０２とステップＳ１０４の処理は、以下のようにして実行することが可能である。
【０１００】
すなわち、以下のような処理を逐次的に行うことで、直流成分の除去と信号の絶対値成分をとる処理を行うことができる。
【０１０１】
【数７】
2014155998.tif 000009
ここで、γは、スムージング係数であり、ｓt´は、ハイパスフィルタ処理される。また、サンプリング周期は２５０[Hz]である。
【０１０２】
さらに、信号は、２次のFIR バターワース（butter-worth）フィルタをかけられ（Ｓ１０６）、さらにFIR butter-worthフィルタ出力と矩形窓をかけた以下のような積分により、信号ｕi（ｔ）が得られる（Ｓ１０８）。
【０１０３】
【数８】
2014155998.tif 000010
このような筋電位フィルタを通した出力ｕi（ｔ）により、筋の活性度qiは次式より定義される（Ｓ１１０）。
【０１０４】
【数９】
2014155998.tif 000011
（関節角駆動トルク）
（線形トルク推定モデル）
本実施の形態では、関節トルク算出部２１０６の実行するＥＭＧ信号からの関節駆動トルク推定は、以下のようにして行う。
【０１０５】
すなわち、推定筋張力ζから関節角駆動トルクτemgは、関節の仮想的なプーリー径を仮定し、線形変換ＢとＣを考え、次のＥを最小化するＷを最小二乗法で求める。
【０１０６】
【数１０】
2014155998.tif 000012
τi（ｉ＝１，…）はキャリブレーションの実験データにおいてヒトの筋骨格モデルより算出される関節駆動トルクである。
【０１０７】
以上のようなモデルについて、オフラインでキャリブレーション時に、パラメータＢ，Ｃを設定する。
（インタラクションを考慮したフィードフォワードのアシストトルク調整）
式（３）では筋電位から推定された関節角駆動トルクに加え、インタラクショントルクにより発生するアシストトルクでフォードバック量を調整している。本実施の形態では人およびロボットの運動学的なモデルよりフィードフォワード的に扱う。
【０１０８】
【数１１】
2014155998.tif 000013
ここで、Ｊrはロボットのトルクからエンドエフェクタ（一般には、ロボットから見た使用者への相互作用部分であり、この場合、外骨格が使用者に力を作用させる部分）のインタラクション力（アシスト力）へ変換するヤコビアンである。Ｊhはインタラクション力から関節駆動トルクへの影響を表す静力学的なモデルであり、ヒトの関節角変化に対するインタラクション部分の位置変化を表すヤコビアンである。τr＊はロボットの目標出力トルクである。
【０１０９】
ここでの「インタラクション部分」とは、エンドエフェクタのことであり、具体的には装具・バンドが拘束している点で、ヒト筋骨格とロボットの接点のことである。ただし、厳密には、ヒト筋骨格とロボットとは、面（分布力や圧力）で接触することになり，力だけでなくモーメントも伝達されるものの、ここでは簡単化して点で接触して集中力を発生していると近似的に考えている。
【０１１０】
もちろん、インタラクション力を触覚センサなどのセンサで直接測定できれば、それを使用することも可能であるが、式（８）のようにおくことで、直接測定の困難さを回避することができる。
［１軸アシスト装置についての実験］
以下では、本実施の形態の制御方式の有用性を確認するために、１軸アシスト装置による筋電位に基づくトルクフォードバックのシミュレーションと実機により従来手法（式（１））と本実施の形態の手法（式（８））を比較した結果について説明する。
（１軸アシスト装置のシミュレーション）
図７は、シミュレーションでの関節角度、人の筋活性、１軸アシスト装置によるアシスト力、およびロボットと人のトルクを比較して示す図である。
【０１１１】
図７（ａ）は式（１）を筋電位フィードバックとして実装し、図７（ｂ）は式（８）を実装した場合のシミュレーション結果である。
【０１１２】
パラメータはそれぞれ（α1，α2）=（０．４５，０）および（α1，α2）=（０．４５，０．４５）である。シミュレータは、実機と同様にエアマッスルの駆動力をプーリーによってトルクに変換するアシスト装置である。人の筋の活性度のコントローラおよびエアマッスルのダイナミクスは以下のように１次遅れのある系として定義する。
【０１１３】
【数１２】
2014155998.tif 000014
ここで、Ｃfはエアマッスルの発生力の時定数であり、Ｃsは人の筋活性度の時定数である。Ｃfは次節で述べるエアマッスルの時定数に合わせて０．０８[sec]とし、Ｃsは運動野からαニューロンへの信号伝達時間に合わせて０．０２[sec]とした。
【０１１４】
シミュレーションは水平位置より４５[deg]の角度を初期位置としてヒトとロボットのトルクの合計値の目標値τtotal＊（＝τh＊）を正弦波で与えた。
【０１１５】
式（１）と式（８）の違いを強調するために、ｕt＊とτtotal＊ を３．７５[sec] 以降０とし、ｕt をトルクコントローラ入力から切断している。
【０１１６】
シミュレーションの結果（図７（ａ）と図７（ｂ））を比較すると、筋電位の入力が切れた後（３．７５ｓｅｃ以後）のエアマッスルＰＡＭ（Pnuematic Airtificial Muscle）の発生力ｆと、その結果発生するアシスト力τr（上から３番目および４番目のグラフ）の時間変化に、従来手法と本実施の形態の手法の違いが強調されている。図７（ａ）のfは短時間で発生力が０になっており、この変化はＰＡＭの空圧の時間遅れによるものであり、時定数も一致している。図７（ｂ）ではfが発生力が０になるまでに図７（ａ）よりも大きな時定数で変化している。
【０１１７】
これは、本実施の形態でのアシストトルクの調整項がｆ＊にダイナミクスを与える結果となったからである。図７（ａ）ではアシストトルクは約１/２であり、これはパラメータ（α1，α2）=（０．４５，０）および（α1，α2）=（０．４５，０．４５）の違い、すなわち図７（ｂ）においてはインタラクション力の調整項が存在することにより、筋電位に基づくフィードバック（アシストトルク）が増加することがわかる。
【０１１８】
（空電ハイブリッド式アクチュエータ用のトルク制御器）
以下では、空電ハイブリッド式アクチュエータ用のトルク制御器による力制御可能な１軸アシスト装置について説明する。エアマッスルＰＡＭは非常に大きな力を発生するため、高強力・高弾性の化繊ワイヤであっても伸びがモデルの誤差の原因となる。「腱スプリングモデル」を仮定した力制御と小型のモーターが発生するトルクにより、緻密な関節トルク制御が可能である。
【０１１９】
以下に説明するように、このようなトルク制御器（腱スプリングモデル算出部２１２０の実行する機能）は、制御部１３４に対して、機械学習を実行させることにより、実行させる構成とすることが可能である。
【０１２０】
空圧式エアマッスルは軽量であるもの、圧縮空気（または圧縮気体）のエネルギーをゴムチューブにより収縮力に変換することにより大きな力を生み出すことができる。
【０１２１】
エアマッスルが、力を発生させる原理は、空圧式の空気袋が埋め込まれたらせん状のファイバーが、圧縮空気が送り込まれて空気袋が膨張すると、長手方向（縦方向）に収縮するというものである。
【０１２２】
より詳しく説明すると、両端部に栓をしたゴムチューブに対して、動径方向を拘束するよう、表面にらせん状にまかれたファイバーを被せた構造になっている。このゴムチューブの中に空気を送り込むと、空気の圧力でゴムチューブが膨張する。しかし、動径方向はファイバーによって拘束されているため膨張できず、径方向の膨張に引っ張られて縦方向に収縮する。膨らみながら収縮する様が動物の筋肉に似ているところが人工筋肉と呼ばれる。
【０１２３】
アクチュエーター自体が軽く、柔らかい。さらに、ゴムチューブの内面全体がアクチュエーターの収縮に寄与するため、断面積のみで圧力を受ける構造になっている一般的なエアシリンダー等よりも、パワー・ウェイト・レシオを大きく取りやすい。反面、上述したように、一般に空気圧による機器制御は空気の収縮・膨張などによる、制御遅れが大きく、素早い動作は苦手である。
【０１２４】
なお、「空気袋」は、流体により膨張ないし収縮運動をするものであれば、袋中に流入するものは、空気に限られないので、より一般には、「流体袋」と表現する。
【０１２５】
図８は、外骨格型ロボットの構成において、１自由度分の動作を行う空電ハイブリッドアクチュエータの部分を抜き出した外観を示す図である。
【０１２６】
図８において、空圧式エアマッスルのうち上側のエアマッスル３０２は、屈筋(ＰＡＭ１)に相当し、下側のエアマッスル３０４ものは伸筋(ＰＡＭ２)に相当する。
【０１２７】
伸筋ＰＡＭ２は、屈筋ＰＡＭ１とは相反し反対向きの力を生成する。
【０１２８】
プーリー付回転関節３１０には、電動モータからの駆動力も印加され、空圧式エアマッスルからの駆動力と電動モータからの駆動力とが合成される。
【０１２９】
エアマッスル３０２からの駆動力は、ワイヤ３０６によりプーリー付回転関節３１０に伝達され、エアマッスル３０４からの駆動力は、ワイヤ３０８によりプーリー付回転関節３１０に伝達される。このような構成により、駆動力の伝達機構自体を軽量化することが可能となる。
【０１３０】
プーリー付回転関節３１０へ加えられた駆動トルクにより、脚（または腕）３５０が駆動される。
【０１３１】
図９は、図８に示した１自由度分の動作を行う空電ハイブリッドアクチュエータのシステムを制御する構成を説明するための機能ブロック図である。
【０１３２】
図９においては、内部制御装置１０は、マルチファンクションボード(Multi Function board)として構成される。
【０１３３】
外部制御装置２０に接続されたマルチファンクションボード１０は、外部制御装置２０からのコマンドに応じて、アクチュエーターを制御する。具体的には、マルチファンクションボード１０は、空圧式エアマッスル３０２および３０４の収縮を制御するためのバルブ３０１とバルブ３０２、ならびに、電動モータ３１２を制御するためのモータドライバ３１１を制御する。さらに、マルチファンクションボード１０は、関節角度θを検知する角度エンコーダ３２４およびエアマッスルからの駆動力を検知するロードセル、能動関節に加わるトルクを検知するトルクセンサからの計測データを読み取ることで、これらに基づき、以下に説明するような制御を実行する。
【０１３４】
角度エンコーダ３２４は、たとえば、直交エンコーダを使用することが可能である。
【０１３５】
エアマッスル３０２および３０４ならびに電動モータ３１２からの駆動力は、プーリー付回転関節３１０において合成され、脚３５０にトルクτを与える。
【０１３６】
図１０は、各アクチュエーターが、関節を動かす態様を説明するための概念図である。
【０１３７】
以下では、まず、図１０（ａ）に示されるように、腱ワイヤの伸長を無視できる場合を考える。
【０１３８】
空圧式エアマッスルフォースからの駆動力fPAMi（ｉ＝１，２，…）は、ワイヤとプーリにより、腕（または脚）にトルクとして転送され/変換される。このプーリーを介して、電動モータの駆動力も伝達される。空圧式エアマッスルフォースからの駆動力fPAMiは、空圧式エアーマッスルと駆動力をプーリーに伝達するワイヤーとの接合部に設けられたロードセルにより検知される。
【０１３９】
以下のモデルの説明は、拮抗筋と電動モータの組合せの場合に限らず、たとえば、伸筋と電動モータの組合せにも広く適用することが可能である。ただし、以下では、説明の簡単のために、単純な２つの拮抗筋のモデルを考える。
【０１４０】
トルクτPAMは、駆動力fPAMiにより、以下の式ように表される。
【０１４１】
【数１３】
2014155998.tif 000015
τPAMは、空圧式エアマッスルの駆動力によるトルクであり、一方、r0はプーリー半径で、この単純なモデルでは、r0は定数とする。
【０１４２】
モーターの駆動力によるトルクは、空圧式エアマッスルの駆動力と並行して伝達される。例えば、小さなトルクが機械的なベルトによりプーリー経由で伝達される。
【０１４３】
空圧式エアマッスルは、直流的または低周波のトルクの生成において優れている。また、電動モータによる駆動力は、迅速で周波数の高いトルクとして、空圧式エアマッスルの駆動力によるトルクτPAMの誤差をカバーするので、電動式のトルクτmotorの値は、空圧式エアマッスルによるトルクよりも小さな値でよい。
【０１４４】
トータルのトルクは、以下の式で表される。
【０１４５】
τ＝τPAM＋τmotor
したがって、特に限定されないが、たとえば、内部制御装置１０が、第一の閾値よりも（「より」は「以上」も含む、とする）高周波であるトルクに対して電動モータ３１２を追加動作させることは好適である。 なお、追加動作とは、エアマッスル３０２および３０４に加えて電動モータ３１２を動作させることである。ここで、第一の閾値は、例えば、３Hzである。
【０１４６】
空圧式エアマッスルから、ワイヤがプーリーへ駆動力を伝達し、指定された一方向へ関節を駆動する。
【０１４７】
ワイヤとしては、金属線、チェーンあるいは機械的なベルトなどの方式と比較して、軽量で強く、柔軟であるので、液晶ポリマー繊維を採用することができる。液晶ポリマー繊維としては、たとえば、クラレ社製ベクトラン（登録商標）を用いることができ、これは、高分子繊維であるにも関わらず、高強力高弾性率を有する素材として知られている。
【０１４８】
（空圧式エアマッスルの力学モデル）
空圧式エアマッスル、ワイヤおよびプーリを含む空圧式アクチュエーターは、人間の筋肉と多くの共通点を持っている。
【０１４９】
外骨格型ロボットが人間の重さを支持する際に、ワイヤは、その特性において、人間の腱に類似している。
【０１５０】
そこで、以下に説明するように、平衡の仮定を置いて、空圧式エアマッスルのためのトルク制御器を導出する。
【０１５１】
空圧式エアマッスルの圧力の制御には、たとえば、比例的に圧力を調整するバルブを使用し、圧力ｐはクローズドフィードバックループにより制御され、十分に安定であるものとする。
【０１５２】
過渡状態においては、バルブ圧力と空圧式エアマッスルの圧力との間には、空気力学的な運動が存在するものの、一定の時定数の後には、ロボットと人間の重量のような外部の運動上の制約条件に駆動力が釣り合うまで、空圧式エアマッスルは収縮する。
【０１５３】
したがって、空気回路の動力学の影響は小さく、準静的な動作では無視することができる。
【０１５４】
この均衡点で、駆動力の生成は、内圧および収縮割合に依存し、空圧式エアマッスルの駆動力モデルは、以下のように与えられる：
【０１５５】
【数１４】
2014155998.tif 000016
ここで、εは、収縮割合であり、Ｄ０とψ０とは、常圧における空圧式エアマッスルの径と、空圧式エアマッスルにおいて、空圧式の空気袋が埋め込まれたらせん状のファイバーの収縮方向に直交する方向に対する巻方向の傾きの角度である。
【０１５６】
このような「空圧式エアマッスルの駆動力モデル」は、たとえば、以下の文献に開示されている。
【０１５７】
公知文献５：K. Inoue. Rubbertuators and applications for robots. In Proceedings of the 4th international symposium on Robotics Research, pp. 57-63. MIT Press, 1988.
公知文献６：D.G. Caldwell, A. Razak, and MJ Goodwin. Braided pneumatic muscle actuators. In Proceedings of the IFAC Conference on Intelligent Autonomous Vehicles, pp. 507-512, 1993.
空気シリンダーと異なり、関節角が変化することによってトルクは非線形に変化する。
運動上の制約条件が常に不変であるという仮定の下では、空圧式エアマッスルの圧力は常に同じ均衡点での収縮割合ε（ｐ）を与えることになる。すなわち、収縮割合ε（ｐ）は、圧力ｐの関数となる。
【０１５８】
このとき、このような「空圧式エアマッスルの駆動力モデル」の逆モデルは、機械学習により学習することができる。すなわち、所望の駆動力ｆを得るための圧力ｐを求めることに相当する。内部制御装置１０は、このような圧力ｐをエアマッスルに与えるように、バルブを制御することで、所望の駆動力を出力させることになる。
【０１５９】
例えば、ハートマン等は、動力学的な学習について、以下の文献で提案している。
【０１６０】
公知文献７：Christoph Hartmann, Joschka Boedecker, Oliver Obst, Shuhei Ikemoto, and Minoru Asada. Real-time inverse dynamics learning for musculoskeletal robots based on echo state gaussian process regression. In Accepted at RSS 2012, 2012.
しかしながら、一般的には、運動上の制約条件をダイナミックに変更し、異なる外力Ｆと釣り合うので、この仮定は、上述したような外骨格型ロボットにおいては、厳密には成り立たない。
【０１６１】
したがって、この場合の均衡点での収縮割合εは、圧力と外力の関数として、ε（ｐ，Ｆ）と表現される。
【０１６２】
（腱スプリング均衡モデル）
上述したような液晶ポリマー繊維製のワイヤは、高強力高弾性率を有するので、通常は、張力に対して長さが変化しないとの近似がよく成り立つ。高強度・高弾性率繊維としては、ベクトラン（登録商標）のようなポリアリレート繊維の他に、たとえば、超高分子量ポリエチレン繊維、ＰＢＯ繊維などがある。
【０１６３】
しかし、人間（たとえば、体重は６０ｋｇある)を保持するためには、空圧式エアマッスルは、張力として、典型的には、３０００Ｎを生成する（最大値としては、たとえば、５０００Ｎが必要になる）。
【０１６４】
このような大きな力が加わると、たとえ、液晶ポリマー繊維製のワイヤ（腱ワイヤ）であっても、その伸長量が無視できず、均衡点を変化させてしまう。
【０１６５】
以下では、図１０（ｂ）に示すように、腱ワイヤの伸長の影響を考慮したモデルを検討する。
【０１６６】
そこで、液晶ポリマー繊維製の腱の力学モデルとして、以下に説明するような、線形の腱スプリング・モデルを導入する。
【０１６７】
ｆ＝ｋΔε    （Ｓ２）
すなわち、このモデルでは、腱ワイヤ（ポリアリレート繊維などの液晶ポリマー繊維製のワイヤ）は、バネに近似され、Δεは、力によって引き起こされた余分なエアマッスルの収縮であり、ｋはばね定数である。つまり、腱が力によって伸びる分、空圧式エアマッスルは、余分に収縮することが必要になる。このようにバネに近似される腱ワイヤのことを「腱スプリング」と呼ぶことにする。
【０１６８】
このような腱に相当するワイヤを、腱スプリングとして近似しない場合、空圧式エアマッスルの駆動力モデルg()は３つのパラメーター(そのうちの２つは依存関係にある)を備えた以下の二次式（３）で表される：
【０１６９】
【数１５】
2014155998.tif 000017
ここで、
２次式の係数には、以下の関係が成り立つ。
【０１７０】
【数１６】
2014155998.tif 000018
均衡点で、駆動力ｆ＊は、腱の伸長によって、以下のように減少する。
【０１７１】
【数１７】
2014155998.tif 000019
ここで、g´（）が腱スプリングの仮定の下での駆動力モデルであり、ε＊は、見かけの収縮割合であり、εは、空圧式エアマッスルの実際の収縮割合であり、Δεは、空圧式エアマッスルの実際の収縮割合への付加項である。
【０１７２】
式（Ｓ２）から、Δε＝ｆ＊／ｋ  となる。
【０１７３】
実際の収縮割合εは、上述したような外骨格型ロボットの検出機構１２４では、直接測定するのが困難である。したがって、その代りに、推定収縮割合εest（＝ε＊−Δε）を使用する。収縮割合ε＊は、以下の式のように、関節角から計算される。
【０１７４】
【数１８】
2014155998.tif 000020
ここで、Ｌ０は、関節角θ＝０のときの空圧式エアマッスルの初期長さである。つまり、「見かけの収縮割合」とは、アクチュエータにより駆動される関節の角度変化から導かれる収縮割合であって、角度センサから検知できる量に基づく収縮割合であるのに対して、「実際の収縮割合」とは、駆動力ｆ＊を生成するために、エアマッスルが現実に収縮している割合のことである。
【０１７５】
収縮割合ε＊において、駆動力ｆ＊が必要とされる場合、逆モデルｇ−１（）から、必要な圧力ｐ＊が、以下のようにして導出される。
【０１７６】
【数１９】
2014155998.tif 000021
最小二乗アルゴリズムを使用して、既知の空圧式エアマッスルのデータから、空圧式エアマッスルパラメーターａ、ｂおよびｃを評価できる。
【０１７７】
さらに、たとえば、キャリブレーション実験で、１自由度システムによって均衡点駆動力を測定する際に、バネ定数ｋについては、手動で調整することで最適値を実験的に求めることができる。
【０１７８】
運動をアシストする主要なトルク(所定の周波数以下の成分のトルク、たとえば、直流成分のトルク)は、空圧式エアマッスルによって供給され、実際のトルクτPAMsは、ロードセルにより測定することができる。上述したように、高周波トルク(所定の周波数を超える成分のトルク：交流成分のトルク)は電動モータによって生成される。
【０１７９】
したがって、モータのトルクは以下の式（Ｓ７）のように表現される。
【０１８０】
【数２０】
2014155998.tif 000022
以上のようにして、腱スプリングモデルを腱ワイヤに適用することにより、関節角と腱ワイヤに印加される駆動力を検知する構成において、空圧式エアマッスルの駆動力を正しく反映した逆モデルを、たとえば、機械学習により、構築することが可能となる。その結果、空電アクチュエータを制御する際に、その周波数応答性を向上させることが可能となる。
（１軸アシストシステムによる実験）
図１１は、１軸アシスト装置の実験セットアップの外観を示す図である。
【０１８１】
図１１（ａ）では被験者（健常者）が装着しており、２つの筋電位電極と拘束部が示されている。図１１（ｂ）では拮抗で配置されている２つのＰＡＭ（PAM１とPAM２）が式（Ｓ６）による重力補償トルク制御により、おもりを静止させている様子を示している。
【０１８２】
図１２は、筋電位計測のための体表電極の位置の例を示す図である。
【０１８３】
図１２に示すように、それぞれ上腕二頭筋と上腕三頭筋を計測対象としている。アンプの倍率は１００００倍に設定し、−１０[V]から１０[V]で計測可能なＡＤ変換器により、筋電位を計測する。
【０１８４】
図１３は、筋電位計測により推定された被験者の関節駆動トルクτhと１軸アシスト装置の駆動トルクτrを示す図である。
【０１８５】
０．５［Ｈｚ］程度でおもりを上下させることを被験者のタスクとした。図１３（ａ）はパラメータ（α1，α2）=（０．９，０）による結果であり、式（３）の第２項の調整項がない場合の実験結果であり、従来手法に相当する。
【０１８６】
図１３（ｂ）は、パラメータ（α1，α2）=（０．４５，０．４５）による結果であり、本実施の形態の手法の実装である。
【０１８７】
図１３（ａ）に示すように、上述のパラメータ設定ではヒト−ロボット間のインタラクションは、図１３（ａ）の７［ｓｅｃ］以降では、図１３（ｂ）に比べて不安定である。これは、運動の方向が変化する際に発生するＰＡＭのモデル誤差やベルトの緩みなどが原因となり、筋電位より推定されるトルクフィードバックとそれに伴うヒトの筋電位の変動が不安定化したものと考えられる。図１３（ｂ）のコントローラの内部にはフィードフォワード的にインタラクション力を考慮するアシスト力の調整項が入っており、アシスト装置が発生するトルクの急激な変動を押さえる結果、安定化された。
【０１８８】
すなわち、α1とα2はアシストトルクの応答とシステムの安定化を調整するパラメータであり、これらはトレードオフの関係にある。α1に比べてα2を大きく調整することで、アシストのインタラクションの安定化を計ることができ、逆にα2に比べてα1を大きくすることで筋の活性度に対するアシストトルクの応答性が向上する。
【０１８９】
以上説明したように、本実施の形態では、筋電位フィードバックに基づくアシスト力のフィードバック制御において、インタラクションにより人の筋張力の減少が起こることに注目し、これを調整するフィードバックコントローラを用いる。
【０１９０】
一例として、１軸アシスト装置により従来手法と本実施の形態の手法を比較した結果により、実験的にも、本実施の形態の制御方式の有効性が示された。
【０１９１】
本実施の形態で示されたように、フィードフォワード的にインタラクションの発生力を考慮し、それを調整項としてフィードバックすることで、筋電位フィードバックが安定化できる。
【０１９２】
なお、インタラクションを現実に計測し、フィードバックすることで、インタラクションを考慮した運動アシストを行ってもよい。
【０１９３】
今回開示された実施の形態は、本発明を具体的に実施するための構成の例示であって、本発明の技術的範囲を制限するものではない。本発明の技術的範囲は、実施の形態の説明ではなく、特許請求の範囲によって示されるものであり、特許請求の範囲の文言上の範囲および均等の意味の範囲内での変更が含まれることが意図される。
【符号の説明】
【０１９４】
１  外骨格型ロボット、１０  内部制御装置、１１ Ｉ／Ｆ部、１２  外骨格、２０  外部制御装置、１２１  ベース、１２２  下半身、１２３  能動関節、１２４  検出機構、１３１  記録処理部、１３２  記憶装置、１３３  計測装置、１３４  制御部、１３５  出力装置、３０２，３０４  エアマッスル、３１０  プーリー付回転関節、３１２  電動モータ。
【図５】
2014155998.tif 000023
【図６】
2014155998.tif 000024
【図９】
2014155998.tif 000025
【図１】
2014155998.tif 000026
【図２】
2014155998.tif 000027
【図３】
2014155998.tif 000028
【図４】
2014155998.tif 000029
【図７】
2014155998.tif 000030
【図８】
2014155998.tif 000031
【図１０】
2014155998.tif 000032
【図１１】
2014155998.tif 000033
【図１２】
2014155998.tif 000034
【図１３】
2014155998.tif 000035