支持面揺動時の姿勢制御に対する視覚の影響

システム工学群
知能メカトロダイナミクス研究室

2.

周波数応答実験

2.1 周波数応答実験概要
本研究では，応答の再現性確保のために支持面の水平揺動に対する人の受動的な応答の計測を目的とし床反力計を搭載した揺動装置を用いて周波数応答実験を行った．人体の応答は，反射マーカを頭，肩，腸骨稜，大転子，膝，くるぶし，支持面に左右対称に計 14 個設置し，モーションキャプチャで計測した．実験条件として被験者は 20 代男性 7 名，実験回数は 12 回とし，スタンス幅は 230mm とし，視線はまっすぐ遠くを見るように指示した．人が支持面の動きを予測できないように 3 パターンの揺動波形を用意した．図 1 のように人体の足首を回転軸とした 1 自由度の剛体とみなし，マーカの位置座標からその姿勢角と支持面の加速度の関係から周波数応答線図を導出した．
2.2 周波数応答線図
支持面の加速度を入力とし，人の姿勢角を出力とした周波数応答線図を図 2 に示す．被験者 1 名分の実験回数 12 回の平均値を白丸，標準偏差をエラーバーとして記す．
図 2 を見ると，1rad/s 付近で位相が 90°程度遅れていることがわかる．この特性は被験者 G を除く全被験者に共通してみられており，先行研究でも同じような傾向が確認されている．低周波帯域では視覚による影響が現れると考えられており，位相特性が視覚の影響によって励起されたと仮定したモデルを検討する．
2.3 モデルの構築
本研究では，図 1 のように人体を 1 自由度の剛体リンクと仮定してモデルの定義を行う．人のバランス制御則については現在確立されていないが一般的には随意制御として遅延を含む状態フィードバックと扱われることが多い(1)．しかし，この考え方では低周波帯域の位相特性が再現できない．そこ

京平で，本研究では一般的な状態フィードバック制御に低周波位相特性を再現するための制御則を付加することを検討する．
モデルの検討にあたり，ベースとなる図 2 の 1 自由度剛体リンクモデルの運動方程式は次のようになる．
(𝐽𝑏 + 𝑚𝑏 𝑙𝑏2 )𝜃̈𝑏 − 𝑚𝑏 𝑔𝑙𝑏 = 𝑚𝑏 𝑙𝑏 𝑋̈ + 𝜏𝑏
(1)
θb は鉛直軸からの角度（左回りが正），Xs は支持面の絶対変位，lb は回転軸から質量中心までの長さ，mb は身体部の質量
Jb は身体部の質量中心まわりの慣性モーメントを表す．式（1）の右辺第 1 項は支持面の加速度によって生じる慣性力の影響を表し，第 2 項は人体の姿勢制御トルクを表す．一般的な随意制御の制御則として，以下のような遅延を含む状態フィードバック制御とおく．
𝜏𝑏 = −𝑘𝑝 𝜃𝑏 (𝑡 − 𝑡𝑑 ) − 𝑘𝑑 𝜃̇𝑏 (𝑡 − 𝑡𝑑 )(2)
k p ,k d は姿勢角と角速度に応じたフィードバックの係数を表し，t d は応答の遅延時間を表す．
上記の仮定したシステムでは低周波帯域の位相特性を再現することができないため，本研究では，上記のシステムに支持面の速度フィードバックを加えたモデルを検討する．このモデルは以下のように仮定する．
𝜏𝑏 = −𝑘𝑝 {𝜃𝑏 (𝑡 − 𝑡𝑑 ) − 𝜃̅𝑏 (𝑡)} − 𝑘𝑑 𝜃̇𝑏 (𝑡 − 𝑡𝑑 )
(3)
1𝛼̇̅
̅
̇
𝜃𝑏 + 𝜃𝑏 = 𝑋𝑠𝑐𝑐

式(3)は支持面の速度に起因して生じる目標姿勢角度を加えたものである．𝜃𝑏̅ は 1 次遅れ系とみなして，cはその時定数，αは比例定数を表す．このモデルに対して支持面の加速度を入力，姿勢角を出力として伝達関数を求めると以下のようになる．
𝐺(𝑠)𝑚𝑏 𝑙𝑏 𝑐𝑠 2 + 𝑚𝑏 𝑙𝑏 𝑠 + 𝛼𝑘𝑝=𝑠(𝑐𝑠 + 1){(𝐽𝑏 + 𝑚𝑏 𝑙𝑏2 )𝑠 2 + 𝑘𝑑 𝑒 −𝑡𝑑 𝑠 𝑠 + 𝑘𝑝 𝑒 −𝑡𝑑 𝑠 − 𝑚𝑏 𝑔𝑙𝑏 }(4)
式(4)のように定義した伝達関数が実験結果に一致するようにフィッティングを行い，未知の制御パラメータであるk p ,k d ,t d ,c,αの 5 つを同定した．この推定したシステムと実験結果の比較を図 2 に示す．このグラフから，仮定したモデルと実験結果はほぼ一致していることがわかり，支持面速度フィードバックを加えたモデルは安定なシステムであると言える．

1

Experimental result
Estimated system

Gain

