A 3D lower limb musculoskeletal model for simultaneous estimation of musculo-tendon, joint contact, ligament and bone forces during gait

Automatic generation of knee kinematic models from medical imaging

BACKGROUND AND OBJECTIVE

Three-dimensional spatial mechanisms have been used to accurately predict passive knee kinematics, and have shown potential to be used in optimized multibody kinematic models. Such multi-body models are anatomically consistent and can accurately predict passive knee kinematics, but require extensive medical image processing and thus are not widely adopted. This study aimed to automate the generation of kinematic models of tibiofemoral (TFJ) and patellofemoral (PFJ) joints from segmented magnetic resonance imaging (MRI) and compare them against a corresponding manual pipeline.

METHODS

From segmented MRI of eight healthy participants (four females; aged 14.0 ± 2.6 years), geometric parameters (i.e., articular surfaces, ligament attachments) were determined both automatically and manually, and then assembled into TFJ and PFJ kinematic models to predict passive kinematics. The TFJ model was a six-link mechanism with deformable ligamentous constraints, whereas PFJ was a modified hinge. The ligament length changes through TFJ flexion were prescribed to literature strain profile. The geometric parameters were optimized to ensure physiological kinematic predictions through a Multiple Objective Particle Swarm Optimization.

RESULTS

Geometric parameters showed strong agreement between automatic and manual pipelines (median error of 2.8 mm for anatomical landmarks and 1.5 mm for ligament lengths). Predicted TFJ and PFJ kinematics from the two pipelines were not statistically different, except for tibial superior/inferior translation near terminal TFJ extension. The TFJ kinematics predicted from the automatic pipeline had mean errors of 3.6° and 12.4° for adduction/abduction and internal/external rotation, respectively, and <7 mm mean translational error compared to the manual pipeline. Predicted PFJ had <9° mean rotational errors and <6 mm mean translational errors.

CONCLUSIONS

The automatic pipeline developed and presented here can predict passive knee kinematics comparable to a manual pipeline, but removes laborious manual processing and provides a systematic approach to model creation.

Shape Memory Alloys Patches to Mimic Rolling, Sliding, and Spinning Movements of the Knee

Every year, almost 4 million patients received medical care for knee osteoarthritis. Osteoarthritis involves progressive deterioration or degenerative changes in the cartilage, leading to inflammation and pain as the bones and ligaments are affected. To enhance treatment and surgical outcomes, various studies analyzing the biomechanics of the human skeletal system by fabricating simulated bones, particularly those reflecting the characteristics of patients with knee osteoarthritis, are underway. In this study, we fabricated replicated bones that mirror the bone characteristics of patients with knee osteoarthritis and developed a skeletal model that mimics the actual movement of the knee. To create patient-specific replicated bones, models were extracted from computerized tomography (CT) scans of knee osteoarthritis patients. Utilizing 3D printing technology, we replicated the femur and tibia, which bear the weight of the body and support movement, and manufactured cartilage capable of absorbing and dispersing the impact of knee joint loads using flexible polymers. Furthermore, to implement knee movement in the skeletal model, we developed artificial muscles based on shape memory alloys (SMAs) and used them to mimic the rolling, sliding, and spinning motions of knee flexion. The knee movement was investigated by changing the SMA spring's position, the number of coils, and the applied voltage. Additionally, we developed a knee-joint-mimicking system to analyze the movement of the femur. The proposed artificial-skeletal-model-based knee-joint-mimicking system appears to be applicable for analyzing skeletal models of knee patients and developing surgical simulation equipment for artificial joint replacement surgery.

Development of a three-dimensional muscle-driven lower limb model developed using an improved CFD-FE method

Analysis of the musculoskeletal movements (gait analysis) is needed in many scenarios. The in vivo method has some difficulties. For example, recruiting human subjects for the gait analysis is challenging due to many issues. In addition, when plenty of subjects are required, the follow-up experiments take a long period and the dropout of subjects always occurs. An efficient and reliable in silico simulation platform for gait analysis has been desired for a long time. Therefore, a technique using three-dimensional (3D) muscle modeling to drive the 3D musculoskeletal model was developed and the application of the technique in the simulation of lower limb movements was demonstrated. A finite element model of the lower limb with anatomically high fidelity was developed from the MRI data, where the main muscles, the bones, the subcutaneous tissues, and the skin were reconstructed. To simulate the active behavior of 3D muscles, an active, fiber-reinforced hyperelastic muscle model was developed using the user-defined material (VUMAT) model. Two typical movements, that is, hip abduction and knee lifting, were simulated by activating the responsible muscles. The results show that it is reasonable to use the improved CFD-FE method proposed in the present study to simulate the active contraction of the muscle, and it is feasible to simulate the movements by activating the relevant muscles. The results from the present technique closely match the physiological scenario and thus the technique developed has a great potential to be used in the in silico human simulation platform for many purposes.

Estimation of two wear factors for total hip arthroplasty: A simulation study based on musculoskeletal modelling

BACKGROUND

Primary causes of surgical revision after total hip arthroplasty are polyethylene wear and implant loosening. These factors are particularly related to joint friction and thus patients' physical activity. Assessing implant wear over time according to patients' morphology and physical activity level is key to improve follow-up and patients' quality of life.

METHODS

An approach initially proposed for tibiofemoral prosthetic wear estimation was adapted to compute two wear factors (force-velocity, directional wear intensity) using a musculoskeletal model. It was applied on 17 participants with total hip arthroplasty to compute joint angular velocity, contact force, sliding velocity, and wear factors during common daily living activities.

FINDINGS

Differences were observed between gait, sitting down, and standing up tasks. An incremental increase of both global wear factors (time-integral) was observed during gait from slow to fast speeds (p ≤ 0.01). Interestingly, these two wear factors did not result in same trend for sitting down and standing up tasks. Compared to gait, one cycle of sitting down or standing up tends to induce higher friction-related wear but lower cross-shear-related wear. Depending on the wear factor, significant differences can be found between sitting down and gait at slow speed (p ≤ 0.05), and between sitting down (p ≤ 0.05) or standing up (p ≤ 0.05) and gait at fast speed. Furthermore, depending on the activity, wear can be fostered by joint contact force and/or sliding velocity.

