Subject-specific hip geometry affects predicted hip joint contact forces during gait

Hip loading affects bone remodeling and implant fixation. In this study, we have analyzed the effect of subject-specific modeling of hip geometry on muscle activation patterns and hip contact forces during gait, using musculoskeletal modeling, inverse dynamic analysis and static optimization. We first used sensitivity analysis to analyze the effect of isolated changes in femoral neck-length (NL) and neck-shaft angle (NSA) on calculated muscle activations and hip contact force during the stance phase of gait. A deformable generic musculoskeletal model was adjusted incrementally to adopt a physiological range of NL and NSA. In a second similar analysis, we adjusted hip geometry to the measurements from digitized radiographs of 20 subjects with primary hip osteoarthrosis. Finally, we studied the effect of hip abductor weakness on muscle activation patterns and hip contact force. This analysis showed that differences in NL (41-74 mm) and NSA (113-140 degrees ) affect the muscle activation of the hip abductors during stance phase and hence hip contact force by up to three times body weight. In conclusion, the results from both the sensitivity and subject-specific analysis showed that at the moment of peak contact force, altered NSA has only a minor effect on the loading configuration of the hip. Increased NL, however, results in an increase of the three hip contact-force components and a reduced vertical loading. The results of these analyses are essential to understand modified hip joint loading, and for planning hip surgery for patients with osteoarthrosis.     

Subject-specific hip geometry and hip joint centre location affects calculated contact forces at the hip during gait

Hip loading affects the development of hip osteoarthritis, bone remodelling and osseointegration of implants. In this study, we analyzed the effect of subject-specific modelling of hip geometry and hip joint centre (HJC) location on the quantification of hip joint moments, muscle moments and hip contact forces during gait, using musculoskeletal modelling, inverse dynamic analysis and static optimization. For 10 subjects, hip joint moments, muscle moments and hip loading in terms of magnitude and orientation were quantified using three different model types, each including a different amount of subject-specific detail: (1) a generic scaled musculoskeletal model, (2) a generic scaled musculoskeletal model with subject-specific hip geometry (femoral anteversion, neck-length and neck-shaft angle) and (3) a generic scaled musculoskeletal model with subject-specific hip geometry including HJC location. Subject-specific geometry and HJC location were derived from CT. Significant differences were found between the three model types in HJC location, hip flexion–extension moment and inclination angle of the total contact force in the frontal plane. No model agreement was found between the three model types for the calculation of contact forces in terms of magnitude and orientations, and muscle moments. Therefore, we suggest that personalized models with individualized hip joint geometry and HJC location should be used for the quantification of hip loading. For biomechanical analyses aiming to understand modified hip joint loading, and planning hip surgery in patients with osteoarthritis, the amount of subject-specific detail, related to bone geometry and joint centre location in the musculoskeletal models used, needs to be considered.     

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.     

Sensitivity analysis of hip joint centre estimation based on three-dimensional CT scans

In morphological analysis of the femur, the hip joint centre (HJC) is generally determined using a 3D model of the femoral head based on medical images. However, the portion of the image selected to represent the femoral head may influence the HJC. We determined if this influence invalidates the results of three HJC calculation methods, one of which we introduce here.

To isolate femoral heads in cadaver CT images, thresholds were applied to the distance between femur and acetabulum models. The sensitivity of the HJC to these thresholds and the differences between methods were quantified.

For thresholds between 6 and 9 mm and healthy hips, differences between methods were below 1 mm and all methods were insensitive to threshold changes. For higher thresholds, the fovea capitis femoris disturbed the HJC. In two deformed hips, the new method performed superiorly. We conclude that for normal hips all methods produce valid results.     

Influence of the method of TM joint total replacement implantation on the loading of the joint on the opposite side

The temporomandibular (TM) joint is one of the most used joints in the human body, and any defect in this joint has a significant influence on quality of life. The objective of this study was to create a parametric numerical finite element (FE) analysis to compare the effect of surgical techniques used for total TM joint replacement implantation on loading the TM joint on the other side. Our hypothesis is that for the optimal function of all total TM joint replacements used in clinical practice it is crucial to devise a minimally invasive surgical technique, whereby there is minimum resection of masticatory muscles. This factor is more important than the design of the usually used total TM joint replacements. The extent of muscle resection influences the mechanical loading of the whole system. In the parametric FE analyses, the magnitude of the TM joint loading was compared for four different ranges of muscle resections during bite, using an anatomical model. The results obtained from all FE analyses support our hypothesis that an increasing extent of the muscle resection increased the magnitude of the TM joint overloading on the opposite side. The magnitude of the TM joint overloading increased depending on the muscle resection to 235% for bite on an incisor and up to 491% for bite on molars. Our study leads to a recommendation that muscle resection be minimised during replacement implantation and to a proposal that the attachment of the condylar part of the TM joint replacement be modified.     

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.     

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.     

