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

Predicting Knee Joint Contact Forces During Normal Walking Using Kinematic Inputs With a Long-Short Term Neural Network

Knee joint contact forces are commonly estimated via surrogate measures (i.e., external knee adduction moments or musculoskeletal modeling). Despite its capabilities, modeling is not optimal for clinicians or persons with limited experience. The purpose of this study was to design a novel prediction method for knee joint contact forces that is simplistic in terms of required inputs. This study included marker trajectories and instrumented knee forces during normal walking from the "Grand Challenge" (n = 6) and "CAMS" (n = 2) datasets. Inverse kinematics were used to derive stance phase hip (sagittal, frontal, transverse), knee (sagittal, frontal), ankle (sagittal), and trunk (frontal) kinematics. A long-short term memory network (LSTM) was created using matlab to predict medial and lateral knee force waveforms using combinations of the kinematics. The Grand Challenge and CAMS datasets trained and tested the network, respectively. Musculoskeletal modeling forces were derived using static optimization and joint reaction tools in OpenSim. Waveform accuracy was determined as the proportion of variance and root-mean-square error between network predictions and in vivo data. The LSTM network was highly accurate for medial forces (R2 = 0.77, RMSE = 0.27 BW) and required only frontal hip and knee and sagittal hip and ankle kinematics. Modeled medial force predictions were excellent (R2 = 0.77, RMSE = 0.33 BW). Lateral force predictions were poor for both methods (LSTM R2 = 0.18, RMSE = 0.08 BW; modeling R2 = 0.21, RMSE = 0.54 BW). The designed LSTM network outperformed most reports of musculoskeletal modeling, including those reached in this study, revealing knee joint forces can accurately be predicted by using only kinematic input variables.

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

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

Pre-operative Planning of High Tibial Osteotomy With ChatGPT: Are We There Yet?

INTRODUCTION

ChatGPT (Chat Generative Pre-trained Transformer), developed by OpenAI (San Francisco, CA, USA), has gained attention in the medical field. It has the potential to enhance and simplify tasks, such as preoperative planning in orthopedic surgery. We aimed to test ChatGPT's accuracy in measuring the angle of correction for high tibial osteotomy for cases planned and performed at a tertiary teaching hospital in Singapore.

MATERIALS AND METHODS

Peri-operative angular parameters from 114 consecutive patients who underwent medial opening wedge high tibial osteotomy (MOWHTO) were used to query ChatGPT 3.0. First ChatGPT 3.0 was queried on what information it required to plan a MOWHTO. Based on its response, pre-operative medial proximal tibial angle (MPTA) and joint line congruence angle (JLCA) were provided. ChatGPT 3.0 then responded with its recommended angle of correction. This was compared against the manually planned surgical correction by our fellowship-trained surgeon. A root mean square analysis was then performed to compare ChatGPT 3.0 and manual planning.

RESULTS

The root mean square error (RMSE) of ChatGPT 3.0 in predicting correction angle in MWHTO was 2.96, suggesting a very poor model fit.

CONCLUSION

Although ChatGPT 3.0 represents a significant breakthrough in large language models with extensive capabilities, it is not currently optimized to effectively perform complex pre-operative planning in orthopedic surgery, specifically in the context of MOWHTO. Further refinement and consideration of specific factors are necessary to enhance its accuracy and suitability for such applications.

Effect of subscapularis repair on joint contact forces based on degree of posterior-superior rotator cuff tear severity in reverse shoulder arthroplasty

Massive irreparable rotator cuff tears (RCTs) affect the clinical outcomes of reverse shoulder arthroplasty (RSA). However, the effects of subscapularis repair on the outcomes of RSA, based on the degree of posterior-superior RCTs, are unclear. This study aimed to examine the effect of subscapularis repair on three-dimensional joint contact forces (JCFs) based on the degree of posterior-superior RCT severity in RSA. Ten human in vivo experimental data were used as input to the musculoskeletal model. A six-degrees-of-freedom (DOF) anatomical shoulder model was developed and validated against three-dimensional JCFs. The 6-DOF musculoskeletal shoulder model of RSA was then developed by importing the reverse shoulder implant into the validated anatomical shoulder model. Based on the various types of posterior-superior RCT severity, inverse dynamic simulations of subscapularis-torn and subscapularis-repaired models of RSA were performed: from isolated supraspinatus tears to partial or massive tears of the infraspinatus and teres minor. The intact rotator cuff model of RSA was also simulated for comparison with the different types of models. Our results showed that the more posterior-superior RCTs progressed in RSA, the more superior JCFs were observed at 90°, 105°, and 120° abduction in the subscapularis-torn model. However, subscapularis repair decreased the superior JCF at those angles sufficiently. In addition, the teres minor muscle-tendon force increased as infraspinatus bundle tears progressed in both the subscapularis-torn and -repaired models, in order to compensate for the reduced force during abduction. However, the teres minor muscle-tendon force was not as high as that of the infraspinatus muscle-tendon, which could result in muscle force imbalance between repaired subscapularis and teres minor. Therefore, our results suggest that repairing the subscapularis and the repairable infraspinatus during RSA can improve glenohumeral joint stability in the superior-inferior direction by restoring muscle force balance between the anterior cuff (i.e., subscapularis) and posterior cuff (i.e., infraspinatus and teres minor). The findings of this study can help clinician decide whether to repair the rotator cuff during RSA to enhance joint stability.

Influence of knee flexion on early femoral fracture healing: A combined analysis of musculoskeletal dynamics and finite elements

BACKGROUND AND OBJECTIVES

Knee flexion causes a certain amount of misalignment and relative movement of the fractured ends of the femur fracture, and if the flexion angle is too large it will affect the stability of the fracture and the healing process, making it challenging to design a safe range of flexion. However, due to a lack of basic understanding of the effect of knee flexion on the mechanical environment at the fracture site, clinicians are often unable to provide an objective and safe range of motion in flexion based on subjective experience. The aim of this study was to evaluate the effect of knee flexion on plate and fracture healing using finite element analysis (FEA).

METHODS

A human musculoskeletal model was constructed based on CT scan data, and the mechanical properties of the fracture site were changed by adjusting the knee flexion angle. The joint forces, muscle forces and moments acting on the femur were obtained by inverse dynamics analysis, and the biomechanical properties of the fracture-plate system were analyzed using finite elements. A finite element model of the fracture-plate system without muscle loading was also constructed. The effect of knee flexion on the safety of plate fixation and fracture healing was evaluated in terms of the biomechanical properties of the plate and the interfragmentary motion of the fracture.

RESULTS

As the knee flexion angle increases, the von Mises stress of the locked compression plate (LCP) first increases, then decreases, then increases again. The deformation from compression bending to tension twisting occurs simultaneously. At 30° of flexion, shear interfragmentary motion (SIM) was dominant and inhibited fracture healing; at more than 45° of flexion, the plate was twisted and deformed to the lateral side of the body, and the fracture site underwent greater misalignment and relative motion, with destructive effects on bone scabs and healing tissues. If muscle loading is not taken into account, the plate will undergo predominantly bending deformation and will overestimate the interfragmentary strain in the far and near cortex.

