Forward dynamic simulation of bipedal walking in the Japanese macaque: investigation of causal relationships among limb kinematics, speed, and energetics of bipedal locomotion in a nonhuman primate

Japanese macaques that have been trained for monkey performances exhibit a remarkable ability to walk bipedally. In this study, we dynamically reconstructed bipedal walking of the Japanese macaque to investigate causal relationships among limb kinematics, speed, and energetics, with a view to understanding the mechanisms underlying the evolution of human bipedalism. We constructed a two-dimensional macaque musculoskeletal model consisting of nine rigid links and eight principal muscles. To generate locomotion, we used a trajectory-tracking control law, the reference trajectories of which were obtained experimentally. Using this framework, we evaluated the effects of changes in cycle duration and gait kinematics on locomotor efficiency. The energetic cost of locomotion was estimated based on the calculation of mechanical energy generated by muscles. Our results demonstrated that the mass-specific metabolic cost of transport decreased as speed increased in bipedal walking of the Japanese macaque. Furthermore, the cost of transport in bipedal walking was reduced when vertical displacement of the hip joint was virtually modified in the simulation to be more humanlike. Human vertical fluctuations in the body's center of mass actually contributed to energy savings via an inverted pendulum mechanism.     

Energetic costs of bipedal and quadrupedal walking in Japanese macaques

We investigated the energetic costs of quadrupedal and bipedal walking in two Japanese macaques. The subjects were engaged in traditional bipedal performance for years, and are extremely adept bipeds. The experiment was conducted in an airtight chamber with a gas analyzer. The subjects walked quadrupedally and bipedally at fixed velocities (<5 km/hr) on a treadmill in the chamber for 2.5-6 min. We estimated energy consumption from carbon dioxide (CO2) production. While walking bipedally, energetic expenditure increased by 30% relative to quadrupedalism in one subject, and by 20% in another younger subject. Energetic costs increased linearly with velocity in quadrupedalism and bipedalism, with bipedal/quadrupedal ratios remaining almost constant. Our experiments were relatively short in duration, and thus the observed locomotor costs may include presteady-state high values. However, there was no difference in experimental duration between bipedal and quadrupedal trials. Thus, the issue of steady state cannot cancel the difference in energetic costs. Furthermore, we observed that switching of locomotor mode (quadrupedalism to bipedalism) during a session resulted in a significant increase of CO2 production. Taylor and Rowntree ([1973] Science 179:186-187) noted that the energetic costs for bipedal and quadrupedal walking were the same in chimpanzees and capuchin monkeys. Although the reason for this inconsistency is not clear, species-specific differences should be considered regarding bipedal locomotor energetics among nonhuman primates. Extra costs for bipedalism may not be great in these macaques. Indeed, it is known that suspensory locomotion in Ateles consumes 1.3-1.4 times as much energy relative to quadrupedal progression. This excess ratio surpasses the bipedal/quadrupedal energetic ratios in these macaques.     

Dynamic optimization of human walking

A three-dimensional, neuromusculoskeletal model of the body was combined with dynamic optimization theory to simulate normal walking on level ground. The body was modeled as a 23 degree-of-freedom mechanical linkage, actuated by 54 muscles. The dynamic optimization problem was to calculate the muscle excitation histories, muscle forces, and limb motions subject to minimum metabolic energy expenditure per unit distance traveled. Muscle metabolic energy was calculated by slimming five terms: the basal or resting heat, activation heat, maintenance heat, shortening heat, and the mechanical work done by all the muscles in the model. The gait cycle was assumed to be symmetric; that is, the muscle excitations for the right and left legs and the initial and terminal states in the model were assumed to be equal. Importantly, a tracking problem was not solved. Rather only a set of terminal constraints was placed on the states of the model to enforce repeatability of the gait cycle. Quantitative comparisons of the model predictions with patterns of body-segmental displacements, ground-reaction forces, and muscle activations obtained from experiment show that the simulation reproduces the salient features of normal gait. The simulation results suggest that minimum metabolic energy per unit distance traveled is a valid measure of walking performance.     

Effect of posture and locomotion on energy expenditure

