(57)【要約】      （修正有）【課題】簡単かつ迅速に人間の状態を推定できるようにした状態推定システムを提供する。【解決手段】四隅にロードセルを備えたフォースプレート２を用いて、起立する被検体のフォースプレート２の平面方向に作用する剪断力と、圧力中心を第一演算部５で算出する。算出した値を用いて、事前に用意された被検体のパラメーターを用いて、矢状面における１リンクモデルを用いて、足関節のモーメントや角加速度、角度などを算出し、２リンクモデルを用いて下半身と上半身の質量中心における角加速度、角度、股関節のモーメントを算出する。それら算出された値と前額面のモデルを用いて、各部の角加速度、角度、モーメントを算出する。【選択図】図２
2020069078.tif 000002
【特許請求の範囲】
【請求項１】
複数のロードセルを備えたフォースプレートと、
当該ロードセルからの出力値に基づいて、当該フォースプレートの平面に沿った被検体からの剪断力と被検体の立位状態における圧力中心を算出する第一演算部と、
当該第一演算部で算出された剪断力と圧力中心を用いて、被検体の各部位の角度、角速度、角加速度、モーメントの少なくとも一つを算出する第二演算部と、
当該第二演算部で算出された値を出力する出力部と、
を備えてなる状態推定システム。
【請求項２】
前記第二演算部が、あらかじめ足部質量、身体部質量、フォースプレートから各関節までの高さ、足部の質量中心までの高さ、身体部の質量中心の高さをパラメーターとして演算されるものである請求項１に記載の状態推定システム。
【請求項３】
前記第二演算部が、人間の矢状面における足関節のモーメントを算出するものである請求項２に記載の状態推定システム。
【請求項４】
第二演算部が、矢状面における足関節のモーメントを算出した後に、矢状面における下半身と上半身の２つの質量中心の加速度と変位、速度を算出していくものである請求項２に記載の状態推定システム。
【請求項５】
前記第二演算部が、矢状面における下半身と上半身の２つの質量中心の加速度と変位、速度を算出した後に、股関節のモーメントを算出していくものである請求項２に記載の状態推定システム。
【請求項６】
前記第二演算部が、前額面における足関節のモーメントを算出するものである請求項２に記載の状態推定システム。
【請求項７】
第二演算部が、前額面における足関節のモーメントを算出した後に、前額面における下半身と上半身の２つの質量中心の加速度と変位、速度を算出していくものである請求項２に記載の状態推定システム。
【請求項８】
前記第二演算部が、この前額面における下半身と上半身の２つの質量中心の加速度と変位、速度を算出した後、股関節のモーメントを推定するようにしたものである請求項２に記載の状態推定システム。
【請求項９】
前記フォースプレートか゛、両足に独立して設けられたものて゛ある請求項1に記載の状態推定システム。
【請求項１０】
複数のロードセルを備えたフォースプレートを用いて、フォースプレート上に立脚した被検体の状態を推定する状態推定方法において、
前記ロードセルの値を用いて、当該フォースプレートの平面に沿った被検体からの剪断力と被検体の圧力中心を算出する第一演算ステップと、
当該第一演算部で算出された剪断力と圧力中心を用いて、被検体の各部位の角度、角速度、角加速度、モーメントの少なくとも一つを算出する第二演算ステップと、
当該第二演算部で算出された値を出力するステップと、
を備えた状態推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、人間の各部位の動き、姿勢、各関節回りのモーメントなどを推定できるようにした状態推定システムに関するものであって、より具体的には、フォースプレートに起立した状態で下半身、上半身の動きや姿勢、足関節や股関節のモーメントなどを推定できるようにした状態推定システムに関するものである。
【背景技術】
【０００２】
従来より、起立している人間の各部位の角度や位置などを推定できるようにしたシステムとして、下記の特許文献１に記載されるようなシステムなどが存在している。
【０００３】
この特許文献１に記載されるシステムは、人間の左右の少なくとも一方の大腿部に装着された第一加速度センサーと、この人間の上半身に装着された第二加速度センサーとを備えてなるものであって、これら第一加速度センサーや第二加速度センサーで計測された加速度を用いて、人間の下半身の腰や膝などの関節の屈曲状態などを推定できるようにしたものである。
【０００４】
また、別の従来例として、人間の姿勢を推定する場合、人間の各部位にマーキングを施し、そのマーキングをモーションキャプチャーで取得することによって、マークの変位情報から各関節の位置や各部位の動きなどを検出できるようにしたシステムもある。
【０００５】
このようなシステムを用いれば、人間の足や脚、上半身などの角度や動きなどを推定することができるため、例えば、これを患者に適用した場合に、その患者のリハビリ時における回復状況などを客観的に判断することができるようになる。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２０１８−０１５０２３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、上述のようなシステムを用いて人間の各部における位置や動きなどの状態を推定する場合、次のような問題がある。
【０００８】
すなわち、上述のようなシステムでは、人間の複数の部位に加速度センサーやモーションキャプチャーのためのマークを付さなければならず、加速度センサーを装着するための準備や、マークを付するための準備などに手間がかかるといった問題があった。また、加速度センサーやモーションキャプチャーを用いた場合は、身体表面の位置や加速度しか検出することができず、質量中心の変位を推定する場合、どうしても誤差を生じてしまうといった問題があった。
【０００９】
そこで、本発明は、上記課題に着目して、簡単かつ迅速に被検体の各部位における位置や動き、モーメントなどの状態を推定できるようにした状態推定システムを提供することを目的とする。
【課題を解決するための手段】
【００１０】
すなわち、本発明は上記課題を解決するために、複数のロードセルを備えたフォースプレートと、当該ロードセルからの出力値に基づいて、当該フォースプレートの平面に沿った被験体からの剪断力と被検体の立位状態における圧力中心を算出する第一演算部と、当該第一演算部で算出された剪断力と圧力中心を用いて、被検体の各部位の角度、角速度、角加速度、モーメントの少なくとも一つを算出する第二演算部と、当該第二演算部で算出された値を出力する出力部とを備えるようにしたものである。
【００１１】
また、このような演算を行う場合、あらかじめ、足部質量、身体部質量、フォースプレートから各関節までの高さ、足部の質量中心までの高さ、身体部の質量中心の高さをパラメーターとして演算できるようにする。
【００１２】
そして、前記第二演算部で演算を行う場合、人間の矢状面における足関節のモーメントを算出する。
【００１３】
そして、このように足関節のモーメントを算出した後に、矢状面における下半身と上半身の２つの質量中心の加速度と変位、速度を算出していくようにする。
【００１４】
同様にして、このように矢状面における下半身と上半身の２つの質量中心の加速度と変位、速度を算出した後、股関節のモーメントを推定していくようにする。
【００１５】
また、別の発明では、前記第二演算部で、前額面における足関節のモーメントを算出する。
【００１６】
このとき、前額面における足関節のモーメントを算出した後に、前額面における下半身と上半身の２つの質量中心の加速度と変位、速度を算出していくようにする。
【００１７】
そして、この前額面における下半身と上半身の２つの質量中心の加速度と変位、速度を算出した後、股関節のモーメントを推定していくようにする。
【００１８】
また、前記フォースプレートを左右の足に対応して独立して設けることによって、下半身のモーメントを別々に求めるようにすることもできる。
【発明の効果】
【００１９】
このように、本発明によれば、フォースプレーの剪断力と圧力中心だけで被検体の各部位の角度や角加速度、モーメントなどの状態を検出することができるため、従来のようにモーションキャプチャーのためのマークを施すことなく、簡単に被検体の状態を推定していくことができるようになる。
【図面の簡単な説明】
【００２０】
【図１】本発明の一実施の形態におけるフォースプレートに起立した状態を示す図
【図２】同形態における機能ブロック図
【図３】同形態の矢状面における矢状面における１リンクモデルを示す図
【図４】同形態の矢状面における矢状面における２リンクモデルを示す図
【図５】同形態の前額面における状態を検出する場合の図
【図６】同形態の推定方法を示すフローチャート
【図７】従来例と本発明の矢状面における１リンクモデルの比較例
【図８】従来例と本発明の矢状面における２リンクモデルの比較例
【図９】従来例と本発明の矢状面における股関節モーメントの比較例
【図１０】従来例と本発明のフォースプレートを揺動させた場合における、矢状面における１リンクモデルの比較例
【図１１】従来例と本発明のフォースプレートを揺動させた場合における、矢状面における２リンクモデルの比較例
【図１２】従来例と本発明のフォースプレートを揺動させた場合における、矢状面における股関節モーメントの比較例
【図１３】従来例と本発明の前額面における１リンクモデルの比較例
【図１４】従来例と本発明の前額面における２リンクモデルの比較例
【図１５】従来例と本発明の前額面におけるモーメントの比較例
【発明を実施するための最良の形態】
【００２１】
以下、本発明の一実施の形態について図面を参照しながら説明する。
【００２２】
この実施の形態における状態推定システム１は、フォースプレート２上に起立する人間の各部位の位置や速度、加速度、モーメントなどの状態（以下、「状態」と略する）を推定できるようにしたものであって、図１や図２に示すように、複数のロードセル４を備えたフォースプレート２と、このフォースプレート２に起立した被検体である人間からの水平面方向に沿った剪断力と、フォースプレート２上の圧力中心を算出する第一演算部５と、この第一演算部５で演算された剪断力と圧力中心を用いて、被検体の各部の角度、加速度、角加速度、モーメントなどの状態を算出する第二演算部６とを備えるようにしたものである。以下、この実施の形態における状態推定システム１について詳細に説明する。なお、この実施の形態において、説明の関係上、図１に示すように、フォースプレート２の前後方向をX軸方向、左右方向をY軸方向、上下方向をZ軸方向として説明し、被検体にX軸方向に向かって起立してもらうものとする。また、矢状面とはXZ平面、前額面とはXY平面を示すものとする。
【００２３】
まず、フォースプレート２は、矩形状のプレート３の四隅下方にロードセル４を備えて構成される。このロードセル４は、図示しないひずみゲージを起歪体に貼り付けて構成されるものであって、その起歪体の変形をひずみゲージを用いて検出して、XYZ方向に作用する荷重や、そのXYZ軸回りのモーメントなどを検出できるようになっている。そして、これらのロードセル４を用いて、プレート３に作用する被検体からのXY平面に沿った剪断力や、XY平面に作用するZ軸方向からの荷重に対する圧力中心などを算出できるようになっている。また、この実施の形態では、人間を揺動させて各部位の状態の変化を検出できるようになっており、フォースプレート２の動きを検出する加速度センサーや、そのフォースプレートの移動状態を検出するためのロータリーエンコーダーなどの揺動検出手段８が設けられている。なお、ここでは、加速度センサーやロータリーエンコーダーで揺動検出手段８を構成するようにしているが、モーターでフォースプレート２を前後揺動させる場合は、このモーターへの電圧入力値を揺動検出手段８としてもよい。
【００２４】
このようなフォースプレート８に起立された被検体の重心は、動作を行う瞬間において、前後方向に移動し、これに伴って足からフォースプレート２に剪断力が作用するようになる。第一演算部５では、この剪断力を算出するとともに、その被検体のZ軸方向の荷重などによって圧力中心を算出できるようにしている。このとき、剪断力を算出する場合は、四隅に設けられたロードセル４のX軸方向の荷重の合計を算出し、サンプリング時間における前後の荷重合計値の差分を剪断力として算出する。一方、圧力中心を算出する場合は、四隅に設けられたロードセル４のZ軸方向の荷重を検出し、圧力中心の座標を未知数として、モーメントの釣り合い式を用いて圧力中心の座標を算出していく。そして、このように算出された剪断力や圧力中心を第二演算部６に出力する。
【００２５】
第二演算部６では、この剪断力や圧力中心を用いて、まず、矢状面における被検体の足関節の角度や角加速度、足関節に作用するモーメントを算出する。これらの値を算出するに際して、事前に被検体の足部の質量、身体部の質量、身体部の質量中心まわりの慣性モーメント、足関節の地面からの高さ、足部の質量中心の高さ、足関節から身体部質量中心までの長さをパラメーターとして入力しておくものとする。これらのパラメーターは、被検体から直接測定してもよいし、あるいは被検体の身長や体重から過去の統計値を用いて推測してもよい。
【００２６】
これらのパラメーターや剪断力、圧力中心を用いて、矢状面における足と脚の関節の角度や角加速度、足関節に作用するモーメントを算出する場合、図３の１リンクモデルを用いて、次のようにして算出する。
【００２７】
＜矢状面における状態＞
【００２８】
まず、足部におけるX方向、Z方向、および、足関節まわりのモーメントは、次の数１によってそれぞれ算出される。なお、パラメーターとして、mfを足部質量、mbを身体部質量、Jbを身体部質量中心まわりの慣性モーメント、Lfを足関節の高さ、lfを足部質量中心の高さ、l1を足関節から下半身質量中心までの長さ、Mを被検体の体重とする。
【００２９】
【数１】
2020069078.tif 000003
【００３０】
一方、身体部におけるX方向、Z方向、および、足関節まわりのモーメントは、次の数２によってそれぞれ算出される。
【００３１】
【数２】
2020069078.tif 000004
【００３２】
そして、上記数１と数２の第一式（X軸方向の釣り合い式）の和より、次の数３が得られる。
【００３３】
【数３】
2020069078.tif 000005
【００３４】
また、上記数１と数２の第三式（モーメントの釣り合い式）の和より、次の数４が得られる。
【００３５】
【数４】
2020069078.tif 000006
【００３６】
これらの数３と数４の右辺はすべてパラメーターや計測値であり、左辺は算出対象となる足関節の角度と角加速度の連立１次方程式になる。そして、これを解くことにより、数５や数６のように足関節の角加速度と角度を得ることができる。
【００３７】
【数５】
2020069078.tif 000007
【００３８】
【数６】
2020069078.tif 000008
【００３９】
また、足関節に作用するモーメントについては、数１の第三式（モーメントの釣り合い式）より、床反力に作用する力（Z軸方向の力）のみから数７のようにして求められる。
【００４０】
【数７】
2020069078.tif 000009
【００４１】
次に、このように算出された足関節の角度、角加速度、モーメントから、今度は、図４に示す２リンクモデルを用いて被検体の下半身、上半身の状態を算出していく。
【００４２】
人間の下半身、上半身のそれぞれの角度や角加速度、あるいは、下半身と上半身の関節部分である股関節回りのモーメントを算出する場合、次のようにして算出する。なお、ここでは、パラメーターとして、mfを足部質量、m1、m2をそれぞれ下半身および上半身の質量、J1、J2を下半身および上半身の各質量中心まわりの慣性モーメント、Lfを足関節の高さ、lfを足部質量中心の高さ、l1を足関節から下半身質量中心までの長さ、L1を足関節から股関節の長さ、l2を股関節から上半身質量中心までの長さとする。これらのパラメーターは、同様に被検体の身長や体重などを用い、過去の統計値などから推定されるものとする。
【００４３】
パラメーターや剪断力、圧力中心、および、数５から数７を用いて算出された足関節の角度、角加速度、足関節まわりのモーメントを用い、２リンクモデルである図４のように、下半身のX軸方向の釣り合い式、Z軸方向の釣り合い式、足関節まわりのモーメントは、次の数８のように成立する。
【００４４】
【数８】
2020069078.tif 000010
【００４５】
一方、上半身におけるX軸方向の釣り合い式、Z軸方向の釣り合い式、足関節まわりのモーメントは、次の数９のように成立する。
【００４６】
【数９】
2020069078.tif 000011
【００４７】
そして、数１や数８、数９の第一式（X軸方向の釣り合い式）の和より、X軸方向の運動方程式が数１０のように得られる。
【００４８】
【数１０】
2020069078.tif 000012
【００４９】
一方、数１や数８、数９の第二式（Z軸方向の釣り合い式）の和より、Z軸方向の運動方程式が数１１のように得られる。
【００５０】
【数１１】
2020069078.tif 000013
【００５１】
ここで、Mは被検体の体重を示すものである。そして、数１や、数８、数９の第三式（モーメントの釣り合い式）に数１１を代入すると、回転に関する運動方程式が、数１２のように得られる。
【００５２】
【数１２】
2020069078.tif 000014
【００５３】
ここで、図３と図４のリンクモデルの質量中心は、次の数１３のように表される。
【００５４】
【数１３】
2020069078.tif 000015
【００５５】
この関係を数１２の左辺第３項と第４項に代入すると、次の数１４のようになる。
【００５６】
【数１４】
2020069078.tif 000016
【００５７】
この数１４の右辺第１項は数６で推定していると仮定すると、数１０と数１４の右辺はすべて既知であり、これらは、下半身の角加速度と上半身の角加速度の２元連立１次方程式となる。そして、これを解くことにより、下半身と上半身の角加速度が得られ、このように得られた角加速度を積分することによって下半身と上半身の角度がそれぞれ計算される。
【００５８】
このように図４の２リンクモデルの角加速度と角度が算出できた場合、股関節まわりのモーメントは、数１と数８の第三式の和に数１と数８の第一式から算出したRhxを代入することにより、次の数１５で表される。
【００５９】
【数１５】
2020069078.tif 000017
【００６０】
このように、上半身と下半身の角加速度および角度を推定することができれば、数１５により自動的に股関節まわりのモーメントNhを得ることができるようになる。
【００６１】
そして、このように算出された、足関節や下半身、上半身の角加速度、角速度、モーメントなどを出力し、被検体の状態を判断できるようにする。なお、このような値を算出する場合、フォースプレートが加速している場合は、揺動検出手段８で検出された方向の加速度などを考慮して算出していく。
【００６２】
＜前額面における状態＞
【００６３】
次に、前額面（XY平面）におけるモデルを図５に示す。このモデルでは、身体を両足、両脚部、骨盤部、上半身の６つのセグメントに構成し、両足関節、両股関節、腰部に５つのジョイントを持つモデルとしている。ここで、（Xf,r,Zf,r）、(Xf,l,Zf,l)を左右足部の質量中心、(Xl,r,Zl,r)、(Xl,l,Zl,l)を左右脚部の質量中心、(Xp,Zp)を骨盤部の質量中心、(Xu,Zu)を上半身の質量中心の座標とする。身体のパラメーターとして、mfを足部の質量、mlを脚部質量、mwを骨盤部の質量、muを上半身の質量、Jlを脚部質量中心まわりの慣性モーメント、Juを上半身質量中心まわりの慣性モーメントとする。また、Lfを足関節の高さ、lfを足部質量中心の高さ、llを足関節から脚部質量中心までの長さ、Llを足関節から股関節までの長さ、lpを股関節から骨盤部質量中心までの長さ、Lpを股関節から腰部ジョイントまでの長さ、luを腰部ジョイントから上半身質量中心までの長さとする。そして、θl,r、θl,lは鉛直軸からみた両脚部の姿勢角度、θuを上半身の角度とする。
【００６４】
両足のスタンス幅は、股関節幅wと一致すると仮定し、両脚部は常に平行で角度は等しいと仮定する。これにより、θl=θl,r＝θl,lとなる。また、膝関節による上下動はなく、姿勢角の変化は微小であるとする。これにより、線形システムと考えることができ、各関節ジョイントには、図５に示すような拘束力とモーメントが作用していると考える。
【００６５】
このような前提のもとで、左右の足、左右の脚、骨盤、上半身は次の数１６の関係式が成立する。
【００６６】
【数１６】
2020069078.tif 000018
【００６７】
次に、各セグメントの運動方程式をフリーボディダイアグラムに基づいて求める。ただし、以下の式では、セグメントの水平加速度は、数１６より支持面加速度と各セグメントの角加速度に変換する。また、左右の脚部姿勢角は常に等しいと仮定しているため、角度および角加速度も同じとして計算する。
【００６８】
これにより、右足部の水平方向、鉛直方向、モーメントの式は次の数１７でそれぞれ示される。
【００６９】
【数１７】
2020069078.tif 000019
【００７０】
同様に、左足部の水平方向、鉛直方向、モーメントの式は次の数１８でそれぞれ示される。
【００７１】
【数１８】
2020069078.tif 000020
【００７２】
同様に、右脚部の水平方向、鉛直方向、モーメントの式は次の数１９でそれぞれ示される。
【００７３】
【数１９】
2020069078.tif 000021
【００７４】
同様に、左脚部の水平方向、鉛直方向、モーメントの式は次の数２０でそれぞれ示される。
【００７５】
【数２０】
2020069078.tif 000022
【００７６】
同様に、骨盤部の水平方向、鉛直方向、モーメントの式は次の数２１でそれぞれ示される。
【００７７】
【数２１】
2020069078.tif 000023
【００７８】
同様に、上半身の水平方向、鉛直方向、モーメントの式は次の数２２でそれぞれ示される。
【００７９】
【数２２】
2020069078.tif 000024
【００８０】
上記式に基づいて、床反力と姿勢の関係式を順に示す。まず、数１７から数２２の第一式（水平方向の関係式）より、Y軸に関する並進運動の運動方程式は次のようになる。
【００８１】
【数２３】
2020069078.tif 000025
【００８２】
ここで、Mは体重を示し、Ry=Ry,r=Ry,lとする。次に、数１７から数２２の第二式（鉛直方向の関係式）より、Z軸に関する力の釣り合いは次のようになる。
【００８３】
【数２４】
2020069078.tif 000026
【００８４】
次に、回転に関する運動方程式を示す。式の整理のために、Rhy,r、Rhy,l、Rhz,r、Rhz,l、が必要となるので、まず、数１７と数１９の第一式より、Rhy,rを得る。
【００８５】
【数２５】
2020069078.tif 000027
【００８６】
同様に、数１８と数２０の第二式より、Rhy,lを得る。
【００８７】
【数２６】
2020069078.tif 000028
【００８８】
数１７と数１９の第二式より、Rhz,rを得る。
【００８９】
【数２７】
2020069078.tif 000029
【００９０】
数１８と数２０の第二式より、Rhz,lを得る。
【００９１】
【数２８】
2020069078.tif 000030
【００９２】
また、Rwyも必要になるので、数１７から数２１の第一式より次式を得る。
【００９３】
【数２９】
2020069078.tif 000031
【００９４】
回転に関する式の総和として、数１７から数２２の第三式の総和を計算し、Rhy,r、Rhy,l、Rhz,r、Rhz,lの部分に数２５から数２８を代入すると、次式が得られる。
【００９５】
【数３０】
2020069078.tif 000032
【００９６】
ここで、mb=M-2mfは身体部の質量である。１リンクモデルと２リンクモデルの質量中心の関係より、次式が成り立つ。
【００９７】
【数３１】
2020069078.tif 000033
【００９８】
そして、これを数３０に代入すると、次式が得られる。
【００９９】
【数３２】
2020069078.tif 000034
【０１００】
これにより、図３の１リンクモデルの姿勢角θbが既知の場合、数２３と数３２の連立１次方程式を解くことで上半身および下半身の角速度を求めることができる。また、身体が床面から受けるモーメントはNsであり、プレートが一枚である場合は、計測されたX軸まわりのモーメントMxを反転させたものになり、Ns=-Mxとなる。
【０１０１】
一方、フォースプレートが２枚の場合は、次式で表される。
【０１０２】
【数３３】
2020069078.tif 000035
【０１０３】
ここで、wはスタンス幅、Rz,r、Rz,lは圧力中心を左右の足の仮想的な中心一から相対変位で表したものである。Yp1、Yp2はそれぞれのフォースプレートの圧力中心、Rz1、Rz2は垂直反力、hは２枚のフォースプレートの原点間の距離である。この添字1は右足、添字2は左足のフォースプレートの出力値を示す。
【０１０４】
数２３と数３０がこの２リンクモデルの基礎式となる。これを１リンクモデルとするために拘束条件を与えると次式が得られる。
【０１０５】
【数３４】
2020069078.tif 000036
【０１０６】
【数３５】
2020069078.tif 000037
【０１０７】
以上の通り、前額面においても矢状面と同じ理論で展開することができる。すなわち、矢状面の１リンクモデルの運動方程式である数３と数４が、前額面の１リンクモデルの運動方程式である数３４と数３５に対応しており、矢状面における２リンクモデルの運動方程式である数１０と数１２が、前額面の２リンクモデルの運動方程式である数２３と数３０に対応する。
【０１０８】
＜前額面の関節モーメントの推定＞
【０１０９】
次に、フォースプレート２に基づく姿勢推定の結果を用いて、前額面の関節モーメントの推定を行う。フォースプレート２が一枚の場合は、「下半身モーメント（両股関節・両足関節モーメントの和）」と「腰部まわりのモーメント」が推定できる。一方、片足ずつ２枚のフォースプレートで計測する場合には、「左足関節モーメント」「右足関節モーメント」「左股関節モーメント」「右股関節モーメント」「腰部まわりのモーメント」の５つのモーメントが推定できる。
【０１１０】
まず、左右の足関節モーメントは、数１７と数１８のそれぞれの第三式の和から次のように求められる。
【０１１１】
【数３６】
2020069078.tif 000038
【０１１２】
【数３７】
2020069078.tif 000039
【０１１３】
また、右股関節モーメントは、数１７と数１９のそれぞれの第三式の和を計算し、Rhy,r、Rhz,rの部分に数２５と数２９を代入すると、次のようになる。
【０１１４】
【数３８】
2020069078.tif 000040
【０１１５】
左股関節モーメントは、数１８と数２０のそれぞれの第三式の和で計算できる。
【０１１６】
【数３９】
2020069078.tif 000041
【０１１７】
このように２枚のフォースプレート２があれば、左右の関節モーメントを独立に推定することができる。
【０１１８】
次に、左右の足関節トルクの和は、数３６と数３７の和から次のようになる。
【０１１９】
【数４０】
2020069078.tif 000042
【０１２０】
また、左右の股関節トルクの和は、数３８と数３９の和から次のようになる。
【０１２１】
【数４１】
2020069078.tif 000043
【０１２２】
数４０と数４１では、左右のフォースプレート２の垂直反力の差Rz,r-Rz,lが必要となるため、股関節と足関節を区別して計算することができない。一方、下半身全体のモーメントτlは数４０と数４１の差より、
【０１２３】
【数４２】
2020069078.tif 000044
となるので、こちらは１枚のフォースプレート２の計測から推定することができる。
【０１２４】
最後に腰部まわりのモーメントを求める。この計算では、数１７から数２１の第一式の和で次式が得られる。
【０１２５】
【数４３】
2020069078.tif 000045
【０１２６】
数１７から数２１の第三式の計算をしたうえで、数４３を代入すると、次式のようになる。
【０１２７】
【数４４】
2020069078.tif 000046
【０１２８】
こちらも１枚のフォースプレート２の計測値によって計算することができるようになる。
【０１２９】
次に、このように構成された状態推定システム１を用いて、人間の各部における角度、角加速度、モーメントなどを推定していく方法について、図６のフローチャートを用いて説明する。
【０１３０】
まず、事前に被検体である人間の身長や体重などを計測しておき、演算に必要な各パラメーターを算出しておく。そして、その状態でフォースプレート２に起立してもらい、プレート３の四隅に設けられたロードセル４を用いて、プレート３の平面に沿った剪断力や、プレート３にかかる圧力中心を第一演算部５を用いて算出する（ステップＳ１）。
【０１３１】
そして、このように既知数であるパラメーターや剪断力を用いて、矢状面である１リンクモデルの足関節の角加速度や角度、足関節に作用するモーメントを数５、数６、数７を用いて求める（ステップＳ２）。
【０１３２】
次に、このように足関節の角加速度や角度、モーメントが算出できた状態で、次に２リンクモデルを用いて、下半身や上半身の角加速度、角度、股関節のモーメントなどを数１０、数１４、数１５から求める（ステップＳ３）。
【０１３３】
今度は、このように算出された１リンクモデルの角度θbを用い、数２３や数３２を用いて前額面における角加速度を求め、これを積分することにより角度を求める。また、身体が床面から受けるモーメントは、圧力中心におけるX軸まわりのモーメントであるため、Ns=-Mxとなり、フォースプレート２が２枚の場合は、数３３からモーメントが求められる（ステップＳ４）。
【０１３４】
そして、最後に前額面の関節モーメントを数３９、数３９、数４２、数４４を用いて算出する（ステップＳ５）。このように矢状面の１リンクモデルから２リンクモデル、また、前額面におけるモデルの順に計算をしていき、各部の角加速度、角度、モーメントを算出して出力する（ステップＳ６）。
【０１３５】
このように上記実施の形態によれば、複数のロードセル３を備えたフォースプレート２と、当該ロードセル３からの出力値に基づいて、当該フォースプレート２の平面に沿った被検体からの剪断力と被検体の立位状態における圧力中心を算出する第一演算部５と、当該第一演算部５で算出された剪断力と圧力中心を用いて、被検体の各部位の角度、角速度、角加速度、モーメントの少なくとも一つを算出する第二演算部６と、当該第二演算部６で算出された値を出力する出力部７とを備えるようにしたので、フォースプレート２の剪断力と圧力中心だけから人間の各部の状態を検出することができ、従来のようにモーションキャプチャーのマークなどを施すことなく、簡単に被検体の状態を推定していくことができるようになる。
【実施例１】
【０１３６】
上述のフォースプレートのみを用いた状態方法による結果と、従来例であるモーションキャプチャーとフォースプレートを用いて検出した結果の比較例を示す。実証実験においては、30秒間計測で行い、モーションキャプチャーで計測された値を太線、上記実施の形態で推定された値を細線で示す。サンプリング周波数は、ともに100Hzとし、モーションキャプチャーで計測された角速度および角加速度の算出は、数値微分により行った。被検体として、身長180cm、体重73kgの男性に静止床面に起立してもらった。
【０１３７】
図７の上から順に矢状面における１リンクモデルの角度、角速度、角加速度を示す。図７からわかるようにモーションキャプチャーで計算された値と近似していることがわかる。
【０１３８】
また、２リンクモデルの角度、角速度、角加速度を図８に示す。同様に、図８においても、モーションキャプチャーで計算された値と近似していることがわかる。
【０１３９】
さらに、股関節まわりのモーメントを図９に示す。図９においても、モーションキャプチャーで計算された値と近似していることがわかる。
【０１４０】
次に、フォースプレートに振動装置を備えて揺動実験を行った場合の実証結果を示す。フォースプレートの揺動は、0.1Hzから1.0Hzの間で複数の揺動波形を重ね合わせて不規則に揺動させるようにした。
【０１４１】
この図１０に示すように、１リンクモデルにおいても、従来のモーションキャプチャーを用いた手法と近似していることが分かり、２リンクモデル（図１１）や股関節モーメント（図１２）においても非常に近似していることが分かった。
【０１４２】
一方、前額面の実証実験を行った。この場合も同様にモーションキャプチャーとの比較例を示す。１リンクモデルによる値を図１３、２リンクモデルによる値を図１４、関節モーメントの値を図１５に示す。これらの値から分かるように、フォースプレートのみを用いた推定方法によっても、モーションキャプチャーによる推定と非常に近似していることが分かった。
【符号の説明】
【０１４３】
１・・・状態推定システム
２・・・フォースプレート
３・・・プレート
４・・・ロードセル
５・・・第一演算部
６・・・第二演算部
７・・・出力部
８・・・揺動検出手段
【図１】
2020069078.tif 000047
【図２】
2020069078.tif 000048
【図３】
2020069078.tif 000049
【図４】
2020069078.tif 000050
【図５】
2020069078.tif 000051
【図６】
2020069078.tif 000052
【図７】
2020069078.tif 000053
【図８】
2020069078.tif 000054
【図９】
2020069078.tif 000055
【図１０】
2020069078.tif 000056
【図１１】
2020069078.tif 000057
【図１２】
2020069078.tif 000058
【図１３】
2020069078.tif 000059
【図１４】
2020069078.tif 000060
【図１５】
2020069078.tif 000061