Refining muscle geometry and wrapping in the TLEM 2 model for improved hip contact force prediction

Sensitivity of knee cartilage biomechanics in finite element analysis to selected Musculoskeletal models

Knee joint kinematics and kinetics analyzed by musculoskeletal (MS) modeling are often utilized in finite element (FE) models, estimating tissue-level mechanical responses. We compared knee cartilage stresses, strains, and centers of pressure of FE models driven by two widely used MS models, implemented in AnyBody and OpenSim. Minor discrepancies in the results were observed between the models. AnyBody-driven FE models showed slightly higher stresses in the medial tibial cartilage, while OpenSim-driven FE models estimated more anterior and lateral center of pressure. Recognizing these differences in the MS-FE models is important to ensure reliable analysis of cartilage mechanics and failure and simulation of rehabilitation.

Cost function criteria using muscle synergies: Exploring the potential of muscle synergy hypothesis

BACKGROUND AND OBJECTIVE

Solving the redundant optimization problem for human muscles depends on the cost function. Choosing the appropriate cost function helps to address a specific problem. Muscle synergies are currently limited to those obtained by electromyography. Furthermore, debate continues regarding whether muscle synergy is derived or real. This study proposes new cost functions based on the muscle synergy hypothesis for solving the optimal muscle force output problem through musculoskeletal modeling.

METHODS

We propose two new computational cost functions involving muscle synergies, which are extracted from muscle activations predicted by musculoskeletal modelling rather than electromyography. In this study, we constructed a musculoskeletal model for simulation using the "Grand Challenge Competition to Predict In Vivo Knee Loads" dataset. Muscle synergies were obtained using non-negative matrix factorization. Two cost functions with muscle synergies were constructed by integrating the polynomial and min/max criterion. Two new functions were verified and validated in normal, smooth, and bouncy gaits.

RESULTS

The muscle synergies based on normal gaits were classified into four modules. The cosine similarities of the first three modules were all >0.9. In the normal and smooth gaits, the forces in most muscles predicted using the two new functions were within three standard deviations of the root mean square error for electromyographic comparisons. Predicted muscle force curves using the four methods as well as characteristic points (i.e., time points in the gait cycle when the significant difference was observed between normal and bouncy gaits) were obtained to validate their predictive capabilities.

CONCLUSIONS

This study constructed two new cost functions involving muscle synergies, verified and validated the ability, and explored the potential of muscle synergy hypothesis.

Virtual reality skateboarding training for balance and functional performance in degenerative lumbar spine disease

BACKGROUND

Degenerative lumbar spine disease (DLD) is a prevalent condition in middle-aged and elderly individuals. DLD frequently results in pain, muscle weakness, and motor impairment, which affect postural stability and functional performance in daily activities. Simulated skateboarding training could enable patients with DLD to engage in exercise with less pain and focus on single-leg weight-bearing. The purpose of this study was to investigate the effects of virtual reality (VR) skateboarding training on balance and functional performance in patients with DLD.

METHODS

Fourteen patients with DLD and 21 age-matched healthy individuals completed a 6-week program of VR skateboarding training. The motion capture and force platform systems were synchronized to collect data during a single-leg stance test (SLST). Musculoskeletal simulation was utilized to calculate muscle force based on the data. Four functional performance tests were conducted to evaluate the improvement after the training. A Visual Analogue Scale (VAS) was also employed for pain assessment.

RESULTS

After the training, pain intensity significantly decreased in patients with DLD (p = 0.024). Before the training, patients with DLD took longer than healthy individuals on the five times sit-to-stand test (p = 0.024). After the training, no significant between-group differences were observed in any of the functional performance tests (p > 0.05). In balance, patients with DLD were similar to healthy individuals after the training, except that the mean frequency (p = 0.014) was higher. Patients with DLD initially had higher biceps femoris force demands (p = 0.028) but shifted to increased gluteus maximus demand after the training (p = 0.037). Gluteus medius strength significantly improved in patients with DLD (p = 0.039), while healthy individuals showed consistent muscle force (p > 0.05).

CONCLUSION

This is the first study to apply the novel VR skateboarding training to patients with DLD. VR skateboarding training enabled patients with DLD to achieve the training effects in a posture that relieves lumbar spine pressure. The results also emphasized the significant benefits to patients with DLD, such as reduced pain, enhanced balance, and improved muscle performance.

Hip joint contact forces are lower in people with femoroacetabular impingement syndrome during squat tasks

It remains unknown if hip joint forces during squat tasks are altered in people with femoroacetabular impingement syndrome (FAIS). The aim of this study is to compare hip joint forces between people with FAIS and healthy controls during double leg squat and single leg squat tasks and within limbs during a single leg squat task in people with FAIS. Kinematic and kinetic data were collected in eight people with FAIS and eight healthy matched controls using 3D motion capture and force plates. AnyBody Modeling System was used to perform musculoskeletal simulations to estimate hip joint angles, forces, and moments for all participants. Estimates were postprocessed with AnyPyTools and converted into normalized time series to be compared using a 1D statistical nonparametric mapping (SnPM) approach. SnPM with an independent samples t-test model was used to compare people with FAIS to controls, while a paired samples model was used to compare involved to uninvolved limb in people with FAIS. Patients demonstrated lower proximodistal force compared to controls (p < 0.01) and compared to the uninvolved side (p = 0.01) for single leg squat. The smaller joint contact forces in people with FAIS compared to controls could represent a strategy of reduced muscle forces to avoid pain and symptoms during this high demand task. These findings when combined with imaging data could help assess the severity of FAIS on hip related function during higher demand tasks.

Investigation of stair ascending and descending activities on the lifespan of hip implants

Total hip arthroplasty (THA) surgeries among young patients are on the increase, so it is crucial to predict the lifespan of hip implants correctly and produce solutions to improve longevity. Current implants are designed and tested against walking conditions to predict the wear rates. However, it would be reasonable to include the additional effects of other daily life activities on wear rates to predict convergent results to clinical outputs. In this study, 14 participants are recruited to perform stair ascending (AS), descending (DS), and walking activities to obtain kinematic and kinetic data for each cycle using marker based Qualisys motion capture (MOCAP) system. AnyBody Modeling System using the Calibrated Anatomical System Technique (CAST) full body marker set are performed Multibody simulations. The 3D generic musculoskeletal model used in this study is a marker-based full-body motion capture model (AMMR,2.3.1 MoCapModel) consisting of the upper extremity and the Twente Lower Extremity Model (TLEM2). The dynamic wear prediction model detailing the intermittent and overall wear rates for CoCr-on-XLPE bearing couple is developed to investigate the wear mechanism under 3D loading for AS, DS, and walking activities over 5 million cycles (Mc) by using finite element modelling technique. The volumetric wear rates of XLPE liner under AS, DS, and walking activities over 5-Mc are predicted as 27.43, 23.22, and 18.84 mm3/Mc respectively. Additionally, the wear rate was predicted by combining stair activities and gait cycles based on the walk-to-stair ratio. By adding the effect of stair activities, the volumetric wear rate of XLPE is predicted as 22.02 mm3/Mc which is equivalent to 19.41% of walking. In conclusion, in this study, the effect of including other daily life activities is demonstrated and evidence is provided by matching them to the clinical data as opposed to simulator test results of implants under ISO 14242 boundary conditions.

Load-induced blood marker kinetics in patients with medial knee compartment osteoarthritis are associated with accumulated load and patient reported outcome measures

Background

