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

Finite element analysis of endoscopic cross-overtop decompression for single-segment lumbar spinal stenosis based on real clinical cases

Introduction: For severe degenerative lumbar spinal stenosis (DLSS), the conventional percutaneous endoscopic translaminar decompression (PEID) has some limitations. The modified PEID, Cross-Overtop decompression, ensures sufficient decompression without excessive damage to the facet joints and posterior complex integrity. Objectives: To evaluate the biomechanical properties of Cross-Overtop and provide practical case validation for final decision-making in severe DLSS treatment. Methods: A finite element (FE) model of L4-L5 (M0) was established, and the validity was verified against prior studies. Endo-ULBD (M1), Endo-LOVE (M2), and Cross-Overtop (M3) models were derived from M0 using the experimental protocol. L4-L5 segments in each model were evaluated for the range of motion (ROM) and disc Von Mises stress extremum. The real clinical Cross-Overtop model was constructed based on clinical CT images, disregarding paraspinal muscle influence. Subsequent validation using actual FE analysis results enhances the credibility of the preceding virtual FE analysis. Results: Compared with M0, ROM in surgical models were less than 10°, and the growth rate of ROM ranged from 0.10% to 11.56%, while those of disc stress ranged from 0% to 15.75%. Compared with preoperative, the growth rate of ROM and disc stress were 2.66%-11.38% and 1.38%-9.51%, respectively. The ROM values in both virtual and actual models were less than 10°, verifying the affected segment stability after Cross-Overtop decompression. Conclusion: Cross-Overtop, designed for fully expanding the central canal and contralateral recess, maximizing the integrity of the facet joints and posterior complex, does no significant effect on the affected segmental biomechanics and can be recommended as an effective endoscopic treatment for severe DLSS.

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

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

Sensitivity analysis of paediatric knee kinematics to the graft surgical parameters during anterior cruciate ligament reconstruction: A sequentially linked neuromusculoskeletal-finite element analysis

BACKGROUND AND OBJECTIVE

Incidence of paediatric anterior cruciate ligament (ACL) rupture has increased substantially over recent decades. Following ACL rupture, ACL reconstruction (ACLR) surgery is typically performed to restore passive knee stability. This surgery involves replacing the failed ACL with a graft, however, surgeons must select from range of surgical parameters (e.g., type, size, insertion, and pre-tension) with no robust evidence guiding these decisions. This study presents a systemmatic computational approach to study effects of surgical parameter variation on kinematics of paediatric knees.

METHODS

This study used sequentially-linked neuromusculoskeletal (NMSK) finite element (FE) models of three paediatric knees to estimate the: (i) sensitivity of post-operative knee kinematics to four surgical parameters (type, size, insertion, and pre-tension) through multi-input multi-output sensitivity analysis; (ii) influence of motion and loading conditions throughout stance phase of walking gait on sensitivity indices; and (iii) influence of subject-specific anatomy (i.e., knee size) on sensitivivty indices. A previously validated FE model of the intact knee for each subject served as a reference against which ACLR knee kinematics were compared.

RESULTS

Sensitivity analyses revealed significant influences of surgical parameters on ACLR knee kinematics, albeit without discernible trend favouring any one parameter. Graft size and pre-tension were primary drivers of variation in knee translations and rotations, however, their effects fluctuated across stance indicating motion and loading conditions affect system sensitivity to surgical parameters. Importantly, the sensitivity of knee kinematics to surgical parameter varied across subjects, indicating geometry (i.e., knee size) influenced system sensitivity. Notably, alterations in graft parameters yielded substantial effects on kinematics (normalized root-mean-square-error > 10 %) compared to intact knee models, indicating surgical parameters vary post-operative knee kinematics.

CONCLUSIONS

Overall, this initial study highlights the importance of surgical parameter selection on post-operative kinematics in the paediatric ACLR knee, and provides evidence of the need for personalized surgical planning to ultimately enhance patient outcomes.

Osteoarthritis year in review 2023: Biomechanics

