Contact stress distributions on the femoral head of the emu (Dromaius novaehollandiae)

Osteonecrosis of the femoral head remains a challenging orthopaedic problem. The disease frequently progresses to femoral head collapse, leading to debilitating osteoarthritis in the affected hip(s). Since a major goal of pre-collapse interventions is to forestall the need for hip arthroplasty, it is important that any animal models used to develop or study such interventions also have a natural history of progression to femoral head collapse. The emu (Dromaius novaehollandiae), a large flightless bird native to Australia, consistently progresses to femoral head collapse when osteonecrosis is experimentally induced cryogenically. Full biomechanical characterization of the demands this animal places on its hip is an important consideration in future usage of this model. This study reports in vitro measurement of the contact stress distributions on the emu femoral head during stance phase of the gait cycle, using Fuji pressure-sensitive film. Applied hip loadings were based upon ground reaction forces and hip flexion angles recorded in vivo. The contact stress data showed reasonable homology with the human hip, both in terms of stress magnitude and sites of habitual loading on the femoral head.     

Contributions of individual muscles to hip joint contact force in normal walking

The human hip joint withstands high contact forces during daily activity and is therefore susceptible to injury and structural deterioration over time. Knowledge of muscle-force contributions to hip joint loading may assist in the development of strategies to prevent and manage conditions such as osteoarthritis, femoro-acetabular impingement and fracture. The main aim of this study was to determine the contributions of individual muscles to hip contact force in normal walking. Muscle contributions to hip contact force were calculated based on a previously published dynamic optimization solution for normal walking, which provided the time histories of joint motion, ground reaction forces, and muscle forces during the stance and swing phases of gait. The force developed by each muscle plus its contribution to the ground reaction force were used to determine the muscle’s contribution to hip contact force. Muscles were the major contributors to hip contact force, with gravitational and centrifugal forces combined contributing less than 5% of the total contact force. Four muscles that span the hip – gluteus medius, gluteus maximus, iliopsoas, and hamstrings – contributed most significantly to the three components of the hip contact force and hip contact impulse (integral of hip contact force over time). Three muscles that do not span the hip – vasti, soleus, and gastrocnemius – also contributed substantially to hip joint loading. These results provide additional insight into lower-limb muscle function during walking and may also be relevant to studies of cartilage degeneration and bone remodelling at the hip.     

Musculoskeletal model choice influences hip joint load estimations during gait

The prevalence of musculoskeletal modeling studies investigating hip contact forces and the number of models used to conduct such investigations has increased in recent years. However, the consistency between models remain unknown and differences in model predicted hip contact forces between studies are difficult to distinguish from natural inter-individual differences. The purpose of this study was therefore to evaluate differences in hip joint contact forces during gait between four OpenSim models. These models included the generic models gait2392 and the Arnold Lower Limb Model, as well as the hip specific models hip2372 and London Lower Limb Model. Data from four individuals who have had a total hip replacement with instrumented hip implants performing slow, normal, and fast walking trials were taken from the HIP98 database to evaluate the various models effectiveness at estimating hip loads. Muscle forces were estimated using static optimization and hip contact forces were calculated using the JointReaction analysis in OpenSim. Results indicated that, for gait, the hip specific London Lower Limb Model consistently predicted peak push-off hip joint contact forces with lower magnitude and timing errors compared to the other models. Likewise, root mean square error values were lowest and correlation coefficients were highest for the London Lower Limb Model. These results suggest that the London Lower Limb Model is the most appropriate model for investigations focused on hip joint loading.     

Walking with increased ankle pushoff decreases hip muscle moments