This study aimed to quantify the mechanoresponse of 10 blood marker candidates for joint metabolism to a walking stress test in patients with knee osteoarthritis and to determine the association among marker kinetics and with accumulated load and patient reported outcomes.

Methods

24 patients with knee osteoarthritis completed questionnaires, and a 30-minute walking stress test with six blood serum samples and gait analysis. Concentrations of cartilage oligomeric matrix protein (COMP), matrix metalloproteinases (MMP)-1, -3, and -9, epitope resulting from cleavage of type II collagen by collagenases (C2C), type II procollagen (CPII), interleukin (IL)-6, proteoglycan (PRG)-4, A disintegrin and metalloproteinase with thrombospondin motifs (ADAMTS)-4, and resistin were determined by enzyme-linked immunosorbent assays, Joint load (moments and compartmental forces) was estimated using musculoskeletal modeling using gait analysis data.

Results

COMP and MMP-3 showed an immediate increase after the walking stress followed by a decrease. MMP-9 and resistin showed a delayed decrease below pre-stress levels. ∆COMP correlated with ∆MMP-3 for most time points. ∆MMP-9 correlated with ∆resistin for most time points. The load-induced increase in blood marker levels correlated among blood markers and time points. C2C and resistin correlated positively and C2C/CPII and MMP2 correlated negatively with load during gait. Immediate relative ∆CPII and ∆MMP1 and delayed relative ∆COMP, ∆IL6, ∆C2C, ∆CPII, ∆MMP1 and ∆MMP3 correlated with the load accumulated during the walking stress. Baseline C2C levels correlated with Knee Osteoarthritis Outcome Score (KOOS) subscales and load-induced changes in MMP-3 with KOOS and Short Form 36 quality of life subscores (P<0.05).

Conclusions

The distinct and differentiated physiological response to the walking stress depends on accumulated load and appears relevant for patient reported osteoarthritis outcome and quality of life and warrants further investigation in the context of disease progression.ClinicalTrials.gov registration: NCT02622204.

EMG-Driven Musculoskeletal Model Calibration With Wrapping Surface Personalization

Muscle forces and joint moments estimated by electromyography (EMG)-driven musculoskeletal models are sensitive to the wrapping surface geometry defining muscle-tendon lengths and moment arms. Despite this sensitivity, wrapping surface properties are typically not personalized to subject movement data. This study developed a novel method for personalizing OpenSim cylindrical wrapping surfaces during EMG-driven model calibration. To avoid the high computational cost of repeated OpenSim muscle analyses, the method uses two-level polynomial surrogate models. Outer-level models specify time-varying muscle-tendon lengths and moment arms as functions of joint angles, while inner-level models specify time-invariant outer-level polynomial coefficients as functions of wrapping surface parameters. To evaluate the method, we used walking data collected from two individuals post-stroke and performed four variations of EMG-driven lower extremity model calibration: 1) no calibration of scaled generic wrapping surfaces (NGA), 2) calibration of outer-level polynomial coefficients for all muscles (SGA), 3) calibration of outer-level polynomial coefficients only for muscles with wrapping surfaces (LSGA), and 4) calibration of cylindrical wrapping surface parameters for muscles with wrapping surfaces (PGA). On average compared to NGA, SGA reduced lower extremity joint moment matching errors by 31%, LSGA by 24%, and PGA by 12%, with the largest reductions occurring at the hip. Furthermore, PGA reduced peak hip joint contact force by 47% bodyweight, which was the most consistent with published in vivo measurements. The proposed method for EMG-driven model calibration with wrapping surface personalization produces physically realistic OpenSim models that reduce joint moment matching errors while improving prediction of hip joint contact force.

The effect of pathological shoulder rhythm on muscle and joint forces after reverse shoulder arthroplasty, a numerical analysis

BACKGROUND

Compromised abduction ability after reverse shoulder arthroplasty is primarily linked to limited glenohumeral range of motion while scapulothoracic mobility can typically be maintained. Glenohumeral joint forces strongly depend on the resulting scapulohumeral rhythm, however, an association between the acting muscle and joint forces and the subject-specific scapulohumeral rhythm after reverse shoulder arthroplasty has not been established.

METHODS

Eleven reverse shoulder arthroplasty patients were divided into groups of poor and excellent abduction ability. Subject-specific models were developed and scaled for each patient using existing motion capture data in AnyBody™. Shoulder muscle and joint forces were obtained using inverse dynamics calculations during shoulder abduction to 100° in the scapula plane. The scapulohumeral rhythm, the resting abduction angle and internal body forces between the outcome groups were compared using a Mann Whitney U test.

FINDINGS

The mean glenohumeral and scapulothoracic contribution to overall shoulder abduction for the excellent group was on average 9.7% higher and 21.4% lower, respectively, compared to the mean of the poor group. For shoulder abduction angles between 30° and 60°, the excellent group demonstrated on average 25% higher muscle forces in the anterior deltoid which was significantly higher compared to the poor outcome patients. Scapulothoracic muscle activity did not differ significantly between the two functional groups.

INTERPRETATION

Accordingly, rehabilitation strategies focusing on strengthening the anterior part of the deltoid in particular may improve clinical outcomes.

Reverse dynamics analysis of contact force and muscle activities during the golf swing after total hip arthroplasty

There are no reports on hip kinetics including contact forces and muscle activities during the golf swing after total hip arthroplasty (THA). The aim of this study was to identify the characteristics of three-dimensional dynamics during the golf swing. Ten unilateral primary THA patients participated in motion capture test of their driver golf swing. The driver swing produced approximately 20-30° of rotation in both lead and trail replaced hips. The mean hip contact forces (HCFs) of lead and trail replaced hips were 5.1 and 6.6 × body weight, respectively. Left and right THAs showed similar HCFs of lead and trail hips. More than 60% of the Percent maximum voluntary isometric contraction was found in bilateral iliopsoas muscles in all unilateral THA. Three factors [female sex, lower modified Harris Hip Score, and higher HCF of surgical side] were associated with the golf-related replacement hip pain. Golf is an admissible sport after THA because driver swings do not contribute excessive rotation or contact forces to hip prostheses. HCF could be reduced through swing adjustments, which may allow patients with golf-related replacement hip pain to develop a comfortable golf game free from pain.

The Non-Affected Muscle Volume Compensates for the Partial Loss of Strength after Injection of Botulinum Toxin A

Local botulinum toxin (BTX-A, Botox®) injection in overactive muscles is a standard treatment in patients with cerebral palsy. The effect is markedly reduced in children above the age of 6 to 7. One possible reason for this is the muscle volume affected by the drug. Nine patients (aged 11.5; 8.7–14.5 years) with cerebral palsy GMFCS I were treated with BTX-A for equinus gait at the gastrocnemii and soleus muscles. BTX-A was administered at one or two injection sites per muscle belly and with a maximum of 50 U per injection site. Physical examination, instrumented gait analysis, and musculoskeletal modelling were used to assess standard muscle parameters, kinematics, and kinetics during gait. Magnetic resonance imaging (MRI) was used to detect the affected muscle volume. All the measurements were carried out pre-, 6 weeks post-, and 12 weeks post-BTX-A. Between 9 and 15% of the muscle volume was affected by BTX-A. There was no effect on gait kinematics and kinetics after BTX-A injection, indicating that the overall kinetic demand placed on the plantar flexor muscles remained unchanged. BTX-A is an effective drug for inducing muscle weakness. However, in our patient cohort, the volume of the affected muscle section was limited, and the remaining non-affected parts were able to compensate for the weakened part of the muscle by taking over the kinetic demands associated with gait, thus not enabling a net functional effect in older children. We recommend distributing the drug over the whole muscle belly through multiple injection sites.