INTERPRETATION

This study demonstrated the potential of wear estimation to highlight activities inducing a higher risk of implant wear after total hip arthroplasty from motion capture data.

Can a Musculoskeletal Model Adapted to Knee Implant Geometry Improve Prediction of 3D Contact Forces and Moments?

Tibiofemoral contact loads are crucial parameters in the onset and progression of osteoarthrosis. While contact loads are frequently estimated from musculoskeletal models, their customization is often limited to scaling musculoskeletal geometry or adapting muscle lines. Moreover, studies have usually focused on superior-inferior contact force without investigating three-dimensional contact loads. Using experimental data from six patients with instrumented total knee arthroplasty (TKA), this study customized a lower limb musculoskeletal model to consider the positioning and the geometry of the implant at knee level. Static optimization was performed to estimate tibiofemoral contact forces and contact moments as well as musculotendinous forces. Predictions from both a generic and a customized model were compared to the instrumented implant measurements. Both models accurately predict superior-inferior (SI) force and abduction-adduction (AA) moment. Notably, the customization improves prediction of medial-lateral (ML) force and flexion-extension (FE) moments. However, there is subject-dependent variability in the prediction of anterior-posterior (AP) force. The customized models presented here predict loads on all joint axes and in most cases improve prediction. Unexpectedly, this improvement was more limited for patients with more rotated implants, suggesting a need for further model adaptations such as muscle wrapping or redefinition of hip and ankle joint centers and axes.

A Two-Degree-of-Freedom Knee Model Predicts Full Three-Dimensional Tibiofemoral and Patellofemoral Joint Motion During Functional Activity

Six kinematic parameters are needed to fully describe three-dimensional (3D) bone motion at a joint. At the knee, the relative movements of the femur and tibia are often represented by a 1-degree-of-freedom (1-DOF) model with a single flexion-extension axis or a 2-DOF model comprising a flexion-extension axis and an internal-external rotation axis. The primary aim of this study was to determine the accuracy with which 1-DOF and 2-DOF models predict the 3D movements of the femur, tibia and patella during daily activities. Each model was created by fitting polynomial functions to 3D tibiofemoral (TF) and patellofemoral (PF) kinematic data recorded from 10 healthy individuals performing 6 functional activities. Model cross-validation analyses showed that the 2-DOF model predicted 3D knee kinematics more accurately than the 1-DOF model. At the TF joint, mean root-mean-square (RMS) errors across all activities and all participants were 3.4°|mm (deg or mm) for the 1-DOF model and 2.4°|mm for the 2-DOF model. At the PF joint, mean RMS errors were 4.0°|mm and 3.9°|mm for the 1-DOF and 2-DOF models, respectively. These results indicate that a 2-DOF model with two rotations as inputs may be used with confidence to predict the full 3D motion of the knee-joint complex.

A subject-specific musculoskeletal model to predict the tibiofemoral contact forces during daily living activities

Accurate prediction of tibiofemoral contact force (TFCF) during daily living activities is significant for understanding the initiation, progression, and treatment of knee osteoarthritis (KOA). However, the diversity of target activities, prediction accuracy, and computational efficiency of the current musculoskeletal simulations need to be further improved. In this study, a subject-specific musculoskeletal model considered the tibiofemoral alignment, medial-lateral contact locations, secondary tibiofemoral and all patellofemoral motions, and knee ligaments was proposed to predict the TFCFs during the five daily activities (normal walking, sit-to-stand, stand-to-sit, stair ascent, and stair descent) in OpenSim software. The standing lower-limbs-full-length radiograph, local radiograph of knee joint, motion capture data, and force plate data of eighteen subjects were acquired as the input data of the musculoskeletal model. The results showed good agreements of TFCFs between the predictions based on our proposed musculoskeletal model and the in-vivo measurements based on instrumented knee implants during the five daily activities (RMSE: 0.16 ∼ 0.31 BW, R²: 0.88 ∼ 0.97, M: -0.11 ∼ -0.02, P: 0.03 ∼ 0.10, and C: 0.04 ∼ 0.14). Additionally, the order of the peak total and lateral TFCFs from low to high was normal walking, stair ascent and stand-to-sit, and stair descent and sit-to-stand (P < 0.05), and the peak medial TFCF was stand-to-sit, sit-to-stand, normal walking, stair ascent and stair descent (P < 0.05). The outcomes of this study are valuable for further understanding the knee biomechanics during daily living activities and providing theoretical guidance for the treatments of KOA.

Functional Resistance Training Differentially Alters Gait Kinetics After Anterior Cruciate Ligament Reconstruction: A Pilot Study

BACKGROUND

Quadriceps weakness is common after anterior cruciate ligament (ACL) reconstruction and can alter gait mechanics. Functional resistance training (FRT) is a novel approach to retraining strength after injury, but it is unclear how it alters gait mechanics. Therefore, we tested how 3 different types of FRT devices: a knee brace resisting extension (unidirectional brace), a knee brace resisting extension and flexion (bidirectional brace), and an elastic band pulling backwards on the ankle (elastic band)-acutely alter gait kinetics in this population.

HYPOTHESIS

The type of FRT device will affect ground-reaction forces (GRFs) during and after the training. Specifically, the uni- and bidirectional braces will increase GRFs when compared with the elastic band.

STUDY DESIGN

Crossover study.

LEVEL OF EVIDENCE

Level 2.

METHODS

A total of 15 individuals with ACL reconstruction received FRT with each device over 3 separate randomized sessions. During training, participants walked on a treadmill while performing a tracking task with visual feedback. Sessions contained 5 training trials (180 seconds each) with rest between. Vertical and anterior-posterior GRFs were assessed on the ACL-reconstructed leg before, during, and after training. Changes in GRFs were compared across devices using 1-dimensional statistical parametric mapping.

RESULTS