Energy expenditure for human adults and infants and for dogs was measured in resting (supine or lateral) posture, in bipedal posture and locomotion, and in quadrupedal posture and locomotion. Variations in respiratory and heart rate and in body temperature were utilized in this comparative study. Oxygen consumption was also measured in human adults. In human adults, bipedal posture and locomotion were shown to be much less energy-consuming than corresponding quadrupedal posture and locomotion. The opposite was observed in adult dogs, where bipedalism was shown to be much more energy-consuming than quadrupedalism. In addition, this study demonstrated, for human adults in their natural erect posture, an energy expenditure barely higher than in supine or lateral resting posture, while the dogs in their natural quadrupedal stance, the energy expenditure is much higher than in their resting posture. With respect to energy, therefore, humans are more adapted to bipedalism than dogs to quadrupedalism. Human children, at the transitional stage between quadrupedalism and bipedalism, have high and almost equal requirements for all postures and locomotions. This demonstrates, in term of energy, their incomplete adaptation to erect behavior.     

The morphology of the gluteus maximus during human evolution: Prerequisite or consequence of the upright bipedal posture?

The human gluteus maximus differs from that of the other hominoids because of its size and bony attachments. These differences raise questions concerning their sequence of appearance in human evolution. Given that humans practice a unique locomotor style, one wonders if the human gluteus maximus morphology is a prerequisite or a consequence of upright bipedal locomotion. This question is addressed using a computer model that evaluates muscle leverage in a variety of locomotor postures. In this model, the human-like, or ape-like, muscular pattern is imposed upon a representative hindlimb of each of the five extant hominoids. Shapes of the skeletal elements (i.e. ilium and ischium lengths) are adjusted in the computer to simulate an evolutionary progression from an ape to a human skeletal morphology. Changes in the leverage of different parts of the gluteus maximus (measured as moment arms) are monitored during this transition. The results show how the mechanical leverages of the gluteus maximus would have changed in a variety of hypothetical evolutionary sequences that describe an ape to human transition. Although the hominoid models exhibit minor differences in these simulations, they all show that the postural and locomotor functions of the gluteus maximus would become more difficult if musculoskeletal morphology changed to the human-like pattern before erect bipedal posture was adopted. Conversely, small adjustments in the ape-like musculoskeletal condition support an erect bipedal posture. These results suggest that a human like posture would have preceded the appearance of the human-like musculoskeletal morphology. Human gluteal morphology, therefore, is a consequence and not a prerequisite of the upright bipedal posture.     

Force, work and power output of lower limb muscles during human maximal-effort countermovement jumping.

The purpose of this study was to simulate human maximal-effort countermovement jumping with a three-dimensional neuromusculoskeletal model. The specific aim was to investigate muscle force, work and power output of major lower limb muscles during the motion. A neuromusculoskeletal model that has nine rigid body segments, 20 degrees of freedom, 32 Hill-type lower limb muscles was developed. The neural activation input signal was represented by a series of step functions with step duration of 0.05 s. The excitation-contraction dynamics of the contractile element, the tissues around the joints to limit the joint range of motion, as well as the foot-ground interaction were implemented. A simulation was started from a standing posture. Optimal pattern of the activation input signal was searched through numerical optimization with a goal of maximizing the height reached by the mass center of body after jumping up. As a result, feasible kinematics, ground reaction force profile and muscle excitation profile were generated. It was found that monoarticular muscles had major contributions of mechanical work and power output, whereas biarticular muscles had minor contributions. Hip adductors, abductors and external rotator muscles were vigorously activated, although their mechanical work and power output was minor because of their limited length change during the motion. Joint flexor muscles such as m. iliopsoas, m. biceps femoris short head and m. tibialis anterior were activated in the beginning of the motion with an effect of facilitating the generation of a countermovement.     

Speed modulation in hylobatid bipedalism: a kinematic analysis

Gibbons are highly arboreal apes, and it is expected that their bipedal locomotion will show some particularities related to the arboreal environment. Previous research has shown that, during hylobatid bipedalism, unsupported phases are rare and stride frequencies are relatively low. This study confirms previous findings, and we suggest that low stride frequencies and the absence of unsupported phases are ways to reduce disadvantageous branch oscillations during arboreal travel. Despite these restrictions, gibbons are able to locomote at a wide range of speeds, implying that they likely exploit other mechanisms to modulate their locomotor speed. To investigate this possibility, we collected video images of a large number of spontaneous bipedal bouts of four untrained white-handed gibbons by using an instrumented walkway with four synchronized cameras. These video images were digitized to obtain a quantification of the 3D kinematics of hylobatid bipedalism. We defined a large number of spatiotemporal and kinematic gait variables, and the relationship between these gait variables and (dimensionless) speed was statistically tested. It was found that gibbons mainly increase stride length to increase their locomotor speed; the main speed-modulating mechanisms are hip and ankle excursion and coupled knee and ankle extension at toe-off. Although aerial phases are rare, gibbons generally adopt a bipedal bouncing gait at most speeds and a clear-cut gait transition, as seen in human locomotion, is absent. Comparison with human and bonobo bipedalism showed that the variability of the 3D joint angles of the hind limb are comparable during human and gibbon bipedalism, and much lower than during bonobo bipedalism. The low variability found in gibbons might be related to constraints imposed by the arboreal environment. These arboreal constraints clearly affect the bipedal gait characteristics of gibbons, but do not constrain the ability to adopt a bipedal bouncing gait during terrestrial locomotion.     