Prediction of ground reaction forces and moments during walking in children with cerebral palsy

Introduction

Gait analysis is increasingly used to support clinical decision-making regarding diagnosis and treatment planning for movement disorders. As a key part of gait analysis, inverse dynamics can be applied to estimate internal loading conditions during movement, which is essential for understanding pathological gait patterns. The inverse dynamics calculation uses external kinetic information, normally collected using force plates. However, collection of external ground reaction forces (GRFs) and moments (GRMs) can be challenging, especially in subjects with movement disorders. In recent years, a musculoskeletal modeling-based approach has been developed to predict external kinetics from kinematic data, but its performance has not yet been evaluated for altered locomotor patterns such as toe-walking. Therefore, the goal of this study was to investigate how well this prediction method performs for gait in children with cerebral palsy.

Methods

The method was applied to 25 subjects with various forms of hemiplegic spastic locomotor patterns. Predicted GRFs and GRMs, in addition to associated joint kinetics derived using inverse dynamics, were statistically compared against those based on force plate measurements.

Results

The results showed that the performance of the predictive method was similar for the affected and unaffected limbs, with Pearson correlation coefficients between predicted and measured GRFs of 0.71–0.96, similar to those previously reported for healthy adults, despite the motor pathology and the inclusion of toes-walkers within our cohort. However, errors were amplified when calculating the resulting joint moments to an extent that could influence clinical interpretation.

Conclusion

To conclude, the musculoskeletal modeling-based approach for estimating external kinetics is promising for pathological gait, offering the possibility of estimating GRFs and GRMs without the need for force plate data. However, further development is needed before implementation within clinical settings becomes possible.

Subject-specific computational modeling of acromioclavicular and coracoclavicular ligaments

BACKGROUND

Disruption of the acromioclavicular joint is a common injury. Despite the different surgical procedures described for treating this injury, complications such as loss of reduction and failure of fixation remain unacceptably high. We developed a computer model of the acromioclavicular joint to better understand the biomechanical contributions of the ligaments that are typically injured.

METHODS

Six fresh frozen human cadaveric shoulders were tested on an AMTI VIVO 6-degree of freedom test platform to measure force-displacement in inferior translation, anteroposterior translation, and internal rotation before and after sequentially transecting the coracoclavicular and acromioclavicular ligaments. These data were used to construct computer models of each specimen. Three-dimensional computed tomographic scans were used to generate a rigid-body dynamics model using the AnyBody Modeling System. The scapula and clavicle were connected by the acromioclavicular joint capsule, the conoid ligament, and the trapezoid ligament. Subject-specific ligament properties were calculated by matching computer predictions to experimental force-displacement data.

RESULTS

The calculated free lengths of the conoid, trapezoid, and acromioclavicular ligaments were 13.5 (±3.2), 11.8 (±2.4), and 11.0 (±2.7) mm, respectively. The calculated stiffnesses of the conoid, trapezoid, and acromioclavicular ligaments were 34.3 (±6.3), 28.4 (±3.2), and 33.8 (±8.2) N/mm, respectively. Root mean square deviation (RMSD) of predicted force-displacement curves relative to experimental force-displacement curves (during inferior and anteroposterior translation) was less than 1 mm. For validation of subject-specific models, after ligament properties were calculated, the RMSD of the predicted torque over 15° of internal rotation was 12% of maximum rotational torque (average for 6 specimens).

DISCUSSION AND CONCLUSION

Acromioclavicular disruption results in multidirectional instability, which requires careful consideration of the individual contributions of the injured ligaments. In addition, variations in patient anatomy can significantly affect the biomechanical stability of the reconstruction. Subject-specific models can enhance our understanding of the individual and collective biomechanical contributions of the injured soft tissues to the multiaxial stability of the acromioclavicular joint. These models may also be useful for analyzing and assessing biomechanical stability after various types of surgical reconstruction.

High tibial osteotomy effectively redistributes compressive knee loads during walking

The objectives of this study were to estimate pre- and postoperative lower limb kinematics and kinetics and knee intra-articular forces during gait using musculoskeletal modeling in a cohort of patients with knee osteoarthritis (OA) undergoing high tibial osteotomy (HTO), compare these to controls, and determine correlations between changes in these parameters and Knee Injury and Osteoarthritis Outcome Score (KOOS) subscores after HTO. Sixteen patients with isolated, symptomatic medial compartment knee OA completed pre- and postoperative gait analysis (mean follow-up time: 8.6 months). Sixteen age- and sex-matched asymptomatic volunteers participated as controls. Musculoskeletal modeling was used to evaluate lower limb joint moments and knee contact forces during gait. While HTO had limited influence on sagittal plane kinematics and moments, significant changes in the load distribution at the knee after HTO were observed with a lower postoperative compressive load on the medial compartment during midstance and a higher compressive load on the lateral compartment during early and late stance. Moreover, the lateral shear force in midstance was significantly lower after HTO. Changes in the external knee adduction moment (KAM) did not always coincide with reductions in the knee compressive force in the medial compartment. Biomechanical changes did not correlate with improvements in KOOS subscores. Hence, HTO effectively unloaded the medial compartment by redistributing part of the overall compressive force to the lateral compartment during gait with limited influence on gait function. The KAM may not adequately describe compartmental load magnitude or changes induced by interventions at the compartment level. Clinical trial registration: ClinicalTrials. gov Identifier-NCT02622204. Clinical significance: This study provides important evidence for changes in joint level loads after corrective osteotomy as joint preserving surgery and emphasizes the need for additional biomechanical outcomes of such interventions.

The functional role of hip muscles during gait in patients with increased femoral anteversion

BACKGROUND

Femoral anteversion affects the lever arm and moment-generating capacity of the hip abductors, while an increased hip internal rotation during walking was proposed to be a compensatory mechanism to restore the abductive lever arm. Children with isolated increased femoral anteversion, however, do not always present a deficit in the net hip abduction moment during gait, suggesting that a more comprehensive understanding of the effect of morphology and motion on muscle forces and moments is needed to aid clinical decision making.

RESEARCH QUESTION

Are muscle contributions to hip joint moments and muscle forces altered in patients with increased femoral anteversion and internally rotated gait pattern compared to a control group of typically developing children? And how would the functional role of the muscle be altered if the patients walked straight?

METHODS

This follow-up study compared patients with increased femoral anteversion (n = 42, 12.8 ± 1.9 years, femoral anteversion: 39.6 ± 6.9°) to controls (n = 9, 12.0 ± 3.0 years, femoral anteversion: 18.7 ± 4.1°). Muscle forces and moment contributions were calculated using personalized musculoskeletal models. Additionally, a hypothetical scenario, in which the gait of the controls was modelled with an anteverted femoral morphology, was used to understand what would happen if the patients walked straight.

RESULTS

Gluteus medius abductive contribution was lower in patients compared to controls, despite a comparable net abduction moment around the hip. Patients presented lower muscle forces. However, if modelled to walk straight, they would require higher forces as well as a larger co-contraction of both hip internal and external rotators in the transversal plane.

SIGNIFICANCE

This study suggests that patients with increased femoral anteversion walking with an internally rotated gait pattern present lower muscle forces, but when modelled to walk straight muscle forces increase. The current results provide important information to better understand this condition and improve treatment recommendations in these patients.

On the estimation of hip joint loads through musculoskeletal modeling

