Comparison between kinetic and kinetic-kinematic driven knee joint finite element models

Finite element analysis of meniscus contact mechanical behavior based on kinematic simulation of abnormal gait

The meniscus plays a crucial role in the proper functioning of the knee joint, and when it becomes damaged, partial removal or replacement is necessary to restore proper function. Understanding the stress and deformation of the meniscus during various movements is essential for developing effective materials for meniscus repair. However, accurately estimating the contact mechanics of the knee joint can be challenging due to its complex shape and the dynamic changes it undergoes during movement. To address this issue, the open-source software SCONE can be used to establish a kinematics model that monitors the different states of the knee joint during human motion and obtains relevant gait kinematics data. To evaluate the stress and deformation of the meniscus during normal human movement, values of different states in the movement gait can be selected for finite element analysis (FEA) of the knee joint. This analysis enables researchers to assess changes in the meniscus. To evaluate meniscus damage, it is necessary to obtain changes in its mechanical behavior during abnormal movements. This information can serve as a reference for designing and optimizing the mechanical performance of materials used in meniscus repair and replacement.

Analysis of cartilage loading and injury correlation in knee varus deformity

Knee varus (KV) deformity leads to abnormal forces in the different compartments of the joint cavity and abnormal mechanical loading thus leading to knee osteoarthritis (KOA). This study used computer-aided design to create 3-dimensional simulation models of KOA with varying varus angles to analyze stress distribution within the knee joint cavity using finite element analysis for different varus KOA models and to compare intra-articular loads among these models. Additionally, we developed a cartilage loading model of static KV deformity to correlate with dynamic clinical cases of cartilage injury. Different KV angle models were accurately simulated with computer-aided design, and the KV angles were divided into (0°, 3°, 6°, 9°, 12°, 15°, and 18°) 7 knee models, and then processed with finite element software, and the Von-Mises stress distribution and peak values of the cartilage of the femoral condyles, medial tibial plateau, and lateral plateau were obtained by simulating the human body weight in axial loading while performing the static extension position. Finally, intraoperative endoscopy visualization of cartilage injuries in clinical cases corresponding to KV deformity subgroups was combined to find cartilage loading and injury correlations. With increasing varus angle, there was a significant increase in lower limb mechanical axial inward excursion and peak Von-Mises stress in the medial interstitial compartment. Analysis of patients' clinical data demonstrated a significant correlation between varus deformity angle and cartilage damage in the knee, medial plateau, and patellofemoral intercompartment. Larger varus deformity angles could be associated with higher medial cartilage stress loads and increased cartilage damage in the corresponding peak stress area. When the varus angle exceeds 6°, there is an increased risk of cartilage damage, emphasizing the importance of early surgical correction to prevent further deformity and restore knee function.

Integration of Swin UNETR and statistical shape modeling for a semi-automated segmentation of the knee and biomechanical modeling of articular cartilage

