Subject-specific modeling of muscle force and knee contact in total knee arthroplasty

Kinematic-kinetic compliant acetabular cup positioning based on preoperative motion tracking and musculoskeletal modeling for total hip arthroplasty

The invention of the surgical robot enabled accurate component implantation during total hip arthroplasty (THA). However, a preoperative surgical planning methodology is still lacking to determine the acetabular cup alignment considering the patient-specific hip functions during daily activities such as walking. To simultaneously avoid implant edgeloading and impingement, this study established a kinematic-kinetic compliant (KKC) acetabular cup positioning method based on preoperative gait kinematics measurement and musculoskeletal modeling. Computed tomography images around the hip joint and their biomechanical data during gait, including motion tracking and foot-ground reaction forces, were collected. Using the reconstructed pelvic and femur geometries, the patient-specific hip muscle insertions were located in the lower limb musculoskeletal model via point cloud registration. The designed cup orientation has to be within the patient-specific safe zone to prevent implant impingement, and the optimized value selected based on the time-dependent hip joint reaction force to minimize the risk of edgeloading. As a validation of the proposed musculoskeletal model, the predicted lower limb muscle activations for seven patients were correlated with their surface electromyographic measurements, and the computed hip contact force was also in quantitative agreement with data from the literature. However, the designed cup orientations were not always within the well-known Lewinnek safe zone, highlighting the importance of KKC surgical planning based on patient-specific biomechanical evaluations.

Breaking the Saturation of Sensitivity for Ultrawide Range Flexible Pressure Sensors by Soft-Strain Effect

The flexible pressure sensors with a broad pressure range and unsaturated sensitivity are highly desired in practical applications. However, pressure sensors by piezoresistive effect are always limited by the compressibility of sensing layers, resulting in a theoretically decreasing sensitivity of less than 100%. Here, a unique strategy is proposed that utilizes the strain effect, simultaneously achieving a trade-off between a wider pressure detection range and unsaturated sensitivity. Ascribed to the strain effect of sensing layers induced by interlaced microdomes, the sensors possess an increased sensitivity (5.22-70 MPa-1) over an ultrawide pressure range (45 Pa-4.1 MPa), a high-pressure resolution (5 Pa), fast response/recovery time (30/45 ms), and a robust response under a high-pressure loading of 3.5 MPa for more than 5000 cycles. These superior sensing performances allow the sensor to monitor large pressure. The flexible pressure sensor array can assist doctors in restoring the neutral mechanical axis, tracking knee flexion angles, and extracting gait features. Moreover, the flexible sensing array can be integrated into the joint motion surveillance system to map the balance medial-lateral contact forces on the metal compartments in real time, demonstrating the potential for further development into precise medical human-machine interfaces during total knee replacement surgery.

Simultaneous Evaluation of Tibiofemoral and Patellofemoral Mechanics in Total Knee Arthroplasty: A Combined Experimental and Computational Approach

Contemporary total knee arthroplasty (TKA) has not fully restored natural patellofemoral (P-F) mechanics across the patient population. Previous experimental simulations have been limited in their ability to create dynamic, unconstrained, muscle-driven P-F articulation while simultaneously controlling tibiofemoral (T-F) contact mechanics. The purpose of this study was to develop a novel experimental simulation and validate a corresponding finite element model to evaluate T-F and P-F mechanics. A commercially available wear simulator was retrofitted with custom fixturing to evaluate whole-knee TKA mechanics with varying patella heights during a simulated deep knee bend. A corresponding dynamic finite element model was developed to validate kinematic and kinetic predictions against experimental measurements. Patella alta reduced P-F reaction forces in early and midflexion, corresponding with an increase in T-F forces that indicated an increase in extensor mechanism efficiency. Due to reduced wrapping of the extensor mechanism in deeper flexion for the alta condition, peak P-F forces in flexion increased from 101% to 135% of the applied quadriceps load for the baja and alta conditions, respectively. Strong agreement was observed between the experiment and model predictions with root-mean-square errors (RMSE) for P-F kinematics ranging from 0.8 deg to 3.3 deg and 0.7 mm to 1.4 mm. RMSE for P-F forces ranged from 7.4 N to 53.6 N. By simultaneously controlling dynamic, physiological loading of the T-F and P-F joint, this novel experimental simulation and validated model will be a valuable tool for investigation of future TKA designs and surgical techniques.

Development of Multibundle Virtual Ligaments to Simulate Knee Mechanics After Total Knee Arthroplasty