Noninvasive estimation of joint loads is still an open challenge in biomechanics. Although musculoskeletal modeling represents a solid resource, multiple improvements are still necessary to obtain accurate predictions of joint loads and to translate such potential into practical utility. The present study, focused on the hip joint, is aimed at reviewing the state-of-the-art literature on the estimation of hip joint reaction forces through musculoskeletal modeling. Our literature inspection, based on well-defined selection criteria, returned seventeen works, which were compared in terms of methods and results. Deviations between predicted and in vivo measured hip joint loads, taken from the OrthoLoad database, were assessed through quantitative deviation indices. Despite the numerous modeling and computational improvements made over the last two decades, predicted hip joint loads still deviate from their experimental counterparts and typically overestimate them. Several critical aspects have emerged that affect muscle force estimation, hence joint loads. Among them, the physical fidelity of the musculoskeletal model, with its parameters and geometry, plays a crucial role. Also, predicted joint loads are markedly affected by the selected muscle recruitment strategy, which reflects the underlying motor control policy. Practical guidelines for researchers interested in noninvasive estimation of hip joint loads are also provided.

Effect of seat tilting on biomechanics of L4 in tractor drivers

The tilting of the tractors' seat during deep tillage operations affects the drivers' sitting position, which can lead to lumbar spine injuries. To investigate the effect of seat tilting on the driver's L4 biomechanics, we built a tractor-driver musculoskeletal model in AnyBody™. The maximum activity of the driver's erector spinae at different tilted angles were measured by sEMG and compared with the simulation results to validate the model. The spatial position of the driver's spine at different tilted angles were obtained by 3 D motion capture. The model simulated the driver's spine posture during the actual tilt and investigated the effects of different tilted angles and vibration on the biomechanics of the driver's L4 . The results showed that as the tilt angle of the tractor increased, the load on the driver's L4 also increased, especially the shear force increased at a faster rate than the axial and normal forces, with the shear force on the driver's L4 increasing from 0 N to 138.7 N when the tractor was tilted from 0° to 15°. When vibration was applied to the musculoskeletal model, the maximum value of the shear force on the driver became progressively greater as the angle of tilt of the tractor increased. Overall, tilting the tractor can have a large impact on the biomechanics of the driver's L4, and tilting the tractor may be an important cause of lumbar spine injuries for tractor drivers.

Muscle coordination retraining inspired by musculoskeletal simulations reduces knee contact force

Humans typically coordinate their muscles to meet movement objectives like minimizing energy expenditure. In the presence of pathology, new objectives gain importance, like reducing loading in an osteoarthritic joint, but people often do not change their muscle coordination patterns to meet these new objectives. Here we use musculoskeletal simulations to identify simple changes in coordination that can be taught using electromyographic biofeedback, achieving the therapeutic goal of reducing joint loading. Our simulations predicted that changing the relative activation of two redundant ankle plantarflexor muscles—the gastrocnemius and soleus—could reduce knee contact force during walking, but it was unclear whether humans could re-coordinate redundant muscles during a complex task like walking. Our experiments showed that after a single session of walking with biofeedback of summary measures of plantarflexor muscle activation, healthy individuals reduced the ratio of gastrocnemius-to-soleus muscle activation by 25 ± 15% (p = 0.004, paired t test, n = 10). Participants who walked with this “gastrocnemius avoidance” gait pattern reduced late-stance knee contact force by 12 ± 12% (p = 0.029, paired t test, n = 8). Simulation-informed coordination retraining could be a promising treatment for knee osteoarthritis and a powerful tool for optimizing coordination for a variety of rehabilitation and performance applications.

Comparison of kinematic parameters of children gait obtained by inverse and direct models

The purpose of this study is to compare differences between kinematic parameters of pediatric gait obtained by direct kinematics (DK) (Plug-in-Gait) and inverse kinematics (IK) (AnyBody) models. Seventeen healthy children participated in this study. Both lower extremities were examined using a Vicon 8-camera motion capture system and a force plate. Angles of the hip, knee, and ankle joints were obtained based on DK and IK models, and ranges of motion (ROMs) were identified from them. The standard error of measurement, root-mean-squared error, correlation r, and magnitude-phase (MP) metrics were calculated to compare differences between the models’ outcomes. The determined standard error of measurement between ROMs from the DK and IK models ranged from 0.34° to 0.58°. A significant difference was found in the ROMs with the exception of the left hip’s internal/external rotation. The mean RMSE of all joints’ amplitudes exceeded the clinical significance limit and was 13.6 ± 4.0°. The best curve angles matching nature were found in the sagittal plane, where r was 0.79 to 0.83 and MP metrics were 0.05 to 0.30. The kinematic parameters of pediatric gait obtained by IK and DK differ significantly. Preferably, all of the results obtained by DK must be validated/verified by IK, in order to achieve a more accurate functional assessment of the individual. Furthermore, the use of IK expands the capabilities of gait analysis and allows for kinetic characterisation.

Increased Femoral Anteversion Does Not Lead to Increased Joint Forces During Gait in a Cohort of Adolescent Patients

Orthopedic complications were previously reported for patients with increased femoral anteversion. A more comprehensive analysis of the influence of increased femoral anteversion on joint loading in these patients is required to better understand the pathology and its clinical management. Therefore, the aim was to investigate lower-limb kinematics, joint moments and forces during gait in adolescent patients with increased, isolated femoral anteversion compared to typically developing controls. Secondly, relationships between the joint loads experienced by the patients and different morphological and kinematic features were investigated. Patients with increased femoral anteversion (n = 42, 12.8 ± 1.9 years, femoral anteversion: 39.6 ± 6.9°) were compared to typically developing controls (n = 9, 12.0 ± 3.0 years, femoral anteversion: 18.7 ± 4.1°). Hip and knee joint kinematics and kinetics were calculated using subject-specific musculoskeletal models. Differences between patients and controls in the investigated outcome variables (joint kinematics, moments, and forces) were evaluated through statistical parametric mapping with Hotelling T2 and t-tests (α = 0.05). Canonical correlation analyses (CCAs) and regression analyses were used to evaluate within the patients’ cohort the effect of different morphological and kinematic predictors on the outcome variables. Predicted compressive proximo-distal loads in both hip and knee joints were significantly reduced in patients compared to controls. A gait pattern characterized by increased knee flexion during terminal stance (KneeFlex_tSt) was significantly correlated with hip and knee forces, as well as with the resultant force exerted by the quadriceps on the patella. On the other hand, hip internal rotation and in-toeing, did not affect the loads in the joints. Based on the finding of the CCAs and linear regression analyses, patients were further divided into two subgroups based KneeFlex_tSt. Patients with excessive KneeFlex_tSt presented a significantly higher femoral anteversion than those with normal KneeFlex_tSt. Patients with excessive KneeFlex_tSt presented significantly larger quadriceps forces on the patella and a larger posteriorly-oriented shear force at the knee, compared to patients with normal KneeFlex_tSt, but both patients’ subgroups presented only limited differences in terms of joint loading compared to controls. This study showed that an altered femoral morphology does not necessarily lead to an increased risk of joint overloading, but instead patient-specific kinematics should be considered.

Multiscale modelling for investigating the long-term time-dependent biphasic behaviour of the articular cartilage in the natural hip joint