Resistance applied via bidirectional brace acutely increased gait kinetics during terminal stance/pre-swing (ie, push-off), while resistance applied via elastic band acutely increased gait kinetics during initial contact/loading (ie, braking). Both braces behaved similarly, but the unidirectional brace was less effective for increasing push-off GRFs.

CONCLUSION

FRT after ACL reconstruction can acutely alter gait kinetics during training. Devices can be applied to selectively alter gait kinetics. However, the long-term effects of FRT after ACL reconstruction with these devices are still unknown.

CLINICAL RELEVANCE

FRT may be applied to alter gait kinetics of the involved limb after ACL reconstruction, depending on the device used.

An Efficient One-Step Moment Balancing Algorithm for Computing Medial and Lateral Knee Compartment Contact Forces

The knee adduction moment is associated with the progression of knee osteoarthritis (OA). The adduction moment reflects the net effect of muscles, passive tissues and bone-on-bone contact forces. Medial compartment OA is more common than lateral and therefore our ability to correctly partition bone-on-bones forces across the medial and lateral compartments is key to understanding mechanical factors associated with the onset and progression of knee OA. In this technical brief we present an efficient one-step moment balancing algorithm linking the sagittal and frontal planes in the determination of musculotendon forces. Novel to the one-step approach is the introduction of a penalty function limiting total compressive force from acting in the lateral compartment when the internal moment is net abduction (i.e., external knee adduction). Medial and lateral knee contact forces were computed using the one-step moment balancing algorithm for 10 subjects walking at a self-selected pace and compared to values determined using a well-established two-step frontal moment balancing approach. Overall, average peak differences in magnitude and timing were small and the ensemble-averaged contact force profiles were similar between methods. The only statistical difference was slightly larger (0.2 BWs) peak medial contact force for the one-step method during the first half of stance, however these differences are small relative to peak values and would not likely alter interpretation of the data. The 1-step moment balancing method is a more efficient methodology for computing medial and lateral knee contact forces that can be used in place of two-step frontal plane moment balancing.

Risk factors for throwing elbow injuries during pitching analyzed by simulation using human musculoskeletal model in youth baseball pitchers

BACKGROUND

Pitching mechanics are believed to be risk factors for throwing elbow injury. Thus, a prospective study of abnormal mechanics in youth baseball players is needed. This study aimed to analyze the ulnar collateral ligament during normal pitching using SIMM (Software for Interactive Musculoskeletal Modeling) for analysis and investigate the risk parameters of throwing elbow injuries in youth baseball players. We hypothesized that excessive ulnar collateral ligament force during pitching would be a risk factor for throwing elbow injuries in this population.

METHODS

In this cohort study, youth baseball pitchers (aged 9-11 years) were instructed to throw a ball into a netted target. Using a SIMM musculoskeletal model, we analyzed the force of the anterior band of the anterior oblique ligament, posterior band of the anterior oblique ligament (AOL_PB), and elbow varus moment during pitching (foot contact to ball release). We calculated the integral of each force of the anterior band of the anterior oblique ligament and AOL_PB during pitching and summarized these data to establish an impulse at the medial epicondyle. Each participant was followed up for 12 months to assess the occurrence of throwing elbow injury.

RESULTS

During the 12-month follow-up period, 18 pitchers (28.1%) reported throwing elbow injuries in the throwing arm. The results of this study showed that the maximum AOL_PB force and the impulse at the medial epicondyle were risk factors for throwing elbow injuries. The maximum AOL_PB force was significantly higher in the throwing elbow injury group than in the uninjured group (59.4 ± 17.8 N vs. 47.1 ± 17.5 N, P = .014). The impulse at the medial epicondyle was also significantly different (11.1 ± 4.0 N ･ s in the throwing elbow injury group vs. 8.3 ± 4.4 N ･ s in the uninjured group, P = .025).

CONCLUSIONS

Increasing the AOL_PB force or the impulse at the medial epicondyle may increase the risk of throwing elbow injuries in youth baseball pitchers. It may be possible to reduce injury risk by focusing on ways to decrease AOL_PB load and cumulative stress on the medial epicondyle throughout the throwing motion while still maintaining high levels of ball velocity.

Knee loading in OA subjects is correlated to flexion and adduction moments and to contact point locations

This study evaluated the association of contact point locations with the knee medial and lateral contact force (F_med, F_lat) alterations in OA and healthy subjects. A musculoskeletal model of the lower limb with subject-specific tibiofemoral contact point trajectories was used to estimate the F_med and F_lat in ten healthy and twelve OA subjects during treadmill gait. Regression analyses were performed to evaluate the correlation of the contact point locations, knee adduction moment (KAM), knee flexion moment (KFM), frontal plane alignment, and gait speed with the F_med and F_lat. Medial contact point locations in the medial–lateral direction showed a poor correlation with the F_med in OA (R² = 0.13, p = 0.01) and healthy (R² = 0.24, p = 0.001) subjects. Anterior–posterior location of the contact points also showed a poor correlation with the F_med of OA subjects (R² = 0.32, p < 0.001). Across all subjects, KAM and KFM remained the best predictors of the F_med and F_lat, respectively (R² between 0.62 and 0.69). Results suggest different mechanisms of contact force distribution in OA joints. The variations in the location of the contact points participate partially to explains the F_med variations in OA subjects together with the KFM and KAM.

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.

Predictions of Anterior Cruciate Ligament Dynamics From Subject-Specific Musculoskeletal Models and Dynamic Biplane Radiography