CONCLUSIONS

Knee flexion causes the plate to deform from compression bending to extension and torsion, which has an important impact on the safety and healing process of the fracture, and this study provides a biomechanical basis to guide the clinician in the post-operative rehabilitation of femoral fractures in the clinical setting.

A neuromusculoskeletal modelling approach to bilateral hip mechanics due to unexpected lateral perturbations during overground walking

BACKGROUND

Current studies on how external perturbations impact gait dynamics have primarily focused on the changes in the body's center of mass (CoM) during treadmill walking. The biomechanical responses, in particular to the multi-planar hip joint coordination, following perturbations in overground walking conditions are not completely known.

METHODS

In this study, a customized gait-perturbing device was designed to impose controlled lateral forces onto the subject's pelvis during overground walking. The biomechanical responses of bilateral hips were simulated by subject-specific neuromusculoskeletal models (NMS) driven by in-vivo motion data, which were further evaluated by statistical parameter mapping (SPM) and muscle coactivation index (CI) analysis. The validity of the subject-specific NMS was confirmed through comparison with measured surface electromyographic signals.

RESULTS

Following perturbations, the sagittal-plane hip motions were reduced for the leading leg by 18.39° and for the trailing leg by 8.23°, while motions in the frontal and transverse plane were increased, with increased hip abduction for the leading leg by 10.71° and external rotation by 9.06°, respectively. For the hip kinetics, both the bilateral hip joints showed increased abductor moments during midstance (20%-30% gait cycle) and decreased values during terminal stance (38%-48%). Muscle CI in both sagittal and frontal planes was significantly decreased for perturbed walking (p < 0.05), except for the leading leg in the sagittal plane.

CONCLUSION

The distinctive phase-dependent biomechanical response of the hip demonstrated its coordinated control strategy for balance recovery due to gait perturbations. And the changes in muscle CI suggested a potential mechanism for rapid and precise control of foot placement through modulation of joint stiffness properties. These findings obtained during actual overground perturbation conditions could have implications for the improved design of wearable robotic devices for balance assistance.

Subject-specific tribo-contact conditions in total knee replacements: a simulation framework across scales

Fundamental knowledge about in vivo kinematics and contact conditions at the articulating interfaces of total knee replacements are essential for predicting and optimizing their behavior and durability. However, the prevailing motions and contact stresses in total knee replacements cannot be precisely determined using conventional in vivo measurement methods. In silico modeling, in turn, allows for a prediction of the loads, velocities, deformations, stress, and lubrication conditions across the scales during gait. Within the scope of this paper, we therefore combine musculoskeletal modeling with tribo-contact modeling. In the first step, we compute contact forces and sliding velocities by means of inverse dynamics approach and force-dependent kinematic solver based upon experimental gait data, revealing contact forces during healthy/physiological gait of young subjects. In a second step, the derived data are employed as input data for an elastohydrodynamic model based upon the finite element method full-system approach taking into account elastic deformation, the synovial fluid's hydrodynamics as well as mixed lubrication to predict and discuss the subject-specific pressure and lubrication conditions.

Comparison of joint load, motions and contact stress and bone-implant interface micromotion of three implant designs for total ankle arthroplasty

BACKGROUND AND OBJECTIVE

Loosening and wear are still the main problems for the failure of total ankle arthroplasty, which are closely related to the micromotion at the bone-implant interface and the contact stress and joint motions at the articular surfaces. Implant design is a key factor to influence the ankle force, motions, contact stress, and bone-implant interface micromotion. The purpose of this study is to evaluate the differences in these parameters of INBONE II, INFINITY, and a new anatomic ankle implant under the physiological walking gait of three patients.

METHODS

This was achieved by using an in-silico simulation framework combining patient-specific musculoskeletal multibody dynamics and finite element analysis. Each implant was implanted into the musculoskeletal multibody dynamics model, respectively, which was driven by the gait data to calculate ankle forces and motions. These were then used as the boundary conditions for the finite element model, and the contact stress and the bone-implant interface micromotions were calculated.

RESULTS

The total ankle contact forces were not significantly affected by articular surface geometries of ankle implants. The range of motion of the ankle joint implanted with INFINITY was a little larger than that with INBONE II. The anatomic ankle implant design produced a greater range of motion than INBONE II, especially the internal-external rotation. The fixation design of INFINITY achieved lower bone-implant interface micromotion compared with INBONE II. The anatomic ankle implant design produced smaller contact stress with no evident edge contact and a smaller tibia-implant interface micromotion. In addition, significant differences in the magnitudes and tendencies of total ankle contact forces and motions among different patients were found.

CONCLUSIONS

The articular surface geometry of ankle implants not only affected the ankle motions and contact stress distribution but also affected the bone-implant interface micromotions. The anatomic ankle implant had good performance in recovering ankle joint motion, equalizing contact stress, and reducing bone-implant interface micromotion. INFINITY's fixation design could achieve smaller bone-implant interface micromotion than INBONE II.

[Movement analysis and musculoskeletal simulation in non-union treatment-Experiences and first clinical results]

BACKGROUND

The mechanical boundary conditions of the non-union and osteosynthetic construct are a key determinant of fracture healing after revision surgery. Aim of this study was to introduce a movement analysis and simulation workflow to determine the mechanical conditions during non-union healing in vivo.

MATERIAL AND METHODS

On an individual case basis after non-union revision surgery we performed an accelerometry-based movement analysis. The results were then used as input for a musculoskeletal simulation of the non-union, osteosynthetic construct as well as adjacent joints mechanical boundary conditions.

RESULTS

A total of 13 patients were analyzed with our new workflow. The introduced protocol allows an in vivo determination of the mechanical boundary conditions. On clinical follow-up all patients showed radiographic consolidation of the non-union.

CONCLUSION

The introduced workflow allows a clinically applicable determination of the mechanical boundary conditions of fracture and non-union healing. Further studies can now determine the effect of the introduced technique for mechanically optimized postoperative aftercare regimes as well as biomechanically adapted surgical treatment.

A Computer-Assisted Approach Regarding the Optimization of the Geometrical Planning of Medial Opening Wedge High Tibial Osteotomy

Opening wedge high tibial osteotomy (OWHTO) is a surgical procedure often used to eliminate the effects of knee osteoarthritis, a disease that is becoming more widespread worldwide. Optimizing the geometric planning of this operation is a very important preparatory step for the success of the intervention and rapid postoperative recovery. This optimization is performed in two main directions. The first direction evaluates the intraoperative behavior of the tibia during the osteotomy by optimizing four geometric parameters that characterize geometric planning. The second direction aims at a postoperative evaluation of the flat tibia-osteosynthesis assembly taking into account the optimal position on the medial–lateral articular line through which the corrected mechanical axis of the tongue passes and implicitly offloads the transfer from the medial area to the side of the knee. The research methods used are exclusively computer-assisted such as: computer-aided design (hereinafter CAD) for geometric modeling of the tibia taking into account the real bone structure, the finite element method (hereinafter FEM) for performing numerical analyses and design of the experiment (hereinafter DOE) for the design of the research. The results obtained are eloquent and clearly presented and can be important elements for orthopedic doctors at the geometric planning stage of the OWHTO.

