明細書 :人体運動評価装置、方法、およびプログラム
| 発行国 | 日本国特許庁(JP) | ||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 公報種別 | 特許公報(B2) | ||||||||||||
| 特許番号 | 特許第5991532号 (P5991532) | ||||||||||||
| 公開番号 | 特開2014-113225 (P2014-113225A) | ||||||||||||
| 登録日 | 平成28年8月26日(2016.8.26) | ||||||||||||
| 発行日 | 平成28年9月14日(2016.9.14) | ||||||||||||
| 公開日 | 平成26年6月26日(2014.6.26) | ||||||||||||
| 発明の名称または考案の名称 | 人体運動評価装置、方法、およびプログラム | ||||||||||||
| 国際特許分類 | A61B 5/11 (2006.01) A61B 5/22 (2006.01) G06F 17/50 (2006.01) |
||||||||||||
| ＦＩ | A61B 5/10 310A A61B 5/22 B G06F 17/50 612L |
||||||||||||
| 請求項の数または発明の数 | 9 | ||||||||||||
| 全頁数 | 13 | ||||||||||||
| 出願番号 | 特願2012-267838 (P2012-267838) | ||||||||||||
| 出願日 | 平成24年12月7日(2012.12.7) | ||||||||||||
| 審査請求日 | 平成27年9月1日(2015.9.1) | ||||||||||||
| 特許権者または実用新案権者 | 【識別番号】504136568 【氏名又は名称】国立大学法人広島大学 |
||||||||||||
| 発明者または考案者 | 【氏名】栗田 雄一 【氏名】辻 敏夫 【氏名】櫻田 浩平 |
||||||||||||
| 個別代理人の代理人 | 【識別番号】110001427、【氏名又は名称】特許業務法人前田特許事務所 | ||||||||||||
| 審査官 | 【審査官】▲高▼ 芳徳 | ||||||||||||
| 参考文献・文献 | 特開２００５－６３３７５（ＪＰ，Ａ） 特開２０１０－２００７９７（ＪＰ，Ａ） 特表２００７－５３３０４３（ＪＰ，Ａ） 国際公開第２０１０／０９５６３６（ＷＯ，Ａ１） 特開２０１０－２９３４０（ＪＰ，Ａ） 特開２００６－１１０２１７（ＪＰ，Ａ） |
||||||||||||
| 調査した分野 | Ａ６１Ｂ ５／１１ － ５／１１３ Ａ６１Ｂ ５／２２ Ｇ０６Ｆ １７／５０ |
||||||||||||
| 特許請求の範囲 | 【請求項1】 人体の骨格、関節、および骨格筋をモデル化した筋骨格モデルを用いて人体運動を評価する装置であって、 前記筋骨格モデルの肢が到達運動をするときの前記肢の運動軌道を受けて前記筋骨格モデルにおける各関節の経時的な角度変化を算出する関節角度算出部と、 前記各関節の経時的な角度変化から各関節の経時的なトルク変化を算出する関節トルク算出部と、 前記各関節の経時的なトルク変化から前記筋骨格モデルにおける各骨格筋の経時的な筋力変化を算出する筋力算出部と、 前記各骨格筋の経時的な筋力変化から各骨格筋の経時的な筋活性度変化を算出する筋活性度算出部と、 前記各骨格筋の経時的な筋活性度変化に運動指令依存ノイズが重畳された状態で前記肢が前記到達運動をするときの到達運動終点のばらつきを評価する終点ばらつき評価部と、 前記肢の到達運動をさまざまに変えて、各到達運動に関して前記終点ばらつき評価部によって評価された到達運動終点のばらつきに基づいて、人間にとって位置決め精度を高くする到達運動を判定する到達運動判定部とを備えている ことを特徴とする人体運動評価装置。 【請求項2】 請求項１に記載の人体運動評価装置において、 前記終点ばらつき評価部は、 前記各骨格筋の経時的な筋活性度変化に前記運動指令依存ノイズを重畳して前記肢に複数回にわたって前記到達運動を試行させて各回の到達運動終点を算出し、 前記各回の到達運動終点と前記到達運動の基準終点との誤差を算出し、 前記誤差の平均値を算出することで前記到達運動終点のばらつきを評価する ことを特徴とする人体運動評価装置。 【請求項3】 請求項１に記載の人体運動評価装置において、 前記終点ばらつき評価部は、前記運動指令依存ノイズに関して定義された分散の時間積分値を前記筋骨格モデルにおける骨格筋ごとに算出し、これら時間積分値の総和を算出することで前記到達運動終点のばらつきを評価する ことを特徴とする人体運動評価装置。 【請求項4】 人体の骨格、関節、および骨格筋をモデル化した筋骨格モデルを用いて人体運動を評価する方法であって、 前記筋骨格モデルの肢が到達運動をするときの前記肢の運動軌道を受けて前記筋骨格モデルにおける各関節の経時的な角度変化を算出し、 前記各関節の角度変化から各関節の経時的なトルク変化を算出し、 前記各関節の経時的なトルク変化から前記筋骨格モデルにおける各骨格筋の経時的な筋力変化を算出し、 前記各骨格筋の経時的な筋力変化から各骨格筋の経時的な筋活性度変化を算出し、 前記各骨格筋の経時的な筋活性度変化に運動指令依存ノイズが重畳された状態で前記肢が前記到達運動をするときの到達運動終点のばらつきを評価し、 前記肢の到達運動をさまざまに変えて、各到達運動に関して前記評価した到達運動終点のばらつきに基づいて、人間にとって位置決め精度を高くする到達運動を判定する ことを特徴とする人体運動評価方法。 【請求項5】 請求項４に記載の人体運動評価方法において、 前記各骨格筋の経時的な筋活性度変化に前記運動指令依存ノイズを重畳して前記肢に複数回にわたって前記到達運動を試行させて各回の到達運動終点を算出し、 前記各回の到達運動終点と前記到達運動の基準終点との誤差を算出し、 前記誤差の平均値を算出することで前記到達運動終点のばらつきを評価する ことを特徴とする人体運動評価方法。 【請求項6】 請求項４に記載の人体運動評価方法において、 前記運動指令依存ノイズに関して定義された分散の時間積分値を前記筋骨格モデルにおける骨格筋ごとに算出し、これら時間積分値の総和を算出することで前記到達運動終点のばらつきを評価する ことを特徴とする人体運動評価方法。 【請求項7】 人体の骨格、関節、および骨格筋をモデル化した筋骨格モデルを用いてコンピュータに人体運動を評価させるプログラムであって、 前記筋骨格モデルの肢が到達運動をするときの前記肢の運動軌道を受けて前記筋骨格モデルにおける各関節の経時的な角度変化を算出する関節角度算出手段、 前記各関節の経時的な角度変化から各関節の経時的なトルク変化を算出する関節トルク算出手段、 前記各関節の経時的なトルク変化から前記筋骨格モデルにおける各骨格筋の経時的な筋力変化を算出する筋力算出手段、 前記各骨格筋の経時的な筋力変化から各骨格筋の経時的な筋活性度変化を算出する筋活性度算出手段、 前記各骨格筋の経時的な筋活性度変化に運動指令依存ノイズが重畳された状態で前記肢が前記到達運動をするときの到達運動終点のばらつきを評価する終点ばらつき評価手段、および 前記肢の到達運動をさまざまに変えて、各到達運動に関して前記終点ばらつき評価手段によって評価された到達運動終点のばらつきに基づいて、人間にとって位置決め精度を高くする到達運動を判定する到達運動判定手段として、コンピュータを機能させる ことを特徴とする人体運動評価プログラム。 【請求項8】 請求項７に記載の人体運動評価プログラムにおいて、 前記終点ばらつき評価手段は、 前記各骨格筋の経時的な筋活性度変化に前記運動指令依存ノイズを重畳して前記肢に複数回にわたって前記到達運動を試行させて各回の到達運動終点を算出し、 前記各回の到達運動終点と前記到達運動の基準終点との誤差を算出し、 前記誤差の平均値を算出することで前記到達運動終点のばらつきを評価する ことを特徴とする人体運動評価プログラム。 【請求項9】 請求項７に記載の人体運動評価プログラムにおいて、 前記終点ばらつき評価手段は、前記運動指令依存ノイズに関して定義された分散の時間積分値を前記筋骨格モデルにおける骨格筋ごとに算出し、これら時間積分値の総和を算出することで前記到達運動終点のばらつきを評価する ことを特徴とする人体運動評価プログラム。 |
||||||||||||
| 発明の詳細な説明 | 【技術分野】 【0001】 本発明は、人体の運動を評価する技術に関し、特に、筋骨格モデルを用いて人間にとって位置決め精度を高くする到達運動を判定する技術に関する。 【背景技術】 【0002】 人間の身体的特徴や生理学的特徴を工学的に利用しようとする人間工学が各種工業製品のヒューマンインタフェースの開発などに盛んに応用されている。人間工学に基づくデザインは人間にとって使いやすいだけではなく、人間が起こしそうなミスを未然に防いだりするのにも役立っている。 【0003】 人間工学では人体がさまざまなモデルで表され、コンピュータ上で人体のさまざまな動きがモデルによってシミュレートされる。モデルの例として、もっとも複雑なものとして有限要素モデルがあり、簡易なものとして人体の骨格、関節、および骨格筋をモデル化した筋骨格モデルがある。例えば、筋骨格モデルに基づく人体運動評価として、技能や感性を定量的に評価する評価指標を備え、人間の姿勢に近い姿勢を自動的に計算できる評価システムが、本願発明者らによって提案されている（例えば、特許文献１参照）。 【先行技術文献】 【0004】 【特許文献1】特開２０１１－１４１７０６号公報 【発明の概要】 【発明が解決しようとする課題】 【0005】 筋骨格モデルに基づく従来の人体運動評価では、筋骨格モデルにさまざまな到達運動をさせたときの各筋肉の筋活性度の時間積分値およびその総和が算出され、当該時間積分値の総和が小さいほどその到達運動の運動効率が高いと評価される。このような消費エネルギー最小規範に基づくと、人間がより疲れにくい動き、すなわち、エネルギー効率に優れた運動を判定することができるという利点がある。 【0006】 一方、脳が骨格筋に対して運動を指令するとき、運動強度に応じたノイズ（運動指令依存ノイズ）が当該指令に重畳されることが知られている。この知見に基づいて、人間は運動指令依存ノイズが最小となるような運動軌道を選択しているという終点分散最小規範が提唱されている。終点分散最小規範に基づくと、消費エネルギー最小規範では評価できなかったような、人間にとって自然で動きやすい運動の評価ができることが期待される。 【0007】 上記問題に鑑み、本発明は、人体の運動を位置決め精度の観点から評価できるようにすることを課題とする。 【課題を解決するための手段】 【0008】 本発明の一局面に従った人体運動評価装置は、人体の骨格、関節、および骨格筋をモデル化した筋骨格モデルを用いて人体運動を評価する装置であって、前記筋骨格モデルの肢が到達運動をするときの前記肢の運動軌道を受けて前記筋骨格モデルにおける各関節の経時的な角度変化を算出する関節角度算出部と、前記各関節の経時的な角度変化から各関節の経時的なトルク変化を算出する関節トルク算出部と、前記各関節の経時的なトルク変化から前記筋骨格モデルにおける各骨格筋の経時的な筋力変化を算出する筋力算出部と、前記各骨格筋の経時的な筋力変化から各骨格筋の経時的な筋活性度変化を算出する筋活性度算出部と、前記各骨格筋の経時的な筋活性度変化に運動指令依存ノイズが重畳された状態で前記肢が前記到達運動をするときの到達運動終点のばらつきを評価する終点ばらつき評価部と、前記肢の到達運動をさまざまに変えて、各到達運動に関して前記終点ばらつき評価部によって評価された到達運動終点のばらつきに基づいて、人間にとって位置決め精度を高くする到達運動を判定する到達運動判定部とを備えている。 【0009】 また、本発明の一局面に従った人体運動評価方法は、人体の骨格、関節、および骨格筋をモデル化した筋骨格モデルを用いて人体運動を評価する方法であって、前記筋骨格モデルの肢が到達運動をするときの前記肢の運動軌道を受けて前記筋骨格モデルにおける各関節の経時的な角度変化を算出し、前記各関節の角度変化から各関節の経時的なトルク変化を算出し、前記各関節の経時的なトルク変化から前記筋骨格モデルにおける各骨格筋の経時的な筋力変化を算出し、前記各骨格筋の経時的な筋力変化から各骨格筋の経時的な筋活性度変化を算出し、前記各骨格筋の経時的な筋活性度変化に運動指令依存ノイズが重畳された状態で前記肢が前記到達運動をするときの到達運動終点のばらつきを評価し、前記肢の到達運動をさまざまに変えて、各到達運動に関して前記評価した到達運動終点のばらつきに基づいて、人間にとって位置決め精度を高くする到達運動を判定する。 【0010】 また、本発明の一局面に従った人体運動評価プログラムは、人体の骨格、関節、および骨格筋をモデル化した筋骨格モデルを用いてコンピュータに人体運動を評価させるプログラムであって、前記筋骨格モデルの肢が到達運動をするときの前記肢の運動軌道を受けて前記筋骨格モデルにおける各関節の経時的な角度変化を算出する関節角度算出手段、前記各関節の経時的な角度変化から各関節の経時的なトルク変化を算出する関節トルク算出手段、前記各関節の経時的なトルク変化から前記筋骨格モデルにおける各骨格筋の経時的な筋力変化を算出する筋力算出手段、前記各骨格筋の経時的な筋力変化から各骨格筋の経時的な筋活性度変化を算出する筋活性度算出手段、前記各骨格筋の経時的な筋活性度変化に運動指令依存ノイズが重畳された状態で前記肢が前記到達運動をするときの到達運動終点のばらつきを評価する終点ばらつき評価手段、および前記肢の到達運動をさまざまに変えて、各到達運動に関して前記終点ばらつき評価手段によって評価された到達運動終点のばらつきに基づいて、人間にとって位置決め精度を高くする到達運動を判定する到達運動判定手段として、コンピュータを機能させる。 【0011】 これら人体運動評価装置、方法、およびプログラムによると、筋骨格モデルに基づく肢のさまざまな到達運動に関して到達運動終点のばらつきが評価され、当該評価された到達運動終点のばらつきに基づいて、人間にとって位置決め精度を高くする到達運動が判定される。なお、肢とは、上腕、下肢、手先、指先、足先などのことを指す。これにより、人体の運動を位置決め精度の観点から評価することができる。 【0012】 上記人体運動評価装置、方法、およびプログラムのそれぞれにおいて、前記各骨格筋の経時的な筋活性度変化に前記運動指令依存ノイズを重畳して前記肢に複数回にわたって前記到達運動を試行させて各回の到達運動終点を算出し、前記各回の到達運動終点と前記到達運動の基準終点との誤差を算出し、前記誤差の平均値を算出することで前記到達運動終点のばらつきを評価することができる。 【0013】 あるいは、上記人体運動評価装置、方法、およびプログラムのそれぞれにおいて、前記運動指令依存ノイズに関して定義された分散の時間積分値を前記筋骨格モデルにおける骨格筋ごとに算出し、これら時間積分値の総和を算出することで前記到達運動終点のばらつきを評価することができる。 【発明の効果】 【0014】 本発明によると、人体の運動を位置決め精度の観点から評価することができる。 【図面の簡単な説明】 【0015】 【図1】本発明の一実施形態に係る人体運動評価装置の主要部の機能ブロック図 【図2】一例に係る筋骨格モデル（２関節６筋上肢モデル）の模式図 【図3】図２に示した２関節６筋上肢モデルにある到達運動をさせたときの手先の運動軌道および各骨格筋の筋力変化を示す図 【図4】運動指令依存ノイズが重畳されたときのある到達運動（人間にとって自然な運動軌道）の終点のばらつきを示す図 【図5】運動指令依存ノイズが重畳されたときの別の到達運動（人間にとって不自然な運動軌道）の終点のばらつきを示す図 【図6】到達運動判定部による判定対象の到達運動を示す模式図 【図7】図６に示した各到達運動の終点のばらつきを示す図 【図8】図７に示した各到達運動の終点の誤差平均値を示す箱ひげ図 【図9】人の到達運動の実験概要を示す模式図 【図10】図９に示した各到達運動の終点のばらつきを示す図 【図11】図１０に示した各到達運動の終点の誤差平均値を示す箱ひげ図 【発明を実施するための形態】 【0016】 以下、図面を参照しながら本発明を実施するための形態について説明する。なお、本発明は、以下の実施形態に限定されるものではない。 【0017】 図１は、本実施形態に係る人体運動評価装置１０の主要部の機能ブロック図である。人体運動評価装置１０は、関節角度算出部１１、関節トルク算出部１２、筋力算出部１３、筋活性度算出部１４、終点ばらつき評価部１５、および到達運動判定部１６を備えている。なお、人体運動評価装置１０は、上記各構成要素を半導体集積回路などで構成した専用ハードウェアとして実施することができる。あるいは、上記各構成要素をコンピュータプログラムで記述して、ＰＣなどの汎用コンピュータに当該コンピュータプログラムを実行させることで人体運動評価装置１０を汎用コンピュータ上で実施することもできる。さらに、人体運動評価装置１０はハードウェアとソフトウェアとを組み合わせて実施することもできる。 【0018】 人体運動評価装置１０は、筋骨格モデルを用いて人間にとって位置決め精度を高くする到達運動を判定する。なお、説明の便宜のため、本実施形態では筋骨格モデルとして図２に示すような２関節６筋上肢モデルを用いるが、本発明で利用可能な筋骨格モデルはこれに限定されない。 【0019】 関節角度算出部１１は、２関節６筋上肢モデル（筋骨格モデル）の手先（肢）が到達運動をするときの手先の運動軌道を受けて、２関節６筋上肢モデルにおける各関節の経時的な角度変化を算出する。当該２関節６筋上肢モデルにおいて手先は平面上を運動するため、関節角度算出部１１に手先の到達運動として手先点１０１の座標位置（Ｘ，Ｙ）が入力される。そして、関節角度算出部１１は、入力された座標位置（Ｘ，Ｙ）から、図２に示した関節Ｊ１，Ｊ２の角度ｑ１，ｑ２を算出する。関節角度ｑ１，ｑ２は、手先点１０１の座標位置（Ｘ，Ｙ）の変化に応じて経時的に変化するパラメータである。 【0020】 なお、手先点１０１の座標位置（Ｘ，Ｙ）は、人体運動評価装置１０に数値入力してもよいし、人体に装着されたモーションキャプチャ（図示せず）などを人体運動評価装置１０に接続して当該モーションキャプチャなどから取得することもできる。 【0021】 関節トルク算出部１２は、関節Ｊ１，Ｊ２の経時的な角度変化（関節角度ｑ１，ｑ２）から関節Ｊ１，Ｊ２の経時的なトルク変化を算出する。関節トルク算出部１２は、例えば、手先点１０１が外力Ｆの運動をするときの運動方程式に基づいて、図２に示した関節Ｊ１，Ｊ２のトルクｔ１，ｔ２を算出する。関節トルクｔ１，ｔ２は、手先点１０１の座標位置（Ｘ，Ｙ）の変化に応じて経時的に変化するパラメータである。 【0022】 筋力算出部１３は、関節Ｊ１，Ｊ２の経時的なトルク変化（関節トルクｔ１，ｔ２）から筋骨格モデルにおける各骨格筋の経時的な筋力変化を算出する。図２に示したように２関節６筋上肢モデルには６つの筋力Ｆ１～Ｆ６があり、筋力算出部１３は、関節トルクｔ１，ｔ２を生み出すような各骨格筋の運動として筋力Ｆ１～Ｆ６を算出する。筋力Ｆ１～Ｆ６は、手先点１０１の座標位置（Ｘ，Ｙ）の変化に応じて経時的に変化するパラメータである。 【0023】 図３は、図２に示した２関節６筋上肢モデルにある到達運動をさせたときの手先の運動軌道および各骨格筋の筋力変化を示す。図３（ａ）は、２関節６筋上肢モデルをｘｙ平面上に投影した図である。関節Ｊ１はｘｙ平面の原点に固定されており、手先点１０１は運動軌道１０２に沿って始点位置（０，０．２）から終点位置（０．３，０．２）まで到達運動をする。図３（ｂ）は、２関節６筋上肢モデルに上記到達運動をさせたときの、筋力算出部１３が算出した筋力Ｆ１～Ｆ６をプロットしたグラフである。 【0024】 図１に戻り、筋活性度算出部１４は、各骨格筋の経時的な筋力変化（筋力Ｆ１～Ｆ６）から各骨格筋の経時的な筋活性度変化を算出する。筋活性度とは、骨格筋を構成する運動単位が活性している割合のことである。したがって、筋活性度は０から１までの値をとる。筋活性度と筋力との間には、例えば式（１）で表されるＨａｔｚｅの筋収縮モデルがある。筋活性度算出部１４は、例えば、式（１）に基づいて各骨格筋の筋活性度α１～α６を算出する。筋活性度α１～α６は、手先点１０１の座標位置（Ｘ，Ｙ）の変化に応じて経時的に変化するパラメータである。 Ｆｉ＝Ｔｍａｘ＿ｉ・（αｉ・ｋ（ｘ）・ｈ（ｘ’）＋ｆｐ（ｘ）） … （１） ただし、Ｔｍａｘ＿ｉはｉ番目の骨格筋の最大筋力、すなわち、筋力Ｆｉの最大値、ｋ（ｘ），ｈ（ｘ’），ｆｐ（ｘ）は筋の長さｘと速度ｘ’によって決まる筋特性である。 【0025】 ここで筋力には、その出力が大きくなるほど、出力の変動の度合いが大きくなる、という性質があることが知られている（信号強度依存ノイズ）。このような知見に基づいて、Ｈａｒｒｉｓらは運動指令にその大きさに比例したノイズ（運動指令依存ノイズ）が重畳されると仮定し、運動軌道は運動指令依存ノイズが重畳された到達運動の終点の揺らぎ（分散）が最小になるように計画されるという終点分散最小規範を提唱した。Ｈａｒｒｉｓらによる運動指令依存ノイズは式（２）のように定式化される。 σｗ２＝ｋｕ２ … （２） ただし、ｗは運動指令に重畳するノイズ、σｗ２はｗの分散、ｕは運動指令、ｋはノイズ係数である。 【0026】 終点ばらつき評価部１５は、各骨格筋の経時的な筋活性度変化（筋活性度α１～α６）に運動指令依存ノイズ１００が重畳された状態で手先１０１が到達運動をするときの到達運動終点のばらつきを評価する。運動指令依存ノイズ１００は、式（２）における運動指令に重畳するノイズｗに該当する。筋活性度は運動指令とみなすことができるため、式（２）より、筋活性度αｉに重畳される運動指令依存ノイズｗｉの分散σｗｉ２は式（３）のように表される。 σｗｉ２＝ｋαｉ２ … （３） すなわち、筋活性度αｉに重畳される運動指令依存ノイズは、式（３）を満たす確率変数をランダムに発生させることで擬似することができる。このとき、運動指令異存ノイズが付加された筋活性度αｉ’は式（４）で表される。 αｉ’＝αｉ＋ｗｉ … （４） 【0027】 終点ばらつき評価部１５は、筋活性度αｉについて式（３）を満たすランダムな確率変数、すなわち、運動指令依存ノイズ１００を生成する。そして、終点ばらつき評価部１５は、筋活性度αｉに運動指令依存ノイズ１００を重畳して手先１０１に複数回（数１０～数１００回程度）、到達運動を試行させて各回の到達運動終点を算出し、各回の到達運動終点と到達運動の基準終点との誤差を算出し、さらに誤差の平均値を算出することで到達運動終点のばらつきを評価する。誤差は、到達運動の基準終点からの到達運動終点のずれ、すなわち、到達運動終点と到達運動の基準終点との間の距離である。 【0028】 図４は、運動指令依存ノイズが重畳されたときのある到達運動の終点のばらつきを示す。当該到達運動の運動軌道は緩やかな弧を描いており、人間にとって自然な運動軌道であると考えられる。一方、図５は、運動指令依存ノイズが重畳された別の到達運動の終点のばらつきを示す。当該到達運動の運動軌道は蛇行しており、人間にとって不自然な運動軌道であると考えられる。図４（ａ）および図５（ａ）は、２関節６筋上肢モデルをｘｙ平面上に投影した図である。関節Ｊ１はｘｙ平面の原点に固定されており、手先点１０１は運動軌道１０２に沿って始点位置（－０．２７，０．３）から終点位置（０．２７，０．３５）まで到達運動をする。当該終点位置が上記の到達運動の基準終点に該当する。図４（ｂ）および図５（ｂ）は、２関節６筋上肢モデルを用いて上記到達運動を複数回試行させたときの手先１０１の移動軌跡をプロットした図である。便宜のため、２関節６筋上肢モデルに指示した運動軌道を白の破線で示している。図４および図５から、人間にとって不自然と考えられる運動軌道で到達運動が行われると到達運動終点が大きくばらつくことがわかる。図４および図５の例で到達運動終点の誤差平均値を算出した場合、図４の例よりも図５の例の方が誤差平均値が大きく算出される。 【0029】 あるいは、終点ばらつき評価部１５は、到達運動を試行することに代えて、例えば、式（３）で定義された運動指令依存ノイズｗｉの分散σｗｉ２に対して所定の演算処理を行うことで到達運動終点のばらつきを評価することもできる。具体的には、終点ばらつき評価部１５は、筋活性度レベルでの運動指令依存ノイズの分散の時間積分値を骨格筋ごとに算出し、これら時間積分値の総和を算出することで到達運動終点のばらつきを評価することができる（式（５）参照）。このような演算処理による到達運動終点のばらつき評価はあくまでも簡易な評価方法であるため、上記のような到達運動の試行によって求めたばらつきとは必ずしも一致しない可能性があるが、到達運動の試行による到達運動終点のばらつき評価の代替となり得る。 Σｉ＝１Ｎ∫０Ｔσｗｉ２ｄｔ（＝Σｉ＝１Ｎ∫０Ｔｋαｉ２ｄｔ） … （５） ただし、Ｎは骨格筋の総数であり、本例ではＮ＝６、Ｔは到達運動の終了時刻である。 【0030】 また、終点ばらつき評価部１５は、筋力について運動指令依存ノイズの分散の時間積分値を骨格筋ごとに算出し、これら時間積分値の総和を算出することで到達運動終点のばらつきを評価することができる（式（６）参照）。 Σｉ＝１Ｎ∫０ＴσＦｉ２ｄｔ … （６） ただし、Ｎは骨格筋の総数であり、本例ではＮ＝６、Ｔは到達運動の終了時刻、σＦｉ２は式（１）に基づいて計算される筋力Ｆｉの分散である。 【0031】 なお、計算式の提示は省略するが、終点ばらつき評価部１５は、関節トルク、関節角度、手先位置のそれぞれについて運動指令依存ノイズの分散の時間積分値を骨格筋ごとに算出し、これら時間積分値の総和を算出することで到達運動終点のばらつきを評価することもできる。 【0032】 図１に戻り、到達運動判定部１６は、手先の到達運動をさまざまに変えて、各到達運動に関して終点ばらつき評価部１５によって評価された到達運動終点のばらつきに基づいて、人間にとって位置決め精度を高くする到達運動を判定する。図６は、到達運動判定部１６による判定対象の到達運動を模式的に示す。ここでは評価対象の到達運動が２つあり、一方の到達運動の運動軌道１０２ｉｎは、２関節６筋上肢モデルの手先が身体の内側に向けて動くときのものであり、手先点１０１は始点位置（０，０．３）から終点位置（－０．２，０．３）まで移動する。他方の到達運動の運動軌道１０２ｏｕｔは、同じく２関節６筋上肢モデルの手先が身体の外側に向けて動くときのものであり、手先点１０１は始点位置（０，０．３）から終点位置（０．２，０．３）まで移動する。すなわち、手先点１０１の移動距離は運動軌道１０２ｉｎおよび運動軌道１０２ｏｕｔで同じである。なお、関節Ｊ１は上記と同様にｘｙ平面の原点に固定されている。 【0033】 図６に示したように、筋肉の使用率は運動軌道１０２ｉｎおよび運動軌道１０２ｏｕｔでほぼ同じである。したがって、従来の筋発揮効率性の観点から運動軌道１０２ｉｎおよび運動軌道１０２ｏｕｔを評価した場合、両者に差はないと判定される。 【0034】 一方、図７は、図６に示した各到達運動について、終点ばらつき評価部１５によって評価された到達運動終点のばらつきを示す。図７（ａ）は、運動軌道１０２ｉｎの到達運動終点のばらつきを示し、図７（ｂ）は、運動軌道１０２ｏｕｔの到達運動終点のばらつきを示す。図７（ａ）と図７（ｂ）を比較してわかるように、２関節６筋上肢モデルの手先を身体の内側に移動させるよりも外側に移動させる方が到達運動終点のばらつきが小さくなる傾向にある。 【0035】 図８は、終点ばらつき評価部１５によって算出された、図７に示した各到達運動の終点の誤差平均値を箱ひげ図で示したものである。運動軌道１０２ｉｎの到達運動終点の誤差平均値は約０．００８であるのに対して、運動軌道１０２ｏｕｔの到達運動終点の誤差平均値は約０．００７である。Ｐ値は０．００１未満である。したがって、到達運動判定部１６は、この数値評価の結果から、運動軌道１０２ｉｎの到達運動および運動軌道１０２ｏｕｔの到達運動のうち、後者の方が人間にとって位置決め精度を高くする到達運動であると判定する。 【0036】 次に、人体運動評価装置１０の評価結果の妥当性を検証するために実際に人に到達運動をさせたときの実験結果について説明する。図９は、実験概要を模式的に示す。被験者３名（２２±１歳）はいずれもアイマスクを着用して、右手の手先を右肩前方の位置から左側、すなわち、身体内側に移動させる運動軌道１０２ｉｎ、および右手の手先を右肩前方の位置から右側、すなわち、身体外側に移動させる運動軌道１０２ｏｕｔで１０回ずつ到達運動を行った。運動軌道１０２ｉｎは、右肩の位置を原点に固定して、右手の手先点１０１が始点位置（０，０．２）から終点位置（－０．２，０．２）まで移動する運動軌道である。運動軌道１０２ｏｕｔは、右肩の位置を原点に固定して、右手の手先点１０１が始点位置（０，０．２）から終点位置（０．３，０．２）まで移動する運動軌道である。なお、いずれの運動軌道も運動時間は０．５秒であり、実験前に２０回の練習運動を行った。手先点１０１の計測にはＳＰＩＣＥ社製の光学式トラッキングシステムＯｐｔｉＴｒａｃｋ Ｖ１２０を使用した。 【0037】 図１０は、図９に示した各到達運動の到達運動終点のばらつきを示す。図１０（ａ）は、運動軌道１０２ｉｎの到達運動終点のばらつきを示し、図１０（ｂ）は、運動軌道１０２ｏｕｔの到達運動終点のばらつきを示す。 【0038】 図１１は、図１０に示した各到達運動の終点の誤差平均値を箱ひげ図で示したものである。運動軌道１０２ｉｎの到達運動終点の誤差平均値は約０．０４であるのに対して、運動軌道１０２ｏｕｔの到達運動終点の誤差平均値は約０．０２である。Ｐ値は０．００１未満である。したがって、この数値評価の結果から、運動軌道１０２ｉｎの到達運動および運動軌道１０２ｏｕｔの到達運動のうち、後者の方が実際に人間にとって位置決め精度を高くする到達運動であることが示される。この実験結果は人体運動評価装置１０の評価結果と合致しており、人体運動評価装置１０による評価結果が妥当であることを実証するものである。 【0039】 以上のように本実施形態によると、人体の運動を位置決め精度の観点から評価することができる。これにより、各種工業製品のヒューマンインタフェースにおける各種部品（例えば、自動車のパワーウィンドウのボタンなど）の配置に関して、人間が位置決めしやすい位置を見つけることができる。 【産業上の利用可能性】 【0040】 本発明に係る人体運動評価技術は、人体の運動を位置決め精度の観点から評価することができるため、各種工業製品のヒューマンインタフェースの開発用途などに有用である。 【符号の説明】 【0041】 １０ 人体運動評価装置 １１ 関節角度算出部 １２ 関節トルク算出部 １３ 筋力算出部 １４ 筋活性度算出部 １５ 終点ばらつき評価部 １６ 到達運動判定部 １００ 運動指令依存ノイズ |
||||||||||||
| 図面 |
|