Three-dimensional knee joint contact forces during walking in unilateral transtibial amputees

Individuals with unilateral transtibial amputations have greater prevalence of osteoarthritis in the intact knee joint relative to the residual leg and non-amputees, but the cause of this greater prevalence is unclear. The purpose of this study was to compare knee joint contact forces and the muscles contributing to these forces between amputees and non-amputees during walking using forward dynamics simulations. We predicted that the intact knee contact forces would be higher than those of the residual leg and non-amputees. In the axial and mediolateral directions, the intact and non-amputee legs had greater peak tibio-femoral contact forces and impulses relative to the residual leg. The peak axial contact force was greater in the intact leg relative to the non-amputee leg, but the stance phase impulse was greater in the non-amputee leg. The vasti and hamstrings muscles in early stance and gastrocnemius in late stance were the largest contributors to the joint contact forces in the non-amputee and intact legs. Through dynamic coupling, the soleus and gluteus medius also had large contributions, even though they do not span the knee joint. In the residual leg, the prosthesis had large contributions to the joint forces, similar to the soleus in the intact and non-amputee legs. These results identify the muscles that contribute to knee joint contact forces during transtibial amputee walking and suggest that the peak knee contact forces may be more important than the knee contact impulses in explaining the high prevalence of intact leg osteoarthritis.     

The effect on conventional gait model kinematics and kinetics of hip joint centre equations in adult healthy gait

The Conventional Gait Model (CGM) needs to benefit from large investigations on localization of the hip joint centre (HJC). Incorrect positions from the native equations were demonstrated (Sangeux et al., 2014; Harrington et al., 2007). More accurate equations were proposed but their impact on kinematics and kinetic CGM outputs was never evaluated. This short communication aims at examining if adoption of new HJC equations would alter standard CGM outputs. Sixteen able bodied participants underwent a full 3-D optoelectronic gait analysis followed by a 3-D ultrasound localization of their hips. Data were processed through the open source python package pyCGM2 replicating kinematic and kinetic processing of the native CGM. Compared with 3D ultrasound location, Hara equations improved the accuracy of sagittal plane kinematics (0.6°) and kinetics (0.02 N m kg⁻¹) for the hip. The worst case participant exhibited Harrington’s equations reached a deviation of 3° for the sagittal kinematics. In the coronal plane, Hara and Harrington equations presented similar differences (1°) for the hip whilst Davis equations had the largest deviation for hip abduction (2.7°) and hip abductor moment (0.10 N m kg⁻¹).Both Harrington and Hara equations improved the CGM location of the HJC. Hara equations improved results in the sagittal plane, plus utilise a single anthropometrics measurement, leg length, that may be more robust. However, neither set of equations had significant effect on kinematics. We reported some effects on kinetics, particularly in the coronal plane, which warrant caution in interpreting outputs using different sets of equations.     

Subject-specific knee joint geometry improves predictions of medial tibiofemoral contact forces

Estimating tibiofemoral joint contact forces is important for understanding the initiation and progression of knee osteoarthritis. However, tibiofemoral contact force predictions are influenced by many factors including muscle forces and anatomical representations of the knee joint. This study aimed to investigate the influence of subject-specific geometry and knee joint kinematics on the prediction of tibiofemoral contact forces using a calibrated EMG-driven neuromusculoskeletal model of the knee. One participant fitted with an instrumented total knee replacement walked at a self-selected speed while medial and lateral tibiofemoral contact forces, ground reaction forces, whole-body kinematics, and lower-limb muscle activity were simultaneously measured. The combination of generic and subject-specific knee joint geometry and kinematics resulted in four different OpenSim models used to estimate muscle–tendon lengths and moment arms. The subject-specific geometric model was created from CT scans and the subject-specific knee joint kinematics representing the translation of the tibia relative to the femur was obtained from fluoroscopy. The EMG-driven model was calibrated using one walking trial, but with three different cost functions that tracked the knee flexion/extension moments with and without constraint over the estimated joint contact forces. The calibrated models then predicted the medial and lateral tibiofemoral contact forces for five other different walking trials. The use of subject-specific models with minimization of the peak tibiofemoral contact forces improved the accuracy of medial contact forces by 47% and lateral contact forces by 7%, respectively compared with the use of generic musculoskeletal model.     

Simulated hip abductor strengthening reduces peak joint contact forces in patients with total hip arthroplasty

