The development of a segment-based musculoskeletal model of the lower limb: introducing FreeBody

An arm swing enhances the proximal-to-distal delay in joint extension during a countermovement jump

An abundance of degrees of freedom (DOF) exist when executing a countermovement jump (CMJ). This research aims to simplify the understanding of this complex system by comparing jump performance and independent functional DOF (fDOF) present in CMJs without (CMJNoArms) and with (CMJArms) an arm swing. Principal component analysis was used on 39 muscle forces and 15 3-dimensional joint contact forces obtained from kinematic and kinetic data, analyzed in FreeBody (a segment-based musculoskeletal model). Jump performance was greater in CMJArms with the increased ground contact time resulting in higher external (p = 0.012), hip (p < 0.001) and ankle (p = 0.009) vertical impulses, and slower hip extension enhancing the proximal-to-distal joint extension strategy. This allowed the hip muscles to generate higher forces and greater time-normalized hip vertical impulse (p = 0.006). Three fDOF were found for the muscle forces and 3-dimensional joint contact forces during CMJNoArms, while four fDOF were present for CMJArms. This suggests that the underlying anatomy provides mechanical constraints during a CMJ, reducing the demand on the control system. The additional fDOF present in CMJArms suggests that the arms are not mechanically coupled with the lower extremity, resulting in additional variation within individual motor strategies.

The Estimation of Knee Medial Force with Substitution Parameters during Walking and Turning

PURPOSE

Knee adduction, flexion moment, and adduction angle are often used as surrogate parameters of knee medial force. To verify whether these parameters are suitable as surrogates under different walking states, we investigated the correlation between knee medial loading with the surrogates during walking and turning.

METHODS

Sixteen healthy subjects were recruited to complete straight walk (SW), step turn (ST), and crossover turn (CT). Knee joint moments were obtained using inverse dynamics, and knee medial force was computed using a previously validated musculoskeletal model, Freebody. Linear regression was used to predict the peak of knee medial force with the peaks of the surrogate parameters and walking speed.

RESULTS

There was no significant difference in walking speed among these three tasks. The peak knee adduction moment (pKAM) was a significant predictor of the peak knee medial force (pKMF) for SW, ST, and CT (p < 0.001), while the peak knee flexion moment (pKFM) was only a significant predictor of the pKMF for SW (p = 0.034). The statistical analysis showed that the pKMF increased, while the pKFM and the peak knee adduction angle (pKAA) decreased significantly during CT compared to those of SW and ST (p < 0.001). The correlation analysis indicated that the knee parameters during SW and ST were quite similar.

CONCLUSIONS

This study investigated the relationship between knee medial force and some surrogate parameters during walking and turning. KAM was still the best surrogate parameter for SW, ST, and CT. It is necessary to consider the type of movement when comparing the surrogate predictors of knee medial force, as the prediction equations differ significantly among movement types.

Smartphone videos-driven musculoskeletal multibody dynamics modelling workflow to estimate the lower limb joint contact forces and ground reaction forces

The estimation of joint contact forces in musculoskeletal multibody dynamics models typically requires the use of expensive and time-consuming technologies, such as reflective marker-based motion capture (Mocap) system. In this study, we aim to propose a more accessible and cost-effective solution that utilizes the dual smartphone videos (SPV)-driven musculoskeletal multibody dynamics modeling workflow to estimate the lower limb mechanics. Twelve participants were recruited to collect marker trajectory data, force plate data, and motion videos during walking and running. The smartphone videos were initially analyzed using the OpenCap platform to identify key joint points and anatomical markers. The markers were used as inputs for the musculoskeletal multibody dynamics model to calculate the lower limb joint kinematics, joint contact forces, and ground reaction forces, which were then evaluated by the Mocap-based workflow. The root mean square error (RMSE), mean absolute deviation (MAD), and Pearson correlation coefficient (ρ) were adopted to evaluate the results. Excellent or strong Pearson correlations were observed in most lower limb joint angles (ρ = 0.74 ~ 0.94). The averaged MADs and RMSEs for the joint angles were 1.93 ~ 6.56° and 2.14 ~ 7.08°, respectively. Excellent or strong Pearson correlations were observed in most lower limb joint contact forces and ground reaction forces (ρ = 0.78 ~ 0.92). The averaged MADs and RMSEs for the joint lower limb joint contact forces were 0.18 ~ 1.07 bodyweight (BW) and 0.28 ~ 1.32 BW, respectively. Overall, the proposed smartphone video-driven musculoskeletal multibody dynamics simulation workflow demonstrated reliable accuracy in predicting lower limb mechanics and ground reaction forces, which has the potential to expedite gait dynamics analysis in a clinical setting.

Determining a musculoskeletal system's pre-stretched state using continuum-mechanical forward modelling and joint range optimization

The subject-specific range of motion (RoM) of a musculoskeletal joint system is balanced by pre-tension levels of individual muscles, which affects their contraction capability. Such an inherent pre-tension or pre-stretch of muscles is not measureable with in vivo experiments. Using a 3D continuum mechanical forward simulation approach for motion analysis of the musculoskeletal system of the forearm with 3 flexor and 2 extensor muscles, we developed an optimization process to determine the muscle fibre pre-stretches for an initial arm position, which is given human dataset. We used RoM values of a healthy person to balance the motion in extension and flexion. The performed sensitivity study shows that the fibre pre-stretches of the m. brachialis, m. biceps brachii and m. triceps brachii with 91 % dominate the objective flexion ratio, while m. brachiradialis and m. anconeus amount 7.8 % and 1.2 % . Within the multi-dimensional space of the surrogate model, 3D sub-spaces of primary variables, namely the dominant muscles and the global objective, flexion ratio, exhibit a path of optimal solutions. Within this optimal path, the muscle fibre pre-stretch of two flexors demonstrate a negative correlation, while, in contrast, the primary extensor, m. triceps brachii correlates positively to each of the flexors. Comparing the global optimum with four other designs along the optimal path, we saw large deviations, e.g., up to 15∘ in motion and up to 40% in muscle force. This underlines the importance of accurate determination of fibre pre-stretch in muscles, especially, their role in pathological muscular disorders and surgical applications such as free muscle or tendon transfer.