A better understanding of the time-dependent biomechanical behaviour of the biphasic hip articular cartilage (AC) under physiological loadings is important to understand the onset of joint pathology and guide the clinical treatment. Current computational studies for the biphasic hip AC were usually limited to short-term duration or using elaborate loading. The present study aimed to develop a multiscale computational modelling to investigate the long-term biphasic behaviour of the hip AC under physiological loadings over multiple gait cycles. Two-scale computational modelling including a musculoskeletal model and a finite element model of the natural hip was created. These two models were then combined and used to investigate the biphasic behaviour of hip AC over 80 gait cycles. The results showed that the interstitial fluid pressure in the AC supported over 89% of the loading during gait. When the contact area was located at the AC centre, the contact pressure and fluid pressure increased over time from the first cycle to the 80th cycle, while when the contact area approached the edge, these pressures decreased first dramatically and then slowly over time. The peak stresses and strains in the solid matrix of the AC remained at a low level and increased over time from the first cycle to the 80th cycle. This study demonstrated that the long-term temporal variations of the biphasic behaviour of hip AC under physiological loadings are significant. The methodology has potentially important implications in the biomechanical studies of human cartilage and supporting the development of cartilage substitution.

Altered Muscle Contributions are Required to Support the Stance Limb During Voluntary Toe-Walking

Toe-walking characterizes several neuromuscular conditions and is associated with a reduction in gait stability and efficiency, as well as in life quality. The optimal choice of treatment depends on a correct understanding of the underlying pathology and on the individual biomechanics of walking. The objective of this study was to describe gait deviations occurring in a cohort of healthy adult subjects when mimicking a unilateral toe-walking pattern compared to their normal heel-to-toe gait pattern. The focus was to characterize the functional adaptations of the major lower-limb muscles which are required in order to toe walk. Musculoskeletal modeling was used to estimate the required muscle contributions to the joint sagittal moments. The support moment, defined as the sum of the sagittal extensive moments at the ankle, knee, and hip joints, was used to evaluate the overall muscular effort necessary to maintain stance limb stability and prevent the collapse of the knee. Compared to a normal heel-to-toe gait pattern, toe-walking was characterized by significantly different lower-limb kinematics and kinetics. The altered kinetic demands at each joint translated into different necessary moment contributions from most muscles. In particular, an earlier and prolonged ankle plantarflexion contribution was required from the soleus and gastrocnemius during most of the stance phase. The hip extensors had to provide a higher extensive moment during loading response, while a significantly higher knee extension contribution from the vasti was necessary during mid-stance. Compensatory muscular activations are therefore functionally required at every joint level in order to toe walk. A higher support moment during toe-walking indicates an overall higher muscular effort necessary to maintain stance limb stability and prevent the collapse of the knee. Higher muscular demands during gait may lead to fatigue, pain, and reduced quality of life. Toe-walking is indeed associated with significantly larger muscle forces exerted by the quadriceps to the patella and prolonged force transmission through the Achilles tendon during stance phase. Optimal treatment options should therefore account for muscular demands and potential overloads associated with specific compensatory mechanisms.

Multiscale biomechanics of the biphasic articular cartilage in the natural hip joint during routine activities

BACKGROUND AND OBJECTIVE

The investigation of the biomechanical behaviour of the articular cartilage (AC) under physiological loading is important to understand the joint function and onset of pathologies. This study aimed to develop a multiscale computational modelling approach and apply the approach to investigate the time-dependant biphasic behaviour of the AC in the natural hip joint under repetitive physiological loading over 80 cycles amongst six routine activities.

METHODS

A subject-specific musculoskeletal multibody dynamics (MBD) model was developed based on the anthropometry and motion capture data collected for a male subject. A corresponding FE model of the natural hip joint with biphasic AC was created based on the bone geometries exported from the MBD model. A multiscale computational modelling was then developed to couple the MBD model and the FE model and used to investigate the time-dependant biphasic behaviour of the AC under subject-specific physiological loading over 80 cycles amongst six routine activities.

RESULTS

The results showed that for all the activities considered, the interstitial fluid pressure in the AC supported over 80% of the loading. The maximum values of the peak contact pressure and peak fluid pressure for the whole cycle increased firstly and then remained stable over time from the 1st cycle to the 80th cycle. At these instants, the contact areas were located at the centre region of the AC. By contrast, when the contact area was located at the edge of the AC, these peak pressures were found to increase over time for some of the activities (squat, ascending stairs, descending stairs) but decrease for the other activities (normal walking, standing up, sitting down).

CONCLUSION

This study for the first time developed a multiscale computational modelling approach to couple a musculoskeletal MBD model of the body and a detailed FE model of the natural hip joint with biphasic AC, which enabled the evaluation of time-dependant biphasic behaviour of the AC under realistic physiological loading conditions. The study may have important implications in biomechanical studies of human cartilage to understand the joint function and biomechanical factors related to joint disease, and to support the development of cartilage substitution.

Sensitivity of musculoskeletal models to variation in muscle architecture parameters

Musculoskeletal models, like all theoretical models of physical processes, depend on the assumptions needed to construct the model. For musculoskeletal models, these assumptions include, among other things, the kinematic data, the kinetic data and the muscle parameters. The former (dynamic) data can be acquired relatively easily from living subjects, but the latter are usually based on limited information, frequently determined from cadaver studies performed on elderly individuals. Previously, we determined the sensitivity of forces to dynamic differences among 10 humans walking on a straight path. Here, we assess the sensitivity of the muscle and joint reaction forces developed in human walking to variable muscle parameters obtained from 10 living adults, whose data were recently reported, and compared the results with the values from a standard model that depends on cadaveric data. We found that, while the force patterns across the stance cycle were similar among muscle parameter models, differences of as much as 15% in the force magnitude were produced. Whether or not the variation between the standard model and other muscle parameters is important depends on why the forces are required.

Association between global sagittal malalignment and increasing hip joint contact force, analyzed by a novel musculoskeletal modeling system

Patients with adult spinal deformity have various standing postures. Although several studies have reported a relationship between sagittal alignment and exacerbation of hip osteoarthritis, information is limited regarding how spinopelvic sagittal alignment changes affect hip joint loading. This study aimed to investigate the relationship between sagittal spinopelvic-lower limb alignment and the hip joint contact force (HCF) using a novel musculoskeletal model. We enrolled 20 women (78.3±6.7 years) from a single institution. Standing lateral radiographs were acquired to measure thoracic kyphosis, lumbar lordosis, the pelvic tilt, sacral slope, sagittal vertical axis (SVA), femur obliquity angle, and knee flexion angle. In the model simulation, the Anybody Modeling System was used, which alters muscle pathways using magnetic resonance imaging data. Each patient’s alignment was entered into the model; the HCF and hip moment in the standing posture were calculated using inverse dynamics analysis. The relationship between the HCF and each parameter was examined using Spearman’s correlation coefficient (r). The patients were divided into low SVA and high SVA groups, with a cutoff value of 50 mm for the SVA. The HCF was 168.2±60.1 N (%BW) and positively correlated with the SVA (r = 0.6343, p<0.01) and femur obliquity angle (r = 0.4670, p = 0.03). The HCF were 122.2 and 214.1 N (75.2% difference) in the low SVA and high SVA groups, respectively (p<0.01). The flexion moment was also increased in the high SVA group compared with that in the low SVA group (p = 0.03). The SVA and femur obliquity angle are factors related to the HCF, suggesting an association between adult spinal deformity and the exacerbation of hip osteoarthritis. Future studies will need to assess the relationship between the hip joint load and sagittal spinopelvic parameters in dynamic conditions.

Influence of stance width and toe direction on medial knee contact force during bodyweight squats