In vivo knee ligament forces are important to consider for informing rehabilitation or clinical interventions. However, they are difficult to directly measure during functional activities. Musculoskeletal models and simulations have become the primary methods by which to estimate in vivo ligament loading. Previous estimates of anterior cruciate ligament (ACL) forces range widely, suggesting that individualized anatomy may have an impact on these predictions. Using ten subject-specific (SS) lower limb musculoskeletal models, which include individualized musculoskeletal geometry, muscle architecture, and six degree-of-freedom knee joint kinematics from dynamic biplane radiography (DBR), this study provides SS estimates of ACL force (anteromedial-aACL; and posterolateral-pACL bundles) during the full gait cycle of treadmill walking. These forces are compared to estimates from scaled-generic (SG) musculoskeletal models to assess the effect of musculoskeletal knee joint anatomy on predicted forces and the benefit of SS modeling in this context. On average, the SS models demonstrated a double force peak during stance (0.39-0.43 xBW per bundle), while only a single force peak during stance was observed in the SG aACL. No significant differences were observed between continuous SG and SS ACL forces; however, root mean-squared differences between SS and SG predictions ranged from 0.08 xBW to 0.27 xBW, suggesting SG models do not reliably reflect forces predicted by SS models. Force predictions were also found to be highly sensitive to ligament resting length, with ±10% variations resulting in force differences of up to 84%. Overall, this study demonstrates the sensitivity of ACL force predictions to SS anatomy, specifically musculoskeletal joint geometry and ligament resting lengths, as well as the feasibility for generating SS musculoskeletal models for a group of subjects to predict in vivo tissue loading during functional activities.

Accuracy of the tibiofemoral contact forces estimated by a subject-specific musculoskeletal model with fluoroscopy-based contact point trajectories

Accurate estimation of the tibiofemoral contact forces relies on exact kinematics and joint geometry. Subject-specific kinematic constraints representing contact point trajectories derived from fluoroscopic measurements during lunge are introduced in a musculoskeletal model of the lower limb and compared to generic kinematic constraints. The medial, lateral, and total contact forces during gait and squat are validated using the data of four patients with an instrumented prosthesis. The accuracy of the estimated contact forces (both with subject-specific and generic kinematic constraints) remains close to the level reported in the literature. The mean root mean square errors range from 0.32 to 0.52 body weights for gait and from 0.27 to 0.72 body weights for squat. The impact of the subject-specific contact point trajectories is not found substantial or consistent between patients and tasks. Indeed, the kinematics of the total knee prostheses remains close to the kinematics of a hinge joint and the contact point locations remain generally centred at 20 mm from the tibia centreline (close to the constant value defined in the generic constraints). The contact point trajectories are also suspected to differ between tasks (lunge vs. gait and squat). While the contact point trajectories have been reported to be sensitive model parameters, no clear improvement of the contact force accuracy is demonstrated on patients with instrumented prosthesis. The introduction (as kinematic constraints) of fluoroscopy-based contact point trajectories may be considered in cases where these trajectories are significantly altered, as reported for osteoarthritis patients.

Strain Assessment of Deep Fascia of the Thigh During Leg Movement: An in situ Study

Fascia is a fibrous connective tissue present all over the body. At the lower limb level, the deep fascia that is overlying muscles of the outer thigh and sheathing them (fascia lata) is involved in various pathologies. However, the understanding and quantification of the mechanisms involved in these sheathing effects are still unclear. The aim of this study is to observe and quantify the strain field of the fascia lata, including the iliotibial tract (ITT), during a passive movement of the knee. Three fresh postmortem human subjects were studied. To measure hip and knee angles during knee flexion-extension, passive movements from 0° to around 120° were recorded with a motion analysis system and strain fields of the fascia were acquired using digital image correlation. Strains were computed for three areas of the fascia lata: anterior fascia, lateral fascia, and ITT. Mean principal strains showed different strain mechanisms depending on location on the fascia and knee angle. For the ITT, two strain mechanisms were observed depending on knee movement: compression is observed when the knee is extended relative to the reference position of 47°, however, tension and pure shear can be observed when the knee is flexed. For the anterior and lateral fascia, in most cases, minor strain is higher than major strain in absolute value, suggesting high tissue compression probably due to microstructural fiber rearrangements. This in situ study is the first attempt to quantify the superficial strain field of fascia lata during passive leg movement. The study presents some limitations but provides a step in understanding strain mechanism of the fascia lata during passive knee movement.

An Anatomical-Based Subject-Specific Model of In-Vivo Knee Joint 3D Kinematics From Medical Imaging

Biomechanical models of the knee joint allow the development of accurate procedures as well as novel devices to restore the joint natural motion. They are also used within musculoskeletal models to perform clinical gait analysis on patients. Among relevant knee models in the literature, the anatomy-based spatial parallel mechanisms represent the joint motion using rigid links for the ligaments’ isometric fibres and point contacts for the articular surfaces. To customize analyses, therapies and devices, there is the need to define subject-specific models, but relevant procedures and their accuracy are still questioned. A procedure is here proposed and validated to define a customized knee model based on a spatial parallel mechanism. Computed tomography, magnetic resonance and 3D-video-fluoroscopy were performed on a healthy volunteer to define the personalized model geometry. The model was then validated by comparing the measured and the replicated joint motion. The model showed mean absolute difference and standard deviations in translations and rotations, respectively of 0.98 ± 0.40 mm and 0.68 ± 0.29 ° for the tibia–femur motion, and of 0.77 ± 0.15 mm and 2.09 ± 0.69 ° for the patella–femur motion. These results show that accurate personalized spatial models of knee kinematics can be obtained from in-vivo imaging.

Evaluation of predicted patellofemoral joint kinematics with a moving-axis joint model

The main objectives of this study were to expand the moving-axis joint model concept to the patellofemoral joint and evaluate the patellar motion against experimental patellofemoral kinematics. The experimental data was obtained through 2D-to-3D bone reconstruction of EOS images and segmented MRI data utilizing an iterative closest point optimization technique. Six knee model variations were developed using the AnyBody Modeling System and subject-specific bone geometries. These models consisted of various combinations of tibiofemoral (hinge, moving-axis, and interpolated) and patellofemoral (hinge and moving-axis) joint types. The newly introduced interpolated tibiofemoral joint is calibrated from the five EOS quasi-static lunge positions. The patellofemoral axis of the hinge model was defined by performing surface fits to the patellofemoral contact area; and the moving-axis model was defined based upon the position of the patellofemoral joint at 0° and 90° tibiofemoral-flexion. In between these angles, the patellofemoral axis moved linearly as a function of tibiofemoral-flexion, while outside these angles, the axis remained fixed. When using a moving-axis tibiofemoral joint, a hinge patellofemoral joint offers (-5.12 ± 1.23 mm, 5.81 ± 0.97 mm, 14.98 ± 2.30°, -4.35 ± 1.95°) mean differences (compared to EOS) while a moving-axis patellofemoral model provides (-2.69 ± 1.04 mm, 1.13 ± 0.80 mm, 12.63 ± 2.03°, 1.74 ± 1.46°) in terms of lateral-shift, superior translation, patellofemoral-flexion, and patellar-rotation, respectively. Furthermore, the model predictive capabilities increased as a direct result of adding more calibrated positions to the tibiofemoral model (hinge-1, moving-axis-2, and interpolated-5). Overall, a novel subject-specific moving-axis patellofemoral model has been established; that produces realistic patellar motion and is computationally fast enough for clinical applications.