Hip joint and muscle loading for persons with bilateral transfemoral/through-knee amputations: biomechanical differences between full-length articulated and foreshortened non-articulated prostheses

BACKGROUND

Currently, there is little available in-depth analysis of the biomechanical effect of different prostheses on the musculoskeletal system function and residual limb internal loading for persons with bilateral transfemoral/through-knee amputations (BTF). Commercially available prostheses for BTF include full-length articulated prostheses (microprocessor-controlled prosthetic knees with dynamic response prosthetic feet) and foreshortened non-articulated stubby prostheses. This study aims to assess and compare the BTF musculoskeletal function and loading during gait with these two types of prostheses.

METHODS

Gait data were collected from four male traumatic military BTF and four able-bodied (AB) matched controls using a 10-camera motion capture system with two force plates. BTF completed level-ground walking trials with full-length articulated and foreshortened non-articulated stubby prostheses. Inverse kinematics, inverse dynamics and musculoskeletal modelling simulations were conducted.

RESULTS

Full-length articulated prostheses introduced larger stride length (by 0.5 m) and walking speed (by 0.3 m/s) than stubbies. BTF with articulated prostheses showed larger peak hip extension angles (by 10.1°), flexion moment (by 1.0 Nm/kg) and second peak hip contact force (by 3.8 bodyweight) than stubbies. There was no difference in the hip joint loading profile between BTF with stubbies and AB for one gait cycle. Full-length articulated prostheses introduced higher hip flexor muscle force impulse than stubbies.

CONCLUSIONS

Compared to stubbies, BTF with full-length articulated prostheses can achieve similar activity levels to persons without limb loss, but this may introduce detrimental muscle and hip joint loading, which may lead to reduced muscular endurance and joint degeneration. This study provides beneficial guidance in making informed decisions for prosthesis choice.

Bilateral Asymmetry in Knee and Hip Musculoskeletal Loading During Stair Ascending/Descending in Individuals with Unilateral Mild-to-Moderate Medial Knee Osteoarthritis

Most cases of unilateral knee osteoarthritis (OA) progress to bilateral OA within 10 years. Biomechanical asymmetries have been implicated in contralateral OA development; however, gait analysis alone does not consistently detect asymmetries in OA patient gait. Stair ambulation is a more demanding activity that may be more suited to reveal between-leg asymmetries in OA patients. The objective of this study was to investigate the between-leg biomechanical differences in patients with unilateral mild-to-moderate knee OA. Sixteen unilateral mild-to-moderate medial knee OA patients and 16 healthy individuals underwent kinematic and kinetic analysis of stair ascent and descent. Stair ascent produced higher loading and muscle forces in the unaffected limb compared to the OA limb, and stair descent produced lower loading on the OA limb compared to healthy subjects. These biomechanical differences were apparent in the ankle, knee, and hip joints. The implications of these findings are that OA patients rely more heavily on their unaffected sides than the affected side in stair ascent, a strategy that may be detrimental to the unaffected joint health. The reduction in affected limb loading in stair descent is thought to be related to minimizing pain.

The effect of muscle atrophy in people with unilateral transtibial amputation for three activities: Gait alone does not tell the whole story

Amputation imposes significant challenges in locomotion to millions of people with limb loss worldwide. The decline in the use of the residual limb results in muscle atrophy that affects musculoskeletal dynamics in daily activities. The aim of this study was to quantify the lower limb muscle volume discrepancy based on magnetic resonance (MR) imaging and to combine this with motion analysis and musculoskeletal modelling to quantify the effects in the dynamics of key activities of daily living. Eight male participants with traumatic unilateral transtibial amputation were recruited who were at least six months after receiving their definitive prostheses. The muscle volume discrepancies were found to be largest at the knee extensors (35 %, p = 0.008), followed by the hip abductors (17 %, p = 0.008). Daily activities (level walking, standing up from a chair and ascending one step) were measured in a motion analysis laboratory and muscle and joint forces quantified using a detailed musculoskeletal model for people with unilateral transtibial amputation which was calibrated in terms of the muscle volume discrepancies post-amputation at a subject-specific level. Knee extensor muscle forces were lower at the residual limb than the intact limb for all activities (p ≤ 0.008); residual limb muscle forces of the hip abductors (p ≤ 0.031) and adductors (p ≤ 0.031) were lower for standing-up and ascending one step. While the reduced knee extensor force has been reported by other studies, our results suggest a new biomechanically-based mitigation strategy to improve functional mobility, which could be achieved through strengthening of the hip abd/adductor muscles.

Characterization of muscle recruitment during gait of bilateral transfemoral and through-knee persons with limb loss