In a simple bipedal walking model, an impulsive push along the trailing limb (similar to ankle plantar flexion) or a torque at the hip can power level walking. This suggests a tradeoff between ankle and hip muscle requirements during human gait. People with anterior hip pain may benefit from walking with increased ankle pushoff if it reduces hip muscle forces. The purpose of our study was to determine if simple instructions to alter ankle pushoff can modify gait dynamics and if resulting changes in ankle pushoff have an effect on hip muscle requirements during gait. We hypothesized that changes in ankle kinetics would be inversely related to hip muscle kinetics. Ten healthy subjects walked on a custom split-belt force-measuring treadmill at 1.25m/s. We recorded ground reaction forces and lower extremity kinematic data to calculate joint angles and internal muscle moments, powers and angular impulses. Subjects walked under three conditions: natural pushoff, decreased pushoff and increased pushoff. For the decreased pushoff condition, subjects were instructed to push less with their feet as they walked. Conversely, for the increased pushoff condition, subjects were instructed to push more with their feet. As predicted, walking with increased ankle pushoff resulted in lower peak hip flexion moment, power and angular impulse as well as lower peak hip extension moment and angular impulse (p<0.05). Our results emphasize the interchange between hip and ankle kinetics in human walking and suggest that increased ankle pushoff during gait may help to compensate for hip muscle weakness or injury and reduce hip joint forces.     

Musculo-skeletal loading conditions at the hip during walking and stair climbing.

Musculo-skeletal loading plays an important role in the primary stability of joint replacements and in the biological processes involved in fracture healing. However, current knowledge of musculo-skeletal loading is still limited. In the past, a number of musculo-skeletal models have been developed to estimate loading conditions at the hip. So far, a cycle-to-cycle validation of predicted musculo-skeletal loading by in vivo measurements has not been possible. The aim of this study was to determine the musculo-skeletal loading conditions during walking and climbing stairs for a number of patients and compare these findings to in vivo data.Following total hip arthroplasty, four patients underwent gait analysis during walking and stair climbing. An instrumented femoral prosthesis enabled simultaneous measurement of in vivo hip contact forces. On the basis of CT and X-ray data, individual musculo-skeletal models of the lower extremity were developed for each patient. Muscle and joint contact forces were calculated using an optimization algorithm. The calculated peak hip contact forces both over- and under-estimated the measured forces. They differed by a mean of 12% during walking and 14% during stair climbing.For the first time, a cycle-to-cycle validation of predicted musculo-skeletal loading was possible for walking and climbing stairs in several patients. In all cases, the comparison of in vivo measured and calculated hip contact forces showed good agreement.Thus, the authors consider the presented approach as a useful means to determine valid conditions for the analysis of prosthesis loading, bone modeling or remodeling processes around implants and fracture stability following internal fixation.     

Influence of altered gait patterns on the hip joint contact forces

Children who exhibit gait deviations often present a range of bone deformities, particularly at the proximal femur. Altered gait may affect bone growth and lead to deformities by exerting abnormal stresses on the developing bones. The objective of this study was to calculate variations in the hip joint contact forces with different gait patterns. Muscle and hip joint contact forces of four children with different walking characteristics were calculated using an inverse dynamic analysis and a static optimisation algorithm. Kinematic and kinetic analyses were based on a generic musculoskeletal model scaled down to accommodate the dimensions of each child. Results showed that for all the children with altered gaits both the orientation and magnitude of the hip joint contact force deviated from normal. The child with the most severe gait deviations had hip joint contact forces 30% greater than normal, most likely due to the increase in muscle forces required to sustain his crouched stance. Determining how altered gait affects joint loading may help in planning treatment strategies to preserve correct loading on the bone from a young age.     

Hip contact forces and gait patterns from routine activities.

In vivo loads acting at the hip joint have so far only been measured in few patients and without detailed documentation of gait data. Such information is required to test and improve wear, strength and fixation stability of hip implants. Measurements of hip contact forces with instrumented implants and synchronous analyses of gait patterns and ground reaction forces were performed in four patients during the most frequent activities of daily living. From the individual data sets an average was calculated. The paper focuses on the loading of the femoral implant component but complete data are additionally stored on an associated compact disc. It contains complete gait and hip contact force data as well as calculated muscle activities during walking and stair climbing and the frequencies of daily activities observed in hip patients. The mechanical loading and function of the hip joint and proximal femur is thereby completely documented. The average patient loaded his hip joint with 238% BW (percent of body weight) when walking at about 4 km/h and with slightly less when standing on one leg. This is below the levels previously reported for two other patients (Bergmann et al., Clinical Biomechanics 26 (1993) 969-990). When climbing upstairs the joint contact force is 251% BW which is less than 260% BW when going downstairs. Inwards torsion of the implant is probably critical for the stem fixation. On average it is 23% larger when going upstairs than during normal level walking. The inter- and intra-individual variations during stair climbing are large and the highest torque values are 83% larger than during normal walking. Because the hip joint loading during all other common activities of most hip patients are comparably small (except during stumbling), implants should mainly be tested with loading conditions that mimic walking and stair climbing.     