Squats are frequently performed to strengthen the quadriceps (Quad) and gluteus maximus (GM) in sports and clinical fields. Since the squat itself produces a large knee contact force, clarifying the relationship between the squat techniques and the knee contact force is important. However, the influence of different squat techniques on the medial knee contact force (KCFmed), which would result in knee disease, remains unclear. This study aimed to investigate the influence of various squat techniques on KCFmed during bodyweight squats. Since muscle strengthening by the squat is inevitable, we additionally aimed to explore the effect of a different squat technique on the quadriceps (Quad) and gluteus maximus (GM) forces. Twelve healthy adults performed squats with different stance widths (narrow stance, NS; middle stance, MS; and wide stance, WS) and different toe directions (0° of forefoot abduction - NEUT and 30°forefoot abduction - OUT). The KCFmed, Quad force, and GM force were computed using a musculoskeletal model with marker trajectories and ground reaction forces. The KCFmed in NS was significantly larger than that in MS and WS, and KCFmed in OUT was significantly larger than that in NEUT. The Quad force in OUT was significantly larger than that in the NEUT, and the GM force significantly became larger as the stance width became narrower. These findings suggest that squats in MS and NEUT may be suitable for reducing KCFmed while maintaining the Quad and GM forces.

A review of musculoskeletal modelling of human locomotion

Locomotion through the environment is important because movement provides access to key resources, including food, shelter and mates. Central to many locomotion-focused questions is the need to understand internal forces, particularly muscle forces and joint reactions. Musculoskeletal modelling, which typically harnesses the power of inverse dynamics, unites experimental data that are collected on living subjects with virtual models of their morphology. The inputs required for producing good musculoskeletal models include body geometry, muscle parameters, motion variables and ground reaction forces. This methodological approach is critically informed by both biological anthropology, with its focus on variation in human form and function, and mechanical engineering, with a focus on the application of Newtonian mechanics to current problems. Here, we demonstrate the application of a musculoskeletal modelling approach to human walking using the data of a single male subject. Furthermore, we discuss the decisions required to build the model, including how to customize the musculoskeletal model, and suggest cautions that both biological anthropologists and engineers who are interested in this topic should consider.

Effects of load mass and position on the dynamic loading of the knees, shoulders and lumbar spine during lifting: a musculoskeletal modelling approach

Musculoskeletal models may enhance our understanding of the dynamic loading of the joints during manual material handling. This study used state-of-the-art musculoskeletal models to determine the effects of load mass, asymmetry angle, horizontal location and deposit height on the dynamic loading of the knees, shoulders and lumbar spine during lifting. Recommended weight limits and lifting indices were also calculated using the NIOSH lifting equation. Based on 1832 lifts from 22 subjects, we found that load mass had the most substantial effect on L5-S1 compression. Increments in asymmetry led to large increases in mediolateral shear, while load mass and asymmetry had significant effects on anteroposterior shear. Increased deposit height led to higher shoulder forces, while the horizontal location mostly affected the forces in the knees and shoulders. These results generally support the findings of previous research, but notable differences in the trends and magnitudes of the estimated forces were observed.

Quantifying Muscle Forces and Joint Loading During Hip Exercises Performed With and Without an Elastic Resistance Band

An increase in hip joint contact forces (HJCFs) is one of the main contributing mechanical causes of hip joint pathologies, such as hip osteoarthritis, and its progression. The strengthening of the surrounding muscles of the joint is a way to increase joint stability, which results in the reduction of HJCF. Most of the exercise recommendations are based on expert opinions instead of evidence-based facts. This study aimed to quantify muscle forces and joint loading during rehabilitative exercises using an elastic resistance band (ERB). Hip exercise movements of 16 healthy volunteers were recorded using a three-dimensional motion capture system and two force plates. All exercises were performed without and with an ERB and two execution velocities. Hip joint kinematics, kinetics, muscle forces, and HJCF were calculated based on the musculoskeletal simulations in OpenSim. Time-normalized waveforms of the different exercise modalities were compared with each other and with reference values found during walking. The results showed that training with an ERB increases both target muscle forces and HJCF. Furthermore, the ERB reduced the hip joint range of motion during the exercises. The type of ERB used (soft vs. stiff ERB) and the execution velocity of the exercise had a minor impact on the peak muscle forces and HJCF. The velocity of exercise execution, however, had an influence on the total required muscle force. Performing hip exercises without an ERB resulted in similar or lower peak HJCF and lower muscle forces than those found during walking. Adding an ERB during hip exercises increased the peak muscle and HJCF but the values remained below those found during walking. Our workflow and findings can be used in conjunction with future studies to support exercise design.

Contact force path in total hip arthroplasty: effect of cup medialisation in a whole-body simulation

BACKGROUND

Cup medialisation down to the true acetabular floor in total hip arthroplasty with a compensatory femoral offset increase seems to be mechanically advantageous for the abductor muscles due to the relocation of the lever arms (body weight lever arm decreased, abductor lever arm increased). However, limited information is currently available about the effects of this reconstruction type at the head cup interface, compared to an anatomical reconstruction that maintains the natural lever arms. Through a whole-body simulation analysis, we compared medialised versus anatomical reconstruction in THA to analyse the effects on: (1) contact force magnitude at the head cup interface; (2) contact force path in the cup; and (3) abductor activity.

METHODS

Musculoskeletal simulations were performed to calculate the above-mentioned parameters using inverse dynamics analysis. The differences between the virtually implanted THAs were calculated to compare the medialised versus anatomical reconstruction.

RESULTS

Cup medialisation with compensatory femoral offset increase led to: (1) a reduction in contact force magnitude at the head cup interface up to 6.6%; (2) a similar contact force path in the cup in terms of sliding distance and aspect ratio; and (3) a reduction in abductor activity up to 17.2% (gluteus medius).

CONCLUSIONS

In our opinion, these potential biomechanical gains do not generally justify a fully medialised reconstruction, especially in younger patients that are more likely to undergo revision surgery in their lifetime. Cup medialisation should be performed until sufficient press fit and bony coverage of a properly sized and oriented cup can be achieved.

Electromyography Based Validation of a Musculoskeletal Hand Model

Regarding the prevention of injuries and rehabilitation of the human hand, musculoskeletal simulations using an inverse dynamics approach allow for insights of the muscle recruitment and thus acting forces on the hand. Currently, several hand models from various research groups are in use, which are mainly validated by the comparison of numerical and anatomical moment arms. In contrast to this validation and model-building technique by cadaver studies, the aim of the present study is to further validate a recently published hand model [1] by analyzing numerically calculated muscle activities in comparison to experimentally measured electromyographical signals of the muscles. Therefore, the electromyographical signals of 10 hand muscles of five test subjects performing seven different hand movements were measured. The kinematics of these tasks were used as input for the hand model, and the numerical muscle activities were computed. To analyze the relationship between simulated and measured activities, the time difference of the muscle on- and off-set points were calculated, which resulted in a mean on- and off-set time difference of 0.58 s between the experimental data and the model. The largest differences were detected for movements that mainly addressed the wrist. One major issue comparing simulated and measured muscle activities of the hand is cross-talk. Nevertheless, the results show that the hand model fits the experiment quite accurately despite some limitations and is a further step towards patient-specific modelling of the upper extremity.

The influence of shaft stiffness on joint kinematics and kinetics during hiking