Preclinical evaluation of total knee arthroplasty (TKA) components is essential to understanding their mechanical behavior and developing strategies for improving joint stability. While preclinical testing of TKA components has been useful in quantifying their effectiveness, such testing can be criticized for lacking clinical relevance, as the important contributions of surrounding soft tissues are either neglected or greatly simplified. The purpose of our study was to develop and determine if subject-specific virtual ligaments reproduce a similar behavior as native ligaments surrounding TKA joints. Six TKA knees were mounted to a motion simulator. Each was subjected to tests of anterior-posterior (AP), internal-external (IE), and varus-valgus (VV) laxity. The forces transmitted through major ligaments were measured using a sequential resection technique. By tuning the measured ligament forces and elongations to a generic nonlinear elastic ligament model, virtual ligaments were designed and used to simulate the soft tissue envelope around isolated TKA components. The average root-mean-square error (RMSE) between the laxity results of TKA joints with native versus virtual ligaments was 3.5 ± 1.8 mm during AP translation, 7.5 ± 4.2 deg during IE rotations, and 2.0 ± 1.2 deg during VV rotations. Interclass correlation coefficients (ICCs) indicated a good level of reliability for AP and IE laxity (0.85 and 0.84). To conclude, the advancement of virtual ligament envelopes as a more realistic representation of soft tissue constraint around TKA joints is a valuable approach for obtaining clinically relevant kinematics when testing TKA components on joint motion simulators.

Restoring Anatomical Features in Primary Total Knee Arthroplasty

Today, the number of people affected by gonarthrosis symptoms is increasing proportionally. Total knee arthroplasty (TKA) is a successful intervention that aims to reduce pain and restore knee function. However, studies have shown that active young patients still have limitations in performing activities such as skiing, golfing, surfing, and dancing. Over the last few years, total knee arthroplasty has undergone significant changes. Most of the modern TKA implants are designed to reproduce the normal biomechanics of the knee joint, mimicking the physiological pattern with greater compliance in the medial compartment between the tibial insert and femoral condyle and less congruence on the lateral side. Unfortunately, functional outcomes are compromised in approximately half of TKA patients. This loss may be caused by the abnormal kinematics and inherent instability of many contemporary implants. The proper alignment of the femoral component during TKA is a crucial step that influences postoperative results. The position of the femoral component in the axial plane is responsible for flexion stability, knee joint kinematics, flexion alignment, and patellar tracking. The main goal when choosing a type of prosthesis is to achieve an adequate recovery that leads to an improvement in mobility and an increase in the efficiency of the quadriceps.

A novel method to accurately recreate in vivo loads and kinematics in computational models of the knee

Despite availability of in vivo knee loads and kinematics data, conventional load- and displacement-controlled configurations still can't accurately predict tibiofemoral loads from kinematics or vice versa. We propose a combined load- and displacement-control method for joint-level simulations of the knee to reliably reproduce in vivo contact mechanics. Prediction errors of the new approach were compared to those of conventional purely load- or displacement-controlled models using in vivo implant loads and kinematics for multiple subjects and activities (CAMS-Knee dataset). Our method reproduced both loads and kinematics more closely than conventional models and thus demonstrates clear advantages for investigating tibiofemoral contact or wear.

Using Finite Element Modeling in Bone Mechanoadaptation

PURPOSE OF THE REVIEW

Bone adapts structure and material properties in response to its mechanical environment, a process called mechanoadpatation. For the past 50 years, finite element modeling has been used to investigate the relationships between bone geometry, material properties, and mechanical loading conditions. This review examines how we use finite element modeling in the context of bone mechanoadpatation.

RECENT FINDINGS

Finite element models estimate complex mechanical stimuli at the tissue and cellular levels, help explain experimental results, and inform the design of loading protocols and prosthetics. FE modeling is a powerful tool to study bone adaptation as it complements experimental approaches. Before using FE models, researchers should determine whether simulation results will provide complementary information to experimental or clinical observations and should establish the level of complexity required. As imaging technics and computational capacity continue increasing, we expect FE models to help in designing treatments of bone pathologies that take advantage of mechanoadaptation of bone.

Using a finite element model of the thumb to study Trapeziometacarpal joint contact during lateral pinch

BACKGROUND