Introduction: Due to loss in musculoskeletal capacity, there is an increased burden on the residual limbs of bilateral transfemoral and through-knee persons with limb loss. This reduced capacity is associated with an increased cost of walking that is detrimental to functionality. Compensatory gait strategies are adopted by this population. However, how these strategies relate to specific muscle recruitment is not known. The primary aim of this study is to characterize muscle recruitment during gait of this population. The secondary aim is to assess whether the measured kinematics can be actuated when the endurance of specific muscles is reduced and if this is the case, which alternative muscles facilitate this. Methods: 3D gait data and high-resolution magnetic resonance images were acquired from six bilateral transfemoral and through-knee persons with limb loss. Subject-specific anatomical muscle models were developed for each participant, and a validated musculoskeletal model was used to quantify muscle forces in two conditions: during normal gait (baseline) and when muscles, which were identified as functioning above a "healthy" level at baseline, have a reduced magnitude of maximum force capacity (reduced endurance simulation). To test the hypothesis that there are differences in muscle forces between the baseline trials and the simulations with reduced muscular endurance, a Bonferroni corrected two-way ANOVA with repeated measures was completed between the two states. Results: The baseline analysis showed that the hip flexors experience relatively high muscle activations during gait. The reduced endurance simulation found two scenarios. First, for 5 out of the 12 simulations, the baseline kinematics could not be reproduced with the reduced muscular capacity. Second, for 7 out of 12 cases where the baseline kinematics were achieved, this was possible with compensatory increased activation of some muscles with similar functions (p ≤ 0.003). Discussion: Evidently, due to the loss of the ankle plantar flexors, gait imposes a high demand on the flexor muscle group of the residual limb. This study highlights how the elevated cost of gait in this population manifests in muscle recruitment. To enhance functionality, it is critical to consider the mechanical demand on the hip flexors and to develop rehabilitation interventions accordingly.

Prediction of in vivo hip contact forces during common activities of daily living using a segment-based musculoskeletal model

Background: Quantifying in vivo hip muscle and contact forces during activities of daily living (ADL) provides valuable information for diagnosis and treatment of hip-related disorders. The objective of this study was to utilize Freebody, a segment-based musculoskeletal model, for the prediction of hip contact forces using a novel objective function during seven common ADLs and validate its performance against the publicly available HIP98 dataset. Methods: Marker data, ground reaction forces, and hip contact forces during slow, normal, and fast walking, stair ascent and descent, and standing up and sitting down were extracted for 3 subjects from the HIP98 dataset. A musculoskeletal anatomical dataset was scaled to match the dimensions of each subject, and muscle and hip contact forces were estimated by minimizing a novel objective function, which was the summation of the muscle stresses squared and body weight-normalised hip contact force. The accuracy of predictions were quantified using several metrics, and muscle forces were qualitatively compared to experimental EMGs in the literature. Results: FreeBody predicted the hip contact forces during the ADLs with encouraging accuracy: The root mean squared error of predictions were 44.0 ± 8.5, 47.4 ± 6.5, and 59.8 ± 7.1% BW during slow, normal, and fast walking, 44.2 ± 16.8% and 53.3 ± 12.2% BW for stair ascent and descent, and 31.8 ± 8.2% and 17.1 ± 5.0% BW for standing up and sitting down, respectively. The error in prediction of peak hip contact forces were 14-18%, 24-28%, 17-35% for slow, normal, and fast walking, 7-25% and 15-32% in stair ascent and descent, and around 10% for standing up and sitting down. The coefficient of determination was larger than 0.90 in all activities except in standing up (0.86 ± 0.08). Conclusion: This study has implemented a novel objective function in a segment-based musculoskeletal model, FreeBody, for the prediction of hip contact forces during a large range of ADLs. The model outputs compare favourably for all ADLs and are the best in standing up and sitting down, while muscle activation patterns are consistent with experimental EMGs from literature. This new objective function addresses one of the major limitations associated with musculoskeletal models in the literature, namely the high non-physiological predicted hip joint contact forces.

Modelling and in vivo evaluation of tendon forces and strain in dynamic rehabilitation exercises: a scoping review

OBJECTIVES

Although exercise is considered the preferred approach for tendinopathies, the actual load that acts on the tendon in loading programmes is usually unknown. The objective of this study was to review the techniques that have been applied in vivo to estimate the forces and strain that act on the human tendon in dynamic exercises used during rehabilitation.

DESIGN

Scoping review.

DATA SOURCES

Embase, PubMed, Web of Science and Google Scholar were searched from database inception to February 2021.

ELIGIBILITY CRITERIA

Cross-sectional studies available in English or Spanish language were included if they focused on evaluating the forces or strain of human tendons in vivo during dynamic exercises. Studies were excluded if they did not evaluate tendon forces or strain; if they evaluated running, walking, jumping, landing or no dynamic exercise at all; and if they were conference proceedings or book chapters.

DATA EXTRACTION AND SYNTHESIS

Data extracted included year of publication, study setting, study population characteristics, technique used and exercises evaluated. The studies were grouped by the types of techniques and the tendon location.

RESULTS

Twenty-one studies were included. Fourteen studies used an indirect methodology based on inverse dynamics, nine of them in the Achilles and five in the patellar tendon. Six studies implemented force transducers for measuring tendon forces in open carpal tunnel release surgery patients. One study applied an optic fibre technique to detect forces in the patellar tendon. Four studies measured strain using ultrasound-based techniques.

CONCLUSIONS

There is a predominant use of inverse dynamics, but force transducers, optic fibre and estimations from strain data are also used. Although these tools may be used to make general estimates of tendon forces and strains, the invasiveness of some methods and the loss of immediacy of others make it difficult to provide immediate feedback to the individuals.

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.

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.

Maximum dorsiflexion increases Achilles tendon force during exercise for midportion Achilles tendinopathy

