Lorin P. Maletsky

Dynamic finite element knee simulation for evaluation of knee replacement mechanics

In vitro pre-clinical testing of total knee replacement (TKR) devices is a necessary step in the evaluation of new implant designs. Whole joint knee simulators, like the Kansas knee simulator (KKS), provide a controlled and repeatable loading environment for comparative evaluation of component designs or surgical alignment under dynamic conditions. Experimental testing, however, is time and cost prohibitive for design-phase evaluation of tens or hundreds of design variations. Experimentally-verified computational models provide an efficient platform for analysis of multiple components, sizes, and alignment conditions. The purpose of the current study was to develop and verify a computational model of a dynamic, whole joint knee simulator. Experimental internal-external and valgus-varus laxity tests, followed by dynamic deep knee bend and gait simulations in the KKS were performed on three cadaveric specimens. Specimen-specific finite element (FE) models of posterior-stabilized TKR were created from magnetic resonance images and CAD geometry. The laxity data was used to optimize mechanical properties of tibiofemoral soft-tissue structures on a specimen-specific basis. Each specimen was subsequently analyzed in a computational model of the experimental KKS, simulating both dynamic activities. The computational model represented all joints and actuators in the experimental setup, including a proportional-integral-derivative (PID) controller to drive quadriceps actuation. The computational model was verified against six degree-of-freedom patellofemoral (PF) and tibiofemoral (TF) kinematics and actuator loading during both deep knee bend and gait activities, with good agreement in trends and magnitudes between model predictions and experimental kinematics; differences were less than 1.8 mm and 2.2° for PF and TF translations and rotations. The whole joint FE simulator described in this study can be applied to investigate a wide range of clinical and research questions.

Simulating Dynamic Activities Using a Five-Axis Knee Simulator

This work describes the design and capabilities of the Purdue Knee Simulator: Mark II and a sagittal-plane model of the machine. This five-axis simulator was designed and constructed to simulate dynamic loading activities on either cadaveric knee specimens or total knee prostheses mounted on fixtures. The purpose of the machine was to provide a consistent, realistic loading of the knee joint, allowing the kinematics and specific loading of the structures of the knee to be determined based on condition, articular geometry, and simulated activity. The sagittal-plane model of the knee simulator was developed both to predict the loading at the knee from arbitrary inputs and to generate the necessary inputs required to duplicate specified joint loading. Measured tibio-femoral compressive force and quadriceps tension were shown to be in good agreement with the predicted loads from the model. A controlled moment about the ankle-flexion axis was also shown to change the loading on the quadriceps.

Verification of predicted specimen-specific natural and implanted patellofemoral kinematics during simulated deep knee bend.

Verified computational models represent an efficient method for studying the relationship between articular geometry, soft-tissue constraint, and patellofemoral (PF) mechanics. The current study was performed to evaluate an explicit finite element (FE) modeling approach for predicting PF kinematics in the natural and implanted knee. Experimental three-dimensional kinematic data were collected on four healthy cadaver specimens in their natural state and after total knee replacement in the Kansas knee simulator during a simulated deep knee bend activity. Specimen-specific FE models were created from medical images and CAD implant geometry, and included soft-tissue structures representing medial-lateral PF ligaments and the quadriceps tendon. Measured quadriceps loads and prescribed tibiofemoral kinematics were used to predict dynamic kinematics of an isolated PF joint between 10 degrees and 110 degrees femoral flexion. Model sensitivity analyses were performed to determine the effect of rigid or deformable patellar representations and perturbed PF ligament mechanical properties (pre-tension and stiffness) on model predictions and computational efficiency. Predicted PF kinematics from the deformable analyses showed average root mean square (RMS) differences for the natural and implanted states of less than 3.1 degrees and 1.7 mm for all rotations and translations. Kinematic predictions with rigid bodies increased average RMS values slightly to 3.7 degrees and 1.9 mm with a five-fold decrease in computational time. Two-fold increases and decreases in PF ligament initial strain and linear stiffness were found to most adversely affect kinematic predictions for flexion, internal-external tilt and inferior-superior translation in both natural and implanted states. The verified models could be used to further investigate the effects of component alignment or soft-tissue variability on natural and implant PF mechanics.