A new parametric approach for modeling hip forces during gait

An analytical parametric model was developed to estimate the natural biological variations in muscle forces and their effect on the hip forces subject only to physiological constraints and not predefined optimization criterion. Force predictions are based on the joint kinematics and kinetics of each subject, a previously published muscle model, and physiological constraints on the muscle force distributions. The model was used to determine the hip contact forces throughout the stance phase of gait of a subject with a total hip replacement (THR). The parametrically modeled peak hip force without antagonistic muscle activity varied from 2.7 to 3.2 Body Weights (mean 2.9 Body Weights), which agreed well with published in vivo measurements from instrumented THRs in other subjects. For every 10% increase in antagonistic activity, the mean peak hip force increased by 0.2 Body Weights. The parametric model allows one to examine the effect of specific muscle weaknesses or increased antagonistic muscle activity on the hip forces. The model also provides a tool for studying the effect of gait adaptations on hip forces, as predictions are based on each individual's gait data. Differences in peak forces between subjects can then be evaluated relative to the uncertainty in not knowing the precise muscle force distributions.     

Assessment of the influence of foot orthoses in the hip loading conditions during walking: a single case study

Despite their large clinical application, the understanding of the effects of foot orthoses on the lower limb kinematics and kinetics is limited. In this context, we propose an advanced musculoskeletal model to assess the influence of foot orthoses in the loading conditions within an osteoarthritic hip joint during gait. Experimental data are collected for a single pathological subject presenting a coxarthrosis (with and without orthoses), and a healthy subject during walking. An inverse dynamic approach coupled with an optimisation method evaluates the forces developed by 14 muscles and the hip contact reaction force. Contact reaction and muscular force magnitudes are closed whether the patient is walking with or without foot orthoses. Nevertheless, contact reaction amplitudes and orientations show differences in relation to those calculated for the healthy subject. The results obtained allow us to formulate some assumptions concerning the causes of coxarthrosis evolution and treatment.     

Center of mass mechanics of chimpanzee bipedal walking

Center of mass (CoM) oscillations were documented for 81 bipedal walking strides of three chimpanzees. Full‐stride ground reaction forces were recorded as well as kinematic data to synchronize force to gait events and to determine speed. Despite being a bent‐hip, bent‐knee (BHBK) gait, chimpanzee walking uses pendulum‐like motion with vertical oscillations of the CoM that are similar in pattern and relative magnitude to those of humans. Maximum height is achieved during single support and minimum height during double support. The mediolateral oscillations of the CoM are more pronounced relative to stature than in human walking when compared at the same Froude speed. Despite the pendular nature of chimpanzee bipedalism, energy recoveries from exchanges of kinetic and potential energies are low on average and highly variable. This variability is probably related to the poor phasic coordination of energy fluctuations in these facultatively bipedal animals. The work on the CoM per unit mass and distance (mechanical cost of transport) is higher than that in humans, but lower than that in bipedally walking monkeys and gibbons. The pronounced side sway is not passive, but constitutes 10% of the total work of lifting and accelerating the CoM. CoM oscillations of bipedally walking chimpanzees are distinctly different from those of BHBK gait of humans with a flat trajectory, but this is often described as “chimpanzee‐like” walking. Human BHBK gait is a poor model for chimpanzee bipedal walking and offers limited insights for reconstructing early hominin gait evolution. Am J Phys Anthropol 156:422–433, 2015. © 2014 Wiley Periodicals, Inc.     

An open source lower limb model: Hip joint validation