Finite element (FE) analysis is widely used in different fields of orthopaedic surgery, however, its application to the trapeziometacarpal joint has been limited due to the small size, complex biconcave-convex joint geometry, and complex musculature. The goal of this study was to improve upon existing models by creating a muscle-driven FE thumb model and use the model to simulate the biomechanical effect of hand therapy exercises and ligament reconstructive surgeries.

METHODS

Bone and cartilage geometry were based on a CT dataset of a subject performing a static lateral pinch task. A previously validated musculoskeletal model was utilized to extract electromyography (EMG)-driven muscle forces. Five ligaments with biomechanical significance were modeled as springs using literature values and attached according to their anatomical landmarks.

FINDINGS

The biomechanical consequence of various interventions was proxied as a change in the maximum cartilage stress. The result shows tightening the dorsal ligament complex (dorsal radial ligament, dorsal central ligament, posterior oblique ligament) is the most effective, achieving a stress reduction of 4.8%. Five exercises used in hand therapies were modeled, among which thenar eminence strengthening showed the most prominent stress reduction of 4.0%. Four ligament reconstructive surgeries were modeled, with Eaton-Littler reconstruction showed the most significant stress reduction of 25.0%.

INTERPRETATION

Among the routinely utilized treatment options for early thumb osteoarthritis, we found that three methods: dorsal ligament imbrication, thenar eminence exercise, and the Eaton-Littler method may confer biomechanical advantages cartilage loading. These advantages align with the clinically observed favorable outcomes.

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

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

Superimposition of ground reaction force on tibial-plateau supporting diagnostics and post-operative evaluations in high-tibial osteotomy. A novel methodology

BACKGROUND

A fully personalised combination of Gait Analysis (GA), including Ground Reaction Force (GRF), and patient-specific knee joint morphology has not yet been reported. This can provide valuable biomechanical insight in normal and pathological conditions. Abnormal knee varus results in medial knee condylar hyper-compression and osteoarthritis, which can be prevented by restoring proper condylar load distribution via High Tibial Osteotomy (HTO).

RESEARCH QUESTION

This study was aimed at reporting on an original methodology, merging GA, GRF and Computer-Tomography (CT) to depict a patient-specific representation of the knee mechanical condition during locomotion. It was hypothesised that HTO results in a lateralized pattern of GRF with respect to the tibial plateau.

METHODS

Four patients selected for HTO received clinical, radiological and instrumental examinations, pre- and post-operatively at 6-month follow-up. GA was performed during level walking and more demanding motor tasks using a 9-camera motion-capture system, combined with two force platforms, and an established protocol. Additional skin markers were positioned around the tibial-plateau rim. Weight-bearing CT scans of the knee were collected while still wearing these markers. Proximal tibial and marker morphological models were reconstructed. The markers from CT reconstruction were then registered to the corresponding trajectories as tracked by GA data. Resulting registration matrices were used to report GRF vectors on the plane best matching the tibial-plateau model and the intersection paths were calculated.

RESULTS AND SIGNIFICANCE

The registration procedure was successfully executed, with a max registration error of about 3 mm. GRF intersection paths were found medially to the tibial plateau pre-op, and lateralized post-op, thus much closer to the knee centre, as expected after HTO. The exploitation of the present methodology offers personalised quantification of the original mechanical misalignment and of the effect of surgical correction which could enhance diagnostics and planning of HTO as well as other knee treatments.

Kinematics and kinetics comparison of ultra-congruent versus medial-pivot designs for total knee arthroplasty by multibody analysis

Nowadays, several configurations of total knee arthroplasty (TKA) implants are commercially available whose designs resulted from clinical and biomechanical considerations. Previous research activities led to the development of the so-called medial-pivot (MP) design. However, the actual benefits of the MP, with respect to other prosthesis designs, are still not well understood. The present work compares the impact of two insert geometries, namely the ultra-congruent (UC) and medial-pivot (MP), on the biomechanical behaviour of a bicondylar total knee endoprosthesis. For this purpose, a multibody model of a lower limb was created alternatively integrating the two implants having the insert geometry discretized. Joint dynamics and contact pressure distributions were evaluated by simulating a squat motion. Results showed a similar tibial internal rotation range of about 3.5°, but an early rotation occurs for the MP design. Furthermore, the discretization of the insert geometry allowed to efficiently derive the contact pressure distributions, directly within the multibody simulation framework, reporting peak pressure values of 33 MPa and 20 MPa for the UC and MP, respectively. Clinically, the presented findings confirm the possibility, through a MP design, to achieve a more natural joint kinematics, consequently improving the post-operative patient satisfaction and potentially reducing the occurrence of phenomena leading to the insert loosening.

