# Do horizontal propulsive forces influence the nonlinear structure of locomotion?

- Max J Kurz
^{1}Email author and - Nicholas Stergiou
^{2}

**4**:30

**DOI: **10.1186/1743-0003-4-30

© Kurz and Stergiou; licensee BioMed Central Ltd. 2007

**Received: **13 October 2006

**Accepted: **15 August 2007

**Published: **15 August 2007

## Abstract

### Background

Several investigations have suggested that changes in the nonlinear gait dynamics are related to the neural control of locomotion. However, no investigations have provided insight on how neural control of the locomotive pattern may be directly reflected in changes in the nonlinear gait dynamics. Our simulations with a passive dynamic walking model predicted that toe-off impulses that assist the forward motion of the center of mass influence the nonlinear gait dynamics. Here we tested this prediction in humans as they walked on the treadmill while the forward progression of the center of mass was assisted by a custom built mechanical horizontal actuator.

### Methods

Nineteen participants walked for two minutes on a motorized treadmill as a horizontal actuator assisted the forward translation of the center of mass during the stance phase. All subjects walked at a self-select speed that had a medium-high velocity. The actuator provided assistive forces equal to 0, 3, 6 and 9 percent of the participant's body weight. The largest Lyapunov exponent, which measures the nonlinear structure, was calculated for the hip, knee and ankle joint time series. A repeated measures one-way analysis of variance with a t-test post hoc was used to determine significant differences in the nonlinear gait dynamics.

### Results

The magnitude of the largest Lyapunov exponent systematically increased as the percent assistance provided by the mechanical actuator was increased.

### Conclusion

These results support our model's prediction that control of the forward progression of the center of mass influences the nonlinear gait dynamics. The inability to control the forward progression of the center of mass during the stance phase may be the reason the nonlinear gait dynamics are altered in pathological populations. However, these conclusions need to be further explored at a range of walking speeds.

## Background

Human and animal locomotion is typically described as having a periodic movement pattern. For example, it can be readily observed that the legs oscillate to-and-fro with a limit cycle behavior that is similar to the pendulum motions of a clock [1, 2]. Any variations from this periodic pattern have traditionally been considered to be "noise" within the neuromuscular system [3, 4]. However, recent investigations have confirmed that the step-to-step variations that are present in gait may not be strictly noise. Rather these variations may have a deterministic structure [3, 5–12]. Several authors have noted that the structure of the nonlinear gait dynamics is influenced by the health of the neuromuscular system. These results imply that the observed changes in the nonlinear gait dynamics may be related to the organization of the nervous system for functional and stable gait [3, 5–9]. Although this seems plausible, no efforts have been made to explore what neural control strategies can govern the nonlinear gait dynamics. Such insight may lead to new clinical methods for assessing the health of the neuromuscular system, and may lead to new metrics that can be used to guide the rehabilitation of the neuromuscular system.

The control of human locomotion can be globally divided into the stance and swing phases. A major determinant of the stance phase is the ability of the neuromuscular system to redirect the center of mass forward and over the support limb for each step of the gait cycle [13, 14]. Proper neuromuscular control of the center of mass allows for the locomotive system to take advantage of the energy exchange that is associated with the inverted pendulum dynamics [13]. Possibly the neural control strategies that dictate the forward progression of the center of mass may also influence the nonlinear structure of human locomotion. In this investigation, we explored if the control of the forward progression of the center of mass during the stance phase may govern the nonlinear gait dynamics.

To test the prediction of our model, we built a mechanical horizontal actuator that assisted the forward motion of the center of mass during the stance phase. Based on the results of our simulations, we hypothesized that the magnitude of the largest Lyapunov exponent is dependent on neural control of the forward progression of the center of mass during the stance phase of gait.

## Methods

### Nonlinear Analysis Techniques

The following analysis techniques were used to quantify the nonlinear structure of the gait patterns of the computer simulations, and the complementary human experiment conducted in this investigation.

From the original time series (*i.e*., knee angle), the state space was reconstructed based on Taken's embedding theorem [22, 23]. The reconstruction process involved creating time-lagged copies of the original time series. Equation 1 presents the reconstructed state vector where **y**(t) was the reconstructed state vector, x(t) was the original time series data, and x(t-T_{i}) was time delay copies of x(t).

**y**(t) = [x(t), x(t-T_{1}), x(t-T_{2}),...].