Effects of the Anterolateral Ligament and Anterior Cruciate Ligament on Knee Joint Mechanics: A Biomechanical Study Using Computational Modeling

Background:

Recent studies on lateral knee anatomy have reported the presence of a true ligament structure, the anterolateral ligament (ALL), in the anterolateral region of the knee joint. However, its biomechanical effects have not been fully elucidated.

Purpose:

To investigate, by using computer simulation, the association between the ALL and anterior cruciate ligament (ACL) under dynamic loading conditions.

Study Design:

Descriptive laboratory study; Level of evidence, 5.

Methods:

The authors combined medical imaging from 5 healthy participants with motion capture to create participant-specific knee models that simulated the entire 12 degrees of freedom of tibiofemoral (TF) and patellofemoral (PF) joint behaviors. These dynamic computational models were validated using electromyographic data, muscle activation data, and data from previous experimental studies. Forces exerted on the ALL with ACL deficiency and on the ACL with ALL deficiency, as well as TF and PF contact forces with deficiencies of the ACL, ALL, and the entire ligament structure, were evaluated under gait and squat loading. A single gait cycle and squat cycle were divided into 11 time points (periods 0.0-1.0). Simulated ligament forces and contact forces were compared using nonparametric repeated-measures Friedman tests.

Results:

Force exerted on the ALL significantly increased with ACL deficiency under both gait- and squat-loading conditions. With ACL deficiency, the mean force on the ALL increased by 129.7% under gait loading in the 0.4 period ( P < .05) and increased by 189% under high flexion during the entire cycle of squat loading ( P < .05). A similar trend of significantly increased force on the ACL was observed with ALL deficiency. Contact forces on the TF and PF joints with deficiencies of the ACL, ALL, and entire ligament structure showed a complicated pattern. However, contact force exerted on TF and PF joints with respect to deficiencies of ACL and ALL significantly increased under both gait- and squat-loading conditions.

Conclusion:

The results of this computer simulation study indicate that the ACL and the ALL of the lateral knee joint act as secondary stabilizers to each other under dynamic load conditions.

Geometrical Planning of the Medial Opening Wedge High Tibial Osteotomy—An Experimental Approach

This article presents an experimental approach to the geometrical planning of the medial opening wedge high tibial osteotomy surgery which, as it is known, is an efficient surgical strategy quite widely used in treating knee osteoarthritis. While most of the published papers focus on analyzing this surgery from a medical point of view, we suggest a postoperative experimental evaluation of the intervention from a biomechanical point of view. The geometrical planning and, more specifically, the determination of the point of intersection between the corrected mechanical axis and the medial-lateral articular line of the knee, is a problem quite often debated in literature. This paper aims to experimentally investigate the behavior of the tibia with an open wedge osteotomy fixed with a locking plate, TomoFix (DE Puy Synthes), taking into account two positions of the mechanical axis of the leg on the width of the tibial plateau, measured from medial to lateral at 50% and 62.5% (Fujisawa point), respectively. The variations of the force relative to the deformation, strains, and displacements resulting from the progressive loading of the tibial plateau are studied. The research results reveal that using the Fujisawa point is better for conducting the correction not only for medical reasons, but also from a mechanical point of view.

Dynamic effect of three locking plates fixated to humeral fracture based on multibody musculoskeletal model

OBJECTIVE

This study attempts to analyse the biomechanical effect of internal fixation (plated in parallel or plated vertically) on the basis of distal humeral fractures on musculoskeletal multibody dynamics using AnyBody in Finite Element Method.

METHOD

Humeral 3D models were reconstructed by MIMICS after volunteers' CT image input in *.dicom format, and processed by Geomagic Studio for surfaces, while locking plates and screws were then designed by Pro-E. A humeral model of T-type fracture was created and assembled in Hypermesh, to integrate fixtures (e.g., MPL/PML/ML), to grid the mesh and then assign materials. A musculoskeletal model of the upper limb was established by AnyBody to simulate elbow flexion and extension. They were finally imported to Abaqus for boundary conditions and dynamic analysis.

RESULT

In terms of Von Mises stress, its maximum increased and then decreased gradually during the joint motion, but p > 0.05 in SPSS suggests no significant difference for all three fixtures. In terms of displacement, when the elbow was at 90°, each motional pattern reached its peak as follows: ML180° = 0.28 mm, MPL90° = 0.49 mm & PML90° = 0.54 mm during flexion; ML180° = 0.073 mm, MPL90° = 0.10 mm & PML90° = 0.12 mm during extension. p < 0.05 suggests a significant difference for the displacements of all three fixations. p = 0.007 < 0.01667 suggests the significant difference between the two fixations, for example, PML90° and ML180°, indicating that the peak displacement of ML180° is less than that of PML90°.

CONCLUSION

After generally analysed in musculoskeletal dynamics, the biomechanical property of the fixtures was presented as follows: the displacement of the parallel plate was less than that of the vertical, and the parallel plate may optimise the clinical reduction anatomically.

Dynamics analysis for flexion and extension of elbow joint motion based on musculoskeletal model of Anybody

PURPOSE

Little is known about how biomechanics governs the biological nature for humeral motion dynamically. Elbow motion ought to be investigated based on a musculoskeletal model and evidence the physiologic principle of upper limbs.

METHOD

A humeral model was reconstructed by MIMICS after CT images input in *.dicom format, it was processed by Geomagic Studio for Surfaces, then gridded mesh and assigned materials by Hypermesh. On the other hand, a musculoskeletal model was built by AnyBody, physical motions were then simulated to export boundary condition and myodynamia during flexion and extension. Finally, all the humeral model and boundary were imported to Abaqus for finite element analysis.

RESULT

During the simulative motion of flexion, the primary muscles are brachii biceps, brachialis anticus and teretipronator, their myodynamia increased and then decreased gradually, and reached its peak value at 30°; During extension, the main muscles are triceps brachii and brachialis anticus, their myodynamia increased and then decreased gradually too, and reached peak at 50°; In these two cases, their strain and displacement distributed at the middle of humerus.

CONCLUSION

AnyBody is a novel modelling system to simulate physical motion, for example flexion and extension. Biceps brachii and brachialis anticus are functional for flexion, and triceps brachii plays a key role in extension critically. This simulation confirms the physiologic rule for sport event, humeral fixation and postoperative healing with clinical significance that minimizing joint forces from injury onset may promote pain-free ways.

Effect of Rotator Cuff Deficiencies on Muscle Forces and Glenohumeral Contact Force After Anatomic Total Shoulder Arthroplasty Using Musculoskeletal Multibody Dynamics Simulation