Simulation studies, such as finite element (FE) modeling, provide insight into knee joint mechanics without patient involvement. Generic FE models mimic the biomechanical behavior of the tissue, but overlook variations in geometry, loading, and material properties of a population. Conversely, subject-specific models include these factors, resulting in enhanced predictive precision, but are laborious and time intensive. The present study aimed to enhance subject-specific knee joint FE modeling by incorporating a semi-automated segmentation algorithm using a 3D Swin UNETR for an initial segmentation of the femur and tibia, followed by a statistical shape model (SSM) adjustment to improve surface roughness and continuity. For comparison, a manual FE model was developed through manual segmentation (i.e., the de-facto standard approach). Both FE models were subjected to gait loading and the predicted mechanical response was compared. The semi-automated segmentation achieved a Dice similarity coefficient (DSC) of over 98% for both the femur and tibia. Hausdorff distance (mm) between the semi-automated and manual segmentation was 1.4 mm. The mechanical results (max principal stress and strain, fluid pressure, fibril strain, and contact area) showed no significant differences between the manual and semi-automated FE models, indicating the effectiveness of the proposed semi-automated segmentation in creating accurate knee joint FE models. We have made our semi-automated models publicly accessible to support and facilitate biomechanical modeling and medical image segmentation efforts ( https://data.mendeley.com/datasets/k5hdc9cz7w/1 ).

Effects of gait modifications on tissue-level knee mechanics in individuals with medial tibiofemoral osteoarthritis: A proof-of-concept study towards personalized interventions

Gait modification is a common nonsurgical approach to alter the mediolateral distribution of knee contact forces, intending to decelerate or postpone the progression of mechanically induced knee osteoarthritis (KOA). Nevertheless, the success rate of these approaches is controversial, with no studies conducted to assess alterations in tissue-level knee mechanics governing cartilage degradation response in KOA patients undertaking gait modifications. Thus, here we investigated the effect of different conventional gait conditions and modifications on tissue-level knee mechanics previously suggested as indicators of collagen network damage, cell death, and loss of proteoglycans in knee cartilage. Five participants with medial KOA were recruited and musculoskeletal finite element analyses were conducted to estimate subject-specific tissue mechanics of knee cartilages during two gait conditions (i.e., barefoot and shod) and six gait modifications (i.e., 0°, 5°, and 10° lateral wedge insoles, toe-in, toe-out, and wide stance). Based on our results, the optimal gait modification varied across the participants. Overall, toe-in, toe-out, and wide stance showed the greatest reduction in tissue mechanics within medial tibial and femoral cartilages. Gait modifications could effectually alter maximum principal stress (~20 ± 7%) and shear strain (~9 ± 4%) within the medial tibial cartilage. Nevertheless, lateral wedge insoles did not reduce joint- and tissue-level mechanics considerably. Significance: This proof-of-concept study emphasizes the importance of the personalized design of gait modifications to account for biomechanical risk factors associated with cartilage degradation.

Validation and evaluation of subject-specific finite element models of the pediatric knee

Finite element (FE) models have been widely used to investigate knee joint biomechanics. Most of these models have been developed to study adult knees, neglecting pediatric populations. In this study, an atlas-based approach was employed to develop subject-specific FE models of the knee for eight typically developing pediatric individuals. Initially, validation simulations were performed at four passive tibiofemoral joint (TFJ) flexion angles, and the resulting TFJ and patellofemoral joint (PFJ) kinematics were compared to corresponding patient-matched measurements derived from magnetic resonance imaging (MRI). A neuromusculoskeletal-(NMSK)-FE pipeline was then used to simulate knee biomechanics during stance phase of walking gait for each participant to evaluate model simulation of a common motor task. Validation simulations demonstrated minimal error and strong correlations between FE-predicted and MRI-measured TFJ and PFJ kinematics (ensemble average of root mean square errors < 5 mm for translations and < 4.1° for rotations). The FE-predicted kinematics were strongly correlated with published reports (ensemble average of Pearson's correlation coefficients (ρ) > 0.9 for translations and ρ > 0.8 for rotations), except for TFJ mediolateral translation and abduction/adduction rotation. For walking gait, NMSK-FE model-predicted knee kinematics, contact areas, and contact pressures were consistent with experimental reports from literature. The strong agreement between model predictions and experimental reports underscores the capability of sequentially linked NMSK-FE models to accurately predict pediatric knee kinematics, as well as complex contact pressure distributions across the TFJ articulations. These models hold promise as effective tools for parametric analyses, population-based clinical studies, and enhancing our understanding of various pediatric knee injury mechanisms. They also support intervention design and prediction of surgical outcomes in pediatric populations.

Towards a Transferable Modeling Method of the Knee to Distinguish Between Future Healthy Joints from Osteoarthritic Joints: Data from the Osteoarthritis Initiative

Computational models can be used to predict the onset and progression of knee osteoarthritis. Ensuring the transferability of these approaches among computational frameworks is urgent for their reliability. In this work, we assessed the transferability of a template-based modeling strategy, based on the finite element (FE) method, by implementing it on two different FE softwares and comparing their results and conclusions. For that, we simulated the knee joint cartilage biomechanics of 154 knees using healthy baseline conditions and predicted the degeneration that occurred after 8 years of follow-up. For comparisons, we grouped the knees using their Kellgren-Lawrence grade at the 8-year follow-up time and the simulated volume of cartilage tissue that exceeded age-dependent thresholds of maximum principal stress. We considered the medial compartment of the knee in the FE models and used ABAQUS and FEBio FE softwares for simulations. The two FE softwares detected different volumes of overstressed tissue in corresponding knee samples (p < 0.01). However, both programs correctly distinguished between the joints that remained healthy and those that developed severe osteoarthritis after the follow-up (AUC = 0.73). These results indicate that different software implementations of a template-based modeling method similarly classify future knee osteoarthritis grades, motivating further evaluations using simpler cartilage constitutive models and additional studies on the reproducibility of these modeling strategies.

Signatures of disease progression in knee osteoarthritis: insights from an integrated multi-scale modeling approach, a proof of concept

Introduction: Knee osteoarthritis (KOA) is characterized by articular cartilage degeneration. It has been widely accepted that the mechanical joint environment plays a significant role in the onset and progression of this disease. In silico models have been used to study the interplay between mechanical loading and cartilage degeneration, hereby relying mainly on two key mechanoregulatory factors indicative of collagen degradation and proteoglycans depletion. These factors are the strain in collagen fibril direction (SFD) and maximum shear strain (MSS) respectively. Methods: In this study, a multi-scale in silico modeling approach was used based on a synergy between musculoskeletal and finite element modeling to evaluate the SFD and MSS. These strains were evaluated during gait based on subject-specific gait analysis data collected at baseline (before a 2-year follow-up) for a healthy and progressive early-stage KOA subject with similar demographics. Results: The results show that both SFD and MSS factors allowed distinguishing between a healthy subject and a KOA subject, showing progression at 2 years follow-up, at the instance of peak contact force as well as during the stance phase of the gait cycle. At the peak of the stance phase, the SFD were found to be more elevated in the KOA patient with the median being 0.82% higher in the lateral and 0.4% higher in the medial compartment of the tibial cartilage compared to the healthy subject. Similarly, for the MSS, the median strains were found to be 3.6% higher in the lateral and 0.7% higher in the medial tibial compartment of the KOA patient compared to the healthy subject. Based on these intersubject SFD and MSS differences, we were additionally able to identify that the tibial compartment of the KOA subject at risk of progression. Conclusion/discussion: We confirmed the mechanoregulatory factors as potential biomarkers to discriminate patients at risk of disease progression. Future studies should evaluate the sensitivity of the mechanoregulatory factors calculated based on this multi-scale modeling workflow in larger patient and control cohorts.

Effect of Whole Body Parameters on Knee Joint Biomechanics: Implications for ACL Injury Prevention During Single-Leg Landings

BACKGROUND

Previous studies have examined the effect of whole body (WB) parameters on anterior cruciate ligament (ACL) strain and loads, as well as knee joint kinetics and kinematics. However, articular cartilage damage occurs in relation to ACL failure, and the effect of WB parameters on ACL strain and articular cartilage biomechanics during dynamic tasks is unclear.

PURPOSES

(1) To investigate the effect of WB parameters on ACL strain, as well as articular cartilage stress and contact force, during a single-leg cross drop (SLCD) and single-leg drop (SLD). (2) To identify WB parameters predictive of high ACL strain during these tasks.

STUDY DESIGN

Descriptive laboratory study.

METHODS

Three-dimensional motion analysis data from 14 physically active men and women were recorded during an SLCD and SLD. OpenSim was used to obtain their kinematics, kinetics, and muscle forces for the WB model. Using these data in kinetically driven finite element simulations of the knee joint produced outputs of ACL strains and articular cartilage stresses and contact forces. Spearman correlation coefficients were used to assess relationships between WB parameters and ACL strain and cartilage biomechanics. Moreover, receiver operating characteristic curve analyses and multivariate binary logistic regressions were used to find the WB parameters that could discriminate high from low ACL strain trials.

RESULTS

Correlations showed that more lumbar rotation away from the stance limb at peak ACL strain had the strongest overall association (ρ = 0.877) with peak ACL strain. Higher knee anterior shear force (ρ = 0.895) and lower gluteus maximus muscle force (ρ = 0.89) at peak ACL strain demonstrated the strongest associations with peak articular cartilage stress or contact force in ≥1 of the analyzed tasks. The regression model that used muscle forces to predict high ACL strain trials during the dominant limb SLD yielded the highest accuracy (93.5%), sensitivity (0.881), and specificity (0.952) among all regression models.

CONCLUSION

WB parameters that were most consistently associated with and predictive of high ACL strain and poor articular cartilage biomechanics during the SLCD and SLD tasks included greater knee abduction angle at initial contact and higher anterior shear force at peak ACL strain, as well as lower gracilis, gluteus maximus, and medial gastrocnemius muscle forces.

CLINICAL RELEVANCE

Knowledge of which landing postures create a high risk for ACL or cartilage injury may help reduce injuries in athletes by avoiding those postures and practicing the tasks with reduced high-risk motions, as well as by strengthening the muscles that protect the knee during single-leg landings.

Computational assessment on the impact of collagen fiber orientation in cartilages on healthy and arthritic knee kinetics/kinematics

BACKGROUND

The inhomogeneous distribution of collagen fiber in cartilage can substantially influence the knee kinematics. This becomes vital for understanding the mechanical response of soft tissues, and cartilage deterioration including osteoarthritis (OA). Though the conventional computational models consider geometrical heterogeneity along with fiber reinforcements in the cartilage model as material heterogeneity, the influence of fiber orientation on knee kinetics and kinematics is not fully explored. This work examines how the collagen fiber orientation in the cartilage affects the healthy (intact knee) and arthritic knee response over multiple gait activities like running and walking.

METHODS

A 3D finite element knee joint model is used to compute the articular cartilage response during the gait cycle. A fiber-reinforced porous hyper elastic (FRPHE) material is used to model the soft tissue. A split-line pattern is used to implement the fiber orientation in femoral and tibial cartilage. Four distinct intact cartilage models and three OA models are simulated to assess the impact of the orientation of collagen fibers in a depth wise direction. The cartilage models with fibers oriented in parallel, perpendicular, and inclined to the articular surface are investigated for multiple knee kinematics and kinetics.

FINDINGS

The comparison of models with fiber orientation parallel to articulating surface for walking and running gait has the highest elastic stress and fluid pressure compared with inclined and perpendicular fiber-oriented models. Also, the maximum contact pressure is observed to be higher in the case of intact models during the walking cycle than for OA models. In contrast, the maximum contact pressure is higher during running in OA models than in intact models. Additionally, parallel-oriented models produce higher maximum stresses and fluid pressure for walking and running gait than proximal-distal-oriented models. Interestingly, during the walking cycle, the maximum contact pressure with intact models is approximately three times higher than on OA models. In contrast, the OA models exhibit higher contact pressure during the running cycle.

INTERPRETATION

Overall, the study indicates that collagen orientation is crucial for tissue responsiveness. This investigation provides insights into the development of tailored implants.

Monitoring Knee Contact Force with Force-Sensing Insoles

Numerous applications exist for monitoring knee contact force (KCF) throughout activities of daily living. However, the ability to estimate these forces is restricted to a laboratory setting. The purposes of this study are to develop KCF metric estimation models and explore the feasibility of monitoring KCF metrics via surrogate measures derived from force-sensing insole data. Nine healthy subjects (3F, age 27 ± 5 years, mass 74.8 ± 11.8 kg, height 1.7 ± 0.08 m) walked at multiple speeds (0.8-1.6 m/s) on an instrumented treadmill. Thirteen insole force features were calculated as potential predictors of peak KCF and KCF impulse per step, estimated with musculoskeletal modeling. The error was calculated with median symmetric accuracy. Pearson product-moment correlation coefficients defined the relationship between variables. Models develop per-limb demonstrated lower prediction error than those developed per-subject (KCF impulse: 2.2% vs 3.4%; peak KCF: 3.50% vs. 6.5%, respectively). Many insole features are moderately to strongly associated with peak KCF, but not KCF impulse across the group. We present methods to directly estimate and monitor changes in KCF using instrumented insoles. Our results carry promising implications for internal tissue loads monitoring outside of a laboratory with wearable sensors.

Robust automatic hexahedral cartilage meshing framework enables population-based computational studies of the knee

Osteoarthritis of the knee is increasingly prevalent as our population ages, representing an increasing financial burden, and severely impacting quality of life. The invasiveness of in vivo procedures and the high cost of cadaveric studies has left computational tools uniquely suited to study knee biomechanics. Developments in deep learning have great potential for efficiently generating large-scale datasets to enable researchers to perform population-sized investigations, but the time and effort associated with producing robust hexahedral meshes has been a limiting factor in expanding finite element studies to encompass a population. Here we developed a fully automated pipeline capable of taking magnetic resonance knee images and producing a working finite element simulation. We trained an encoder-decoder convolutional neural network to perform semantic image segmentation on the Imorphics dataset provided through the Osteoarthritis Initiative. The Imorphics dataset contained 176 image sequences with varying levels of cartilage degradation. Starting from an open-source swept-extrusion meshing algorithm, we further developed this algorithm until it could produce high quality meshes for every sequence and we applied a template-mapping procedure to automatically place soft-tissue attachment points. The meshing algorithm produced simulation-ready meshes for all 176 sequences, regardless of the use of provided (manually reconstructed) or predicted (automatically generated) segmentation labels. The average time to mesh all bones and cartilage tissues was less than 2 min per knee on an AMD Ryzen 5600X processor, using a parallel pool of three workers for bone meshing, followed by a pool of four workers meshing the four cartilage tissues. Of the 176 sequences with provided segmentation labels, 86% of the resulting meshes completed a simulated flexion-extension activity. We used a reserved testing dataset of 28 sequences unseen during network training to produce simulations derived from predicted labels. We compared tibiofemoral contact mechanics between manual and automated reconstructions for the 24 pairs of successful finite element simulations from this set, resulting in mean root-mean-squared differences under 20% of their respective min-max norms. In combination with further advancements in deep learning, this framework represents a feasible pipeline to produce population sized finite element studies of the natural knee from subject-specific models.

Computer-aided Design of Distal Femoral Osteotomy for the Valgus Knee and Effect of Correction Angle on Joint Loading by Finite Element Analysis

Objective

Lateral open‐wedge distal femoral osteotomy (DFO) has been used to treat valgus deformity of the knee, with good clinical outcomes. However, there is a lack of biomechanical studies regarding the angle of correction. The objective of this study was to apply computer‐aided design (CAD) for osteotomy planning in a three‐dimensional (3D) anatomical model and to assess the biomechanical differences among the varying correction angles on joint loading by finite element analysis (FEA).

Methods

To model different angles of lateral open‐wedge DFO correction, the CAD software package Mimics 21.0 was used to accurately simulate the operated knee. The femur was cut to 0°, 2°, 4°, 6°, 8°, and 10° of varus (equivalent to hip‐knee‐ankle angles of 180°, 178°, 176°, 174°, 172°, and 170°, respectively). The original knee model and the corrected models were processed by FE software. Then, the FE models were subjected to an axial force to obtain the von Mises stress (VMS) and shear stress distributions within the femoral cartilages and menisci.

Results

Under a compressive load of 740 N, the highest VMS in lateral and medial compartments of the intact knee model was 3.418 and 3.303 MPa. The maximum value of both the VMS and the shear stress in the lateral compartment decreased as the varus angle increased, but the corresponding values in the medial compartment increased. When the hip‐knee‐ankle (HKA) angle was 180°, the VMS in the lateral and medial compartments was balanced (3.418 and 3.303 MPa, respectively). Meanwhile, when the HKA angle was 178° (3.488 and 3.625 MPa, respectively), the shear stress in the lateral and medial compartments was balanced. In addition, the magnitude of change in the stress was significantly higher in the medial compartment (90.9%) than in the lateral compartment (19.3%).

Conclusion

The optimal correction angle of the valgus knee is close to neutral alignment or slightly varus (0° ‐ 2°). Overcorrection is not recommended, as it can result in a steep increase of the stress within the medial compartment and may accelerate the process of medial compartment OA.

Influence of material heterogeneity on the mechanical response of articulated cartilages in a knee joint

Structurally, the articular cartilages are heterogeneous owing to nonuniform distribution and orientation of its constituents. The oversimplification of this soft tissue as a homogeneous material is generally considered in the simulation domain to estimate contact pressure along with other physical responses. Hence, there is a need for investigating knee cartilages for their actual response to external stimuli. In this article, impact of material and geometrical heterogeneity of the cartilage is resolved using well known material models. The findings are compared with conventional homogeneous models. The results indicate vital differences in contact pressure distribution and tissue deformation. Further, this study paves way for standardizing material models to extract maximum information possible for investigating knee mechanics with variable geometry and case specific parameters.

Subject-specific biomechanical analysis to estimate locations susceptible to osteoarthritis-Finite element modeling and MRI follow-up of ACL reconstructed patients

The aims of this case-control study were to: (1) Identify cartilage locations and volumes at risk of osteoarthritis (OA) using subject-specific finite element (FE) models; (2) Quantify the relationships between the simulated biomechanical parameters and T₂ and T_1ρ relaxation times of magnetic resonance imaging (MRI). We created subject-specific FE models for seven patients with anterior cruciate ligament (ACL) reconstruction and six controls based on a previous proof-of-concept study. We identified locations and cartilage volumes susceptible to OA, based on maximum principal stresses and absolute maximum shear strains in cartilage exceeding thresholds of 7 MPa and 32%, respectively. The locations and volumes susceptible to OA were compared qualitatively and quantitatively against 2-year longitudinal changes in T₂ and T_1ρ relaxation times. The degeneration volumes predicted by the FE models, based on excessive maximum principal stresses, were significantly correlated (r = 0.711, p < 0.001) with the degeneration volumes determined from T₂ relaxation times. There was also a significant correlation between the predicted stress values and changes in T₂ relaxation time (r = 0.649, p < 0.001). Absolute maximum shear strains and changes in T_1ρ relaxation time were not significantly correlated. Five out of seven patients with ACL reconstruction showed excessive maximum principal stresses in either one or both tibial cartilage compartments, in agreement with follow-up information from MRI. Expectedly, for controls, the FE models and follow-up information showed no degenerative signs. Our results suggest that the presented modelling methodology could be applied to prospectively identify ACL reconstructed patients at risk of biomechanically driven OA, particularly by the analysis of maximum principal stresses of cartilage.

A musculoskeletal finite element model of rat knee joint for evaluating cartilage biomechanics during gait

Abnormal loading of the knee due to injuries or obesity is thought to contribute to the development of osteoarthritis (OA). Small animal models have been used for studying OA progression mechanisms. However, numerical models to study cartilage responses under dynamic loading in preclinical animal models have not been developed. Here we present a musculoskeletal finite element model of a rat knee joint to evaluate cartilage biomechanical responses during a gait cycle. The rat knee joint geometries were obtained from a 3-D MRI dataset and the boundary conditions regarding loading in the joint were extracted from a musculoskeletal model of the rat hindlimb. The fibril-reinforced poroelastic (FRPE) properties of the rat cartilage were derived from data of mechanical indentation tests. Our numerical results showed the relevance of simulating anatomical and locomotion characteristics in the rat knee joint for estimating tissue responses such as contact pressures, stresses, strains, and fluid pressures. We found that the contact pressure and maximum principal strain were virtually constant in the medial compartment whereas they showed the highest values at the beginning of the gait cycle in the lateral compartment. Furthermore, we found that the maximum principal stress increased during the stance phase of gait, with the greatest values at midstance. We anticipate that our approach serves as a first step towards investigating the effects of gait abnormalities on the adaptation and degeneration of rat knee joint tissues and could be used to evaluate biomechanically-driven mechanisms of the progression of OA as a consequence of joint injury or obesity.

Osteoarthritis is a disease of the musculoskeletal system which is characterized by the degradation of articular cartilage. Changes in the knee loading after injuries or obesity contribute to the development of cartilage degeneration. Since injured cartilage cannot be reversed back to intact conditions, small animal models have been widely used for investigating osteoarthritis progression mechanisms. Moreover, experimental studies have been complemented with numerical models to overcome inherent limitations such as cost, difficulties to obtain accurate measures and replicate degenerative situations in the knee joint. However, computational models to study articular cartilage responses under dynamic loading in small animal models have not been developed. Thus, here we present a musculoskeletal finite element model (MSFE) of a rat knee joint to evaluate cartilage biomechanical responses during gait. Our computational model considers both the anatomical and locomotion characteristics of the rat knee joint for estimating mechanical responses in the articular cartilage. We suggest that our approach can be used to investigate tissue adaptations based on the mechanobiological responses of the cartilage to prevent the progression of osteoarthritis.

Biomechanics of the lateral meniscus: evidences from narrative review

Lateral meniscus plays a crucial role in the knee stability and function. Although complex settings are required for lateral meniscus biomechanical assessment, such characteristics have been investigated over time to optimize the daily clinical practice surgical procedure. The aim of the present study was to provide an overview of literature knowledge regarding the biomechanics of lateral meniscus and give further insights about novel experimental analyses. Studies regarding lateral meniscus mobility, effect on knee laxity, and contact mechanics were included in the review. The effect of meniscal lesion, meniscal repair, partial meniscectomy, and meniscus allograft transplantation were retrieved for either in vivo, cadaveric, or in-silico computational settings. Knee laxity was evaluated both in presence of isolated meniscal tears and in anterior cruciate ligament (ACL) deficient knees. Surgical navigation systems and robots were mainly used for the investigations in vivo and in cadaveric studies, while accelerometers emerged as an alternative for ambulatory assessments. Contact mechanics was only assessed in cadaveric and computational studies. Great effort has been put into exploring lateral meniscus biomechanics from multiple perspectives. Strong evidence emerged regarding the importance of lateral meniscus repair: increased knee stability at high degrees of knee flexion and reduction of peak contact pressures on tibial cartilage were the most reported benefits. The meniscus allograft transplantation also emerged as a concrete solution for irreparable tears and further studies are needed to investigate its long-term influence on knee stability.

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.

An EMG-assisted muscle-force driven finite element analysis pipeline to investigate joint- and tissue-level mechanical responses in functional activities: towards a rapid assessment toolbox

Joint tissue mechanics (e.g., stress and strain) are believed to have a major involvement in the onset and progression of musculoskeletal disorders, e.g., knee osteoarthritis (KOA). Accordingly, considerable efforts have been made to develop musculoskeletal finite element (MS-FE) models to estimate highly detailed tissue mechanics that predict cartilage degeneration. However, creating such models is time-consuming and requires advanced expertise. This limits these complex, yet promising MS-FE models to research applications with few participants and makes the models impractical for clinical assessments. Also, these previously developed MS-FE models have not been used to assess activities other than gait. This study introduces and verifies a semi-automated rapid state-of-the-art MS-FE modeling and simulation toolbox incorporating an electromyography- (EMG) assisted MS model and a muscle-force driven FE model of the knee with fibril-reinforced poro(visco)elastic cartilages and menisci. To showcase the usability of the pipeline, we estimated joint- and tissue-level knee mechanics in 15 KOA individuals performing different daily activities. The pipeline was verified by comparing the estimated muscle activations and joint mechanics to existing experimental data. To determine the importance of EMG-assisted MS approach, results were compared to those from the same FE models but driven by static-optimization-based MS models. The EMG-assisted MS-FE pipeline bore a closer resemblance to experiments compared to the static-optimization-based MS-FE pipeline. Importantly, the developed pipeline showed great potential as a rapid MS-FE analysis toolbox to investigate multiscale knee mechanics during different activities of individuals with KOA. The template FE model of the study is freely available here.

Sealing mechanism study of laryngeal mask airways via 3D modelling and finite element analysis

Proper sealing of laryngeal mask airways (LMAs) is critical for airway management in clinical use. Understanding the sealing mechanism can significantly help front-line anaesthetists to reduce the incidence of adverse events. However, anaesthetists, who may not have the most substantial engineering backgrounds, lack intuitive ways to develop an understanding of the LMA sealing mechanism effectively. The paper aims to study the LMA-pharynx sealing mechanisms from the perspective of front-line anaesthetists. We use a computer-aided 3D modelling technique to visualise the LMA—pharynx interactions, which helps anaesthetists identify the critical areas of complications. Furthermore, we conduct a quantitative pressure distribution analysis of the LMA-pharynx contacting surface using the finite element analysis technique, which helps further understand the sealing mechanics in those areas. We present two cases studies based on one male volunteer, aged 50, inserted with a ProSeal LMA. In the first case, a relatively low cuff pressure (CP) was applied to simulate the clinical circumstances in which complications related to air leakage are most likely to happen; in the second case, we increase the CP to a relatively high value to simulate the scenarios with an increased risk of complications related to high mucosal pressure. The experiments suggest the follows: (1) Sore throat complications related to high mucosal pressure is most likely to occur in the hypopharynx with a high CP setting, particularly in the areas where the cricoid cartilage presses the mucosa. (2) The narrow hyoid bone super horn width likely causes LMA insertion difficulties. (3) Insufficient CP will significantly increase the risk of air leakage in the oropharynx. A complete sealing pressure line in the contacting surface will be formed with sufficient CP, thereby preventing the air leakage into the oral.

Magnetic Resonance Imaging-based biomechanical simulation of cartilage: A systematic review

MRI-based mathematical and computational modeling studies can contribute to a better understanding of the mechanisms governing cartilage's mechanical performance and cartilage disease. In addition, distinct modeling of cartilage is needed to optimize artificial cartilage production. These studies have opened up the prospect of further deepening our understanding of cartilage function. Furthermore, these studies reveal the initiation of an engineering-level approach to how cartilage disease affects material properties and cartilage function. Aimed at researchers in the field of MRI-based cartilage simulation, research articles pertinent to MRI-based cartilage modeling were identified, reviewed, and summarized systematically. Various MRI applications for cartilage modeling are highlighted, and the limitations of different constitutive models used are addressed. In addition, the clinical application of simulations and studied diseases are discussed. The paper's quality, based on the developed questionnaire, was assessed, and out of 79 reviewed papers, 34 papers were determined as high-quality. Due to the lack of the best constitutive models for various clinical conditions, researchers may consider the effect of constitutive material models on the cartilage disease simulation. In the future, research groups may incorporate various aspects of machine learning into constitutive models and MRI data extraction to further refine the study methodology. Moreover, researchers should strive for further reproducibility and rigorous model validation and verification, such as gait analysis.

"Fishing net" suture augmenting tension-band wiring fixation in the treatment of inferior pole fracture of the patella

INTRODUCTION

Inferior pole fracture of the patella (IPFP) is difficult to repair and stabilize clinically. Although various fixation techniques have been developed, fixation strength and mobility remain daunting challenges to orthopedic surgeons. The goal of this research is to evaluate the biomechanical strength and clinical outcomes of a novel "fishing net" suture fixation procedure.

MATERIALS AND METHODS

Four finite element models, modified tension-band wiring fixation, anchor suture fixation, basket plate fixation and "fishing net" suture fixation were built to compare the fixing efficacy of "fishing net" suture fixation with three other fixation methods during IPFP fixation. From January 2018 to February 2019, 17 patients who suffered IPFP and treated by "fishing net" suture (FNS) fixation were compared with 20 patients treated by tension-band wiring (TBW) fixation in database and the two groups were evaluated postoperatively using the modified Cincinnati knee rating system.

RESULTS

Biomechanical evaluation showed that the relative displacement values of proximal patella measured by three pairs of points on both sides of the fracture line were the lowest using the "fishing net" suture fixation, while fixation using tension-band wiring and basket plate showed similar levels of stability that were less desirable than the "fishing net" method. As to clinical outcomes, there were 17 (100%) patients exhibited excellent or good results with no internal fixation failures in the FNS group compared to three internal fixation failures in the TBW group.

CONCLUSION

The biomechanical and clinical results suggest that the "fishing net" suture fixation is a viable candidate for fixation of IPFP.

Biomechanics of the Femoral Head Cartilage and Subchondral Trabecular Bone in Osteoporotic and Osteopenic Fractures

This study aimed to investigate the relationship between the micro structural properties of the subchondral trabecular bone (STB) and the macro mechanical properties of the articular cartilage (AC) in patients with osteoporotic (OP) and osteopenic (OPE) fractures. Sixteen femoral head samples (OP;OPE, n = 8 each) were obtained from female patients who underwent hip hemiarthroplasty. STB and AC specimens were harvested from those heads. Bone specimens were scanned using µ-CT to determine the micro structural properties. In-situ nondestructive compressive tests were performed for the cartilages to obtain elastic properties. The finite element technique was implemented on STB models created from µ-CT data to compute apparent elastic modulus. In addition, dynamic cyclic destructive tests were performed on STB and AC specimens to assess failure cycles. The results demonstrated that STB specimens in OPE group have more interconnected structure and higher cyclic dynamic strength than those in OP group. Furthermore, bone mineral density, failure cycle, and trabecular number of STB were positively correlated with the cartilage failure cycle, which indicates that STB alteration may affect the macroscopic mechanical properties of AC. The findings suggest that STB loss correlates with a decrease in cartilage strength and that improving of bone quality may prevent cartilage weakness.

Prediction of local fixed charge density loss in cartilage following ACL injury and reconstruction: A computational proof-of-concept study with MRI follow-up

The purpose of this proof-of-concept study was to develop three-dimensional patient-specific mechanobiological knee joint models to simulate alterations in the fixed charged density (FCD) around cartilage lesions during the stance phase of the walking gait. Two patients with anterior cruciate ligament (ACL) reconstructed knees were imaged at 1 and 3 years after surgery. The magnetic resonance imaging (MRI) data were used for segmenting the knee geometries, including the cartilage lesions. Based on these geometries, finite element (FE) models were developed. The gait of the patients was obtained using a motion capture system. Musculoskeletal modeling was utilized to calculate knee joint contact and lower extremity muscle forces for the FE models. Finally, a cartilage adaptation algorithm was implemented in both FE models. In the algorithm, it was assumed that excessive maximum shear and deviatoric strains (calculated as the combination of principal strains), and fluid velocity, are responsible for the FCD loss. Changes in the longitudinal T_1ρ and T₂ relaxation times were postulated to be related to changes in the cartilage composition and were compared with the numerical predictions. In patient 1 model, both the excessive fluid velocity and strain caused the FCD loss primarily near the cartilage lesion. T_1ρ and T₂ relaxation times increased during the follow-up in the same location. In contrast, in patient 2 model, only the excessive fluid velocity led to a slight FCD loss near the lesion, where MRI parameters did not show evidence of alterations. Significance: This novel proof-of-concept study suggests mechanisms through which a local FCD loss might occur near cartilage lesions. In order to obtain statistical evidence for these findings, the method should be investigated with a larger cohort of subjects.

Development and evaluation of a new procedure for subject-specific tensioning of finite element knee ligaments

Subject-specific tensioning of ligaments is essential for the stability of the knee joint and represents a challenging aspect in the development of finite element models. We aimed to introduce and evaluate a new procedure for the quantification of ligament prestrains from biplanar X-ray and CT data. Subject-specific model evaluation was performed by comparing predicted femorotibial kinematics with the in vitro response of six cadaveric specimens. The differences obtained using personalized models were comparable to those reported in similar studies in the literature. This study is the first step toward the use of simplified, personalized knee FE models in clinical context such as ligament balancing.

Rapid CT-based Estimation of Articular Cartilage Biomechanics in the Knee Joint Without Cartilage Segmentation

Knee osteoarthritis (OA) is a painful joint disease, causing disabilities in daily activities. However, there is no known cure for OA, and the best treatment strategy might be prevention. Finite element (FE) modeling has demonstrated potential for evaluating personalized risks for the progression of OA. Current FE modeling approaches use primarily magnetic resonance imaging (MRI) to construct personalized knee joint models. However, MRI is expensive and has lower resolution than computed tomography (CT). In this study, we extend a previously presented atlas-based FE modeling framework for automatic model generation and simulation of knee joint tissue responses using contrast agent-free CT. In this method, based on certain anatomical dimensions measured from bone surfaces, an optimal template is selected and scaled to generate a personalized FE model. We compared the simulated tissue responses of the CT-based models with those of the MRI-based models. We show that the CT-based models are capable of producing similar tensile stresses, fibril strains, and fluid pressures of knee joint cartilage compared to those of the MRI-based models. This study provides a new methodology for the analysis of knee joint and cartilage mechanics based on measurement of bone dimensions from native CT scans.

Identification of locations susceptible to osteoarthritis in patients with anterior cruciate ligament reconstruction: Combining knee joint computational modelling with follow-up T1ρ and T2 imaging

BACKGROUND

Finite element modelling can be used to evaluate altered loading conditions and failure locations in knee joint tissues. One limitation of this modelling approach has been experimental comparison. The aims of this proof-of-concept study were: 1) identify areas susceptible to osteoarthritis progression in anterior cruciate ligament reconstructed patients using finite element modelling; 2) compare the identified areas against changes in T₂ and T_1ρ values between 1-year and 3-year follow-up timepoints.