_{i}) for creating the state vector was determined by estimating when information about the state of the dynamic system at x(t) was different from the information contained in its time-delayed copy using an average mutual information algorithm [22]. Equation 2 presents the average mutual information algorithm used in this investigation where T was the time delay, x(t) was the original data, x(t+T) was the time delay data, P(x(t), x(t+T)) was the joint probability for measurement of x(t) and x(t+T), P(x(t)) was the probability for measurement of x(t), and P(x(t+T)) was the probability for measurement of x(t+T).

The number of embedding dimensions of the data time series was calculated to unfold the dynamics of the system in an appropriate state space. An inappropriate number of embedding dimensions may result in a projection of the dynamics of the system that has orbital crossings in the state space that are due to false neighbors, and not the actual dynamics of the system [22, 23]. To unfold the state space, we systematically inspected x(t), and its neighbors in various dimensions (*e.g*., dimension = 1, 2, 3,...etc.). The appropriate embedding dimension was identified when the neighbors of the x(t) stopped being un-projected by the addition of further dimensions of the state vector. For example, the global false nearest neighbors algorithm compares the points in the attractor at a given dimension d_{E}

Y(t) = [x(t), x(t + T), x(t + 2T), ... x(t + (d_{E}-1) T)].

**y**^{NN}(t) = [x^{NN}(t), x^{NN}(t + T), x^{NN}(t + 2T), ... x^{NN}(t + (d_{E}-1) T)].

**y**(t) is the current point being considered, and

**y**

^{NN}(t) is the nearest neighbor. If the distance between the points at the next dimension (e.g., d

_{E+1}) is greater than the distance calculated at the current dimension (e.g., d

_{E}), then the point is considered a false neighbor and further embeddings are necessary to unfold the attractor. The percentage of false nearest neighbors was calculated at higher dimensions until the percent nearest neighbors dropped to zero (Figure 4). The embedding dimension that had zero percent false nearest neighbors was used to re-construct the attractor in an appropriate state space. Equation 5 presents a reconstructed state vector where d

_{E}was the number of embedding dimensions,

**y**(t) was the d

_{E}-1 dimensional state vector, x(t) was the original data, and T was the time delay.

**y**(t) = [x(t), x(t + T), x(t + 2T), ... x(t + (d_{E}-1) T)].

The Applied Nonlinear Dynamics software was used to calculate the time lags and embedding dimensions for the computer simulations and complementary human experiments.

*Chaos Data Analyzer*(American Institute of Physics) to numerically calculate the largest Lyapunov exponent.

### Walking Model

where *θ* was the angle of the stance leg, *φ* was the angle of the swing leg and $\dot{\theta}$, $\ddot{\theta}$ and $\ddot{\phi}$ were the respective time derivatives, *γ* is the angle of the walking surface, and t is time. Equation 6 represents the stance leg and equation 7 represents the swing leg. Derivations of the equations of motion for the walking model are detailed in Garcia *et al*. [17].

The governing equations were integrated using a modified version of *Matlab's* (MathWorks, Natick, MA) ODE45. The ODE45 was modified to integrate the equations of motion with a tolerance of 10^{-11}, and to stop integrating when the angle of the swing leg angle was twice as large as the stance leg angle (Equation 8).

*φ* - 2*θ* = 0

where "+" indicated the behavior of the model just after heel-contact, "-" indicated the behavior of the model just before heel-contact, and J represents an instantaneous toe-off impulse that is directed toward the center of mass. J was dimensionless and had a normalization factor M(g l)^{1/2}. Further details on the derivation of the transition equation are found in Kuo [24]. A toe-off impulse was included in our model because several experimental investigations with humans have demonstrated that the ankle joint is a major contributor for the forward progression of the center of mass during locomotion [25, 26].

Analyses of the locomotive patterns of the model were performed from 3,000 footfalls with the first 500 footfalls removed to be certain that the model converged to the given attractor. The step time intervals were used to classify the gait pattern of the walking model. The influence of the toe-off impulse on the nonlinear structure of the model's gait was explored by systematically increasing J in the transition equation while the ramp angle remained constant in the governing equations. For each simulation, the largest Lyapunov exponent was calculated for the respective step time interval time series to determine how the altered toe-off impulses influenced the model's nonlinear gait dynamics using a time lag of one and an embedding dimension of three [19].

### Experimental Procedures

Nineteen subjects (14 Females, 5 Males; Age = 25.89 ± 5 years; Weight = 665.8 ± 79.8 N; Height = 1.68 ± 0.06 m) volunteered to participate in this investigation. All subjects were in good health and free from any musculoskeletal injuries and disorders. The experimental protocol used in this investigation was approved by the University's Internal Review Board and all subjects provided written informed consent. All subjects had treadmill walking experience prior to participating in the experiment.