Lower extremity muscle strength training is a focus of rehabilitation following total hip arthroplasty (THA). Strength of the hip abductor muscle group is a predictor of overall function following THA. The purpose of this study was to investigate the effects of hip abductor strengthening following rehabilitation on joint contact forces (JCFs) in the lower extremity and low back during a high demand step down task. Five THA patients performed lower extremity maximum isometric strength tests and a stair descent task. Patient-specific musculoskeletal models were created in OpenSim and maximum isometric strength parameters were scaled to reproduce measured pre-operative joint torques. A pre-operative forward dynamic simulation of each patient performing the stair descent was constructed using their corresponding patient-specific model to predict JCFs at the ankle, knee, hip, and low back. The hip abductor muscles were strengthened with clinically supported increases (0–30%) above pre-operative values in a probabilistic framework to predict the effects on peak JCFs (99% confidence bounds). Simulated hip abductor strengthening resulted in lower peak JCFs relative to pre-operative for all five patients at the hip (18.9–23.8 ± 16.5%) and knee (20.5–23.8 ± 11.2%). Four of the five patients had reductions at the ankle (7.1–8.5 ± 11.3%) and low back (3.5–7.0 ± 5.3%) with one patient demonstrating no change. The reduction in JCF at the hip joint and at joints other than the hip with hip abductor strengthening demonstrates the dynamic and mechanical interdependencies of the knee, hip and spine that can be targeted in early THA rehabilitation to improve overall patient function.     

Alterations of musculoskeletal models for a more accurate estimation of lower limb joint contact forces during normal gait: A systematic review

Musculoskeletal modelling is a methodology used to investigate joint contact forces during a movement. High accuracy in the estimation of the hip or knee joint contact forces can be obtained with subject-specific models. However, construction of subject-specific models remains time consuming and expensive. The purpose of this systematic review of the literature was to identify what alterations can be made on generic (i.e. literature-based, without any subject-specific measurement other than body size and weight) musculoskeletal models to obtain a better estimation of the joint contact forces. The impact of these alterations on the accuracy of the estimated joint contact forces were appraised.The systematic search yielded to 141 articles and 24 papers were included in the review. Different strategies of alterations were found: skeletal and joint model (e.g. number of degrees of freedom, knee alignment), muscle model (e.g. Hill-type muscle parameters, level of muscular redundancy), and optimisation problem (e.g. objective function, design variables, constraints). All these alterations had an impact on joint contact force accuracy, so demonstrating the potential for improving the model predictions without necessarily involving costly and time consuming medical images. However, due to discrepancies in the reported evidence about this impact and despite a high quality of the reviewed studies, it was not possible to highlight any trend defining which alteration had the largest impact.     

In vivo hip joint loads and pedal forces during ergometer cycling

The rising prevalence of osteoarthritis and an increase in total hip replacements calls for attention to potential therapeutic activities. Cycling is considered as a low impact exercise for the hip joint and hence recommended. However, there are limited data about hip joint loading to support this claim. The aim of this study was to measure synchronously the in vivo hip joint loads and pedal forces during cycling. The in vivo hip joint loads were measured in 5 patients with instrumented hip implants. Data were collected at several combinations of power and cadence, at two saddle heights.Joint loads and pedal forces showed strong linear correlation with power. So the relationship between the external pedal forces and internal joint forces was shown. While cycling at different cadences the minimum joint loads were acquired at 60 RPM. The lower saddle height configuration results in an approximately 15% increase compared to normal saddle height.The results offered new insights into the actual effects of cycling on the hip joint and can serve as useful tools while developing an optimum cycling regimen for individuals with coxarthrosis or following total hip arthroplasty. Due to the relatively low contact forces, cycling at a moderate power level of 90 W at a normal saddle height is suitable for patients.     

Influence of the 3D inverse dynamic method on the joint forces and moments during gait

The joint forces and moments are commonly used in gait analysis. They can be computed by four different 3D inverse dynamic methods proposed in the literature, either based on vectors and Euler angles, wrenches and quaternions, homogeneous matrices, or generalized coordinates and forces. In order to analyze the influence of the inverse dynamic method, the joint forces and moments were computed during gait on nine healthy subjects. A ratio was computed between the relative dispersions (due to the method) and the absolute amplitudes of the gait curves. The influence of the inverse dynamic method was negligible at the ankle (2%) but major at the knee and the hip joints (40%). This influence seems to be due to the dynamic computation rather than the kinematic computation. Compared to the influence of the joint center location, the body segment inertial parameter estimation, and more, the influence of the inverse dynamic method is at least of equivalent importance. This point should be confirmed with other subjects, possibly pathologic, and other movements.     

Joint contact loading in forefoot and rearfoot strike patterns during running

