(57)【特許請求の範囲】
【請求項１】
垂直荷重からCOPを計測する重心計を用い、Δｔ間の立位状態における重心計上の力の作用点（COP）の変位量ΔCOPを測定するステップと、重心位置の移動速度（ＶＣＯＰ）をＶＣＯＰ＝ΔCOP／Δｔで、重心移動加速度（αＣＯＰ）をαＣＯＰ＝ΔＶＣＯＰ／Δｔで算出するステップと、ｔ時の重心に作用する仮想の力Ｆｘｔ，ＦｙｔとＦｚｔを次式で求めるステップとから成る重心計を用いて身体重心に作用する仮想の力を算出する方法。
【数１】
0005921892.tif 000012
ただし、Ｇは重力加速度、Ｆｔは時刻ｔの荷重である。
【請求項２】
重心計には左足を載せるプレートＬと右足を載せるプレートＲを組み合わせたものを用い、請求項１に記載の方法によってそれぞれ仮想床反力を求めるステップと、左右のプレートＬ，Ｒ上の力の作用点（COP）ＰＬとＰＲを結ぶ線分を左右のプレートの荷重比で内分して合成COPを算出するステップと、前記求めた左右の仮想床反力を前記の線分に対する合成COPの位置比率でＸ，Ｙ成分ごとに算出するステップとからなる重心計を用いて身体重心に作用する仮想の力を算出する方法。
【請求項３】
垂直荷重からCOPを計測する重心計を用い、人間モデルから身体重心（COG）の床面からの高さｌと軸（足関節）の粘性摩擦係数Ｃ＝［kgm2/S］ 及び慣性モーメントを参照し、立位状態を振り子モデルとして演算するステップと、体重をＷ，時刻ｔ の垂直力成分を Ｆｔ，時刻ｔ の前後力成分を Ｆy、左右力成分をＦｘとして、次の振り子の式からθの値を算出し、極座標（ｌ，θ）から身体重心に作用する仮想の力とCOGの空間の座標値ｚ、ｙを求めるステップと、同様の手法でCOGの空間の座標値ｚ、ｘを求めるステップと、からなる重心計を用いて身体重心位置の近似値を求める方法。
【数３】
0005921892.tif 000013
【請求項４】
垂直荷重からCOPを計測する重心計を用い、人間モデルから身体重心（COG）の床面からの高さｌと軸（足関節）の粘性摩擦係数Ｃ＝［kgm2/S］ 及び慣性モーメントを参照し、立位状態からしゃがみ込んだ状態を振り子モデルとして演算するステップと、体重をＷ，時刻ｔ のしゃがみ込み時にかかる垂直力成分を Ｆｔ’，時刻ｔ の前後力成分を Ｆy、左右力成分をＦｘとして、次の振り子の式からθの値を算出し、極座標（ｌ，θ）から身体重心に作用する仮想の力とCOGの空間の座標値ｚ、ｙを求めるステップと、同様の手法でCOGの空間の座標値ｚ、ｘを求めるステップと、からなる重心計を用いてしゃがみ込んだ状態でのCOGの三次元空間位置（ｘ、ｙ、ｚ）の近似値を得る方法。
【数４】
0005921892.tif 000014
【請求項５】
左足を載せるプレートＬと右足を載せるプレートＲを隣接させた重心計と、該重心計が検出したデータを基に人間モデルから身体重心（COG）の床面からの高さｌと軸（足関節）の粘性摩擦係数Ｃ＝［kgm2/S］ 及び慣性モーメントを参照し、立位状態を振り子モデルとして演算する手段と、体重をＷ，時刻ｔ の垂直力成分を Ｆｔ，時刻ｔ の前後力成分を Ｆy、左右力成分をＦｘとして、次の振り子の式からθの値を算出し、極座標（ｌ，θ）から身体重心に作用する仮想の力とCOGの空間の座標値ｚ、ｙを求める手段と、同様の手法でCOGの空間の座標値ｚ、ｘを求める手段とからなる重心計を用いて身体重心位置の近似値を求めるシステム。
【数５】
0005921892.tif 000015
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は重心計を用いて身体重心に作用する力と身体重心位置の近似値を求める技術に関する。
【背景技術】
【０００２】
身体重心とその点に作用する力（直交する３方向成分：Fz、Fx、Fy）を計測する方法としては、現在複数台のCCDカメラで三次元位置計測装置（以下モーションキャプチャという）と人間モデル（年齢、性別、身長、体重と各体節の重心位置及び質量と慣性モーメント等のデータベース）から演算する方法が広く用いられているが、据え置き型で場所の制約をうけること、校正に手間がかかること、非常に高価であることなどがあり、研究機関や大病院でのみ用いられているのが現状である。
そこで立位状態において重心計を使い簡便に身体重心位置とFz、Fx、Fyを計測可能にする方法を開発することを指向した。重心計は安価であり病院に普及しているので重心計でFz、Fx、Fy及び身体重心が計測できれば多くの患者に有益になると考えられるからである。
【０００３】
現在重心計で算出するプレート上の力の作用点（COP：center of pressure）は下記の方法で求められる。
図６の様に点Ｐの荷重Ｗにより、支点Ａ、ＢにＷＡ、ＷＢの荷重がかかる場合、点Ｐを中心とするモーメント（回転力）は、Ｐ、Ａ間の距離をｘとしＡ、Ｂ間の距離をＬとすると
ｘ・ＷＡ＝（Ｌ−ｘ）・ＷＢ  
の関係となるから、ｘは次式（１）で表せる。
ｘ＝Ｌ・ＷＢ／(ＷＡ＋ＷＢ)     ‥‥‥‥‥‥‥‥（１）
従って、点Ａ、Ｂでのそれぞれの垂直荷重を検出すれば、Ａ、Ｂ間の距離Ｌを用いて支点Ｐの位置を示すｘを求めることができる。
ただし、重心計の算出する点Ｐの位置は身体の重心（COG：Center Of Gravity）の位置ではなく、プレート上の力の作用点であるCOPの位置であり、このCOPと身体重心の位置は近似するが一致しない。
【０００４】
COPと身体重心COGとの関係については、トロント大学の政二慶氏らは、非特許文献１において「COPは身体重心の動揺と強く相関するものの、同一ではない。静止立位時の身体を１リンクの倒立振子で近似すると、倒立振子の運動方程式は
【数１】
0005921892.tif 000002
と表され、ここでＩ は倒立振子の慣性モーメント、θ は身体重心変位角、ｍ は身体質量、ｈ は重心高、ｇ は重力加速度、そしてＴ は足関節トルクである。身体重心COG の挙動を表すθ はトルクＴによって制御され、トルクＴ は神経中枢系により制御される筋収縮および足関節周りの粘弾性により供給される。ここで、COP 変位と足関節トルクの間にＴ≒ｍｇｕ が成り立つ（ｕは足関節を基点としたCOP変位）。」としている。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】政二  慶・他「高齢者の立体制御モデル」日本バイオメカニクス学会誌 Vol.11 No.4 2007年
【発明の概要】
【発明が解決しようとする課題】
【０００６】
上記の政二論文等から、COPは足関節トルクつまり足関節の力のモーメントの係数であり、COPの変位によって身体体重が左右されることを意味し、身体重心COGの変位はCOPの変位及び、足関節モーメントの変位に比例し、足関節モーメントは床反力の変位に相関することが明らかにされた。しかし、重心計のCOP値と身体重心COGの一義的な対応関係は明らかにされていない。
本発明の目的は重心計を用いて重心に作用する仮想の力Fx、Fyを求める方法と、本来は、三次元位置計測装置と人間モデル（年齢、性別、身長、体重と各体節の重心位置及び質量と慣性モーメント等のデータベース）を基に演算する身体重心の近似値を垂直荷重からCOPを計測する重心計を用いて求める方法を提供することにある。
【課題を解決するための手段】
【０００７】
本発明の重心計を用いて身体重心に作用する仮想の力（三次元成分：Ｆｘｔ，Ｆｙｔ，Ｆｚｔ）を求める方法は、垂直荷重からCOPを計測する重心計を用い、Δｔ間（微小時間）の立位状態における重心計上の力の作用点（COP）の変位量ΔCOPを測定するステップと、重心位置の移動速度（ＶＣＯＰ）をＶＣＯＰ＝ΔCOP／Δｔで、重心移動加速度（αＣＯＰ）をαＣＯＰ＝ΔＶＣＯＰ／Δｔで算出するステップと、ｔ時の重心に作用する仮想の力Ｆｘｔ，ＦｙｔとＦｚｔを次式で求めるステップとから成るものとした。
【数２】
0005921892.tif 000003
ただし、Ｇは重力加速度、Ｆｔは時刻ｔの荷重である。
また、左右の足をそれぞれに載せる二枚のプレートからなる重心計を用いて身体重心に作用する力を求める方法は、左足を載せるプレートＬ、右足を載せるプレートＲに対して先の方法によってそれぞれ仮想床反力を求めるステップと、左右のプレートＰ，Ｌ上の力の作用点（COP）ＰＬとＰＲを結ぶ線分を左右のプレートの荷重比で内分して合成COPを算出するステップと、前記求めた左右の仮想床反力を前記の線分に対する合成COPの位置比率でＸ，Ｙ成分ごとに算出するステップとからなる重心計を用いて身体重心に作用する力を求めるものとした。
【０００８】
本発明の重心計を用いて身体重心位置の近似値を求める方法は、人間モデルから身体重心（COG）の床面からの高さｌと軸（足関節）の粘性摩擦係数Ｃ＝［kgm2/S］及び慣性モーメントを参照し、立位状態を振り子モデルとして演算するステップと、体重を Ｗ，時刻ｔ の垂直力成分を Ｆｔ，時刻ｔ の前後力成分をＦy、左右力成分をＦｘとして、次の振り子の式からθの値を算出し、極座標（ｌ，θ）から身体重心に作用する仮想の力とCOGの空間の座標値ｚ、ｙを求めるステップと、同様の手法でCOGの空間の座標値ｚ、ｘを求めるステップと、からなるものとした。
【数３】
0005921892.tif 000004
【０００９】
重心計上でしゃがみ込んだ状態でのCOGの三次元空間位置（ｘ、ｙ、ｚ）の近似値を得る本発明の方法は、人間モデルから身体重心（COG）の床面からの高さｌと軸（足関節）の粘性摩擦係数Ｃ＝［kgm2/S］及び慣性モーメントを参照し、立位状態からしゃがみ込んだ状態を振り子モデルとして演算するステップと、体重を Ｗ，時刻ｔ のしゃがみ込み時にかかる垂直力成分を Ｆｔ’，時刻ｔの前後力成分を Ｆy、左右力成分をＦｘとして、次の振り子の式からθの値を算出し、極座標（ｌ，θ）から身体重心に作用する仮想の力とCOGの空間の座標値ｚ、ｙを求めるステップと、同様の手法でCOGの空間の座標値ｚ、ｘを求めるステップとからなるものとした。
【数４】
0005921892.tif 000005
【００１０】
本発明の重心計を用いて身体重心位置の近似値を求めるシステムは、左足を載せるプレートＬと右足を載せるプレートＲを隣接させた重心計と、該重心計が検出したデータを基に人間モデルから身体重心（COG）の床面からの高さｌと軸（足関節）の粘性摩擦係数Ｃ＝［kgm2/S］ 及び慣性モーメントを参照し、立位状態を振り子モデルとして演算する手段と、体重をＷ，時刻ｔ の垂直力成分を Ｆｔ，時刻ｔ の前後力成分を Ｆy、左右力成分をＦｘとして、次の振り子の式からθの値を算出し、極座標（ｌ，θ）から身体重心に作用する仮想の力とCOGの空間の座標値ｚ、ｙを求める手段と、同様の手法でCOGの空間の座標値ｚ、ｘを求める手段とからなる。
【数５】
0005921892.tif 000006
【発明の効果】
【００１１】
本発明の重心計を用いて身体重心に作用する力と身体重心位置の近似値を求める方法は、重心計と所定の演算処理を行う装置（ＰＣ）とから、モーションキャプチャを用いて求めた計測値と同様な三次元情報を得ることができるものである。したがって、モーションキャプチャを備えていない所でも容易に身体重心に作用する力と身体重心（COG）の近似値を得ることができる。
【００１２】
本発明の重心計を用いて身体重心に作用する力と身体重心位置の近似値を求めるシステムは、高価で大がかりなモーションキャプチャを用いることなく、単に重心計と所定の演算処理を行う装置（ＰＣ）とからなる簡素なシステムを備えるだけで、モーションキャプチャを用いて求めた計測値と同様なデータを得ることができるものである。
【図面の簡単な説明】
【００１３】
【図１】一枚の三角プレートの場合の床反力演算を説明する図である。
【図２】点ＰLと点Ｐ間の距離のＸ成分をRxm、Ｙ成分をRxm、点ＰRと点Ｐ間の距離のＸ成分をRxn、Ｙ成分をRxnとして求める方法を説明する図である。
【図３】立位状態を振り子モデルとして演算することを説明する図である。
【図４】本発明の方法を搭載したシステムのディスプレイ表示画面である。
【図５】従来の重心計システムでの可能な表示を比較例として示した図である。
【図６】バーの支点Ａ，Ｂ間に荷重Ｗがかかると、支点Ａ、ＢにＷＡ、ＷＢの分荷重がかかるることを説明する図である。
【図７】一枚の三角プレート上のCOP演算を説明する図である。
【図８】二枚の三角プレート上のCOP演算を説明する図である。
【発明を実施するための形態】
【００１４】
三角プレート１枚形態の重心計のCOP演算について考察する。図７に示すように正三角形のプレートの各頂点近傍には３つの荷重センサー（Ｐ１、Ｐ２、Ｐ３）が配置されている。プレートの各頂点から底辺に下ろした垂線の交点Ｃを原点としたプレート上の点Ｐ（ｘ，ｙ）に荷重Ｗの物体が存在するとき、３つの荷重センサー（Ｐ１、Ｐ２、Ｐ３）が、それぞれ荷重Ｗ１、Ｗ２、Ｗ３（Ｗ＝Ｗ１＋Ｗ２＋Ｗ３）を示した場合、原点Ｃを中心にＸ方向のモーメントとＹ方向のモーメントの釣り合い式は
Ｗ・ｘ＝ｍ（Ｗ１−Ｗ２）                 ‥‥‥‥‥‥‥‥（３）
Ｗ・ｙ＝ｎ２・Ｗ３−ｎ１（Ｗ１＋Ｗ２）   ‥‥‥‥‥‥‥‥（４）
となる。点Ｐの座標（ｘ，ｙ）は、次式（５）（６）によって求めることができる。すなわち、
ｘ＝ｍ（Ｗ１−Ｗ２）／Ｗ                  ‥‥‥‥‥‥‥‥（５）
ｙ＝｛ｎ２Ｗ３−ｎ１（Ｗ１−Ｗ２）｝／Ｗ  ‥‥‥‥‥‥‥‥（６）
ここで、ｍはＰ１、Ｐ２の中間点をＱとしたときＰ１、Ｑ間の長さ、ｎ１はＣ点からＱ点までの距離、ｎ２はＣ点からＰ３点までの距離を示している。
【００１５】
三角プレートを２枚使用した場合のCOP演算について考察する。上記の三角プレートを、互いの１辺を図８に示すように中央で突き合わせて配置し、中央の点Ｃを原点とする。左プレートのCOPであるＰＬと右プレートのCOPであるＰＲとの合成COPであるＰ（Ｘ，Ｙ）の座標値はＰＬとＰＲを結ぶ線分を各プレートの荷重比で内分し、次式（７）（８）のように表現できる。
Ｘ＝（ＸＬ・ＷＬ−ＸＲ・ＷＲ）／Ｗ          ‥‥‥‥‥‥‥‥（７）
Ｙ＝（ＹＬ・ＷＬ−ＹＲ・ＷＲ）／Ｗ          ‥‥‥‥‥‥‥‥（８）
そして荷重Ｗは左右の荷重の和（Ｗ＝ ＷＬ＋ＷＲ）となる。
また、プレートが四角形でも同様に荷重比でCOPを求められることは容易に理解されよう。
【００１６】
上記の要領で得られた重力計のCOP値を用いて、重心に作用する仮想の力Fx、Fyを求める計算方法を検討する。
まず、プレート一枚の場合（図１参照）をベクトルによる床反力演算手法によって解析する。現在の時刻ｔの重心位置をCOPｔ 、荷重をＦｔ とし、時刻ｔ−Δｔ の重心位置をCOPｔ−Δｔ、荷重をＦｔ−Δｔ とすると、重心位置の移動量ΔCOP は以下の式で表される。
ΔCOP＝COPｔ−COPｔ−Δｔ
＝（COPｔ・ｘ−COPｔ−Δｔ・ｘ，COPｔ・ｙ−COPｔ−Δｔ・ｙ）‥‥（９）
この時の重心位置移動速度ＶＣＯＰ と重心移動加速度αＣＯＰは、以下の式（10）（11）で表される。
ＶＣＯＰ＝ΔCOP／Δｔ         ‥‥‥‥‥‥‥‥（10）
αＣＯＰ＝ΔＶＣＯＰ／Δｔ    ‥‥‥‥‥‥‥‥（11）
これにより、質量ｍ の物体が、COPｔ−ΔｔからCOPｔまで移動するのに必要な力ｆＣＯＰは、以下の式（13）で表すことができる。
【数６】
0005921892.tif 000007
ここで、質量ｍは  ｍ＝Ｆｔ／Ｇなる関係にある。ただし、Ｇ は重力加速度である。
この力ｆＣＯＰをプレートにかかるｘ、ｙ 方向の力とすると、仮想床反力ｆ は以下の式（14）で表すことができる。
【数７】
0005921892.tif 000008
ここで、ΔＶＣＯＰ・x は重心位置移動速度ＶＣＯＰのＸ方向成分を、ΔＶＣＯＰ・yは重心位置移動速度ＶＣＯＰのＹ方向成分を表わしている。このｆをプレート１枚のときの仮想床反力と定義し、そのＸ，Ｙ，Ｚ方向成分をＦＸｔ，ＦＹｔ，ＦＺｔとすれば明示的に次式のように
【数８】
0005921892.tif 000009
で示されることになる。
【００１７】
プレート二枚の場合も一枚と同様に左プレートＬ、右プレートＲでそれぞれ仮想床反力を求め明示的に
ＦＬ＝（ＦLxt ×RatioL1，ＦLyt ×RatioL2，ＦLzt ×Ratio）
ＦＲ＝（ＦRxt ×RatioR1，ＦRyt ×RatioR2，ＦRzt ×Ratio）
と置く。
RatioL1，RatioL2，RatioR1，RatioR2は、図８の点ＰＬと点ＰＲで構成する直線を内分する点Ｐの比率で成分Ｘ，Ｙごとに求める。図２に示すように点ＰLと点Ｐ間の距離のＸ成分をRxm、Ｙ成分をRxm、点ＰRと点Ｐ間の距離のＸ成分をRxn、Ｙ成分をRxnとすると次のようになる。
RatioL1＝Rxn／（Rxm＋Rxn）
RatioL2＝Ryn／（Rym＋Ryn）
RatioR1＝Rxm／（Rxm＋Rxn）
RatioR2＝Rym／（Rym＋Ryn）
Ratio は計測対象が立位状態であることから、Ratio＝１とする。当然ながら立位状態ではない歩行などの動作の場合は１とはならない。
【００１８】
次に、上記の解析を基に身体重心の求め方を提示する。従来は、三次元位置計測装置（モーションキャプチャ）と人間モデル（年齢、性別、身長、体重と各体節の重心位置と質量及び慣性モーメント等のデータベース）から演算していた身体重心を近似値として以下の方法で求める。
人間モデルから身体重心（COG）の床面からの高さｌと軸（足関節）の粘性摩擦係数Ｃ＝［kgm2/S］ 及び慣性モーメント参照し、図３に示すように立位状態を振り子モデルとして演算する。（立位状態でしか使えない。）
【００１９】
振り子の式は
【数９】
0005921892.tif 000010
上記式でＧは重力加速度である。この上記式を解き、θの値を算出し、（ｌ，θ）からCOGの空間の座標値ｚ、ｙを求める。
Fxについても同様に解くことができCOGの三次元空間位置（ｘ、ｙ、ｚ）が得られる。
【００２０】
体重を Ｗ，時刻ｔ のしゃがみ込み時にかかる垂直力成分をＦｔ’，時刻ｔ の前後力成分を Ｆy  として
Ｆz＝Ｗ−Ｆｔ’      ‥‥‥‥‥‥‥‥（17）
この式（17）により立位状態からのしゃがみ込みにも対応し身体重心COGの座標位置が計算できる。
【数１０】
0005921892.tif 000011
上記式を解き、θの値を算出し、（ｌ，θ）からCOGの空間の座標値ｚ、ｙを求める。
Fxについても同様に解くことができしゃがみ込んだ状態でのCOGの三次元空間位置（ｘ、ｙ、ｚ）が得られる。
【実施例】
【００２１】
本発明の重心計を用いて身体重心に作用する力及び身体重心位置の近似値を求める方法を組込んだシステムによって、従来複数台のCCDカメラからなる三次元位置計測装置（モーションキャプチャ）と人間モデル（年齢、性別、身長、体重と各体節の重心位置及び質量と慣性モーメント等のデータベース）から演算する方法によって得られたデータと同様な鉛直方向（Ｚ方向）位置を含んだ三次元データが得られたことを示す。
【００２２】
図４は本発明の方法を搭載したシステムのディスプレイ表示画面である。重心計には三角プレートのものを二枚左右に隣接させ、それぞれに左右の足を乗せて立位状態で25秒間計測し、そのデータをディスプレイ上に表示させた。
左側のメインの表示はプレートの上方からの画面ではなく斜視図形態として身体重心位置を立体感覚で表示されるようにしている。プレート上には左右のCOPの軌跡と合成COPの軌跡及び左右のプレートが検出した身体重心変位角θの瞬時値がベクトル表示される。この左右のベクトル値を合成したものがその時点の身体重心となって、表示されている。画面右側上段のグラフは荷重のデータであり、左右のプレートの検出値が振幅の大きな振動波形として重なり、その上方に合成荷重が示されている。画面右側中段の表示は左右のプレートのCOPの軌跡と合成COPが水平面（ＸＹ平面）で表示されている。画面右側下段の表示は計測した基礎データの表表示である。荷重と、COPが総合値、左足、右足情報として表示されている。なお、荷重欄は被検者の体重であるが、直接表示することが個人情報として嫌がられるため、一般にある値を体重で割って％表示することが当業界では慣用されている。
【００２３】
従来の重心計システムでは鉛直方向（Ｚ方向）を含む三次元データを得ることはできないので、可能なシステム表示は図５に示すようになる。すなわち、左側のメインの表示は二次元データであるからプレートの上方からの画面となり、身体重心位置を立体感覚で表示させることはできない。プレート上には本発明と同様に左右のCOPの軌跡と合成COPの軌跡を表示できるが、合成COPは正確な位置ではない。そして、左右のプレートが検出する身体重心変位角θの瞬時値を得ることはできないのでベクトル表示させることは不可能である。画面右側上段のグラフと画面右側下段の検出データの表示は全く同じものとすることができるが、画面右側中段の表示では左右のプレートのCOPの軌跡は同じものが表示できるものの、合成COPは三次元位置としての身体重心を投影したものではないので、近似して見えても不正確な軌跡となっている。
【００２４】
以上のように、本発明の方法を実施したシステムでは、モーションキャプチャを備えていなくても、従来の重心計で検出できるデータを基に、身体重心の三次元位置を測定することが可能であることを実証できた。そして、重心計とパソコンの組み合わせからなる従来の同機種のものが測定できなかった身体重心位置とそこの作用する力をデータとして得ることができるという画期的なシステムとなっている。
【産業上の利用可能性】
【００２５】
本発明の方法を組込んだシステムは、モーションキャプチャを備えた大規模システムと同様な計測をハードとしては従来の重心計とパソコンの組み合わせからなる重心計測システムで行えることから、設備費用は格段に安価となるため、小さな病院や診療所、保健所、個人医院、学校の保健室、企業等の健康管理室など広い部署で設置可能となり、その普及が見込まれる。
【図１】
0005921892.tif 000016
【図２】
0005921892.tif 000017
【図３】
0005921892.tif 000018
【図６】
0005921892.tif 000019
【図７】
0005921892.tif 000020
【図８】
0005921892.tif 000021
【図４】
0005921892.tif 000022
【図５】
0005921892.tif 000023