Verification of Predicted Knee Replacement Kinematics During Simulated Gait in the Kansas Knee Simulator

Evaluating total knee replacement kinematics and contact pressure distributions is an important element of preclinical assessment of implant designs. Although physical testing is essential in the evaluation process, validated computational models can augment these experiments and efficiently evaluate perturbations of the design or surgical variables. The objective of the present study was to perform an initial kinematic verification of a dynamic finite element model of the Kansas knee simulator by comparing predicted tibio- and patellofemoral kinematics with experimental measurements during force-controlled gait simulation. A current semiconstrained, cruciate-retaining, fixed-bearing implant mounted in aluminum fixtures was utilized. An explicit finite element model of the simulator was developed from measured physical properties of the machine, and loading conditions were created from the measured experimental feedback data. The explicit finite element model allows both rigid body and fully deformable solutions to be chosen based on the application of interest. Six degrees-of-freedom kinematics were compared for both tibio- and patellofemoral joints during gait loading, with an average root mean square (rms) translational error of 1.1 mm and rotational rms error of 1.3 deg. Model sensitivity to interface friction and damping present in the experimental joints was also evaluated and served as a secondary goal of this paper. Modifying the metal-polyethylene coefficient of friction from 0.1 to 0.01 varied the patellar flexion-extension and tibiofemoral anterior-posterior predictions by 7 deg and 2 mm, respectively, while other kinematic outputs were largely insensitive.

The influence of total knee arthroplasty geometry on mid-flexion stability: An experimental and finite element study

Fluoroscopic evaluation of total knee arthroplasty (TKA) has reported sudden anterior translation of the femur relative to the tibia (paradoxical anterior motion) for some cruciate-retaining designs. This motion may be tied to abrupt changes in the femoral sagittal radius of curvature characteristic of traditional TKA designs, as the geometry transitions from a large load-bearing distal radius to a smaller posterior radius which can accommodate femoral rollback. It was hypothesized that a gradually reducing radius may attenuate sudden changes in anterior-posterior motion that occur in mid-flexion with traditional discrete-radius designs. A combined experimental and computational approach was employed to test this hypothesis. A previously developed finite element (FE) model of the Kansas knee simulator (KKS), virtually implanted with multiple implant designs, was used to predict the amount of paradoxical anterior femoral slide during a simulated deep knee bend. The model predicted kinematics demonstrated that incorporating a gradually reducing radius in mid-flexion reduced the magnitude of paradoxical anterior translation between 21% and 68%, depending on the conformity of the tibial insert. Subsequently, both a dual-radius design and a modified design incorporating gradually reducing radii were tested in vitro in the KKS for verification. The model-predicted and experimentally observed kinematics exhibited good agreement, while the average experimental kinematics demonstrated an 81% reduction in anterior translation with the modified design. The FE model demonstrated sufficient sensitivity to appropriately differentiate kinematic changes due to subtle changes in implant design, and served as a useful pre-clinical design-phase tool to improve implant kinematics.

Computational Knee Ligament Modeling Using Experimentally Determined Zero-Load Lengths

This study presents a subject-specific method of determining the zero-load lengths of the cruciate and collateral ligaments in computational knee modeling. Three cadaver knees were tested in a dynamic knee simulator. The cadaver knees also underwent manual envelope of motion testing to find their passive range of motion in order to determine the zero-load lengths for each ligament bundle. Computational multibody knee models were created for each knee and model kinematics were compared to experimental kinematics for a simulated walk cycle. One-dimensional non-linear spring damper elements were used to represent cruciate and collateral ligament bundles in the knee models. This study found that knee kinematics were highly sensitive to altering of the zero-load length. The results also suggest optimal methods for defining each of the ligament bundle zero-load lengths, regardless of the subject. These results verify the importance of the zero-load length when modeling the knee joint and verify that manual envelope of motion measurements can be used to determine the passive range of motion of the knee joint. It is also believed that the method described here for determining zero-load length can be used for in vitro or in vivo subject-specific computational models.

Evaluating knee replacement mechanics during ADL with PID-controlled dynamic finite element analysis