Musculoskeletal lower limb models have been shown to be able to predict hip contact forces (HCFs) that are comparable to in vivo measurements obtained from instrumented prostheses. However, the muscle recruitment predicted by these models does not necessarily compare well to measured electromyographic (EMG) signals. In order to verify if it is possible to accurately estimate HCFs from muscle force patterns consistent with EMG measurements, a lower limb model based on a published anatomical dataset (Klein Horsman et al., 2007. Clinical Biomechanics. 22, 239-247) has been implemented in the open source software OpenSim. A cycle-to-cycle hip joint validation was conducted against HCFs recorded during gait and stair climbing trials of four arthroplasty patients (Bergmann et al., 2001. Journal of Biomechanics. 34, 859-871). Hip joint muscle tensions were estimated by minimizing a polynomial function of the muscle forces. The resulting muscle activation patterns obtained by assessing multiple powers of the objective function were compared against EMG profiles from the literature. Calculated HCFs denoted a tendency to monotonically increase their magnitude when raising the power of the objective function; the best estimation obtained from muscle forces consistent with experimental EMG profiles was found when a quadratic objective function was minimized (average overestimation at experimental peak frame: 10.1% for walking, 7.8% for stair climbing). The lower limb model can produce appropriate balanced sets of muscle forces and joint contact forces that can be used in a range of applications requiring accurate quantification of both. The developed model is available at the website https://simtk.org/home/low_limb_london.     

Model of loadings acting on the femoral bone during gait

An evaluation of the model of loadings acting on the femoral bone during the whole gait cycle was the main aim of the paper. A computer simulation of the musculoskeletal system based on the gait data collected during gait was used to determine the muscle forces as well as the hip joint reaction. Kinematic parameters as well as the ground reaction force for ninety-nine healthy persons of both sexes (18–36 years old) who had no history of musculoskeletal disease were registered during normal gait with preferred speed and used as inputs for musculoskeletal modelling and numerical simulation with the use of the AnyBody software. Time waveforms of the values of force generated by 21 muscles having attachments on the femoral bone as well as the hip joint reaction force were obtained. Directions of particular forces were presented using a femoral coordinate system. Attachment points for all muscle forces were obtained on the basis of the unscaled standard model with the length of the femur equal to 0.41 m. The presented model of loadings acting on the femoral bone element can be useful for the biomechanical analysis of bone development and remodelling as well as for the optimisation of implant or bone stabilizer design and pre-clinical testing.     

Long-term hip loading in unilateral total hip replacement patients is no different between limbs or compared to healthy controls at similar walking speeds

Variation in hip joint contact forces directly influences the performance of total hip replacements (THRs). Measurement and calculation of contact forces in THR patients has been limited by small sample sizes, wide variation in patient and surgical factors, and short-term follow-up. This study hypothesised that, at long-term follow-up, unilateral THR patients have similar calculated hip contact forces compared to controls walking at similar (self-selected) speeds and, in contrast, THR patients walking at slower (self-selected) speeds have reduced hip contact forces. It was further hypothesised that there is no difference in calculated hip contact forces between operated and non-operated limbs at long-term follow-up for both faster and slower patients. Gait analysis data for THR patients walking at faster (walking speed: 1.29 ± 0.12 m/s; n = 11) and slower (walking speed: 0.72 ± 0.09 m/s; n = 11) speeds were used. Healthy subjects constituted the control group (walking speed: 1.36 ± 0.12 m/s; n = 10). Hip contact forces were calculated using static optimisation. There was no significant difference (p > 0.31) in hip contact forces between faster and control groups. Conversely, force was reduced at heel strike by 19% (p = 0.002), toe-off by 31% (p < 0.001) and increased at mid-stance by 15% (p = 0.02) for the slower group compared to controls. There were no differences between operated and non-operated limbs for the slower group or the faster group, suggesting good biomechanical recovery at long-term follow-up. Loading, at different walking speeds, presented here can improve the relevance of preclinical testing methods.     

Effect of arm swinging on lumbar spine and hip joint forces