Comparison of Kinematics and Contact Mechanics in Normal Knee and Total Knee Replacements: A Computational Investigation

An objective of total knee replacement (TKR) is to restore the mechanical function of a normal knee. Joint kinematics and contact mechanics performance are two of the primary indices that indicate the success of TKR devices. The aim of this study was to compare the kinematics and contact mechanics of TKR and normal knee joints. An experimentally evaluated finite-element (FE) knee model was developed and used to investigate the performance of four TKR designs (fixed cruciate-retaining (CR), mobile CR, posterior-stabilized (PS), medial pivot design (MP)) and the normal knee joint during a gait cycle. The predicted kinematic results showed that the MP design presented similar kinematics to those of the normal knee joint and did not demonstrate paradoxical motion of the femur. A considerably larger contact area and lower contact pressure were found on the normal knee joint (1315 mm², and 14.8 MPa, respectively) than on the TKRs, which was consistent with the previous in-vivo fluoroscopic investigation. The mobile CR and PS designs exhibited the smallest and greatest contact pressures of the four TKR designs, respectively. The results of the present study help to understand the kinematics and contact mechanics in the TKR during the gait cycle, and provide comprehensive information about the performance of the normal knee joint.

Machine learning for rapid estimation of lower extremity muscle and joint loading during activities of daily living

Joint contact and muscle forces estimated with musculoskeletal modeling techniques offer useful metrics describing movement quality that benefit multiple research and clinical applications. The expensive processing of laboratory data associated with generating these outputs presents challenges to researchers and clinicians, including significant time and expertise requirements that limit the number of subjects typically evaluated. The objective of the current study was to develop and compare machine learning techniques for rapid, data-driven estimation of musculoskeletal metrics from derived gait lab data. OpenSim estimates of patient joint and muscle forces during activities of daily living were simulated using laboratory data from 70 total knee replacement patients and used to develop 4 different machine learning algorithms. Trained machine learning models predicted both trend and magnitude of estimated joint contact (mean correlation coefficients ranging from 0.93 to 0.94 during gait) and muscle forces (mean correlation coefficients ranging from 0.83 to 0.91 during gait) based on anthropometrics, ground reaction forces, and joint angle data. Patient mechanics were accurately predicted by recurrent neural networks, even after removing dependence on key subsets of predictor features. The ability to quickly estimate patient mechanics from derived measurements of movement has the potential to broaden the impact of musculoskeletal modeling by enabling faster assessment in both clinical and research settings.

Total Knee Replacement: Subject-Specific Modeling, Finite Element Analysis, and Evaluation of Dynamic Activities

This study presents a semi-automatic framework to create subject-specific total knee replacement finite element models, which can be used to analyze locomotion patterns and evaluate knee dynamics. In recent years, much scientific attention was attracted to pre-clinical optimization of customized total knee replacement operations through computational modeling to minimize post-operational adverse effects. However, the time-consuming and laborious process of developing a subject-specific finite element model poses an obstacle to the latter. One of this work's main goals is to automate the finite element model development process, which speeds up the proposed framework and makes it viable for practical applications. This pipeline's reliability was ratified by developing and validating a subject-specific total knee replacement model based on the 6th SimTK Grand Challenge data set. The model was validated by analyzing contact pressures on the tibial insert in relation to the patient's gait and analysis of tibial contact forces, which were found to be in accordance with the ones provided by the Grand Challenge data set. Subsequently, a sensitivity analysis was carried out to assess the influence of modeling choices on tibial insert's contact pressures and determine possible uncertainties on the models produced by the framework. Parameters, such as the position of ligament origin points, ligament stiffness, reference strain, and implant-bone alignment were used for the sensitivity study. Notably, it was found that changes in the alignment of the femoral component in reference to the knee bones significantly affect the load distribution at the tibiofemoral joint, with an increase of 206.48% to be observed at contact pressures during 5° internal rotation. Overall, the models produced by this pipeline can be further used to optimize and personalize surgery by evaluating the best surgical parameters in a simulated manner before the actual surgery.

A musculoskeletal modelling approach to explain sit-to-stand difficulties in older people due to changes in muscle recruitment and movement strategies