Theories of bipedal walking: an odyssey

In this paper six theories of bipedal walking, and the evidence in support of the theories, are reviewed. They include: evolution, minimising energy consumption, maturation in children, central pattern generators, linking control and effect, and robots on two legs. Specifically, the six theories posit that: (1) bipedalism is the fundamental evolutionary adaptation that sets hominids--and therefore humans--apart from other primates; (2) locomotion is the translation of the centre of gravity along a pathway requiring the least expenditure of energy; (3) when a young child takes its first few halting steps, his or her biomechanical strategy is to minimise the risk of falling; (4) a dedicated network of interneurons in the spinal cord generates the rhythm and cyclic pattern of electromyographic signals that give rise to bipedal gait; (5) bipedal locomotion is generated through global entrainment of the neural system on the one hand, and the musculoskeletal system plus environment on the other; and (6) powered dynamic gait in a bipedal robot can be realised only through a strategy which is based on stability and real-time feedback control. The published record suggests that each of the theories has some measure of support. However, it is important to note that there are other important theories of locomotion which have not been covered in this review. Despite such omissions, this odyssey has explored the wide spectrum of bipedal walking, from its origins through to the integration of the nervous, muscular and skeletal systems.     

Evaluating functional roles of phase resetting in generation of adaptive human bipedal walking with a physiologically based model of the spinal pattern generator

The central pattern generators (CPGs) in the spinal cord strongly contribute to locomotor behavior. To achieve adaptive locomotion, locomotor rhythm generated by the CPGs is suggested to be functionally modulated by phase resetting based on sensory afferent or perturbations. Although phase resetting has been investigated during fictive locomotion in cats, its functional roles in actual locomotion have not been clarified. Recently, simulation studies have been conducted to examine the roles of phase resetting during human bipedal walking, assuming that locomotion is generated based on prescribed kinematics and feedback control. However, such kinematically based modeling cannot be used to fully elucidate the mechanisms of adaptation. In this article we proposed a more physiologically based mathematical model of the neural system for locomotion and investigated the functional roles of phase resetting. We constructed a locomotor CPG model based on a two-layered hierarchical network model of the rhythm generator (RG) and pattern formation (PF) networks. The RG model produces rhythm information using phase oscillators and regulates it by phase resetting based on foot-contact information. The PF model creates feedforward command signals based on rhythm information, which consists of the combination of five rectangular pulses based on previous analyses of muscle synergy. Simulation results showed that our model establishes adaptive walking against perturbing forces and variations in the environment, with phase resetting playing important roles in increasing the robustness of responses, suggesting that this mechanism of regulation may contribute to the generation of adaptive human bipedal locomotion.     

Muscle force production during bent-knee,bent-hip walking in humans

Researchers have long debated the locomotor posture used by the earliest bipeds. While many agree that by 3–4 Ma (millions of years ago), hominins walked with an extended-limb human style of bipedalism, researchers are still divided over whether the earliest bipeds walked like modern humans, or walked with a more bent-knee, bent-hip (BKBH) ape-like form of locomotion. Since more flexed postures are associated with higher energy costs, reconstructing early bipedal mechanics has implications for the selection pressures that led to upright walking. The purpose of this study is to determine how modern human anatomy functions in BKBH walking to clarify the links between morphology and energy costs in different mechanical regimes. Using inverse dynamics, we calculated muscle force production at the major limb joints in humans walking in two modes, both with extended limbs and BKBH. We found that in BKBH walking, humans must produce large muscle forces at the knee to support body weight, leading to higher estimated energy costs. However, muscle forces at the hip remained similar in BKBH and extended limb walking, suggesting that anatomical adaptations for hip extension in humans do not necessarily diminish the effective mechanical advantage at the hip in more flexed postures. We conclude that the key adaptations for economical walking, regardless of joint posture, seem to center on maintaining low muscle forces at the hip, primarily by keeping low external moments at the hip. We explore the implications of these results for interpreting locomotor energetics in early hominins, including australopithecines and Ardipithecus ramidus.     

