Quantitative analysis of human ankle characteristics at different gait phases and speeds for utilizing in ankle-foot prosthetic design
BioMedical Engineering OnLine volume 13, Article number: 19 (2014)
Ankle characteristics vary in terms of gait phase and speed change. This study aimed to quantify the components of ankle characteristics, including quasi-stiffness and work in different gait phases and at various speeds.
The kinetic and kinematic data of 20 healthy participants were collected during normal gait at four speeds. Stance moment-angle curves were divided into three sub-phases including controlled plantarflexion, controlled dorsiflexion and powered plantarflexion. The slope of the moment-angle curves was quantified as quasi-stiffness. The area under the curves was defined as work.
The lowest quasi-stiffness was observed in the controlled plantarflexion. The fitted line to moment-angle curves showed R2 > 0.8 at controlled dorsiflexion and powered plantarflexion. Quasi-stiffness was significantly different at different speeds (P = 0.00). In the controlled dorsiflexion, the ankle absorbed energy; by comparison, energy was generated in the powered plantarflexion. A negative work value was recorded at slower speeds and a positive value was observed at faster speeds. Ankle peak powers were increased with walking speed (P = 0.00).
Our findings suggested that the quasi-stiffness and work of the ankle joint can be regulated at different phases and speeds. These findings may be clinically applicable in the design and development of ankle prosthetic devices that can naturally replicate human walking at various gait speeds.
The overall function of a human ankle during the stance phase of walking can be explored in terms of moment versus angle relation . This relation reveals an almost linear loop-shaped curve for the majority of stance phase at slow and normal walking speeds [1–4]. The curve shows a clockwise hysteresis loop at slower speeds, whereas a counter-clockwise loop is displayed at higher speeds . Considering the moment-angle curve, researchers suggested that human ankle can be replaced with a passive rotational spring-damper system at slow and normal speeds, but an augmented active mechanism is necessary at fast speeds [3, 5, 6].
The conventional passive prosthetic feet are limited to relatively rigid ankle mechanisms. This condition has prompted researchers to develop passive elastic prosthetic feet with energy storing and returning capabilities . Although energy-storing prostheses exhibit more advantages than conventional feet , both devices show the same biomechanical behaviors which deviate from normal ankle function [1, 5, 8]. These limitations can be decreased by designing the mechanical parameters of the prosthetic ankle based on the human ankle characteristics [9–11].
Human ankle characteristics, such as quasi-stiffness and work derived from the moment-angle curve, are keys to design the prosthetic feet [1, 3, 5, 6, 12]. Quasi-stiffness refers to the slope of the curve and defined as overall ankle resistance to motion [3, 13–15]. Previous studies indicated that ankle quasi-stiffness is an important design aspect which should be tuned based on gait sub-phases and speeds in a successful prosthetic design. The values of ankle quasi-stiffness can be used as a guideline for researchers to select and adjust the spring coefficient of the prostheses [3, 5, 16]. Ankle work is calculated from the area under the curve and quantifies the amount of energy absorbed or produced by the joint structures [2, 5]. In addition, the ankle displays negative work in mid-stance and positive work at terminal stance. The ankle work and quasi-stiffness varies throughout a gait cycle or in terms of gait speed [2, 5, 8, 17]. A successful prosthetic design should incorporate adaptive characteristics similar to those of the human ankle in a gait cycle or at different speeds [1, 3, 5, 6].
Previous studies explored the ankle moment-angle relation at different speeds and showed some aspects of ankle function [3, 5, 15, 18]. However, the findings of these studies cannot be directly applied in the design of ankle prostheses because a range of ankle characteristics at different gait phases and speeds have not been presented.
In this study, the ankle moment-angle curve was quantitatively analyzed using a motion capture system. Useful information in terms of quasi-stiffness and work was also analyzed and presented in distinct sub-phases of stance at a range of walking speeds that can be applied for prosthetic designs.
Participants and protocol
Twenty healthy volunteers, 14 males and 6 females, without history of foot or ankle and lower extremity orthopedic or neurological pathologies participated in this study. On average, the participants were 23.4 (SD 3.7) years old, 67.6 (SD 10.2) kg in weight, and 174.2 (SD 7.7) cm in height. All of the participants provided written informed consent approved by the Ethics Committee of the University of Social Welfare and Rehabilitation Sciences.
Analyses were performed in a standard gait laboratory which is equipped with five infrared cameras (Vicon 460, Vicon Motion System Ltd., UK) and two force plates (Kistler Instrument AG, Switzerland). Data were collected at a rate of 100 Hz. Reflective markers were placed on the anatomical landmarks according to the recommendation of Vicon ‘Plug-In-Gait’ marker set . This set included the ankle (lateral malleolus), toe (dorsum of the foot between first and second metatarsals), heel, tibia (one third distal), knee (lateral femoral condoyle), femur (one third distal) and ASIS (anterior superior iliac spines). A marker was placed on the sacrum to compute the average walking speed . Each participant was instructed to walk barefoot at four walking speeds including normal, slow, very slow and fast speeds. At first, they were instructed to walk at their self-selected normal speed. Then they were asked to walk slower and faster than their normal speeds. The trials were accepted when the foot was completely in contact with the platform .
The marker positions and ground reaction forces were processed using the Vicon Plug-In-Gait model (Workstation version 4.6) to derive the ankle angle (degree), moment (N m), and power (W). Ankle moments were computed by an inverse dynamics approach. Ankle power was defined as multiplying the ankle moment by ankle angular velocity. Moment and power were then normalized by each individual’s body mass (N m kg−1, W kg−1). Average walking speed was defined as the total displacement of the sacrum marker divided by time at a defined distance . Kinetic and kinematic data were low-pass filtered using a zero-lag, sixth-order Butterworth filter with a cut-off frequency of 10 Hz . Linear interpolation was applied to the data points to establish equal lengths of data sets . Zero ankle position was defined as the point where the foot segment was perpendicular to the tibia segment. The dorsiflexion angle and moment of the ankle were considered positive.
Each moment-angle loop in the stance phase was divided into three sub-phases: controlled plantarflexion, controlled dorsiflexion; and powered plantarflexion [12, 21]. The beginning of the controlled plantarflexion phase was characterized by a heel strike and ended at the maximum plantarflexion in the early stance. The subsequent controlled dorsiflexion phase was identified from the end of controlled plantarflexion to the maximum ankle dorsiflexion in mid-stance, in which the ankle angle increased and power remained negative. Powered plantarflexion was defined as the phase from the end of controlled dorsiflexion when ankle power changed to positive values to the toe-off. In this sub-phase, the ankle obtained the maximum plantarflexion and power remained positive. The quasi-stiffnesses or K of the ankle (N m kg−1 rad−1) in controlled plantarflexion (KCP), controlled dorsiflexion (KCD), and powered plantarflexion (KPP) were estimated as the slopes of the corresponding linear regression lines to the moment-angle data at each sub-phase (Equation 1) (Figure 1) [1–3, 14, 15, 18].
Where M is the ankle moment, θ is the joint angle and K stands for the quasi-stiffness.
The values of work exerted by the ankle (N m rad kg−1) in controlled plantarflexion (WCP), controlled dorsiflexion (WCD), and powered plantarflexion (WPP) were calculated as the area under moment-angle curve at each sub-phase on the basis of the trapezoidal approximation approach. The area within the loop was considered as the total work exerted by the ankle in the stance phase (Wtotal); this area was computed by subtracting the absolute areas under powered plantarflexion (|WPP|) and controlled dorsiflexion (|WCD|) curves [2, 3, 5].
The normal distribution of data was examined by Kolmogorov-Smirnov test. Linear regression models using least squares approach were applied to examine the relationship between moment and angle. The goodness of fit of the regression models was evaluated using the coefficients of determination (R2). Repeated measures ANOVA with Bonferroni-adjusted post-hoc test was conducted to compare the means at different sub-phases and speeds. Statistical analyses were performed using SPSS version 19 (SPSS Inc., Chicago, IL, USA). P < 0.05 was considered as the significance level of all of the tests.
The average walking speeds at the four categories were 0.95 (0.07), 1.06 (SD 0.09), 1.25 (SD 0.13), and 1.58 (SD 0.12) m s−1. The mean differences among the walking speeds in the four categories were statistically significant (P < 0.05).
The regression lines fitted to the controlled plantarflexion moment-angle data showed 0.49 ≤ R2 ≤ 0.75 (Table 1). No significant differences were found between KCPs at all of the speeds used in this study (P > 0.05). Moment curves versus angle curves in the controlled dorsiflexion phase were observed in the regression models, but R2 decreased as speed increased (0.83 ≤ R2 ≤ 0.96, Table 1). Statistically significant differences were found in KCD in terms of walking speeds. KCD at fast and normal speeds was significantly higher than that at slow and very slow speeds (P = 0.00, Table 1). The moment-angle curves of powered plantarflexion showed R2 > 0.95 (Table 1). The results showed no significant differences between KPPs at all of the walking speeds (Table 1).
Statistically significant differences were found in Wtotal in all speeds (P = 0.00, Table 2). Moreover, in early stance (controlled plantarflexion phase) and mid-stance to terminal stance (controlled dorsiflexion phase), power was negative, indicating energy absorption. In the pre-swing (powered plantarflexion) phase, power was positive, showing energy generation in this period (Figure 2). The mean ± SD of stance power at slow to fast speeds were −0.12 (0.13), 0.03 (0.15), 0.17 (0.18), 0.46 (0.22) W kg−1. The mean power differences were significant (P < 0.05). Maximum ankle power was increased with walking speed (Figure 3).
In the present study, a range of ankle quasi-stiffness and work was investigated in terms of three distinct stance sub-phases at four walking speeds. This data set can be beneficial for designers who define the mechanical requirements of prosthetic ankle-foot devices.
Findings showed that the general pattern of ankle behavior were consistent with previous studies, demonstrating that human ankle changes from a passive to an active system in response to increased walking speed [5, 18, 22].
In the controlled plantarflexion phase, findings showed lower values of quasi-stiffness and work than other sub-phases. In this period, ankle begins to plantarflex under the eccentric contraction of dorsiflexor muscles to provide shock absorption, and to control the rate of floor impact . Regression lines fitted to the moment versus angle data showed a great variability in this period. This might be due to the number of the kinetics and kinematics data in the controlled plantarflexion (about 5 data points at the frequency of 100 Hz) as this phase occurs in a short period of time. However, fits with high R2s values were seen in the majority of the participants. Therefore, spring-like behavior of ankle joint can be considered in the prosthetic design.
In the controlled dorsiflexion phase, ankle showed moment-angle relationship with acceptable R2 values. In this interval, tibia moves forward over the stationary foot, ankle begins and continues to dorsiflex with increasing plantarflexor moment that peaks at the end of the phase . Moreover, negative work is done by ankle plantarflexors to decelerate the rate of leg forward movement [3, 18, 21]. Thus, it could be simply assumed that ankle shows spring-like behavior with considerable amount of energy absorption in the controlled dorsiflexion phase. Considering the walking speed, results indicated that ankle becomes stiffer at higher speeds in the controlled dorsiflexion. Similarly, Shamaei et al.  concluded that movement of the ankle joint becomes more nonlinear in controlled dorsiflexion phase at higher speeds. It can be inferred that ankle acts possibly more linearly elastic in slower speeds, like a linear spring, but at higher speeds more complicated and augmented systems are added. Thus, a spring with an adjustable stiffness modulated according to the walking speed can be considered in the prosthetic design in this period.
During the powered plantarflexion interval, body weight abruptly transfers to the contra lateral limb. Ankle joint plantarflexes under a declining moment while plantarflexors produces a positive peak power to push the body forward . The results at the powered plantarflexion phase showed positive peak power at all speeds that were increased with walking speed. Hence, results suggest that in powered plantarflexion interval, ankle functions like a spring–damper structure at very slow speed. However, the damping effect gradually vanishes at slow speeds. At higher walking speeds, an active source of energy progressively turns up. This suggests the demands for implementation of augmented systems (such as motor) in combination with passive spring-damper elements at fast walking speeds.
This study was performed on 20 healthy young participants at different ranges of age, mobility, and walking speed. Therefore, the results of the analyses cannot be applied in other populations, such as those with lower limb amputees or diabetic patients. The kinetic and kinematic data of the ankle were derived using a single-segment Vicon Plug-In-Gait model; as such, this factor should be considered when our results are applied since power and work may be overestimated .
This work provides a quantitative understanding of ankle characteristics aimed to be assistive for designers in developing new ankle-foot prostheses. The values of the ankle quasi-stiffness and work presented in this study can be used as a guideline for the selection and adjustment of the spring coefficient of the prostheses. However, variations in the parameters exist at different sub-phases and walking speeds. These adjustment strategies should be considered in the design of ankle prosthetic systems to attain more natural behavior.
Quasi-stiffness of the ankle in controlled plantarflexion
Quasi-stiffness of the ankle in controlled dorsiflexion
Quasi-stiffness of the ankle in powered plantarflexion
Work exerted by the ankle in controlled plantarflexion
Work exerted by the ankle in controlled dorsiflexion
Work exerted by the ankle in powered plantarflexion.
Versluys R, Beyl P, Van Damme M, Desomer A, Van Ham R, Lefeber D: Prosthetic feet: State-of-the-art review and the importance of mimicking human ankle-foot biomechanics. Disabil Rehabil Assist Technol 2009,4(2):65–75. 10.1080/17483100802715092
Crenna P, Frigo C: Dynamics of the ankle joint analyzed through moment–angle loops during human walking: Gender and age effects. Hum Mov Sci 2011,30(6):1185–1198. 10.1016/j.humov.2011.02.009
Shamaei K, Sawicki GS, Dollar AM: Estimation of quasi-stiffness and propulsive work of the human ankle in the stance phase of walking. PLoS One 2013,8(3):e59935. 10.1371/journal.pone.0059935
Frigo C, Crenna P, Jensen LM: Moment-angle relationship at lower limb joints during human walking at different velocities. J Electromyogr Kinesiol 1996,6(3):177–190. 10.1016/1050-6411(96)00030-2
Hansen AH, Childress DS, Miff SC, Gard SA, Mesplay KP: The human ankle during walking: implications for design of biomimetic ankle prostheses. J Biomech 2004,37(10):1467–1474. 10.1016/j.jbiomech.2004.01.017
Takahashi KZ, Stanhope SJ: Mechanical energy profiles of the combined ankle–foot system in normal gait: insights for prosthetic designs. Gait Posture 2013,38(4):818–823. 10.1016/j.gaitpost.2013.04.002
Hsu M-J, Nielsen DH, Lin-Chan S-J, Shurr D: The effects of prosthetic foot design on physiologic measurements, self-selected walking velocity, and physical activity in people with transtibial amputation. Arch Phys Med Rehabil 2006,87(1):123–129. 10.1016/j.apmr.2005.07.310
Hansen AH, Miff SC, Childress DS, Gard SA, Meier MR: Net external energy of the biologic and prosthetic ankle during gait initiation. Gait Posture 2010,31(1):13–17. 10.1016/j.gaitpost.2009.08.237
Singer E, Ishai G, Kimmel E: Parameter estimation for a prosthetic ankle. Ann Biomed Eng 1995,23(5):691–696. 10.1007/BF02584466
Safaeepour Z, Esteki A, Tabatabai Ghomshe F, Mousavai ME: Design and development of a novel viscoelastic ankle-foot prosthesis based on the human ankle biomechanics. Prosthet Orthot Int 2014. in press
Eshraghi A, Abu Osman NA, Gholizadeh H, Ali S, Shadgan B: 100 top-cited scientific papers in limb prosthetics. Biomed Eng Online 2013,12(1):1–12. 10.1186/1475-925X-12-1
Au S, Berniker M, Herr H: Powered ankle-foot prosthesis to assist level-ground and stair-descent gaits. Neural Netw 2008,21(4):654–666. 10.1016/j.neunet.2008.03.006
Latash ML, Zatsiorsky VM: Joint stiffness: myth or reality? Hum Mov Sci 1993,12(6):653–92. 10.1016/0167-9457(93)90010-M
Butler RJ, Crowell Iii HP, Davis IM: Lower extremity stiffness: implications for performance and injury. Clin Biomech 2003,18(6):511–517. 10.1016/S0268-0033(03)00071-8
Gabriel RC, Abrantes J, Granata K, Bulas-Cruz J, Melo-Pinto P, Filipe V: Dynamic joint stiffness of the ankle during walking: gender-related differences. Phys Ther Sport 2008,9(1):16–24. 10.1016/j.ptsp.2007.08.002
Safaeepour Z, Esteki A, Mousavi M, Tabatabaei F: Quantitative analysis of the human ankle viscoelastic behavior at different gait speeds. In 5th Kuala Lumpur International Conference on Biomedical Engineering. Edited by: Osman N, Abas W, Wahab A, Ting H-N. Berlin Heidelberg: Springer; 2011:200–202.
Winter DA: Energy generation and absorption at the ankle and knee during fast, natural, and slow cadences. Clin Orthop Relat Res 1983, 175: 147–54.
Davis RB, DeLuca PA: Gait characterization via dynamic joint stiffness. Gait Posture 1996,4(3):224–231. 10.1016/0966-6362(95)01045-9
Davis RB, Õunpuu S, Tyburski D, Gage JR: A gait analysis data collection and reduction technique. Hum Mov Sci 1991,10(5):575–587. 10.1016/0167-9457(91)90046-Z
Winter DA: Biomechanics and motor control of human movement. 3rd edition. Toronto, Canada: Wiley; 2005.
Perry J, Burnfield JM: Gait Analysis: Normal and Pathological Function. 2nd edition. Thorofare, USA: Slack Inc; 2010.
Orendurff MS, Segal AD, Aiona MD, Dorociak RD: Triceps surae force, length and velocity during walking. Gait Posture 2005,21(2):157–163. 10.1016/j.gaitpost.2004.01.009
Bruening DA, Cooney KM, Buczek FL: Analysis of a kinetic multi-segment foot model part II: kinetics and clinical implications. Gait Posture 2012,35(4):535–540. 10.1016/j.gaitpost.2011.11.012
The authors would like to thank biomechanics laboratory at Department of Ergonomics, University of Social Welfare and Rehabilitation Sciences, Tehran, Iran. Contribution of Hoda Nabavi and Sara Soltani in data collection is appreciated. This study was supported by UM/MOHE/HIR Project No. D000014-16001.
The authors declare no competing interests.
ZS carried out the research, participated in the study design, data analysis and drafted the manuscript. AE and FTG involved in study design, revising the manuscript and analysis of data. NAAO revised it critically for important intellectual content. All authors read and approved the final manuscript.
About this article
Cite this article
Safaeepour, Z., Esteki, A., Ghomshe, F.T. et al. Quantitative analysis of human ankle characteristics at different gait phases and speeds for utilizing in ankle-foot prosthetic design. BioMed Eng OnLine 13, 19 (2014). https://doi.org/10.1186/1475-925X-13-19
- Ankle biomechanics
- Ankle-foot mechanisms