By 2050 the proportion of over 65s is predicted to be 20% of the population. The consequences of an age-related reduction in muscle mass have not been fully investigated and, therefore, the aim of the present study was to quantify the muscle and joint contact forces using musculoskeletal modelling, during a sit-to-stand activity, to better explain difficulties in performing everyday activities for older people. A sit-to-stand activity with and without the use of arm rests was observed in ninety-five male participants, placed into groups of young (aged 18-35 years), middle-aged (aged 40-60 years) or older adults (aged 65 years and over). Older participants demonstrated significantly lower knee extensor and joint forces than the young when not using arm rests, compensating through elevated hip extensor and ankle plantarflexor muscle activity. The older group were also found to have higher shoulder joint contact forces whilst using arm rests. This tendency to reorganise muscle recruitment to include neighbouring groups or other parts of the body could make everyday activities more susceptible to age-related functional decline. Reductions in leg strength, via age- or atrophy- related means, creates increased reliance on the upper body and may result in further lower limb atrophy through disuse. The eventual decline of upper body function reduces strength reserves, leading to increased vulnerability, dependence on others and risk of institutionalisation.

Simulated hip abductor strengthening reduces peak joint contact forces in patients with total hip arthroplasty

Lower extremity muscle strength training is a focus of rehabilitation following total hip arthroplasty (THA). Strength of the hip abductor muscle group is a predictor of overall function following THA. The purpose of this study was to investigate the effects of hip abductor strengthening following rehabilitation on joint contact forces (JCFs) in the lower extremity and low back during a high demand step down task. Five THA patients performed lower extremity maximum isometric strength tests and a stair descent task. Patient-specific musculoskeletal models were created in OpenSim and maximum isometric strength parameters were scaled to reproduce measured pre-operative joint torques. A pre-operative forward dynamic simulation of each patient performing the stair descent was constructed using their corresponding patient-specific model to predict JCFs at the ankle, knee, hip, and low back. The hip abductor muscles were strengthened with clinically supported increases (0-30%) above pre-operative values in a probabilistic framework to predict the effects on peak JCFs (99% confidence bounds). Simulated hip abductor strengthening resulted in lower peak JCFs relative to pre-operative for all five patients at the hip (18.9-23.8 ± 16.5%) and knee (20.5-23.8 ± 11.2%). Four of the five patients had reductions at the ankle (7.1-8.5 ± 11.3%) and low back (3.5-7.0 ± 5.3%) with one patient demonstrating no change. The reduction in JCF at the hip joint and at joints other than the hip with hip abductor strengthening demonstrates the dynamic and mechanical interdependencies of the knee, hip and spine that can be targeted in early THA rehabilitation to improve overall patient function.

A lower extremity model for muscle-driven simulation of activity using explicit finite element modeling

A key strength of computational modeling is that it can provide estimates of muscle, ligament, and joint loads, stresses, and strains through non-invasive means. However, simulations that can predict the forces in the muscles during activity while maintaining sufficient complexity to realistically represent the muscles and joint structures can be computationally challenging. For this reason, the current state of the art is to apply separate rigid-body dynamic and finite-element (FE) analyses in series. However, the use of two or more disconnected models often fails to capture key interactions between the joint-level and whole-body scales. Single framework MSFE models have the potential to overcome the limitations associated with disconnected models in series. The objectives of the current study were to create a multi-scale FE model of the human lower extremity that combines optimization, dynamic muscle modeling, and structural FE analysis in a single framework and to apply this framework to evaluate the mechanics of healthy knee specimens during two activities. Two subject-specific FE models (Model 1, Model 2) of the lower extremity were developed in ABAQUS/Explicit including detailed representations of the muscles. Muscle forces, knee joint loading, and articular contact were calculated for two activities using an inverse dynamics approach and static optimization. Quadriceps muscle forces peaked at the onset of chair rise (2174 N, 1962 N) and in early stance phase (510 N, 525 N), while gait saw peak forces in the hamstrings (851 N, 868 N) in midstance. Joint forces were similar in magnitude to available telemetric patient data. This study demonstrates the feasibility of detailed quasi-static, muscle-driven simulations in an FE framework.

A computationally efficient strategy to estimate muscle forces in a finite element musculoskeletal model of the lower limb