Anatomic total shoulder arthroplasty (ATSA) is widely used to treat the diseases of the glenohumeral (GH) joint. However, the incidence of rotator cuff tears after ATSA increases during follow-up. The effects of rotator cuff deficiencies after ATSA on the biomechanics of the GH joint are to be investigated. In this study, a musculoskeletal multibody dynamics model of ATSA was established using a force-dependent kinematics (FDK) method. The biomechanical effects were predicted during arm abduction under different rotator cuff deficiencies. The deltoid forces were increased under the rotator cuff deficiencies, the maximum deltoid forces were increased by 36% under the subscapularis deficiency and by 53% under the supraspinatus, infraspinatus, subscapularis, and teres minor deficiencies. The maximum GH contact forces were decreased by 11.3% under supraspinatus and infraspinatus deficiencies but increased by 24.8% under subscapularis deficiency. The maximum subscapularis force was decreased by 17% under only infraspinatus tear during arm abduction. The results suggested that the changes in the biomechanics of the GH joint induced by rotator cuff deficiencies after ATSA increase the deltoid muscle energy expenditure and joint instability, which result in postoperative less satisfactory clinical outcomes. The changes in rotator cuff muscle forces deserve more attention for understanding the evolution of rotator cuff tear after ATSA.

Development and validation of an in-silico virtual testing rig for analyzing total knee arthroplasty performance during passive deep flexion: A feasibility study

The use of in-silico finite element (FE) models has become more common in orthopedic applications and in the design of biomedical devices, since they can provide results comparable to in vitro experiments while maintaining lower cost. The main downside of this kind of analysis is the high computing time, as it can reach hours or even days to complete; this limitation makes it then not suitable for time-sensitive applications, such as probabilistic analyses or helping clinicians in surgical pre-planning or intra-operative setting. In-silico multibody (MB) simulations, on the other hand, are significantly faster than FE simulations (considering each component of the model as a rigid body); although deformability of each model component is a necessary feature in some applications (e.g. simulation of implant-bone micromotions), several outputs of interest in orthopedic applications, such as implant kinematics and contact forces, do not require a fully deformable model. Therefore, this feasibility study aimed to develop a MB model of a human knee joint implanted with a Total Knee Arthroplasty; a 10 second flexion movement up to 105° was then simulated and the results compared with validated FE models results (under similar boundary conditions) from literature, to perform a preliminary validation in terms of kinematic and kinetic results between the two methods. The agreement and relatively low computing time obtaining with this approach represent a promising starting point for subsequent studies and applications of such techniques in the clinical field.

Predicting ground reaction and tibiotalar contact forces after total ankle arthroplasty during walking

The motion capture and force plates data are essential inputs for musculoskeletal multibody dynamics models to predict in vivo tibiotalar contact forces. However, it could be almost impossible to obtain valid force plates data in old patients undergoing total ankle arthroplasty under some circumstances, such as smaller gait strides and inconsistent walking speeds during gait analysis. To remove the dependence of force plates, this study has established a patient-specific musculoskeletal multibody dynamics model with total ankle arthroplasty by combining a foot-ground contact model based on elastic contact elements. And the established model could predict ground reaction forces, ground reaction moments and tibiotalar contact forces simultaneously. Three patients' motion capture and force plates data during their normal walking were used to establish the patient-specific musculoskeletal models and evaluate the predicted ground reaction forces and ground reaction moments. Reasonable accuracies were achieved for the predicted and measured ground reaction forces and ground reaction moments. The predicted tibiotalar contact forces for all patients using the foot-ground contact model had good consistency with those using force plates data. These findings suggested that the foot-ground contact model could take the place of the force plates data for predicting the tibiotalar contact forces in other total ankle arthroplasty patients, thus providing a simplified and valid platform for further study of the patient-specific prosthetic designs and clinical problems of total ankle arthroplasty in the absence of force plates data.

Patellofemoral Mechanics: a Review of Pathomechanics and Research Approaches

PURPOSE OF REVIEW

The patellofemoral joint is a complicated articulation of the patella and femur that is prone to pathologies. The purpose of this review is to report on the current methods of investigating patellofemoral mechanics, factors that affect joint function, and future directions in patellofemoral joint research with emerging technologies and techniques.

RECENT FINDINGS

While previous hypotheses have suggested that the patella is only a moment arm extender, recent literature has suggested that the patella influences the control of knee moments and forces acting on the tibia as well as contributes to various aspects of patellar function with minimal neural input. With advancements in simulating a six-degrees-of-freedom patellofemoral joint, we have gained a better understanding of patella motion and have shown that geometry and muscle activations impact patella mechanics. Research into influences on patella mechanics from other joints such as the hip and foot has become more prevalent. In this review, we report current in vivo, in vitro, and in silico approaches to studying the patellofemoral joint. Kinematic and anatomical factors that affect patellofemoral joint function such as patella alta and tilt or bone morphology and ligaments are discussed. Moving forward, we suggest that advanced in vivo dynamic imaging methods coupled to musculoskeletal simulation will provide further understanding of patellofemoral pathomechanics and allow engineers and clinicians to design interventions to mitigate or prevent patellofemoral pathologies.

Immediate Effects of Medially Posted Insoles on Lower Limb Joint Contact Forces in Adult Acquired Flatfoot: A Pilot Study

Flatfoot is linked to secondary lower limb joint problems, such as patellofemoral pain. This study aimed to investigate the influence of medial posting insoles on the joint mechanics of the lower extremity in adults with flatfoot. Gait analysis was performed on fifteen young adults with flatfoot under two conditions: walking with shoes and foot orthoses (WSFO), and walking with shoes (WS) in random order. The data collected by a vicon system were used to drive the musculoskeletal model to estimate the hip, patellofemoral, ankle, medial and lateral tibiofemoral joint contact forces. The joint contact forces in WSFO and WS conditions were compared. Compared to the WS group, the second peak patellofemoral contact force (p < 0.05) and the peak ankle contact force (p < 0.05) were significantly lower in the WSFO group by 10.2% and 6.8%, respectively. The foot orthosis significantly reduced the peak ankle eversion angle (p < 0.05) and ankle eversion moment (p < 0.05); however, the peak knee adduction moment increased (p < 0.05). The reduction in the patellofemoral joint force and ankle contact force could potentially inhibit flatfoot-induced lower limb joint problems, despite a greater knee adduction moment.

Evaluation and validation of 2D biomechanical models of the knee for radiograph-based preoperative planning in total knee arthroplasty

