3D active-passive response of human knee joint in gait is markedly altered when simulated as a planar 2D joint

Computational biomechanics of human knee joint in stair ascent: Muscle-ligament-contact forces and comparison with level walking

About a third of knee joint disorders originate from the patellofemoral (PF) site that makes stair ascent a difficult activity for patients. A detailed finite element model of the knee joint is coupled to a lower extremity musculoskeletal model to simulate the stance phase of stair ascent. It is driven by the mean of measurements on the hip-knee-ankle moments-angles as well as ground reaction forces reported in healthy individuals. Predicted muscle activities compare well to the recorded electromyography data. Peak forces in quadriceps (3.87 BW, body weight, at 20% instance in our 607 N subject), medial hamstrings (0.77 BW at 20%), and gastrocnemii (1.21 BW at 80%) are estimated. Due to much greater flexion angles-moments in the first half of stance, large PF contact forces (peak of 3.1 BW at 20% stance) and stresses (peak of 4.83 MPa at 20% stance) are estimated that exceed their peaks in level walking by fourfold and twofold, respectively. Compared with level walking, ACL forces diminish in the first half of stance but substantially increase later in the second half (peak of 0.76 BW at 75% stance). Under nearly similar contact forces at 20% of stance, the contact stress on the tibiofemoral (TF) medial plateau reaches a peak (9.68 MPa) twice that on the PF joint suggesting the vulnerability of both joints. Compared with walking, stair ascent increases peak ACL force and both peak TF and PF contact stresses. Reductions in the knee flexion moment and/or angle appear as a viable strategy to mitigate internal loads and pain.

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.

Towards Tailored Rehabilitation by Implementation of a Novel Musculoskeletal Finite Element Analysis Pipeline

Tissue-level mechanics (e.g., stress and strain) are important factors governing tissue remodeling and development of knee osteoarthritis (KOA), and hence, the success of physical rehabilitation. To date, no clinically feasible analysis toolbox has been introduced and used to inform clinical decision making with subject-specific in-depth joint mechanics of different activities. Herein, we utilized a rapid state-of-the-art electromyography-assisted musculoskeletal finite element analysis toolbox with fibril-reinforced poro(visco)elastic cartilages and menisci to investigate knee mechanics in different activities. Tissue mechanical responses, believed to govern collagen damage, cell death, and fixed charge density loss of proteoglycans, were characterized within 15 patients with KOA while various daily activities and rehabilitation exercises were performed. Results showed more inter-participant variation in joint mechanics during rehabilitation exercises compared to daily activities. Accordingly, the devised workflow may be used for designing subject-specific rehabilitation protocols. Further, results showed the potential to tailor rehabilitation exercises, or assess capacity for daily activity modifications, to optimally load knee tissue, especially when mechanically-induced cartilage degeneration and adaptation are of interest.

Biomechanical Characteristics of the Knee Joint during Gait in Obese versus Normal Subjects

Knee osteoarthritis (OA) is a growing source of pain and disability. Obesity is the most important avoidable risk factor underlying knee OA. The processes by which obesity impacts osteoarthritis are of tremendous interest to osteoarthritis researchers and physicians, where the joint mechanical load is one of the pathways generally thought to cause or intensify the disease process. In the current work, we developed a hybrid framework that simultaneously incorporates a detailed finite element model of the knee joint within a musculoskeletal model to compute lower extremity muscle forces and knee joint stresses in normal-weight (N) and obese (OB) subjects during the stance phase gait. This model accounts for the synergy between the active musculature and passive structures. In comparing OB subjects and normal ones, forces significantly increased in all muscle groups at most instances of stance. Mainly, much higher activation was computed with lateral hamstrings and medial gastrocnemius. Cartilage contact average pressure was mostly supported by the medial plateau and increased by 22%, with a larger portion of the load transmitted via menisci. This medial compartment experienced larger relative movement and cartilage stresses in the normal subjects and continued to do so with a higher level in the obese subjects. Finally, the developed bioengineering frame and the examined parameters during this investigation might be useful clinically in evaluating the initiation and propagation of knee OA.

Alterations in the Functional Knee Alignment Are Not an Effective Strategy to Modify the Mediolateral Distribution of Knee Forces During Closed Kinetic Chain Exercises

Pain felt while performing rehabilitation exercises could be a reason for the low adherence of knee osteoarthritis patients to physical rehabilitation. Reducing compressive forces on the most affected knee regions may help to mitigate the pain. Knee frontal plane positioning with respect to pelvis and foot (functional knee alignment) has been shown to modify the mediolateral distribution of the tibiofemoral joint contact force in walking. Hence, different functional knee alignments could be potentially used to modify joint loading during rehabilitation exercises. The aim was to understand whether utilizing different alignments is an effective strategy to unload specific knee areas while performing rehabilitation exercises. Eight healthy volunteers performed 5 exercises with neutral, medial, and lateral knee alignment. A musculoskeletal model was modified for improved prediction of tibiofemoral contact forces and used to evaluate knee joint kinematics, moments, and contact forces. Functional knee alignment had only a small and inconsistent effect on the mediolateral distribution joint contact force. Moreover, the magnitude of tibiofemoral and patellofemoral contact forces, knee moments, and measured muscle activities was not significantly affected by the alignment. Our results suggest that altering the functional knee alignment is not an effective strategy to unload specific knee regions in physical rehabilitation.

A Computationally Efficient Lower Limb Finite Element Musculoskeletal Framework Directly Driven Solely by Inertial Measurement Unit Sensors

Finite element musculoskeletal (FEMS) approaches using concurrent musculoskeletal and finite element models driven by motion data such as marker-based motion trajectory can provide insight into the interactions between the knee joint secondary kinematics, contact mechanics, and muscle forces in subject-specific biomechanical investigations. However, these data-driven FEMS systems have two major disadvantages that make them challenging to apply in clinical environments: they are computationally expensive and they require expensive and inconvenient equipment for data acquisition. In this study, we developed an FEMS model of the lower limb driven solely by inertial measurement unit sensors that includes the tissue geometries of the entire knee joint and combines muscle modeling and elastic foundation theory-based contact analysis of knee into a single framework. The model requires only the angular velocities and accelerations measured by the sensors as input, and the target outputs (knee contact mechanics, secondary kinematics, and muscle forces) are predicted from the convergence results of iterative calculations of muscle force optimization and knee contact mechanics. To evaluate its accuracy, the model was compared with in vivo experimental data during gait. The maximum contact pressure (12.6 MPa) in the rigid body contact analysis occurred on the medial side of the cartilage at the maximum loading response. The proposed computationally efficient framework drastically reduced the computational time (97.5% reduction) in comparison with the conventional deformable finite element analysis. The developed framework combines measurement convenience and computational efficiency and shows promise for clinical applications.

A novel coupled musculoskeletal finite element model of the spine - Critical evaluation of trunk models in some tasks

Spine musculoskeletal (MS) models make simplifying assumptions on the intervertebral joint degrees-of-freedom (rotational and/or translational), representation (spherical or beam-like joints), and properties (linear or nonlinear). They also generally neglect the realistic structure of the joints with disc nuclei/annuli, facets, and ligaments. We aim to develop a novel MS model where trunk muscles are incorporated into a detailed finite element (FE) model of the ligamentous T12-S1 spine thus constructing a gold standard coupled MS-FE model. Model predictions are compared under some tasks with those of our earlier spherical joints, beam joints, and hybrid (uncoupled) MS-FE models. The coupled model predicted L4-L5 intradiscal pressures (R² ≅ 0.97, RMSE ≅ 0.27 MPa) and L1-S1 centers of rotation (CoRs) in agreement to in vivo data. Differences in model predictions grew at larger trunk flexion angles; at the peak (80°) flexion the coupled model predicted, compared to the hybrid model, much smaller global/local muscle forces (~38%), segmental (~44%) and disc (~22%) compression forces but larger segmental (~9%) and disc (~17%) shear loads, ligament forces at the lower lumbar levels (by up to 57%) and facet forces at all levels. The spherical/beam joints models predicted much greater muscle forces and segmental loads under larger flexion angles. Unlike the spherical joints model with fixed CoRs, the beam joints model predicted CoRs closer (RMSE = 2.3 mm in flexion tasks) to those of the coupled model. The coupled model offers a great potential for future studies towards improvement of surgical techniques, management of musculoskeletal injuries and subject-specific simulations.

Increased Q-factor increases medial compartment knee joint contact force during cycling

As Q-Factor (QF: inter-pedal distance) is increased, the internal knee abduction moment (KAbM) also increases, however it is unknown if this increased KAbM is associated with increased medial compartment knee joint contact force in cycling. In the absence of in vivo measurement, musculoskeletal modeling simulations may provide a viable option for estimating knee joint contact forces in cycling. The primary purpose of this study was to investigate the effect of increasing QF on knee joint total (TCF), and medial (MCF) compartment contact force during ergometer cycling. The secondary purpose was to evaluate whether KAbM and knee extension moment are accurate predictors of MCF in cycling. Musculoskeletal simulations were performed to estimate TCF and MCF for sixteen participants cycling at an original QF (150 mm), and wide QF (276 mm), at 80 W and 80 rotations per minute. Paired samples t-tests were used to detect differences between QF conditions. MCF increased significantly, however, TCF did not change at wide QF. Peak knee extensor muscle force did not change at wide QF. Peak knee flexor muscle force was significantly reduced with wide QF. Regression analyses showed KAbM and knee extension moments explained 87.4% of the variance in MCF when considered alongside QF. The increase of MCF may be attributed to increased frontal-plane pedal reaction force moment arm. Future research may seek to implement QF modulation as a part of rehabilitation or training procedures utilizing cycling in cases where medial compartment joint loading is of importance.

A patient-specific 3D model of the knee to compare the femoral rollback before and after total knee arthroplasty (TKA)

Purpose

Total knee arthroplasty (TKA) is nowadays performed as a standard procedure on a large number of patients suffering from arthrosis. Replacing the knee joint causes changes in the geometry and kinematics of the knee, which are unique to each individual. This research focuses on the method to detect these changes after TKA and on the impact on the knee movement. This approach could reduce complications in patients with post-operative pain and reduce the number of revisions.

Methods

A 3D model of a patient’s knee was made by measuring the movement with a medically certified infrared stereo camera. This measurement was combined with the 3D model of the patient’s bones, previously segmented from the CT scan. This model is printed in 3D, one part being the mechanism that follows the movement of the patient, and the other part being the 3D copy of the femur and tibia bones. The knee replacement operation is performed directly on the model and the resulting rollback is being measured before and after TKA.

Results

We observe a difference in the rollback before and after TKA on the 3D printed model. The variation in size and shape of the femoral implant compared to the natural femur condyles is one of the reasons for the changes in the rollback effect. The rollback is half as large after the prosthesis insertion, which confirms the fact that the femoral prosthesis geometry influences the knee kinematics.

Conclusions

In this study, a first 3D model combining the patient-specific kinematic and the geometry of his bones has been constructed. This model allows the surgeon to validate the plan of the operation, but also to understand the problems and consequences generated by the prosthesis insertion. The rollback is one of the most important motion of the knee joint and this behavior could be quantified, providing comparative analysis of the knee joint before and after the operation. As a future study, the model could be used to analyse more parameters of the TKA such as the impact of different implantation methods.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40634-020-00319-6.

12 Degrees of Freedom Muscle Force Driven Fibril-Reinforced Poroviscoelastic Finite Element Model of the Knee Joint

Accurate knowledge of the joint kinematics, kinetics, and soft tissue mechanical responses is essential in the evaluation of musculoskeletal (MS) disorders. Since in vivo measurement of these quantities requires invasive methods, musculoskeletal finite element (MSFE) models are widely used for simulations. There are, however, limitations in the current approaches. Sequentially linked MSFE models benefit from complex MS and FE models; however, MS model's outputs are independent of the FE model calculations. On the other hand, due to the computational burden, embedded (concurrent) MSFE models are limited to simple material models and cannot estimate detailed responses of the soft tissue. Thus, first we developed a MSFE model of the knee with a subject-specific MS model utilizing an embedded 12 degrees of freedom (DoFs) knee joint with elastic cartilages in which included both secondary kinematic and soft tissue deformations in the muscle force estimation (inverse dynamics). Then, a muscle-force-driven FE model with fibril-reinforced poroviscoelastic cartilages and fibril-reinforced poroelastic menisci was used in series to calculate detailed tissue mechanical responses (forward dynamics). Second, to demonstrate that our workflow improves the simulation results, outputs were compared to results from the same FE models which were driven by conventional MS models with a 1 DoF knee, with and without electromyography (EMG) assistance. The FE model driven by both the embedded and the EMG-assisted MS models estimated similar results and consistent with experiments from literature, compared to the results estimated by the FE model driven by the MS model with 1 DoF knee without EMG assistance.

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

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

