JP4805818B2 - 移動ロボットの歩容生成装置 - Google Patents
移動ロボットの歩容生成装置 Download PDFInfo
- Publication number
- JP4805818B2 JP4805818B2 JP2006510511A JP2006510511A JP4805818B2 JP 4805818 B2 JP4805818 B2 JP 4805818B2 JP 2006510511 A JP2006510511 A JP 2006510511A JP 2006510511 A JP2006510511 A JP 2006510511A JP 4805818 B2 JP4805818 B2 JP 4805818B2
- Authority
- JP
- Japan
- Prior art keywords
- reaction force
- floor reaction
- motion
- error
- force component
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 230000005021 gait Effects 0.000 title claims description 756
- 238000006243 chemical reaction Methods 0.000 claims description 1407
- 230000033001 locomotion Effects 0.000 claims description 1091
- 238000012937 correction Methods 0.000 claims description 334
- 238000000034 method Methods 0.000 claims description 281
- 230000008569 process Effects 0.000 claims description 265
- 238000012545 processing Methods 0.000 claims description 188
- 230000008859 change Effects 0.000 claims description 148
- 238000011156 evaluation Methods 0.000 claims description 46
- 238000013519 translation Methods 0.000 claims description 46
- 238000004364 calculation method Methods 0.000 claims description 40
- 230000009471 action Effects 0.000 claims description 16
- 230000003252 repetitive effect Effects 0.000 claims description 5
- 238000012804 iterative process Methods 0.000 claims description 4
- 210000002683 foot Anatomy 0.000 description 101
- 210000002414 leg Anatomy 0.000 description 98
- 230000001133 acceleration Effects 0.000 description 73
- 230000005484 gravity Effects 0.000 description 34
- 238000006073 displacement reaction Methods 0.000 description 19
- 210000001503 joint Anatomy 0.000 description 17
- 239000002131 composite material Substances 0.000 description 12
- 230000007246 mechanism Effects 0.000 description 12
- 230000014509 gene expression Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 210000000544 articulatio talocruralis Anatomy 0.000 description 6
- 238000001514 detection method Methods 0.000 description 6
- 230000000737 periodic effect Effects 0.000 description 6
- 210000000629 knee joint Anatomy 0.000 description 5
- 239000000470 constituent Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 210000002310 elbow joint Anatomy 0.000 description 4
- 210000003857 wrist joint Anatomy 0.000 description 4
- 210000003423 ankle Anatomy 0.000 description 3
- 210000004394 hip joint Anatomy 0.000 description 3
- 210000000323 shoulder joint Anatomy 0.000 description 3
- 230000006641 stabilisation Effects 0.000 description 3
- 238000011105 stabilization Methods 0.000 description 3
- 206010048669 Terminal state Diseases 0.000 description 2
- 238000012854 evaluation process Methods 0.000 description 2
- 210000003128 head Anatomy 0.000 description 2
- 210000001624 hip Anatomy 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000002945 steepest descent method Methods 0.000 description 2
- 210000000689 upper leg Anatomy 0.000 description 2
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 239000013013 elastic material Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 210000003127 knee Anatomy 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 210000001699 lower leg Anatomy 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 230000003094 perturbing effect Effects 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J13/00—Controls for manipulators
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B62—LAND VEHICLES FOR TRAVELLING OTHERWISE THAN ON RAILS
- B62D—MOTOR VEHICLES; TRAILERS
- B62D57/00—Vehicles characterised by having other propulsion or other ground- engaging means than wheels or endless track, alone or in addition to wheels or endless track
- B62D57/02—Vehicles characterised by having other propulsion or other ground- engaging means than wheels or endless track, alone or in addition to wheels or endless track with ground-engaging propulsion means, e.g. walking members
- B62D57/032—Vehicles characterised by having other propulsion or other ground- engaging means than wheels or endless track, alone or in addition to wheels or endless track with ground-engaging propulsion means, e.g. walking members with alternately or sequentially lifted supporting base and legs; with alternately or sequentially lifted feet or skid
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J5/00—Manipulators mounted on wheels or on carriages
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Transportation (AREA)
- Robotics (AREA)
- Manipulator (AREA)
Description
前記移動ロボットの所定の第1動力学モデルを用いて該移動ロボットの運動を作成するようにあらかじめ構成された所定のアルゴリズムを実行することにより、前記床反力成分が前記許容範囲を満たすように、前記目標運動の仮値である仮運動を作成する仮目標運動作成手段と、
前記第1動力学モデルと該第1動力学モデルよりも動力学的精度が高い所定の第2動力学モデルとを用い、前記移動ロボットの任意の運動に対して、該運動によって前記第2動力学モデル上で発生する前記床反力成分と該運動によって前記第1動力学モデル上で発生する前記床反力成分との差である床反力成分誤差を求める床反力成分誤差算出手段と、
前記作成された仮運動から前記床反力成分誤差算出手段によって求められる床反力成分誤差が所定の許容誤差範囲に収まっているか否かを評価する評価手段と、
前記評価手段の評価において前記仮運動に対応する床反力成分誤差が前記許容誤差範囲に収まっている場合には、前記仮運動を前記目標運動として決定すると共に、前記仮運動に対応する床反力成分誤差が前記許容誤差範囲から逸脱している場合には、前記仮運動を少なくとも1回以上修正して前記目標運動を決定する目標運動決定手段とを備え、
前記仮運動に対応する床反力成分誤差が前記許容誤差範囲から逸脱している場合において、前記目標運動決定手段が前記仮運動の第n回目(n:n≧1である整数)の修正によって求める修正後運動を第n修正後運動とし、該第n修正後運動から前記床反力成分誤差算出手段により求められる床反力成分誤差を第n床反力成分誤差とし、前記仮運動を第0修正後運動とし、前記仮運動に対応する床反力成分誤差を第0床反力成分誤差としたとき、
前記目標運動決定手段は、第n修正後運動によって前記第1動力学モデル上で発生する前記床反力成分に、第n−1床反力成分誤差と少なくとも該第n−1床反力成分誤差に応じて決定した床反力補正量とのいずれかを加えたものが前記許容範囲を満たすように前記所定のアルゴリズムを実行することにより、第n修正後運動を決定する修正後運動決定処理と、前記決定した第n修正後運動に対応する第n床反力成分誤差と第n−1修正後運動に対応する第n−1床反力成分誤差との差、および、前記第n床反力成分誤差と少なくとも前記第n−1床反力成分誤差に応じて決定した前記床反力補正量との差のいずれか一方として定義される床反力成分誤差変化量ΔFMが0に収束したか否かを判断する収束判断処理とを少なくとも該床反力成分誤差変化量ΔFMが0に収束したと判断されるまで繰り返し、この繰り返し処理における最後の修正後運動決定処理で決定した修正後運動を前記目標運動として決定することを特徴とするものである。
前記目標運動の仮値である仮運動を作成する仮目標運動作成手段と、
前記移動ロボットの所定の第1の動力学モデルと該第1動力学モデルよりも動力学的精度が高い所定の第2動力学モデルとを用い、前記移動ロボットの任意の運動に対して、該運動によって前記第2動力学モデル上で発生する前記床反力成分と該運動によって前記第1動力学モデル上で発生する前記床反力成分との差である床反力成分誤差を求める床反力成分誤差算出手段と、
前記仮運動を少なくとも1回以上修正して前記目標運動を決定する目標運動決定手段とを備え、
前記目標運動決定手段が前記仮運動の第n回目(n:n≧1である整数)の修正によって求める修正後運動を第n修正後運動とし、該第n修正後運動から前記床反力成分誤差算出手段により求められる床反力成分誤差を第n床反力成分誤差とし、前記仮運動を第0修正後運動とし、前記仮運動から前記床反力成分誤差算出手段により求められる床反力成分誤差を第0床反力成分誤差としたとき、
前記目標運動決定手段は、前記第1動力学モデルを用いて該移動ロボットの運動を作成するようにあらかじめ構成された所定のアルゴリズムを実行することにより、第n修正後運動によって前記第1動力学モデル上で発生する前記床反力成分に、第n−1床反力成分誤差と少なくとも該第n−1床反力成分誤差に応じて決定した床反力補正量とのいずれかを加えたものが前記許容範囲を満たすように第n修正後運動を決定する修正後運動決定処理と、前記決定した第n修正後運動に対応する第n床反力成分誤差と第n−1修正後運動に対応する第n−1床反力成分誤差との差、および、前記第n床反力成分誤差と少なくとも前記第n−1床反力成分誤差に応じて決定した前記床反力補正量との差のいずれか一方として定義される床反力成分誤差変化量ΔFMが0に収束したか否かを判断する収束判断処理とを少なくとも該床反力成分誤差変化量ΔFMが0に収束したと判断されるまで繰り返し、この繰り返し処理における最後の修正後運動決定処理で決定した修正後運動を前記目標運動として決定することを特徴とするものである。
前記移動ロボットの所定の第1動力学モデルを用いて該移動ロボットの運動を作成するようにあらかじめ構成された所定のアルゴリズムを実行することにより、前記第1床反力成分が前記目標値に一致し、且つ、前記第2床反力成分が前記許容範囲を満たすように、前記目標運動の仮値である仮運動を作成する仮目標運動作成手段と、
前記第1動力学モデルと該第1動力学モデルよりも動力学的精度が高い所定の第2動力学モデルとを用い、前記移動ロボットの任意の運動に対して、該運動によって前記第2動力学モデル上で発生する前記第1床反力成分と該運動によって前記第1動力学モデル上で発生する前記第1床反力成分との差である床反力成分誤差Aerrを求めると共に、該運動によって前記第2動力学モデル上で発生する前記第2床反力成分と該運動によって前記第1動力学モデル上で発生する前記第2床反力成分との差である床反力成分誤差Berrを求める床反力成分誤差算出手段と、
前記作成された仮運動から前記床反力成分誤差算出手段によって求められる床反力成分誤差Aerr,Berrのうちの床反力成分誤差Aerrが所定の第1許容誤差範囲に収まっているか否かを評価すると共に、床反力成分誤差Berrが所定の第2許容誤差範囲に収まっているか否かを評価する評価手段と、
前記評価手段の評価において前記仮運動に対応する床反力成分誤差Aerr,Berrの両者がそれぞれに対応する前記第1許容誤差範囲および第2許容誤差範囲に収まっている場合には、前記仮運動を前記目標運動として決定すると共に、前記仮運動に対応する床反力成分誤差Aerr,Berrの少なくともいずれか一方がそれに対応する前記第1許容誤差範囲または第2許容誤差範囲から逸脱している場合には、前記仮運動を少なくとも1回以上修正して前記目標運動を決定する目標運動決定手段とを備え、
前記仮運動に対応する床反力成分誤差Aerr,Berrの少なくともいずれか一方がそれに対応する前記第1許容誤差範囲または第2許容誤差範囲から逸脱している場合において、前記目標運動決定手段が前記仮運動の第n回目(n:n≧1である整数)の修正によって求める修正後運動を第n修正後運動とし、該第n修正後運動から前記床反力成分誤差算出手段により求められる床反力成分誤差Aerr,Berrをそれぞれ第n床反力成分誤差Aerr(n),Berr(n)とし、前記仮運動を第0修正後運動とし、前記仮運動に対応する床反力成分誤差Aerr,Berrをそれぞれ第0床反力成分誤差Aerr(0),Berr(0)としたとき、
前記目標運動決定手段は、第n修正後運動によって前記第1動力学モデル上で発生する前記第1床反力成分に、第n−1床反力成分誤差Aerr(n-1)と少なくとも該第n−1床反力成分誤差Aerr(n-1)に応じて決定した第1床反力補正量とのいずれかを加えたものが前記目標値に一致し、且つ、該第n修正後運動によって前記第1動力学モデル上で発生する前記第2床反力成分に、第n−1床反力成分誤差Berr(n-1)と少なくとも該第n−1床反力成分誤差Berr(n-1)に応じて決定した第2床反力補正量とのいずれかを加えたものが前記許容範囲を満たすように前記所定のアルゴリズムを実行することにより、第n修正後運動を決定する修正後運動決定処理と、前記決定した第n修正後運動に対応する第n床反力成分誤差Aerr(n)と第n−1修正後運動に対応する第n−1床反力成分誤差Aerr(n-1)との差、および、前記第n床反力成分誤差Aerr(n)と少なくとも前記第n−1床反力成分誤差Aerr(n-1)に応じて決定した前記第1床反力補正量との差のいずれか一方として定義される床反力成分誤差変化量ΔAerrと、前記決定した第n修正後運動に対応する第n床反力成分誤差Berr(n)と第n−1修正後運動に対応する第n−1床反力成分誤差Berr(n-1)との差、および、前記第n床反力成分誤差Berr(n)と少なくとも前記第n−1床反力成分誤差Berr(n-1)に応じて決定した前記第2床反力補正量との差のいずれか一方として定義される床反力成分誤差変化量ΔBerrとがそれぞれ0に収束したか否かを判断する収束判断処理とを少なくとも該床反力成分誤差変化量ΔAerr,ΔBerrの両者が0に収束したと判断されるまで繰り返し、この繰り返し処理における最後の前記修正後運動決定処理で決定した修正後運動を前記目標運動として決定することを特徴とするものである。
前記目標運動の仮値である仮運動を作成する仮目標運動作成手段と、
前記移動ロボットの所定の第1動力学モデルと該第1動力学モデルよりも動力学的精度が高い所定の第2動力学モデルとを用い、前記移動ロボットの任意の運動に対して、該運動によって前記第2動力学モデル上で発生する前記第1床反力成分と該運動によって前記第1動力学モデル上で発生する前記第1床反力成分との差である床反力成分誤差Aerrを求めると共に、該運動によって前記第2動力学モデル上で発生する前記第2床反力成分と該運動によって前記第1動力学モデル上で発生する前記第2床反力成分との差である床反力成分誤差Berrを求める床反力成分誤差算出手段と、
前記仮運動を少なくとも1回以上修正して前記目標運動を決定する目標運動決定手段とを備え、
前記目標運動決定手段が前記仮運動の第n回目(n:n≧1である整数)の修正によって求める修正後運動を第n修正後運動とし、該第n修正後運動から前記床反力成分誤差算出手段により求められる床反力成分誤差Aerr,Berrをそれぞれ第n床反力成分誤差Aerr(n),Berr(n)とし、前記仮運動を第0修正後運動とし、前記仮運動から前記床反力成分誤差算出手段により求められる床反力成分誤差Aerr,Berrをそれぞれ第0床反力成分誤差Aerr(0),Berr(0)としたとき、
前記目標運動決定手段は、前記第1動力学モデルを用いて該移動ロボットの運動を作成するようにあらかじめ構成された所定のアルゴリズムを実行することにより、第n修正後運動によって前記第1動力学モデル上で発生する前記第1床反力成分に、第n−1床反力成分誤差Aerr(n-1)と少なくとも該第n−1床反力成分誤差Aerr(n-1)に応じて決定した第1床反力補正量とのいずれかを加えたものが前記目標値に一致し、且つ、該第n修正後運動によって前記第1動力学モデル上で発生する前記第2床反力成分に、第n−1床反力成分誤差Berr(n-1)と少なくとも該第n−1床反力成分誤差Berr(n-1)に応じて決定した第2床反力補正量とのいずれかを加えたものが前記許容範囲を満たすように第n修正後運動を決定する修正後運動決定処理と、前記決定した第n修正後運動に対応する第n床反力成分誤差Aerr(n)と第n−1修正後運動に対応する第n−1床反力成分誤差Aerr(n-1)との差、および、前記第n床反力成分誤差Aerr(n)と少なくとも前記第n−1床反力成分誤差Aerr(n-1)に応じて決定した前記第1床反力補正量との差のいずれか一方として定義される床反力成分誤差変化量ΔAerrと、前記決定した第n修正後運動に対応する第n床反力成分誤差Berr(n)と第n−1修正後運動に対応する第n−1床反力成分誤差Berr(n-1)との差、および、前記第n床反力成分誤差Berr(n)と少なくとも前記第n−1床反力成分誤差Berr(n-1)に応じて決定した前記第2床反力補正量との差のいずれか一方として定義される床反力成分誤差変化量ΔBerrとがそれぞれ0に収束したか否かを判断する収束判断処理とを少なくとも該床反力成分誤差変化量ΔAerr,ΔBerrの両者が0に収束したと判断されるまで繰り返し、この繰り返し処理における最後の前記修正後運動決定処理で決定した修正後運動を前記目標運動として決定することを特徴とするものである。
前記移動ロボットの所定の第1動力学モデルを用いて該移動ロボットの運動を作成するようにあらかじめ構成された所定のアルゴリズムを実行することにより、前記目標ZMPを満足し、且つ、前記並進床反力水平成分が前記許容範囲を満たすように、前記目標運動の仮値である仮運動を作成する仮目標運動作成手段と、
前記第1動力学モデルと該第1動力学モデルよりも動力学的精度が高い所定の第2動力学モデルとを用い、前記移動ロボットの任意の運動に対して、該運動によって前記第2動力学モデル上で算出されるZMPと該運動によって前記第1動力学モデル上で算出されるZMPとの差であるZMP誤差ZMPerrを求めると共に、該運動によって前記第2動力学モデル上で発生する前記並進床反力水平成分と該運動によって前記第1動力学モデル上で発生する前記並進床反力水平成分との差である並進床反力水平成分誤差Ferrを求める誤差算出手段と、
前記作成された仮運動から前記誤差算出手段によって求められるZMP誤差ZMPerrおよび並進床反力水平成分誤差FerrのうちのZMP誤差ZMPerrが所定の第1許容誤差範囲に収まっているか否かを評価すると共に、並進床反力水平成分誤差Ferrが所定の第2許容誤差範囲に収まっているか否かを評価する評価手段と、
前記評価手段の評価において前記仮運動に対応するZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferrの両者がそれぞれに対応する前記第1許容誤差範囲および第2許容誤差範囲に収まっている場合には、前記仮運動を前記目標運動として決定すると共に、前記仮運動に対応するZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferrの少なくともいずれか一方がそれに対応する前記第1許容誤差範囲または第2許容誤差範囲から逸脱している場合には、前記仮運動を少なくとも1回以上修正して前記目標運動を決定する目標運動決定手段とを備え、
前記仮運動に対応するZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferrの少なくともいずれか一方がそれに対応する前記第1許容誤差範囲または第2許容誤差範囲から逸脱している場合において、前記目標運動決定手段が前記仮運動の第n回目(n:n≧1である整数)の修正によって求める修正後運動を第n修正後運動とし、該第n修正後運動から前記床反力成分誤差算出手段により求められるZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferrをそれぞれ第nZMP誤差ZMPerr(n)および第n並進床反力水平成分誤差Ferr(n)とし、前記仮運動を第0修正後運動とし、前記仮運動に対応するZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferrをそれぞれ第0ZMP誤差ZMPerr(0)および第0並進床反力水平成分誤差Ferr(0)としたとき、
前記目標運動決定手段は、第n修正後運動によって前記第1動力学モデル上で算出されるZMPに、第n−1ZMP誤差ZMPerr(n-1)と少なくとも該第n−1ZMP誤差ZMPerr(n-1)に応じて決定したZMP補正量とのいずれかを加えたものが前記目標ZMPに一致し、且つ、該第n修正後運動によって前記第1動力学モデル上で発生する前記並進床反力水平成分に、第n−1並進床反力水平成分誤差Ferr(n-1)と少なくとも該第n−1並進床反力水平成分誤差Ferr(n-1)に応じて決定した床反力補正量とのいずれかを加えたものが前記許容範囲を満たすように前記所定のアルゴリズムを実行することにより、第n修正後運動を決定する修正後運動決定処理と、前記決定した第n修正後運動に対応する第nZMP誤差ZMPerr(n)と第n−1修正後運動に対応する第n−1ZMP誤差ZMPerr(n-1)との差、および、前記第nZMP誤差ZMPerr(n)と少なくとも前記第n−1ZMP誤差ZMPerr(n-1)に応じて決定した前記ZMP補正量との差のいずれか一方として定義されるZMP誤差変化量ΔZMPerrと、前記決定した第n修正後運動に対応する第n並進床反力水平成分誤差Ferr(n)と第n−1修正後運動に対応する第n−1並進床反力水平成分誤差Ferr(n-1)との差、および、第n並進床反力水平成分誤差Ferr(n)と少なくとも前記第n−1並進床反力水平成分誤差Ferr(n-1)に応じて決定した前記床反力補正量との差のいずれか一方として定義される並進床反力水平成分誤差変化量ΔFerrとがそれぞれ0に収束したか否かを判断する収束判断処理とを少なくとも該ZMP誤差変化量ΔZMPerrおよび並進床反力水平成分誤差変化量ΔFerrの両者が0に収束したと判断されるまで繰り返し、この繰り返し処理における最後の前記修正後運動決定処理で決定した修正後運動を前記目標運動として決定することを特徴とするものである。
前記目標運動の仮値である仮運動を作成する仮目標運動作成手段と、
前記移動ロボットの所定の第1動力学モデルと該第1動力学モデルよりも動力学的精度が高い所定の第2動力学モデルとを用い、前記移動ロボットの任意の運動に対して、該運動によって前記第2動力学モデル上で算出されるZMPと該運動によって前記第1動力学モデル上で算出されるZMPとの差であるZMP誤差ZMPerrを求めると共に、該運動によって前記第2動力学モデル上で発生する前記並進床反力水平成分と該運動によって前記第1動力学モデル上で発生する前記並進床反力水平成分との差である並進床反力水平成分誤差Ferrを求める誤差算出手段と、
前記仮運動を少なくとも1回以上修正して前記目標運動を決定する目標運動決定手段とを備え、
前記目標運動決定手段が前記仮運動の第n回目(n:n≧1である整数)の修正によって求める修正後運動を第n修正後運動とし、該第n修正後運動から前記床反力成分誤差算出手段により求められるZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferrをそれぞれ第nZMP誤差ZMPerr(n)および第n並進床反力水平成分誤差Ferr(n)とし、前記仮運動を第0修正後運動とし、前記仮運動から前記床反力成分誤差算出手段により求められるZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferrをそれぞれ第0ZMP誤差ZMPerr(0)および第0並進床反力水平成分誤差Ferr(0)としたとき、
前記目標運動決定手段は、前記第1動力学モデルを用いて該移動ロボットの運動を作成するようにあらかじめ構成された所定のアルゴリズムを実行することにより、第n修正後運動によって前記第1動力学モデル上で算出されるZMPに、第n−1ZMP誤差ZMPerr(n-1)と少なくとも該第n−1ZMP誤差ZMPerr(n-1)に応じて決定したZMP補正量とのいずれかを加えたものが前記目標ZMPに一致し、且つ、該第n修正後運動によって前記第1動力学モデル上で発生する前記並進床反力水平成分に、第n−1並進床反力水平成分誤差Ferr(n-1)と少なくとも該第n−1並進床反力水平成分誤差Ferr(n-1)に応じて決定した床反力補正量とのいずれかを加えたものが前記許容範囲を満たすように第n修正後運動を決定する修正後運動決定処理と、前記決定した第n修正後運動に対応する第nZMP誤差ZMPerr(n)と第n−1修正後運動に対応する第n−1ZMP誤差ZMPerr(n-1)との差、および、前記第nZMP誤差ZMPerr(n)と少なくとも前記第n−1ZMP誤差ZMPerr(n-1)に応じて決定した前記ZMP補正量との差のいずれか一方として定義されるZMP誤差変化量ΔZMPerrと、前記決定した第n修正後運動に対応する第n並進床反力水平成分誤差Ferr(n)と第n−1修正後運動に対応する第n−1並進床反力水平成分誤差Ferr(n-1)との差、および、第n並進床反力水平成分誤差Ferr(n)と少なくとも前記第n−1並進床反力水平成分誤差Ferr(n-1)に応じて決定した前記床反力補正量との差のいずれか一方として定義される並進床反力水平成分誤差変化量ΔFerrとがそれぞれ0に収束したか否かを判断する収束判断処理とを少なくとも該ZMP誤差変化量ΔZMPerrおよび並進床反力水平成分誤差変化量ΔFerrの両者が0に収束したと判断されるまで繰り返し、この繰り返し処理における最後の前記修正後運動決定処理で決定した修正後運動を前記目標運動として決定することを特徴とするものである。
前記床反力成分誤差算出手段により求められる前記第m床反力成分誤差Aerr,Berr(m:m≧0である整数)のうちの第m床反力成分誤差Aerrは、第m修正後運動によって該第m修正後運動の各時刻に前記第2動力学モデル上で発生する前記第1床反力成分の瞬時値と該第m修正後運動によって該時刻に前記第1動力学モデル上で発生する前記第1床反力成分の瞬時値との差の前記所定期間における時系列から構成されると共に、前記第m床反力成分誤差Berrは、第m修正後運動の各時刻に前記第2動力学モデル上で発生する前記第2床反力成分の瞬時値と該第m修正後運動によって該時刻に前記第1動力学モデル上で発生する前記第2床反力成分の瞬時値との差の前記所定期間における時系列から構成され、前記評価手段の前記所定の第1許容誤差範囲は、前記第0床反力成分誤差Aerr(0)を構成する時系列パターンのうちの所定の第1特徴量に対する許容誤差範囲であると共に前記所定の第2許容誤差範囲は、前記第0床反力成分誤差Berr(0)を構成する時系列のパターンのうちの所定の第2特徴量に対する許容誤差範囲であり、前記目標運動決定手段の前記修正後運動決定処理は、前記第n修正後運動の各時刻tにおいて、該第n修正後運動によって時刻tに前記第1動力学モデル上で発生する前記第1床反力成分の瞬時値に、前記第n−1床反力成分誤差Aerr(n-1)の時刻tでの値と少なくとも該値に応じて決定した前記第1床反力補正量の値とのいずれかを加えたものが該時刻tでの前記目標値を満足し、且つ、該第n修正後運動によって該時刻tに前記第1動力学モデル上で発生する前記第2床反力成分の瞬時値に、前記第n−1床反力成分誤差Berr(n-1)の時刻tでの値と少なくとも該値に応じて決定した前記第2床反力補正量の値とのいずれかを加えたものが該時刻tでの前記許容範囲を満たすように該時刻tにおける第n修正後運動の瞬時値を決定する処理であり、前記床反力成分誤差変化量ΔAerrは、前記第n床反力成分誤差Aerr(n)の各時刻での値と第n−1床反力成分誤差Aerr(n-1)の該時刻での値との差、および、前記第n床反力成分誤差Aerr(n)の各時刻での値と少なくとも第n−1床反力成分誤差Aerr(n-1)の該時刻での値に応じて決定した前記第1床反力補正量の値との差のいずれか一方の前記所定期間における時系列から構成されると共に、前記床反力成分誤差変化量ΔBerrは、前記第n床反力成分誤差Berr(n)の各時刻での値と第n−1床反力成分誤差Berr(n-1)の該時刻での値との差、および、前記第n床反力成分誤差Berr(n)の各時刻での値と少なくとも第n−1床反力成分誤差Berr(n-1)の該時刻での値に応じて決定した前記第2床反力補正量の値との差のいずれか一方の前記所定期間における時系列から構成され、
前記目標運動決定手段の前記収束判断処理は、前記床反力成分誤差変化量ΔAerrを構成する時系列のパターンのうちの所定の第3特徴量が該第3特徴量に対する所定の許容変化量範囲に収まったときに該床反力成分誤差変化量ΔAerrが0に収束したと判断すると共に、前記床反力成分誤差変化量ΔBerrを構成する時系列のパターンのうちの所定の第4特徴量が該第4特徴量に対する所定の許容変化量範囲に収まったときに該床反力成分誤差変化量ΔBerrが0に収束したと判断する処理であることが好ましい(第15発明)。
D1)単純化モデルを用いて生成した歩容(単純化モデル歩容)を変位次元補正用モデルを用いて修正してなる歩容(以降、変位次元補正歩容と呼ぶ)よりも高い精度で動力学的平衡条件を満足する。
D2)真のZMP(目標ZMPのまわりに補正目標床反力モーメントを発生させることによって修正された本来の定義を満足するZMP)は、ZMP許容範囲(安定余裕が十分維持できる許容範囲)に存在する。
D3)床反力水平成分は床反力水平成分許容範囲内になる。
という条件を満足するように、
E1)単純化モデル歩容の上体位置姿勢を補正する。
あるいは
E2)単純化モデル歩容の上体位置姿勢を補正すると共に、目標ZMPまわりの補正目標床反力モーメントを出力する(目標床反力を補正する)。
Zsup:支持脚質点鉛直位置、Zswg:遊脚質点鉛直位置、Zb:上体質点鉛直位置、ZGtotal:全体重心鉛直位置、Xsup:支持脚質点水平位置、Xswg:遊脚質点水平位置、Xb:上体質点水平位置、XGtotal:全体重心水平位置、θby:鉛直方向に対するY軸回りの上体姿勢角(傾斜角)、mb:上体質点質量、msup:支持脚質点質量、mswg:遊脚質点質量、mtotal:ロボット総質量(=mb+msup+mswg)、J:上体慣性モーメント(上体回転モードにおける等価慣性モーメント)、Fx:床反力水平成分(詳しくは並進床反力の前後方向(X軸)成分)、Fz:床反力鉛直成分(詳しくは並進床反力の鉛直方向(Z軸)成分)、My:目標ZMPまわりの床反力モーメント(詳しくは床反力モーメントの左右軸(Y軸)まわり成分)、g:重力加速度。なお、任意の変数Xに対してd2X/dt2は、Xの2階微分値を意味するものとする。
Fz=mb*(g+d2Zb/dt2)+msup*(g+d2Zsup/dt2)
+mswg*(g+d2Zswg/dt2) ……式1a
Fx=mb*d2Xb/dt2+msup*d2Xsup/dt2+mswg*d2Xswg/dt2
……式1b
My=−mb*(Xb−Xzmp)+(g+d2Zb/dt2)+mb*(Zb−Zzmp)*d2Xb/dt2
−msup*(Xsup−Xzmp)*(g+d2Zsup/dt2)
+msup*(Zsup−Zzmp)*d2Xsup/dt2
−mswg*(Xswg−Xzmp)*(g+d2Zswg/dt2)
+mswg*(Zswg−Zzmp)*(d2Xswg/dt2)+J*d2θby/dt2
……式1c
なお、ロボット全体重心位置には次の関係式が成立する。
ZGtotal=(mb*Zb+msup*Zsup+mswg*Zswg)/mtotal …式1d
XGtotal=(mb*Xb+msup*Xsup+mswg*Xswg)/mtotal …式1e
この単純化モデルでは、脚体2,2の動力学(各質点2m,2mの動力学)と上体3の動力学(質点24m及びフライホイールFHの動力学)とが相互に非干渉に構成されると共に、ロボット1全体の動力学は、それらの線形結合で表される。また、さらに、上体3の運動と床反力との関係は、上体3の並進運動(上体並進モード)と床反力との関係、並びに上体3の回転運動(上体回転モード)と床反力との関係に分離される。具体的には、上体質点24mの水平運動によって発生する床反力は、上体3の水平方向並進運動(上体並進モード)によって発生する床反力に対応し、フライホイールの回転運動によって発生する床反力は、上体3の回転運動(上体回転モード)、すなわち上体3の姿勢変化運動によって発生する床反力に対応する。
発散成分=上体質点水平位置+上体質点水平速度/ω0 …式2
この式2のω0は所定の値である。なお、本実施形態では、歩容における発散成分は、式2の上体質点水平位置、上体質点水平速度の代わりに、それぞれ上体水平位置、上体水平速度を用いて求められる。
前記S302の次に、S304を経て、定常歩容の初期時刻Tsから終端時刻Ts+Tcycまでの各時刻に対し、S306〜S312の処理が実行される。S306では、S302で初期化した誤差Merr_p(k)およびFerr_p(k)の時系列を用い、単純化モデル上で時刻kの目標ZMP(S022の中で決定されたZMP軌道パラメータにより規定される目標ZMP)まわりに−Merr_p(k)が発生し、且つ、並進床反力水平成分(単純化モデル上での歩容の全体重心の並進加速度水平成分に伴う慣性力に釣り合う並進床反力)にFerr_p(k)を加えたものが時刻kの床反力水平成分許容範囲(S022の中で決定されたパラメータにより規定される床反力水平成分許容範囲)に収まるように定常歩容の瞬時値の時系列が作成される。この場合、S306では、Merr_p(k)=0且つFerr_p(k)=0であるから、S306で作成される(仮作成される)歩容は、結局、単純化モデル上で目標ZMPを満足し、且つ、床反力水平成分許容範囲を満足する歩容となる。S306のより具体的な処理は後述するが、S306の処理の中では、各時刻kにおいて、作成中の歩容の運動によってセミフルモデル上で発生する目標ZMPまわりの床反力モーメント水平成分と該運動によってセミフルモデル上で発生する目標ZMPまわりの床反力モーメント水平成分との差Merr(k)が求められると共に、該運動によってセミフルモデル上で発生する並進床反力水平成分と該運動によって単純化モデル上で発生する並進床反力水平成分との差Ferr(k)が求められる。
Msmpl(k)=0−Merr_p(k) ……式3a
Fsmpl(k)=Fx'(k)−Ferr_p(k) ……式3b
すなわち、S410〜S414までの処理で作成する歩容は、目標ZMPまわりに−Merr_p(k)(より正確には、目標ZMPまわりの本来の目標床反力モーメント水平成分(=0)からMerr_p(k)を差し引いたもの)が発生するように作成されるので、その歩容の運動によって単純化モデル上で発生する目標ZMPまわりの床反力モーメントMsmpl(k)は、−Merr_p(k)に一致もしくはほぼ一致する。従って、Msmpl(k)は、前記式3aにより求めてもよい。なお、当然ながら、式3a中の「0」は省略してもよい。
Merr(k)=Msemifull(k)−Msmpl(k) ……式3c
Ferr(k)=Fsemifull(k)−Fsmpl(k) ……式3d
補足すると、本実施形態では、図13のS306およびS322のサブルーチン処理を同じにしたが、S306およびS322のサブルーチン処理を各別に構築してもよい。この場合、S306の処理では、当初からMerr_p(k)=0、Ferr_p(k)=0を前提としてアルゴリズムを構築しておけばよく、その場合には、図13のS302の処理は不要である。また、S306の処理の中のS415において、前記式3aを適用する場合には常にMsmpl(k)=0となるので、S322の処理の中では、S415の処理を省略してもよい。以上説明した補足事項は、後述する今回歩容の作成処理(後述のS028の処理)においても同様である。
Merr_p(k)=Merr_p(k)+p*ΔMerr(k) ……式10a
Ferr_p(k)=Ferr_p(k)+p*ΔFerr(k) ……式10b
これらの式10a、10bにおいて、pは1以下、0以上の所定の係数である。
また、図13のS300〜S312の処理の終了後、図13のS318〜S328の処理を実行する前の最初のS316の判断処理を省略してもよい。すなわち、図13のS300〜S312の処理の後、必ず1回は、S318〜S328の処理を実行するようにしてもよい。同様に、図18のS800〜S812の処理の終了後、図18のS818〜S828の処理を実行する前の最初のS816の判断処理を省略し、S818〜S828の処理が必ず1回は実行されるようにしてもよい。このようにした実施形態(以下、第2実施形態という)は、前記第2発明、第4発明、第5発明、第6発明、第10発明、第14発明、第16発明の実施形態となる。この場合のこれらの発明と第2実施形態との対応関係は、評価手段を除いて(第2実施形態では評価手段は備えられない)、前述の実施形態(本明細書で最初に説明した実施形態(以降、第1実施形態という))と、第1発明もしくは第3発明などとの対応関係と同様である。
Claims (18)
- 移動ロボットの運動によって発生する床反力の所定の床反力成分の許容範囲を設定し、その許容範囲を満足するように、移動ロボットの目標運動を少なくとも含む目標歩容を生成する歩容生成装置において、
前記移動ロボットの所定の第1動力学モデルを用いて該移動ロボットの運動を作成するようにあらかじめ構成された所定のアルゴリズムを実行することにより、前記床反力成分が前記許容範囲を満たすように、前記目標運動の仮値である仮運動を作成する仮目標運動作成手段と、
前記第1動力学モデルと該第1動力学モデルよりも動力学的精度が高い所定の第2動力学モデルとを用い、前記移動ロボットの任意の運動に対して、該運動によって前記第2動力学モデル上で発生する前記床反力成分と該運動によって前記第1動力学モデル上で発生する前記床反力成分との差である床反力成分誤差を求める床反力成分誤差算出手段と、
前記作成された仮運動から前記床反力成分誤差算出手段によって求められる床反力成分誤差が所定の許容誤差範囲に収まっているか否かを評価する評価手段と、
前記評価手段の評価において前記仮運動に対応する床反力成分誤差が前記許容誤差範囲に収まっている場合には、前記仮運動を前記目標運動として決定すると共に、前記仮運動に対応する床反力成分誤差が前記許容誤差範囲から逸脱している場合には、前記仮運動を少なくとも1回以上修正して前記目標運動を決定する目標運動決定手段とを備え、
前記仮運動に対応する床反力成分誤差が前記許容誤差範囲から逸脱している場合において、前記目標運動決定手段が前記仮運動の第n回目(n:n≧1である整数)の修正によって求める修正後運動を第n修正後運動とし、該第n修正後運動から前記床反力成分誤差算出手段により求められる床反力成分誤差を第n床反力成分誤差とし、前記仮運動を第0修正後運動とし、前記仮運動に対応する床反力成分誤差を第0床反力成分誤差としたとき、
前記目標運動決定手段は、第n修正後運動によって前記第1動力学モデル上で発生する前記床反力成分に、第n−1床反力成分誤差と少なくとも該第n−1床反力成分誤差に応じて決定した床反力補正量とのいずれかを加えたものが前記許容範囲を満たすように前記所定のアルゴリズムを実行することにより、第n修正後運動を決定する修正後運動決定処理と、前記決定した第n修正後運動に対応する第n床反力成分誤差と第n−1修正後運動に対応する第n−1床反力成分誤差との差、および、前記第n床反力成分誤差と少なくとも前記第n−1床反力成分誤差に応じて決定した前記床反力補正量との差のいずれか一方として定義される床反力成分誤差変化量ΔFMが0に収束したか否かを判断する収束判断処理とを少なくとも該床反力成分誤差変化量ΔFMが0に収束したと判断されるまで繰り返し、この繰り返し処理における最後の修正後運動決定処理で決定した修正後運動を前記目標運動として決定することを特徴とする移動ロボットの歩容生成装置。 - 移動ロボットの運動によって発生する床反力の所定の床反力成分の許容範囲を設定し、その許容範囲を満足するように、移動ロボットの目標運動を少なくとも含む目標歩容を生成する歩容生成装置において、
前記目標運動の仮値である仮運動を作成する仮目標運動作成手段と、
前記移動ロボットの所定の第1の動力学モデルと該第1動力学モデルよりも動力学的精度が高い所定の第2動力学モデルとを用い、前記移動ロボットの任意の運動に対して、該運動によって前記第2動力学モデル上で発生する前記床反力成分と該運動によって前記第1動力学モデル上で発生する前記床反力成分との差である床反力成分誤差を求める床反力成分誤差算出手段と、
前記仮運動を少なくとも1回以上修正して前記目標運動を決定する目標運動決定手段とを備え、
前記目標運動決定手段が前記仮運動の第n回目(n:n≧1である整数)の修正によって求める修正後運動を第n修正後運動とし、該第n修正後運動から前記床反力成分誤差算出手段により求められる床反力成分誤差を第n床反力成分誤差とし、前記仮運動を第0修正後運動とし、前記仮運動から前記床反力成分誤差算出手段により求められる床反力成分誤差を第0床反力成分誤差としたとき、
前記目標運動決定手段は、前記第1動力学モデルを用いて該移動ロボットの運動を作成するようにあらかじめ構成された所定のアルゴリズムを実行することにより、第n修正後運動によって前記第1動力学モデル上で発生する前記床反力成分に、第n−1床反力成分誤差と少なくとも該第n−1床反力成分誤差に応じて決定した床反力補正量とのいずれかを加えたものが前記許容範囲を満たすように第n修正後運動を決定する修正後運動決定処理と、前記決定した第n修正後運動に対応する第n床反力成分誤差と第n−1修正後運動に対応する第n−1床反力成分誤差との差、および、前記第n床反力成分誤差と少なくとも前記第n−1床反力成分誤差に応じて決定した前記床反力補正量との差のいずれか一方として定義される床反力成分誤差変化量ΔFMが0に収束したか否かを判断する収束判断処理とを少なくとも該床反力成分誤差変化量ΔFMが0に収束したと判断されるまで繰り返し、この繰り返し処理における最後の修正後運動決定処理で決定した修正後運動を前記目標運動として決定することを特徴とする移動ロボットの歩容生成装置。 - 移動ロボットの運動によって発生する床反力の所定の第1床反力成分の目標値と該第1床反力成分と異なる所定の第2床反力成分の許容範囲とを設定し、該第1床反力成分の目標値と第2床反力成分の許容範囲とを満足するように、移動ロボットの目標運動を少なくとも含む目標歩容を生成する歩容生成装置において、
前記移動ロボットの所定の第1動力学モデルを用いて該移動ロボットの運動を作成するようにあらかじめ構成された所定のアルゴリズムを実行することにより、前記第1床反力成分が前記目標値に一致し、且つ、前記第2床反力成分が前記許容範囲を満たすように、前記目標運動の仮値である仮運動を作成する仮目標運動作成手段と、
前記第1動力学モデルと該第1動力学モデルよりも動力学的精度が高い所定の第2動力学モデルとを用い、前記移動ロボットの任意の運動に対して、該運動によって前記第2動力学モデル上で発生する前記第1床反力成分と該運動によって前記第1動力学モデル上で発生する前記第1床反力成分との差である床反力成分誤差Aerrを求めると共に、該運動によって前記第2動力学モデル上で発生する前記第2床反力成分と該運動によって前記第1動力学モデル上で発生する前記第2床反力成分との差である床反力成分誤差Berrを求める床反力成分誤差算出手段と、
前記作成された仮運動から前記床反力成分誤差算出手段によって求められる床反力成分誤差Aerr,Berrのうちの床反力成分誤差Aerrが所定の第1許容誤差範囲に収まっているか否かを評価すると共に、床反力成分誤差Berrが所定の第2許容誤差範囲に収まっているか否かを評価する評価手段と、
前記評価手段の評価において前記仮運動に対応する床反力成分誤差Aerr,Berrの両者がそれぞれに対応する前記第1許容誤差範囲および第2許容誤差範囲に収まっている場合には、前記仮運動を前記目標運動として決定すると共に、前記仮運動に対応する床反力成分誤差Aerr,Berrの少なくともいずれか一方がそれに対応する前記第1許容誤差範囲または第2許容誤差範囲から逸脱している場合には、前記仮運動を少なくとも1回以上修正して前記目標運動を決定する目標運動決定手段とを備え、
前記仮運動に対応する床反力成分誤差Aerr,Berrの少なくともいずれか一方がそれに対応する前記第1許容誤差範囲または第2許容誤差範囲から逸脱している場合において、前記目標運動決定手段が前記仮運動の第n回目(n:n≧1である整数)の修正によって求める修正後運動を第n修正後運動とし、該第n修正後運動から前記床反力成分誤差算出手段により求められる床反力成分誤差Aerr,Berrをそれぞれ第n床反力成分誤差Aerr(n),Berr(n)とし、前記仮運動を第0修正後運動とし、前記仮運動に対応する床反力成分誤差Aerr,Berrをそれぞれ第0床反力成分誤差Aerr(0),Berr(0)としたとき、
前記目標運動決定手段は、第n修正後運動によって前記第1動力学モデル上で発生する前記第1床反力成分に、第n−1床反力成分誤差Aerr(n-1)と少なくとも該第n−1床反力成分誤差Aerr(n-1)に応じて決定した第1床反力補正量とのいずれかを加えたものが前記目標値に一致し、且つ、該第n修正後運動によって前記第1動力学モデル上で発生する前記第2床反力成分に、第n−1床反力成分誤差Berr(n-1)と少なくとも該第n−1床反力成分誤差Berr(n-1)に応じて決定した第2床反力補正量とのいずれかを加えたものが前記許容範囲を満たすように前記所定のアルゴリズムを実行することにより、第n修正後運動を決定する修正後運動決定処理と、前記決定した第n修正後運動に対応する第n床反力成分誤差Aerr(n)と第n−1修正後運動に対応する第n−1床反力成分誤差Aerr(n-1)との差、および、前記第n床反力成分誤差Aerr(n)と少なくとも前記第n−1床反力成分誤差Aerr(n-1)に応じて決定した前記第1床反力補正量との差のいずれか一方として定義される床反力成分誤差変化量ΔAerrと、前記決定した第n修正後運動に対応する第n床反力成分誤差Berr(n)と第n−1修正後運動に対応する第n−1床反力成分誤差Berr(n-1)との差、および、前記第n床反力成分誤差Berr(n)と少なくとも前記第n−1床反力成分誤差Berr(n-1)に応じて決定した前記第2床反力補正量との差のいずれか一方として定義される床反力成分誤差変化量ΔBerrとがそれぞれ0に収束したか否かを判断する収束判断処理とを少なくとも該床反力成分誤差変化量ΔAerr,ΔBerrの両者が0に収束したと判断されるまで繰り返し、この繰り返し処理における最後の前記修正後運動決定処理で決定した修正後運動を前記目標運動として決定することを特徴とする移動ロボットの歩容生成装置。 - 移動ロボットの運動によって発生する床反力の所定の第1床反力成分の目標値と該第1床反力成分と異なる所定の第2床反力成分の許容範囲とを設定し、該第1床反力成分の目標値と第2床反力成分の許容範囲とを満足するように、移動ロボットの目標運動を少なくとも含む目標歩容を生成する歩容生成装置において、
前記目標運動の仮値である仮運動を作成する仮目標運動作成手段と、
前記移動ロボットの所定の第1動力学モデルと該第1動力学モデルよりも動力学的精度が高い所定の第2動力学モデルとを用い、前記移動ロボットの任意の運動に対して、該運動によって前記第2動力学モデル上で発生する前記第1床反力成分と該運動によって前記第1動力学モデル上で発生する前記第1床反力成分との差である床反力成分誤差Aerrを求めると共に、該運動によって前記第2動力学モデル上で発生する前記第2床反力成分と該運動によって前記第1動力学モデル上で発生する前記第2床反力成分との差である床反力成分誤差Berrを求める床反力成分誤差算出手段と、
前記仮運動を少なくとも1回以上修正して前記目標運動を決定する目標運動決定手段とを備え、
前記目標運動決定手段が前記仮運動の第n回目(n:n≧1である整数)の修正によって求める修正後運動を第n修正後運動とし、該第n修正後運動から前記床反力成分誤差算出手段により求められる床反力成分誤差Aerr,Berrをそれぞれ第n床反力成分誤差Aerr(n),Berr(n)とし、前記仮運動を第0修正後運動とし、前記仮運動から前記床反力成分誤差算出手段により求められる床反力成分誤差Aerr,Berrをそれぞれ第0床反力成分誤差Aerr(0),Berr(0)としたとき、
前記目標運動決定手段は、前記第1動力学モデルを用いて該移動ロボットの運動を作成するようにあらかじめ構成された所定のアルゴリズムを実行することにより、第n修正後運動によって前記第1動力学モデル上で発生する前記第1床反力成分に、第n−1床反力成分誤差Aerr(n-1)と少なくとも該第n−1床反力成分誤差Aerr(n-1)に応じて決定した第1床反力補正量とのいずれかを加えたものが前記目標値に一致し、且つ、該第n修正後運動によって前記第1動力学モデル上で発生する前記第2床反力成分に、第n−1床反力成分誤差Berr(n-1)と少なくとも該第n−1床反力成分誤差Berr(n-1)に応じて決定した第2床反力補正量とのいずれかを加えたものが前記許容範囲を満たすように第n修正後運動を決定する修正後運動決定処理と、前記決定した第n修正後運動に対応する第n床反力成分誤差Aerr(n)と第n−1修正後運動に対応する第n−1床反力成分誤差Aerr(n-1)との差、および、前記第n床反力成分誤差Aerr(n)と少なくとも前記第n−1床反力成分誤差Aerr(n-1)に応じて決定した前記第1床反力補正量との差のいずれか一方として定義される床反力成分誤差変化量ΔAerrと、前記決定した第n修正後運動に対応する第n床反力成分誤差Berr(n)と第n−1修正後運動に対応する第n−1床反力成分誤差Berr(n-1)との差、および、前記第n床反力成分誤差Berr(n)と少なくとも前記第n−1床反力成分誤差Berr(n-1)に応じて決定した前記第2床反力補正量との差のいずれか一方として定義される床反力成分誤差変化量ΔBerrとがそれぞれ0に収束したか否かを判断する収束判断処理とを少なくとも該床反力成分誤差変化量ΔAerr,ΔBerrの両者が0に収束したと判断されるまで繰り返し、この繰り返し処理における最後の前記修正後運動決定処理で決定した修正後運動を前記目標運動として決定することを特徴とする移動ロボットの歩容生成装置。 - 前記床反力成分は、前記移動ロボットに作用する床反力の並進床反力水平成分であることを特徴とする請求項1または2記載の移動ロボットの歩容生成装置。
- 前記第1床反力成分は、前記移動ロボットに作用する床反力の、所定の作用点まわりの床反力モーメント水平成分であり、前記第2床反力成分は、前記移動ロボットに作用する床反力の並進床反力水平成分であることを特徴とする請求項3または4記載の移動ロボットの歩容生成装置。
- 移動ロボットの目標ZMPと該移動ロボットの運動によって発生する並進床反力水平成分の許容範囲とを設定し、該目標ZMPと並進床反力水平成分の許容範囲とを満足するように、移動ロボットの目標運動を少なくとも含む目標歩容を生成する歩容生成装置において、
前記移動ロボットの所定の第1動力学モデルを用いて該移動ロボットの運動を作成するようにあらかじめ構成された所定のアルゴリズムを実行することにより、前記目標ZMPを満足し、且つ、前記並進床反力水平成分が前記許容範囲を満たすように、前記目標運動の仮値である仮運動を作成する仮目標運動作成手段と、
前記第1動力学モデルと該第1動力学モデルよりも動力学的精度が高い所定の第2動力学モデルとを用い、前記移動ロボットの任意の運動に対して、該運動によって前記第2動力学モデル上で算出されるZMPと該運動によって前記第1動力学モデル上で算出されるZMPとの差であるZMP誤差ZMPerrを求めると共に、該運動によって前記第2動力学モデル上で発生する前記並進床反力水平成分と該運動によって前記第1動力学モデル上で発生する前記並進床反力水平成分との差である並進床反力水平成分誤差Ferrを求める誤差算出手段と、
前記作成された仮運動から前記誤差算出手段によって求められるZMP誤差ZMPerrおよび並進床反力水平成分誤差FerrのうちのZMP誤差ZMPerrが所定の第1許容誤差範囲に収まっているか否かを評価すると共に、並進床反力水平成分誤差Ferrが所定の第2許容誤差範囲に収まっているか否かを評価する評価手段と、
前記評価手段の評価において前記仮運動に対応するZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferrの両者がそれぞれに対応する前記第1許容誤差範囲および第2許容誤差範囲に収まっている場合には、前記仮運動を前記目標運動として決定すると共に、前記仮運動に対応するZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferrの少なくともいずれか一方がそれに対応する前記第1許容誤差範囲または第2許容誤差範囲から逸脱している場合には、前記仮運動を少なくとも1回以上修正して前記目標運動を決定する目標運動決定手段とを備え、
前記仮運動に対応するZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferrの少なくともいずれか一方がそれに対応する前記第1許容誤差範囲または第2許容誤差範囲から逸脱している場合において、前記目標運動決定手段が前記仮運動の第n回目(n:n≧1である整数)の修正によって求める修正後運動を第n修正後運動とし、該第n修正後運動から前記床反力成分誤差算出手段により求められるZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferrをそれぞれ第nZMP誤差ZMPerr(n)および第n並進床反力水平成分誤差Ferr(n)とし、前記仮運動を第0修正後運動とし、前記仮運動に対応するZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferrをそれぞれ第0ZMP誤差ZMPerr(0)および第0並進床反力水平成分誤差Ferr(0)としたとき、
前記目標運動決定手段は、第n修正後運動によって前記第1動力学モデル上で算出されるZMPに、第n−1ZMP誤差ZMPerr(n-1)と少なくとも該第n−1ZMP誤差ZMPerr(n-1)に応じて決定したZMP補正量とのいずれかを加えたものが前記目標ZMPに一致し、且つ、該第n修正後運動によって前記第1動力学モデル上で発生する前記並進床反力水平成分に、第n−1並進床反力水平成分誤差Ferr(n-1)と少なくとも該第n−1並進床反力水平成分誤差Ferr(n-1)に応じて決定した床反力補正量とのいずれかを加えたものが前記許容範囲を満たすように前記所定のアルゴリズムを実行することにより、第n修正後運動を決定する修正後運動決定処理と、前記決定した第n修正後運動に対応する第nZMP誤差ZMPerr(n)と第n−1修正後運動に対応する第n−1ZMP誤差ZMPerr(n-1)との差、および、前記第nZMP誤差ZMPerr(n)と少なくとも前記第n−1ZMP誤差ZMPerr(n-1)に応じて決定した前記ZMP補正量との差のいずれか一方として定義されるZMP誤差変化量ΔZMPerrと、前記決定した第n修正後運動に対応する第n並進床反力水平成分誤差Ferr(n)と第n−1修正後運動に対応する第n−1並進床反力水平成分誤差Ferr(n-1)との差、および、第n並進床反力水平成分誤差Ferr(n)と少なくとも前記第n−1並進床反力水平成分誤差Ferr(n-1)に応じて決定した前記床反力補正量との差のいずれか一方として定義される並進床反力水平成分誤差変化量ΔFerrとがそれぞれ0に収束したか否かを判断する収束判断処理とを少なくとも該ZMP誤差変化量ΔZMPerrおよび並進床反力水平成分誤差変化量ΔFerrの両者が0に収束したと判断されるまで繰り返し、この繰り返し処理における最後の前記修正後運動決定処理で決定した修正後運動を前記目標運動として決定することを特徴とする移動ロボットの歩容生成装置。 - 移動ロボットの目標ZMPと該移動ロボットの運動によって発生する並進床反力水平成分の許容範囲とを設定し、該目標ZMPと並進床反力水平成分の許容範囲とを満足するように、移動ロボットの目標運動を少なくとも含む目標歩容を生成する歩容生成装置において、
前記目標運動の仮値である仮運動を作成する仮目標運動作成手段と、
前記移動ロボットの所定の第1動力学モデルと該第1動力学モデルよりも動力学的精度が高い所定の第2動力学モデルとを用い、前記移動ロボットの任意の運動に対して、該運動によって前記第2動力学モデル上で算出されるZMPと該運動によって前記第1動力学モデル上で算出されるZMPとの差であるZMP誤差ZMPerrを求めると共に、該運動によって前記第2動力学モデル上で発生する前記並進床反力水平成分と該運動によって前記第1動力学モデル上で発生する前記並進床反力水平成分との差である並進床反力水平成分誤差Ferrを求める誤差算出手段と、
前記仮運動を少なくとも1回以上修正して前記目標運動を決定する目標運動決定手段とを備え、
前記目標運動決定手段が前記仮運動の第n回目(n:n≧1である整数)の修正によって求める修正後運動を第n修正後運動とし、該第n修正後運動から前記床反力成分誤差算出手段により求められるZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferrをそれぞれ第nZMP誤差ZMPerr(n)および第n並進床反力水平成分誤差Ferr(n)とし、前記仮運動を第0修正後運動とし、前記仮運動から前記床反力成分誤差算出手段により求められるZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferrをそれぞれ第0ZMP誤差ZMPerr(0)および第0並進床反力水平成分誤差Ferr(0)としたとき、
前記目標運動決定手段は、前記第1動力学モデルを用いて該移動ロボットの運動を作成するようにあらかじめ構成された所定のアルゴリズムを実行することにより、第n修正後運動によって前記第1動力学モデル上で算出されるZMPに、第n−1ZMP誤差ZMPerr(n-1)と少なくとも該第n−1ZMP誤差ZMPerr(n-1)に応じて決定したZMP補正量とのいずれかを加えたものが前記目標ZMPに一致し、且つ、該第n修正後運動によって前記第1動力学モデル上で発生する前記並進床反力水平成分に、第n−1並進床反力水平成分誤差Ferr(n-1)と少なくとも該第n−1並進床反力水平成分誤差Ferr(n-1)に応じて決定した床反力補正量とのいずれかを加えたものが前記許容範囲を満たすように第n修正後運動を決定する修正後運動決定処理と、前記決定した第n修正後運動に対応する第nZMP誤差ZMPerr(n)と第n−1修正後運動に対応する第n−1ZMP誤差ZMPerr(n-1)との差、および、前記第nZMP誤差ZMPerr(n)と少なくとも前記第n−1ZMP誤差ZMPerr(n-1)に応じて決定した前記ZMP補正量との差のいずれか一方として定義されるZMP誤差変化量ΔZMPerrと、前記決定した第n修正後運動に対応する第n並進床反力水平成分誤差Ferr(n)と第n−1修正後運動に対応する第n−1並進床反力水平成分誤差Ferr(n-1)との差、および、第n並進床反力水平成分誤差Ferr(n)と少なくとも前記第n−1並進床反力水平成分誤差Ferr(n-1)に応じて決定した前記床反力補正量との差のいずれか一方として定義される並進床反力水平成分誤差変化量ΔFerrとがそれぞれ0に収束したか否かを判断する収束判断処理とを少なくとも該ZMP誤差変化量ΔZMPerrおよび並進床反力水平成分誤差変化量ΔFerrの両者が0に収束したと判断されるまで繰り返し、この繰り返し処理における最後の前記修正後運動決定処理で決定した修正後運動を前記目標運動として決定することを特徴とする移動ロボットの歩容生成装置。 - 前記床反力成分誤差算出手段は、前記仮運動に対しては、該仮運動によって前記第1動力学モデル上で発生する前記第1床反力成分として、前記目標値を用いることにより該仮運動に対応する前記床反力成分誤差Aerrを算出し、前記仮運動以外の第n修正後運動に対しては、該第n修正後運動によって前記第1動力学モデル上で発生する前記第1床反力成分として、前記目標値から前記第n−1床反力成分誤差Aerr(n-1)と少なくとも該第n−1床反力成分誤差Aerr(n-1)に応じて決定した前記第1床反力補正量とのいずれかを差し引いたものを用いることにより該第n修正後運動に対応する前記床反力成分誤差Aerrを算出することを特徴とする請求項3記載の移動ロボットの歩容生成装置。
- 前記仮運動作成手段は、少なくとも前記第1動力学モデル上で前記第1床反力成分の目標値を満足するように前記仮運動を作成する手段であり、
前記床反力成分誤差算出手段は、前記仮運動に対しては、該仮運動によって前記第1動力学モデル上で発生する前記第1床反力成分として、前記目標値を用いることにより該仮運動に対応する前記床反力成分誤差Aerrを算出し、前記仮運動以外の第n修正後運動に対しては、該第n修正後運動によって前記第1動力学モデル上で発生する前記第1床反力成分として、前記目標値から前記第n−1床反力成分誤差Aerr(n-1)と少なくとも該第n−1床反力成分誤差Aerr(n-1)に応じて決定した前記第1床反力補正量とのいずれかを差し引いたものを用いることにより該第n修正後運動に対応する前記床反力成分誤差Aerrを算出することを特徴とする請求項4記載の移動ロボットの歩容生成装置。 - 前記床反力成分誤差算出手段は、前記仮運動に対しては、該仮運動によって前記第1動力学モデル上で算出されるZMPとして、前記目標ZMPを用いることにより該仮運動に対応する前記ZMP誤差ZMPerrを算出し、前記仮運動以外の第n修正後運動に対しては、該第n修正後運動によって前記第1動力学モデル上で算出されるZMPとして、前記目標ZMPから前記第n−1ZMP誤差ZMPerr(n-1)と少なくとも該第n−1ZMP誤差ZMPerr(n-1)に応じて決定した前記ZMP補正量とのいずれかを差し引いたものを用いることにより該第n修正後運動に対応するZMP誤差ZMPerrを算出することを特徴とする請求項7記載の移動ロボットの歩容生成装置。
- 前記仮運動作成手段は、少なくとも前記第1動力学モデル上で前記目標ZMPを満足するように前記仮運動を作成する手段であり、
前記床反力成分誤差算出手段は、前記仮運動に対しては、該仮運動によって前記第1動力学モデル上で算出されるZMPとして、前記目標ZMPを用いることにより該仮運動に対応する前記ZMP誤差ZMPerrを算出し、前記仮運動以外の第n修正後運動に対しては、該第n修正後運動によって前記第1動力学モデル上で算出されるZMPとして、前記目標ZMPから前記第n−1ZMP誤差ZMPerr(n-1)と少なくとも該第n−1ZMP誤差ZMPerr(n-1)に応じて決定した前記ZMP補正量とのいずれかを差し引いたものを用いることにより該第n修正後運動に対応するZMP誤差ZMPerrを算出することを特徴とする請求項8記載の移動ロボットの歩容生成装置。 - 前記目標運動、仮運動および修正後運動は、所定期間における前記移動ロボットの運動の瞬時値の時系列から構成されるものであり、
前記床反力成分誤差算出手段により求められる前記第m床反力成分誤差(m:m≧0である整数)は、第m修正後運動によって該第m修正後運動の各時刻に前記第2動力学モデル上で発生する前記床反力成分の瞬時値と該第m修正後運動によって該時刻に前記第1動力学モデル上で発生する前記床反力成分の瞬時値との差の前記所定期間における時系列から構成され、
前記評価手段の前記所定の許容誤差範囲は、前記第0床反力成分誤差を構成する時系列のパターンのうちの所定の第1特徴量に対する許容誤差範囲であり、
前記目標運動決定手段の前記修正後運動決定処理は、前記第n修正後運動の各時刻tにおいて、該第n修正後運動によって時刻tに前記第1動力学モデル上で発生する前記床反力成分の瞬時値FM(t)に、前記第n−1床反力成分誤差の該時刻tでの値と少なくとも該値に応じて決定した床反力補正量の値とのいずれかを加えたものが該時刻tでの前記許容範囲を満たすように該時刻tにおける第n修正後運動の瞬時値を決定する処理であり、
前記床反力成分誤差変化量ΔFMは、前記第n床反力成分誤差の各時刻での値と第n−1床反力成分誤差の該時刻での値との差、および、前記第n床反力成分誤差の各時刻での値と少なくとも第n−1床反力成分誤差の該時刻での値に応じて決定した前記床反力補正量の値との差のいずれか一方の前記所定期間における時系列から構成され、
前記目標運動決定手段の前記収束判断処理は、前記床反力成分誤差変化量ΔFMを構成する時系列のパターンのうちの所定の第2特徴量が所定の許容変化量範囲に収まったときに該床反力成分誤差変化量ΔFMが0に収束したと判断する処理であることを特徴とする請求項1記載の移動ロボットの歩容生成装置。 - 前記目標運動、仮運動および修正後運動は、所定期間における前記移動ロボットの運動の瞬時値の時系列から構成されるものであり、
前記床反力成分誤差算出手段により求められる前記第m床反力成分誤差(m:m≧0である整数)は、第m修正後運動によって該第m修正後運動の各時刻に前記第2動力学モデル上で発生する前記床反力成分の瞬時値と該第m修正後運動によって該時刻に前記第1動力学モデル上で発生する前記床反力成分の瞬時値との差の前記所定期間における時系列から構成され、
前記目標運動決定手段の前記修正後運動決定処理は、前記第n修正後運動の各時刻tにおいて、該第n修正後運動によって時刻tに前記第1動力学モデル上で発生する前記床反力成分の瞬時値FM(t)に、前記第n−1床反力成分誤差の該時刻tでの値と少なくとも該値に応じて決定した前記床反力補正量の値とのいずれかを加えたものが該時刻tでの前記許容範囲を満たすように該時刻tにおける第n修正後運動の瞬時値を決定する処理であり、
前記床反力成分誤差変化量ΔFMは、前記第n床反力成分誤差の各時刻での値と第n−1床反力成分誤差の該時刻での値との差、および、前記第n床反力成分誤差の各時刻での値と少なくとも第n−1床反力成分誤差の該時刻での値に応じて決定した前記床反力補正量の値との差のいずれか一方の前記所定期間における時系列から構成され、
前記目標運動決定手段の前記収束判断処理は、前記床反力成分誤差変化量ΔFMを構成する時系列のパターンのうちの所定の特徴量が所定の許容変化量範囲に収まったときに該床反力成分誤差変化量ΔFMが0に収束したと判断する処理であることを特徴とする請求項2記載の移動ロボットの歩容生成装置。 - 前記目標運動、仮運動および修正後運動は、所定期間における前記移動ロボットの運動の瞬時値の時系列から構成されるものであり、
前記床反力成分誤差算出手段により求められる前記第m床反力成分誤差Aerr,Berr(m:m≧0である整数)のうちの第m床反力成分誤差Aerrは、第m修正後運動によって該第m修正後運動の各時刻に前記第2動力学モデル上で発生する前記第1床反力成分の瞬時値と該第m修正後運動によって該時刻に前記第1動力学モデル上で発生する前記第1床反力成分の瞬時値との差の前記所定期間における時系列から構成されると共に、前記第m床反力成分誤差Berrは、第m修正後運動の各時刻に前記第2動力学モデル上で発生する前記第2床反力成分の瞬時値と該第m修正後運動によって該時刻に前記第1動力学モデル上で発生する前記第2床反力成分の瞬時値との差の前記所定期間における時系列から構成され、
前記評価手段の前記所定の第1許容誤差範囲は、前記第0床反力成分誤差Aerr(0)を構成する時系列パターンのうちの所定の第1特徴量に対する許容誤差範囲であると共に前記所定の第2許容誤差範囲は、前記第0床反力成分誤差Berr(0)を構成する時系列のパターンのうちの所定の第2特徴量に対する許容誤差範囲であり、
前記目標運動決定手段の前記修正後運動決定処理は、前記第n修正後運動の各時刻tにおいて、該第n修正後運動によって時刻tに前記第1動力学モデル上で発生する前記第1床反力成分の瞬時値に、前記第n−1床反力成分誤差Aerr(n-1)の時刻tでの値と少なくとも該値に応じて決定した前記第1床反力補正量の値とのいずれかを加えたものが該時刻tでの前記目標値を満足し、且つ、該第n修正後運動によって該時刻tに前記第1動力学モデル上で発生する前記第2床反力成分の瞬時値に、前記第n−1床反力成分誤差Berr(n-1)の時刻tでの値と少なくとも該値に応じて決定した前記第2床反力補正量の値とのいずれかを加えたものが該時刻tでの前記許容範囲を満たすように該時刻tにおける第n修正後運動の瞬時値を決定する処理であり、
前記床反力成分誤差変化量ΔAerrは、前記第n床反力成分誤差Aerr(n)の各時刻での値と第n−1床反力成分誤差Aerr(n-1)の該時刻での値との差、および、前記第n床反力成分誤差Aerr(n)の各時刻での値と少なくとも第n−1床反力成分誤差Aerr(n-1)の該時刻での値に応じて決定した前記第1床反力補正量の値との差のいずれか一方の前記所定期間における時系列から構成されると共に、前記床反力成分誤差変化量ΔBerrは、前記第n床反力成分誤差Berr(n)の各時刻での値と第n−1床反力成分誤差Berr(n-1)の該時刻での値との差、および、前記第n床反力成分誤差Berr(n)の各時刻での値と少なくとも第n−1床反力成分誤差Berr(n-1)の該時刻での値に応じて決定した前記第2床反力補正量の値との差のいずれか一方の前記所定期間における時系列から構成され、
前記目標運動決定手段の前記収束判断処理は、前記床反力成分誤差変化量ΔAerrを構成する時系列のパターンのうちの所定の第3特徴量が該第3特徴量に対する所定の許容変化量範囲に収まったときに該床反力成分誤差変化量ΔAerrが0に収束したと判断すると共に、前記床反力成分誤差変化量ΔBerrを構成する時系列のパターンのうちの所定の第4特徴量が該第4特徴量に対する所定の許容変化量範囲に収まったときに該床反力成分誤差変化量ΔBerrが0に収束したと判断する処理であることを特徴とする請求項3記載の移動ロボットの歩容生成装置。 - 前記目標運動、仮運動および修正後運動は、所定期間における前記移動ロボットの運動の瞬時値の時系列から構成されるものであり、
前記床反力成分誤差算出手段により求められる前記第m床反力成分誤差Aerr,Berr(m:m≧0である整数)のうちの第m床反力成分誤差Aerrは、第m修正後運動によって該第m修正後運動の各時刻に前記第2動力学モデル上で発生する前記第1床反力成分の瞬時値と該第m修正後運動によって該時刻に前記第1動力学モデル上で発生する前記第1床反力成分の瞬時値との差の前記所定期間における時系列から構成されると共に、前記第m床反力成分誤差Berrは、第m修正後運動の各時刻に前記第2動力学モデル上で発生する前記第2床反力成分の瞬時値と該第m修正後運動によって該時刻に前記第1動力学モデル上で発生する前記第2床反力成分の瞬時値との差の前記所定期間における時系列から構成され、
前記目標運動決定手段の前記修正後運動決定処理は、前記第n修正後運動の各時刻tにおいて、該第n修正後運動によって時刻tに前記第1動力学モデル上で発生する前記第1床反力成分の瞬時値に、前記第n−1床反力成分誤差Aerr(n-1)の時刻tでの値と少なくとも該値に応じて決定した前記第1床反力補正量の値とのいずれかを加えたものが該時刻tでの前記目標値を満足し、且つ、該第n修正後運動によって該時刻tに前記第1動力学モデル上で発生する前記第2床反力成分の瞬時値に、前記第n−1床反力成分誤差Berr(n-1)の時刻tでの値と少なくとも該値に応じて決定した前記第2床反力補正量の値とのいずれかを加えたものが該時刻tでの前記許容範囲を満たすように該時刻tにおける第n修正後運動の瞬時値を決定する処理であり、
前記床反力成分誤差変化量ΔAerrは、前記第n床反力成分誤差Aerr(n)の各時刻での値と第n−1床反力成分誤差Aerr(n-1)の該時刻での値との差、および、前記第n床反力成分誤差Aerr(n)の各時刻での値と少なくとも第n−1床反力成分誤差Aerr(n-1)の該時刻での値に応じて決定した前記第1床反力補正量の値との差のいずれか一方の前記所定期間における時系列から構成されると共に、前記床反力成分誤差変化量ΔBerrは、前記第n床反力成分誤差Berr(n)の各時刻での値と第n−1床反力成分誤差Berr(n-1)の該時刻での値との差、および、前記第n床反力成分誤差Berr(n)の各時刻での値と少なくとも第n−1床反力成分誤差Berr(n-1)の該時刻での値に応じて決定した前記第2床反力補正量の値との差のいずれか一方の前記所定期間における時系列から構成され、
前記目標運動決定手段の前記収束判断処理は、前記床反力成分誤差変化量ΔAerrを構成する時系列のパターンのうちの所定の第3特徴量が該第3特徴量に対する所定の許容変化量範囲に収まったときに該床反力成分誤差変化量ΔAerrが0に収束したと判断すると共に、前記床反力成分誤差変化量ΔBerrを構成する時系列のパターンのうちの前記第4特徴量が該第4特徴量に対する所定の許容変化量範囲に収まったときに該床反力成分誤差変化量ΔBerrが0に収束したと判断する処理であることを特徴とする請求項3記載の移動ロボットの歩容生成装置。 - 前記目標運動、仮運動および修正後運動は、所定期間における前記移動ロボットの運動の瞬時値の時系列から構成されるものであり、
前記床反力成分誤差算出手段により求められる前記第mZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferr(m:m≧0である整数)のうちの第mZMP誤差ZMPerrは、第m修正後運動によって該第m修正後運動の各時刻に前記第2動力学モデル上で算出されるZMPの瞬時値と該第m修正後運動によって該時刻に前記第1動力学モデル上で算出されるZMPの瞬時値との差の前記所定期間における時系列から構成されると共に、前記第m並進床反力水平成分誤差Berrは、第m修正後運動の各時刻に前記第2動力学モデル上で発生する並進床反力水平成分の瞬時値と該第m修正後運動によって該時刻に前記第1動力学モデル上で発生する並進床反力水平成分の瞬時値との差の前記所定期間における時系列から構成され、
前記評価手段の前記所定の第1許容誤差範囲は、前記第0ZMP誤差ZMPerr(0)を構成する時系列のパターンのうちの所定の第1特徴量に対する許容誤差範囲であると共に、前記所定の第2許容誤差範囲は、前記第0並進床反力水平成分誤差Ferr(0)を構成する時系列のパターンのうちの所定の第2特徴量に対する許容誤差範囲であり、
前記目標運動決定手段の前記修正後運動決定処理は、前記第n修正後運動の各時刻tにおいて、該第n修正後運動によって時刻tに前記第1動力学モデル上で算出されるZMPの瞬時値に、前記第n−1ZMP誤差ZMPerr(n-1)の時刻tでの値と少なくとも該値に応じて決定した前記ZMP補正量の値とのいずれかを加えたものが該時刻tでの前記目標ZMPを満足し、且つ、該第n修正後運動によって該時刻tに前記第1動力学モデル上で発生する並進床反力水平成分の瞬時値に、前記第n−1並進床反力水平成分誤差Ferr(n-1)の時刻tでの値と少なくとも該値に応じて決定した前記床反力補正量の値とのいずれかを加えたものが該時刻tでの前記許容範囲を満たすように該時刻tにおける第n修正後運動の瞬時値を決定する処理であり、
前記ZMP誤差変化量ΔZMPerrは、前記第nZMP誤差ZMPerr(n)の各時刻での値と第n−1ZMP誤差ZMPerr(n-1)の該時刻での値との差、および、前記第nZMP誤差ZMPerr(n)の各時刻での値と少なくとも第n−1ZMP誤差ZMPerr(n-1)の該時刻での値に応じて決定した前記ZMP補正量の値との差のいずれか一方の前記所定期間における時系列から構成されると共に、前記並進床反力水平成分誤差変化量ΔFerrは、前記第n並進床反力水平成分誤差Ferr(n)の各時刻での値と第n−1並進床反力水平成分誤差Ferr(n-1)の該時刻での値との差、および、前記第n並進床反力水平成分誤差Ferr(n)の各時刻での値と少なくとも第n−1並進床反力水平成分誤差Ferr(n-1)の該時刻での値に応じて決定した前記床反力補正量の値との差のいずれか一方の前記所定期間における時系列から構成され、
前記目標運動決定手段の前記収束判断処理は、前記ZMP誤差変化量ΔZMPerrを構成する時系列のパターンのうちの所定の第3特徴量が該第3特徴量に対する所定の許容変化量範囲に収まったときに該ZMP成分誤差変化量ΔZMPerrが0に収束したと判断すると共に、前記並進床反力水平成分誤差変化量ΔFerrを構成する時系列のパターンのうちの所定の第4特徴量が該第4特徴量に対する所定の許容変化量範囲に収まったときに該並進床反力水平成分誤差変化量ΔFerrが0に収束したと判断する処理であることを特徴とする請求項7記載の移動ロボットの歩容生成装置。 - 前記目標運動、仮運動および修正後運動は、所定期間における前記移動ロボットの運動の瞬時値の時系列から構成されるものであり、
前記床反力成分誤差算出手段により求められる前記第mZMP誤差ZMPerrおよび並進床反力水平成分誤差Ferr(m:m≧0である整数)のうちの第mZMP誤差ZMPerrは、第m修正後運動によって該第m修正後運動の各時刻に前記第2動力学モデル上で算出されるZMPの瞬時値と該第m修正後運動によって該時刻に前記第1動力学モデル上で算出されるZMPの瞬時値との差の前記所定期間における時系列から構成されると共に、前記第m並進床反力水平成分誤差Berrは、第m修正後運動の各時刻に前記第2動力学モデル上で発生する並進床反力水平成分の瞬時値と該第m修正後運動によって該時刻に前記第1動力学モデル上で発生する並進床反力水平成分の瞬時値との差の前記所定期間における時系列から構成され、
前記目標運動決定手段の前記修正後運動決定処理は、前記第n修正後運動の各時刻tにおいて、該第n修正後運動によって時刻tに前記第1動力学モデル上で算出されるZMPの瞬時値に、前記第n−1ZMP誤差ZMPerr(n-1)の時刻tでの値と少なくとも該値に応じて決定した前記ZMP補正量の値とのいずれかを加えたものが該時刻tでの前記目標ZMPを満足し、且つ、該第n修正後運動によって該時刻tに前記第1動力学モデル上で発生する並進床反力水平成分の瞬時値に、前記第n−1並進床反力水平成分誤差Ferr(n-1)の時刻tでの値と少なくとも該値に応じて決定した前記床反力補正量の値とのいずれかを加えたものが該時刻tでの前記許容範囲を満たすように該時刻tにおける第n修正後運動の瞬時値を決定する処理であり、
前記ZMP誤差変化量ΔZMPerrは、前記第nZMP誤差ZMPerr(n)の各時刻での値と第n−1ZMP誤差ZMPerr(n-1)の該時刻での値との差、および、前記第nZMP誤差ZMPerr(n)の各時刻での値と少なくとも第n−1ZMP誤差ZMPerr(n-1)の該時刻での値に応じて決定した前記ZMP補正量の値との差のいずれか一方の前記所定期間における時系列から構成されると共に、前記並進床反力水平成分誤差変化量ΔFerrは、前記第n並進床反力水平成分誤差Ferr(n)の各時刻での値と第n−1並進床反力水平成分誤差Ferr(n-1)の該時刻での値との差、および、前記第n並進床反力水平成分誤差Ferr(n)の各時刻での値と少なくとも第n−1並進床反力水平成分誤差Ferr(n-1)の該時刻での値に応じて決定した前記床反力補正量の値との差のいずれか一方の前記所定期間における時系列から構成され、
前記目標運動決定手段の前記収束判断処理は、前記ZMP誤差変化量ΔZMPerrを構成する時系列のパターンのうちの所定の第3特徴量が該第3特徴量に対する所定の許容変化量範囲に収まったときに該ZMP成分誤差変化量ΔZMPerrが0に収束したと判断すると共に、前記並進床反力水平成分誤差変化量ΔFerrを構成する時系列のパターンのうちの所定の第4特徴量が該第4特徴量に対する所定の許容変化量範囲に収まったときに該並進床反力水平成分誤差変化量ΔFerrが0に収束したと判断する処理であることを特徴とする請求項8記載の移動ロボットの歩容生成装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006510511A JP4805818B2 (ja) | 2004-02-27 | 2005-02-28 | 移動ロボットの歩容生成装置 |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004052823 | 2004-02-27 | ||
JP2004052823 | 2004-02-27 | ||
JP2006510511A JP4805818B2 (ja) | 2004-02-27 | 2005-02-28 | 移動ロボットの歩容生成装置 |
PCT/JP2005/003346 WO2005082582A1 (ja) | 2004-02-27 | 2005-02-28 | 移動ロボットの歩容生成装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2005082582A1 JPWO2005082582A1 (ja) | 2007-10-25 |
JP4805818B2 true JP4805818B2 (ja) | 2011-11-02 |
Family
ID=34908714
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2006510510A Expired - Fee Related JP4805817B2 (ja) | 2004-02-27 | 2005-02-28 | 移動ロボットの歩容生成装置 |
JP2006510511A Expired - Fee Related JP4805818B2 (ja) | 2004-02-27 | 2005-02-28 | 移動ロボットの歩容生成装置 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2006510510A Expired - Fee Related JP4805817B2 (ja) | 2004-02-27 | 2005-02-28 | 移動ロボットの歩容生成装置 |
Country Status (5)
Country | Link |
---|---|
US (2) | US7715945B2 (ja) |
EP (2) | EP1738879B1 (ja) |
JP (2) | JP4805817B2 (ja) |
KR (2) | KR101140812B1 (ja) |
WO (2) | WO2005082583A1 (ja) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7649536B1 (en) * | 2006-06-16 | 2010-01-19 | Nvidia Corporation | System, method, and computer program product for utilizing natural motions of a user to display intuitively correlated reactions |
KR100835354B1 (ko) * | 2006-07-05 | 2008-06-04 | 삼성전자주식회사 | 보행로봇 및 그의 제어방법 |
US7943874B2 (en) * | 2007-09-18 | 2011-05-17 | Honeywell International Inc. | Ground contact switch for personal navigation system |
JP5251253B2 (ja) * | 2008-05-20 | 2013-07-31 | 株式会社安川電機 | 脚式歩行ロボットの安定化制御装置 |
JP5232124B2 (ja) * | 2009-10-28 | 2013-07-10 | 本田技研工業株式会社 | 脚式移動ロボットの制御装置 |
KR101687629B1 (ko) | 2010-01-18 | 2016-12-20 | 삼성전자주식회사 | 인간형 로봇 및 그 보행 제어방법 |
US8437993B2 (en) | 2010-08-12 | 2013-05-07 | International Business Machines Corporation | Dynamically adjusting simulation fidelity based on changes in activity levels of components |
US8428921B2 (en) | 2010-08-12 | 2013-04-23 | International Business Machines Corporation | Dynamically adjusting simulation fidelity based on checkpointed fidelity state |
US8468005B2 (en) | 2010-08-12 | 2013-06-18 | International Business Machines Corporation | Determining simulation fidelity in a self-optimized simulation of a complex system |
US8407035B2 (en) * | 2010-08-12 | 2013-03-26 | International Business Machines Corporation | Dynamically adjusting simulation fidelity in a self-optimized simulation of a complex system |
WO2016185877A1 (ja) * | 2015-05-20 | 2016-11-24 | ボッシュ株式会社 | ウェアラブルロボットの制御装置、ウェアラブルロボット、及び、ウェアラブルロボットの制御方法 |
US9586316B1 (en) * | 2015-09-15 | 2017-03-07 | Google Inc. | Determination of robotic step path |
KR102556924B1 (ko) * | 2016-09-05 | 2023-07-18 | 삼성전자주식회사 | 보행 보조 방법 및 이를 수행하는 장치 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002326173A (ja) * | 2001-04-27 | 2002-11-12 | Honda Motor Co Ltd | 脚式移動ロボットの動作生成装置 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3167404B2 (ja) * | 1992-02-26 | 2001-05-21 | 本田技研工業株式会社 | ロボットの関節駆動制御装置 |
JP3148827B2 (ja) | 1992-04-30 | 2001-03-26 | 本田技研工業株式会社 | 脚式移動ロボットの歩行制御装置 |
JP3273443B2 (ja) * | 1992-05-22 | 2002-04-08 | 本田技研工業株式会社 | ロボットのリンクなどの軌道生成方法及び装置 |
US5872893A (en) * | 1996-07-25 | 1999-02-16 | Honda Giken Kogyo Kabushiki Kaisha | Gait generation system of legged mobile robot |
JP3663034B2 (ja) | 1996-07-25 | 2005-06-22 | 本田技研工業株式会社 | 脚式移動ロボットの歩容生成装置 |
JP3629133B2 (ja) | 1997-01-31 | 2005-03-16 | 本田技研工業株式会社 | 脚式移動ロボットの制御装置 |
JP3672406B2 (ja) * | 1997-01-31 | 2005-07-20 | 本田技研工業株式会社 | 脚式移動ロボットの歩容生成装置 |
JP3443077B2 (ja) | 1999-09-20 | 2003-09-02 | ソニー株式会社 | ロボットの運動パターン生成装置及び運動パターン生成方法、並びにロボット |
JP3615702B2 (ja) * | 1999-11-25 | 2005-02-02 | ソニー株式会社 | 脚式移動ロボットの動作制御装置及び動作制御方法、並びに、脚式移動ロボット |
US6876903B2 (en) * | 2000-11-17 | 2005-04-05 | Honda Giken Kogyo Kabushiki Kaisha | Gait pattern generating device for legged mobile robot |
JP3679105B2 (ja) * | 2001-12-28 | 2005-08-03 | 本田技研工業株式会社 | 脚式移動ロボットの歩容生成装置 |
DE60236457D1 (de) * | 2001-12-28 | 2010-07-01 | Honda Motor Co Ltd | Gangerzeugungsvorrichtung für beweglichen roboter mit beinen und steuervorrichtung |
WO2003061917A1 (fr) | 2002-01-18 | 2003-07-31 | Honda Giken Kogyo Kabushiki Kaisha | Dispositif de commande pour robot bipede |
US6999851B2 (en) * | 2002-08-30 | 2006-02-14 | Sony Corporation | Robot apparatus and motion controlling method therefor |
US7236852B2 (en) * | 2002-10-11 | 2007-06-26 | Sony Corporation | Motion controlling apparatus and method and remote controlling apparatus and method for legged mobile robot |
US7054718B2 (en) * | 2002-10-11 | 2006-05-30 | Sony Corporation | Motion editing apparatus and method for legged mobile robot and computer program |
-
2005
- 2005-02-28 EP EP05719663A patent/EP1738879B1/en not_active Not-in-force
- 2005-02-28 KR KR1020067017634A patent/KR101140812B1/ko not_active IP Right Cessation
- 2005-02-28 JP JP2006510510A patent/JP4805817B2/ja not_active Expired - Fee Related
- 2005-02-28 KR KR1020067017772A patent/KR101160161B1/ko active IP Right Grant
- 2005-02-28 WO PCT/JP2005/003345 patent/WO2005082583A1/ja active Application Filing
- 2005-02-28 JP JP2006510511A patent/JP4805818B2/ja not_active Expired - Fee Related
- 2005-02-28 EP EP05719662A patent/EP1736286B1/en not_active Not-in-force
- 2005-02-28 US US10/597,931 patent/US7715945B2/en active Active
- 2005-02-28 WO PCT/JP2005/003346 patent/WO2005082582A1/ja active Application Filing
- 2005-02-28 US US10/597,933 patent/US7774098B2/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002326173A (ja) * | 2001-04-27 | 2002-11-12 | Honda Motor Co Ltd | 脚式移動ロボットの動作生成装置 |
Also Published As
Publication number | Publication date |
---|---|
WO2005082583A1 (ja) | 2005-09-09 |
WO2005082582A1 (ja) | 2005-09-09 |
US20070156284A1 (en) | 2007-07-05 |
EP1738879A1 (en) | 2007-01-03 |
KR20060129423A (ko) | 2006-12-15 |
US20070061038A1 (en) | 2007-03-15 |
JP4805817B2 (ja) | 2011-11-02 |
JPWO2005082583A1 (ja) | 2007-10-25 |
KR101140812B1 (ko) | 2012-05-03 |
EP1738879B1 (en) | 2011-12-21 |
US7715945B2 (en) | 2010-05-11 |
KR101160161B1 (ko) | 2012-06-27 |
EP1736286A4 (en) | 2009-12-09 |
EP1736286B1 (en) | 2011-11-30 |
KR20060126803A (ko) | 2006-12-08 |
US7774098B2 (en) | 2010-08-10 |
EP1738879A4 (en) | 2009-12-09 |
EP1736286A1 (en) | 2006-12-27 |
JPWO2005082582A1 (ja) | 2007-10-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4805818B2 (ja) | 移動ロボットの歩容生成装置 | |
JP5221688B2 (ja) | 脚式移動ロボットの歩容生成装置 | |
JP4800038B2 (ja) | 移動ロボットの歩容生成装置 | |
EP1825967B1 (en) | Gait generator of legged mobile robot | |
JP4199236B2 (ja) | 脚式移動ロボットの歩容生成装置 | |
JP5465074B2 (ja) | 2足移動ロボットの制御装置 | |
JP5468974B2 (ja) | 2足移動ロボットの制御装置及び歩容生成装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20071213 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20101207 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110204 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20110809 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20110811 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4805818 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140819 Year of fee payment: 3 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |