Stair Climbing Stabilization of the HRP-4 Humanoid Robot . Stphane - PowerPoint PPT Presentation

Stair Climbing Stabilization of the HRP-4 Humanoid Robot . Stéphane Caron December 11 and 14, 2018 JRL Seminar, CNRS-AIST Joint Robotics Laboratory, Tsukuba (December 11) Seminar at the Department of Mechano-Informatics, the University of Tokyo (December 14)

context . COMANOID project - https://comanoid.cnrs.fr 1

demonstrator scenario . 2

stair climbing part . https://www.youtube.com/watch?v=vFCFKAunsYM 3

system overview . 4 Whole-body Walking Whole-body Kinematic Pattern Admittance Generation Control Control DCM DCM Control Observer

walking pattern generation . 5 Desired Kinematic Footstep Whole-body Walking Whole-body Targets Locations Kinematic Pattern Admittance Generation Control Control DCM DCM Control Observer

linear inverted pendulum mode m IEEE/RSJ International Conference on Intelligent Robots and Systems . 2001. 1. Shuuji Kajita, Fumio Kanehiro, Kenji Kaneko, Kazuhito Yokoi et Hirohisa Hirukawa. « The 3D 1 . 6 • Equation of motion : • Floating base dynamics : M ¨ q + N = S T τ + J T f ¨ c = 1 ∑ i f i ˙ L c = ∑ i ( p i − c ) × f i • Angular momentum ˙ L c = 0 : c = ω 2 ( c − z ) ¨ with ω 2 = g / h and z the ZMP Linear Inverted Pendulum Mode : A simple modeling for a biped walking pattern generation ». In :

divergent component of motion . rence on Intelligent Robots and Systems . 2011. 2. Johannes Englsberger, Christian Ott, Maximo Roa, Alin Albu-Schäffer, Gerhard Hirzinger et 2 • Unstable dynamics : 7 c • Divergent Component of Motion : • LIPM equation of motion : c = ω 2 ( c − z ) ¨ ξ = c + ˙ ω ˙ ξ = ω ( ξ − z ) al. « Bipedal walking control based on capture point dynamics ». In : IEEE/RSJ International Confe-

walking pattern generation . Generate a CoM-ZMP trajectory that is : Consistent Feasible • ZMP belongs to support area • Contact force within friction cone Viable Not falling. For this system, same as Figure adapted from [Gri+17]. 8 c ( t ) = ω 2 ( c ( t ) − z ( t )) ∀ t > 0 , ¨ bounded : ∃ M > 0 , ∀ t > 0 , ∥ c ( t ) ∥ ≤ M

walking pattern generation . So far we have tested three methods : • Linear Model Predictive Control [Wie06] • Foot-guided Agile Control through ZMP Manipulation [SY17] • Capturability of Variable-Height Inverted Pendulum [Car+18] 9

linear model predictive control . 2006. 3. Pierre-Brice Wieber. « Trajectory free linear model predictive control for stable walking in the 3 • Guarantee of recursive feasibility : [CWF17] • Variable step timings : [BW17] • Variable CoM height : [Bra+15] Allows a number of extensions : • Viability : terminal DCM • Feasibility : ZMP in support area • Consistency : state equation Constraints • Minimize CoM jerk • Track desired CoM velocity • Track desired ZMP reference Cost function Formulate preview control [Kaj+03] as a Quadratic Program (QP) : 10 presence of strong perturbations ». In : IEEE-RAS International Conference on Humanoid Robots .

foot-guided agile control through zmp manipulation . tems . 2017. 4. Tomomichi Sugihara et Takanobu Yamamoto. « Foot-guided Agile Control of a Biped Robot 4 • Call many times to adapt step timings • Analytical solution : • Finite horizon, continuous time dynamics 11 subject to Predictive control with ZMP as input minimize ∫ T ( z ( t ) − z d ) 2 d t ⇒ Feasibility (best effort) z ( t ) 0 c = ω 2 ( c − z ) ¨ ⇒ Consistency c ( T ) + ˙ c ( T ) = ξ d ⇒ Viability ω z ∗ (0) = z d + 2( ξ (0) − z d ) − ( ξ d − z d ) e − ω T 1 − e − 2 ω T through ZMP Manipulation ». In : IEEE/RSJ International Conference on Intelligent Robots and Sys-

variable-height inverted pendulum . ouvertes.fr/hal-01689331/document . 5. Stéphane Caron, Adrien Escande, Leonardo Lanari et Bastien Mallein. « Capturability-based 5 • Call many times to adapt step timings 12 • New input λ > 0 for height variations : ¨ c = λ ( c − z ) + g • Viability ⇒ boundedness condition : ∫ ∞ ( λ ( t ) r ( t ) − g ) e − Ω( t ) d t ξ (0) = 0 • Solve : tailored optimization (30-50 µ s) Analysis, Optimization and Control of 3D Bipedal Walking ». 2018. url : https://hal.archives-

visualization . Visualization of stair climbing pattern in mc_rtc 13