During level walking, arm swing plays a key role in improving dynamic stability. In vivo investigations with a telemeterized vertebral body replacement showed that spinal loads can be affected by differences in arm positions during sitting and standing. However, little is known about how arm swing could influence the lumbar spine and hip joint forces and motions during walking. The present study aims to provide better understanding of the contribution of the upper limbs to human gait, investigating ranges of motion and joint reaction forces.A three-dimensional motion analysis was carried out via a motion capturing system on six healthy males and five patients with hip instrumented implant. Each subject performed walking with different arm swing amplitudes (small, normal, and large) and arm positions (bound to the body, and folded across the chest). The motion data were imported in a commercial musculoskeletal analysis software for kinematic and inverse dynamic investigation.The range of motion of the thorax with respect to the pelvis and of the pelvis with respect to the ground in the transversal plane were significantly associated with arm position and swing amplitude during gait. The hip external-internal rotation range of motion statistically varied only for non-dominant limb. Unlike hip joint reaction forces, predicted peak spinal loads at T12-L1 and L5-S1 showed significant differences at approximately the time of contralateral toe off and contralateral heel strike.Therefore, arm position and swing amplitude have a relevant effect on kinematic variables and spinal loads, but not on hip loads during walking.     

Influence of weak hip abductor muscles on joint contact forces during normal walking: probabilistic modeling analysis

The weakness of hip abductor muscles is related to lower-limb joint osteoarthritis, and joint overloading may increase the risk for disease progression. The relationship between muscle strength, structural joint deterioration and joint loading makes the latter an important parameter in the study of onset and follow-up of the disease. Since the relationship between hip abductor weakness and joint loading still remains an open question, the purpose of this study was to adopt a probabilistic modeling approach to give insights into how the weakness of hip abductor muscles, in the extent to which normal gait could be unaltered, affects ipsilateral joint contact forces. A generic musculoskeletal model was scaled to each healthy subject included in the study, and the maximum force-generating capacity of each hip abductor muscle in the model was perturbed to evaluate how all physiologically possible configurations of hip abductor weakness affected the joint contact forces during walking. In general, the muscular system was able to compensate for abductor weakness. The reduced force-generating capacity of the abductor muscles affected joint contact forces to a mild extent, with 50th percentile mean differences up to 0.5 BW (maximum 1.7 BW). There were greater increases in the peak knee joint loads than in loads at the hip or ankle. Gluteus medius, particularly the anterior compartment, was the abductor muscle with the most influence on hip and knee loads. Further studies should assess if these increases in joint loading may affect initiation and progression of osteoarthritis.     

Subject-specific geometrical detail rather than cost function formulation affects hip loading calculation*

This study assessed the relative importance of introducing an increasing level of medical image-based subject-specific detail in bone and muscle geometry in the musculoskeletal model, on calculated hip contact forces during gait. These forces were compared to introducing minimization of hip contact forces in the optimization criterion. With an increasing level of subject-specific detail, specifically MRI-based geometry and wrapping surfaces representing the hip capsule, hip contact forces decreased and were more comparable to contact forces measured using instrumented prostheses (average difference of 0.69 BW at the first peak compared to 1.04 BW for the generic model). Inclusion of subject-specific wrapping surfaces in the model had a greater effect than altering the cost function definition.     

Contributions of muscles to terminal-swing knee motions vary with walking speed

Many children with cerebral palsy walk with diminished knee extension during terminal swing, at speeds much slower than unimpaired children. Treatment of these gait abnormalities is challenging because the factors that extend the knee during normal walking, over a range of speeds, are not well understood. This study analyzed a series of three-dimensional, muscle-driven dynamic simulations to determine whether the relative contributions of individual muscles and other factors to angular motions of the swing-limb knee vary with walking speed. Simulations were developed that reproduced the measured gait dynamics of seven unimpaired children walking at self-selected, fast, slow, and very slow speeds (7 subjects×4 speeds=28 simulations). In mid-swing, muscles on the stance limb made the largest net contribution to extension of the swing-limb knee at all speeds examined. The stance-limb hip abductors, in particular, accelerated the pelvis upward, inducing reaction forces at the swing-limb hip that powerfully extended the knee. Velocity-related forces (i.e., Coriolis and centrifugal forces) also contributed to knee extension in mid-swing, though these contributions were diminished at slower speeds. In terminal swing, the hip flexors and other muscles on the swing-limb decelerated knee extension at the subjects’ self-selected, slow, and very slow speeds, but had only a minimal net effect on knee motions at the fastest speeds. Muscles on the stance limb helped brake knee extension at the subjects’ fastest speeds, but induced a net knee extension acceleration at the slowest speeds. These data—which show that the contributions of muscular and velocity-related forces to terminal-swing knee motions vary systematically with walking speed—emphasize the need for speed-matched control subjects when attempting to determine the causes of a patient's abnormal gait.     