METHODS

Two patient-specific finite element models of knee joints with anterior cruciate ligament reconstruction were created. The knee geometry was based on clinical magnetic resonance imaging and joint loading was obtained via motion capture. We evaluated biomechanical parameters linked with cartilage degeneration and compared the identified risk areas against T₂ and T_1ρ maps.

FINDINGS

The risk areas identified by the finite element models matched the follow-up magnetic resonance imaging findings. For Patient 1, excessive values of maximum principal stresses and shear strains were observed in the posterior side of the lateral tibial and femoral cartilage. For Patient 2, high values of maximum principal stresses and shear strains of cartilage were observed in the posterior side of the medial joint compartment. For both patients, increased T₂ and T_1ρ values between the follow-up times were observed in the same areas.

INTERPRETATION

Finite element models with patient-specific geometries and motions and relatively simple material models of tissues were able to identify areas susceptible to post-traumatic knee osteoarthritis. We suggest that the methodology presented here may be applied in large cohort studies.

Influence of Gait Cycle Loads on Stress Distribution at The Residual Limb/Socket Interface of Transfemoral Amputees: A Finite Element Analysis

A Finite Element Analysis (FEA) was performed to evaluate the interaction between residual limb and socket when considering the dynamic loads of the gait cycle. Fourteen transfemoral amputees participated in this study, where their residual limbs (i.e., soft tissues and bone), and their sockets were reconstructed. The socket and the femur were defined as elastic materials, while the bulk soft tissues were defined as a hyperelastic material. Each model included the donning, standing, and gait cycle phase, with load and boundary conditions applied accordingly. The influence of adding the dynamic loads related to the gait cycle were compared against the modelling of the static load equivalent to the standing position resulting in changes of 23% ± 19% in the maximum values and in an increase in the size of the regions where they were located. Additionally, the possible correspondence between comfort and the location of peak loadbearing at the residual-limb/socket interface was explored. Consequently, the comfort perceived by the patient could be estimated based on the locations of the maximum stresses (i.e., if they coincide with the pressure tolerant or sensitive regions of the residual limb).