Concurrent multiscale simulation strategies are required in computational biomechanics to study the interdependence between body scales. However, detailed finite element models rarely include muscle recruitment due to the computational burden of both the finite element method and the optimization strategies widely used to estimate muscle forces. The aim of this study was twofold: first, to develop a computationally efficient muscle force prediction strategy based on proportional-integral-derivative (PID) controllers to track gait and chair rise experimental joint motion with a finite element musculoskeletal model of the lower limb, including a deformable knee representation with 12 degrees of freedom; and, second, to demonstrate that the inclusion of joint-level deformability affects muscle force estimation by using two different knee models and comparing muscle forces between the two solutions. The PID control strategy tracked experimental hip, knee, and ankle flexion/extension with root mean square errors below 1°, and estimated muscle, contact and ligament forces in good agreement with previous results and electromyography signals. Differences up to 11% and 20% in the vasti and biceps femoris forces, respectively, were observed between the two knee models, which might be attributed to a combination of differing joint contact geometry, ligament behavior, joint kinematics, and muscle moment arms. The tracking strategy developed in this study addressed the inevitable tradeoff between computational cost and model detail in musculoskeletal simulations and can be used with finite element musculoskeletal models to efficiently estimate the interdependence between muscle forces and tissue deformation.

Evaluation of the accuracy of musculoskeletal simulation during squats by means of instrumented knee prostheses

Standard musculoskeletal simulation tools now offer widespread access to internal loading conditions for use in improving rehabilitation concepts or training programmes. However, despite broad reliance on their outcome, the accuracy of such loading estimations, specifically in deep knee flexion, remains generally unknown. The aim of this study was to evaluate the error of tibio-femoral joint contact force (JCF) calculations using musculoskeletal simulation compared to in vivo measured JCFs in subjects with instrumented total knee endoprostheses during squat exercises. Using the early but common "Gait2392_simbody" (OpenSim) scaled musculoskeletal models, tibio-femoral JCFs were calculated in 6 subjects for 5 repetitions of squats. Tibio-femoral JCFs of 0.8-3.2 times bodyweight (BW) were measured. While the musculoskeletal simulations underestimated the measured knee JCFs at low flexion angles, an average error of less than 20% was achieved between approximately 25°-60° knee flexion. With an average error that behaved almost linearly with knee flexion angle, an overestimation of approximately 60% was observed at deep flexion (ca. 80°), with an absolute maximum error of ca. 1.9BW. Our data indicate that loading estimations from early musculoskeletal gait models at both high and low knee joint flexion angles should be interpreted carefully.

Influence of varying stem and metaphyseal sleeve size on the primary stability of cementless revision tibial trays used to reconstruct AORI IIA defects. A simulation study

Traditionally, diaphyseal stems have been utilized to augment the stability of revision total knee replacement (rTKR) implants. More recently metaphyseal augments, such as sleeves, have been introduced to further augment component fixation. The effect of augments such as stems and sleeves have on the primary stability of a rTKR implant is poorly understood, however it has important implications on the complexity, costs and survivorship of the procedure. Finite element analysis was used to investigate the primary stability and strain distribution of various size stems and sleeves used in conjunction with a cementless revision tibial tray. The model was built from computer tomography images of a single healthy tibia obtained from an 81-year-old patient to which an Anderson Orthopaedic Research Institute (AORI) IIA defect was virtually added. The influences of varying body mass index (BMI) and bone modulus were also investigated. Stemless sleeves were found to provided adequate primary implant stability (average implant micro-motion <50 μm) for the studied defect. Addition of a stem did not enhance the primary stability. Furthermore, this study found that varying BMI and bone modulus had a considerable effect on strain distribution but negligible effect on micro-motion in the sleeve area. In conclusion, the addition of diaphyseal stem to a metaphyseal sleeve had little benefit in enhancing the primary stability of tibial trays augmented when simulating reconstructions of AORI IIA tibial defects. Additional studies are required to determine the relative benefit of the diaphyseal stem when using metaphyseal sleeves defects with more extensive bone loss. © 2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:1876-1886, 2018.

In vivo comparison of medialized dome and anatomic patellofemoral geometries using subject-specific computational modeling