Rehabilitation is an important treatment for non-insertional Achilles tendinopathy. To date, eccentric loading exercises (ECC) have been the predominant choice; however, mechanical evidence underlying their use remains unclear. Other protocols, such as heavy slow resistance loading (HSR), have shown comparable outcomes, but with less training time. This study aims to identify the effect of external loading and other variables that influence Achilles tendon (AT) force in ECC and HSR. Ground reaction force and kinematic data during ECC and HSR were collected from 18 healthy participants for four loading conditions. The moment arms of the AT were estimated from MRIs of each participant. AT force then was calculated using the ankle torque obtained from inverse dynamics. In the eccentric phase, the AT force was not larger than in the concentric phase in both ECC and HSR. Under the same external load, the force through the AT was larger in ECC with the knee bent than in HSR with the knee straight due to increased dorsiflexion angle of the ankle. Multivariate regression analysis showed that external load and maximum dorsiflexion angle were significant predictors of peak AT force in both standing and seated positions. Therefore, to increase the effectiveness of loading the AT, exercises should apply adequate external load and reach maximum dorsiflexion during the movement. Peak dorsiflexion angle affected the AT force in a standing position at twice the rate of a seated position, suggesting standing could prove more effective for the same external loading and peak dorsiflexion angle.

Higher knee contact forces might underlie increased osteoarthritis rates in high functioning amputees: A pilot study

High functioning military transtibial amputees (TTAs) with well-fitted state of the art prosthetics have gait that is indistinguishable from healthy individuals, yet they are more likely to develop knee osteoarthritis (OA) of their intact limbs. This contrasts with the information at the knees of the amputated limbs that have been shown to be at a significantly reduced risk of pain and OA. The hypothesis of this study is that biomechanics can explain the difference in knee OA risk. Eleven military unilateral TTAs and eleven matched healthy controls underwent gait analysis. Muscle forces and joint contact forces at the knee were quantified using musculoskeletal modeling, validated using electromyography measurements. Peak knee contact forces for the intact limbs on both the medial and lateral compartments were significantly greater than the healthy controls (P ≤ .006). Additionally, the intact limbs had greater peak semimembranosus (P = .001) and gastrocnemius (P ≤ .001) muscle forces compared to the controls. This study has for the first time provided robust evidence of increased force on the non-affected knees of high functioning TTAs that supports the mechanically based hypothesis to explain the documented higher risk of knee OA in this patient group. The results suggest several protentional strategies to mitigate knee OA of the intact limbs, which may include the improvements of the prosthetic foot control, socket design, and strengthening of the amputated muscles.

Review of musculoskeletal modelling in a clinical setting: Current use in rehabilitation design, surgical decision making and healthcare interventions

BACKGROUND

Musculoskeletal modelling is a common means by which to non-invasively analyse movement. Such models have largely been used to observe function in both healthy and patient populations. However, utility in a clinical environment is largely unknown. The aim of this review was to explore existing uses of musculoskeletal models as a clinical intervention, or decision-making, tool.

METHODS

A literature search was performed using PubMed and Scopus to find articles published since 2010 and relating to musculoskeletal modelling and joint and muscle forces.

FINDINGS

4662 abstracts were found, of which 39 relevant articles were reviewed. Journal articles were categorised into 5 distinct groups: non-surgical treatment, orthoses assessment, surgical decision making, surgical intervention assessment and rehabilitation regime assessment. All reviewed articles were authored by collaborations between clinicians and engineers/modellers. Current uses included insight into the development of osteoarthritis, identifying candidates for hamstring lengthening surgery, and the assessment of exercise programmes to reduce joint damage.

INTERPRETATION

There is little evidence showing the use of musculoskeletal modelling as a tool for patient care, despite the ability to assess long-term joint loading and muscle overuse during functional activities, as well as clinical decision making to avoid unfavourable treatment outcomes. Continued collaboration between model developers should aim to create clinically-friendly models which can be used with minimal input and experience by healthcare professionals to determine surgical necessity and suitability for rehabilitation regimes, and in the assessment of orthotic devices.

Differences in Motor Control Strategies of Jumping Tasks, as Revealed by Group and Individual Analysis

The aim of this study was to investigate the motor control strategies adopted when performing two jumping tasks with different task demands when analysed at an individual and group level. Twenty-two healthy individuals performed two jumping tasks: jumping without the use of an arm swing (CMJnas) and jumping starting in a plantar flexed position with the use of an arm swing (PF). Principal component analysis (PCA) was performed using hip, knee and ankle joint moment data on individual (PCAi) and group data (PCAc). The results demonstrate that a greater number of PCs are required to explain the majority of variance within the dataset in the PF condition at both an individual and group level, compared to CMJnas condition. Although common control strategies were observed between the two jumping conditions, differences in the organisation of the movement (PC loading coefficients) were observed. Results from the group analysis did not completely reflect the individual strategies used to perform each jumping task and highlight the value in performing individual analysis to determine emergent control strategies.

Task Demand Changes Motor Control Strategies in Vertical Jumping

The purpose of this study was to examine the motor control strategies employed to control the degrees of freedom when performing a lower limb task with constraints applied at the hip, knee, and ankle. Thirty-five individuals performed vertical jumping tasks: hip flexed, no knee bend, and plantar flexed. Joint moment data from the hip, knee, and ankle were analyzed using principal component analysis (PCA). In all PCA performed, a minimum of two and maximum of six principal components (PC) were required to describe the movements. Similar reductions in dimensionality were observed in the hip flexed and no knee bend conditions (3PCs), compared to the plantar flexed condition (5PCs). A proximal to distal reduction in variability was observed for the hip flexed and no knee bend conditions but not for the plantar flexed condition. Collectively, the results suggest a reduction in the dimensionality of the movement occurs despite the constraints imposed within each condition and would suggest that dimensionality reduction and motor control strategies are a function of the task demands.

A Femoral Clamp to Reduce Soft Tissue Artifact: Accuracy and Reliability in Measuring Three-Dimensional Knee Kinematics During Gait