EMG-Assisted Muscle Force Driven Finite Element Model of the Knee Joint with Fibril-Reinforced Poroelastic Cartilages and Menisci

Abnormal mechanical loading is essential in the onset and progression of knee osteoarthritis. Combined musculoskeletal (MS) and finite element (FE) modeling is a typical method to estimate load distribution and tissue responses in the knee joint. However, earlier combined models mostly utilize static-optimization based MS models and muscle force driven FE models typically use elastic materials for soft tissues or analyze specific time points of gait. Therefore, here we develop an electromyography-assisted muscle force driven FE model with fibril-reinforced poro(visco)elastic cartilages and menisci to analyze knee joint loading during the stance phase of gait. Moreover, since ligament pre-strains are one of the important uncertainties in joint modeling, we conducted a sensitivity analysis on the pre-strains of anterior and posterior cruciate ligaments (ACL and PCL) as well as medial and lateral collateral ligaments (MCL and LCL). The model produced kinematics and kinetics consistent with previous experimental data. Joint contact forces and contact areas were highly sensitive to ACL and PCL pre-strains, while those changed less cartilage stresses, fibril strains, and fluid pressures. The presented workflow could be used in a wide range of applications related to the aetiology of cartilage degeneration, optimization of rehabilitation exercises, and simulation of knee surgeries.

EMG-Assisted Muscle Force Driven Finite Element Model of the Knee Joint with Fibril-Reinforced Poroelastic Cartilages and Menisci

Abnormal mechanical loading is essential in the onset and progression of knee osteoarthritis. Combined musculoskeletal (MS) and finite element (FE) modeling is a typical method to estimate load distribution and tissue responses in the knee joint. However, earlier combined models mostly utilize static-optimization based MS models and muscle force driven FE models typically use elastic materials for soft tissues or analyze specific time points of gait. Therefore, here we develop an electromyography-assisted muscle force driven FE model with fibril-reinforced poro(visco)elastic cartilages and menisci to analyze knee joint loading during the stance phase of gait. Moreover, since ligament pre-strains are one of the important uncertainties in joint modeling, we conducted a sensitivity analysis on the pre-strains of anterior and posterior cruciate ligaments (ACL and PCL) as well as medial and lateral collateral ligaments (MCL and LCL). The model produced kinematics and kinetics consistent with previous experimental data. Joint contact forces and contact areas were highly sensitive to ACL and PCL pre-strains, while those changed less cartilage stresses, fibril strains, and fluid pressures. The presented workflow could be used in a wide range of applications related to the aetiology of cartilage degeneration, optimization of rehabilitation exercises, and simulation of knee surgeries.

3D MRI Modeling of Thin and Spatially Complex Soft Tissue Structures without Shrinkage: Lamprey Myosepta as an Example

3D imaging techniques enable the nondestructive analysis and modeling of complex structures. Among these, MRI exhibits good soft tissue contrast, but is currently less commonly used for nonclinical research than X-ray CT, even though the latter requires contrast-staining that shrinks and distorts soft tissues. When the objective is the creation of a realistic and complete 3D model of soft tissue structures, MRI data are more demanding to acquire and visualize and require extensive post-processing because they comprise noncubic voxels with dimensions that represent a trade-off between tissue contrast and image resolution. Therefore, thin soft tissue structures with complex spatial configurations are not always visible in a single MRI dataset, so that standard segmentation techniques are not sufficient for their complete visualization. By using the example of the thin and spatially complex connective tissue myosepta in lampreys, we developed a workflow protocol for the selection of the appropriate parameters for the acquisition of MRI data and for the visualization and 3D modeling of soft tissue structures. This protocol includes a novel recursive segmentation technique for supplementing missing data in one dataset with data from another dataset to produce realistic and complete 3D models. Such 3D models are needed for the modeling of dynamic processes, such as the biomechanics of fish locomotion. However, our methodology is applicable to the visualization of any thin soft tissue structures with complex spatial configurations, such as fasciae, aponeuroses, and small blood vessels and nerves, for clinical research and the further exploration of tensegrity. Anat Rec, 301:1745-1763, 2018. © 2018 Wiley Periodicals, Inc.