Thorough preoperative planning in total knee arthroplasty is essential to reduce implant failure by proper implant sizing and alignment. The “gold standard” in conventional preoperative planning is based on anterior-posterior long-leg radiographs. However, the coronal component alignment is still an open discussion in literature, since studies have reported contradictory outcomes on survivorship, indicating that optimal individual alignment goals still need to be defined. Two-dimensional biomechanical models of the knee have the potential to predict joint forces and, therefore, objectify therapy planning. Previously published two-dimensional biomechanical models were evaluated and validated for the first time in this study by comparison of model predictions to corresponding in vivo measurements obtained from telemetric implants for a one- and two-leg stance. Model input parameters were acquired from weight-bearing anterior-posterior long-leg radiographs and statistical assumptions for patient-specific model adaptation. The overall time from initialization to load prediction was in the range of 7–8 minutes per patient for all models. However, no model could accurately predict the correct trend of knee joint forces over patients. Two dimensional biomechanical models of the knee have the potential to improve preoperative planning in total knee arthroplasty by providing additional individual biomechanical information to the surgeon. Although integration into the clinical workflow might be performed with acceptable costs, the models’ accuracy is insufficient for the moment. Future work is needed for model optimization and more sophisticated modelling approaches.

The anterolateral ligament is a secondary stabilizer in the knee joint: A validated computational model of the biomechanical effects of a deficient anterior cruciate ligament and anterolateral ligament on knee joint kinematics

Objectives

The aim of this study was to investigate the biomechanical effect of the anterolateral ligament (ALL), anterior cruciate ligament (ACL), or both ALL and ACL on kinematics under dynamic loading conditions using dynamic simulation subject-specific knee models.

Methods

Five subject-specific musculoskeletal models were validated with computationally predicted muscle activation, electromyography data, and previous experimental data to analyze effects of the ALL and ACL on knee kinematics under gait and squat loading conditions.

Results

Anterior translation (AT) significantly increased with deficiency of the ACL, ALL, or both structures under gait cycle loading. Internal rotation (IR) significantly increased with deficiency of both the ACL and ALL under gait and squat loading conditions. However, the deficiency of ALL was not significant in the increase of AT, but it was significant in the increase of IR under the squat loading condition.

Conclusion

The results of this study confirm that the ALL is an important lateral knee structure for knee joint stability. The ALL is a secondary stabilizer relative to the ACL under simulated gait and squat loading conditions.Cite this article: Bone Joint Res 2019;8:509-517.

Musculoskeletal modeling of total ankle arthroplasty using force-dependent kinematics for predicting in vivo joint mechanics

In vivo load and motion in the ankle joint play a key role in the understanding of the failure mechanism and function outcomes of total ankle arthroplasty. However, a thorough understanding of the biomechanics of the ankle joint in daily activities is lacking. The objective of this study was to develop a novel lower extremity musculoskeletal multibody dynamics model with total ankle arthroplasty considering the 6 degrees of freedom of the ankle joint motions and the deformable contact mechanics of the implant, based on force-dependent kinematics method. A patient who underwent total ankle arthroplasty surgery was considered. The walking gait data of the patient was measured in a gait laboratory and used as the input for the patient-specific musculoskeletal modeling. The predictions from the musculoskeletal model of total ankle arthroplasty included dorsiflexion-plantar flexion, inversion-eversion, internal-external rotation, anterior-posterior translation, inferior-superior translation, and medial-lateral translation of the tibiotalar joint, the ankle contact forces, the muscle activations, and the ligament forces. The magnitudes and tendencies of the predicted results were all within reasonable ranges, as compared with the data available in the literature. The predicted peak total ankle contact force was 6.55 body weight. In addition, the peak contact forces of the lateral and medial compartments were 4.22 body weight and 2.59 body weight, respectively. This study provides a potential new platform for the design of a better ankle prosthesis, the improvement of the operation techniques of the clinicians, and the accelerated postoperative recovery of the patients.

High Tibial Osteotomy: Review of Techniques and Biomechanics

High tibial osteotomy becomes increasingly important in the treatment of cartilage damage or osteoarthritis of the medial compartment with concurrent varus deformity. HTO produces a postoperative valgus limb alignment with shifting the load-bearing axis of the lower limb laterally. However, maximizing procedural success and postoperative knee function still possess many difficulties. The key to improve the postoperative satisfaction and long-term survival is the understanding of the vital biomechanics of HTO in essence. This review article discussed the alignment principles, surgical technique, and fixation plate of HTO as well as the postoperative gait, musculoskeletal dynamics, and contact mechanics of the knee joint. We aimed to highlight the recent findings and progresses on the biomechanics of HTO. The biomechanical studies on HTO are still insufficient in the areas of gait analysis, joint kinematics, and joint contact mechanics. Combining musculoskeletal dynamics modelling and finite element analysis will help comprehensively understand in vivo patient-specific biomechanics after HTO.

Insert conformity variation affects kinematics and wear performance of total knee replacements

BACKGROUND

The insert conformity is a critical factor for successful total knee replacement which must be considered in design of the implant. However, the effects of conformity on knee kinematics and wear under physiological environment are often neglected in previous studies. The present study involved evaluating the biomechanics and wear performance with regard to different insert conformity in total knee replacement.

METHODS

Different tibial inserts with different sagittal and coronal conformity levels were created and analyzed using a previously developed wear prediction framework, coupling a patient-specific musculoskeletal multibody dynamics simulation, finite element and wear analysis. The contact mechanics, kinematics, and wear performance were compared during 10 million cycles of wear simulation.

FINDINGS

The findings revealed that the knee kinematics was affected by sagittal conformity design variables, which further influenced the wear of insert bearing surface. Additionally, kinematics and wear of artificial knee joint were much more sensitive to sagittal than coronal conformity of tibial insert. The lower sagittal conformity designs had lower wear rates, worn area and contact area. In turn, the wear of insert bearing surface also changed insert conformity, and further impacted on knee kinematics.

INTERPRETATION

The present study indicated that the sagittal conformity design of insert surface played a crucial role to improve contact mechanics and kinematics performance and minimize wear of total knee replacement. The optimization of insert conformity should be considered carefully in implant design and surgical procedures.

Biomechanical evaluation of the influence of posterolateral corner structures on cruciate ligaments forces during simulated gait and squatting

Posterolateral corner (PLC) structures of the knee joint comprise complex anatomical soft tissues that support static and dynamic functional movements of the knee. Most previous studies analyzed posterolateral stability in vitro under static loading conditions. This study aimed to evaluate the contributions of the lateral (fibular) collateral ligament (LCL), popliteofibular ligament (PFL), and popliteus tendon (PT) to cruciate ligament forces under simulated dynamic loading conditions by using selective individual resection. We combined medical imaging and motion capture of healthy subjects (four males and one female) to develop subject-specific knee models that simulated the 12 degrees of freedom of tibiofemoral and patellofemoral joint behaviors. These computational models were validated by comparing electromyographic (EMG) data with muscle activation data and were based on previous experimental studies. A rigid multi-body dynamics simulation using a lower extremity musculoskeletal model was performed to incorporate intact and selective resection of ligaments, based on a novel force-dependent kinematics method, during gait (walking) and squatting. Deficiency of the PLC structures resulted in increased loading on the posterior cruciate ligament and anterior cruciate ligament. Among PLC structures, the PT is the most influential on cruciate ligament forces under dynamic loading conditions.