Predictive Simulations of Neuromuscular Coordination and Joint-Contact Loading in Human Gait

We implemented direct collocation on a full-body neuromusculoskeletal model to calculate muscle forces, ground reaction forces and knee contact loading simultaneously for one cycle of human gait. A data-tracking collocation problem was solved for walking at the normal speed to establish the practicality of incorporating a 3D model of articular contact and a model of foot-ground interaction explicitly in a dynamic optimization simulation. The data-tracking solution then was used as an initial guess to solve predictive collocation problems, where novel patterns of movement were generated for walking at slow and fast speeds, independent of experimental data. The data-tracking solutions accurately reproduced joint motion, ground forces and knee contact loads measured for two total knee arthroplasty patients walking at their preferred speeds. RMS errors in joint kinematics were < 2.0° for rotations and < 0.3 cm for translations while errors in the model-computed ground-reaction and knee-contact forces were < 0.07 BW and < 0.4 BW, respectively. The predictive solutions were also consistent with joint kinematics, ground forces, knee contact loads and muscle activation patterns measured for slow and fast walking. The results demonstrate the feasibility of performing computationally-efficient, predictive, dynamic optimization simulations of movement using full-body, muscle-actuated models with realistic representations of joint function.

Musculoskeletal multibody dynamics simulation of the contact mechanics and kinematics of a natural knee joint during a walking cycle

Detailed knowledge of the in vivo loading and kinematics in the knee joint is essential to understand its normal functions and the aetiology of osteoarthritis. Computer models provide a viable non-invasive solution for estimating joint loading and kinematics during different physiological activities. However, the joint loading and kinematics of the tibiofemoral and patellofemoral joints during a gait cycle were not typically investigated concurrently in previous computational simulations. In this study, a natural knee architecture was incorporated into a lower extremity musculoskeletal multibody dynamics model based on a force-dependent kinematics approach to investigate the contact mechanics and kinematics of a natural knee joint during a walking cycle. Specifically, the contact forces between the femoral/tibial articular cartilages and menisci and between the femoral and tibial/patellar articular cartilages were quantified. The contact forces and kinematics of the tibiofemoral and patellofemoral joints and the muscle activations and ligament forces were predicted simultaneously with a reasonable level of accuracy. The developed musculoskeletal multibody dynamics model with a natural knee architecture can serve as a potential platform for assisting clinical decision-making and postoperative rehabilitation planning.

Biomechanical comparison of fixation systems in posterior wall fracture of acetabular by finite element analysis

BACKGROUND

The use of reconstruction plates and lag screws has been recommended for fractures to the posterior wall of the acetabulum. However, little information about the rigidity of recommended forms of fracture fixation is available. This study aimed to evaluate the biomechanical difference among the fixation systems.

METHODS

A posterior wall fracture, which is represented by softer elements with lower elastic modulus, was created along an arc of 40-90° about the acetabular rim. Three different fixation systems: screws alone, reconstruction plate with screws, reconstruction plate with T-shaped plates were used to fix the posterior wall fractures to the acetabulum.

RESULTS

All three fixation system can be used to obtain good functional outcomes. The reconstruction plate with T-shaped plates was beneficial to increasing the effective stiffness, decreasing the stress concentration and enhancing the rigidity of fracture fixation. So this fixation system served an ideal result in the analysis.

CONCLUSION

Theoretically, the reconstruction plate with T-shaped plates system may reduce many of the risks and limitations compared to the other fixation systems. This fixation system may result in a clinical benefit.

A computational framework for simultaneous estimation of muscle and joint contact forces and body motion using optimization and surrogate modeling

Concurrent estimation of muscle activations, joint contact forces, and joint kinematics by means of gradient-based optimization of musculoskeletal models is hindered by computationally expensive and non-smooth joint contact and muscle wrapping algorithms. We present a framework that simultaneously speeds up computation and removes sources of non-smoothness from muscle force optimizations using a combination of parallelization and surrogate modeling, with special emphasis on a novel method for modeling joint contact as a surrogate model of a static analysis. The approach allows one to efficiently introduce elastic joint contact models within static and dynamic optimizations of human motion. We demonstrate the approach by performing two optimizations, one static and one dynamic, using a pelvis-leg musculoskeletal model undergoing a gait cycle. We observed convergence on the order of seconds for a static optimization time frame and on the order of minutes for an entire dynamic optimization. The presented framework may facilitate model-based efforts to predict how planned surgical or rehabilitation interventions will affect post-treatment joint and muscle function.

Knee medial and lateral contact forces in a musculoskeletal model with subject-specific contact point trajectories

Contact point (CP) trajectory is a crucial parameter in estimating medial/lateral tibio-femoral contact forces from the musculoskeletal (MSK) models. The objective of the present study was to develop a method to incorporate the subject-specific CP trajectories into the MSK model. Ten healthy subjects performed 45 s treadmill gait trials. The subject-specific CP trajectories were constructed on the tibia and femur as a function of extension-flexion using low-dose bi-plane X-ray images during a quasi-static squat. At each extension-flexion position, the tibia and femur CPs were superimposed in the three directions on the medial side, and in the anterior-posterior and proximal-distal directions on the lateral side to form the five kinematic constraints of the knee joint. The Lagrange multipliers associated to these constraints directly yielded the medial/lateral contact forces. The results from the personalized CP trajectory model were compared against the linear CP trajectory and sphere-on-plane CP trajectory models which were adapted from the commonly used MSK models. Changing the CP trajectory had a remarkable impact on the knee kinematics and changed the medial and lateral contact forces by 1.03 BW and 0.65 BW respectively, in certain subjects. The direction and magnitude of the medial/lateral contact force were highly variable among the subjects and the medial-lateral shift of the CPs alone could not determine the increase/decrease pattern of the contact forces. The suggested kinematic constraints are adaptable to the CP trajectories derived from a variety of joint models and those experimentally measured from the 3D imaging techniques.