Successful outcome following total knee arthroplasty (TKA) with patella resurfacing is partly determined by the restoration of patellofemoral (PF) function and recovery of the quadriceps mechanism. The current study compared two patellar TKA geometries (medialized dome and anatomic) to determine their impact on PF mechanics and quadriceps function. In-vivo, subject-specific patellar mechanics were evaluated using a sequential experimental and modeling approach. First, stereo radiography, marker-based motion capture, and force plate data were collected for TKA patients (10 dome, 10 anatomic) performing a knee extension and lunge. Second, subject-specific, whole-body, musculoskeletal models, including 6 degrees-of-freedom (DOF) knee joint kinematics, were created for each subject and activity to predict quadriceps forces. Last, finite element models of each subject and activity were created to predict PF kinematics, patellar loading, moment arm, and patellar tendon angle. Differences in mechanics between dome and anatomic patients were highlighted during load-bearing (lunge) activity. Anatomic subjects demonstrated greater PF flexion angles (avg. 11 ± 3°) compared to dome subjects during lunge. Similar to the natural knee, contact locations on the patella migrated inferior to superior as the knee flexed in anatomic subjects, but remained relatively superior in dome subjects. Differences in kinematics and contact location likely contributed to altered mechanics with anatomic subjects presenting greater load transfer from the quadriceps to the patellar tendon in deep flexion (>75°), and dome subjects demonstrating larger contact forces during lunge. Although there was substantial patient variability, evaluations of PF mechanics suggested improved quadriceps function and more natural kinematics in the anatomic design. © 2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:1910-1918, 2018.

The interaction of muscle moment arm, knee laxity, and torque in a multi-scale musculoskeletal model of the lower limb

INTRODUCTION

Musculoskeletal modeling allows insight into the interaction of muscle force and knee joint kinematics that cannot be measured in the laboratory. However, musculoskeletal models of the lower extremity commonly use simplified representations of the knee that may limit analyses of the interaction between muscle forces and joint kinematics. The goal of this research was to demonstrate how muscle forces alter knee kinematics and consequently muscle moment arms and joint torque in a musculoskeletal model of the lower limb that includes a deformable representation of the knee.

METHODS

Two musculoskeletal models of the lower limb including specimen-specific articular geometries and ligament deformability at the knee were built in a finite element framework and calibrated to match mean isometric torque data collected from 12 healthy subjects. Muscle moment arms were compared between simulations of passive knee flexion and maximum isometric knee extension and flexion. In addition, isometric torque results were compared with predictions using simplified knee models in which the deformability of the knee was removed and the kinematics at the joint were prescribed for all degrees of freedom.

RESULTS

Peak isometric torque estimated with a deformable knee representation occurred between 45° and 60° in extension, and 45° in flexion. The maximum isometric flexion torques generated by the models with deformable ligaments were 14.6% and 17.9% larger than those generated by the models with prescribed kinematics; by contrast, the maximum isometric extension torques generated by the models were similar. The change in hamstrings moment arms during isometric flexion was greater than that of the quadriceps during isometric extension (a mean RMS difference of 9.8 mm compared to 2.9 mm, respectively).

DISCUSSION

The large changes in the moment arms of the hamstrings, when activated in a model with deformable ligaments, resulted in changes to flexion torque. When simulating human motion, the inclusion of a deformable joint in a multi-scale musculoskeletal finite element model of the lower limb may preserve the realistic interaction of muscle force with knee kinematics and torque.

Comparison of Marker-Based and Stereo Radiography Knee Kinematics in Activities of Daily Living

Movement of the marker positions relative to the body segments obscures in vivo joint level motion. Alternatively, tracking bones from radiography images can provide precise motion of the bones at the knee but is impracticable for measurement of body segment motion. Consequently, researchers have combined marker-based knee flexion with kinematic splines to approximate the translations and rotations of the tibia relative to the femur. Yet, the accuracy of predicting six degree-of-freedom joint kinematics using kinematic splines has not been evaluated. The objectives of this study were to (1) compare knee kinematics measured with a marker-based motion capture system to kinematics acquired with high speed stereo radiography (HSSR) and describe the accuracy of marker-based motion to improve interpretation of results from these methods, and (2) use HSSR to define and evaluate a new set of knee joint kinematic splines based on the in vivo kinematics of a knee extension activity. Simultaneous measurements were recorded from eight healthy subjects using HSSR and marker-based motion capture. The marker positions were applied to three models of the lower extremity to calculate tibiofemoral kinematics and compared to kinematics acquired with HSSR. As demonstrated by normalized RMSE above 1.0, varus-valgus rotation (1.26), medial-lateral (1.26), anterior-posterior (2.03), and superior-inferior translations (4.39) were not accurately measured. Using kinematic splines improved predictions in varus-valgus (0.81) rotation, and medial-lateral (0.73), anterior-posterior (0.69), and superior-inferior (0.49) translations. Using splines to predict tibiofemoral kinematics as a function knee flexion can lead to improved accuracy over marker-based motion capture alone, however this technique was limited in reproducing subject-specific kinematics.