Hiking boots provide an interface for walking in challenging environments, typically equipped with a shaft to provide ankle joint stability in rough terrains. Currently it is unclear if the ankle joint is stabilized to an extent that protects against ankle injuries, and if so, to what degree this added ankle stability sacrifices ankle mobility and hence decreases efficient gait propulsion. The aim of the present study was to compare the effect of shaft construction and stiffness on lower extremity kinematics and kinetics during level and step-down walking to simulate hiking conditions. Thirteen healthy males walked in one low-cut and three shafted commercially available hiking shoes with varying shaft stiffness. Lower extremity kinematics and ground reaction forces were recorded simultaneously. During level walking, ankle plantar-dorsiflexion range of motion was significantly reduced for the stiffest shaft hiking shoe compared to the low-cut shoe. A reduction in the muscle contribution to ankle joint work was found for all shafted shoes compared to the low-cut shoe. The reduced ankle joint work for the shafted shoes conversely increased eccentric knee joint work. Kinematic and kinetic differences between shoes diminished during box step-down walking. The present study shows that shaft height and stiffness can influence ankle joint range of motion, and ankle and knee joint work, with the high-shaft shoes redistributing load from the ankle to the knee joint. This may have implications for gait efficiency and increase the risk of knee joint loading or injuries.

Subject-Specific Modeling of Femoral Torsion Influences the Prediction of Hip Loading During Gait in Asymptomatic Adults

Hip osteoarthritis may be caused by increased or abnormal intra-articular forces, which are known to be related to structural articular cartilage damage. Femoral torsional deformities have previously been correlated with hip pain and labral damage, and they may contribute to the onset of hip osteoarthritis by exacerbating the effects of existing pathoanatomies, such as cam and pincer morphologies. A comprehensive understanding of the influence of femoral morphotypes on hip joint loading requires subject-specific morphometric and biomechanical data on the movement characteristics of individuals exhibiting varying degrees of femoral torsion. The aim of this study was to evaluate hip kinematics and kinetics as well as muscle and joint loads during gait in a group of adult subjects presenting a heterogeneous range of femoral torsion by means of personalized musculoskeletal models. Thirty-seven healthy volunteers underwent a 3D gait analysis at a self-selected walking speed. Femoral torsion was evaluated with low-dosage biplanar radiography. The collected motion capture data were used as input for an inverse dynamics analysis. Personalized musculoskeletal models were created by including femoral geometries that matched each subject’s radiographically measured femoral torsion. Correlations between femoral torsion and hip kinematics and kinetics, hip contact forces (HCFs), and muscle forces were analyzed. Within the investigated cohort, higher femoral antetorsion led to significantly higher anteromedial HCFs during gait (medial during loaded stance phase and anterior during swing phase). Most of the loads during gait are transmitted through the anterior/superolateral quadrant of the acetabulum. Correlations with hip kinematics and muscle forces were also observed. Femoral antetorsion, through altered kinematic strategies and different muscle activations and forces, may therefore lead to altered joint mechanics and pose a risk for articular damage. The method proposed in this study, which accounts for both morphological and kinematic characteristics, might help in identifying in a clinical setting patients who, as a consequence of altered femoral torsional alignment, present more severe functional impairments and altered joint mechanics and are therefore at a higher risk for cartilage damage and early onset of hip osteoarthritis.

Muscle forces and the demands of human walking

Reconstructing the locomotor behavior of extinct animals depends on elucidating the principles that link behavior, function, and morphology, which can only be done using extant animals. Within the human lineage, the evolution of bipedalism represents a critical transition, and evaluating fossil hominins depends on understanding the relationship between lower limb forces and skeletal morphology in living humans. As a step toward that goal, here we use a musculoskeletal model to estimate forces in the lower limb muscles of ten individuals during walking. The purpose is to quantify the consistency, timing, and magnitude of these muscle forces during the stance phase of walking. We find that muscles which act to support or propel the body during walking demonstrate the greatest force magnitudes as well as the highest consistency in the shape of force curves among individuals. Muscles that generate moments in the same direction as, or orthogonal to, the ground reaction force show lower forces of greater variability. These data can be used to define the envelope of load cases that need to be examined in order to understand human lower limb skeletal load bearing.

Summary: A musculoskeletal model of human walking reveals the consistency, timing, and magnitude of lower limb muscle forces across the stance phase.

The impact of anatomical uncertainties on the predictions of a musculoskeletal hand model - a sensitivity study

Outputs of musculoskeletal models should be considered probabilistic rather than deterministic as they are affected by inaccuracies and estimations associated with the development of the model. One of these uncertainties being critical for modeling arises from the determination of the muscles' line of action and the physiological cross-sectional area. Therefore, the aim of this study was to evaluate the outcome sensitivity of model predictions from a musculoskeletal hand model in comparison to the uncertainty of these input parameters. For this purpose, the kinematics and muscle activities of different hand movements (abduction of the fingers, abduction of the thumb, and flexion of the thumb) were recorded. One thousand simulations were calculated for each movement using the Latin hypercube sampling method with a corresponding variation of the muscle origin/insertion points and the cross-sectional area. Comparing the standard hand to simulations incorporating uncertainties of input parameters shows no major deviations in on- and off-set time point of muscle activities. About 60% of simulations are located within a ± 30% interval around the standard model concerning joint reaction forces. The comparison with the variation of the input data leads to the conclusion that the standard hand model is able to provide not over-scattered outcomes and, therefore, can be considered relatively stable. These results are of practical importance to the personalization of a musculoskeletal model with subject-specific bone geometries and hence changed muscle line of action.

Influence of simulated hip muscle weakness on hip joint forces during deep squatting

This study aimed to determine the effects of simulated hip muscle weakness on changes in hip joint forces during deep squat motion. Ten healthy individuals performed squat motion at three different positions (0° foot angle [N-squat], 10° toe-in [IN-squat], and 30° toe-out [OUT-squat]). A scaled musculoskeletal model for each participant was used to calculate the muscle and hip joint forces. For each hip muscle, models of full strength, mild muscle weakness (15% decrease), and severe muscle weakness (30% decrease) were created. The muscles affecting the hip joint forces were identified, and the rate of change in the joint forces was compared among the three squat conditions. The anterior hip joint force was increased in the muscle weakness models of the inferior gluteus maximus (iGlutMax) and iGlutMax+deep external rotator (ExtRot) muscles. With 30% muscle weakness of these muscles, statistically significant differences in the rate of increase in the anterior joint force were observed in the following order: IN-squat (iGlutMax, 29.5%; iGlutMax+ExtRot, 41.4%), N-squat (iGlutMax, 18.3%; iGlutMax+ExtRot, 27.8%), and OUT-squat (iGlutMax, 5.6%; iGlutMax+ExtRot, 9.3%). OUT-squat may be recommended to minimize the increase in hip joint forces if accompanied by hip muscle weakness.

Effect of the underlying cadaver data and patient-specific adaptation of the femur and pelvis on the prediction of the hip joint force estimated using static models

The prediction of the hip joint force (HJF) is a fundamental factor for the prevention of edge loading in total hip arthroplasty. Naturally, the loading of the liner of the acetabular component depends on the HJF acting on the artificial joint. In contrast to dynamic musculoskeletal models, static models for HJF prediction do not require motion analysis of the patient. However, patient-specific adaptability and validity of static models have to be scrutinized. In this study, a modular framework for HJF prediction using static models is introduced to compare the results of different cadaver templates that are the basis of most static and dynamic models, and different scaling laws for the patient-specific adaptation with in vivo HJF of ten patients for one-leg stance and level walking. The results revealed the significant effect of the underlying cadaver template used for the prediction of the HJF (p < 0.01). A higher degree of patient-specific scaling of the cadaver template often did not significantly reduce the prediction error. Three static models with the lowest prediction errors were compared to results of dynamic models from literature. The prediction error of the peak HJF of the static models (median absolute errors below 15% body weight in magnitude and below 5° in direction) was similar in magnitude and even smaller in direction compared to dynamic models. The necessary reduction of a load-based target zone for the prevention of edge loading due to the uncertainty of the HJF prediction has to be considered in the preoperative planning. The framework for HJF prediction is openly accessible at https://github.com/RWTHmediTEC/HipJointForceModel.