Validated computational knee simulations are valuable tools for design phase development of knee replacement devices. Recently, a dynamic finite element (FE) model of the Kansas knee simulator was kinematically validated during gait and deep flexion cycles. In order to operate the computational simulator in the same manner as the experiment, a proportional–integral–derivative (PID) controller was interfaced with the FE model to control the quadriceps actuator excursion and produce a target flexion profile regardless of implant geometry or alignment conditions. The controller was also expanded to operate multiple actuators simultaneously in order to produce in vivo loading conditions at the joint during dynamic activities. Subsequently, the fidelity of the computational model was improved through additional muscle representation and inclusion of relative hip–ankle anterior–posterior (A–P) motion. The PID-controlled model was able to successfully recreate in vivo loading conditions (flexion angle, compressive joint load, medial–lateral load distribution or varus-valgus torque, internal–external torque, A–P force) for deep knee bend, chair rise, stance-phase gait and step-down activities.

Computational modeling of a dynamic knee simulator for reproduction of knee loading.

As a first step towards reproducing desired three-dimensional joint loading and motion on a dynamic knee simulator, the goal of this study was to develop and verify a three-dimensional computational model that generated control profiles for the simulator using desired knee loading and motion as model inputs. The developed model was verified by predicting tibio-femoral loading on an instrumented analog knee for given actuator forces and the ability to generate simulator control profiles was demonstrated using a three-dimensional walking profile. The model predicted axial tibia loading for a sagittal-plane dual-limb squat within 1% of measured peak loading. Adding out-of-sagittal-plane forces decreased the accuracy of load prediction. The model generated control profiles to the simulator that produced axial tibia loading within 16% of desired for walking. Discrepancies in predicted and measured quadriceps forces influenced the accuracy of the generated control profiles. Future work will replace the analog knee in both the model and machine with a prosthetic knee.

Mechanics of post‐cam engagement during simulated dynamic activity

Posterior‐stabilized (PS) total knee arthroplasty (TKA) components employ a tibial post and femoral cam mechanism to guide anteroposterior knee motion in lieu of the posterior cruciate ligament. Some PS TKA patients report a clicking sensation when the post and cam engage, while severe wear and fracture of the post; we hypothesize that these complications are associated with excessive impact velocity at engagement. We evaluated the effect of implant design on engagement dynamics of the post‐cam mechanism and resulting polyethylene stresses during dynamic activity. In vitro simulation of a knee bend activity was performed for four cadaveric specimens implanted with PS TKA components. Post‐cam engagement velocity and flexion angle at initial contact were determined. The experimental data were used to validate computational predictions of PS mechanics using the same loading conditions. A lower limb model was subsequently utilized to compare engagement mechanics of eight TKA designs, relating differences between implants to geometric design features. Flexion angle and post‐cam velocity at engagement demonstrated considerable ranges among designs (23°–89°, and 0.05–0.22 mm/°, respectively). Post‐cam velocity was correlated (r = 0.89) with tibiofemoral condylar design features. Condylar geometry, in addition to post‐cam geometry, played a significant role in minimizing engagement velocity and forces and stresses in the post. This analysis guides selection and design of PS implants that facilitate smooth post‐cam engagement and reduce edge loading of the post.

Statistical Modeling to Characterize Relationships Between Knee Anatomy and Kinematics

The mechanics of the knee are complex and dependent on the shape of the articular surfaces and their relative alignment. Insight into how anatomy relates to kinematics can establish biomechanical norms, support the diagnosis and treatment of various pathologies (e.g., patellar maltracking) and inform implant design. Prior studies have used correlations to identify anatomical measures related to specific motions. The objective of this study was to describe relationships between knee anatomy and tibiofemoral (TF) and patellofemoral (PF) kinematics using a statistical shape and function modeling approach. A principal component (PC) analysis was performed on a 20‐specimen dataset consisting of shape of the bone and cartilage for the femur, tibia and patella derived from imaging and six‐degree‐of‐freedom TF and PF kinematics from cadaveric testing during a simulated squat. The PC modes characterized links between anatomy and kinematics; the first mode captured scaling and shape changes in the condylar radii and their influence on TF anterior–posterior translation, internal‐external rotation, and the location of the femoral lowest point. Subsequent modes described relations in patella shape and alta/baja alignment impacting PF kinematics. The complex interactions described with the data‐driven statistical approach provide insight into knee mechanics that is useful clinically and in implant design.