The accurate measurement of full six degrees-of-freedom (6DOFs) knee joint kinematics is prohibited by soft tissue artifact (STA), which remains the greatest source of error. The purpose of this study was to present and assess a new femoral clamp to reduce STA at the thigh. It was hypothesized that the device can preserve the natural knee joint kinematics pattern and outperform a conventional marker mounted rigid cluster during gait. Six healthy subjects were asked to walk barefoot on level ground with a cluster marker set (cluster gait) followed by a cluster-clamp-merged marker set (clamp gait) and their kinematics was measured using the cluster method in cluster gait and the cluster and clamp methods simultaneously in clamp gait. Two operators performed the gait measurement. A 6DOFs knee joint model was developed to enable comparison with the gold standard knee joint kinematics measured using a dual fluoroscopic imaging technique. One-dimensional (1D) paired t-tests were used to compare the knee joint kinematics waveforms between cluster gait and clamp gait. The accuracy was assessed in terms of the root-mean-square error (RMSE), coefficient of determination, and Bland-Altman plots. Interoperator reliability was assessed using the intraclass correlation coefficient (ICC). The result showed that the femoral clamp did not change the walking speed and knee joint kinematics waveforms. Additionally, clamp gait reduced the rotation and translation errors in the transverse plane and improved the interoperator reliability when compared to the rigid cluster method, suggesting a more accurate and reliable measurement of knee joint kinematics.

Reconstruction of the lower limb bones from digitised anatomical landmarks using statistical shape modelling

BACKGROUND

Bone shapes strongly influence force and moment predictions of kinematic and musculoskeletal models used in motion analysis. The precise determination of joint reference frames is essential for accurate predictions. Since clinical motion analysis typically does not include medical imaging, from which bone shapes may be obtained, scaling methods using reference subjects to create subject-specific bone geometries are widely used.

RESEARCH QUESTION

This study investigated if lower limb bone shape predictions from skin-based measurements, utilising an underlying statistical shape model (SSM) that corrects for soft tissue artefacts in digitisation, can be used to improve conventional linear scaling methods of bone geometries.

METHODS

SSMs created from 35 healthy adult femurs and tibiae/fibulae were used to reconstruct bone shapes by minimising the distance between anatomical landmarks on the models and those digitised in the motion laboratory or on medical images. Soft tissue artefacts were quantified from magnetic resonance images and then used to predict distances between landmarks digitised on the skin surface and bone. Reconstruction results were compared to linearly scaled models by measuring root mean squared distances to segmented surfaces, calculating differences of commonly used anatomical measures and the errors in the prediction of the hip joint centre.

RESULTS

SSM reconstructed surface predictions from varying landmark sets from skin and bone landmarks were more accurate compared to linear scaling methods (2.60-2.95 mm vs. 3.66-3.87 mm median error; p < 0.05). No significant differences were found between SSM reconstructions from bony landmarks and SSM reconstructions from digitised landmarks obtained in the motion lab and therefore reconstructions using skin landmarks are as accurate as reconstructions from landmarks obtained from medical images.

SIGNIFICANCE

These results indicate that SSM reconstructions can be used to increase the accuracy in obtaining bone shapes from surface digitised experimental data acquired in motion lab environments.

A musculoskeletal modelling approach to explain sit-to-stand difficulties in older people due to changes in muscle recruitment and movement strategies

By 2050 the proportion of over 65s is predicted to be 20% of the population. The consequences of an age-related reduction in muscle mass have not been fully investigated and, therefore, the aim of the present study was to quantify the muscle and joint contact forces using musculoskeletal modelling, during a sit-to-stand activity, to better explain difficulties in performing everyday activities for older people. A sit-to-stand activity with and without the use of arm rests was observed in ninety-five male participants, placed into groups of young (aged 18-35 years), middle-aged (aged 40-60 years) or older adults (aged 65 years and over). Older participants demonstrated significantly lower knee extensor and joint forces than the young when not using arm rests, compensating through elevated hip extensor and ankle plantarflexor muscle activity. The older group were also found to have higher shoulder joint contact forces whilst using arm rests. This tendency to reorganise muscle recruitment to include neighbouring groups or other parts of the body could make everyday activities more susceptible to age-related functional decline. Reductions in leg strength, via age- or atrophy- related means, creates increased reliance on the upper body and may result in further lower limb atrophy through disuse. The eventual decline of upper body function reduces strength reserves, leading to increased vulnerability, dependence on others and risk of institutionalisation.

Micromotion and Push-Out Evaluation of an Additive Manufactured Implant for Above-the-Knee Amputees

In comparison to through-knee amputees the outcomes for above-the-knee amputees are relatively poor; based on this novel techniques have been developed. Most current percutaneous implant-based solutions for transfemoral amputees make use of high stiffness intramedullary rods for skeletal fixation, which can have risks including infection, femoral fractures, and bone resorption due to stress shielding. This work details the cadaveric testing of a short, cortical bone stiffness-matched subcutaneous implant, produced using additive manufacture, to determine bone implant micromotion and push-out load. The results for the micromotions were all <20 μm and the mean push-out load was 2,099 Newtons. In comparison to a solid control, the stiffness-matched implant exhibited significantly higher micromotion distributions and no significant difference in terms of push-out load. These results suggest that, for the stiffness-matched implant at time zero, osseointegration would be facilitated and that the implant would be securely anchored. For these metrics, this provides justification for the use of a short-stem implant for transfemoral amputees in this subcutaneous application. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 37:2104-2111, 2019.

Knee Forces During Landing in Men and Women