The subjects walked on the treadmill for two minutes at a self-selected pace while the mechanical horizontal actuator assisted the forward motion of the center of mass. The average walking speed that was used for all conditions was 1.01 ± 0.2 ms^{-1} (Cadence = 116 ± 8 steps/min). The horizontal force actuator supplied a force equal to 0%, 3%, 6% and 9% of the subject's body weight. These percentages were selected based on our pilot data, where we determined that they provided the minimal distortion of the normal gait pattern. A high-speed digital four camera motion capture system (Motion Analysis, Santa Rosa, California) was used to capture the three dimensional positions of reflective markers placed on the lower extremity at 60 Hz. Triangulations of markers were placed on the thigh, shank and foot segments. A standing calibration was used to correct for misalignment of the markers with the local coordinate system of each of the lower extremity segments. This was accomplished by having the subjects stand in a calibration fixture that was aligned with the global reference system. Custom software was used to calculate the three-dimensional segment and joint angles consistent with Vaughan *et al*. [29] from the corrected positions of the segment markers. The joint angle time series were analyzed unfiltered in order to get a more accurate representation of the variability within the system [30]. Previous investigations have indicated that filtering the data may eliminate important information and provide a skewed view of the system's inherent variability [31]. Using the nonlinear analysis techniques discussed in Section A of the methods, the largest Lyapunov exponents for the respective joint angle time series were numerically calculated with an embedding dimension of six.

A one-way analysis of variances (ANOVA) with repeated-measures design was performed for each joint to determine statistical significance between the means of the respective assistance conditions. Furthermore, we used dependent t-tests with a Bonferroni adjustment as a post-hoc test to analyze if the respective assistance conditions were different from the no assistance condition. The alpha level was defined as *P* < 0.05. A linear trend analysis was performed if statistical differences were found. The trend analysis allowed us to infer if the nonlinear structure of the human gait pattern scaled in a similar fashion as the passive dynamic walking model computer simulations.

## Results and Discussion

### Simulation Results

Our simulations indicated that systematic increase in the toe-off impulse (J>0) resulted in the largest Lyapunov exponent to have a greater magnitude. For example, at a ramp angle of 0.0185 radians the largest Lyapunov exponent for the model's nonlinear gait pattern was 0.285 when J = 0. However, if a toe-off impulse was used to assist the forward progression of the center of mass (J = 0.001), the largest Lyapunov exponent of the model's gait pattern increased to a value of 0.363. These results are further detailed in Figure 2 where it is apparent that the largest Lyapunov exponent's magnitude linearly increased as a greater toe-off impulse was used to assist the forward progression of the model's center of mass. Therefore, the simulations predict that an increase in the propulsive forces that govern the forward translation of the center of mass during the stance phase will result in a linear increase in the magnitude of the largest Lyapunov exponent in a human's gait pattern.

### Experimental Results

The experimental results are consistent with the hypothesis that the nonlinear structure of gait is dependent on the neural control of the forward progression of the center of mass during the stance phase of gait. As the mechanical actuator increased the amount of assistance supplied to the center of mass, the magnitude of the largest Lyapunov exponent systematically increased for the hip and ankle joints. These results imply that the performance of the hip and ankle joints during the stance phase may be related to the changes in the nonlinear structure noted in previous investigations [6, 7, 9]. This is consistent with previous experimental studies where it has been concluded that the stance phase dynamics are dependent on the ankle and hip joints' control properties. The ankle joint supplies a large amount of power for the forward progression of the center of mass [25, 26] and the hip joint stabilizes the trunk during the early and late portions of the stance phase [32, 33]. However, it cannot be completely concluded if the changes in the nonlinear structure of the hip joint were a result of normal torso control during the stance phase. Since the mechanical horizontal actuator was attached at the waist of the subject, it may have artificially created instabilities in the torso which required an altered control strategy at the hip joint that would not have been present if the center of mass was actuated purely by a toe-off impulse.

The nonlinear structure of the knee joint during the horizontal assistance conditions was not significantly different from normal walking. The lack of clear results for the knee joint may be related to its functional role during gait. The behavior of the knee joint is largely attributed to maintaining the inverted pendulum during stance and limb clearance during the swing [32]. Hence, the knee joint has less influence on the forward progression of the center of mass [14]. However, further inspection of Figure 6 indicates that with the exception of 0%, the knee follows the same increasing linear trend as the ankle and hip joint. Possibly, the knee joint's nonlinear behavior may be also sensitive to the assistive force provided during the stance phase. However, further exploration of this notion is necessary before we can make this conclusion. Possibly, by altering the walking velocity of the subject, the linear trend at the knee joint may be further magnified.