The impact of hip implant alignment on muscle and joint loading during dynamic activities

BACKGROUND

Component alignment is an important consideration in total hip arthroplasty. The impact of changes in alignment on muscle forces and joint contact forces during dynamic tasks are not well understood, and have the potential to influence surgical decision making. The objectives of this study were to assess the impact of femoral head/stem and cup component placement on hip muscle and joint contact forces during tasks of daily living and to identify which alignment parameters have the greatest impact on joint loading.

METHODS

Using a series of strength-calibrated, subject-specific musculoskeletal models of patients performing gait, sit-to-stand and step down tasks, component alignments were perturbed and joint contact and muscle forces evaluated.

FINDINGS

Based on the range of alignments reported clinically, variation in head/stem anteversion-retroversion had the largest impact of any degree of freedom throughout all three tasks; average contact forces 413.5 (319.1) N during gait, 262.7 (256.4) N during sit to stand, and 572.7 (228.1) N during the step down task. The sensitivity of contact force to anteversion-retroversion of the head/stem was 31.5 N/° for gait, which was similar in magnitude to anterior-posterior position of the cup (34.6 N/m for gait). Additionally, superior-inferior cup alignment resulted in 16.4 (4.9)° of variation in the direction of the hip joint contact force across the three tasks, with the most inferior cup placements moving the force vector towards the cup equator at the point of peak joint contact force.

INTERPRETATION

A quantitative understanding of the impact and potential tradeoffs when altering component alignment is valuable in supporting surgical decision making.

Kinematics of mobile-bearing unicompartmental knee arthroplasty compared to native: results from an in vitro study

INTRODUCTION

Fixed-bearing unicompartmental knee arthroplasty (UKA) closely replicates native knee kinematics. As few studies have assessed kinematics following mobile-bearing (MB) UKA, the current study aimed to investigate whether MB UKA preserves natural knee kinematics.

MATERIALS AND METHODS

Seven fresh-frozen full-leg cadaver specimens were prepared and mounted in a kinematic rig that allowed all degrees of freedom at the knee. Three motion patterns, passive flexion-extension (0°-110° flexion), open-chain extension (5°-70° flexion) and squatting (30°-100° flexion), were performed pre- and post-implantation of a medial MB UKA and compared in terms of rotational and translational knee joint kinematics in the different anatomical planes, respectively.

RESULTS

In terms of frontal plane rotational kinematics, MB UKA specimens were in a more valgus orientation for all motion patterns. In the axial plane, internal rotation of the tibia before and after UKA was consistent, regardless of motion task, with no significant differences. In terms of frontal plane, i.e., inferior-superior, translations, the FMCC was significantly higher in UKA knees in all flexion angles and motor tasks, except in early flexion during passive motion. In terms of axial plane, i.e., anteroposterior (AP), translations, during open-chain activities, the femoral medial condyle center (FMCC) tended to be more posterior following UKA relative to the native knee in mid-flexion and above. AP excursions of the FMCC were small in all tested motions, however. There was substantial AP translation of the femoral lateral condyle center during passive motion before and after UKA, which was significantly different for flexion angles > 38°.

CONCLUSIONS

Our study data demonstrate that the kinematics of the unloaded knee following MB UKA closely resemble those of the native knee while relative medial overstuffing with UKA will result in the joint being more valgus. However, replacing the conforming and rigidly fixed medial meniscus with a mobile inlay may successfully prevent aberrant posterior translation of the medial femoral compartment during passive motion and squatting motion.

How Can We Use Computational Modeling to Improve Total Knee Arthroplasty? Modeling Stability and Mobility in the Implanted Knee

Validated computational models promise a virtual platform to create optimal articular surfaces that best achieve desired implant characteristics. Today, designers can parametrically define the primary geometric features of an implant, and automatically modify design variables until stability/mobility performance objectives are best achieved. This preclinical, virtual design iteration minimizes the development cycle compared with testing physical prototypes and, by evaluating a broader scope of design concepts, likely improves the clinical performance of the final product. However, the scenario described is not without shortcomings and requires thorough understanding of the capabilities and the limitations of the models used. Although models typically represent the articular interface well, the interaction with the patient and the surgical process includes significant variability and increase in complexity. We present current modeling capabilities for the estimation of implant stability/mobility, with further suggestions for answering the difficult question of how an implant might perform throughout the population.