Sex differences in biomechanics may provide one explanation for the greater incidence of knee injuries in women, but few studies have compared internal forces. In this study, a musculoskeletal model was used to compare male and female, bilateral and unilateral landings based on motion capture and force plate data. Participants were classified as landing medially or laterally loaded based upon the mediolateral load share at the knee (bilateral: p < 0.001, η²=0.452; unilateral: p < 0.001, η² = 0.444). Knee kinematics and ground reaction forces were not different between the two groups (p > 0.05, η² = 0.001 - 0.059), but there were differences in muscular recruitment. Landing strategy did not appear to be dependent on sex. However, for both medially and laterally loaded bilateral landings men had greater gluteal (p = 0.017, η² = 0.085) and hamstrings forces (p < 0.001, η² = 0.183), whereas women had greater quadriceps forces (p = 0.004, η² = 0.116). This study demonstrates an association between muscular recruitment and medially loaded landings. Landing strategy seems to be a function of skill not sex; however, within a particular landing strategy there may be sex differences in muscular activation that contribute to the difference in injury rates.

Principal Component Analysis Reveals the Proximal to Distal Pattern in Vertical Jumping Is Governed by Two Functional Degrees of Freedom

The successful completion of motor tasks requires effective control of multiple degrees of freedom (DOF), with adaptations occurring as a function of varying performance constraints. In this study we sought to compare the emergent coordination strategies employed in vertical jumping under different task constraints [countermovement jump (CMJ) with arm swing-CMJas and no arm swing-CMJnas]. In order to achieve this, principal component analysis (PCA) was conducted on joint moment waveform data from the hip, knee and ankle. This statistical approach has the advantage of analyzing the whole movement within a time series and reduces multidimensional datasets to lower dimensions for analysis. Both individual and group analyses were conducted. For individual analysis, PCA was conducted on combined hip, knee, and ankle joint moment data for each individual across both CMJnas (thirty-eight participants), and CMJas (twenty-two participants) conditions. PCA was also performed comparing all data from each individual across CMJnas and CMJas conditions. The results revealed a maximum of three principal components (PC) explained over 90% of the variance in the data sets for both conditions and within individual and group analyses. For individual analysis, no more than 2PCs were required for both conditions. For group analysis, CMJas required 3PCs to explain over 90% of the variance within the dataset and CMJnas only required 2PCs. Reconstruction of the original NJM waveforms from the PCA output demonstrates a greater loading of hip and knee joint moments to PC1, with PC2 showing a greater loading to ankle joint moment. The reduction in dimensions of the original data shows the proximal to distal extension pattern in the sagittal plane, typical of vertical jumping tasks, is governed by only 2 functional DOF, at both a group, and individual level, rather than the typically reported 3 mechanical DOF in some forms of jumping.

Improving Musculoskeletal Model Scaling Using an Anatomical Atlas: The Importance of Gender and Anthropometric Similarity to Quantify Joint Reaction Forces

OBJECTIVE

The accuracy of a musculoskeletal model relies heavily on the implementation of the underlying anatomical dataset. Linear scaling of a generic model, despite being time and cost efficient, produces substantial errors as it does not account for gender differences and inter-individual anatomical variations. The hypothesis of this study is that linear scaling to a musculoskeletal model with gender and anthropometric similarity to the individual subject produces similar results to the ones that can be obtained from a subject-specific model.

METHODS

A lower limb musculoskeletal anatomical atlas was developed consisting of ten datasets derived from magnetic resonance imaging of healthy subjects and an additional generic dataset from the literature. Predicted muscle activation and joint reaction force were compared with electromyography and literature data. Regressions based on gender and anthropometry were used to identify the use of atlas.

RESULTS

Primary predictors of differences for the joint reaction force predictions were mass difference for the ankle (p < 0.001) and length difference for the knee and hip (p ≤ 0.017). Gender difference accounted for an additional 3% of the variance (p ≤ 0.039). Joint reaction force differences at the ankle, knee, and hip were reduced by between 50% and 67% (p = 0.005) when using a musculoskeletal model with the same gender and similar anthropometry in comparison with a generic model.

CONCLUSION

Linear scaling with gender and anthropometric similarity can improve joint reaction force predictions in comparison with a scaled generic model.

SIGNIFICANCE

The presented scaling approach and atlas can improve the fidelity and utility of musculoskeletal models for subject-specific applications.

Extra excitation of biceps femoris during neuromuscular electrical stimulation reduces knee medial loading

Medial knee joint osteoarthritis (OA) is a debilitating and prevalent condition. Surgical treatment consists of redistributing the forces from the medial to the lateral compartment through osteotomy, or replacing the joint surfaces. As the mediolateral load distribution is related to the action of the musculature around the knee, the aim of this study was to devise a technique to redistribute these forces non-surgically through changes in muscle excitation. Eight healthy subjects participated in the experiment, and neuromuscular electrical stimulation was used to change the muscle forces around the knee. A musculoskeletal model was used to quantify the loading on the medial compartment of the knee, and a novel algorithm devised and implemented to simulate neuromuscular electrical stimulation. The forces and moments at the knee, ground reaction forces, walking velocity and step length were quantified before and after stimulation. Stimulation of the biceps femoris resulted in a significant decrease in the second peak of the medial knee joint loading by up to 0.17 body weight (p = 0.016). Kinematic parameters were not significantly affected. Neuromuscular electrical stimulation can decrease the peak loads on the medial compartment of the knee, and thus offers a promising therapy for medial knee joint OA.

Deep Learning for Musculoskeletal Force Prediction