The experiments conducted here were only performed at a medium-high walking velocity. This walking velocity may not be representative of the walking velocity that a disabled subject may select. Since we did not test the influence of horizontal assistance at a wide range of speeds we cannot generalize our results to all populations. Future investigations should explore how the interactive effect of walking speed and forward progression of the center of mass on the nonlinear structure of gait. These insights may lead to new insights on the nature of nonlinear gait patterns and may guide the development of rehabilitative protocols that are aimed at restoring a healthy nonlinear gait.

The passive dynamic walking model was able to predict the changes in the nonlinear structure of human locomotion as the forward progression of the center of mass was assisted. Although this model is highly simplified compared to the human locomotive system, it appears that it provides a well suited template for modeling the control properties of nonlinear gait dynamics. The additions of more life-like properties to this model may prove fruitful for the future research that is directed toward understanding how the neuromuscular properties influence the nonlinear structure of human locomotion. Such simulations and models will provide further insight on what neuromechanical variables influence the nonlinear gait dynamics.

## Conclusion

Horizontal propulsive forces that are applied during the stance phase influence the nonlinear structure of human locomotion. The experimental results presented here infer that the changes in the nonlinear structure may be related to the proper utilization the hip and ankle joint musculature to control the forward progression of the center of mass. Future investigation should determine if the results presented here can be extended to individuals with altered nonlinear gait patterns (*i.e*., elderly, Parkinson's disease). The initial step toward making this connection should be directed towards determining if the results presented here are consistent for different walking speeds. This scientific information will provide further insight on which neuromechanical components that are responsible for changes in the nonlinear structure of gait, and may lead to a better understanding of why the nonlinear gait pattern is altered in pathological populations.

## Declarations

### Acknowledgements

Funding was provided by the Nebraska Research Initiative Grant awarded to NS and the Texas Learning and Computational Center grant awarded to MJK.

## Authors’ Affiliations

## References

- Clark JE, Phillips SJ:
**A longitudinal study of intralimb coordination in the first year of independent walking: a dynamical systems analysis.***Child Devel*1987,**64:**1143-57. 10.2307/1131331View ArticleGoogle Scholar - Kurz MJ, Stergiou N:
**Applied dynamic systems theory for the analysis of movement.**In*Innovative Analyses of Human Movement*. Edited by: Stergiou N. Human Kinetics: Champaign, Illinois; 2004.Google Scholar - Hausdorff JM, Peng CK, Ladin Z, Wei JY, Goldberger AL:
**Iswalking a random walk? Evidence for long-range correlations in stride interval of human gait.***J Appl Physiol*1995,**78**(1):349-58.PubMedGoogle Scholar - Riley MA, Turvey MT:
**Variability and determinism in motor behavior.***J Mot Behav*2002,**34**(2):99-125.View ArticlePubMedGoogle Scholar - Dingwell JB, Cusumano JP:
**Nonlinear time series analysis of normal and pathological human walking.***Chaos*2000,**10**(4):848-863. 10.1063/1.1324008View ArticlePubMedGoogle Scholar - Buzzi UH, Stergiou N, Kurz MJ, Hageman PA, Heidel J:
**Nonlinear dynamics indicates aging affects variability during gait.***Clin Biomech*2003,**18:**435-43. 10.1016/S0268-0033(03)00029-9View ArticleGoogle Scholar - Hausdorff JM, Mitchell SL, Firtion R, Peng CK, Cudkowicz ME, Wei JY, Goldberger AL:
**Altered fractal dynamics of gait: reduced stride-interval correlations with aging and Huntington's disease.***J Appl Phyiol*1997,**82**(1):262-269.Google Scholar - Miller DJ, Stergiou N, Kurz MJ:
**An improved surrogate method for detecting the presence of chaos in gait.***J Biomech*2006,**39**(15):2873-6. 10.1016/j.jbiomech.2005.10.019View ArticlePubMedGoogle Scholar - Stergiou N, Buzzi UH, Kurz MJ, Heidel J:
**Nonlinear tools in human movement.**In*Innovative Analysis of Human Movement*. Edited by: Stergiou N. Champaign, Illinois: Human Kinetics; 2004.Google Scholar - Stergiou N, Moraiti C, Giakas G, Ristanis S, Georgoulis AD:
**The effect of the walking speed on the stability of the anterior cruciate ligament deficient knee.***Clin Biomech*2004,**19**(9):957-63. 10.1016/j.clinbiomech.2004.06.008View ArticleGoogle Scholar - West BJ, Griffin L:
**Allometric control of human gait.***Fractals*1998,**6**(2):101-108. 10.1142/S0218348X98000122View ArticleGoogle Scholar - West BL, Griffin L:
**Allometric control, inverse power laws and human gait.***Chaos, Solitons and Fractals*1999,**10**(9):1519-1527. 10.1016/S0960-0779(98)00149-0View ArticleGoogle Scholar - Cavagna GA, Heglund NC, Taylor RC:
**Mechanical work in terrestrial locomotion: two basic mechanisms for minimizing energy expenditure.***Am J Physiol*1977,**233**(5):R243-R261.PubMedGoogle Scholar - Kuo AD, Donelan JM, Ruina A:
**Energetic consequences of walking like an inverted pendulum: step-to-step transitions.***Exerc Sport Sci Rev*2005,**33**(2):88-97. 10.1097/00003677-200504000-00006View ArticlePubMedGoogle Scholar - Glass L:
**Synchoronization and rhythmic processes in physiology.***Nature*2001,**410:**277-284. 10.1038/35065745View ArticlePubMedGoogle Scholar - Full RJ, Koditschek DE:
**Templates and anchors: neuromechanical hypotheses of legged locomotion on land.***J Exp Biol*1999,**202:**3325-3332.PubMedGoogle Scholar - Garcia M, Chatterjee A, Ruina A, Coleman M:
**The simplest walking model: stability, complexity, and scaling.***ASME J Biomech Eng*1998,**120**(2):281-288.View ArticleGoogle Scholar - Goswami A, Thuilot B, Espiau B: A study of the passivegait of a compass-like biped robot: symmetry and chaos. Int J Robotics Res 1998.,17(12):Google Scholar
- Kurz MJ, Stergiou N, Heidel J, Foster T:
**A template for the exploration of chaotic locomotive patterns.***Chaos, Solitons and Fractals*2005,**23:**485-493. 10.1016/j.chaos.2004.04.034View ArticleGoogle Scholar - Kurz MJ, Stergiou N:
**An artificial neural network that utilizes hip joint actuations to control bifurcations and chaos in a passive dynamic bipedal walking model.***Biol Cybern*2005,**93**(3):213-21. 10.1007/s00422-005-0579-6View ArticlePubMedGoogle Scholar - Kurz MJ, Stergiou N:
**Hip actuations can be used to control bifurcations and chaos in a passive dynamic walking model.***ASME J Biomech Eng*2007,**129**(2):216-22. 10.1115/1.2486008View ArticleGoogle Scholar - Kantz H, Schreiber T:
*Nonlinear Time Series Analysis*. second edition. Cambridge University Press: New York; 2004.Google Scholar - Abarbanel HDI:
*Analysis of Observed Chaotic Data*. New York: Springer; 1996.View ArticleGoogle Scholar - Kuo AD:
**Energetics of actively powered locomotion using the simplest walking model.***ASME J Biomech Eng*2002,**124:**113-120. 10.1115/1.1427703View ArticleGoogle Scholar - Gottschall JS, Kram R:
**Energy cost and muscular activity required for leg swing during walking.***J Appl Physiol*2004,**99:**23-30. 10.1152/japplphysiol.01190.2004View ArticleGoogle Scholar - Winter DA:
**Energy generation and absorption at the ankle and knee during fast, natural, slow cadences.***Clinical Orthopedics*1983,**175:**147-154.Google Scholar - Chang YH, Kram R:
**Metabolic cost of generating horizontal forces during human running.***J Appl Physiol*1999,**86**(5):1657-1662.PubMedGoogle Scholar - Chang YH, Hamerski CM, Kram R:
**Applied horizontal force increases impact loading in reduced-gravity running.***J Biomech*2001,**34:**679-685. 10.1016/S0021-9290(00)00196-2View ArticlePubMedGoogle Scholar - Vaughan CL, Davis BL, O'Connor JC:
*Dynamics of Human Gait*. Champaign, Illinois: Human Kinetics; 1992.Google Scholar - Mees AI, Judd K:
**Dangers of geometric filtering.***Physica D Nonlinear Phenomena*1993,**68:**427-436. 10.1016/0167-2789(93)90135-NView ArticleGoogle Scholar - Rapp PE:
**A guide to dynamical analysis.***Integr Physiol Behav Sci*1994,**29**(3):311-327.View ArticlePubMedGoogle Scholar - Perry J:
*Gait Analysis Normal and Pathological Function*. SLACK Incorporated: Thorofare, New Jersey; 1992.Google Scholar - Winter DA:
**Kinematic and kinetic patterns in human gait: variability and compensating effects.***Human Mov Sci*1984,**3:**51-76. 10.1016/0167-9457(84)90005-8View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.