The effects of knee pain on knee contact force and external knee adduction moment in patients with knee osteoarthritis

Knee osteoarthritis (OA) is a major cause of knee pain, leading to physical dysfunction. External knee adduction moment (KAM), a surrogate measure of knee contact force (KCF) in the medial compartment, is related to knee pain, but the association between KCF and pain severity remains unclear. This study aimed to reveal the differences in KCF due to pain severity. Twenty-eight patients with knee OA were evaluated knee symptoms including pain severity via the Knee Society Score. Based on the median symptom score, 17 points in this study, subjects were classified as having Mild symptomatic OA (n = 15) and Severe symptomatic OA (n = 13). Subjects walked three times at a comfortable speed along a six-meter walkway, and we calculated KAM during the stance phase. KCF magnitude and distribution were also computed using the subject-specific musculoskeletal model, considering physical characteristics such as the femorotibial angle measured by X-ray. No differences in physical characteristics such as femorotibial angle and gait speed were found by symptom severity, whereas KAM and medial KCF at minimum and second peak in Severe symptomatic OA patients were significantly greater than those in Mild symptomatic OA. A significant medial shift of KCF in Severe symptomatic OA was also seen at first peak and minimum. Severe symptomatic OA had a greater medial KCF and medial shift of KCF. Detailed evaluations of KCF magnitude and distribution in addition to KAM would provide crucial information on knee contact force in relation to symptom severity.

Development and validation of a finite-element musculoskeletal model incorporating a deformable contact model of the hip joint during gait

Musculoskeletal models provide non-invasive and subject-specific biomechanical investigations of the musculoskeletal system. In a musculoskeletal model, muscle forces contribute to the deformation and kinematics of the joint which in turn would alter moment arms of muscles and ground reaction forces and thus affect the prediction of muscle forces and contact forces and contact mechanics of the joint. By far, deformable contact models of the hip have not been considered in musculoskeletal models, and the role of kinematics and deformation within the hip in muscle forces and hip contact mechanics is unknown. In this study, an FE musculoskeletal model including bones, joints and muscles of the lower extremity was developed. A deformable contact model of the hip joint was incorporated and coupled into the musculoskeletal model. Joint angles and ground reaction forces during gait were used as inputs. Optimization minimizing the sum of muscle stresses squared was performed directly to the FE musculoskeletal model in order to simultaneously solve muscle forces and contact forces and contact stresses of the hip joint within a single framework. The calculated hip contact forces corresponded well to the in vivo measurement data. The maximum hip contact stress was 6.5 MPa and occurred at weight-acceptance. The influence of kinematics and deformation in the hip on muscles forces and hip contact forces was minimal and not sensitive to variations in the thickness and properties of the joint cartilage during gait. This suggests that the uncoupled approach in which the hip contact forces and contact mechanics are simulated in separate frameworks would serve as an effective and efficient alternative for subject-specific modelling of the hip. This study provides guidance for the level of complexity needed for future hip models and can be used to evaluate biomechanical changes of the musculoskeletal system following interventions.

Automated Generation of Three-Dimensional Complex Muscle Geometries for Use in Personalised Musculoskeletal Models

The geometrical representation of muscles in computational models of the musculoskeletal system typically consists of a series of line segments. These muscle anatomies are based on measurements from a limited number of cadaveric studies that recently have been used as atlases for creating subject-specific models from medical images, so potentially restricting the options for personalisation and assessment of muscle geometrical models. To overcome this methodological limitation, we propose a novel, completely automated technique that, from a surface geometry of a skeletal muscle and its attachment areas, can generate an arbitrary number of lines of action (fibres) composed by a user-defined number of straight-line segments. These fibres can be included in standard musculoskeletal models and used in biomechanical simulations. This methodology was applied to the surfaces of four muscles surrounding the hip joint (iliacus, psoas, gluteus maximus and gluteus medius), segmented on magnetic resonance imaging scans from a cadaveric dataset, for which highly discretised muscle representations were created and used to simulate functional tasks. The fibres’ moment arms were validated against measurements and models of the same muscles from the literature with promising outcomes. The proposed approach is expected to improve the anatomical representation of skeletal muscles in personalised biomechanical models and finite element applications.

Current Preclinical Testing of New Hip Arthroplasty Technologies Does Not Reflect Real-World Loadings: Capturing Patient-Specific and Activity-Related Variation in Hip Contact Forces

BACKGROUND

Total hip arthroplasty (THA) implants are routinely tested for their tribological performance through regulatory preclinical wear testing (eg, ISO-14242). The standardized loading conditions defined in these tests consist of simplified waveforms, which do not specifically represent in vivo loads in different groups of patients. The aim of this study is to investigate, through musculoskeletal modeling, patient-specific and activity-related variation in hip contact forces (HCFs) in a large cohort of THA patients during common activities of daily living (ADLs).

METHODS

A total of 132 THA patients participated in a motion-capture analysis while performing different ADLs, including walk, fast walk, stair ascent, and descent (locomotor); sit to stand, stand to sit, squat, and lunge (nonlocomotor). HCFs were then calculated using the AnyBody Modeling System and qualitatively compared across all activities. The influence of gender on HCFs was analyzed through statistical parametric mapping analysis.

RESULTS

Systematic differences were found in HCF magnitudes and individual components in both locomotor and nonlocomotor ADLs. The qualitative analysis of the ADLs revealed a large range and a large variability in forces experienced at the hip during different activities. Significant differences in the 3-dimensional loading patterns were observed between males and females across most activities.

CONCLUSION

THA patients present a large variability in the forces experienced at the hip joint during their daily life. The interpatient variation might partially explain the heterogeneity observed in implant survival rates. A more extensive preclinical implant testing standard under clinically relevant loading conditions has been advocated to better predict and avoid clinical wear problems.

Patient characteristics affect hip contact forces during gait

OBJECTIVE

To examine hip contact force (HCF), calculated through multibody modelling, in a large total hip replacement (THR) cohort stratified by patient characteristics such as body mass index (BMI), age and function.

METHOD

132 THR patients undertook one motion capture session of gait analysis at a self-selected walking speed. HCFs were then calculated using the AnyBody Modelling System. Patients were stratified into three BMI groups, five age groups, and finally three functional groups determined by their self-selected gait speed. By means of statistical parametric mapping (SPM), statistical analyses of the 1-dimensional time series were performed to separately evaluate the influence of age, BMI and functionality on HCF.

RESULTS

The mean predicted HCFs were comparable to HCFs measured with instrumented prostheses reported in the literature. The SPM analysis revealed a statistically significant positive linear correlation between BMI and HCF, indicating that obese patients are more likely to experience higher HCF during most of the stance phase, while a statistically significant negative correlation with age was found only during the late swing-phase. Patients with higher functional ability exhibited significantly increased peak HCF, while patients with lower functional ability demonstrated lower HCFs overall and a pathological flattening of the typical double hump force profile.

CONCLUSION

HCFs experienced at the bearing surface are highly dependent on patient characteristics. BMI and functional ability were determined to have the biggest influence on contact forces. Current preclinical testing standards do not reflect this.

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.