Musculoskeletal models permit the determination of internal forces acting during dynamic movement, which is clinically useful, but traditional methods may suffer from slowness and a need for extensive input data. Recently, there has been interest in the use of supervised learning to build approximate models for computationally demanding processes, with benefits in speed and flexibility. Here, we use a deep neural network to learn the mapping from movement space to muscle space. Trained on a set of kinematic, kinetic and electromyographic measurements from 156 subjects during gait, the network's predictions of internal force magnitudes show good concordance with those derived by musculoskeletal modelling. In a separate set of experiments, training on data from the most widely known benchmarks of modelling performance, the international Grand Challenge competitions, generates predictions that better those of the winning submissions in four of the six competitions. Computational speedup facilitates incorporation into a lab-based system permitting real-time estimation of forces, and interrogation of the trained neural networks provides novel insights into population-level relationships between kinematic and kinetic factors.

An important role of the biarticular hamstrings is to exert internal/external rotation moments on the tibia during vertical jumping

Most research considering biarticular muscle function has tended to focus on the sagittal plane. Instead, the purpose of this study was to evaluate the internal/external rotation moment arms of the biarticular muscles of the knee, and then to explore their function. The FreeBody musculoskeletal model of the lower limb was used to calculate the moment arms and moments that each of the muscles of the knee exerted on the proximal tibia of 12 athletic males during vertical jumping. Biceps femoris and tensor fascia latae were external rotators of the tibia, whereas semimembranosus, semitendinosus, sartorius, gracilis, popliteus and the patellar tendon were internal rotators. The magnitudes of the internal/external rotation and flexion moments exerted on the tibia by the biarticular hamstrings were similar, suggesting that the creation of internal/external rotation is a key aspect of their role. One potential reason is to stabilise the tibia during femoral extension (and it is argued that it may be helpful to characterise the creation of active joint stability as the stabilisation of one segment during the rotation of an adjacent segment). A second explanation may be to mechanically couple hip abduction when the hip is flexed with internal rotation of the tibia.

Effects of an 8-week strength training intervention on tibiofemoral joint loading during landing: a cohort study

Objectives

To use a musculoskeletal model of the lower limb to evaluate the effect of a strength training intervention on the muscle and joint contact forces experienced by untrained women during landing.

Methods

Sixteen untrained women between 18 and 28 years participated in this cohort study, split equally between intervention and control groups. The intervention group trained for 8 weeks targeting improvements in posterior leg strength. The mechanics of bilateral and unilateral drop landings from a 30 cm platform were recorded preintervention and postintervention, as was the isometric strength of the lower limb during a hip extension test. The internal muscle and joint contact forces were calculated using FreeBody, a musculoskeletal model.

Results

The strength of the intervention group increased by an average of 35% (P<0.05; pre: 133±36 n, post: 180±39 n), whereas the control group showed no change (pre: 152±36 n, post: 157±46 n). There were only small changes from pre-test to post-test in the kinematics and ground reaction forces during landing that were not statistically significant. Both groups exhibited a post-test increase in gluteal muscle force during landing and a lateral to medial shift in tibiofemoral joint loading in both landings. However, the magnitude of the increase in gluteal force and lateral to medial shift was significantly greater in the intervention group.

Conclusion

Strength training can promote a lateral to medial shift in tibiofemoral force (mediated by an increase in gluteal force) that is consistent with a reduction in valgus loading. This in turn could help prevent injuries that are due to abnormal knee loading such as anterior cruciate ligament ruptures, patellar dislocation and patellofemoral pain.

Activation of biceps femoris long head reduces tibiofemoral anterior shear force and tibial internal rotation torque in healthy subjects

The anterior cruciate ligament (ACL) provides resistance to tibial internal rotation torque and anterior shear at the knee. ACL deficiency results in knee instability. Optimisation of muscle contraction through functional electrical stimulation (FES) offers the prospect of mitigating the destabilising effects of ACL deficiency. The hypothesis of this study is that activation of the biceps femoris long head (BFLH) reduces the tibial internal rotation torque and the anterior shear force at the knee. Gait data of twelve healthy subjects were measured with and without the application of FES and taken as inputs to a computational musculoskeletal model. The model was used to investigate the optimum levels of BFLH activation during FES gait in reducing the anterior shear force to zero. This study found that FES significantly reduced the tibial internal rotation torque at the knee during the stance phase of gait (p = 0.0322) and the computational musculoskeletal modelling revealed that a mean BFLH activation of 20.8% (±8.4%) could reduce the anterior shear force to zero. At the time frame when the anterior shear force was zero, the internal rotation torque was reduced by 0.023 ± 0.0167 Nm/BW, with a mean 188% reduction across subjects (p = 0.0002). In conclusion, activation of the BFLH is able to reduce the tibial internal rotation torque and the anterior shear force at the knee in healthy control subjects. This should be tested on ACL deficient subject to consider its effect in mitigating instability due to ligament deficiency. In future clinical practice, activating the BFLH may be used to protect ACL reconstructions during post-operative rehabilitation, assist with residual instabilities post reconstruction, and reduce the need for ACL reconstruction surgery in some cases.

Reliability and Minimal Detectable Change Values for Predictions of Knee Forces during Gait and Stair Ascent Derived from the FreeBody Musculoskeletal Model of the Lower Limb