Three-dimensional multifractal analysis of trabecular bone under clinical computed tomography

PURPOSE

An adequate understanding of bone structural properties is critical for predicting fragility conditions caused by diseases such as osteoporosis, and in gauging the success of fracture prevention treatments. In this work we aim to develop multiresolution image analysis techniques to extrapolate high-resolution images predictive power to images taken in clinical conditions.

METHODS

We performed multifractal analysis (MFA) on a set of 17 ex vivo human vertebrae clinical CT scans. The vertebrae failure loads (FFailure) were experimentally measured. We combined bone mineral density (BMD) with different multifractal dimensions, and BMD with multiresolution statistics (e.g., skewness, kurtosis) of MFA curves, to obtain linear models to predict FFailure. Furthermore we obtained short- and long-term precisions from simulated in vivo scans, using a clinical CT scanner. Ground-truth data - high-resolution images - were obtained with a High-Resolution Peripheral Quantitative Computed Tomography (HRpQCT) scanner.

RESULTS

At the same level of detail, BMD combined with traditional multifractal descriptors (Lipschitz-Hölder exponents), and BMD with monofractal features showed similar prediction powers in predicting FFailure (87%, adj. R² ). However, at different levels of details, the prediction power of BMD with multifractal features raises to 92% (adj. R2) of FFailure. Our main finding is that a simpler but slightly less accurate model, combining BMD and the skewness of the resulting multifractal curves, predicts 90% (adj. R2) of FFailure.

CONCLUSIONS

Compared to monofractal and standard bone measures, multifractal analysis captured key insights in the conditions leading to FFailure. Instead of raw multifractal descriptors, the statistics of multifractal curves can be used in several other contexts, facilitating further research.

Tenosuspension of the Reflected Head of the Rectus Femoris in Hip Arthroscopy: Description of a Portal and a Surgical Maneuver

Arthroscopy is a surgical technique useful in the treatment of intra- and extra-articular hip pathologies, including femoroacetabular impingement (FAI). In the arthroscopic treatment of pincer-type FAI, tendon interposition of the reflected head of the rectus femoris frequently hinders acetabular edge resection in the anterior-superior region (acetabuloplasty) and labrum repositioning, even causing tendon injury during the surgical procedure. Many surgeons do not give importance to the reflected head of the rectus femoris, and during the procedure they try to avoid it, giving 30° of flexion to the hip. However, this does not always prevent tendon injury caused by the instruments. Some surgeons even resect part of the tendon to improve visualization or to use it as a labral graft. The objective of this work is to describe a reproducible surgical maneuver that allows to perform the anterior-superior acetabular edge resection and chondrolabral union repair with less damage to the reflected head of the rectus femoris.

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.

Modeling the Human Tibiofemoral Joint Using Ex Vivo Determined Compliance Matrices

Several approaches have been used to devise a model of the human tibiofemoral joint for embedment in lower limb musculoskeletal models. However, no study has considered the use of cadaveric 6 × 6 compliance (or stiffness) matrices to model the tibiofemoral joint under normal or pathological conditions. The aim of this paper is to present a method to determine the compliance matrix of an ex vivo tibiofemoral joint for any given equilibrium pose. Experiments were carried out on a single ex vivo knee, first intact and, then, with the anterior cruciate ligament (ACL) transected. Controlled linear and angular displacements were imposed in single degree-of-freedom (DoF) tests to the specimen, and the resulting forces and moments were measured using an instrumented robotic arm. This was done starting from seven equilibrium poses characterized by the following flexion angles: 0 deg, 15 deg, 30 deg, 45 deg, 60 deg, 75 deg, and 90 deg. A compliance matrix for each of the selected equilibrium poses and for both the intact and ACL-deficient specimen was calculated. The matrix, embedding the experimental load-displacement relationship of the examined DoFs, was calculated using a linear least squares inversion based on a QR decomposition, assuming symmetric and positive-defined matrices. Single compliance matrix terms were in agreement with the literature. Results showed an overall increase of the compliance matrix terms due to the ACL transection (2.6 ratio for rotational terms at full extension) confirming its role in the joint stabilization. Validation experiments were carried out by performing a Lachman test (the tibia is pulled forward) under load control on both the intact and ACL-deficient knee and assessing the difference (error) between measured linear and angular displacements and those estimated using the appropriate compliance matrix. This error increased nonlinearly with respect to the values of the load. In particular, when an incremental posterior-anterior force up to 6 N was applied to the tibia of the intact specimen, the errors on the estimated linear and angular displacements were up to 0.6 mm and 1.5 deg, while for a force up to 18 N, the errors were 1.5 mm and 10.5 deg, respectively. In conclusion, the method used in this study may be a viable alternative to characterize the tibiofemoral load-dependent behavior in several applications.

Influence of the Level of Muscular Redundancy on the Validity of a Musculoskeletal Model

