<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Neurorobot.</journal-id>
<journal-title>Frontiers in Neurorobotics</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Neurorobot.</abbrev-journal-title>
<issn pub-type="epub">1662-5218</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fnbot.2017.00049</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Neuroscience</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Human-Derived Disturbance Estimation and Compensation (DEC) Method Lends Itself to a Modular Sensorimotor Control in a Humanoid Robot</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Lippi</surname> <given-names>Vittorio</given-names></name>
<xref ref-type="author-notes" rid="fn001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/278662/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Mergner</surname> <given-names>Thomas</given-names></name>
<uri xlink:href="http://loop.frontiersin.org/people/2823/overview"/>
</contrib>
</contrib-group>
<aff><institution>Neurology, University Clinics of Freiburg</institution> <country>Freiburg, Germany</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Manish Sreenivasa, Heidelberg University, Germany</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Egidio Falotico, Sant&#x00027;Anna School of Advanced Studies, Italy; Malte Schilling, Bielefeld University, Germany</p></fn>
<fn fn-type="corresp" id="fn001"><p>&#x0002A;Correspondence: Vittorio Lippi <email>vittorio.lippi&#x00040;uniklinik-freiburg.de</email></p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>08</day>
<month>09</month>
<year>2017</year>
</pub-date>
<pub-date pub-type="collection">
<year>2017</year>
</pub-date>
<volume>11</volume>
<elocation-id>49</elocation-id>
<history>
<date date-type="received">
<day>31</day>
<month>01</month>
<year>2017</year>
</date>
<date date-type="accepted">
<day>24</day>
<month>08</month>
<year>2017</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2017 Lippi and Mergner.</copyright-statement>
<copyright-year>2017</copyright-year>
<copyright-holder>Lippi and Mergner</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract><p>The high complexity of the human posture and movement control system represents challenges for diagnosis, therapy, and rehabilitation of neurological patients. We envisage that engineering-inspired, model-based approaches will help to deal with the high complexity of the human posture control system. Since the methods of system identification and parameter estimation are limited to systems with only a few DoF, our laboratory proposes a heuristic approach that step-by-step increases complexity when creating a hypothetical human-derived control systems in humanoid robots. This system is then compared with the human control in the same test bed, a posture control laboratory. The human-derived control builds upon the identified disturbance estimation and compensation (DEC) mechanism, whose main principle is to support execution of commanded poses or movements by compensating for external or self-produced disturbances such as gravity effects. In previous robotic implementation, up to 3 interconnected DEC control modules were used in modular control architectures separately for the sagittal plane or the frontal body plane and successfully passed balancing and movement tests. In this study we hypothesized that conflict-free movement coordination between the robot&#x00027;s sagittal and frontal body planes emerges simply from the physical embodiment, not necessarily requiring a full body control. Experiments were performed in the 14 DoF robot Lucy Posturob (i) demonstrating that the mechanical coupling from the robot&#x00027;s body suffices to coordinate the controls in the two planes when the robot produces movements and balancing responses in the intermediate plane, (ii) providing quantitative characterization of the interaction dynamics between body planes including frequency response functions (FRFs), as they are used in human postural control analysis, and (iii) witnessing postural and control stability when all DoFs are challenged together with the emergence of inter-segmental coordination in squatting movements. These findings represent an important step toward controlling in the robot in future more complex sensorimotor functions such as walking.</p></abstract>
<kwd-group>
<kwd>sensory-motor system</kwd>
<kwd>humans</kwd>
<kwd>neuromechanical modeling</kwd>
<kwd>modular control architecture</kwd>
<kwd>humanoid robot experiments</kwd>
</kwd-group>
<contract-num rid="cn001">600698</contract-num>
<contract-num rid="cn001">610454</contract-num>
<contract-sponsor id="cn001">Seventh Framework Programme<named-content content-type="fundref-id">10.13039/501100004963</named-content></contract-sponsor>
<counts>
<fig-count count="11"/>
<table-count count="1"/>
<equation-count count="17"/>
<ref-count count="66"/>
<page-count count="22"/>
<word-count count="5291"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<sec>
<title>Human posture control modeling and neurorobotics</title>
<p>Human postural control attracts considerable interest in healthcare research worldwide for reasons such as &#x0201C;fall of the elderly&#x0201D; and neurological impairments such as ataxia in cerebellar patients or deficient movement control in Parkinson&#x00027;s disease. Developing <italic>model-based</italic> diagnostics as well as therapeutic and rehabilitative interventions is an important aim of this research. Engineering-inspired approaches to model the human sensorimotor control and its failures have a long tradition (e.g., Nashner, <xref ref-type="bibr" rid="B41">1972</xref>; Hajos and Kirchner, <xref ref-type="bibr" rid="B13">1984</xref>; Johansson et al., <xref ref-type="bibr" rid="B21">1988</xref>; Kuo, <xref ref-type="bibr" rid="B22">1995</xref>; Fitzpatrick et al., <xref ref-type="bibr" rid="B9">1996</xref>). These approaches often address <italic>reactive</italic> postural responses to well-controlled external stimuli, which lend themselves to system identification, whereas the exact input for voluntary movements is generally unknown. Two of the following recent modeling approaches were especially influential, yet their clinical application can often prove to be problematic.</p>
<p>One modeling approach used system state estimation based on multi-sensory integration under noise optimization principles. This approach builds on biological textbook knowledge on human sensors and anthropometrics and, based on this knowledge and engineering principles, analyzes the human posture control using noise optimizing principles (van der Kooij et al., <xref ref-type="bibr" rid="B61">2001</xref>; Kuo, <xref ref-type="bibr" rid="B23">2005</xref>). For example, it predicts that, whenever possible, posture control preferentially uses proprioceptive rather than vestibular information, as the vestibular information is the one containing more noise. These models allow general predictions on human preferences for certain sensory environments and may define which sensory deficit in patients tends to increase the danger of falling. However, poor correspondence of these models with human neurophysiology and anatomy restricts its clinical usefulness for diagnosing more specific posture control problems of an individual patient.</p>
<p>The other modeling approach used time series or frequency domain data gathered from posture control experiments in humans to establish the simplest model compatible with known human physiology and anatomy that would allow reproduction of the data with identified model parameters. Best known is the independent channel (IC) model of Peterka (<xref ref-type="bibr" rid="B48">2002</xref>) for single inverted pendulum (SIP) scenarios. It identifies sensory weights depending on the proper selection and interpretation of results using different stimuli and test conditions. From extensions of this model to more complex scenarios it was concluded that available engineering methods are in principle capable of arriving at multi-segmental control models. However, the increasing model complexity and the proliferation of parameters tend to reduce the chance of unequivocally identifying the control parameters in health and disease (Mergner and Peterka, <xref ref-type="bibr" rid="B36">2017</xref>). The problem is aggravated by non-linearities of the human control system (described as detection thresholds in Maurer et al., <xref ref-type="bibr" rid="B31">2006</xref>).</p>
<p>As such, there currently exists a dilemma as to which methodology can be used to establish model-based diagnostics and therapeutic or rehabilitative interventions for patients with impaired postural control. The goal of the present study is to contribute to a <italic>heuristic</italic> solution. By this we mean a practicable solution that proceeds from an established model of the human control in the SIP scenario and uses plausible arguments and steps for its extension to more degrees of freedom (DoFs). To evaluate the appropriateness of these steps, we test their effects in special-purpose humanoid robots with human-inspired anthropometrics, sensors, and actuators. In this way a &#x0201C;real world&#x0201D; challenge is imposed which accounts for noisy and inaccurate sensors and non-ideal actuation and mechanics. For as much correspondence as possible to the human situation, these tests are performed in the same testbed that is also used for the human subjects, i.e., in a human posture control laboratory. As described below, first steps in this heuristic approach have successfully been performed. Point of departure was the &#x0201C;disturbance estimation and compensation&#x0201D; (DEC) model for the SIP scenario, which shares basic similarities with the IC model (Mergner et al., <xref ref-type="bibr" rid="B35">2003</xref>; Maurer et al., <xref ref-type="bibr" rid="B31">2006</xref>; Mergner, <xref ref-type="bibr" rid="B33">2010</xref>). While several extension steps have been described previously (see section Previous and Current Steps in the Heuristic Approach), here we report on an extension to a 14 DoF robot and examine whether a modular architecture consisting of a net of DEC controls in the robot&#x00027;s sagittal plane cooperates in a conflict-free way with a corresponding but independent net of DEC controls in the frontal plane during both disturbance compensation and commanded (&#x0201C;voluntary&#x0201D;) movements.</p>
<p>The following subsections aim to combine state of the art human sensorimotor issues with recent robotics issues dealing with posture control. Here and in later sections, we include brief descriptions on the current state of our bottom-up implementation of a human-like postural control in robots. First, we briefly review the biological basis of the DEC model (section Main Features of the Human DEC Model), then consider related issues in neuroscience and robotics (section Modular Control Issues in Neuroscience and Robotics), and finally present a list of the previously performed steps in our neurorobotics approach and the new steps taken in this study (section Previous and Current Steps in the Heuristic Approach).</p>
</sec>
<sec>
<title>Main features of the human DEC model</title>
<p>Both the DEC model and the clearly simpler IC model can describe results from the same protocols for human balance control experiments. The IC model is a linear model that analytically describes human sway responses to support surface tilt in the frequency domain based on a control by proprioceptive, vestibular, and visual feedback channels (Peterka, <xref ref-type="bibr" rid="B48">2002</xref>). The model allows for identification of important features of the human postural control system, the most important ones being time delays in the order of 100&#x02013;200 ms associated with low loop gain and, as a consequence, soft mechanical compliance, and low energy consumption. It also identifies <italic>sensory reweighting</italic>, meaning that humans adjust their use of sensory information to changes in perturbation amplitude and modality. While the IC model describes this feature by using different sets of control parameters, the DEC model (Figure <xref ref-type="fig" rid="F1">1</xref>) is able to predictively describe the data from various experimental conditions with one set of control parameters. The DEC model contains <italic>synthetic</italic> and <italic>holistic</italic> features-synthetic in the sense that it uses disturbance estimations inspired from studies on multisensory fusions in human self-motion perception (Mergner et al., <xref ref-type="bibr" rid="B34">1997</xref>; Mergner and Rosemeier, <xref ref-type="bibr" rid="B37">1998</xref>), holistic in the sense that it also integrates the control of <italic>movements</italic> in a single structure.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Schema of DEC control of the ankle joint in the single inverted pendulum, SIP, scenario with support surface tilt forward (this model is the basis for the modular control architexture in multi-DoF systems). Four external disturbances (external force <italic>F</italic><sub><italic>ext</italic></sub>, gravitational force <italic>F</italic><sub><italic>g</italic></sub>, foot-space rotation &#x003B1;<sub><italic>FS</italic></sub> and foot-space linear acceleration &#x01E8D;<sub><italic>FS</italic></sub>) give rise to the joint torques (<italic>T</italic><sub><italic>ext</italic></sub>, <italic>T</italic><sub><italic>grav</italic></sub>, <italic>T</italic><sub><italic>prop</italic></sub>, and <italic>T</italic><sub><italic>in</italic>(&#x0003D;<italic>inertial</italic>)</sub>, respectively). Sensors (VEST, vestibular; PROP, proprioceptive; TORQUE) are actually networks that combine signals from a variety of transducers (e.g., 3 VEST sensors result from vestibular canal-otolith intractions; Mergner et al., <xref ref-type="bibr" rid="B38">2009</xref>) to yield information on the following physical quantities: &#x01E8D;<sub><italic>hx</italic></sub>, head acceleration in x (sagittal forward) direction, and &#x003B1;<sub><italic>hs</italic></sub> and <inline-formula><mml:math id="M1"><mml:mover accent="true"><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mo>&#x002D9;</mml:mo></mml:mover></mml:math></inline-formula><sub><italic>hs</italic></sub>, head-in-space angle and angular velocity (dashes indicate proprioceptive coordinate transformation to lower body segments, see Figure <xref ref-type="fig" rid="F2">2</xref>; here they point out that in the SIP scenario these vestibular signals refer to the whole body above the anke joints, includinging the legs, i.e., &#x003B1;<sub><italic>ls</italic></sub> and <inline-formula><mml:math id="M2"><mml:mover accent="true"><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mo>&#x002D9;</mml:mo></mml:mover></mml:math></inline-formula><sub><italic>ls</italic></sub>); &#x003B1;<sub><italic>lf</italic></sub> and <inline-formula><mml:math id="M3"><mml:mover accent="true"><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mo>&#x002D9;</mml:mo></mml:mover></mml:math></inline-formula><sub><italic>lf</italic></sub> represent leg-to-foot angle and angular velocity, respectively. Colored boxes derive disturbance estimates (indicated by hats) from the reported physical quantities. Variable foot-in-space angular velocity (<inline-formula><mml:math id="M4"><mml:mover accent="true"><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mo>&#x002D9;</mml:mo></mml:mover></mml:math></inline-formula><sub><italic>fs</italic></sub>) is time-integrated leading to the estimate <inline-formula><mml:math id="M5"><mml:mover accent="true"><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mo>^</mml:mo></mml:mover></mml:math></inline-formula><sub><italic>fs</italic></sub> of foot-in-space position; by means of this signal the proprioceptive signal &#x003B1;<sub><italic>lf</italic></sub> is being &#x0201C;upgraded&#x0201D; (transformed) into space coordinates (<inline-formula><mml:math id="M6"><mml:msub><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mrow><mml:mi>l</mml:mi><mml:msup><mml:mrow><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:msub></mml:math></inline-formula>). C, neural controler; &#x00394;t, lumped neural time delay; &#x003B1;<sub><italic>ls</italic></sub><italic>!</italic>, desired leg-space angle. Box 1/<italic>mgh</italic> transforms torque into an equivalent of an angle. Passive stiffness and its modulation (in humans achieved by muscle co-contration) is omited here for simplicity (compare Ott et al., <xref ref-type="bibr" rid="B46">2016</xref> as to its role for control stability). In current versions of the DEC concept, the input may be varied in three ways depending on the task. The task shown in this figure is to reach and maintain a given orientation of the supported body segment, here the leg segment (in fixed alignment with upper body), in space (&#x003B1;<sub><italic>ls</italic></sub><italic>!</italic>). In situations where the body COM changes with the body configuration or load distribution, the task refers to the body COM in space (<italic>bs</italic>; &#x003B1;<sub><italic>bs</italic></sub><italic>!</italic>) and requires that the current COM location is taken into account in the control. This has been experimentally tested and modeled for human responses in ankle and hip joints to support surface tilts in the sagittal plane (Hettich et al., <xref ref-type="bibr" rid="B17">2014</xref>). As an intuitive third possibility, one may want to reach and maintain a given joint angle, which in the scenario of this figure would mean to command the leg-foot angle &#x003B1;<sub><italic>lf</italic></sub><italic>!</italic> as input and to neglect estimate <inline-formula><mml:math id="M7"><mml:mover accent="true"><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mo>^</mml:mo></mml:mover></mml:math></inline-formula><sub><italic>fs</italic></sub>.</p></caption>
<graphic xlink:href="fnbot-11-00049-g0001.tif"/>
</fig>
<p>The DEC model extends upon the engineering concept of the servo control by negative feedback (see Wiener, <xref ref-type="bibr" rid="B64">1948</xref>, for the early developments of the concept). In neurology, Merton (<xref ref-type="bibr" rid="B39">1953</xref>) used this concept to explain the role of the muscle stretch reflex for the control of posture and movements. He posits that a PD-controller adjusts the force of the muscles so as to produce the desired pose or movement. This would be achieved through negative feedback from proprioceptive sensors, i.e., by feeding the controller with the difference between the desired and the sensed joint angle. However, later researchers considered this concept to be problematic. One reason was that the biological time delay in the feedback loop does not allow for stable performance when large external disturbances such as gravity require high loop gains (see McIntyre and Bizzi, <xref ref-type="bibr" rid="B32">1993</xref>). The DEC model overcomes this problem by estimating the external disturbances and commanding the servo to produce the extra force required for their compensation. The underlying principle is known in control theory under different names, e.g., &#x0201C;feed forward disturbance correction&#x0201D; (Roffel and Betlem, <xref ref-type="bibr" rid="B52">2006</xref>); similar (Luecke and McGuire, <xref ref-type="bibr" rid="B29">1968</xref>; Zhong et al., <xref ref-type="bibr" rid="B66">2012</xref>); in German consistently &#x0201C;St&#x000F6;rgr&#x000F6;&#x000DF;enaufschaltung;&#x0201D; Bleisteiner et al., <xref ref-type="bibr" rid="B5">1961</xref>). The DEC mechanisms operate <italic>context</italic>- and <italic>intent</italic>-dependently (Mergner, <xref ref-type="bibr" rid="B33">2010</xref>). In the terminology of sensorimotor physiology, we subsume them in contradistinction to &#x0201C;short latency reflexes&#x0201D; under &#x0201C;long latency reflexes,&#x0201D; which are known to be modifiable by higher brain centers such as the cerebral cortex (see Pruszynski and Scott, <xref ref-type="bibr" rid="B51">2012</xref>). Disturbance compensation in the DEC is thought to arise reactively with unforeseen external disturbances as well as in the form of predicted-sensory estimates, issued by higher brain centers with foreseen external or self-produced disturbances. Such a &#x0201C;proactive&#x0201D; compensation is assumed to be advantageous compared to a &#x0201C;reactive&#x0201D; one due to short central time delays and absence of sensory noise, leading to more stable control (see Maurer et al., <xref ref-type="bibr" rid="B31">2006</xref>). Noticeably, several reactive and proactive compensatory actions may arise simultaneously, even during voluntary movements (i.e., the &#x0201C;superposition law&#x0201D; applies despite the non-linarities in the disturbance estimates).</p>
<p>As previously explained in more detail (Maurer et al., <xref ref-type="bibr" rid="B31">2006</xref>), humans use several sensory inputs (proprioceptive, vestibular, visual, haptic contact, torque, and pressure) for posture control. They combine this information in various ways to estimate physical variables such as joint angle and angular velocity from muscle spindle, skin receptor, tendon organ, and torque inputs, as well as head linear and angular acceleration from vestibular otolith and semicircular canal inputs, for example. From the estimated physical variables, the DEC model then derives estimates of four classes of disturbances that may influence the torque acting on the skeletal joints (colored boxes in Figure <xref ref-type="fig" rid="F1">1</xref>): (1) <italic>Rotation</italic> and (2) <italic>translation</italic> of the support (be this a supporting body segment or an external support), (3) <italic>gravity and other field forces</italic>, and (4) <italic>contact forces</italic> (&#x0201C;external torque&#x0201D;). As depicted in Figure <xref ref-type="fig" rid="F1">1</xref> (which refers to the SIP scenario), in the case of balancing, the input of the model is the desired angle of the leg segment (here for SIP) in space &#x003B1;<sub><italic>ls</italic></sub><italic>!</italic>, a signal that is thought to stem from higher brain centers including the cerebral cortex. The cortex builds a coherent and continuously updated sensorimotor model of the body as a whole from many input sources. This is often referred to as &#x0201C;body schema:&#x0201D; a conception from neurology that has found its way into humanoid robotics (Morasso, <xref ref-type="bibr" rid="B40">2013</xref>). In the DEC model, noticeably, both the cortical movement commands from the body schema and the conscious perception of the produced movement almost exclusively reflect kinematics. Signals related to stereotypically occurring forces from gravity, link inertia, inter-link coupling forces, etc. are hidden, so to speak, because they are compensated by the DEC mechanism (details in Mergner, <xref ref-type="bibr" rid="B33">2010</xref>).</p>
<p>In this form, the DEC concept qualified for a modular control architecture in scenarios with &#x0003E;1 DoF (see Section Previous and Current Steps in the Heuristic Approach). Modular control as a concept has attracted considerable interest earlier in neuroscience and robotics. The following section points out differences to the present concept.</p>
</sec>
<sec>
<title>Modular control issues in neuroscience and robotics</title>
<p>In the DEC concept, a modular control by two DEC modules has been used for the modeling of human ankle joint and hip joint responses to support surface tilt and for mimicking the responses in a robot (Hettich et al., <xref ref-type="bibr" rid="B17">2014</xref>). Further work supported this idea of a modular control, for example when we added a third DEC control in a robot for proactive squatting in the knee joints during reactive balancing (Ott et al., <xref ref-type="bibr" rid="B46">2016</xref> and <ext-link ext-link-type="uri" xlink:href="https://www.youtube.com/watch?v=3ALCTMW3Ei4">https://www.youtube.com/watch?v=3ALCTMW3Ei4</ext-link>).</p>
<p>Generally, modularity is often applied at the level of kinematics, as is the case in neuroscience with &#x0201C;movement synergies,&#x0201D; where each control module commands a subset of the controlled DoFs. Specifically, a kinematic synergy is often defined as a function mapping a scalar value to several DoFs. Such synergies are often used to reduce the dimensionality of the control problem with the result that the number of controlled synergies is smaller than the total number of DoFs. The principle has not only inspired research of human walking (e.g., Ivanenko et al., <xref ref-type="bibr" rid="B20">2003</xref>), but also various fields of humanoid control (Hauser et al., <xref ref-type="bibr" rid="B15">2007</xref>). Noticeably, however, the number of required synergies in complex movements may reach the number of DoFs. This may happen for example when movement control is organized in terms of eigenmovements (Alexandrov et al., <xref ref-type="bibr" rid="B1">2017</xref>) where the aim is to free the control from coupling forces stemming from other movements. Certain synergies studied in humans, such as torque synergies, are so far largely neglected in robotics. Another related concept is that of motor primitives, originally meaning basic kinematic, dynamic, or muscular building blocks of movements arising at neuronal levels (Flash and Hochner, <xref ref-type="bibr" rid="B10">2005</xref>). Generally, motor primitives can be combined in several ways, i.e., to act simultaneously by summing them or to obtain superposition of effects, as is the case with synergies, or serial effects when using them in a sequence, as with predefined trajectories or velocity profiles. The concept of motor primitives has found extensive application in robotics for learning of motor tasks (Schaal et al., <xref ref-type="bibr" rid="B53">2003</xref>). Some implementations of modularity are not only advantageous in that they reduce complexity in control design, but also in increasing control robustness (in case one module fails, remaining modules can take over).</p>
<p>Compared to these modularity concepts, the architecture of combining DEC modules is clearly distinct. Here, each DoF of the human skeletal system is controlled by one DEC module. This even applies to situations where two or more modules do the same job and conceptually may be viewed as one module (e.g., during sagittal body sway about the two ankle joints with aligned axes and the body weight equally distributed on both feet). Each module determines the torque to be applied to the controlled DoF. The desired trajectory is specified as an input to each module. All modules have essentially the same structure and have no internal model of the whole system. Yet, the modules operate not completely independently of each other, because they exchange sensory information through coordinate transformation across the joints that interconnect the body segments (Figure <xref ref-type="fig" rid="F2">2A</xref>). Compared to the above examples of modular control, the DEC model can be defined as a low level control system that takes care of the fundamental task of posture control for undisturbed motor execution and acts at the level of joint kinematics. Coordination between different joints may emerge from the interaction between the modules and the body mechanics under task, even if no kinematic synergy is explicitly specified. For example, hip-ankle coordination in the experiments of Hettich et al. (<xref ref-type="bibr" rid="B17">2014</xref>) emerged from the tasks for the ankle and hip controllers to maintain the body COM over the feet as base of support and the trunk upright, respectively.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p><bold>(A)</bold> Schematic representation of the interconnections between the DEC control modules, in terms of down-channeling and up-channeling of sensory information in a given body plane (here sagittal plane; compare Figure <xref ref-type="fig" rid="F4">4</xref> for analogies in the frontal plane). Shown by solid downward arrows is the distribution of the vestibular signal &#x0201C;head angle with respect to the gravitational vertical&#x0201D; (&#x003B1;<sub><italic>hs</italic></sub>). Analog distributions hold for the angular velocity signal <inline-formula><mml:math id="M8"><mml:mover accent="true"><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mo>&#x002D9;</mml:mo></mml:mover></mml:math></inline-formula><sub><italic>hs</italic></sub> and the head linear acceleration signal in x (sagittal forward) direction, &#x01E8D;<sub><italic>Hx</italic></sub>. The gray upward arrows show the up-channeling of processed signals, representing the causal physical interactions in a stack of superimposed body segments. For example, up channeling of <inline-formula><mml:math id="M9"><mml:mover accent="true"><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mo>^</mml:mo></mml:mover></mml:math></inline-formula><sub><italic>fs</italic></sub> from the ankle module to the knee module and fusing it with the proprioceptive &#x003B1;<sub><italic>lf</italic></sub> signal may be used in the knee module as a predictor of &#x003B1;<sub><italic>hs</italic></sub> (a mechanism that bypases the sensory feedback chain and increases control stability). This up-channeling is especially effective if the support surface is stationary, whereby <inline-formula><mml:math id="M10"><mml:mover accent="true"><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mo>^</mml:mo></mml:mover></mml:math></inline-formula><sub><italic>fs</italic></sub> becomes subthreshold (through a velocity threshold) and with it the effect of the particularly noisy vestibular <inline-formula><mml:math id="M11"><mml:mover accent="true"><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mo>&#x002D9;</mml:mo></mml:mover></mml:math></inline-formula><sub><italic>hs</italic></sub> signal; see Mergner et al., <xref ref-type="bibr" rid="B38">2009</xref>). <bold>(B)</bold> Picture of Lucy Posturob (from right and behind) balancing on the motion platform in the human posture control laboratory.</p></caption>
<graphic xlink:href="fnbot-11-00049-g0002.tif"/>
</fig>
<p>Further aspects will be considered in section Discussion.</p>
</sec>
<sec>
<title>Previous and current steps in the heuristic approach</title>
<p>In previous works, the DEC model was subjected to &#x0201C;real world&#x0201D; tests, separately for the sagittal and the frontal plane. The tests were performed with humanoid robots equipped with human-inspired sensors, actuation, and anthropomorphic properties in the human posture control laboratory. Robot responses to support surface tilt in the sagittal plane were successfully compared to human responses once in a SIP scenario (Mergner et al., <xref ref-type="bibr" rid="B38">2009</xref>) and later in a DIP scenario (Hettich et al., <xref ref-type="bibr" rid="B17">2014</xref>). The latter study also showed that it sufficed to use interconnected DEC control modules for the ankle and the hip joints in order to simulate the human responses in the robot (for simplification, the joints on both sides were mechanically coupled). Interestingly, the ankle-hip coordination emerged from the interconnection of the two modules (see section Modular Control Issues in Neuroscience and Robotics). Subsequent experiments focused on a transfer of DEC to the frontal plane (14 DoF robot; Lippi et al., <xref ref-type="bibr" rid="B28">2016</xref>) and the use of the knee joints to test the control of squatting movements in the sagittal plane. To this end, realization of the control in Simulink (The MathWorks Inc., Natick, USA) and using force controlled actuation allowed us to transfer the DEC control to the DLR robot TORO for comparison with a fully model-based control (Ott et al., <xref ref-type="bibr" rid="B46">2016</xref> and <ext-link ext-link-type="uri" xlink:href="https://www.youtube.com/watch?v=3ALCTMW3Ei4">https://www.youtube.com/watch?v=3ALCTMW3Ei4</ext-link>).</p>
<p>This study investigates whether the networks of DEC modules previously controlling the robots&#x00027; sagittal plane and frontal plane separately would adequately cooperate when combined during movements in intermediate planes. Specifically, we asked whether the physical linkages given by the robot&#x00027;s body suffice to guarantee adequate cooperation between the sensorimotor controls, or whether a supervising full body model or some other form of software linkage between the two networks would be required. The experiments are performed in a robot with 14 DoFs (Lucy Posturob, Figure <xref ref-type="fig" rid="F2">2B</xref>) using one and the same set of control parameters for non-trivial scenarios that draw on the scalability of the DEC modular control. The tests included both voluntary movements and balancing responses to passive body motions in an intermediate plane. Postural and control stability across all DoF of the robot were tested when the robot performed squatting movements in the knee joints, evoking the emergence of inter-segmental coordination in the ankle and hip joints. Further experiments aimed to provide a quantitative characterization of the robot&#x00027;s balancing responses in terms of FRFs (as they are used in human postural control analyses) and tested postural stability when increasing frequency of voluntary movements simultaneously in the sagittal and frontal body planes up to the performance limits.</p>
</sec>
</sec>
<sec id="s2">
<title>Methods&#x02014;robotic platform and experimental procedures</title>
<sec>
<title>Humanoid robot Lucy Posturob</title>
<p>The Posturob III robot, <italic>Lucy Posturob</italic> (details in Figure <xref ref-type="fig" rid="F3">3</xref>), was constructed as a humanoid of 1.5 m body height and &#x0007E;20 kg body weight with anthropometric parameters inspired by Winter (<xref ref-type="bibr" rid="B65">1990</xref>). Its body consists of the upper body (HAT, for head, arms, and trunk), pelvis, the two thighs, shanks (lower leg), and foot segments, all made of aluminum and interconnected by hinge joints. The total of 14 DoF (Figures <xref ref-type="fig" rid="F3">3A,B</xref>) comprise 2 DoF per ankle joint, 1 per knee, 3 per hip, and 2 DoF for the lower vertebral column (&#x0201C;lumped&#x0201D; DoFs across the vertebrae of the back bone).</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Details of the humanoid robot Lucy. <bold>(A)</bold> Scheme of the 14 degrees of freedom. <bold>(B)</bold> Picture of Lucy Posturob (arms represented by weights; vestibular system is located in head; no visual system is currently provided by the eyes; electronics for actuation is contained in the pelvis). <bold>(C)</bold> Details of foot with actuation (a, DC motor; b, spindle; c, enconder; d, load cell as torque sensor). Actuation is the same in all DoFs.</p></caption>
<graphic xlink:href="fnbot-11-00049-g0003.tif"/>
</fig>
<p>Using an analog acquisition board, a PC read from mechatronic sensors the signals joint torque, joint angular position and velocity, and the anterior-posterior and medio-lateral pressure distribution under each foot (not used in the present experiments). The same computer also read via USB the signals from a custom-made human-inspired artificial vestibular system (Mergner et al., <xref ref-type="bibr" rid="B38">2009</xref>). The DEC control model was implemented and executed as a compiled Simulink model (Real-Time Windows Target, The MathWorks Inc., Natick, USA). The control system worked at 200 Hz. The delay of the system as a whole is estimated to be 20 ms (which corresponds approximately to the value to which we adjusted all sensory signals for correct time-matched interactions in the disturbance estimations). This delay is, admittedly, shorter that the lumped delays identified in humans (see section Main Features of the Human DEC Model). However, additional features have to be added to the DEC control modules in the Lucy robot before we can aim to repeat some of our experiments with more human-inspired time delays (as in Hettich et al., <xref ref-type="bibr" rid="B17">2014</xref>).</p>
<p>The vestibular sensor was fixed to the robot&#x00027;s HAT segment. A bio-inspired algorithm (Mergner et al., <xref ref-type="bibr" rid="B38">2009</xref>) fuses the inputs from 3 accelerometers and 3 gyrometers to estimate (i) HAT angular velocity in space, (ii) HAT angle with respect to the gravitational vertical, and (iii) linear acceleration of the upper HAT end representing head position (compare Figure <xref ref-type="fig" rid="F3">3A</xref>). Joint actuation was achieved by DC electric motors (part <italic>a</italic> in Figure <xref ref-type="fig" rid="F3">3C</xref>) that rigidly interconnected the two links by a screw/spindle system (part <italic>b</italic>). The spindle drive transformed the rotational movement of the motors into a lever movement, as measured by a linear potentiometer, part c, for producing joint angle. An inner torque control loop was implemented on the on-board robot electronics, receiving the torque command signal from the higher-level bio-inspired DEC algorithms on the PC. The force sensor, part d in Figure <xref ref-type="fig" rid="F3">3C</xref>, measured the tangential force acting on the joint in order to compute the joint torque signal used (i) in the DEC control loop and (ii) in the on-board control of the torque.</p>
</sec>
<sec>
<title>Disturbance estimation in multiple DoF system</title>
<p>The four estimated disturbances in a DEC control module are not directly available as sensory inputs and are hence reconstructed through inter-sensory interactions. The sensory inputs used in these and our previous robot experiments with the DEC model are from the vestibular system, joint proprioception and joint torque (compare Figure <xref ref-type="fig" rid="F1">1</xref>). Currently unconsidered are visual self-motion cues for balancing (see Assl&#x000E4;nder et al., <xref ref-type="bibr" rid="B3">2015</xref>), and foot pressure cues.</p>
<p>In the multi-segment human body, the sensor fusions for the disturbance estimates use signals from both local sensors (proprioceptive joint angle and angular velocity sensors and torque sensors) and remote (e.g., vestibular) sensors. The vestibular signals in Lucy are conveyed to the control modules through down-channeling from their source in the head, as schematically shown in Figure <xref ref-type="fig" rid="F2">2</xref>. Additional information is required for the disturbance estimations concerning mass and inertia distribution of the supported body segments with respect to the supporting joint, as described for the SIP scenario in Mergner (<xref ref-type="bibr" rid="B33">2010</xref>). To further account for momentary changes of these parameters in the multi-segment system, Lippi et al. (<xref ref-type="bibr" rid="B27">2013</xref>) provided a general description of the down-channeling of the processed sensory information. In particular, for the two disturbance estimators of support rotation and gravitational torque, <inline-formula><mml:math id="M12"><mml:mover accent="true"><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mo>^</mml:mo></mml:mover></mml:math></inline-formula><sub><italic>fs</italic></sub> and <inline-formula><mml:math id="M13"><mml:mover accent="true"><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mo>^</mml:mo></mml:mover></mml:math></inline-formula><sub><italic>grav</italic></sub> (see Figure <xref ref-type="fig" rid="F1">1</xref>), the COM location of all segments above a given supporting joint is calculated step-wise downwards. This allows for treating the COMs of all supported segments in a given joint as if they were the COM of a single-segment body. For example, with sagittal support surface tilts the hip joints compensate the gravity impact from any upper body lean, while the ankle joints are compensating the gravity effect due to the lean of the whole body COM. This example also shows that, while the action of a control module is local, the estimated disturbances represent global effects acting on the body. For the <italic>support surface translation estimates</italic> <inline-formula><mml:math id="M14"><mml:mover accent="true"><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mo>^</mml:mo></mml:mover></mml:math></inline-formula><sub><italic>in</italic></sub>, the simplification by down-channeling concerns the joint torque produced by the combined inertial forces exerted by all upper body segments, while for the <italic>contact force estimates</italic> <inline-formula><mml:math id="M15"><mml:mover accent="true"><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mo>^</mml:mo></mml:mover></mml:math></inline-formula><sub><italic>ext</italic></sub> it is the torque produced by combined forces having impact on the supported segments. The sensor fusions used to reconstruct these global variables are distributed among modules and are based on signal exchanges between modules.</p>
<p>The present study explains the generalized concepts of disturbance estimations with down and up channeling between modules (Figure <xref ref-type="fig" rid="F2">2</xref>) in more detail. It is worth noting that an advantage of using up-channeling is that the information of the physical variables conveyed upward are already processed at a lower level. This is especially relevant for the estimate of support surface tilt in the ankle joint. Running the input velocity signal through a threshold reduces noise of this estimate (see below). The up-channeling of this signal is used in upper segments, instead of the local input signals in these segments, which carry more noise. This up-channeled foot(support)-in-space signal contributes to the computation of the variables controlled by the servo loops of the higher modules. Module inputs and outputs are shown in Table <xref ref-type="table" rid="T1">1</xref>. Included are also the inputs coming directly from the sensory system (e.g., joint angle). These sensor fusions are described here in more detail than before, using a generalized notation that allows for an arbitrary number of DoFs.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>Signals and parameters used in the sensor fusion process.</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th valign="top" align="left"><bold>Signal</bold></th>
<th valign="top" align="left"><bold>Symbol</bold></th>
<th valign="top" align="left"><bold>Description</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Body in space sway, up-channeled</td>
<td valign="top" align="left"><inline-formula><mml:math id="M16"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>b</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x0005E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula></td>
<td valign="top" align="left">Angle sway of the COM of all the segments above the controlled link, obtained using the upchanneled <inline-formula><mml:math id="M17"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x0005E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> signal. It can be used in the servo loop as controlled variable</td>
</tr>
<tr>
<td valign="top" align="left">Body in space sway, down-channeled</td>
<td valign="top" align="left"><inline-formula><mml:math id="M18"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>b</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula></td>
<td valign="top" align="left">Angle sway of the COM of all the segments above the controlled link, obtained using the down-channeled <inline-formula><mml:math id="M19"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>signal. It is used to compute external disturbances.</td>
</tr>
<tr>
<td valign="top" align="left">Desired value for the controlled variable</td>
<td valign="top" align="left"><inline-formula><mml:math id="M20"><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>f</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:mo>!</mml:mo></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M21"><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>b</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:mo>!</mml:mo></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M22"><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:mo>!</mml:mo></mml:mrow></mml:math></inline-formula></td>
<td valign="top" align="left">Reference for joint angle (lf), body COM orientation in space (bs) and link orientation in space (ls)</td>
</tr>
<tr>
<td valign="top" align="left">Proprioceptive input</td>
<td valign="top" align="left"><inline-formula><mml:math id="M23"><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>f</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow></mml:math></inline-formula></td>
<td valign="top" align="left">The angular position of the controlled joint. Used to compute <inline-formula><mml:math id="M24"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M25"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x0005E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>. It can be used in the servo loop as controlled variable.</td>
</tr>
<tr>
<td valign="top" align="left">Foothold link orientation in space, downchanneled</td>
<td valign="top" align="left"><inline-formula><mml:math id="M26"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula></td>
<td valign="top" align="left">Orientation in space of the link supporting the controlled joint, it is used to compute <inline-formula><mml:math id="M27"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>b</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>and it is passed to the underlying module as <inline-formula><mml:math id="M28"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>.</td>
</tr>
<tr>
<td valign="top" align="left">Foothold link orientation in space, up-channeled</td>
<td valign="top" align="left"><inline-formula><mml:math id="M29"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x0005E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula></td>
<td valign="top" align="left">In the module controlling the support joint it is computed using <sup>&#x02323;</sup><inline-formula><mml:math id="M30"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mover accent='true'><mml:mi>&#x003B1;</mml:mi><mml:mo>&#x002D9;</mml:mo></mml:mover><mml:mrow><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> as shown in Equation (1). In the other modules it is up-channeled from the underlying module as <inline-formula><mml:math id="M31"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x0005E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>.</td>
</tr>
<tr>
<td valign="top" align="left">Controlled link in space rotation speed, down-channeled.</td>
<td valign="top" align="left"><inline-formula><mml:math id="M32"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mover accent='true'><mml:mi>&#x003B1;</mml:mi><mml:mo>&#x002D9;</mml:mo></mml:mover><mml:mrow><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula></td>
<td valign="top" align="left">It is down-channeled and used for the computation of <inline-formula><mml:math id="M33"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mn>0</mml:mn></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x0005E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>.</td>
</tr>
<tr>
<td>&#x000A0;</td>
<td>&#x000A0;</td>
<td>&#x000A0;</td>
</tr>
<tr>
<td valign="top" align="left">Controlled link in space orientation down-channeled.</td>
<td valign="top" align="left"><inline-formula><mml:math id="M34"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x0005E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula></td>
<td valign="top" align="left">Angle sway of the COM of all the segments above the controlled link. It is used to compute external disturbances.</td>
</tr>
<tr>
<td valign="top" align="left">Controlled link in space orientation up-channeled.</td>
<td valign="top" align="left"><inline-formula><mml:math id="M35"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x0005E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula></td>
<td valign="top" align="left">It can be used in the servo loop as controlled variable</td>
</tr>
<tr>
<td valign="top" align="left">Center of mass position in the controlled plane</td>
<td valign="top" align="left"><inline-formula><mml:math id="M36"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:mstyle mathvariant='bold'><mml:mi>C</mml:mi><mml:mi>O</mml:mi><mml:msup><mml:mi>M</mml:mi><mml:mi>n</mml:mi></mml:msup></mml:mstyle></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula></td>
<td valign="top" align="left">Position in space of the center of mass of all the links above the <italic>n</italic> th joint. It is down-channeled so that in each module it can be updated to take in account the controlled body segment.</td>
</tr>
<tr>
<td valign="top" align="left">Mass of all the segments above the controlled joint</td>
<td valign="top" align="left"><inline-formula><mml:math id="M37"><mml:msubsup><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula></td>
<td valign="top" align="left">This parameter is down-channeled. Each module is updating it adding the mass of the controlled link, that is an internal parameter</td>
</tr>
<tr>
<td valign="top" align="left">Moment of inertia of all the segments above the controlled joint</td>
<td valign="top" align="left"><inline-formula><mml:math id="M38"><mml:msubsup><mml:mrow><mml:mi>J</mml:mi></mml:mrow><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula></td>
<td valign="top" align="left">The moment of inertia is down-channeled and updated in each module on the basis of the internal parameters describing the controlled link and the configuration of the body.</td>
</tr>
</tbody>
</table>
</table-wrap>
<sec>
<title>Support surface tilt</title>
<p>In the SIP scenario of Figure <xref ref-type="fig" rid="F1">1</xref>, an estimate of support surface tilt is obtained by combining a vestibular derived leg-in-space signal and a local ankle proprioceptive leg-on-foot signal. Specifically, the signal used to compensate the support surface rotation is derived from the rotation speed of the foot in space, when the foot is in firm contact with the support surface <inline-formula><mml:math id="M39"><mml:mover accent="true"><mml:mrow><mml:msubsup><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mrow><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msubsup></mml:mrow><mml:mo>&#x002D9;</mml:mo></mml:mover></mml:math></inline-formula>.</p>
<disp-formula id="E1"><label>(1)</label><mml:math id="M40"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mn>1</mml:mn></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x0005E;</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mstyle displaystyle='true'><mml:mrow><mml:msubsup><mml:mo>&#x0222B;</mml:mo><mml:mn>0</mml:mn><mml:mi>t</mml:mi></mml:msubsup><mml:mrow><mml:mi>&#x003C1;</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:msubsup><mml:mover accent='true'><mml:mi>&#x003B1;</mml:mi><mml:mo>&#x0005E;</mml:mo></mml:mover><mml:mrow><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:msup><mml:mrow></mml:mrow><mml:mrow><mml:msup><mml:mrow></mml:mrow><mml:mo>&#x000B7;</mml:mo></mml:msup></mml:mrow></mml:msup><mml:mn>1</mml:mn></mml:mrow></mml:msubsup><mml:mo stretchy='false'>)</mml:mo></mml:mrow></mml:mrow></mml:mstyle><mml:mo>&#x02212;</mml:mo><mml:mi>k</mml:mi><mml:msubsup><mml:mover accent='true'><mml:mi>&#x003B1;</mml:mi><mml:mo>&#x0005E;</mml:mo></mml:mover><mml:mrow><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mn>1</mml:mn></mml:msubsup><mml:mtext>&#x000A0;</mml:mtext><mml:mi>d</mml:mi><mml:mi>&#x003C4;</mml:mi></mml:mrow></mml:math></disp-formula>
<p>In the generalized notation used, the pedix <italic>fs</italic> stands for <italic>foothold in space</italic>, i.e., the orientation in space of the link under the controlled joint, which represents an extension of the variable <italic>foot in space</italic>. Similarly, in the following the variables <italic>ls</italic> for <italic>link in space</italic> and <italic>lf</italic> for <italic>link to foothold</italic> denote the generalization of the variables <italic>leg in space</italic> and <italic>leg to foot</italic> as used in the SIP and the DIP (double inverted pendulum) case before (Mergner, <xref ref-type="bibr" rid="B33">2010</xref>; Hettich et al., <xref ref-type="bibr" rid="B16">2013</xref>; Lippi et al., <xref ref-type="bibr" rid="B27">2013</xref>). This notation now allows for a description of the sensor fusion process with a generic number of modules, where the index <italic>n</italic> represents the position of the controlled joint (<italic>n</italic> &#x0003D; 1 is the ankle, <italic>n</italic> &#x0003D; 2 the knee, etc.). For specific cases, a notation referring to the names of body segments can be used for simplicity, as shown in Figure <xref ref-type="fig" rid="F2">2A</xref> where <italic>ts</italic> for example means <italic>thigh in space</italic>. All the signals in (1) are estimates, with the hat denoting here the up-channeled estimate, &#x003C1;(&#x000B7;) is the threshold function and <italic>t</italic> the current time. The integration is leaky (modulated by the term <italic>k</italic>).</p>
<p>Given the threshold &#x00398; &#x0003E;0, the function &#x003C1;(&#x000B7;) is defined as</p>
<disp-formula id="E2"><label>(2)</label><mml:math id="M41"><mml:mrow><mml:mi>&#x003C1;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x003B1;</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:mi>&#x003B1;</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mi>&#x003B8;</mml:mi><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mi>&#x003B1;</mml:mi><mml:mo>&#x02264;</mml:mo><mml:mo>&#x02212;</mml:mo><mml:mi>&#x003B8;</mml:mi></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mo>&#x02212;</mml:mo><mml:mi>&#x003B8;</mml:mi><mml:mo>&#x0003C;</mml:mo><mml:mi>&#x003B1;</mml:mi><mml:mo>&#x0003C;</mml:mo><mml:mi>&#x003B8;</mml:mi><mml:mtext>&#x000A0;&#x000A0;</mml:mtext></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mi>&#x003B1;</mml:mi><mml:mo>&#x02212;</mml:mo><mml:mi>&#x003B8;</mml:mi><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mi>&#x003B1;</mml:mi><mml:mo>&#x02265;</mml:mo><mml:mi>&#x003B8;</mml:mi></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mrow></mml:mrow></mml:math></disp-formula>
<p>The value of &#x00398; is a parameter of the control module. The presence of the threshold function &#x003C1;(&#x000B7;) introduces a gain non-linearity in the balance behavior (i.e., larger support surface tilts are more compensated than smaller ones), as previously observed in human experiments (e.g., Hettich et al., <xref ref-type="bibr" rid="B17">2014</xref>).</p>
<p>The estimate of support surface tilt from Equation (1) is used to reconstruct the orientation in space of the link above the controlled joint <inline-formula><mml:math id="M42"><mml:mover accent="true"><mml:mrow><mml:msubsup><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo>^</mml:mo></mml:mover></mml:math></inline-formula> inside each module. The signal is up-channeled and fused with the proprioceptive input throughout all the modules of a given body plane. The resulting value in the nth module is</p>
<disp-formula id="E3"><label>(3)</label><mml:math id="M43"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x0005E;</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mn>1</mml:mn></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x0005E;</mml:mo></mml:mover><mml:mo>&#x0002B;</mml:mo><mml:mstyle displaystyle='true'><mml:msubsup><mml:mo>&#x02211;</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>f</mml:mi></mml:mrow><mml:mi>k</mml:mi></mml:msubsup></mml:mrow></mml:mstyle></mml:mrow></mml:math></disp-formula>
<p>This is different with the other three sensory disturbance estimates (below), in that these are based on the down-channeled position. This scheme reproduces human-like responses in robot experiments and model simulations. However, the interactions between up-channeled and down-channeled signals in the general case are still under research.</p>
</sec>
<sec>
<title>Field forces such as gravity</title>
<p>The gravity torque <italic>T</italic><sub><italic>g</italic></sub> in the general case is calculated by</p>
<disp-formula id="E4"><label>(4)</label><mml:math id="M44"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>g</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msubsup><mml:mi>m</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:mi>g</mml:mi><mml:mi>C</mml:mi><mml:mi>O</mml:mi><mml:msubsup><mml:mi>M</mml:mi><mml:mi>x</mml:mi><mml:mi>n</mml:mi></mml:msubsup></mml:mrow></mml:math></disp-formula>
<p>where <inline-formula><mml:math id="M45"><mml:mi>C</mml:mi><mml:mi>O</mml:mi><mml:msubsup><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> is the horizontal component of the position of the center of mass <inline-formula><mml:math id="M46"><mml:msubsup><mml:mrow><mml:mstyle mathvariant='bold-italic' mathsize='normal'><mml:mi>COM</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mtext>&#x000A0;</mml:mtext></mml:mrow><mml:mrow><mml:mstyle class="text"><mml:mtext mathvariant="bold">n</mml:mtext></mml:mstyle></mml:mrow></mml:msubsup></mml:math></inline-formula> of all the segments above the controlled joint. The estimated <inline-formula><mml:math id="M47"><mml:msubsup><mml:mrow><mml:mstyle mathvariant='bold-italic' mathsize='normal'><mml:mi>COM</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mtext>&#x000A0;</mml:mtext></mml:mrow><mml:mrow><mml:mstyle class="text"><mml:mtext mathvariant="bold">n</mml:mtext></mml:mstyle></mml:mrow></mml:msubsup></mml:math></inline-formula> is computed performing the weighted sum</p>
<disp-formula id="E5"><label>(5)</label><mml:math id="M48"><mml:mrow><mml:msup><mml:mrow><mml:mover accent='true'><mml:mrow><mml:mstyle mathvariant='bold-italic' mathsize='normal'><mml:mi>COM</mml:mi></mml:mstyle></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>n</mml:mi></mml:mstyle></mml:msup><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mover accent='true'><mml:mrow><mml:mstyle mathvariant='bold-italic' mathsize='normal'><mml:mi>COM</mml:mi></mml:mstyle></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>n</mml:mi></mml:mstyle><mml:mo>&#x0002B;</mml:mo><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mn>1</mml:mn></mml:mstyle></mml:mrow></mml:msup><mml:mo>&#x0002B;</mml:mo><mml:msup><mml:mi>L</mml:mi><mml:mi>n</mml:mi></mml:msup><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>cos</mml:mi></mml:mstyle><mml:mo stretchy='false'>(</mml:mo><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mo stretchy='false'>)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>sin</mml:mi></mml:mstyle><mml:mo stretchy='false'>(</mml:mo><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mo stretchy='false'>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:msubsup><mml:mi>m</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:msup><mml:mi>h</mml:mi><mml:mi>n</mml:mi></mml:msup><mml:mi>cos</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mo stretchy='false'>)</mml:mo><mml:msubsup><mml:mi>m</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mrow><mml:msubsup><mml:mi>m</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow></mml:mfrac></mml:mrow></mml:math></disp-formula>
<p>where <italic>L</italic><sup><italic>n</italic></sup> is the length of the link controlled by the joint and <italic>h</italic><sup><italic>n</italic></sup> is the distance of the COM of the nth link from the nth joint, <inline-formula><mml:math id="M49"><mml:msubsup><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mtext>&#x000A0;</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> is the mass of the nth link, and <inline-formula><mml:math id="M50"><mml:msubsup><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> the total mass of all the links above the nth joint. The inverted hat denotes estimators based on down-channeled signals. The modules can be set to use <inline-formula><mml:math id="M51"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mrow><mml:mi>C</mml:mi><mml:mi>O</mml:mi><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> or to apply a small angle approximation as done in previous experiments (Mergner, <xref ref-type="bibr" rid="B33">2010</xref>; Hettich et al., <xref ref-type="bibr" rid="B17">2014</xref>). The estimate then becomes</p>
<disp-formula id="E6"><label>(6)</label><mml:math id="M52"><mml:mrow><mml:msub><mml:mrow><mml:mover accent='true'><mml:mi>T</mml:mi><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mi>g</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msubsup><mml:mi>m</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:mi>g</mml:mi><mml:msup><mml:mover accent='true'><mml:mi>h</mml:mi><mml:mo>&#x002DC;</mml:mo></mml:mover><mml:mi>n</mml:mi></mml:msup><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>b</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:mrow></mml:math></disp-formula>
<p>where the expression <inline-formula><mml:math id="M53"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>b</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> &#x0003D; <italic>atan</italic>2(<inline-formula><mml:math id="M54"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mrow><mml:mi>C</mml:mi><mml:mi>O</mml:mi><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M55"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mrow><mml:mi>C</mml:mi><mml:mi>O</mml:mi><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>) represents the angular sway in space of the center of the body mass above the nth joint, while <inline-formula><mml:math id="M56"><mml:msup><mml:mover accent='true'><mml:mi>h</mml:mi><mml:mo>&#x002DC;</mml:mo></mml:mover><mml:mi>n</mml:mi></mml:msup></mml:math></inline-formula> is the average height of the COM of all the segments above the controlled joint.</p>
</sec>
<sec>
<title>Support surface linear acceleration</title>
<p>In the presence of support surface acceleration described in the reference system of the support, an inertial force on the center of mass of the body arises. The external acceleration is computed for each joint. The part of the vestibular head acceleration signal not explained by trunk rotation at the hip or at any joint below is taken to stem from support surface acceleration. This is expressed as</p>
<disp-formula id="E7"><label>(7)</label><mml:math id="M57"><mml:mrow><mml:msub><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>E</mml:mi><mml:mi>X</mml:mi><mml:mi>T</mml:mi><mml:mi>E</mml:mi><mml:mi>R</mml:mi><mml:mi>N</mml:mi><mml:mi>A</mml:mi><mml:mi>L</mml:mi></mml:mrow></mml:msub><mml:mtext>&#x000A0;</mml:mtext><mml:mo>=</mml:mo><mml:mtext>&#x000A0;&#x000A0;</mml:mtext><mml:msub><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>V</mml:mi><mml:mi>E</mml:mi><mml:mi>S</mml:mi><mml:mi>T</mml:mi><mml:mi>I</mml:mi><mml:mi>B</mml:mi><mml:mi>U</mml:mi><mml:mi>L</mml:mi><mml:mi>A</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:msub><mml:mo>&#x02212;</mml:mo><mml:msub><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>E</mml:mi><mml:mi>L</mml:mi><mml:mi>F</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></disp-formula>
<p>where the acceleration produced by the joint movements is:</p>
<disp-formula id="E8"><label>(8)</label><mml:math id="M58"><mml:mrow><mml:msubsup><mml:mrow><mml:mover accent='true'><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>a</mml:mi></mml:mstyle><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>S</mml:mi><mml:mi>E</mml:mi><mml:mi>L</mml:mi><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:msubsup><mml:mrow><mml:mover accent='true'><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>a</mml:mi></mml:mstyle><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>S</mml:mi><mml:mi>E</mml:mi><mml:mi>L</mml:mi><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mi>L</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mfrac><mml:mrow><mml:msup><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:msup><mml:mi>t</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:mi>s</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:msubsup><mml:mrow><mml:mover accent='true'><mml:mi>&#x003B1;</mml:mi><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup><mml:mo stretchy='false'>)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mi>c</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:msubsup><mml:mrow><mml:mover accent='true'><mml:mi>&#x003B1;</mml:mi><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup><mml:mo stretchy='false'>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow></mml:math></disp-formula>
<p>The disturbance <italic>inertial torque</italic> then results from</p>
<disp-formula id="E9"><label>(9)</label><mml:math id="M59"><mml:mrow><mml:msub><mml:mrow><mml:mover accent='true'><mml:mi>T</mml:mi><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:msub><mml:mrow><mml:mover accent='true'><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>a</mml:mi></mml:mstyle><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>E</mml:mi><mml:mi>X</mml:mi><mml:mi>T</mml:mi><mml:mi>E</mml:mi><mml:mi>R</mml:mi><mml:mi>N</mml:mi><mml:mi>A</mml:mi><mml:mi>L</mml:mi></mml:mrow></mml:msub><mml:msubsup><mml:mrow><mml:mover accent='true'><mml:mrow><mml:mstyle mathvariant='bold-italic' mathsize='normal'><mml:mi>COM</mml:mi></mml:mstyle></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mtext>&#x000A0;</mml:mtext><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>n</mml:mi></mml:mstyle></mml:msubsup><mml:msubsup><mml:mi>m</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow></mml:math></disp-formula>
<p>For simplicity, here only the horizontal translation is considered</p>
<disp-formula id="E10"><label>(10)</label><mml:math id="M60"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mover accent='true'><mml:mi>x</mml:mi><mml:mo>&#x000A8;</mml:mo></mml:mover><mml:mrow><mml:mi>E</mml:mi><mml:mi>X</mml:mi><mml:mi>T</mml:mi><mml:mi>E</mml:mi><mml:mi>R</mml:mi><mml:mi>N</mml:mi><mml:mi>A</mml:mi><mml:mi>L</mml:mi></mml:mrow></mml:msub><mml:mi>C</mml:mi><mml:mi>O</mml:mi><mml:msubsup><mml:mi>M</mml:mi><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>y</mml:mi></mml:mstyle><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>n</mml:mi></mml:mstyle></mml:msubsup><mml:msubsup><mml:mi>m</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow></mml:math></disp-formula>
</sec>
<sec>
<title>Contact force disturbance (such as a push or pull)</title>
<p>Humans may sense amount and location of a force exerting impact on the body directly (locally), and in addition may sense the impact that the force has on their balancing in terms of the evoked center of pressure (COP) shifts under the feet, which is proportional to ankle torque. Recall the SIP case scenario of Figure <xref ref-type="fig" rid="F1">1</xref>, where the external torque <italic>T</italic><sub><italic>ext</italic></sub> affects the sensory measure of the actively produced torque (<italic>T</italic><sub><italic>a</italic></sub>). It contributes together with the gravitational torque, inertial torque, passive torque (from muscle and connective tissues) and total torque (<italic>T</italic><sub><italic>g</italic></sub>, <italic>T</italic><sub><italic>in</italic></sub>, <italic>T</italic><sub><italic>p</italic></sub>, and <italic>T</italic><sub><italic>A</italic></sub>, respectively) to <italic>T</italic><sub><italic>ext</italic></sub>, which thus is given by</p>
<disp-formula id="E11"><label>(11)</label><mml:math id="M61"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mtext>ext</mml:mtext></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:mo>&#x02212;</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mi>g</mml:mi></mml:msub><mml:mo>&#x02212;</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x02212;</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mi>p</mml:mi></mml:msub><mml:mo>&#x02212;</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mi>a</mml:mi></mml:msub></mml:mrow></mml:math></disp-formula>
<p>The term <italic>T</italic><sub><italic>A</italic></sub> is computed using the down-channeled signal <inline-formula><mml:math id="M62"><mml:mover accent="false"><mml:mrow><mml:msubsup><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mrow><mml:mi>b</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo>&#x02323;</mml:mo></mml:mover></mml:math></inline-formula> and, <inline-formula><mml:math id="M63"><mml:msubsup><mml:mrow><mml:mi>J</mml:mi></mml:mrow><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>, the moment of inertia of the supported body segments as</p>
<disp-formula id="E12"><label>(12)</label><mml:math id="M64"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mi>d</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>b</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mtext>&#x000A0;</mml:mtext><mml:msubsup><mml:mi>J</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:math></disp-formula>
<p>This expression takes in account that also <inline-formula><mml:math id="M65"><mml:msubsup><mml:mrow><mml:mi>J</mml:mi></mml:mrow><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> may change with respect to time due to the movements of other DoF. Similarly to <inline-formula><mml:math id="M66"><mml:msubsup><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>, also <inline-formula><mml:math id="M67"><mml:msubsup><mml:mrow><mml:mi>J</mml:mi></mml:mrow><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> is computed down-channeling information through the modules. For simplicity the down-channeled value is the moment of inertia computed as</p>
<disp-formula id="E13"><label>(13)</label><mml:math id="M68"><mml:mtable columnalign='left'><mml:mtr><mml:mtd><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>J</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='false'>&#x0007E;</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>J</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup></mml:mrow><mml:mo stretchy='false'>&#x0007E;</mml:mo></mml:mover><mml:mo>&#x0002B;</mml:mo><mml:msubsup><mml:mi>m</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:msup><mml:mrow><mml:mo>&#x02016;</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mover accent='true'><mml:mrow><mml:mstyle mathvariant='bold-italic' mathsize='normal'><mml:mi>COM</mml:mi></mml:mstyle></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>n</mml:mi></mml:mstyle><mml:mo>&#x0002B;</mml:mo><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mn>1</mml:mn></mml:mstyle></mml:mrow></mml:msup><mml:mo>&#x02212;</mml:mo><mml:msup><mml:mrow><mml:mover accent='true'><mml:mrow><mml:mstyle mathvariant='bold-italic' mathsize='normal'><mml:mi>COM</mml:mi></mml:mstyle></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>n</mml:mi></mml:mstyle></mml:msup></mml:mrow><mml:mo>&#x02016;</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup><mml:mtext>&#x000A0;</mml:mtext><mml:mo>&#x0002B;</mml:mo><mml:msup><mml:mi>J</mml:mi><mml:mi>n</mml:mi></mml:msup></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mo>&#x0002B;</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:msup><mml:mi>m</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msup><mml:msup><mml:mrow><mml:mo>&#x02016;</mml:mo><mml:mrow><mml:mstyle mathvariant='bold-italic' mathsize='normal'><mml:mi>CO</mml:mi></mml:mstyle><mml:msubsup><mml:mstyle mathvariant='bold-italic' mathsize='normal'><mml:mi>M</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>l</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mi>k</mml:mi></mml:mstyle></mml:mrow><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>n</mml:mi></mml:mstyle></mml:msubsup><mml:mo>&#x02212;</mml:mo><mml:msup><mml:mrow><mml:mover accent='true'><mml:mrow><mml:mstyle mathvariant='bold-italic' mathsize='normal'><mml:mi>COM</mml:mi></mml:mstyle></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>n</mml:mi></mml:mstyle></mml:msup></mml:mrow><mml:mo>&#x02016;</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>Where <italic>J</italic><sup><italic>n</italic></sup> and <inline-formula><mml:math id="M69"><mml:msubsup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>COM</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>link</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>n</mml:mtext></mml:mstyle></mml:mrow></mml:msubsup></mml:math></inline-formula> are, respectively, the moment of inertia and the position of the COM of the nth link. The moment of inertia used in Equation (12) is</p>
<disp-formula id="E14"><label>(14)</label><mml:math id="M70"><mml:mrow><mml:msubsup><mml:mi>J</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>J</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:mo>&#x0002B;</mml:mo></mml:mrow><mml:mo stretchy='false'>&#x002DC;</mml:mo></mml:mover><mml:msubsup><mml:mi>m</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:msup><mml:mrow><mml:mrow><mml:mo>&#x02016;</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mover accent='true'><mml:mrow><mml:mstyle mathvariant='bold-italic' mathsize='normal'><mml:mi>COM</mml:mi></mml:mstyle></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>n</mml:mi></mml:mstyle></mml:msup></mml:mrow><mml:mo>&#x02016;</mml:mo></mml:mrow></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:math></disp-formula>
<p>Compensating the estimated <italic>T</italic><sub><italic>ext</italic></sub> may imply that, in a system with neural time delays, positive feedback in <italic>T</italic><sub><italic>a</italic></sub> requires a limitation of the compensation (gain &#x0003C; 1, low-pass filtering), similarly to the case of the translation estimator. Model simulations suggest that humans may deal with these flaws by transiently increasing passive stiffness, which has zero time delay, through co-contraction of antagonistic muscles pairs.</p>
<sec>
<title>Servo loop and compensation of external disturbances</title>
<p>The disturbance compensation is implemented by summing the disturbance estimates with the input of the controller (PD) with negative sign for compensation (compare Figure <xref ref-type="fig" rid="F1">1</xref>). In the case of the support surface tilt this can be seen as a coordinate transformation of the controlled variable from joint coordinates to space coordinates. The other disturbances are normalized by <inline-formula><mml:math id="M71"><mml:msubsup><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup><mml:mi>g</mml:mi><mml:msup><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>h</mml:mi></mml:mrow><mml:mo>&#x0007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula>, which represents an angle equivalent to the torque, i.e., the angle that would produce the torque evoked by gravity during body (or segment) lean in linear approximation. This makes the input of the PD controller homogeneous. The torque commanded by the servo controller in the <italic>n</italic>th module is defined by</p>
<disp-formula id="E15"><label>(15)</label><mml:math id="M72"><mml:mtable columnalign='left'><mml:mtr><mml:mtd><mml:msub><mml:mi>T</mml:mi><mml:mi>a</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mi>p</mml:mi></mml:msub><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mi>&#x003B5;</mml:mi><mml:mo>&#x02212;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>g</mml:mi></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>x</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>/</mml:mo><mml:msubsup><mml:mi>m</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:mi>g</mml:mi><mml:msup><mml:mover accent='true'><mml:mi>h</mml:mi><mml:mo>&#x002DC;</mml:mo></mml:mover><mml:mi>n</mml:mi></mml:msup></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mi>d</mml:mi></mml:msub><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mover accent='true'><mml:mi>&#x003B5;</mml:mi><mml:mo>&#x002D9;</mml:mo></mml:mover><mml:mo>&#x02212;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mover accent='true'><mml:mi>T</mml:mi><mml:mo>&#x002D9;</mml:mo></mml:mover><mml:mi>g</mml:mi></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mover accent='true'><mml:mi>T</mml:mi><mml:mo>&#x002D9;</mml:mo></mml:mover><mml:mrow><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mover accent='true'><mml:mi>T</mml:mi><mml:mo>&#x002D9;</mml:mo></mml:mover><mml:mrow><mml:mi>e</mml:mi><mml:mi>x</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>/</mml:mo><mml:msubsup><mml:mi>m</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>p</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:mi>g</mml:mi><mml:msup><mml:mover accent='true'><mml:mi>h</mml:mi><mml:mo>&#x002DC;</mml:mo></mml:mover><mml:mi>n</mml:mi></mml:msup></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>K</italic><sub><italic>p</italic></sub> and <italic>K</italic><sub><italic>d</italic></sub> are the proportional and the derivative gain, respectively, and &#x003B5; is the error of the controlled variable as computed using the up-channeled information</p>
<disp-formula id="E16"><label>(16)</label><mml:math id="M73"><mml:mrow><mml:mtext>&#x000A0;</mml:mtext><mml:mi>&#x003B5;</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:mo>!</mml:mo><mml:mo>&#x02212;</mml:mo><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x0005E;</mml:mo></mml:mover></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mi>i</mml:mi><mml:mi>f</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>t</mml:mi><mml:mi>h</mml:mi><mml:mi>e</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>c</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi><mml:mi>t</mml:mi><mml:mi>r</mml:mi><mml:mi>o</mml:mi><mml:mi>l</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi><mml:mi>d</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>v</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mi>i</mml:mi><mml:mi>a</mml:mi><mml:mi>b</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:mtext>&#x000A0;</mml:mtext></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mtext>&#x000A0;</mml:mtext><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>b</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:mo>!</mml:mo><mml:mo>&#x02212;</mml:mo><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>b</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x002DC;</mml:mo></mml:mover></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mi>i</mml:mi><mml:mi>f</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>t</mml:mi><mml:mi>h</mml:mi><mml:mi>e</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>c</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi><mml:mi>t</mml:mi><mml:mi>r</mml:mi><mml:mi>o</mml:mi><mml:mi>l</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi><mml:mi>d</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>v</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mi>i</mml:mi><mml:mi>a</mml:mi><mml:mi>b</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>b</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>f</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:mo>!</mml:mo><mml:mo>&#x02212;</mml:mo><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>f</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mi>i</mml:mi><mml:mi>f</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>t</mml:mi><mml:mi>h</mml:mi><mml:mi>e</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>c</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi><mml:mi>t</mml:mi><mml:mi>r</mml:mi><mml:mi>o</mml:mi><mml:mi>l</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi><mml:mi>d</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>v</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mi>i</mml:mi><mml:mi>a</mml:mi><mml:mi>b</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>f</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mrow><mml:mtext>&#x000A0;</mml:mtext></mml:mrow></mml:math></disp-formula>
<p>where <inline-formula><mml:math id="M74"><mml:mover accent="false"><mml:mrow><mml:msubsup><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mrow><mml:mi>b</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo>&#x0007E;</mml:mo></mml:mover></mml:math></inline-formula> is an estimate of the COM sway of all the links above the controlled joint, as computed using the up-channeled variable <inline-formula><mml:math id="M75"><mml:mover accent="false"><mml:mrow><mml:msubsup><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mrow><mml:mi>l</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo>^</mml:mo></mml:mover></mml:math></inline-formula>.</p>
<p>The effect of each disturbance input and of the error signal of the servo loop can be adjusted by gains (specified as control module parameters). In contrast, the relation between proportional (<italic>K</italic><sub><italic>p</italic></sub>) and derivative (<italic>K</italic><sub><italic>d</italic></sub>) gain is fixed for all disturbances. Gravity compensation with lasting body lean represents a special case. Modeling of human responses to <italic>slow</italic> support surface tilts yielded better results when using a PID controller instead of a PD controller (Mergner et al., <xref ref-type="bibr" rid="B35">2003</xref>; Peterka, <xref ref-type="bibr" rid="B49">2003</xref>) For force feedback as an alternative for the I, see Peterka (<xref ref-type="bibr" rid="B50">2009</xref>). The solution used here is a gain elevation in the gravity estimator for low frequencies (Schweigart and Mergner, <xref ref-type="bibr" rid="B56">2008</xref>). As an alternative, Ott et al. (<xref ref-type="bibr" rid="B46">2016</xref>) used a PID controller for the servo in robot experiments and in addition a specific PD controller for each disturbance estimator.</p>
<p>The DEC implementation used in the robot experiments of Ott et al. (<xref ref-type="bibr" rid="B46">2016</xref>) included feedback from passive joint stiffness and damping. In humans, passive stiffness and damping amount to &#x0007E;10% of the active stiffness and damping, stemming mainly from connective tissue properties of muscles and tendons. Having impact with virtual zero latency, they improve control stability in face of the considerable time delays of the reflexive loops (lumped delay &#x02265;100 ms for the ankle joint). Their implementation in the robot experiments of Ott et al. (<xref ref-type="bibr" rid="B46">2016</xref>) helped in stabilizing the control, as was the case in our previous experiments with Posturob II (Hettich et al., <xref ref-type="bibr" rid="B17">2014</xref>). Implementation of passive stiffness and damping in the Lucy robot is still pending. Further, the implementation of predicted-sensory estimates with pro-active movements, whose effectiveness has been shown before in Posturob I (Mergner, <xref ref-type="bibr" rid="B33">2010</xref>), awaits implementation in Lucy.</p>
</sec>
</sec>
</sec>
<sec>
<title>DEC control in the frontal plane and combined in the sagittal and frontal planes</title>
<p>Conceptually, one DEC control module can be used to control both legs during stance control in the sagittal plane, as realized in our robot experiments in Posturob I and II (Mergner et al., <xref ref-type="bibr" rid="B38">2009</xref>; Hettich et al., <xref ref-type="bibr" rid="B17">2014</xref>). This applies when the two legs and feet are aligned in parallel such that the rotations in the two ankle joints occur approximately around a common axis (for special conditions in which humans weight the proprioceptive input into the controls of the two legs independently of each other, see Pasma et al., <xref ref-type="bibr" rid="B47">2012</xref>). Analogous simplifications hold for the hip joints when extending the control to upper body rotations, and to the knee joints when considering vertical body movements. A mechanically special situation (&#x0201C;four bar linkage&#x0201D;) with only one DoF is given with biped balancing in the frontal plane. In this study, stance width with approximately parallel legs is exclusively considered (Figure <xref ref-type="fig" rid="F3">3A</xref>). In particular, the torque produced in the frontal plane around the ankle joints is relatively small, but the legs can overall produce a large torque on the support surface and hence on the body COM for body stabilization and movements.</p>
<p>Inspired by a model interpretation of human balancing in the frontal plane (Goodworth and Peterka, <xref ref-type="bibr" rid="B11">2010</xref>), Lippi et al. (<xref ref-type="bibr" rid="B28">2016</xref>) suggested for the Lucy robot a preliminary generalization of the DEC control to the frontal plane, formalizing the body kinematics as a double inverted pendulum. The four bar linkage system with one degree of freedom for the lower body control applies if the knee movements are negligible (true with moderate disturbances) and when the feet are continuously kept in contact with the support surface. Then, a simplified model in terms of the SIP scenario applies to lower body sways about a virtual ankle joint between the two actual feet and can be connected by a virtual link to the pelvis joint (Figure <xref ref-type="fig" rid="F4">4A</xref>). Body balancing is then achieved in terms of controlling the position of the COM of the whole body in the frontal plane by applying the appropriate torque to the support surface and maintaining the vertical orientation of the HAT segment above the pelvis. The desired torque for the virtual joint is distributed in the robot on the four actuated joints (ankles and hips). Using the control equations of the previous section, two DEC modules suffice for the control in the frontal plane, for simplification, one for the lower body and one for the upper body.</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p><bold>(A)</bold> Balance control in the frontal plane modeled as a double inverted pendulum. <bold>(B)</bold> For balancing during the squatting experiments (see section Experimental Procedures and Testbed (5), Voluntary Squatting Movements and Figure <xref ref-type="fig" rid="F10">10</xref>), the control is extended to employ 6 DEC control modules that control 10 mechanical DOFs in the frontal and sagittal planes. Note that horizontal hip control is passive proportional-derivative (PD). This full body posture control can, in principle, be used to start gait. Low loop gain provides mechanical compliance (and keeps energy consumption low).</p></caption>
<graphic xlink:href="fnbot-11-00049-g0004.tif"/>
</fig>
<p>The frontal and sagittal planes are mechanically coupled in the physical robot, as are the intermediate planes between them. Theoretically, the interaction between the planes could be considered problematic in terms of coupling forces if the two planes are controlled independently of each other. <italic>Our hypothesis in this study for the control system of Lucy is that, in face of the low loop gain for the disturbance compensations used in the control, the coupling between the dynamics in the two planes is not critical and should allow us to control the two planes independently of each other</italic>. Independent control in the frontal and sagittal planes based on kinematic synergies has been successfully applied for posture control in a small position-controlled humanoid in earlier studies (Hauser et al., <xref ref-type="bibr" rid="B15">2007</xref>).</p>
<p>In the experiments reported below, Lucy was freely standing while actively controlling its balance in the sagittal and frontal planes using the two sets of DEC control modules: one set for the sagittal plane and the other set for frontal plane, while the controls of the horizontal-plane in the hip joints were &#x0201C;passive&#x0201D; through local joint angle proportional and derivative (PD) &#x0201C;proprioceptive&#x0201D; feedback (indicated in red in Figure <xref ref-type="fig" rid="F4">4B</xref>).</p>
</sec>
<sec>
<title>Experimental procedures and testbed</title>
<p>Performance of the DEC system was investigated experimentally when it controlled the sensorimotor behavior of the humanoid robot during balancing of upright stance while compensating external disturbances and self-produced disturbances arising from voluntary movements. Voluntary movements were produced by defining a reference trajectory for the specific control variables (see below). During the experiments (Figure <xref ref-type="fig" rid="F2">2B</xref>), the robot was standing on a 6 DoF motion platform (Stewart platform) in a human posture control laboratory (Mergner et al., <xref ref-type="bibr" rid="B35">2003</xref>). Motions of the robot were recorded by capturing motion signals from the robot&#x00027;s internal sensors (same as used for its control) for comparison between desired and actual poses or movements. The responses to external disturbances were recorded using an external opto-electronic device (Optotrak 3020; Northern Digital Inc.; Waterloo, Canada).</p>
<p>List of experiments performed:
<list list-type="order">
<list-item><p><italic>Test</italic>: Balancing of COM during upright stance on periodically tilting support surface (SS)&#x02014;this tested steady state balancing performance.</p>
<p><italic>Disturbance:</italic> Sinusoidal SS tilts of peak-to-peak (pp) 2&#x000B0; at 0.1 and 0.2 Hz were applied in three trials: (a) sagittal plane, (b) frontal plane, and (c) 45&#x000B0; with respect to the sagittal and frontal planes.</p>
<p><italic>Aim:</italic> Demonstration of stable balancing in all three planes, with cooperative effects for <italic>c</italic> from the simultaneous balancing in the sagittal and frontal planes.</p></list-item>
<list-item><p><italic>Test:</italic> Balancing of COM during upright stance as in (1), but with transient tilts of the SS. This experiment allowed us to distinguish between static and dynamic balancing performance.</p>
<p><italic>Disturbance:</italic> SS tilts with raised cosine velocity profile (dominant frequency, 0.2 Hz) of peak-to-peak 4&#x000B0; were applied again in three trials: (a) in the sagittal plan, (b) in the frontal plane, and (c) 45&#x000B0; with respect to frontal and sagittal planes.</p>
<p><italic>Aim:</italic> Evaluation of static and dynamic balancing performances in the sagittal and frontal planes and their cooperative effects occurring in the intermediate plane.</p></list-item>
<list-item><p><italic>Test:</italic> Balancing of COM during upright stance while pseudo-randomly tilting the SS using the PRTS stimulus<sup>1</sup>. This stimulus allowed us to describe the balancing behavior in terms of FRFs.</p>
<p><italic>Disturbance:</italic> Support surface tilts of peak-to-peak (pp) 1&#x000B0;, 2&#x000B0;, 4&#x000B0;, and 8&#x000B0;.</p>
<p><italic>Disturbance waveform:</italic> PRTS<xref ref-type="fn" rid="fn0001"><sup>1</sup></xref> (0.016&#x02013;2.2 Hz), applied again in the sagittal, the frontal, and the 45&#x000B0; intermediate planes.</p>
<p><italic>Analysis:</italic> Spectral analysis of the angular excursions of the body COM in space.</p>
<p><italic>Aim:</italic> Demonstration of stable balancing in the sagittal and frontal planes across a broad spectrum of tilt frequencies and the cooperative effects occurring between the balancing in the two planes. Varying stimulus amplitude would allow us to test for the human-like non-linearity of responses (expected from the velocity threshold in the ankle module, see Equations 1 and 2).</p></list-item>
<list-item><p><italic>Test:</italic> Voluntary rapid full body movements in the intermediate sagittal-frontal plane.</p>
<p><italic>Movement command:</italic> Starting from a leaning body COM position in an intermediate plane with 45&#x000B0; orientation with respect to sagittal and frontal planes, a fast voluntary movement was commanded to reorient the COM into the vertical position above the feet, using step function references for all the commanded DoFs. The movement was repeated six times to observe the variability of the response. A similar movement was also performed in the frontal plane and in the sagittal plane separately. This allowed us to observe how the coupling effects are affecting the dynamic response in the two planes.</p>
<p><italic>Aim:</italic> Demonstrating proactive movements and testing the robustness of the system in face of strong self-produced disturbances including coupling effects between different joints and between the two controlled planes.</p></list-item>
<list-item><p><italic>Test:</italic> Squatting movements (knee bending).</p>
<p><italic>Movement and task commands:</italic> Raised sinusoids with 4&#x000B0; amplitude at 0.17 Hz were used for commanding knee-bending in repetitive cycles from and back to straight. With the instructed tasks of COM balancing in the ankle and hip joints, the commanded knee-bending was associated with reactive compensatory movements of the whole body in the ankle joints and of the upper body in the hip joints. The task of the pelvis-HAT joint was to maintain a vertical HAT orientation in space.</p>
<p><italic>Aim:</italic> This test challenged the robot&#x00027;s postural stabilization in that all DoFs in the sagittal and frontal planes were interacting. Demonstrating postural stability with the &#x0201C;emerging&#x0201D; movements of the whole body in the ankle joints and of the upper body in the hip joints were secondary aims of the test. For corresponding experiments restricted to the sagittal plane in a 3 DoF DEC implementation in the TORO robot of DLR, see: <ext-link ext-link-type="uri" xlink:href="https://www.youtube.com/watch?v=3ALCTMW3Ei4">https://www.youtube.com/watch?v=3ALCTMW3Ei4</ext-link>).</p></list-item>
<list-item><p><italic>Test:</italic> Voluntary full body movement at increasing speed.</p>
<p><italic>Movement and task commands:</italic> Voluntary body sway movements were commanded simultaneously in the sagittal and frontal planes (the result was a combined movement in some intermediate plane). The common reference signal followed a sinusoidal function with linearly increasing frequency (i.e., <italic>chirp</italic> signal). The amplitude of the reference signal was set to 1&#x000B0; in the sagittal plane and to 2&#x000B0; in the frontal plane. The frequency range was &#x0007E;0.2&#x02013;0.7 Hz.</p>
<p><italic>Aim:</italic> This experiment tests the frequency limits of the active body sway in the sagittal and frontal planes of the ankle joints. While the previous tests were performed <italic>within</italic> the margins of the system&#x00027;s stability (to characterize the normal behavior of the system), this test pushes the robot beyond these margins to performance limits.</p></list-item>
</list></p>
</sec>
</sec>
<sec id="s3">
<title>Experimental results</title>
<p>In the robot experiments, we investigated how the DEC control mechanisms in the sagittal and frontal body planes interact through mechanical coupling of the robot&#x00027;s body in producing reactive responses to external disturbances applied in the 45&#x000B0; intermediate plane (sections Responses to Sinusoidal Support Surface Tilts; Responses to Raised Cosine Support Surface Tilts; Responses to PRTS Support Surface Tilts). A further experiment investigated this interaction for rapid voluntary full body lean movements in the intermediate plane, which were generated by commanding combined action in the sagittal and frontal planes (section Commanded Fast Full-Body Movements in the Intermediate Plane). A fifth experiment (section Voluntary Squatting Movements) tested control stability across all DoFs of Lucy during commanded squatting.</p>
<sec>
<title>Responses to sinusoidal support surface tilts</title>
<p>Lucy&#x00027;s steady state postural responses in terms of body COM sway evoked by the sinusoidal support surface tilts are shown in Figure <xref ref-type="fig" rid="F5">5</xref>. The responses to the &#x000B1;2&#x000B0; tilts in the sagittal and frontal plane at 0.2 Hz (Figures <xref ref-type="fig" rid="F5">5A<sub>1</sub>,B<sub>1</sub></xref>) and 0.1 Hz (Figure <xref ref-type="fig" rid="F5">5A<sub>2</sub>,B<sub>2</sub></xref>) are compared with corresponding responses in the intermediate 45&#x000B0; plane (Figures <xref ref-type="fig" rid="F5">5C<sub>1</sub>,C<sub>2</sub></xref>). Note that compensation is similarly stable in all three stimulus planes, with some residual COM lean resulting in the direction of the tilt. This under-compensation is &#x0201C;human-like,&#x0201D; stemming mainly from imperfect support surface tilt estimations in the DEC modules with gain &#x0003C; 1 and velocity threshold (see notches around maxima and minima).</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>Lucy&#x00027;s sway responses are balancing its body COM during sinusoidal support surface tilt stimuli of &#x000B1;2&#x000B0; about horizontal. The stimuli were applied at stimulus frequencies of 0.2 Hz and 0.1 Hz in the sagittal plane <bold>(A<sub>1</sub>,A<sub>2</sub>)</bold>, frontal plane <bold>(B<sub>1</sub>,B<sub>2</sub>)</bold> and in an intermediate (45&#x000B0; diagonal) plane <bold>(C<sub>1</sub>,C<sub>2</sub>)</bold>. All sway responses show some under-compensation (body slightly sways in direction of support surface tilt). The responses show slow fluctuations, stemming mainly from vestibular noise (but never showed lasting drifts that might have led to loss of balance). See text for further details.</p></caption>
<graphic xlink:href="fnbot-11-00049-g0005.tif"/>
</fig>
</sec>
<sec>
<title>Responses to raised cosine support surface tilts</title>
<p>Figure <xref ref-type="fig" rid="F6">6</xref> shows Lucy&#x00027;s COM transient and static sway responses in the sagittal, frontal, and intermediate plane to support surface tilt stimuli with raised cosine velocity profile and amplitude of 4&#x000B0; with respect to the horizontal. Dynamic and static response components reflect under-compensation of the tilt stimulus similarly as observed in Figure <xref ref-type="fig" rid="F6">6</xref>. The responses in the intermediate plane approximately reflect the sum of the responses in the sagittal and frontal planes. The finding of larger dynamic responses in the sagittal as compared to the frontal plane mainly reflects a stronger effect from body inertia. We can further note an absence of static error when the support surface has moved back to the horizontal position, owing to properties of the support surface tilt estimator with velocity threshold and leaky integrator (see Equation 1).</p>
<fig id="F6" position="float">
<label>Figure 6</label>
<caption><p>Lucy is balancing its body COM during support surface tilt stimuli with &#x0201C;raised cosine velocity&#x0201D; waveform (4&#x000B0; excursions from horizontal). The tilt was applied in the sagittal plane <bold>(A)</bold>, frontal plane <bold>(B)</bold>, and in the intermediate diagonal plane (45&#x000B0; degrees from the frontal) <bold>(C)</bold>. The robot&#x00027;s responses are slightly under-compensating the disturbances (compare Figure <xref ref-type="fig" rid="F5">5</xref>). (Waveform resembles that of most human voluntary movements; it is derived from a bell-shaped velocity profile, <italic>v</italic>(t) &#x0003D; &#x02013;<italic>A</italic> &#x000B7; <italic>f</italic> &#x000B7; cos(2&#x003C0;<italic>ft</italic>) &#x0002B; <italic>A</italic> &#x000B7; <italic>f</italic>, where <italic>t</italic> is time, <italic>A</italic> is angular displacement, and <italic>f</italic> is dominant frequency).</p></caption>
<graphic xlink:href="fnbot-11-00049-g0006.tif"/>
</fig>
</sec>
<sec>
<title>Responses to PRTS support surface tilts</title>
<p>Testing the robot with the PRTS stimulus allowed us to more comprehensively characterize the frequency and amplitude behavior of the system. Being composed of several velocity step functions, this stimulus has a power spectrum with significant components over the whole range of the frequencies that are of interest here to describe the system dynamics [compare section Experimental Procedures and Testbed (3)]. The stimulus was applied with different amplitudes (pp 1, 2, 4, and 8&#x000B0;) for the support surface tilts in the robot&#x00027;s sagittal, frontal, and intermediate body plane. Lucy&#x00027;s PRTS responses in terms of time series data are shown in Figure <xref ref-type="fig" rid="F7">7</xref>. They again show under-compensation as in the previous experiments with the sine and raised cosine stimuli. Here, the under-compensation exhibits a non-linearity upon increase in stimulus amplitude. The non-linearity stands out better in Figure <xref ref-type="fig" rid="F8">8</xref> where the responses are expressed in the upper panels of the corresponding FRFs as sway (error) gain (zero with full compensation and unity if body motion equals platform motion). Note that gain in Figures <xref ref-type="fig" rid="F8">8A,B</xref> decreases with increasing stimulus amplitude, being lowest with the pp 8&#x000B0; stimulus&#x02014;again as can be expected from the velocity threshold contained in the support surface tilt estimation. The basic feature in terms of gain, phase, and coherence resemble each other across the tests in the three body plane tested (Figures <xref ref-type="fig" rid="F8">8A&#x02013;C</xref>). Interestingly, we observe a smaller amplitude non-linearity for low frequencies in the intermediated plane (Figure <xref ref-type="fig" rid="F8">8C</xref>; also compare across in Figures <xref ref-type="fig" rid="F7">7B1&#x02013;B3</xref>). Again, we can observe that the coupling between the two planes does no harm to the stability of the system.</p>
<fig id="F7" position="float">
<label>Figure 7</label>
<caption><p>Body COM sway responses of Lucy to PRTS support surface tilts of peak-peak amplitudes of 1&#x02013;8&#x000B0; <bold>(A)</bold> presented as averaged (<italic>n</italic> &#x0003D; 6) time series (<bold>B1</bold> sagittal plane, <bold>B2</bold> frontal plane, and <bold>B3</bold> intermediate plane). Note that the shown responses represent under-compensation of the stimuli, which is relatively more pronounced for small than for large stimuli (compare gain curves in Figure <xref ref-type="fig" rid="F8">8</xref>).</p></caption>
<graphic xlink:href="fnbot-11-00049-g0007.tif"/>
</fig>
<fig id="F8" position="float">
<label>Figure 8</label>
<caption><p>Lucy&#x00027;s sway responses in the sagittal, frontal and intermediate planes <bold>(A&#x02013;C)</bold> in terms of frequency response functions (FRFs) and coherence functions (bottom) from response averages of 6 repetitions of the PRTS stimulus. Shown are gain, phase and coherence curves over frequency for the four indicated peak&#x02013;peak stimulus amplitudes. Gain of zero would indicate ideal tilt compensation, a gain of unity that the evoked body COM excursion equals the tilt excursion. The phase gives the temporal response to stimulus relation. Coherence is a measure of the frequency dependent signal-to-noise ratio. Note non-linearity of gain curves (due to the threshold applied to <inline-formula><mml:math id="M76"><mml:mrow><mml:mover accent='true'><mml:mrow><mml:msubsup><mml:mover accent='true'><mml:mi>&#x003B1;</mml:mi><mml:mo>&#x002D9;</mml:mo></mml:mover><mml:mrow><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy='true'>&#x02323;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>) in <bold>(A,B)</bold>, less so in <bold>(C)</bold>.</p></caption>
<graphic xlink:href="fnbot-11-00049-g0008.tif"/>
</fig>
</sec>
<sec>
<title>Commanded fast full-body movements in the intermediate plane</title>
<p>This test challenged the robustness of the modular control in face of rapid self-produced movements in the sagittal-frontal 45&#x000B0; intermediate plane. The movement is associated with strong disturbances acting mainly through coupling forces arising between most of the robot&#x00027;s controlled DoFs. The test may be critical in a system with distributed modular controls that have time delays. Figure <xref ref-type="fig" rid="F9">9A</xref> shows the path of the COM sway in the two planes from the starting position at the right upper corner back to primary position (coordinates x &#x0003D; 0 cm, y &#x0003D; 0 cm). The complex path of the return reflects differences in the system&#x00027;s control and mechanical compliance in the two planes. The Figure <xref ref-type="fig" rid="F9">9B</xref> reports the temporal relaxation in the frontal plane (red) and the sagittal plane (green) following the step. The panel shows decaying oscillations of the relaxation. In the frontal plane, the rise time amounted to 0.584 s, the settling time to 5.115 s, and the overshoot to 33.2% of the step with a peak time of 1.620 s. The corresponding values for the sagittal plane were 0.621 s rise time, 5.141 s for the settling time, and an overshoot of 60.3%. This experiment shows differences in the dynamics in the two body planes, with larger oscillations in the sagittal plane due to higher compliance in this plane. Overall, it demonstrates stable performance when the control and mechanics in the robot&#x00027;s sagittal and frontal planes are dynamically interacting in rapid movements.</p>
<fig id="F9" position="float">
<label>Figure 9</label>
<caption><p><bold>(A,B)</bold> Return of body COM back to primary position in the course of repeated rapid voluntary movements of the COM in intermediate plane. Starting from an eccentric position (2.3 cm forward, 1.7 cm on the right), the robot was commanded to perform a rapid return of the COM to the upright position. The system response is shown in <bold>(A)</bold> and the COM trajectories are shown in <bold>(B)</bold> (for commands, see profiles of dashed and dotted step functions). <bold>(C,D)</bold> In order to visualize the coupling effects between the two body planes, additional experiments were performed where the robot moved only in the sagittal plane <bold>(C)</bold> and only in the frontal plane <bold>(D)</bold>. Note that all plots are displayed with the same scales to facilitate the comparisons. The control was always active in both planes to keep the system stable.</p></caption>
<graphic xlink:href="fnbot-11-00049-g0009.tif"/>
</fig>
<p>In the additional experiments shown in panels Figures <xref ref-type="fig" rid="F9">9C,D</xref>, the movement was restricted to the sagittal and frontal plane, respectively. They demonstrate that a rapid movement in one plane has only a small effect on the COM in the other plane. The cross talk from the sagittal to the frontal plane is again clearly larger than <italic>vice versa</italic>. Note that the responses in Figures <xref ref-type="fig" rid="F9">9C,D</xref> are similar to those shown in Figure <xref ref-type="fig" rid="F9">9B</xref> with respect to the actively moved component. This suggests that the controllers efficiently compensate the disturbance produced by the coupling between the planes. Note also that in all three experiments (Figures <xref ref-type="fig" rid="F9">9A&#x02013;C</xref>) the control was active in all the DoFs of the sagittal and frontal planes in order to keep the body upright.</p>
</sec>
<sec>
<title>Voluntary squatting movements</title>
<p>In this experiment, Lucy performed repetitive squatting knee-bending movements with reactive balancing movements occurring mainly in the ankle and hip joints [see section Experimental Procedures and Testbed (5)]. This experiment challenged Lucy&#x00027;s movement and balancing performance in a situation where all DoFs in the sagittal and frontal planes were interacting. The robot successfully performed this test (see film sections in Figure <xref ref-type="fig" rid="F10">10</xref> and film in <xref ref-type="supplementary-material" rid="SM1">Supplementary Materials</xref>). Lucy executes the voluntary and reactive movements in the sagittal plane, while the frontal plane controls prevent falling by keeping the body upright. No kinematic synergy is imposed explicitly and the resulting joint configurations are produced by the interactions between the modules.</p>
<fig id="F10" position="float">
<label>Figure 10</label>
<caption><p>Lucy performing squatting movements. The robot was commanded to perform knee bending at 0.17 Hz while balancing to keep the body upright. The observed coordination between knee, ankle, and hip movements (in terms of compensatory bending in hip and ankle joints) is a property emerging from the interaction between the control modules during the squatting movements (see text).</p></caption>
<graphic xlink:href="fnbot-11-00049-g0010.tif"/>
</fig>
</sec>
<sec>
<title>Voluntary full body movement at increasing frequency</title>
<p>In this experiment, Lucy performed voluntary body sway movements in the frontal plane and in the sagittal plane simultaneously. The reference trajectory was a sinusoid with increasing frequency (<italic>chirp</italic> signal) with the amplitude of 1&#x000B0; in the frontal plane and 2&#x000B0; in the sagittal plane. In contrast to the previous experiments, which were aimed to characterize the behavior produced by the DEC control within the margins of stability of the system (frequency range similar to that used in human posture control experiments), this trial was performed to the limit of failure (until the robot&#x00027;s feet lost contact with the support surface). The resulting movements are shown in Figure <xref ref-type="fig" rid="F11">11</xref>. In the low frequency range up to 0.4 Hz, tracking performance was almost accurate in both planes. With further increase in frequency, the responses in the sagittal plane develop a peak and then get smaller, while the responses in the frontal plane remain essentially similar as before and then increase in amplitude, before the response becomes unstable at 0.65 Hz, where the robot lost contact with the support surface. Thus, notably, this task pushed the robot to its stability limits. This owed mainly to limitations of the actuators, which were not designed for high frequencies and torques (which also would exceed human capabilities). Interferences between the sagittal and the frontal plane are relatively small, which supports our notion that the two planes of the robot&#x00027;s body can be controlled by two independent control systems.</p>
<fig id="F11" position="float">
<label>Figure 11</label>
<caption><p>Lucy performing a body sway in the ankle joints in the sagittal and the frontal plane simultaneously at increasing frequency (indicated at top of panels; time at bottom). The responses in the sagittal plane show a peak around 0.42 Hz. For higher frequencies, response amplitude tends to decrease. The response in the frontal plane increases markedly above 0.6 Hz. The recording was finished (right boundaries) when the robot&#x00027;s feet started to lose contact with the support surface due to large oscillations in the frontal plane. The difference in the dynamic responses in the two planes is mainly due to the body mechanics, i.e., the different size of the support base and the different mass distribution (which humans may account for by adjustments in the control, as future work may show).</p></caption>
<graphic xlink:href="fnbot-11-00049-g0011.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>Discussion</title>
<p>The main aim of our neurorobotics approach is to investigate whether the networks of DEC modules so far used separately for the sagittal plane and frontal plane (Lippi et al., <xref ref-type="bibr" rid="B28">2016</xref>; Ott et al., <xref ref-type="bibr" rid="B46">2016</xref>) would adequately cooperate with each other without using a supervising full body model or some other form of software linkage between the two networks. In particular, we asked whether the mechanical coupling between the planes given by the physics of the robot&#x00027;s body would provide postural stability when the robot performs movements in the 45&#x000B0; plane intermediate to the sagittal and frontal planes, or when it balances external perturbations in this plane. We envisaged that the underlying DEC with both voluntary movements and postural reactions follows the rules of vector decomposition from the intermediate into the sagittal and frontal planes. In the experiments 1&#x02013;5, the robot performed well within the stability margins based on the mechanical properties of the robot and the dynamics of the controls in the two planes, as ascertained empirically also in experiment 6 (Figure <xref ref-type="fig" rid="F11">11</xref>).</p>
<p>The described experiments and results considerably contribute to our aim of building a robotic system to further develop our human-derived DEC control with the ultimate goal of mirroring in robot experiments the human sensorimotor functions. This approach is complex and cannot be reached in one step, but requires several steps. The here-described steps comprise the cooperation between control modules in the frontal and sagittal planes of our modular control architecture, the testing of both reactive and proactive movement controls in view of control stability including mechanical aspects, and to provide measures such as frequency responses functions for robot-human comparisons. The current steps are important for further DEC developments that aim to implement in future steps more complex sensorimotor functions such as human-like walking.</p>
<p>The main computational challenge we expected in the present experiments was control stability in face of the feedback time delay of 20 ms, which is still much shorter than the known human time delay. Furthermore, the physical anisotropy of the robot&#x00027;s body in the sagittal and frontal planes represented a challenge. This arises among others from differences in body inertia and length of the base of support between the two planes, as reflected in the different responses to external disturbances and the different dynamic performance in voluntary squatting movements shown in the results section (Figures <xref ref-type="fig" rid="F5">5</xref>&#x02013;<xref ref-type="fig" rid="F8">8</xref>, <xref ref-type="fig" rid="F10">10</xref>). Yet, the robot successfully mastered with the same set of control parameters the disturbance and movement scenarios we applied. Thus, our finding demonstrate that the modular DEC control architecture is able to coordinate the movements across the robot&#x00027;s 14 DoF without signal exchange between the DEC module nets for the sagittal and the frontal planes or some higher order control mechanism. This is reminiscent of the so-called embodied approaches (e.g., Brooks, <xref ref-type="bibr" rid="B6">1991</xref>) according to which the body or other loops through the physical world can mediate interaction effects directly, i.e., without the need of explicit connections at control levels. The frequency responses functions obtained in the present experiments give us a comparison basis for future experiments in robots and humans.</p>
<p>The following discussion firstly considers related issues in humanoid robotics and the relevance for robotic neuro-rehabilitation (section Related Issues in Humanoid Robotics and Relevance for Robotic Neuro-Rehabilitation), then insights for the modeling of the human postural control (section Insights for the Modeling of the Human Postural Control), and finally future steps and developments expected for the extended DEC concept (section Future Steps and Developments Expected for the Extended DEC Concept).</p>
<sec>
<title>Related issues in humanoid robotics and relevance for robotic neuro-rehabilitation</title>
<p>Taking human bipedal control in standing and walking as a basis for comparison is an important research topic in humanoid robotics (Torricelli et al., <xref ref-type="bibr" rid="B60">2014</xref>). One reason is that human postural and movement skills are still considered to be superior to those of robots, being more robust and efficient, and also covering a wider range of external conditions (Nori et al., <xref ref-type="bibr" rid="B45">2014</xref>). Another reason is that humanoids acting in the human sphere may profit from human-like sensorimotor behavior when interacting with humans and their world.</p>
<p>Among the various tasks of sensorimotor control, maintaining balance is a <italic>primary task</italic> in the DEC concept as well as a basic rule that is followed in many fields of humanoid robotics. S<italic>econdary tasks</italic> can be performed in parallel with this primary task exploiting the kinematic redundancy of the robot, e.g., by projecting the secondary task into the null space of the Jacobian of the balancing task (Sentis and Khatib, <xref ref-type="bibr" rid="B57">2005</xref>). Care is generally taken to constrain secondary tasks such that they are not conflicting with the balancing task. In the present work, multiple tasks are achieved simultaneously, defined by different control variables for different modules. For example, in the squatting movement (section Voluntary Squatting Movements), the ankle module balances the body, while the knee module performs the vertical body motion and the other modules balance the COM of the supported body segments. In general, each module controls a joint (mechanical or virtual), along with all other segments supported by it. This works under the assumption that the controlled variable is affected directly by the controlled joint angle in that every rotation of the ankle joint produces the same rotation of the COM of the whole body around the ankle. For upright stance, the relation between COM sway and ankle joint is</p>
<disp-formula id="E17"><label>(17)</label><mml:math id="M77"><mml:mrow><mml:mfrac><mml:mo>&#x02202;</mml:mo><mml:mrow><mml:mo>&#x02202;</mml:mo><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>f</mml:mi></mml:mrow><mml:mn>1</mml:mn></mml:msubsup></mml:mrow></mml:mfrac><mml:msubsup><mml:mi>&#x003B1;</mml:mi><mml:mrow><mml:mi>b</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mn>1</mml:mn></mml:msubsup><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
<p>so that the servo controller can use the ankle control torque to control <inline-formula><mml:math id="M78"><mml:msubsup><mml:mrow><mml:mo>&#x003B1;</mml:mo></mml:mrow><mml:mrow><mml:mi>b</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msubsup></mml:math></inline-formula>. The dynamic effects of joint movements on the supporting links are here neglected. The integration of tasks in which a joint is controlled in order to produce an effect on the supporting links, e.g., using hip movements to control sheer forces under the feet, is currently still an open issue in the DEC concept.</p>
<p>An often-used balance control based on sensory signals is in humanoid robotics the method of the zero moment point (ZMP) criterion (Vukobratovi&#x00107; and Borovac, <xref ref-type="bibr" rid="B63">2004</xref>). It allows for balancing against moderate disturbances, which do not require hand contact or a step. This method has been successfully applied in robots with stiff actuation to adjust actual to desired ZMP, e.g., for walking (Hirai et al., <xref ref-type="bibr" rid="B18">1998</xref>; Sentis and Khatib, <xref ref-type="bibr" rid="B58">2006</xref>). However, many robots nowadays use compliant joints, as is the case with Lucy and the human system. The compliance has advantages for robot-world interactions such as collisions and for robot-human interactions. However, the control of compliant joints is more complex due to higher complexity of the dynamics. On the other hand, an important advantage of compliant actuation based on passive stiffness and damping is its immediate response to impact, starting well before time-consuming sensory feedback mechanisms take action (Haddadin et al., <xref ref-type="bibr" rid="B12">2007</xref>). This is comparable to the immediate passive stiffness &#x0201C;feedback loop&#x0201D; in humans from muscles, tendons, and connective tissues, where previous modeling and robot experiments suggested an improved control stability in face of considerable sensory, neural and muscular time delays (Antritter et al., <xref ref-type="bibr" rid="B2">2014</xref>; Ott et al., <xref ref-type="bibr" rid="B46">2016</xref>). While passive stiffness was implemented in Posturob I and II of this laboratory by using pneumatic muscles and in some experiments springs as tendons (e.g., Mergner et al., <xref ref-type="bibr" rid="B38">2009</xref>), it has not yet been implemented in Lucy. This is one of the aims for future experiments with Lucy (see section Future Steps and Developments Expected for the Extended DEC Concept).</p>
<p>In the DEC-controlled robots including Lucy, compliant behavior and low energy consumption are positive side effects from the low control loop gain, which in humans appears to be mainly related to control stability in face of the biological time delays. Effects of the compliance showed up in the presented experiments in the form of residual body sway following external disturbances as well as overshoot and under-damped dynamics with fast voluntary movements. A particularly important beneficial side effect of using low actuation torques was that the robot&#x00027;s feet did not loose contact with the support surface in the experiments 1&#x02013;5. The effects that passive stiffness has for reducing impact magnitude and improving energy efficiency has been an issue in actuator design (e.g., Ham et al., <xref ref-type="bibr" rid="B14">2009</xref>) and recently has received considerable interest in humanoid robotics. The method has been combined for example with active compliance modulation in the robot COMAN (Li et al., <xref ref-type="bibr" rid="B26">2012</xref>). In this robot, joint stiffness is modulated in response to an external disturbing force by shifting desired body COM such that the resulting center of pressure (COP) shift compensates for the disturbance. This mechanism is employed in combination with other assisting mechanisms that controlled upper body orientation and energy dissipation. Furthermore, flexible robots such a humanoids with compliant &#x0201C;ankle&#x0201D; and &#x0201C;hip&#x0201D; joint equivalents have been controlled with the so-called Reaction Null Space (RNS) formalisms that aim for a reactionless motion control via a feed forward mechanism and an error compensation via feedback (Nenchev, <xref ref-type="bibr" rid="B43">2013</xref>).</p>
<p>Currently, our experiments do not include large perturbations; this topic remains to be addressed in future work (section Future Steps and Developments Expected for the Extended DEC Concept). For example, explicit inspirations from human balancing research for humanoids often address the &#x0201C;ankle strategy&#x0201D; and &#x0201C;hip strategy&#x0201D; (Nishio et al., <xref ref-type="bibr" rid="B44">2006</xref>). In the &#x0201C;ankle strategy,&#x0201D; ankle torque suffices to produce shifts of the body COP to compensate moderate perturbations, whereas the &#x0201C;hip strategy&#x0201D; becomes involved or dominates when the COP is expected to exceed the limits of the support base for balancing, be these limits determined by the feet or a restricted or compliant support surface. Another limiting factor is the maximal ankle torque. The hip movement then generates horizontal ground forces in order to keep the COM over the base of support (Nashner and McCollum, <xref ref-type="bibr" rid="B42">1985</xref>). These bio-inspired dynamic balance mechanisms have been technically realized in a compliant humanoid robot by Hyon et al. (<xref ref-type="bibr" rid="B19">2007</xref>) and other groups. Measures related to COM and COP shifts were found to describe postural stability decision limits with increasing perturbation magnitude in three steps: from (1) &#x0201C;COP balancing&#x0201D; in the ankle joints via (2) &#x0201C;centroidal moment point, CMP, balancing&#x0201D; in the hip joints involving a moment about the COM, up to (3) a rescue step, often resulting in a double support (Stephens, <xref ref-type="bibr" rid="B59">2007</xref>). Simulations showed that human &#x0201C;hip&#x0201D; and &#x0201C;ankle&#x0201D; balancing strategies tend to emerge in the presence of perturbations of different magnitudes from the same optimization criterion (Atkeson and Stephens, <xref ref-type="bibr" rid="B4">2007</xref>). These robotic studies in turn provide inspirations from robotics back to human posture control research.</p>
<p>As stated in section Introduction, an important aspect of the DEC control is its modularity. Modular concepts have been widely applied before to both the analysis of human movement and the control of robots. In general, the concept of modular control implies the subdivision of the control problem into several sub-problems that can be managed by separate control modules. Such modules may interact in different ways, e.g., directly by exchange of signals or through the effects they produce on the system. Control modularity has been defined at different levels of abstraction. For example, in the &#x0201C;behavior-based&#x0201D; control systems several sub-behaviors may operate in concert with each other, organized according to a hierarchical principle (Brooks, <xref ref-type="bibr" rid="B6">1991</xref>). In the context of posture control for biped humanoids, such sub-behaviors can be used to control low-level tasks like joint torque profiles, or to process higher-level issues such as adapting control parameters to contextual situations, or performing specific reactive adjustments. For example, an implementation on a simulated humanoid is shown in Luksch (<xref ref-type="bibr" rid="B30">2010</xref>).</p>
<p>Higher-level functionalities can emerge from the complexity of behavior based systems. Recently, the integration of internal representations has been considered in behavior based system research; this allows for a direct solution of issues such as movement planning. It finds its basis in that several tasks in biological systems are considered to rely on internal representations for inverse modeling, forward models and sensor fusion (Schilling, <xref ref-type="bibr" rid="B54">2011</xref>). Internal representations are different from a behavior in that they are not designed specifically as tasks, but instead have a general validity, as e.g., realized in the principle of <italic>mean of multiple computations</italic> (MMC). This method can represent body kinematics independently from the specific task and be generalized to arbitrary joint configurations (e.g., of a planar arm in Schilling, <xref ref-type="bibr" rid="B54">2011</xref> or of a hexapod in Schilling and Cruse, <xref ref-type="bibr" rid="B55">2008</xref>). The DEC control, on the other hand, is specifically addressing posture control and disturbance rejection. Yet, although its bottom-up control design by the integration of locally acting modules shows interesting emerging properties such as segment coordination, it is not expected in general to cover higher level aspects of motor control. An integration of DEC with the above methods can be imagined, e.g., in a set up where high level processing receives sensory input preprocessed through the DEC sensor fusions and where the high level controller outputs into the DEC posture control, e.g., estimates of disturbances that are self-produced by voluntary movements.</p>
<p>Another modularity concept, the Modular Modality Frame model (Ehrenfeld and Butz, <xref ref-type="bibr" rid="B7">2013</xref>), envisages multiple cross-coordinations between cortical sensorimotor representations for processing states of the body and the body parts and the integration of corresponding multisensory information with Bayesian optimality. Again, this concept addresses probabilistic mechanisms of the body schema for movement planning and commanding, whereas the DEC concept deals with the postural control that, on a lower level, helps movement execution by disturbance compensation. Yet, parallels to the DEC concept possibly exist concerning multisensory processing. Specifically, the DEC concept combines in a first step signals from a variety of transducers to obtain estimates of relevant physical variables. In a second step, from these estimated physical variables estimates then the four disturbances are derived. This two-step transformation of transducer signals into estimates of disturbances occurs without probability weighting. The processing in the DEC involves at least partially known transfer characteristics and noise properties and aims to fulfill the commanded task, using a direct reconstruction of the physical relationship between disturbing forces, joint positions and transducer stimulation.</p>
<p>The uncertainty of the results caused by the sensor noise (mainly vestibular noise; van der Kooij and Peterka, <xref ref-type="bibr" rid="B62">2011</xref>) can partially be mitigated by thresholds (Mergner et al., <xref ref-type="bibr" rid="B38">2009</xref>) and disturbance prediction (Mergner, <xref ref-type="bibr" rid="B33">2010</xref>). Estimated disturbances smaller than a given threshold are excluded from the compensation and the control process. This prevents noise from sensory signals that are not participating in the control from influencing the system. To which extent this approach fulfills optimality criteria remains to be investigated.</p>
<p>Exchange of knowledge and inspirations between the robotic and the human domains is especially relevant when it comes to the use of robotic sensorimotor devices such as exoskeletons for neuro-rehabilitation and assistive devices in the development of new invasive or non-invasive stimulation methods to intervene with human sensorimotor functions and more. Knowing the neural control mechanisms behind human sensorimotor behavior will likely prove to be a key for ensuring an intuitive and safe use of such devices by patients.</p>
<p>A related consideration concerns the observation in the present experiments that the robot was able to balance across a broad range of scenarios and conditions using the same sets of control parameters, while automatically adjusting its motor performance to changes in stimulus modality and amplitude, this in combination with the afore discussed low loop gain, and soft mechanical compliance and low energy consumption. Versatility of the control is instrumental for the sensorimotor behavior in a rich and changing environment, which represents a major and still open challenge in the field of humanoid robotics and robotics-inspired assistive devices. The concept of a modular DEC architecture may provide a realistic solution in face of these challenges.</p>
</sec>
<sec>
<title>Insights for the modeling of the human postural control</title>
<p>The coordination of disturbance compensation observed in our experiments when using support tilts in the intermediate plane resulted from mechanical forces that were produced by the controls for the sagittal and frontal planes. A human skeletal body equipped with sensors, actuators and DEC control modules can be expected to function similarly to Lucy. We take this as evidence that the DEC control can be considered as a valid control concept for human reactive postural control. This even applies if the realization in humans would differ in some respect from that of our robot. Differences may owe to the much larger number of sensory transducer and actuators in humans as compared to the robot, which may reflect the fact that humans are experiencing and producing a richer spectrum of sensorimotor behaviors than currently faced by Lucy.</p>
<p>One may object that here exists a conflict between the DEC concept and the notion of those neuroscientists, who assume that movement planning and commanding at the level of the cortical homunculus (compare section Introduction) entails a full-body control. Note, however, that the cortical motor commands, with few exceptions, do not reach the spinal motor centers directly. Rather, the commands pass through and interact with subcortical structures, mainly the basal ganglia and cerebellum, where they undergo modifications. Thereby, they may better fit to the DEC concept. In this view, one important role of these subcortical structures is to complement the voluntary movement commands by predicted-sensory disturbance estimates. These centrally generated estimates are thought to represent fast and low-noisy predictions of the sensory-derived disturbance estimates evoked by the movements, which they substitute for compensation. This allows for a distinction between self-produced and external disturbances and thus allows compensating both even when they occur in superposition (Maurer et al., <xref ref-type="bibr" rid="B31">2006</xref>; Mergner, <xref ref-type="bibr" rid="B33">2010</xref>).</p>
<p>Although the robot experiments provide no direct insights into human postural control, they may provide valuable support for one or the other hypothetical model when this is tested in robot experiments for &#x0201C;real world&#x0201D; robustness (see section Introduction). This neurorobotics approach will allow neuroscientists and roboticists to compare different models both with respect to specific and general performance criteria (e.g., versatility, failsafe robustness, etc.). For example, a concept of how humans may deal with inter-link coupling forces using Eigenmovements has recently been tested in a robot of this laboratory, where it successfully passed the &#x0201C;real world&#x0201D; test (Alexandrov et al., <xref ref-type="bibr" rid="B1">2017</xref>). The DEC concept provides an alternative solution for the problem of inter-segmental force coupling in terms of a bundle of counter measures, which include the DEC mechanisms for <italic>T</italic><sub><italic>ext</italic></sub> and <italic>T</italic><sub><italic>acc</italic></sub> in Figure <xref ref-type="fig" rid="F1">1</xref>, modulation of passive stiffness and damping, and the use of waveforms for commanded movement trajectories, which are optimized for low acceleration and jerk (compare waveform in Figure <xref ref-type="fig" rid="F6">6</xref>).</p>
<p>In the previous robot experiments of Hettich et al. (<xref ref-type="bibr" rid="B17">2014</xref>), the DEC concept was restricted to the DIP scenario and used for the ankle and hip control estimated parameters, which were directly derived from the modeling of human reactive movements. The control parameters related to mass and COM heights of the robot&#x00027;s segments previously used were adjusted in the present experiments to Lucy accordingly. This does not apply, however, to the lumped time delay. Instead, a delay of 20 ms was used for all joints, which is clearly less than that used previously or and cannot be considered &#x0201C;human-like.&#x0201D; Identified lumped delays for the ankle joint control in our previous experiments were &#x02265;100 ms. We have not yet tested such lumped delays in Lucy, but plan to test them together with a number of yet pending robot implementations of related other issues such as the aforementioned envisaged solution for the inter-segmental coupling forces, the up-channeling of control parameters, and the passive stiffness and damping. A further aim is to replace the concept of a lumped time delay by the biologically more appropriate concept of &#x0201C;short latency reflexes&#x0201D; (local proprioceptive mechanisms with latencies of 20&#x02013;40 ms) and &#x0201C;long latency reflexes&#x0201D; (responses with latencies of 60&#x02013;300 ms, which involve volition and intention; compare Mergner, <xref ref-type="bibr" rid="B33">2010</xref>). With the integration of these mechanisms we expect stabilizing effects also for coupling forces and a better dynamic performance in faces of challenges that may occur during walking. This also may apply when horizontal body segment movements are included into the control, such as rapid head movements during gaze shifts (Falotico et al., <xref ref-type="bibr" rid="B8">2017</xref>).</p>
<p>In future experiments, we also will address in more detail the emergence of conflict-free interactions between control constituents, as we observed it in our previous and present robot experiments. It allowed here the superposition of voluntary movements and reactive compensations of external disturbances, but also may allow superposition of two or more disturbances at a time. Compliance to the superposition law may represent an important criterion when judging a given robot performance as human-like.</p>
</sec>
<sec>
<title>Future steps and developments expected for the extended DEC concept</title>
<p>To develop the robotic model of the human sensorimotor system further, we consider the following steps important:
<list list-type="alpha-upper">
<list-item><p>Installing the ability for an automatic reconfiguration of the control model for certain tasks such as walking. There, it will account for changes in the robot&#x00027;s configuration with the alternation between the double leg support and the single leg support while the other leg swings. Smooth human-like walking in future DEC implementations may furthermore benefit from DEC extensions in terms of whole body coordination, which may includes the upper body, and it will involve integration of an open loop gait pattern (compare Lapeyre et al., <xref ref-type="bibr" rid="B24">2013a</xref>,<xref ref-type="bibr" rid="B25">b</xref>).</p></list-item>
<list-item><p>Predictions of disturbance estimates through both down- and up-channeling of signals that are produced during <italic>voluntary and reactive</italic> movements&#x02014;with an expected improvement of control stability from lower noise and shorter time delay as compared to sensor-derived signals.</p></list-item>
<list-item><p>Further improvement of control stability by using passive stiffness modulation during contact force and support translation disturbances and during rapid voluntary and passive movements.</p></list-item>
<list-item><p>Addition of visual self-motion information to improve the balance performance of Lucy. Improvements of accuracy and noise of vestibular and proprioceptive sensor information through interaction with the visual signal can be expected on the basis of human experiments (e.g., Assl&#x000E4;nder et al., <xref ref-type="bibr" rid="B3">2015</xref>).</p></list-item>
<list-item><p>Superposition of more than one external perturbation (e.g., support surface tilts and contact forces) during quite stance and voluntary movements.</p></list-item>
<list-item><p>Implementation of human-derived lumped time delays.</p></list-item>
<list-item><p>Attempts to mimic sensorimotor impairments of neurological patients using DEC-controlled robots (compare Mergner et al., <xref ref-type="bibr" rid="B38">2009</xref>, for balance model without vestibular function).</p></list-item>
</list></p>
</sec>
</sec>
<sec id="s5">
<title>Author contributions</title>
<p>VL and TM performed the experiment and collected and analyzed the data, performed the computer simulations. Both authors contributed to the interpretation of the data, drafted and critically revised the manuscript and finally approved the manuscript for submission.</p>
<sec>
<title>Conflict of interest statement</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p></sec>
</sec>
</body>
<back>
<sec sec-type="supplementary-material" id="s6">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="http://journal.frontiersin.org/article/10.3389/fnbot.2017.00049/full#supplementary-material">http://journal.frontiersin.org/article/10.3389/fnbot.2017.00049/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Video1.MP4" id="SM1" mimetype="video/mp4" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Alexandrov</surname> <given-names>A. V.</given-names></name> <name><surname>Lippi</surname> <given-names>V.</given-names></name> <name><surname>Mergner</surname> <given-names>T.</given-names></name> <name><surname>Frolov</surname> <given-names>A. A.</given-names></name> <name><surname>Hettich</surname> <given-names>G.</given-names></name> <name><surname>Husek</surname> <given-names>D.</given-names></name></person-group> (<year>2017</year>). <article-title>Human-inspired eigenmovement concept provides coupling-free sensorimotor control in humanoid robot</article-title>. <source>Front. Neurorobot.</source> <volume>11</volume>:<fpage>22</fpage>. <pub-id pub-id-type="doi">10.3389/fnbot.2017.00022</pub-id><pub-id pub-id-type="pmid">28487646</pub-id></citation></ref>
<ref id="B2">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Antritter</surname> <given-names>F.</given-names></name> <name><surname>Scholz</surname> <given-names>F.</given-names></name> <name><surname>Hettich</surname> <given-names>G.</given-names></name> <name><surname>Mergner</surname> <given-names>T.</given-names></name></person-group> (<year>2014</year>). <article-title>Stability analysis of human stance control from the system theoretic point of view</article-title>, in <source>2014 European Control Conference, ECC 2014</source> (<publisher-loc>Strasbourg</publisher-loc>), <fpage>1849</fpage>&#x02013;<lpage>1855</lpage>.</citation></ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Assl&#x000E4;nder</surname> <given-names>L.</given-names></name> <name><surname>Hettich</surname> <given-names>G.</given-names></name> <name><surname>Mergner</surname> <given-names>T.</given-names></name></person-group> (<year>2015</year>). <article-title>Visual contribution to human standing balance during support surface tilts</article-title>. <source>Hum. Mov. Sci.</source> <volume>41</volume>, <fpage>147</fpage>&#x02013;<lpage>164</lpage>. <pub-id pub-id-type="doi">10.1016/j.humov.2015.02.010</pub-id><pub-id pub-id-type="pmid">25816794</pub-id></citation></ref>
<ref id="B4">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Atkeson</surname> <given-names>C. G.</given-names></name> <name><surname>Stephens</surname> <given-names>B.</given-names></name></person-group> (<year>2007</year>). <article-title>Multiple balance strategies from one optimization criterion</article-title>, in <source>7th IEEE-RAS International Conference on Humanoid Robots</source> (<publisher-loc>Pittsburgh, PA</publisher-loc>), <fpage>57</fpage>&#x02013;<lpage>64</lpage>.</citation></ref>
<ref id="B5">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Bleisteiner</surname> <given-names>G. V.</given-names></name> <name><surname>Mangoldt</surname> <given-names>W.</given-names></name> <name><surname>Henninger</surname> <given-names>H.</given-names></name> <name><surname>Oetker</surname> <given-names>R.</given-names></name></person-group> (<year>1961</year>). <source>Handbuch der Regelungstechnik</source>. Berlin; <publisher-loc>Heidelberg</publisher-loc>: <publisher-name>Springer</publisher-name>.</citation></ref>
<ref id="B6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Brooks</surname> <given-names>R.</given-names></name></person-group> (<year>1991</year>). <article-title>New approaches to robotics</article-title>. <source>Science</source> <volume>253</volume>, <fpage>1227</fpage>&#x02013;<lpage>1232</lpage>. <pub-id pub-id-type="doi">10.1126/science.253.5025.1227</pub-id><pub-id pub-id-type="pmid">17831441</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ehrenfeld</surname> <given-names>S.</given-names></name> <name><surname>Butz</surname> <given-names>M. V.</given-names></name></person-group> (<year>2013</year>). <article-title>The modular modality frame model: continuous body state estimation and plausibility-weighted information fusion</article-title>. <source>Biol. Cybern.</source> <volume>107</volume>, <fpage>61</fpage>&#x02013;<lpage>82</lpage>. <pub-id pub-id-type="doi">10.1007/s00422-012-0526-2</pub-id><pub-id pub-id-type="pmid">23090574</pub-id></citation></ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Falotico</surname> <given-names>E.</given-names></name> <name><surname>Cauli</surname> <given-names>N.</given-names></name> <name><surname>Kryczka</surname> <given-names>P.</given-names></name> <name><surname>Hashimoto</surname> <given-names>K.</given-names></name> <name><surname>Berthoz</surname> <given-names>A.</given-names></name> <name><surname>Takanishi</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Head stabilization in a humanoid robot: models and implementations</article-title>. <source>Auton. Robots</source> <volume>41</volume>, <fpage>349</fpage>&#x02013;<lpage>365</lpage>. <pub-id pub-id-type="doi">10.1007/s10514-016-9583-z</pub-id></citation></ref>
<ref id="B9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fitzpatrick</surname> <given-names>R.</given-names></name> <name><surname>Burke</surname> <given-names>D.</given-names></name> <name><surname>Gandevia</surname> <given-names>S. C.</given-names></name></person-group> (<year>1996</year>). <article-title>Loop gain of reflexes controlling human standing measured with the use of postural and vestibular disturbances</article-title>. <source>J. Neurophysiol</source>. <volume>76</volume>, <fpage>3994</fpage>&#x02013;<lpage>4008</lpage>. <pub-id pub-id-type="pmid">8985895</pub-id></citation></ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Flash</surname> <given-names>T.</given-names></name> <name><surname>Hochner</surname> <given-names>B.</given-names></name></person-group> (<year>2005</year>). <article-title>Motor primitives in vertebrates and invertebrates</article-title>. <source>Curr. Opin. Neurobiol.</source> <volume>15</volume>, <fpage>660</fpage>&#x02013;<lpage>666</lpage>. <pub-id pub-id-type="doi">10.1016/j.conb.2005.10.011</pub-id><pub-id pub-id-type="pmid">16275056</pub-id></citation></ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Goodworth</surname> <given-names>A. D.</given-names></name> <name><surname>Peterka</surname> <given-names>R. J.</given-names></name></person-group> (<year>2010</year>). <article-title>Influence of stance width on frontal plane postural dynamics and coordination in human balance control</article-title>. <source>J. Neurophysiol.</source> <volume>104</volume>, <fpage>1103</fpage>&#x02013;<lpage>1118</lpage>. <pub-id pub-id-type="doi">10.1152/jn.00916.2009</pub-id><pub-id pub-id-type="pmid">20427616</pub-id></citation></ref>
<ref id="B12">
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Haddadin</surname> <given-names>S.</given-names></name> <name><surname>Schaffer</surname> <given-names>A.</given-names></name> <name><surname>Hirzinger</surname> <given-names>G.</given-names></name></person-group> (<year>2007</year>). <article-title>Safety evaluation of physical human-robot interaction via crash-testing</article-title>, in <source>Robotics: Science and Systems Conference (RSS 2007)</source>, <fpage>395</fpage>&#x02013;<lpage>407</lpage>. Available online at: <ext-link ext-link-type="uri" xlink:href="http://www.phriends.org/RSS_07b.pdf">http://www.phriends.org/RSS_07b.pdf</ext-link></citation></ref>
<ref id="B13">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Hajos</surname> <given-names>A.</given-names></name> <name><surname>Kirchner</surname> <given-names>W.</given-names></name></person-group> (<year>1984</year>). <article-title>The control of body posture in man during electrical stimulation of the vestibular system</article-title>, in <source>Sensory Experience, Adaptation and Perception</source>, eds <person-group person-group-type="editor"><name><surname>Spillmann</surname> <given-names>L.</given-names></name> <name><surname>Wooten</surname> <given-names>B. R.</given-names></name></person-group> (<publisher-loc>Hillsdale, MI; London</publisher-loc>: <publisher-name>Lawrence Erlbaum</publisher-name>), <fpage>255</fpage>&#x02013;<lpage>280</lpage>.</citation></ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ham</surname> <given-names>R.</given-names></name> <name><surname>Sugar</surname> <given-names>T.</given-names></name> <name><surname>Vanderborght</surname> <given-names>B.</given-names></name> <name><surname>Hollander</surname> <given-names>K.</given-names></name> <name><surname>Lefeber</surname> <given-names>D.</given-names></name></person-group> (<year>2009</year>). <article-title>Compliant actuator designs</article-title>. <source>IEEE Robot. Automat. Mag.</source> <volume>16</volume>, <fpage>81</fpage>&#x02013;<lpage>94</lpage>. <pub-id pub-id-type="doi">10.1109/MRA.2009.933629</pub-id></citation></ref>
<ref id="B15">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Hauser</surname> <given-names>H.</given-names></name> <name><surname>Neumann</surname> <given-names>G.</given-names></name> <name><surname>Ijspeert</surname> <given-names>A. J.</given-names></name> <name><surname>Maass</surname> <given-names>W.</given-names></name></person-group> (<year>2007</year>). <article-title>Biologically inspired kinematic synergies provide a new paradigm for balance control of humanoid robots</article-title>, in <source>Humanoid Robots, 2007 7th IEEE-RAS International Conference</source> (<publisher-loc>Pittsburgh, PA</publisher-loc>), <fpage>73</fpage>&#x02013;<lpage>80</lpage>.</citation></ref>
<ref id="B16">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Hettich</surname> <given-names>G.</given-names></name> <name><surname>Lippi</surname> <given-names>V.</given-names></name> <name><surname>Mergner</surname> <given-names>T.</given-names></name></person-group> (<year>2013</year>). <article-title>Human-like sensor fusion mechanisms in a postural control robot</article-title>, in <source>Proceedings of the International Congress on Neurotechnology, Electronics and Informatics</source> (<publisher-loc>Freiburg</publisher-loc>), <fpage>152</fpage>&#x02013;<lpage>160</lpage>.</citation></ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hettich</surname> <given-names>G.</given-names></name> <name><surname>Assl&#x000E4;nder</surname> <given-names>L.</given-names></name> <name><surname>Gollhofer</surname> <given-names>A.</given-names></name> <name><surname>Mergner</surname> <given-names>T.</given-names></name></person-group> (<year>2014</year>). <article-title>Human hip-ankle coordination emerging from multisensory feedback control</article-title>. <source>Hum. Mov. Sci.</source> <volume>37</volume>, <fpage>123</fpage>&#x02013;<lpage>146</lpage>. <pub-id pub-id-type="doi">10.1016/j.humov.2014.07.004</pub-id><pub-id pub-id-type="pmid">25150802</pub-id></citation></ref>
<ref id="B18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hirai</surname> <given-names>K.</given-names></name> <name><surname>Hirose</surname> <given-names>M.</given-names></name> <name><surname>Haikawa</surname> <given-names>Y.</given-names></name> <name><surname>Takenaka</surname> <given-names>T.</given-names></name></person-group> (<year>1998</year>). <article-title>The development of Honda humanoid robot</article-title>. <source>IEEE Int. Conf. Robot. Automat.</source> <volume>2</volume>, <fpage>1321</fpage>&#x02013;<lpage>1326</lpage>. <pub-id pub-id-type="doi">10.1109/ROBOT.1998.677288</pub-id></citation></ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hyon</surname> <given-names>S. H.</given-names></name> <name><surname>Hale</surname> <given-names>J. G.</given-names></name> <name><surname>Cheng</surname> <given-names>G.</given-names></name></person-group> (<year>2007</year>). <article-title>Full-body compliant human-humanoid interaction: balancing in the presence of unknown external forces</article-title>. <source>IEEE Trans. Robot.</source> <volume>23</volume>, <fpage>884</fpage>&#x02013;<lpage>898</lpage>. <pub-id pub-id-type="doi">10.1109/TRO.2007.904896</pub-id></citation></ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ivanenko</surname> <given-names>Y. P.</given-names></name> <name><surname>Grasso</surname> <given-names>R.</given-names></name> <name><surname>Zago</surname> <given-names>M.</given-names></name> <name><surname>Molinari</surname> <given-names>M.</given-names></name> <name><surname>Scivoletto</surname> <given-names>G.</given-names></name> <name><surname>Castellano</surname> <given-names>V.</given-names></name> <etal/></person-group>. (<year>2003</year>). <article-title>Temporal components of the motor patterns expressed by the human spinal cord reflect foot kinematics</article-title>. <source>J. Neurophysiol.</source> <volume>90</volume>, <fpage>3555</fpage>&#x02013;<lpage>3565</lpage>. <pub-id pub-id-type="doi">10.1152/jn.00223.2003</pub-id><pub-id pub-id-type="pmid">12853436</pub-id></citation></ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Johansson</surname> <given-names>R.</given-names></name> <name><surname>Magnusson</surname> <given-names>M.</given-names></name> <name><surname>Akesson</surname> <given-names>M.</given-names></name></person-group> (<year>1988</year>). <article-title>Identification of human postural dynamics</article-title>. <source>IEEE Trans. Biomed. Eng.</source> <volume>35</volume>, <fpage>858</fpage>&#x02013;<lpage>869</lpage>. <pub-id pub-id-type="doi">10.1109/10.7293</pub-id><pub-id pub-id-type="pmid">3192235</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kuo</surname> <given-names>A. D.</given-names></name></person-group> (<year>1995</year>). <article-title>An optimal control model for analyzing human postural balance</article-title>. <source>IEEE Trans. Biomed. Eng.</source> <volume>42</volume>, <fpage>87</fpage>&#x02013;<lpage>101</lpage>. <pub-id pub-id-type="doi">10.1109/10.362914</pub-id><pub-id pub-id-type="pmid">7851935</pub-id></citation></ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kuo</surname> <given-names>A. D.</given-names></name></person-group> (<year>2005</year>). <article-title>An optimal state estimation model of sensory integration in human postural balance</article-title>. <source>J. Neur. Eng.</source> <volume>2</volume>:<fpage>S235</fpage>. <pub-id pub-id-type="doi">10.1088/1741-2560/2/3/S07</pub-id><pub-id pub-id-type="pmid">16135887</pub-id></citation></ref>
<ref id="B24">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Lapeyre</surname> <given-names>M.</given-names></name> <name><surname>Rouanet</surname> <given-names>P.</given-names></name> <name><surname>Oudeyer</surname> <given-names>P. Y.</given-names></name></person-group> (<year>2013a</year>). <article-title>The poppy humanoid robot: Leg design for biped locomotion</article-title>, in <source>2013 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS)</source> (<publisher-loc>Tokyo</publisher-loc>), <fpage>349</fpage>&#x02013;<lpage>356</lpage>.</citation></ref>
<ref id="B25">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Lapeyre</surname> <given-names>M.</given-names></name> <name><surname>Rouanet</surname> <given-names>P.</given-names></name> <name><surname>Oudeyer</surname> <given-names>P. Y.</given-names></name></person-group> (<year>2013b</year>). <article-title>Poppy humanoid platform: Experimental evaluation of the role of a bio-inspired thigh shape</article-title>, in <source>2013 13th IEEE-RAS International Conference on Humanoid Robots (Humanoids)</source> (<publisher-loc>Tokyo</publisher-loc>), <fpage>376</fpage>&#x02013;<lpage>383</lpage>.</citation></ref>
<ref id="B26">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>Z.</given-names></name> <name><surname>Vanderborght</surname> <given-names>B.</given-names></name> <name><surname>Tsagarakis</surname> <given-names>N. G.</given-names></name> <name><surname>Colasanto</surname> <given-names>L.</given-names></name> <name><surname>Caldwell</surname> <given-names>D. G.</given-names></name></person-group> (<year>2012</year>). <article-title>Stabilization for the compliant humanoid robot COMAN exploiting intrinsic and controlled compliance</article-title>, in <source>IEEE International Conference on Robotics and Automation</source> (<publisher-loc>St. Paul, MN</publisher-loc>), <fpage>2000</fpage>&#x02013;<lpage>2006</lpage>.</citation></ref>
<ref id="B27">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Lippi</surname> <given-names>V.</given-names></name> <name><surname>Mergner</surname> <given-names>T.</given-names></name> <name><surname>Hettich</surname> <given-names>G.</given-names></name></person-group> (<year>2013</year>). <article-title>A Bio-inspired modular system for humanoid posture control</article-title>, in <source>Proceedings of IROS 2013 Workshop on Neuroscience and Robotics &#x0201C;Towards a robot-enabled, neuroscience-guided healthy society,&#x0201D;</source> eds <person-group person-group-type="editor"><name><surname>Ugur</surname> <given-names>E.</given-names></name> <name><surname>Oztop</surname> <given-names>E.</given-names></name> <name><surname>Morimoto</surname> <given-names>J.</given-names></name> <name><surname>Ishii</surname> <given-names>S.</given-names></name></person-group> (<publisher-loc>Tokyo</publisher-loc>).</citation></ref>
<ref id="B28">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Lippi</surname> <given-names>V.</given-names></name> <name><surname>Mergner</surname> <given-names>T.</given-names></name> <name><surname>Szumowski</surname> <given-names>M.</given-names></name> <name><surname>Zurawska</surname> <given-names>M. S.</given-names></name> <name><surname>Zieli&#x00144;ska</surname> <given-names>T.</given-names></name></person-group> (<year>2016</year>). <article-title>Human-inspired humanoid balancing and posture control in frontal plane</article-title>, in <source>ROMANSY 21-Robot Design, Dynamics and Control</source>, eds <person-group person-group-type="editor"><name><surname>Parenti-Castelli</surname> <given-names>V.</given-names></name> <name><surname>Schiehlen</surname> <given-names>W.</given-names></name></person-group> (<publisher-loc>Udine</publisher-loc>: <publisher-name>Springer International Publishing</publisher-name>), <fpage>285</fpage>&#x02013;<lpage>292</lpage>.</citation></ref>
<ref id="B29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Luecke</surname> <given-names>R. H.</given-names></name> <name><surname>McGuire</surname> <given-names>M. L.</given-names></name></person-group> (<year>1968</year>). <article-title>Analysis of optimal composite feedback-feedforward control</article-title>. <source>AIChE J.</source> <volume>14</volume>, <fpage>181</fpage>&#x02013;<lpage>189</lpage>. <pub-id pub-id-type="doi">10.1002/aic.690140131</pub-id></citation></ref>
<ref id="B30">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Luksch</surname> <given-names>T.</given-names></name></person-group> (<year>2010</year>). <source>Human-Like Control of Dynamically Walking Bipedal Robots</source>. <publisher-loc>Munich</publisher-loc>: <publisher-name>Verlag Dr. Hut</publisher-name>.</citation></ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Maurer</surname> <given-names>C.</given-names></name> <name><surname>Mergner</surname> <given-names>T.</given-names></name> <name><surname>Peterka</surname> <given-names>R. J.</given-names></name></person-group> (<year>2006</year>). <article-title>Multisensory control of human upright stance</article-title>. <source>Exp. Brain Res.</source> <volume>171</volume>, <fpage>231</fpage>&#x02013;<lpage>250</lpage>. <pub-id pub-id-type="doi">10.1007/s00221-005-0256-y</pub-id><pub-id pub-id-type="pmid">16307252</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>McIntyre</surname> <given-names>J.</given-names></name> <name><surname>Bizzi</surname> <given-names>E.</given-names></name></person-group> (<year>1993</year>). <article-title>Servo hypotheses for the biological control of movement</article-title>. <source>J. Mot. Behav.</source> <volume>25</volume>, <fpage>193</fpage>&#x02013;<lpage>202</lpage>. <pub-id pub-id-type="pmid">12581989</pub-id></citation></ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mergner</surname> <given-names>T.</given-names></name></person-group> (<year>2010</year>). <article-title>A neurological view on reactive human stance control</article-title>. <source>Annu. Rev. Control</source> <volume>34</volume>, <fpage>177</fpage>&#x02013;<lpage>198</lpage>. <pub-id pub-id-type="doi">10.1016/j.arcontrol.2010.08.001</pub-id></citation></ref>
<ref id="B34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mergner</surname> <given-names>T.</given-names></name> <name><surname>Huber</surname> <given-names>W.</given-names></name> <name><surname>Becker</surname> <given-names>W.</given-names></name></person-group> (<year>1997</year>). <article-title>Vestibular-neck interaction and transformation of sensory coordinates</article-title>. <source>J. Vestib. Res.</source> <volume>7</volume>, <fpage>347</fpage>&#x02013;<lpage>367</lpage>. <pub-id pub-id-type="pmid">9218246</pub-id></citation></ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mergner</surname> <given-names>T.</given-names></name> <name><surname>Maurer</surname> <given-names>C.</given-names></name> <name><surname>Peterka</surname> <given-names>R. J.</given-names></name></person-group> (<year>2003</year>). <article-title>A multisensory posture control model of human upright stance</article-title>. <source>Prog. Brain Res.</source> <volume>142</volume>, <fpage>189</fpage>&#x02013;<lpage>201</lpage>. <pub-id pub-id-type="doi">10.1016/S0079-6123(03)42014-1</pub-id><pub-id pub-id-type="pmid">12693262</pub-id></citation></ref>
<ref id="B36">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Mergner</surname> <given-names>T.</given-names></name> <name><surname>Peterka</surname> <given-names>R. J.</given-names></name></person-group> (<year>2017</year>). <article-title>Human sense of balance</article-title>, in <source>Humanoid Robotics: A Reference</source>, eds <person-group person-group-type="editor"><name><surname>Goswami</surname> <given-names>A.</given-names></name> <name><surname>Vadakkepat</surname> <given-names>P.</given-names></name></person-group> (<publisher-name>Springer</publisher-name>).</citation></ref>
<ref id="B37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mergner</surname> <given-names>T.</given-names></name> <name><surname>Rosemeier</surname> <given-names>T.</given-names></name></person-group> (<year>1998</year>). <article-title>Interaction of vestibular, somatosensory and visual signals for postural control and motion perception under terrestrial and microgravity conditions&#x02014;a conceptual model</article-title>. <source>Brain Res. Rev.</source> <volume>28</volume>, <fpage>118</fpage>&#x02013;<lpage>135</lpage>. <pub-id pub-id-type="doi">10.1016/S0165-0173(98)00032-0</pub-id><pub-id pub-id-type="pmid">9795180</pub-id></citation></ref>
<ref id="B38">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mergner</surname> <given-names>T.</given-names></name> <name><surname>Schweigart</surname> <given-names>G.</given-names></name> <name><surname>Fennell</surname> <given-names>L.</given-names></name></person-group> (<year>2009</year>). <article-title>Vestibular humanoid postural control</article-title>. <source>J. Physiol. Paris</source> <volume>103</volume>, <fpage>178</fpage>&#x02013;<lpage>194</lpage>. <pub-id pub-id-type="doi">10.1016/j.jphysparis.2009.08.002</pub-id><pub-id pub-id-type="pmid">19665555</pub-id></citation></ref>
<ref id="B39">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Merton</surname> <given-names>P. A.</given-names></name></person-group> (<year>1953</year>). <article-title>Speculations on the servo control of movement</article-title>, in <source>The Spinal Cord</source>, eds <person-group person-group-type="editor"><name><surname>Malcom</surname> <given-names>J. L.</given-names></name> <name><surname>Gray</surname> <given-names>J. A. B.</given-names></name> <name><surname>Wolstenholme</surname> <given-names>G. E. W.</given-names></name></person-group> (<publisher-loc>Boston, MA</publisher-loc>: <publisher-name>Little Brown</publisher-name>), <fpage>247</fpage>&#x02013;<lpage>260</lpage>.</citation></ref>
<ref id="B40">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Morasso</surname> <given-names>P.</given-names></name></person-group> (<year>2013</year>). <article-title>What is the use of the body schema for humanoid robots?</article-title> <source>Int. J. Mach. Conscious</source> <volume>5</volume>, <fpage>75</fpage>&#x02013;<lpage>94</lpage>. <pub-id pub-id-type="doi">10.1142/S1793843013400064</pub-id></citation></ref>
<ref id="B41">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nashner</surname> <given-names>L. M.</given-names></name></person-group> (<year>1972</year>). <article-title>Vestibular postural control model</article-title>. <source>Kybernetik</source> <volume>10</volume>, <fpage>106</fpage>&#x02013;<lpage>110</lpage>. <pub-id pub-id-type="doi">10.1007/BF00292236</pub-id><pub-id pub-id-type="pmid">4537349</pub-id></citation></ref>
<ref id="B42">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nashner</surname> <given-names>L. M.</given-names></name> <name><surname>McCollum</surname> <given-names>G.</given-names></name></person-group> (<year>1985</year>). <article-title>The organization of human postural movements: a formal basis and experimental synthesis</article-title>. <source>Behav. Brain Sci.</source> <volume>8</volume>, <fpage>135</fpage>&#x02013;<lpage>150</lpage>. <pub-id pub-id-type="doi">10.1017/S0140525X00020008</pub-id></citation></ref>
<ref id="B43">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nenchev</surname> <given-names>D. N.</given-names></name></person-group> (<year>2013</year>). <article-title>Reaction null space of a multibody system with applications in robotics</article-title>. <source>Mech. Sci.</source> <volume>4</volume>, <fpage>97</fpage>&#x02013;<lpage>112</lpage>. <pub-id pub-id-type="doi">10.5194/ms-4-97-2013</pub-id></citation></ref>
<ref id="B44">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Nishio</surname> <given-names>A.</given-names></name> <name><surname>Takahashi</surname> <given-names>K.</given-names></name> <name><surname>Nenchev</surname> <given-names>D. N.</given-names></name></person-group> (<year>2006</year>). <article-title>Balance control of a humanoid robot based on the reaction null space method</article-title>, in <source>Intelligent Robots and Systems, 2006 IEEE/RSJ International Conference on 1996-2001</source> (<publisher-loc>Beijing</publisher-loc>).</citation></ref>
<ref id="B45">
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Nori</surname> <given-names>F.</given-names></name> <name><surname>Peters</surname> <given-names>J.</given-names></name> <name><surname>Padois</surname> <given-names>V.</given-names></name> <name><surname>Babic</surname> <given-names>J.</given-names></name> <name><surname>Ivaldi</surname> <given-names>S.</given-names></name></person-group> (<year>2014</year>). <article-title>Whole-body motion in humans and humanoids</article-title>, in <source>New Research Frontiers for Intelligent Autonomous Systems</source>, <fpage>81</fpage>&#x02013;<lpage>92</lpage>. Available online at: <ext-link ext-link-type="uri" xlink:href="http://hal.upmc.fr/hal-01053094">http://hal.upmc.fr/hal-01053094</ext-link></citation></ref>
<ref id="B46">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ott</surname> <given-names>C.</given-names></name> <name><surname>Henze</surname> <given-names>B.</given-names></name> <name><surname>Hettich</surname> <given-names>G.</given-names></name> <name><surname>Seyde</surname> <given-names>T. N. T. N.</given-names></name> <name><surname>Roa</surname> <given-names>M. A. M. A.</given-names></name> <name><surname>Lippi</surname> <given-names>V.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Good posture, good balance: Comparison of bioinspired and model-based approaches for posture control of humanoid robots</article-title>. <source>IEEE Rob. Autom. Mag.</source> <volume>23</volume>, <fpage>22</fpage>&#x02013;<lpage>33</lpage>. <pub-id pub-id-type="doi">10.1109/MRA.2015.2507098</pub-id></citation></ref>
<ref id="B47">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pasma</surname> <given-names>J. H.</given-names></name> <name><surname>Boonstra</surname> <given-names>T. A.</given-names></name> <name><surname>Campfens</surname> <given-names>S. F.</given-names></name> <name><surname>Schouten</surname> <given-names>A. C.</given-names></name> <name><surname>Van der Kooij</surname> <given-names>H.</given-names></name></person-group> (<year>2012</year>). <article-title>Sensory reweighting of proprioceptive information of the left and right leg during human balance control</article-title>. <source>J. Neurophysiol.</source> <volume>108</volume>, <fpage>1138</fpage>&#x02013;<lpage>1148</lpage>. <pub-id pub-id-type="doi">10.1152/jn.01008.2011</pub-id><pub-id pub-id-type="pmid">22623486</pub-id></citation></ref>
<ref id="B48">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Peterka</surname> <given-names>R. J.</given-names></name></person-group> (<year>2002</year>). <article-title>Sensorimotor integration in human postural control</article-title>. <source>J. Neurophysiol.</source> <volume>88</volume>, <fpage>1097</fpage>&#x02013;<lpage>1118</lpage>. <pub-id pub-id-type="doi">10.1152/jn.00605.2001</pub-id><pub-id pub-id-type="pmid">12205132</pub-id></citation></ref>
<ref id="B49">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Peterka</surname> <given-names>R. J.</given-names></name></person-group> (<year>2003</year>). <article-title>Simplifying the complexities of maintaining balance</article-title>. <source>Eng. Med. Biol. Mag. IEEE</source> <volume>22</volume>, <fpage>63</fpage>&#x02013;<lpage>68</lpage>. <pub-id pub-id-type="doi">10.1109/MEMB.2003.1195698</pub-id><pub-id pub-id-type="pmid">12733461</pub-id></citation></ref>
<ref id="B50">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Peterka</surname> <given-names>R. J.</given-names></name></person-group> (<year>2009</year>). <article-title>Comparison of human and humanoid robot control of upright stance</article-title>. <source>J. Physiol. Paris</source> <volume>103</volume>, <fpage>149</fpage>&#x02013;<lpage>158</lpage>. <pub-id pub-id-type="doi">10.1016/j.jphysparis.2009.08.001</pub-id><pub-id pub-id-type="pmid">19665564</pub-id></citation></ref>
<ref id="B51">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pruszynski</surname> <given-names>J. A.</given-names></name> <name><surname>Scott</surname> <given-names>S. H.</given-names></name></person-group> (<year>2012</year>). <article-title>Optimal feedback control and the long-latency stretch response</article-title>. <source>Exp. Brain Res.</source> <volume>218</volume>, <fpage>341</fpage>&#x02013;<lpage>359</lpage>. <pub-id pub-id-type="doi">10.1007/s00221-012-3041-8</pub-id><pub-id pub-id-type="pmid">22370742</pub-id></citation></ref>
<ref id="B52">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Roffel</surname> <given-names>B.</given-names></name> <name><surname>Betlem</surname> <given-names>B.</given-names></name></person-group> (<year>2006</year>). <source>Process Dynamics and Control: Modeling for Control and Prediction</source>. <publisher-loc>Hoboken, NJ</publisher-loc>: <publisher-name>John Wiley &#x00026; Sons</publisher-name>.</citation></ref>
<ref id="B53">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Schaal</surname> <given-names>S.</given-names></name> <name><surname>Peters</surname> <given-names>J.</given-names></name> <name><surname>Nakanishi</surname> <given-names>J.</given-names></name> <name><surname>Ijspeert</surname> <given-names>A.</given-names></name></person-group> (<year>2003</year>). <article-title>Control, planning, learning, and imitation with dynamic movement primitives</article-title>, in <source>Workshop on Bilateral Paradigms on Humans and Humanoids, IEEE International Conference on Intelligent Robots and Systems</source> (<publisher-loc>Las Vegas, NV</publisher-loc>), <fpage>1</fpage>&#x02013;<lpage>21</lpage>.</citation></ref>
<ref id="B54">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Schilling</surname> <given-names>M.</given-names></name></person-group> (<year>2011</year>). <article-title>Universally manipulable body models&#x02014;dual quaternion representations in layered and dynamic MMCs</article-title>. <source>Auton. Robots</source> <volume>30</volume>, <fpage>399</fpage>&#x02013;<lpage>425</lpage>. <pub-id pub-id-type="doi">10.1007/s10514-011-9226-3</pub-id></citation></ref>
<ref id="B55">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Schilling</surname> <given-names>M.</given-names></name> <name><surname>Cruse</surname> <given-names>H.</given-names></name></person-group> (<year>2008</year>). <article-title>The evolution of cognition&#x02014;from first order to second order embodiment</article-title>, in <source>Modeling Communication with Robots and Virtual Humans</source>, eds <person-group person-group-type="editor"><name><surname>Wachsmuth</surname> <given-names>I.</given-names></name> <name><surname>Knoblich</surname> <given-names>G.</given-names></name></person-group> (<publisher-loc>Berlin</publisher-loc>, <publisher-name>Springer</publisher-name>), <fpage>77</fpage>&#x02013;<lpage>108</lpage>.</citation></ref>
<ref id="B56">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Schweigart</surname> <given-names>G.</given-names></name> <name><surname>Mergner</surname> <given-names>T.</given-names></name></person-group> (<year>2008</year>). <article-title>Human stance control beyond steady state response and inverted pendulum simplification</article-title>. <source>Exp. Brain Res</source>. <volume>185</volume>, <fpage>635</fpage>&#x02013;<lpage>653</lpage>. <pub-id pub-id-type="doi">10.1007/s00221-007-1189-4</pub-id><pub-id pub-id-type="pmid">18030458</pub-id></citation></ref>
<ref id="B57">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sentis</surname> <given-names>L.</given-names></name> <name><surname>Khatib</surname> <given-names>O.</given-names></name></person-group> (<year>2005</year>). <article-title>Synthesis of whole-body behaviors through hierarchical control of behavioral primitives</article-title>. <source>Int. J. Hum. Robot.</source> <volume>2</volume>, <fpage>505</fpage>&#x02013;<lpage>518</lpage>. <pub-id pub-id-type="doi">10.1142/S0219843605000594</pub-id></citation></ref>
<ref id="B58">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Sentis</surname> <given-names>L.</given-names></name> <name><surname>Khatib</surname> <given-names>O.</given-names></name></person-group> (<year>2006</year>). <article-title>A whole-body control framework for humanoids operating in human environments</article-title>, in <source>Proceedings 2006 IEEE International Conference on Robotics and Automation</source> (<publisher-loc>Orlando, FL</publisher-loc>), <fpage>2641</fpage>&#x02013;<lpage>2648</lpage>.</citation></ref>
<ref id="B59">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Stephens</surname> <given-names>B.</given-names></name></person-group> (<year>2007</year>). <article-title>Humanoid push recovery</article-title>, in <source>7th IEEE-RAS International Conference on Humanoid Robots, HUMANOIDS 2007</source> (<publisher-loc>Pittsburgh, PA</publisher-loc>), <fpage>589</fpage>&#x02013;<lpage>595</lpage>.</citation></ref>
<ref id="B60">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Torricelli</surname> <given-names>D.</given-names></name> <name><surname>Mizanoor</surname> <given-names>R. S. M.</given-names></name> <name><surname>Gonzalez</surname> <given-names>J.</given-names></name> <name><surname>Lippi</surname> <given-names>V.</given-names></name> <name><surname>Hettich</surname> <given-names>G.</given-names></name> <name><surname>Asslaender</surname> <given-names>L.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Benchmarking human-like posture and locomotion of humanoid robots: A preliminary scheme</article-title>, in <source>Conference on Biomimetic and Biohybrid Systems</source> (<publisher-loc>Cham</publisher-loc>: <publisher-name>Springer</publisher-name>), <fpage>320</fpage>&#x02013;<lpage>331</lpage>.</citation></ref>
<ref id="B61">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>van der Kooij</surname> <given-names>H.</given-names></name> <name><surname>Jacobs</surname> <given-names>R.</given-names></name> <name><surname>Koopman</surname> <given-names>B.</given-names></name> <name><surname>van der Helm</surname> <given-names>F.</given-names></name></person-group> (<year>2001</year>). <article-title>An adaptive model of sensory integration in a dynamic environment applied to human stance control</article-title>. <source>Biol. Cybern.</source> <volume>84</volume>, <fpage>103</fpage>&#x02013;<lpage>115</lpage>. <pub-id pub-id-type="doi">10.1007/s004220000196</pub-id><pub-id pub-id-type="pmid">11205347</pub-id></citation></ref>
<ref id="B62">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>van der Kooij</surname> <given-names>H.</given-names></name> <name><surname>Peterka</surname> <given-names>R. J.</given-names></name></person-group> (<year>2011</year>). <article-title>Non-linear stimulus-response behavior of the human stance control system is predicted by optimization of a system with sensory and motor noise</article-title>. <source>J. Comp. Neurosci.</source> <volume>30</volume>, <fpage>759</fpage>&#x02013;<lpage>778</lpage>. <pub-id pub-id-type="doi">10.1007/s10827-010-0291-y</pub-id><pub-id pub-id-type="pmid">21161357</pub-id></citation></ref>
<ref id="B63">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Vukobratovi&#x00107;</surname> <given-names>M.</given-names></name> <name><surname>Borovac</surname> <given-names>B.</given-names></name></person-group> (<year>2004</year>). <article-title>Zero-moment point &#x02014; Thirty five years of its life</article-title>. <source>Int. J. of Humanoid Robot.</source> <volume>1</volume>, <fpage>157</fpage>&#x02013;<lpage>173</lpage>. <pub-id pub-id-type="doi">10.1142/S0219843604000083</pub-id></citation></ref>
<ref id="B64">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Wiener</surname> <given-names>N.</given-names></name></person-group> (<year>1948</year>). <source>Cybernetics: Control and Communication in the Animal and the Machine</source>. <publisher-loc>Cambridge, MA</publisher-loc>: <publisher-name>The Technology Press</publisher-name>; <publisher-loc>New York, NY</publisher-loc>: <publisher-name>John Wiley &#x00026; Sons, Inc</publisher-name>.</citation></ref>
<ref id="B65">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Winter</surname> <given-names>D. A.</given-names></name></person-group> (<year>1990</year>). <source>Biomechanics and Motor Control of Human Movement, 2nd Edn</source>. <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Wiley</publisher-name>.</citation></ref>
<ref id="B66">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Zhong</surname> <given-names>H.</given-names></name> <name><surname>Pao</surname> <given-names>L.</given-names></name> <name><surname>de Callafon</surname> <given-names>R.</given-names></name></person-group> (<year>2012</year>). <article-title>Feedforward control for disturbance rejection: model matching and other methods</article-title>, in <source>IEEE 24th Chinese Control and Decision Conference (CCDC)</source> (<publisher-loc>Taiyuan</publisher-loc>), <fpage>3528</fpage>&#x02013;<lpage>3533</lpage>.</citation></ref>
</ref-list>
<fn-group>
<fn id="fn0001"><p><sup>1</sup>PRTS stands for pseudo-random ternary sequence stimulus (see Peterka, <xref ref-type="bibr" rid="B48">2002</xref>). It allows evaluation of gain of the disturbance-evoked body excursion and its phase and coherence (e.g., with support surface tilt) over a defined frequency range. Data processing comprises a spectral analysis of the stimulus vs. the body (e.g., COM) angular excursions in space using a discrete Fourier transform. The pseudorandom sequence of the tilts makes them unpredictable for humans. See Hettich et al. (<xref ref-type="bibr" rid="B17">2014</xref>) for details of the stimulus used.</p></fn>
</fn-group>
<fn-group>
<fn fn-type="financial-disclosure"><p><bold>Funding.</bold> This research is supported by the by the European Commission (FP7 Grant 600698 H2R and 610454 EMBalance).</p>
</fn>
</fn-group>
</back>
</article>