Research concerning forefoot strike pattern (FFS) versus rearfoot strike pattern (RFS) running has focused on the ground reaction force even though internal joint contact forces are a more direct measure of the loads responsible for injury. The main purpose of this study was to determine the internal loading of the joints for each strike pattern. A secondary purpose was to determine if converted FFS and RFS runners can adequately represent habitual runners with regards to the internal joint loading. Using inverse dynamics to calculate the net joint moments and reaction forces and optimization techniques to estimate muscle forces, we determined the axial compressive loading at the ankle, knee, and hip. Subjects consisted of 15 habitual FFS and 15 habitual RFS competitive runners. Each subject ran at a preferred running velocity with their habitual strike pattern and then converted to the opposite strike pattern. Plantar flexor muscle forces and net ankle joint moments were greater in the FFS running compared to the RFS running during the first half of the stance phase. The average contact forces during this period increased by 41.7% at the ankle and 14.4% at the knee joint during FFS running. Peak ankle joint contact force was 1.5 body weights greater during FFS running (p<0.05). There was no evidence to support a difference between habitual and converted running for joint contact forces. The increased loading at the ankle joint for FFS is an area of concern for individuals considering altering their foot strike pattern.     

Comparison of global and joint-to-joint methods for estimating the hip joint load and the muscle forces during walking

A three-dimensional musculoskeletal model of the lower limb was developed to study the influence of biarticular muscles on the muscle force distribution and joint loads during walking. A complete walking cycle was recorded for 9 healthy subjects using the standard optoelectronic motion tracking system. Ground contact forces were also measured using a 6-axes force plate. Inverse dynamics was used to compute net joint reactions (forces and torques) in the lower limb. A static optimization method was then used to estimate muscle forces. Two different approaches were used: in the first one named global method, the biarticular muscles exerted a torque on the two joints they spanned at the same time, and in the second one called joint-by-joint method, these biarticular muscles were divided into two mono-articular muscles with geometrical (insertion, origin, via points) and physiological properties remained unchanged. The hip joint load during the gait cycle was then calculated taking into account the effect of muscle contractions. The two approaches resulted in different muscle force repartition: the biarticular muscles were favoured over any set of single-joint muscles with the same physiological function when using the global method. While the two approaches yielded only little difference in the resultant hip load, the examination of muscle power showed that biarticular muscles could produce positive work at one joint and negative work at the other, transferring energy between body segments and thus decreasing the metabolic cost of movement.     

Effects of hip joint centre mislocation on gait analysis results

Methods to determine the hip joint centre (HJC) location are necessary in gait analysis. It has been demonstrated that the methods proposed in the literature involve large mislocation errors. The choice should be made according to the extent by which HJC location errors distort the estimates of angles and resultant moments at the hip and knee joints. This study aimed at quantifying how mislocation errors propagate to these gait analysis results. Angles and moments at the hip and knee joint were calculated for five able-bodied subjects during level walking. The nominal position of the HJC was determined as the position of the pivot point of a 3D movement of the thigh relative to the pelvis. Angles and moments were then re-calculated after having added to HJC co-ordinates errors in the range of +/-30 mm. Angles and moments at both hip and knee joints were affected by HJC mislocation. The hip moments showed the largest propagation error: a 30 mm HJC anterior mislocation resulted in a propagated error into flexion/extension component of about -22%. The hip abduction/adduction moment was found the second largest affected quantity: a 30 mm lateral HJC mislocation produced a propagated error of about -15%. Finally, a 30 mm posterior HJC mislocation produced a delay of the flexion-to-extension timing in the order of 25% of the stride duration. HJC estimation methods with minimum antero-posterior error should therefore be preferred.     

Simultaneous prediction of muscle and contact forces in the knee during gait

Musculoskeletal models are currently the primary means for estimating in vivo muscle and contact forces in the knee during gait. These models typically couple a dynamic skeletal model with individual muscle models but rarely include articular contact models due to their high computational cost. This study evaluates a novel method for predicting muscle and contact forces simultaneously in the knee during gait. The method utilizes a 12 degree-of-freedom knee model (femur, tibia, and patella) combining muscle, articular contact, and dynamic skeletal models. Eight static optimization problems were formulated using two cost functions (one based on muscle activations and one based on contact forces) and four constraints sets (each composed of different combinations of inverse dynamic loads). The estimated muscle and contact forces were evaluated using in vivo tibial contact force data collected from a patient with a force-measuring knee implant. When the eight optimization problems were solved with added constraints to match the in vivo contact force measurements, root-mean-square errors in predicted contact forces were less than 10 N. Furthermore, muscle and patellar contact forces predicted by the two cost functions became more similar as more inverse dynamic loads were used as constraints. When the contact force constraints were removed, estimated medial contact forces were similar and lateral contact forces lower in magnitude compared to measured contact forces, with estimated muscle forces being sensitive and estimated patellar contact forces relatively insensitive to the choice of cost function and constraint set. These results suggest that optimization problem formulation coupled with knee model complexity can significantly affect predicted muscle and contact forces in the knee during gait. Further research using a complete lower limb model is needed to assess the importance of this finding to the muscle and contact force estimation process.