Voluntary bipedal walking of infant chimpanzees

The voluntary bipedal walking of infant chimpanzees was studied by the analysis of foot force and by motion analysis. The infants were trained to locomote on a level platform without any restrictions on the locomotor pattern. The voluntary bipedal walking was compared with the other types of locomotion at the same age and with the trained bipedal walking performed by other chimpanzees, including adult chimpanzees. The characteristics of voluntary bipedal walking in the infant until one year of age were: (1) high-speed walking with short cycle duration; (2) short stance phase duration; (3) small braking component of the preceding leg and large acceleration of the following leg; (4) one downward peak in the vertical component; and (5) a relatively small transverse component. Bipedal walking usually continued for less than one second and ended in quadrupedal locomotion. During walking, the preceding foot touched the floor, heel first, as in the case of older chimpanzees and humans. At this age, bipedal walking was similar to high-speed locomotion. The voluntary bipedal walking of the two-year-old and frour-yearold chimpanzees was characterized as follows: (1) slower speed than during quadrupedal locomotion, (2) relatively long periods and distances; (3) well balanced accelerating and braking components; and (4) a vertical component showing two downward peaks and a trough in between during numerous trials. The last characteristic means that the body center of gravity is higher in the single stance phase, just as in the bipedal walkinbg of the adult chimpanzees and humans. The bipedal walking of infant chimpanzees was discussed in comparison with the walking of humans, including infants.     

Locomotion in ornithischian dinosaurs: an assessment using three‐dimensional computational modelling

Ornithischian dinosaurs were primitively bipedal with forelimbs modified for grasping, but quadrupedalism evolved in the clade on at least three occasions independently. Outside of Ornithischia, quadrupedality from bipedal ancestors has only evolved on two other occasions, making this one of the rarest locomotory transitions in tetrapod evolutionary history. The osteological and myological changes associated with these transitions have only recently been documented, and the biomechanical consequences of these changes remain to be examined. Here, we review previous approaches to understanding locomotion in extinct animals, which can be broadly split into form–function approaches using analogy based on extant animals, limb‐bone scaling, and computational approaches. We then carry out the first systematic attempt to quantify changes in locomotor muscle function in bipedal and quadrupedal ornithischian dinosaurs. Using three‐dimensional computational modelling of the major pelvic locomotor muscle moment arms, we examine similarities and differences among individual taxa, between quadrupedal and bipedal taxa, and among taxa representing the three major ornithischian lineages (Thyreophora, Ornithopoda, Marginocephalia). Our results suggest that the ceratopsid Chasmosaurus and the ornithopod Hypsilophodon have relatively low moment arms for most muscles and most functions, perhaps suggesting poor locomotor performance in these taxa. Quadrupeds have higher abductor moment arms than bipeds, which we suggest is due to the overall wider bodies of the quadrupeds modelled. A peak in extensor moment arms at more extended hip angles and lower medial rotator moment arms in quadrupeds than in bipeds may be due to a more columnar hindlimb and loss of medial rotation as a form of lateral limb support in quadrupeds. We are not able to identify trends in moment arm evolution across Ornithischia as a whole, suggesting that the bipedal ancestry of ornithischians did not constrain the development of quadrupedal locomotion via a limited number of functional pathways. Functional anatomy appears to have had a greater effect on moment arms than phylogeny, and the differences identified between individual taxa and individual clades may relate to differences in locomotor performance required for living in different environments or for clade‐specific behaviours.     

Trunk orientation causes asymmetries in leg function in small bird terrestrial locomotion