Computational evaluation of altered biomechanics related to articular cartilage lesions observed in vivo

Chondral lesions provide a potential risk factor for development of osteoarthritis. Despite the variety of in vitro studies on lesion degeneration, in vivo studies that evaluate relation between lesion characteristics and the risk for the possible progression of OA are lacking. Here, we aimed to characterize different lesions and quantify biomechanical responses experienced by surrounding cartilage tissue. We generated computational knee joint models with nine chondral injuries based on clinical in vivo arthrographic computed tomography images. Finite element models with fibril-reinforced poro(visco)elastic cartilage and menisci were constructed to simulate physiological loading. Systematically, the lesions experienced increased peak values of maximum principal strain, maximum shear strain, and minimum principal strain in the surrounding chondral tissue (p < 0.01) compared with intact tissue. Depth, volume, and area of the lesion correlated with the maximum shear strain (p < 0.05, Spearman rank correlation coefficient ρ = 0.733-0.917). Depth and volume of the lesion correlated also with the maximum principal strain (p < 0.05, ρ = 0.767, and ρ = 0.717, respectively). However, the lesion area had non-significant correlation with this strain parameter (p = 0.06, ρ = 0.65). Potentially, the introduced approach could be developed for clinical evaluation of biomechanical risks of a chondral lesion and planning an intervention. Statement of Clinical Relevance: In this study, we computationally characterized different in vivo chondral lesions and evaluated their risk of cartilage degeneration. This information is vital in decision-making for intervention in order to prevent post-traumatic osteoarthritis. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res.