walking stabilization . 14 Desired Kinematic Footstep Whole-body Walking Whole-body Targets Locations Kinematic Pattern Admittance Generation Control Control DCM DCM Control Observer

role of stabilization . Actuated joints converge but unactuated floating base diverges : In walking pattern By robot without stabilization Figure adapted from [Tak+09]. 15

floating base facts . Let's review the facts : • Floating base translation is unactuated • Its dynamics are reduced to : • Only way to control it is via indirect force control of the ZMP z 16 c = ω 2 ( c − z ) ¨

indirect force control . ... but our robot is position-controlled ? Split control into two components : Admittance control Allow position changes to improve force tracking Floating-base control Assuming force control, select reaction force to control the floating base 17

visualization . Standing stabilization under external forces 18

admittance control . 19 Commanded Desired Kinematic Kinematic Commanded Footstep Whole-body Walking Whole-body Targets Targets Joint Angles Locations Kinematic Pattern Admittance Generation Control Control Distributed Measured Foot Wrenches Foot Wrenches DCM DCM Control Observer

strategies . Different strategies for different components of the net contact wrench : • CoP at each contact [Kaj+01b] • Pressure distribution [Kaj+10] • CoM admittance control [Nag99] 20

center-of-pressure control . • Rotate end-effector to move its CoP • Assumes compliance at contact : • Apply damping control : Figure adapted from [Kaj+01b] 6 6. Shuuji Kajita, Kazuhito Yokoi, Muneharu Saigo et Kazuo Tanie. « Balancing a Humanoid Robot Using Backdrive Concerned Torque Control and Direct Angular Momentum Feedback ». In : IEEE International Conference on Robotics and Automation . 2001. 21 τ = K e ( θ − θ e ) ˙ θ = A cop ( τ d − τ ) • Closed-loop behavior has τ → τ d

pressure distribution control . Kaneko, Fumio Kanehiro et Kazuhito Yokoi. « Biped walking stabilization based on linear inverted 7. Shuuji Kajita, Mitsuharu Morisawa, Kanako Miura, Shin'ichiro Nakaoka, Kensuke Harada, Kenji 7 Figure adapted from [Kaj+10] 22 • Apply damping control : • Push down with foot that needs more • Net vertical force compensates pressure, lift the other one gravity ⇒ only need to control : ∆ f z = f Rz − f Lz * d z p * d ctrl p Lz Rz ˙ z ctrl = A z (∆ f zd − ∆ f z ) f f Rz Lz pendulum tracking ». In : IEEE/RSJ International Conference on Intelligent Robots and Systems . 2010.

com admittance control • Closed-loop behavior : analysis is 8. Discussions with Pr T. Sugihara. 9 . only possible with delay or 23 • Accelerate CoM against ZMP error : c = ω 2 ( c − z ) • Recall that ¨ ¨ c = A c ( z − z d ) disturbance observer 8 measured desired 9. Ken'ichiro Nagasaka. « Whole-body Motion Generation for a Humanoid Robot by Dynamics Filters ». In : PhD thesis (1999). The University of Tokyo, in Japanese.

choice of strategies . Which ones to choose ? End-effector strategies • CoP at each contact [Kaj+01b] • Pressure distribution [Kaj+10] ... are sufficient to control the net wrench, yet : CoM admittance control [Nag99] • uses other joints, e.g. hips 10. The effect is similar to Model ZMP Control [Tak+09]. 24 • helps recover from ZMP saturation 10

com admittance in stair climbing . 25 Measured ZMP 0.3 Measured DCM 0.2 Desired DCM 0.1 Sagittal Coordinate (m) 0.0 -0.1 Desized ZMP 0.3 Maximum ZMP 0.2 Minimum ZMP 0.1 0.0 -0.1 Time (s) Figure 1 : Top : no CoM admittance control. Bottom : with A c = 20 [Hz 2 ].

dcm control . 26 Commanded Desired Kinematic Kinematic Commanded Footstep Whole-body Walking Whole-body Targets Targets Joint Angles Locations Kinematic Pattern Admittance Generation Control Control Distributed Measured Foot Wrenches Foot Wrenches Estimated Desired DCM Measured Joint Angles DCM DCM DCM Desired CoM & Contacts Measured IMU Orientation Control Observer

dcm control . mation . 2009. 11. Tomomichi Sugihara. « Standing stabilizability and stepping maneuver in planar bipedalism 11 Compensation [Nag99] • Equivalent to Inverted Pendulum 27 • Apply proportional feedback to it : • yields best CoM-ZMP regulator [Sug09] • requires less control input [Tak+09] c = ω 2 ( c − z ) • Assume control of z in ¨ • Control only the DCM ˙ ξ = ω ( ξ − z ) : ξ = ˙ ˙ ξ d + k p ( ξ d − ξ ) [ ] z = z d − 1 + k p ( ξ d − ξ ) ω based on the best COM-ZMP regulator ». In : IEEE International Conference on Robotics and Auto-