In contrast to the upright trunk in humans, trunk orientation in most birds is almost horizontal (pronograde). It is conceivable that the orientation of the heavy trunk strongly influences the dynamics of bipedal terrestrial locomotion. Here, we analyse for the first time the effects of a pronograde trunk orientation on leg function and stability during bipedal locomotion. For this, we first inferred the leg function and trunk control strategy applied by a generalized small bird during terrestrial locomotion by analysing synchronously recorded kinematic (three-dimensional X-ray videography) and kinetic (three-dimensional force measurement) quail locomotion data. Then, by simulating quail gaits using a simplistic bioinspired numerical model which made use of parameters obtained in in vivo experiments with real quail, we show that the observed asymmetric leg function (left-skewed ground reaction force and longer leg at touchdown than at lift-off) is necessary for pronograde steady-state locomotion. In addition, steady-state locomotion becomes stable for specific morphological parameters. For quail-like parameters, the most common stable solution is grounded running, a gait preferred by quail and most of the other small birds. We hypothesize that stability of bipedal locomotion is a functional demand that, depending on trunk orientation and centre of mass location, constrains basic hind limb morphology and function, such as leg length, leg stiffness and leg damping.     

The impact of locomotor energetics on mammalian foraging

The ecological impact of energy expended on an activity stems from its effect on foraging requirements. For locomotion, the effect of moving each additional unit distance probably depends on the proportional increase in energy expenditure. Other common measures of the cost of locomotion do not reflect the impact of energy expenditure on foraging requirements. In terrestrial mammals, both the effect of body mass and the unit cost itself are very small: moving one kilometre requires less than 2% of all other energy expenditures combined. Total locomotor energy expenditure ranges from 1/2% of daily metabolism for a 10 g mammal to 6% for an elephant. Potential sources of bias in the estimation of these costs include systematic bias in estimates of distance traversed and extra energy required for non-linear locomotion. Because larger mammals can readily locomote at greater speeds, the primary locomotor advantage of large size may not be conservation of energy but of time, which can mean greater safety and more or better food.     

A Wider Pelvis Does Not Increase Locomotor Cost in Humans,with Implications for the Evolution of Childbirth

The shape of the human female pelvis is thought to reflect an evolutionary trade-off between two competing demands: a pelvis wide enough to permit the birth of large-brained infants, and narrow enough for efficient bipedal locomotion. This trade-off, known as the obstetrical dilemma, is invoked to explain the relative difficulty of human childbirth and differences in locomotor performance between men and women. The basis for the obstetrical dilemma is a standard static biomechanical model that predicts wider pelves in females increase the metabolic cost of locomotion by decreasing the effective mechanical advantage of the hip abductor muscles for pelvic stabilization during the single-leg support phase of walking and running, requiring these muscles to produce more force. Here we experimentally test this model against a more accurate dynamic model of hip abductor mechanics in men and women. The results show that pelvic width does not predict hip abductor mechanics or locomotor cost in either women or men, and that women and men are equally efficient at both walking and running. Since a wider birth canal does not increase a woman’s locomotor cost, and because selection for successful birthing must be strong, other factors affecting maternal pelvic and fetal size should be investigated in order to help explain the prevalence of birth complications caused by a neonate too large to fit through the birth canal.     

Three-dimensional data-tracking dynamic optimization simulations of human locomotion generated by direct collocation

The aim of this study was to perform full-body three-dimensional (3D) dynamic optimization simulations of human locomotion by driving a neuromusculoskeletal model toward in vivo measurements of body-segmental kinematics and ground reaction forces. Gait data were recorded from 5 healthy participants who walked at their preferred speeds and ran at 2 m/s. Participant-specific data-tracking dynamic optimization solutions were generated for one stride cycle using direct collocation in tandem with an OpenSim-MATLAB interface. The body was represented as a 12-segment, 21-degree-of-freedom skeleton actuated by 66 muscle-tendon units. Foot-ground interaction was simulated using six contact spheres under each foot. The dynamic optimization problem was to find the set of muscle excitations needed to reproduce 3D measurements of body-segmental motions and ground reaction forces while minimizing the time integral of muscle activations squared. Direct collocation took on average 2.7 ± 1.0 h and 2.2 ± 1.6 h of CPU time, respectively, to solve the optimization problems for walking and running. Model-computed kinematics and foot-ground forces were in good agreement with corresponding experimental data while the calculated muscle excitation patterns were consistent with measured EMG activity. The results demonstrate the feasibility of implementing direct collocation on a detailed neuromusculoskeletal model with foot-ground contact to accurately and efficiently generate 3D data-tracking dynamic optimization simulations of human locomotion. The proposed method offers a viable tool for creating feasible initial guesses needed to perform predictive simulations of movement using dynamic optimization theory. The source code for implementing the model and computational algorithm may be downloaded at http://simtk.org/home/datatracking.