While recent literature has clearly demonstrated that an extensive personalization of the musculoskeletal models was necessary to reach high accuracy, several components of the generic models may be further investigated before defining subject-specific parameters. Among others, the choice in muscular geometry and thus the level of muscular redundancy in the model may have a noticeable influence on the predicted musculotendon and joint contact forces. In this context, the aim of this study was to investigate if the level of muscular redundancy can contribute or not to reduce inaccuracies in tibiofemoral contact forces predictions. For that, the dataset disseminated through the Sixth Grand Challenge Competition to Predict In Vivo Knee Loads was applied to a versatile 3D lower limb musculoskeletal model in which two muscular geometries (i.e., two different levels of muscular redundancy) were implemented. This dataset provides tibiofemoral implant measurements for both medial and lateral compartments and thus allows evaluation of the validity of the model predictions. The results suggest that an increase of the level of muscular redundancy corresponds to a better accuracy of total tibiofemoral contact force whatever the gait pattern investigated. However, the medial and lateral contact forces ratio and accuracy were not necessarily improved when increasing the level of muscular redundancy and may thus be attributed to other parameters such as the location of contact points. To conclude, the muscular geometry, among other components of the generic model, has a noticeable impact on joint contact forces predictions and may thus be correctly chosen even before trying to personalize the model.

Dependence of Muscle Moment Arms on In Vivo Three-Dimensional Kinematics of the Knee

Quantification of muscle moment arms is important for clinical evaluation of muscle pathology and treatment, and for estimating muscle and joint forces in musculoskeletal models. Moment arms estimated with musculoskeletal models often assume a default motion of the knee derived from measurements of passive cadaveric flexion. However, knee kinematics are unique to each person and activity. The objective of this study was to estimate moment arms of the knee muscles with in vivo subject- and activity-specific kinematics from seven healthy subjects performing seated knee extension and single-leg lunge to show changes between subjects and activities. 3D knee motion was measured with a high-speed stereo-radiography system. Moment arms of ten muscles were estimated in OpenSim by replacing the default knee motion with in vivo measurements. Estimated inter-subject moment arm variability was similar to previously reported in vitro measurements. RMS deviations up to 9.0 mm (35.2% of peak value) were observed between moment arms estimated with subject-specific knee extension and passive cadaveric motion. The degrees of freedom that most impacted inter-activity differences were superior/inferior and anterior/posterior translations. Musculoskeletal simulations used to estimate in vivo muscle forces and joint loads may provide significantly different results when subject- and activity-specific kinematics are implemented.

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

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

In Vivo Knee Contact Force Prediction Using Patient-Specific Musculoskeletal Geometry in a Segment-Based Computational Model

Segment-based musculoskeletal models allow the prediction of muscle, ligament, and joint forces without making assumptions regarding joint degrees-of-freedom (DOF). The dataset published for the “Grand Challenge Competition to Predict in vivo Knee Loads” provides directly measured tibiofemoral contact forces for activities of daily living (ADL). For the Sixth Grand Challenge Competition to Predict in vivo Knee Loads, blinded results for “smooth” and “bouncy” gait trials were predicted using a customized patient-specific musculoskeletal model. For an unblinded comparison, the following modifications were made to improve the predictions: further customizations, including modifications to the knee center of rotation; reductions to the maximum allowable muscle forces to represent known loss of strength in knee arthroplasty patients; and a kinematic constraint to the hip joint to address the sensitivity of the segment-based approach to motion tracking artifact. For validation, the improved model was applied to normal gait, squat, and sit-to-stand for three subjects. Comparisons of the predictions with measured contact forces showed that segment-based musculoskeletal models using patient-specific input data can estimate tibiofemoral contact forces with root mean square errors (RMSEs) of 0.48–0.65 times body weight (BW) for normal gait trials. Comparisons between measured and predicted tibiofemoral contact forces yielded an average coefficient of determination of 0.81 and RMSEs of 0.46–1.01 times BW for squatting and 0.70–0.99 times BW for sit-to-stand tasks. This is comparable to the best validations in the literature using alternative models.

A Hertzian Integrated Contact Model of the Total Knee Replacement Implant for the Estimation of Joint Contact Forces

The prediction of lower limb muscle and contact forces may provide useful knowledge to assist the clinicians in the diagnosis as well as in the development of appropriate treatment for musculoskeletal disorders. Research studies have commonly estimated joint contact forces using model-based muscle force estimation due to the lack of a reliable contact model and material properties. The objective of this present study was to develop a Hertzian integrated contact model. Then, in vivo elastic properties of the Total Knee Replacement (TKR) implant were identified using in vivo contact forces leading to providing reliable material properties for modeling purposes. First, a patient specific rigid musculoskeletal model was built. Second, a STL-based implant model was designed to compute the contact area evolutions during gait motions. Finally, a Hertzian integrated contact model was defined for the in vivo identification of elastic properties (Young’s modulus and Poisson coefficient) of the instrumented TKR implant. Our study showed a potential use of a new approach to predict the contact forces without knowledge of muscle forces. Thus, the outcomes may lead to accurate and reliable prediction of human joint contact forces for new case study.

An Open-Source Toolbox for Surrogate Modeling of Joint Contact Mechanics

GOAL

Incorporation of elastic joint contact models into simulations of human movement could facilitate studying the interactions between muscles, ligaments, and bones. Unfortunately, elastic joint contact models are often too expensive computationally to be used within iterative simulation frameworks. This limitation can be overcome by using fast and accurate surrogate contact models that fit or interpolate input-output data sampled from existing elastic contact models. However, construction of surrogate contact models remains an arduous task. The aim of this paper is to introduce an open-source program called Surrogate Contact Modeling Toolbox (SCMT) that facilitates surrogate contact model creation, evaluation, and use.

METHODS

SCMT interacts with the third-party software FEBio to perform elastic contact analyses of finite-element models and uses MATLAB to train neural networks that fit the input-output contact data. SCMT features sample point generation for multiple domains, automated sampling, sample point filtering, and surrogate model training and testing.

RESULTS

An overview of the software is presented along with two example applications. The first example demonstrates creation of surrogate contact models of artificial tibiofemoral and patellofemoral joints and evaluates their computational speed and accuracy, while the second demonstrates the use of surrogate contact models in a forward dynamic simulation of an open-chain leg extension-flexion motion.

CONCLUSION

SCMT facilitates the creation of computationally fast and accurate surrogate contact models. Additionally, it serves as a bridge between FEBio and OpenSim musculoskeletal modeling software.

SIGNIFICANCE

Researchers may now create and deploy surrogate models of elastic joint contact with minimal effort.