Computational analysis of knee joint stability following total knee arthroplasty

The overall objective of this study was to introduce knee joint power as a potential measure to investigate knee joint stability following total knee arthroplasty (TKA). Specific aims were to investigate whether weakened knee joint stabilizers cause abnormal kinematics and how it influences the knee joint kinetic (i.e., power) in response to perturbation. Patient-specific musculoskeletal models were simulated with experimental gait data from six TKA patients (baseline models). Muscle strength and ligament force parameter were reduced by up to 30% to simulate weak knee joint stabilizers (weak models). Two different muscle recruitment criteria were tested to examine whether altered muscle recruitment pattern can mask the influence of weakened stabilizers on the knee joint kinematics and kinetics. Level-walking knee joint kinematics and kinetics were calculated though force-dependent kinematic and inverse dynamic analyses. Bode analysis was then recruited to estimate the knee joint power in response to a simulated perturbation. Weak models resulted in larger anterior-posterior (A-P) displacement and internal-external (I-E) rotation compared to baseline (I-E: 18.4 ± 8.5 vs. 11.6 ± 5.7 (deg), A-P: 9.7 ± 5.6 vs. 5.5 ± 4.1 (mm)). Changes in muscle recruitment criterion however altered the results such that A-P and I-E were not notably different from baseline models. In response to the simulated perturbation, weak models versus baseline models generated a delayed power response with unbounded magnitudes. Perturbed power behavior of the knee remained unaltered regardless of the muscle recruitment criteria. In conclusion, impairment at the knee joint stabilizers may or may not lead to excessive joint motions but it notably affects the knee joint power in response to a perturbation. Whether perturbed knee joint power is associated with the patient-reported outcome requires further investigation.

The role of menisci in knee contact mechanics and secondary kinematics during human walking

BACKGROUND

Meniscectomy is likely to result in an increase of joint loading on articular cartilage and initiates cartilage damages. However, the induced alterations in knee contact mechanics and secondary kinematics after the meniscal removal are still unclear during a walking gait. In this study, the role of menisci on the knee contact mechanics and secondary kinematics were investigated during a walking simulation.

METHODS

Two natural knee models with or without menisci were established using a musculoskeletal multibody dynamics framework. Walking simulation were performed to qualify the knee contact forces and secondary kinematics, and ligament forces in these models.

FINDINGS

After the meniscal removal, the redistributed contact forces on the medial tibial cartilage increased by twofold, while the contact area remained almost the same. The subsequent increase of contact pressure suggested potential cartilage damages. In terms of the kinematic alterations, the tibia moved more posteriorly and internally with respect to the femur. And, the displacement in the medial-lateral direction reversed. In addition, the sharp force increase in the anterior cruciate ligament explained the reason why meniscal and anterior cruciate ligament injuries always happened concurrently. And, the anterior lateral ligament may act as the stabilizer in the tibial posterior displacement and varus rotation.

INTERPRETATION

This study shows that menisci served as the joint load distribution, and also as the kinematics constraints in the posterior and medial-lateral direction, which is beneficial to the rehabilitation plan-making of meniscal injuries.

Biomechanics and wear comparison between mechanical and kinematic alignments in total knee arthroplasty

The uses of mechanical and kinematic alignments in total knee arthroplasty are under debate in recent clinical investigations. In this study, the differences in short-term biomechanics and long-term wear volume between mechanical and kinematic alignments in total knee arthroplasty were investigated, based on a subject-specific musculoskeletal multi-body dynamics model during walking gait simulation. An increase of 8.2% in the peak tibiofemoral medial contact force, a posterior contact translation by maximum 4.7 mm and a decrease of 5.5% in the wear volume after a 10-million-cycle simulation were predicted in the kinematic alignment, compared with the mechanical alignment. Nevertheless, the tibiofemoral contact mechanics, the range of motions and the long-term wear were not markedly different between mechanical and kinematic alignments. Furthermore, the mechanical alignment with a posterior tibial slope similar to that under the kinematic alignment was found to produce similar anterior-posterior translation and the range of motion, and an approximate wear volume, compared with the kinematic alignment. The ligament forces under the kinematic alignment were influenced markedly by as much as 25%, 50% and 77% for the medial collateral ligament, lateral collateral ligament and posterior cruciate ligament forces, respectively. And, a maximum increase of 40% for patellofemoral contact force was predicted under the kinematic alignment. These findings suggest that the kinematic alignment is an alternative alignment principle but no marked advantages in biomechanics and wear to the mechanical alignment. The adverse effects of the kinematic alignment on patella loading and soft tissue forces should be noticed.

Prediction of in vivo lower cervical spinal loading using musculoskeletal multi-body dynamics model during the head flexion/extension, lateral bending and axial rotation

Cervical spine diseases lead to a heavy economic burden to the individuals and societies. Moreover, frequent post-operative complications mean a higher risk of neck pain and revision. At present, controversy still exists for the etiology of spinal diseases and their associated complications. Knowledge of in vivo cervical spinal loading pattern is proposed to be the key to answer these questions. However, direct acquisition of in vivo cervical spinal loading remains challenging. In this study, a previously developed cervical spine musculoskeletal multi-body dynamics model was utilized for spinal loading prediction. The in vivo dynamic segmental contributions to head motion and the out-of-plane coupled motion were both taken into account. First, model validation and sensitivity analysis of different segmental contributions to head motion were performed. For model validation, the predicted intervertebral disk compressive forces were converted into the intradiskal pressures and compared with the published experimental measurements. Significant correlations were found between the predicted values and the experimental results. Thus, the reliability and capability of the cervical spine model was ensured. Meanwhile, the sensitivity analysis indicated that cervical spinal loading is sensitive to different segmental contributions to head motion. Second, the compressive, shear and facet joint forces at C3–C6 disk levels were predicted, during the head flexion/extension, lateral bending and axial rotation. Under the head flexion/extension movement, asymmetric loading patterns of the intervertebral disk were obtained. In comparison, symmetrical typed loading patterns were found for the head lateral bending and axial rotation movements. However, the shear forces were dramatically increased during the head excessive extension and lateral bending. Besides, a nonlinear correlation was seen between the facet joint force and the angular displacement. In conclusion, dynamic cervical spinal loading was both intervertebral disk angle-dependent and level-dependent. Cervical spine musculoskeletal multi-body dynamics model provides an attempt to comprehend the in vivo biomechanical surrounding of the human head-neck system.

Biomechanical influence of deficient posterolateral corner structures on knee joint kinematics: A computational study