FreeBody is a musculoskeletal model of the lower limb used to calculate predictions of muscle and joint contact forces. The validation of FreeBody has been described in a number of publications; however, its reliability has yet to be established. The purpose of this study was, therefore, to establish the test-retest reliability of FreeBody in a population of healthy adults in order to add support to previous and future research using FreeBody that demonstrates differences between cohorts after an intervention. We hypothesized that test-retest estimations of knee contact forces from FreeBody would demonstrate a high intra-class correlation. Kinematic and kinetic data from nine older participants (4 men: mean age = 63 ± 11 years; 5 women: mean age = 49 ± 4 years) performing level walking and stair ascent was collected on consecutive days and then analyzed using FreeBody. There was a good level of intra-session agreement between the waveforms for the individual trials of each activity during testing session 1 (R = 0.79-0.97). Similarly, overall there was a good inter-session agreement within subjects (R = 0.69-0.97) although some subjects showed better agreement than others. There was a high level of agreement between the group mean waveforms of the two sessions for all variables (R = 0.882-0.997). The intra-class correlation coefficients (ICC) were very high for peak tibiofemoral joint contact forces (TFJ) and hamstring forces during gait, for peak patellofemoral joint contact forces and quadriceps forces during stair ascent and for peak lateral TFJ and the proportion of TFJ accounted for by the medial compartment during both tasks (ICC = 0.86-0.96). Minimal detectable change (MDC) of the peak knee forces during gait ranged between 0.43 and 1.53 × body weight (18-170% of the mean peak values). The smallest MDCs were found for medial TFJ share (4.1 and 5.8% for walking and stair ascent, respectively, or 4.8 and 6.7% of the mean peak values). In conclusion, the results of this study support the use of FreeBody to investigate the effect of interventions on muscle and joint contact forces at the cohort level, but care should be taken if using FreeBody at the subject level.

Effect of a gluteal activation warm-up on explosive exercise performance

Objectives

To evaluate the effect of a gluteal activation warm-up on the performance of an explosive exercise (the high hang pull (HHP)).

Methods

Seventeen professional rugby union players performed one set of three HHPs (with 80% of their one repetition maximum load) following both a control and activation warm-up. Peak electrical activity of the gluteus maximus and medius was quantified using electromyography (EMG). In addition, the kinematics and kinetics of nine players was also recorded using force plate and motion capture technology. These data were analysed using a previously described musculoskeletal model of the right lower limb in order to provide estimates of the muscular force expressed during the movement.

Results

The mean peak EMG activity of the gluteus maximus was significantly lower following the activation warm-up as compared with the control (p<0.05, effect size d=0.30). There were no significant differences in the mean peak estimated forces in gluteus maximus and medius, the quadriceps or hamstrings (p=0.053), although there was a trend towards increased force in gluteus maximus and hamstrings following the activation warm-up. There were no differences between the ground reaction forces following the two warm-ups.

Conclusion

This study suggests that a gluteal activation warm-up may facilitate recruitment of the gluteal musculature by potentiating the glutes in such a way that a smaller neural drive evokes the same or greater force production during movement. This could in turn potentially improve movement quality.

RIGID MUSCULOSKELETAL MODELS OF THE HUMAN BODY SYSTEMS: A REVIEW

Advancing the knowledge of the biomechanics of the human body is essential to improve the clinical decision-makings of musculoskeletal disorders in the framework of in silico medicine. An impressive number of research projects focused on the development of rigid-body musculoskeletal models have been conducted over the world thanks to the new research directives. However, the application of these models in clinical practices remains a challenging issue. The objective of this review paper was to present the most current rigid-body musculoskeletal models of the human body systems and to analyze their trends and weaknesses for clinical applications. Then, recommendations were proposed for future researches toward fully clinical decision support. A systematic review process was performed. Well-selected studies related to the most current rigid-body 3D musculoskeletal models for each body system component (jaw, cervical spine, upper limbs, lumbar spine, and lower limbs) were summarized and explored. Trends in rigid musculoskeletal modeling are highlighted as personalization, new imaging techniques for specific joint kinematics, and computational efficiency. Weaknesses are highlighted as modeling assumptions, use of generic model, lack of modeling consensus, model validation, and parameter and model uncertainties. Future directions related to joint and muscle modeling, neuro-musculoskeletal modeling, model validation, data and model uncertainty quantification are recommended.

Functional electrical stimulation of gluteus medius reduces the medial joint reaction force of the knee during level walking

Background

By altering muscular activation patterns, internal forces acting on the human body during dynamic activity may be manipulated. The magnitude of one of these forces, the medial knee joint reaction force (JRF), is associated with disease progression in patients with early osteoarthritis (OA), suggesting utility in its targeted reduction. Increased activation of gluteus medius has been suggested as a means to achieve this.

Methods

Motion capture equipment and force plate transducers were used to obtain kinematic and kinetic data for 15 healthy subjects during level walking, with and without the application of functional electrical stimulation (FES) to gluteus medius. Musculoskeletal modelling was employed to determine the medial knee JRF during stance phase for each trial. A further computer simulation of increased gluteus medius activation was performed using data from normal walking trials by a manipulation of modelling parameters. Relationships between changes in the medial knee JRF, kinematics and ground reaction force were evaluated.

Results

In simulations of increased gluteus medius activity, the total impulse of the medial knee JRF was reduced by 4.2 % (p = 0.003) compared to control. With real-world application of FES to the muscle, the magnitude of this reduction increased to 12.5 % (p < 0.001), with significant inter-subject variation. Across subjects, the magnitude of reduction correlated strongly with kinematic (p < 0.001) and kinetic (p < 0.001) correlates of gluteus medius activity.

Conclusions

The results support a major role for gluteus medius in the protection of the knee for patients with OA, establishing the muscle’s central importance to effective therapeutic regimes. FES may be used to achieve increased activation in order to mitigate distal internal loads, and much of the benefit of this increase can be attributed to resulting changes in kinematic parameters and the ground reaction force. The utility of interventions targeting gluteus medius can be assessed in a relatively straightforward way by determination of the magnitude of reduction in pelvic drop, an easily accessed marker of aberrant loading at the knee.