緒言
人体の立位バランス能力は視覚，前庭感覚，体性感覚の組み合わせのフィードフォワード制御，あるいはフィードバック制御によってなされているとされている(1)．この制御則を明らかにすることができれば，個人のバランス能力の推定やリハビリ患者の回復具合の推定などに応用できることが期待される．先行研究によって周波数応答実験を行った際の人のバランス制御は低周波帯域で位相が遅れることがわかっている．我々は視覚による床面速度のフィードバックによってこの遅れが生じているのだと考えた．この仮説を検証するために，本研究では視覚による影響を明らかにすることを目的としている．
人体に視覚だけの情報を与える方法として VR（virtualreality）を用いる．実際に支持面を揺らす実験の被験者目線の動画を VR で見せることで視覚の情報だけを被験者に与える．視覚による影響が顕著に出やすいのは単一周波数だと考え，VR 実験は単一周波数で行った．同じ被験者で比較を行うため，床面揺動による周波数応答実験の結果と，単一周波数 VR 実験の結果を示す．

00Phase [deg]1.

西村

1180124−90−180−2700246Frequency [rad/s]

Fig. 1 Single rigid bodymodelFig. 2 Frequency response
diagram and system
identification

0.4

0.2

0−0.2

Fig. 4Experimental diagram−0.4

Table 1 Velocity condition of support surface in VR test
ExperimentalⅠⅡⅢⅣVideo
500
333250
Amplitude[mm] 1000
0.10.20.3
0.4
Frequency[Hz]10101010Cycle

3.2 結果モーションキャプチャによって得られた位置座標から姿勢角を導出した結果を図 5 に示す．実験動画Ⅰの 5 周期分を示す．また，姿勢角の後半 6 周期分(動画Ⅰは 3 周期分)をフーリエ変換し，導出した周波数応答線図を図 6 に示す．被験者ごとに 4 種類の動画で場合分けし，6 回実験した平均値を出力している．
白丸は揺動周波数，黒線が 0.1Hz，緑線が 0.2Hz，水色線が 0.3Hz，赤線が 0.4Hz を示している．
すべての被験者で動画の周波数に対して上段のゲインでピークが出ていることから，視覚刺激によって姿勢を傾けていることが確認された．また，全体的な傾向として低周波帯域の方がゲインは大きくなることがわかった．しかし，被験者 G に関しては全く逆の傾向を示している．

400.004

Sway frequency0.1 Hz
0.2 Hz0.3 Hz0.4 Hz

0.002

Gain

Gain
Phase (rad)単一周波数 VR 実験

0220−20.2

0.4

Frequency (Hz)Sway frequency
0.1 Hz0.2 Hz0.3 Hz0.4 Hz0.002

000−200.2

0.4

Frequency (Hz)(a) Subject A(c) Subject C

0.006

Gain

3.1 単一周波数 VR 実験概要
本実験は支持面を前額面方向に単一周波数で揺動した時の，立位バランスの視覚による影響を見ることが目的である．
また，揺動周波数の違いによって応答に差が生じるのかを確認する．
被験者に見てもらう動画は図 3 に示す撮影機器を使い撮影した．電動スケートボードに角柱を取り付け，Vuze cameraを取り付けた三脚を角柱に固定した． Vuze camera は 8 つのレンズが付いており，360°全体を撮影することが可能である．
撮影した映像は専用ソフト Humaneyes VR Studio で VR 動画に変換し本実験で使用した．実験に用いる動画は 0.1，0.2，0.3，0.4Hz の角周波数で揺動する 4 種類あり，動画内容を表1 に示す．電動スケートボードの支持面速度𝑋̇を
𝑋̇ = 𝐴𝜔𝑐𝑜𝑠𝜔𝑡(5)としたときの強度Aωを一定に設定した．Aは支持面の振幅，ωは角速度，tは時間である．
本実験の様子を図 4 に示す．人体の姿勢変化を計測するためにモーションキャプチャと床反力計を使用した．反射マー
カの位置，被験者，スタンス幅，は周波数応答実験と同様である．実験回数は被験者一人につき，4 種類の動画を 6 回ず
つ，計 24 回行った．視線は一点を凝視するのではなくまわりを見るようにと指示した．

20

Fig. 5 Posture angle in VR test

0.004

Sway frequency
0.1 Hz0.2 Hz
0.3 Hz0.4 Hz0.004
0.002
0Phase (rad)3.

0Phase (rad)Fig. 3 PhotographDevice20−200.2

0.4Frequency (Hz)(g) Subject G

Fig. 6 Frequency response diagram in VR test
下段の位相については被験者によってばらばらである．仮定したモデルと一致しているのは被験者 B，C，E であるが，他の被験者の位相はあまり変化しない．これは，単一周波数であるため被験者が揺動を予測することが容易であり，この予測によって結果が変わってきていると考えられる．
4.

結言
本研究では，立位人体の前額面バランスに関して支持面水平揺動による周波数応答実験を行った．実験結果で見られた低周波帯域の位相遅れを説明できるモデルを仮定し，安定なシステムが得られた．仮定したモデルは支持面の速度によって姿勢角が変化すると仮定しており，視覚情報が低周波帯域においては大きく働くと考えられる．
そこで，視覚だけの情報を与えて応答が現れるのか確認するため単一周波数 VR 実験を行った．実験結果から，明らかに揺動周波数で人体の応答が出ていることがわかり，視覚情報が立位バランスに関与していることが確認された．
文献(1) Robertson,D.G.E.,Whittlesey,S.N.,Caldwell,G.E.,Kamen,G.
andHamill,J., 身体運動のバイオメカニクス研究法，大修館書店(2008),p.66