The posterolateral corner (PLC) structures including the popliteofibular ligament (PFL), popliteus tendon (PT) and lateral collateral ligament (LCL) are important soft tissues for posterior translational, external rotational, and varus angulation knee joint instabilities. The purpose of this study was to determine the effects of deficient PLC structures on the kinematics of the knee joint under gait and squat loading conditions. We developed subject-specific computational models with full 12-degree-of-freedom tibiofemoral and patellofemoral joints for four male subjects and one female subject. The subject-specific knee joint models were validated with computationally predicted muscle activation, electromyography data, and experimental data from previous study. According to our results, deficiency of the PFL did not significantly influence knee joint kinematics compared to an intact model under gait loading conditions. Compared with an intact model under gait and squat loading conditions, deficiency of the PT led to significant increases in external rotation and posterior translation, while LCL deficiency increased varus angulation. Deficiency of all PLC structures led to the greatest increases in external rotation, varus angulation, and posterior translation. These results suggest that the PT is an important structure for external rotation and posterior translation, while the LCL is important for varus angulation under dynamic loading conditions. © 2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 9999:1-8, 2018.

Prediction of Cervical Spinal Joint Loading and Secondary Motion Using a Musculoskeletal Multibody Dynamics Model Via Force-Dependent Kinematics Approach

STUDY DESIGN

A cervical spine biomechanical investigation using multibody dynamics.

OBJECTIVE

To develop a comprehensive cervical spine multibody dynamics model incorporated with the force-dependent kinematics (FDK) approach, and to study the influence of soft tissue deformation on the joint loading prediction.

SUMMARY OF BACKGROUND DATA

Musculoskeletal multibody dynamics models have been widely used to analyze joint loading. Current cervical spine musculoskeletal models, however, neglect the joint internal motion caused by soft tissue deformation. A novel FDK approach is introduced, which can predict joint internal motion and spinal joint loading simultaneously.

METHODS

A comprehensive cervical spine musculoskeletal model with the posterior facet joints and essential ligaments was developed. To quantify the influence of soft tissue structures on joint loading prediction, four different models with different features were created. These newly developed models were validated, under flexion-extension movement. The predicted intervertebral disc loads (from C3-C4 to C5-C6) were compared with the published cadaveric experimental results. Moreover, the predicted facet joint forces, ligament forces, and anterior-posterior translations of instantaneous centers of rotation were also studied.

RESULTS

The obtained intervertebral disc loads were varied among different models. Model 3 provided the closest prediction of joint loading to the experimental results. Moreover, the facet joint and ligament forces were in similar range of magnitude as literature findings. The predicted instantaneous centers of rotation translational changes were in accordance with the in vivo kinematics observation.

CONCLUSION

In the present study, a validated cervical spine musculoskeletal model was developed, using multibody dynamics and FDK approach. It can simulate the function of musculature and consider joint internal motion, and thus provides more reliable joint loading prediction. This newly developed cervical model can be used as an efficient tool to study the biomechanical behaviors of human cervical spine, and to understand the fundamental pathologies of spinal pains.

LEVEL OF EVIDENCE

N /A.

Parametric study of patient-specific femoral locking plates based on a combined musculoskeletal multibody dynamics and finite element modeling

A combined musculoskeletal multibody dynamics and finite element modeling was performed to investigate the effects of design parameters on the fracture-healing efficiency and the mechanical property of a patient-specific anatomically adjusted femoral locking plate. Specifically, the screw type, the thickness and material of the locking plate, the gap between two femoral fragments (fracture gap) and the distance between bone and plate (interface gap) were evaluated during a human walking. We found that the patient-specific locking plate possessed greater mechanical strength and more efficient fracture healing than the corresponding traditional plate. An optimal patient-specific femoral locking plate would consist of bicortical locking screws, Ti-6Al-4V material and 4.75-mm plate thickness with a fracture gap of 2 mm and an interface gap of 1 mm. The developed patient-specific femoral locking plate based on the patient-specific musculoskeletal mechanical environment was more beneficial to fracture rehabilitation and healing. The patient-specific design method provides an effective research platform for designing and optimizing the patient-specific femoral locking plate under realistic in vivo walking conditions, which can be extended to the design of other implants as well as to other physiological loading conditions related to various daily activities.

Load application for the contact mechanics analysis and wear prediction of total knee replacement

Tibiofemoral contact forces in total knee replacement have been measured at the medial and lateral sites respectively using an instrumented prosthesis, and predicted from musculoskeletal multibody dynamics models with a reasonable accuracy. However, it is uncommon that the medial and lateral forces are applied separately to replace a total axial load according to the ISO standard in the majority of current finite element analyses. In this study, we quantified the different effects of applying the medial and lateral loads separately versus the traditional total axial load application on contact mechanics and wear prediction of a patient-specific knee prosthesis. The load application position played an important role under the medial-lateral load application. The loading set which produced the closest load distribution to the multibody dynamics model was used to predict the contact mechanics and wear for the prosthesis and compared with the total axial load application. The medial-lateral load distribution using the present method was found to be closer to the multibody dynamics prediction than the traditional total axial load application, and the maximum contact pressure and contact area were consistent with the corresponding load variation. The predicted total volumetric wear rate and area were similar between the two load applications. However, the split of the predicted wear volumes on the medial and the lateral sides was different. The lateral volumetric wear rate was 31.46% smaller than the medial from the traditional load application prediction, while from the medial-lateral load application, the lateral side was only 11.8% smaller than the medial. The medial-lateral load application could provide a new and more accurate method of load application for patient-specific preclinical contact mechanics and wear prediction of knee implants.

Finite Element Analysis of the Biomechanical Effects of 3 Posterolateral Corner Reconstruction Techniques for the Knee Joint

PURPOSE

To compare the forces exerted on the cruciate ligaments and the contact stresses on the tibiofemoral (TF) and patellofemoral (PF) joints with respect to 3 different tibial- and fibular-based posterolateral corner (PLC) reconstructions under dynamic loading conditions.

METHODS

A subject-specific finite element knee model was developed by using 3-dimensional anatomic data from motion captures in gait and squat activities, including in vivo knee joint kinematics and muscle forces for the single subject. Cruciate ligament forces and contact stresses on the TF and PF joints under 3 PLC reconstruction techniques (tibial-based, TBR; modified fibular-based, mFBR; conventional fibular-based, cFBR) and PLC-deficient models were compared with those of the intact model in gait and squat loading conditions.

RESULTS

The cruciate ligament forces in the 3 surgical models differed from those in the intact model. The greatest differences in ligament forces from the intact model were found in the cFBR model, whereas there were no remarkable differences between the TBR and mFBR models in both gait and squat loading conditions. Contact stresses on the lateral TF and PF joints of the 3 surgical models were greater than those of the intact model under the squat loading condition.

CONCLUSIONS

The biomechanical effects achieved using the anatomic reconstruction technique were found to be improved compared with those using nonanatomic reconstruction techniques. However, the ligament forces and contact stresses under normal conditions could not be restored through any of the 3 techniques.

CLINICAL RELEVANCE

Anatomic TBR and FBR for grade III PLC injuries could restore better biomechanics in the knee joint compared with nonanatomic reconstruction. However, discrepancy with the normal condition requires further modification of surgical techniques.

The effects of posterior cruciate ligament deficiency on posterolateral corner structures under gait- and squat-loading conditions: A computational knee model