Muscular coordination of knee motion during the terminal-swing phase of normal gait

Children with cerebral palsy often walk with diminished knee extension during the terminal-swing phase, resulting in a troublesome "crouched" posture at initial contact and a shortened stride. Treatment of this gait abnormality is challenging because the factors that extend the knee during normal walking are not well understood, and because the potential of individual muscles to limit terminal-swing knee extension is unknown. This study analyzed a series of three-dimensional, muscle-driven dynamic simulations to quantify the angular accelerations of the knee induced by muscles and other factors during swing. Simulations were generated that reproduced the measured gait dynamics and muscle excitation patterns of six typically developing children walking at self-selected speeds. The knee was accelerated toward extension in the simulations by velocity-related forces (i.e., Coriolis and centrifugal forces) and by a number of muscles, notably the vasti in mid-swing (passive), the hip extensors in terminal swing, and the stance-limb hip abductors, which accelerated the pelvis upward. Knee extension was slowed in terminal swing by the stance-limb hip flexors, which accelerated the pelvis backward. The hamstrings decelerated the forward motion of the swing-limb shank, but did not contribute substantially to angular motions of the knee. Based on these data, we hypothesize that the diminished knee extension in terminal swing exhibited by children with cerebral palsy may, in part, be caused by weak hip extensors or by impaired hip muscles on the stance limb that result in abnormal accelerations of the pelvis.     

Effect of walking speed and severity of hip osteoarthritis on gait variability

Gait analysis in orthopaedic and neurological examinations is important; however, few studies assess gait variability at different walking speeds in patients with varying degrees of hip osteoarthritis. We aimed to clarify (1) how different controlled speeds and (2) various severities of hip osteoarthritis influence gait variability. Gait variability was described by the standard deviation (SD) of the spatial–temporal and mean standard deviation (MeanSD) of angular parameters. The spatial positions of the anatomical points for calculating gait parameters were determined in 20 healthy elderly controls and 20 patients with moderate and 20 patients with severe hip osteoarthritis with a zebris CMS-HS ultrasound-based motion analysis system at three walking speeds. The SD of the spatial–temporal and MeanSD of angular parameters of gait, which together describe gait variability, significantly depended on speed and osteoarthritis severity. The lowest variability in the gait was found near the self-selected walking speeds. Hip joint degeneration significantly worsened variability on the affected side, with non-affected joints and the pelvis compensating by increasing flexibility and adapting to step-by-step motions. Particular attention must be paid to improving gait stability and the reliability of limb movements in the presence of and increasing severity of osteoarthritis.     

A mass-length scaling law for modeling muscle strength in the lower limb

Musculoskeletal computer models are often used to study muscle function in children with and without impaired mobility. Calculations of muscle forces depend in part on the assumed strength of each muscle, represented by the peak isometric force parameter, which is usually based on measurements obtained from cadavers of adult donors. The aim of the present study was twofold: first, to develop a method for scaling lower-limb peak isometric muscle forces in typically-developing children; and second, to determine the effect of this scaling method on model calculations of muscle forces obtained for normal gait. Muscle volumes were determined from magnetic resonance (MR) images obtained from ten children aged from 7 to 13yr. A new mass-length scaling law was developed based on the assumption that muscle volume and body mass are linearly related, which was confirmed by the obtained volume and body mass data. Two musculoskeletal models were developed for each subject: one in which peak isometric muscle forces were estimated using the mass-length scaling law; and another in which these parameters were determined directly from the MR-derived muscle volumes. Musculoskeletal modeling and quantitative gait analysis were then used to calculate lower-limb muscle forces in normal walking. The patterns of muscle forces predicted by the model with scaled peak isometric force values were similar to those predicted by the MR-based model, implying that assessments of muscle function obtained from these two methods are practically equivalent. These results support the use of mass-length scaling in the development of subject-specific musculoskeletal models of children.