Insertion of a pressure sensing arrayminimally affects hindfoot bone kinematics

BACKGROUND

Understanding the development of ankle osteoarthritis (OA) is of high importance and interest; however its causality is poorly understood and several links to joint loading conditions have been made. One way of quantifying joint loading conditions is by measuring the intra-articular pressure distribution during gait simulations performed by in-vitro experimental set-ups. However the effect of inserting a pressure sensing array in the ankle joint could potentially disturb the proper kinematics and therefore the loading conditions.

METHODS

In this study, we performed in-vitro gait simulations in 7 cadaveric feet, before and after inserting a pressure sensing array and quantified the effect on the joints range of motion (ROM). The gait was simulated with a stance phase duration of one second using a custom build cadaveric gait simulator (CGS).

RESULTS

The results show a limited effect in the ROM for all the joints of the hind foot, not exceeding the variability observed in specimens without a sensor. However, no consistent direction (increase/decrease) can be observed.

CONCLUSION

The results suggest that even though the effect of inserting a pressure sensing array is minimal, it needs to be evaluated against the demands/requirements of the application.

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

Traditional approaches to the biomechanical analysis of movement are joint-based; that is the mechanics of the body are described in terms of the forces and moments acting at the joints, and that muscular forces are considered to create moments about the joints. We have recently shown that segment-based approaches, where the mechanics of the body are described by considering the effect of the muscle, ligament and joint contact forces on the segments themselves, can also prove insightful. We have also previously described a simultaneous, optimization-based, musculoskeletal model of the lower limb. However, this prior model incorporates both joint- and segment-based assumptions. The purpose of this study was therefore to develop an entirely segment-based model of the lower limb and to compare its performance to our previous work. The segment-based model was used to estimate the muscle forces found during vertical jumping, which were in turn compared with the muscular activations that have been found in vertical jumping, by using a Geers' metric to quantify the magnitude and phase errors. The segment-based model was shown to have a similar ability to estimate muscle forces as a model based upon our previous work. In the future, we will evaluate the ability of the segment-based model to be used to provide results with clinical relevance, and compare its performance to joint-based approaches. The segment-based model described in this article is publicly available as a GUI-based Matlab® application and in the original source code (at www.msksoftware.org.uk).

Malrotated tibial component increases medial collateral ligament tension in total knee arthroplasty

Malrotation of the tibial component can lead to complications after total knee arthroplasty (TKA). Despite reports of internal rotation being associated with more severe pain or stiffness than external rotation, the biomechanical reasons remain largely unknown. We used a computer simulation model and evaluated traction forces in the lateral collateral ligament (LCL) and medial collateral ligament (MCL) with a malrotated tibial component during squatting. We also examined tibiofemoral and patellofemoral contact forces and stresses under similar conditions. A dynamic musculoskeletal knee model was simulated in three different constrained tibial geometries with a prototype component. The testing conditions were changed between 15° external and 15° internal rotation of the tibial component. With internal rotation of the tibial component, the MCL force increased progressively; the LCL force also increased, but only up to less than half of the MCL force values. A higher degree of constraint of the tibial component was associated with greater femoral rotational movement and higher MCL forces. The tibiofemoral and patellofemoral contact forces were not influenced by malrotation of the tibial component, but the contact stresses increased because of decreased contact area. This altered loading condition could cause patient complaints and polyethylene problems after TKA.

Multi-objective optimisation for musculoskeletal modelling: application to a planar elbow model

One of the open issues in musculoskeletal modelling remains the choice of the objective function that is used to solve the muscular redundancy problem. Some authors have recently proposed to introduce joint reaction forces in the objective function, and the question of the weights associated with musculo-tendon forces and joint reaction forces arose. This question typically deals with a multi-objective optimisation problem. The aim of this study is to illustrate, on a planar elbow model, the ensemble of optimal solutions (i.e. Pareto front) and the solution of a global objective method that represent different compromises between musculo-tendon forces, joint compression force, and joint shear force. The solutions of the global objective method, based either on the minimisation of the sum of the squared musculo-tendon forces alone or on the minimisation of the squared joint compression force and shear force together, are in the same range. Minimising either the squared joint compression force or shear force alone leads to extreme force values. The exploration of the compromises between these forces illustrates the existence of major interactions between the muscular and joint structures. Indeed, the joint reaction forces relate to the projection of the sum of the musculo-tendon forces. An illustration of these interactions, due to the projection relation, is that the Pareto front is not a large surface, like in a typical three-objective optimisation, but almost a curve. These interactions, and the possibility to take them into account by a multi-objective optimisation, seem essential for the application of musculoskeletal modelling to joint pathologies.

Prediction of in vivo joint mechanics of an artificial knee implant using rigid multi-body dynamics with elastic contacts

Lower extremity musculoskeletal computational models play an important role in predicting joint forces and muscle activation simultaneously and are valuable for investigating functional outcomes of the implants. However, current computational musculoskeletal models of total knee replacement rarely consider the bearing surface geometry of the implant. Therefore, these models lack detailed information about the contact loading and joint motion which are important factors for evaluating clinical performances. This study extended a rigid multi-body dynamics simulation of a lower extremity musculoskeletal model to incorporate an artificial knee joint, based upon a novel force-dependent kinematics method, and to characterize the in vivo joint contact mechanics during gait. The developed musculoskeletal total knee replacement model integrated the rigid skeleton multi-body dynamics and the flexible contact mechanics of the tibiofemoral and patellofemoral joints. The predicted contact forces and muscle activations are compared against those in vivo measurements obtained from a single patient with good agreements for the medial contact force (root-mean-square error = 215 N, ρ = 0.96) and lateral contact force (root-mean-square error = 179 N, ρ = 0.75). Moreover, the developed model also predicted the motion of the tibiofemoral joint in all degrees of freedom. This new model provides an important step toward the development of a realistic dynamic musculoskeletal total knee replacement model to predict in vivo knee joint motion and loading simultaneously. This could offer a better opportunity to establish a robust virtual modeling platform for future pre-clinical assessment of knee prosthesis designs, surgical procedures and post-operation rehabilitation.