Objectives

The aim of the current study was to analyse the effects of posterior cruciate ligament (PCL) deficiency on forces of the posterolateral corner structure and on tibiofemoral (TF) and patellofemoral (PF) contact force under dynamic-loading conditions.

Methods

A subject-specific knee model was validated using a passive flexion experiment, electromyography data, muscle activation, and previous experimental studies. The simulation was performed on the musculoskeletal models with and without PCL deficiency using a novel force-dependent kinematics method under gait- and squat-loading conditions, followed by probabilistic analysis for material uncertain to be considered.

Results

Comparison of predicted passive flexion, posterior drawer kinematics and muscle activation with experimental measurements showed good agreement. Forces of the posterolateral corner structure, and TF and PF contact forces increased with PCL deficiency under gait- and squat-loading conditions. The rate of increase in PF contact force was the greatest during the squat-loading condition. The TF contact forces increased on both medial and lateral compartments during gait-loading conditions. However, during the squat-loading condition, the medial TF contact force tended to increase, while the lateral TF contact forces decreased. The posterolateral corner structure, which showed the greatest increase in force with deficiency of PCL under both gait- and squat-loading conditions, was the popliteus tendon (PT).

Conclusion

PCL deficiency is a factor affecting the variability of force on the PT in dynamic-loading conditions, and it could lead to degeneration of the PF joint.

Cite this article: K-T. Kang, Y-G. Koh, M. Jung, J-H. Nam, J. Son, Y.H. Lee, S-J. Kim, S-H. Kim. The effects of posterior cruciate ligament deficiency on posterolateral corner structures under gait- and squat-loading conditions: A computational knee model. Bone Joint Res 2017;6:31–42. DOI: 10.1302/2046-3758.61.BJR-2016-0184.R1.

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

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

Intra-Articular Knee Contact Force Estimation During Walking Using Force-Reaction Elements and Subject-Specific Joint Model2

Joint contact forces measured with instrumented knee implants have not only revealed general patterns of joint loading but also showed individual variations that could be due to differences in anatomy and joint kinematics. Musculoskeletal human models for dynamic simulation have been utilized to understand body kinetics including joint moments, muscle tension, and knee contact forces. The objectives of this study were to develop a knee contact model which can predict knee contact forces using an inverse dynamics-based optimization solver and to investigate the effect of joint constraints on knee contact force prediction. A knee contact model was developed to include 32 reaction force elements on the surface of a tibial insert of a total knee replacement (TKR), which was embedded in a full-body musculoskeletal model. Various external measurements including motion data and external force data during walking trials of a subject with an instrumented knee implant were provided from the Sixth Grand Challenge Competition to Predict in vivo Knee Loads. Knee contact forces in the medial and lateral portions of the instrumented knee implant were also provided for the same walking trials. A knee contact model with a hinge joint and normal alignment could predict knee contact forces with root mean square errors (RMSEs) of 165 N and 288 N for the medial and lateral portions of the knee, respectively, and coefficients of determination (R2) of 0.70 and −0.63. When the degrees-of-freedom (DOF) of the knee and locations of leg markers were adjusted to account for the valgus lower-limb alignment of the subject, RMSE values improved to 144 N and 179 N, and R2 values improved to 0.77 and 0.37, respectively. The proposed knee contact model with subject-specific joint model could predict in vivo knee contact forces with reasonable accuracy. This model may contribute to the development and improvement of knee arthroplasty.

Effect of component mal-rotation on knee loading in total knee arthroplasty using multi-body dynamics modeling under a simulated walking gait

Mal-rotation of the components in total knee arthorplasty (TKA) is a major cause of postoperative complications, with an increased propensity for implant loosening or wear leading to revision. A musculoskeletal multi-body dynamics model was used to perform a parametric study of the effects of the rotational mal-alignments in TKA on the knee loading under a simulated walking gait. The knee contact forces were found to be more sensitive to variations in the varus-valgus rotation of both the tibial and the femoral components and the internal-external rotation of the femoral component in TKA. The varus-valgus mal-rotation of the tibial or femoral component and the internal-external mal-rotation of the femoral component with a 5° variation were found to affect the peak medial contact force by 17.8-53.1%, the peak lateral contact force by 35.0-88.4% and the peak total contact force by 5.2-18.7%. Our findings support the clinical observations that a greater than 3° internal mal-rotation of the femoral component may lead to unsatisfactory pain levels and a greater than 3° varus mal-rotation of the tibial component may lead to medial bone collapse. These findings determined the quantitative effects of the mal-rotation of the components in TKA on the contact load. The effect of such mal-rotation of the components of TKA on the kinematics would be further addressed in future studies.

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

GOAL

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

METHODS

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

RESULTS

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

CONCLUSION

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

SIGNIFICANCE

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

A Subject-Specific Musculoskeletal Modeling Framework to Predict In Vivo Mechanics of Total Knee Arthroplasty

Musculoskeletal (MS) models should be able to integrate patient-specific MS architecture and undergo thorough validation prior to their introduction into clinical practice. We present a methodology to develop subject-specific models able to simultaneously predict muscle, ligament, and knee joint contact forces along with secondary knee kinematics. The MS architecture of a generic cadaver-based model was scaled using an advanced morphing technique to the subject-specific morphology of a patient implanted with an instrumented total knee arthroplasty (TKA) available in the fifth “grand challenge competition to predict in vivo knee loads” dataset. We implemented two separate knee models, one employing traditional hinge constraints, which was solved using an inverse dynamics technique, and another one using an 11-degree-of-freedom (DOF) representation of the tibiofemoral (TF) and patellofemoral (PF) joints, which was solved using a combined inverse dynamic and quasi-static analysis, called force-dependent kinematics (FDK). TF joint forces for one gait and one right-turn trial and secondary knee kinematics for one unloaded leg-swing trial were predicted and evaluated using experimental data available in the grand challenge dataset. Total compressive TF contact forces were predicted by both hinge and FDK knee models with a root-mean-square error (RMSE) and a coefficient of determination (R2) smaller than 0.3 body weight (BW) and equal to 0.9 in the gait trial simulation and smaller than 0.4 BW and larger than 0.8 in the right-turn trial simulation, respectively. Total, medial, and lateral TF joint contact force predictions were highly similar, regardless of the type of knee model used. Medial (respectively lateral) TF forces were over- (respectively, under-) predicted with a magnitude error of M < 0.2 (respectively > −0.4) in the gait trial, and under- (respectively, over-) predicted with a magnitude error of M > −0.4 (respectively < 0.3) in the right-turn trial. Secondary knee kinematics from the unloaded leg-swing trial were overall better approximated using the FDK model (average Sprague and Geers' combined error C = 0.06) than when using a hinged knee model (C = 0.34). The proposed modeling approach allows detailed subject-specific scaling and personalization and does not contain any nonphysiological parameters. This modeling framework has potential applications in aiding the clinical decision-making in orthopedics procedures and as a tool for virtual implant design.