Biomechanics plays a significant yet complex role in osteoarthritis (OA) onset and progression. Identifying alterations in biomechanical factors and their complex interactions is critical for gaining new insights into OA pathophysiology and identification of clearly defined and modifiable mechanical treatment targets. This review synthesized biomechanics studies from March 2022 to April 2023, from which three themes relating to human gait emerged: (1) new insights into the pathogenesis of OA using computational modeling and machine learning, (2) technology-enhanced biomechanical interventions for OA, and (3) out-of-lab biomechanical assessments of OA. We further highlighted future-focused areas which may continue to advance the field of biomechanics in OA, with a particular emphasis on exploiting technology to understand and treat biomechanical mechanisms of OA outside the laboratory. The breadth of studies included in this review highlights the complex role of biomechanics in OA and showcase numerous innovative and outstanding contributions to the field. Exciting cross-disciplinary efforts integrating computational modeling, mobile sensors, and machine learning methods show great promise for streamlining in vivo multi-scale biomechanics workflows and are expected to underpin future breakthroughs in the understanding and treatment of biomechanics in OA.

Model Properties and Clinical Application in the Finite Element Analysis of Knee Joint: A Review

The knee is the most complex joint in the human body, including bony structures like the femur, tibia, fibula, and patella, and soft tissues like menisci, ligaments, muscles, and tendons. Complex anatomical structures of the knee joint make it difficult to conduct precise biomechanical research and explore the mechanism of movement and injury. The finite element model (FEM), as an important engineering analysis technique, has been widely used in many fields of bioengineering research. The FEM has advantages in the biomechanical analysis of objects with complex structures. Researchers can use this technology to construct a human knee joint model and perform biomechanical analysis on it. At the same time, finite element analysis can effectively evaluate variables such as stress, strain, displacement, and rotation, helping to predict injury mechanisms and optimize surgical techniques, which make up for the shortcomings of traditional biomechanics experimental research. However, few papers introduce what material properties should be selected for each anatomic structure of knee FEM to meet different research purposes. Based on previous finite element studies of the knee joint, this paper summarizes various modeling strategies and applications, serving as a reference for constructing knee joint models and research design.

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.

Biomechanical effects of typical lower limb movements of Chen-style Tai Chi on knee joint

The load and stress distribution on cartilage and meniscus of the knee joint in typical lower limb movements of Chen-style Tai Chi (TC) and deep squat (DS) were analyzed using finite element (FE) analysis. The loadings for this analysis consisted of muscle forces and ground reaction force (GRF), which were calculated through the inverse dynamic approach based on kinematics and force plate measurements obtained from motion capture experiments. Thirteen experienced practitioners performed four typical TC movements, namely, single whip (SW), brush knee and twist step (BKTS), stretch down (SD), and part the wild horse's mane (PWHM), which exhibit lower posture and greater lower limb force compared to other TC styles. The results indicated that TC required greater lower limb muscle strength than DS, resulting in greater knee joint forces. The stress on the medial cartilage in SW and BKTS fell within a range conductive to maintaining the balance between anabolism and catabolism of cartilage matrix. This was due to the fact that SW and BKTS reduce the medial to total tibiofemoral contact force ratios through knee abduction, which may effectively alleviate mild medial knee osteoarthritis (KOA). However, the greater medial contact force ratios observed in SD and PWHM resulted in great contact stresses that may aggravate the pain of patients with KOA. To mitigate these effects, practitioners should consider elevating their postures appropriately to reduce knee flexion angles, especially during the single-leg support phase. This adjustment can decrease the required muscle strength, load and stress on the knee joint.

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.

A Digital Twin Framework for Precision Neuromusculoskeletal Health Care: Extension Upon Industrial Standards

There is a powerful global trend toward deeper integration of digital twins into modern life driven by Industry 4.0 and 5.0. Defense, agriculture, engineering, manufacturing, and urban planning sectors have thoroughly incorporated digital twins to great benefit across their respective product lifecycles. Despite clear benefits, a digital twin framework for health and medical sectors is yet to emerge. This paper proposes a digital twin framework for precision neuromusculoskeletal health care. We build upon the International Standards Organization framework for digital twins for manufacturing by presenting best available computational models within a digital twin framework for clinical application. We map a use case for modeling Achilles tendon mechanobiology, highlighting how current modeling practices align with our proposed digital twin framework. Similarly, we map a use case for advanced neurorehabilitation technology, highlighting the role of a digital twin in control of systems where human and machine are interfaced. Future work must now focus on creating an informatic representation to govern how digital data are passed to, from, and within the digital twin, as well as specific standards to declare which measurement systems and modeling methods are acceptable to move toward widespread use of the digital twin framework for precision neuromusculoskeletal health care.

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.

Computational prediction of muscle synergy using a finite element framework for a musculoskeletal model on lower limb

Previous studies have demonstrated that the central nervous system activates muscles in module patterns to reduce the complexity needed to control each muscle while producing a movement, which is referred to as muscle synergy. In previous musculoskeletal modeling-based muscle synergy analysis studies, as a result of simplification of the joints, a conventional rigid-body link musculoskeletal model failed to represent the physiological interactions of muscle activation and joint kinematics. However, the interaction between the muscle level and joint level that exists in vivo is an important relationship that influences the biomechanics and neurophysiology of the musculoskeletal system. In the present, a lower limb musculoskeletal model coupling a detailed representation of a joint including complex contact behavior and material representations was used for muscle synergy analysis using a decomposition method of non-negative matrix factorization (NMF). The complexity of the representation of a joint in a musculoskeletal system allows for the investigation of the physiological interactions in vivo on the musculoskeletal system, thereby facilitating the decomposition of the muscle synergy. Results indicated that, the activities of the 20 muscles on the lower limb during the stance phase of gait could be controlled by three muscle synergies, and total variance accounted for by synergies was 86.42%. The characterization of muscle synergy and musculoskeletal biomechanics is consistent with the results, thus explaining the formational mechanism of lower limb motions during gait through the reduction of the dimensions of control issues by muscle synergy and the central nervous system.

Finite element analysis of a healthy knee joint at deep squatting for the study of tibiofemoral and patellofemoral contact

Background

In non-western countries, deep squatting is a daily activity, and prolonged deep squatting is common among occupational squatters. Household tasks, taking a bath, socializing, using toilets, and performing religious acts are among the activities frequently carried out while squatting by the Asian population. High knee loading is responsible for a knee injury and osteoarthritis. Finite element analysis is an effective tool to determine stresses on the knee joint.

Methods

Magnetic Resonance Imaging (MRI) and Computed Tomographic (CT) images were acquired of one adult without knee injuries. The CT images were acquired at the fully extended knee and one more set of images was acquired with the knee at a deeply flexed knee position. The MRI was acquired with the fully extended knee. 3-Dimensional models of bones were created using CT and soft tissue using MRI with the help of 3D Slicer software. Kinematics and finite element analysis of the knee was performed for standing and deep squatting posture using Ansys Workbench 2022.

Results

High peak stresses were observed at deep squatting compared to standing along with the reduction in the contact area. Peak von Mises stresses on femoral cartilage, tibial cartilage, patellar cartilage, and meniscus were increased from 3.3 MPa to 19.9 MPa, 2.9 MPa to 12.4 MPa, 1.5 MPa to 16.7 MPa and 15.8 MPa to 32.8 MPa respectively during deep squatting. Posterior translation of 7.01 mm, and 12.58 mm was observed for medial and lateral femoral condyle respectively from full extension to 153° knee flexion.

Conclusions

Increased stresses in the knee joint at deep squat posture may cause cartilage damage. A sustained deep squat posture should be avoided for healthy knee joints. More posterior translations of the medial femoral condyle at higher knee flexion angle warrant further investigation.

Maintaining soldier musculoskeletal health using personalised digital humans, wearables and/or computer vision

OBJECTIVES

The physical demands of military service place soldiers at risk of musculoskeletal injuries and are major concerns for military capability. This paper outlines the development new training technologies to prevent and manage these injuries.

DESIGN

Narrative review.

METHODS

Technologies suitable for integration into next-generation training devices were examined. We considered the capability of technologies to target tissue level mechanics, provide appropriate real-time feedback, and their useability in-the-field.

RESULTS

Musculoskeletal tissues' health depends on their functional mechanical environment experienced in military activities, training and rehabilitation. These environments result from the interactions between tissue motion, loading, biology, and morphology. Maintaining health of and/or repairing joint tissues requires targeting the "ideal" in vivo tissue mechanics (i.e., loading and strain), which may be enabled by real-time biofeedback. Recent research has shown that these biofeedback technologies are possible by integrating a patient's personalised digital twin and wireless wearable devices. Personalised digital twins are personalised neuromusculoskeletal rigid body and finite element models that work in real-time by code optimisation and artificial intelligence. Model personalisation is crucial in obtaining physically and physiologically valid predictions.

CONCLUSIONS

Recent work has shown that laboratory-quality biomechanical measurements and modelling can be performed outside the laboratory with a small number of wearable sensors or computer vision methods. The next stage is to combine these technologies into well-designed easy to use products.

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.