Comparison of deformable and elastic foundation finite element simulations for predicting knee replacement mechanics

Rigid body total knee replacement (TKR) models with tibiofemoral contact based on elastic foundation (EF) theory utilize simple contact pressure-surface overclosure relationships to estimate joint mechanics, and require significantly less computational time than corresponding deformable finite element (FE) methods. However, potential differences in predicted kinematics between these representations are currently not well understood, and it is unclear if the estimates of contact area and pressure are acceptable. Therefore, the objectives of the current study were to develop rigid EF and deformable FE models of tibiofemoral contact, and to compare predicted kinematics and contact mechanics from both representations during gait loading conditions with three different implant designs. Linear and nonlinear contact pressure-surface overclosure relationships based on polyethylene material properties were developed using EF theory. All other variables being equal, rigid body FE models accurately estimated kinematics predicted by fully deformable FE models and required only 2% of the analysis time. As expected, the linear EF contact model sufficiently approximated trends for peak contact pressures, but overestimated the deformable results by up to 30%. The nonlinear EF contact model more accurately reproduced trends and magnitudes of the deformable analysis, with maximum differences of approximately 15% at the peak pressures during the gait cycle. All contact area predictions agreed in trend and magnitude. Using rigid models, edge-loading conditions resulted in substantial overestimation of peak pressure. Optimal nonlinear EF contact relationships were developed for specific TKR designs for use in parametric or repetitive analyses where computational time is paramount. The explicit FE analysis method utilized here provides a unique approach in that both rigid and deformable analyses can be run from the same input file, thus enabling simple selection of the most appropriate representation for the analysis of interest.     

Predicting knee replacement damage in a simulator machine using a computational model with a consistent wear factor

Wear of ultrahigh molecular weight polyethylene remains a primary factor limiting the longevity of total knee replacements (TKRs). However, wear testing on a simulator machine is time consuming and expensive, making it impractical for iterative design purposes. The objectives of this paper were first, to evaluate whether a computational model using a wear factor consistent with the TKR material pair can predict accurate TKR damage measured in a simulator machine, and second, to investigate how choice of surface evolution method (fixed or variable step) and material model (linear or nonlinear) affect the prediction. An iterative computational damage model was constructed for a commercial knee implant in an AMTI simulator machine. The damage model combined a dynamic contact model with a surface evolution model to predict how wear plus creep progressively alter tibial insert geometry over multiple simulations. The computational framework was validated by predicting wear in a cylinder-on-plate system for which an analytical solution was derived. The implant damage model was evaluated for 5 million cycles of simulated gait using damage measurements made on the same implant in an AMTI machine. Using a pin-on-plate wear factor for the same material pair as the implant, the model predicted tibial insert wear volume to within 2% error and damage depths and areas to within 18% and 10% error, respectively. Choice of material model had little influence, while inclusion of surface evolution affected damage depth and area but not wear volume predictions. Surface evolution method was important only during the initial cycles, where variable step was needed to capture rapid geometry changes due to the creep. Overall, our results indicate that accurate TKR damage predictions can be made with a computational model using a constant wear factor obtained from pin-on-plate tests for the same material pair, and furthermore, that surface evolution method matters only during the initial "break in" period of the simulation.     

Computational wear prediction of a total knee replacement from in vivo kinematics

Wear of ultra-high molecular weight polyethylene bearings in total knee replacements remains a major limitation to the longevity of these clinically successful devices. Few design tools are currently available to predict mild wear in implants based on varying kinematics, loads, and material properties. This paper reports the implementation of a computer modeling approach that uses fluoroscopically measured motions as inputs and predicts patient-specific implant damage using computationally efficient dynamic contact and tribological analyses. Multibody dynamic simulations of two activities (gait and stair) with two loading conditions (70-30 and 50-50 medial-lateral load splits) were generated from fluoroscopic data to predict contact pressure and slip velocity time histories for individual elements on the tibial insert surface. These time histories were used in a computational wear analysis to predict the depth of damage due to wear and creep experienced by each element. Predicted damage areas, volumes, and maximum depths were evaluated against a tibial insert retrieved from the same patient who provided the in vivo motions. Overall, the predicted damage was in close agreement with damage observed on the retrieval. The gait and stair simulations separately predicted the correct location of maximum damage on the lateral side, whereas a combination of gait and stair was required to predict the correct location on the medial side. Predicted maximum damage depths were consistent with the retrieval as well. Total computation time for each damage prediction was less than 30 min. Continuing refinement of this approach will provide a robust tool for accurately predicting clinically relevant wear in total knee replacements.     

Explicit finite element modeling of total knee replacement mechanics

Joint kinematics and contact mechanics dictate the success of current total knee replacement (TKR) devices. Efficient computer models present an effective way of evaluating these characteristics. Predicted contact stress and area due to articulations at the tibio-femoral and patello-femoral interfaces indicate potential clinical performance. Previous finite element (FE) knee models have generally been used to predict contact stresses and/or areas during static or quasi-static loading conditions. Explicit dynamic FE analyses have recently been used to efficiently predict TKR kinematics and contact mechanics during dynamic loading conditions. The objective of this study was to develop and experimentally validate an explicit FE TKR model that incorporates tibio-femoral and patello-femoral articulations. For computational efficiency, we developed rigid body analyses that can reasonably reproduce the kinematics, contact pressure distribution, and contact area of a fully deformable system. Results from the deformable model showed that the patello-femoral and tibio-femoral kinematics were in good agreement with experimental knee simulator measurements. Kinematic results from the rigid body analyses were nearly identical to those from the fully deformable model, and the contact pressure and contact area correlation was acceptable given the great reduction in analysis time. Component mesh density studied had little effect on the predicted kinematics, particularly for the patellar component, and small effects on the predicted contact pressure and area. These analyses have shown that, at low computational cost, a force-control dynamic simulation of a gait cycle can yield useful and predictable results.     

Three-dimensional dynamic hip contact area and pressure distribution during activities of daily living

Estimation of the hip joint contact area and pressure distribution during activities of daily living is important in predicting joint degeneration mechanism, prosthetic implant wear, providing biomechanical rationales for preoperative planning and postoperative rehabilitation. These biomechanical data were estimated utilizing a generic hip model, the Discrete Element Analysis technique, and the in vivo hip joint contact force data. The three-dimensional joint potential contact area was obtained from the anteroposterior radiograph of a subject and the actual joint contact area and pressure distribution in eight activities of daily living were calculated. During fast, normal, and slow walking, the peak pressure of moderate magnitude was located at the lateral roof of the acetabulum during mid-stance. In standing up and sitting down, and during knee bending, the peak pressures were located at the edge of the posterior horn and the magnitude of the peak pressure during sitting down was 2.8 times that of normal walking. The peak pressure was found at the lateral roof in climbing up stairs which was higher than that in going down stairs. These results can be used to rationalize rehabilitation protocols, functional restrictions after complex acetabular reconstructions, and prosthetic component wear and fatigue test set up. The same model and analysis can provide further insight to soft tissue loading and pathology such as labral injury. When the pressure distribution on the acetabulum is inverted onto the femoral head, prediction of subchondral bone collapse associated with avascular necrosis can be achieved with improved accuracy.     

Mechanical buckling of artery under pulsatile pressure

Tortuosity that often occurs in carotid and other arteries has been shown to be associated with high blood pressure, atherosclerosis, and other diseases. However the mechanisms of tortuosity development are not clear. Our previous studies have suggested that arteries buckling could be a possible mechanism for the initiation of tortuous shape but artery buckling under pulsatile flow condition has not been fully studied. The objectives of this study were to determine the artery critical buckling pressure under pulsatile pressure both experimentally and theoretically, and to elucidate the relationship of critical pressures under pulsatile flow, steady flow, and static pressure. We first tested the buckling pressures of porcine carotid arteries under these loading conditions, and then proposed a nonlinear elastic artery model to examine the buckling pressures under pulsatile pressure conditions. Experimental results showed that under pulsatile pressure arteries buckled when the peak pressures were approximately equal to the critical buckling pressures under static pressure. This was also confirmed by model simulations at low pulse frequencies. Our results provide an effective tool to predict artery buckling pressure under pulsatile pressure.     

A multibody knee model with discrete cartilage prediction of tibio-femoral contact mechanics

Combining musculoskeletal simulations with anatomical joint models capable of predicting cartilage contact mechanics would provide a valuable tool for studying the relationships between muscle force and cartilage loading. As a step towards producing multibody musculoskeletal models that include representation of cartilage tissue mechanics, this research developed a subject-specific multibody knee model that represented the tibia plateau cartilage as discrete rigid bodies that interacted with the femur through deformable contacts. Parameters for the compliant contact law were derived using three methods: (1) simplified Hertzian contact theory, (2) simplified elastic foundation contact theory and (3) parameter optimisation from a finite element (FE) solution. The contact parameters and contact friction were evaluated during a simulated walk in a virtual dynamic knee simulator, and the resulting kinematics were compared with measured in vitro kinematics. The effects on predicted contact pressures and cartilage–bone interface shear forces during the simulated walk were also evaluated. The compliant contact stiffness parameters had a statistically significant effect on predicted contact pressures as well as all tibio-femoral motions except flexion–extension. The contact friction was not statistically significant to contact pressures, but was statistically significant to medial–lateral translation and all rotations except flexion–extension. The magnitude of kinematic differences between model formulations was relatively small, but contact pressure predictions were sensitive to model formulation. The developed multibody knee model was computationally efficient and had a computation time 283 times faster than a FE simulation using the same geometries and boundary conditions.     

Computational wear simulation of patellofemoral articular cartilage during in vitro testing

Though changes in normal joint motions and loads (e.g., following anterior cruciate ligament injury) contribute to the development of knee osteoarthritis, the precise mechanism by which these changes induce osteoarthritis remains unknown. As a first step toward identifying this mechanism, this study evaluates computational wear simulations of a patellofemoral joint specimen wear tested on a knee simulator machine. A multibody dynamic model of the specimen mounted in the simulator machine was constructed in commercial computer-aided engineering software. A custom elastic foundation contact model was used to calculate contact pressures and wear on the femoral and patellar articular surfaces using geometry created from laser scan and MR data. Two different wear simulation approaches were investigated--one that wore the surface geometries gradually over a sequence of 10 one-cycle dynamic simulations (termed the "progressive" approach), and one that wore the surface geometries abruptly using results from a single one-cycle dynamic simulation (termed the "non-progressive" approach). The progressive approach with laser scan geometry reproduced the experimentally measured wear depths and areas for both the femur and patella. The less costly non-progressive approach predicted deeper wear depths, especially on the patella, but had little influence on predicted wear areas. Use of MR data for creating the articular and subchondral bone geometry altered wear depth and area predictions by at most 13%. These results suggest that MR-derived geometry may be sufficient for simulating articular cartilage wear in vivo and that a progressive simulation approach may be needed for the patella and tibia since both remain in continuous contact with the femur.     

Heterogeneous mechanics of the mouse pulmonary arterial network

Individualized modeling and simulation of blood flow mechanics find applications in both animal research and patient care. Individual animal or patient models for blood vessel mechanics are based on combining measured vascular geometry with a fluid structure model coupling formulations describing dynamics of the fluid and mechanics of the wall. For example, one-dimensional fluid flow modeling requires a constitutive law relating vessel cross-sectional deformation to pressure in the lumen. To investigate means of identifying appropriate constitutive relationships, an automated segmentation algorithm was applied to micro-computerized tomography images from a mouse lung obtained at four different static pressures to identify the static pressure–radius relationship for four generations of vessels in the pulmonary arterial network. A shape-fitting function was parameterized for each vessel in the network to characterize the nonlinear and heterogeneous nature of vessel distensibility in the pulmonary arteries. These data on morphometric and mechanical properties were used to simulate pressure and flow velocity propagation in the network using one-dimensional representations of fluid and vessel wall mechanics. Moreover, wave intensity analysis was used to study effects of wall mechanics on generation and propagation of pressure wave reflections. Simulations were conducted to investigate the role of linear versus nonlinear formulations of wall elasticity and homogeneous versus heterogeneous treatments of vessel wall properties. Accounting for heterogeneity, by parameterizing the pressure/distention equation of state individually for each vessel segment, was found to have little effect on the predicted pressure profiles and wave propagation compared to a homogeneous parameterization based on average behavior. However, substantially different results were obtained using a linear elastic thin-shell model than were obtained using a nonlinear model that has a more physiologically realistic pressure versus radius relationship.     

Finite element modeling for predicting the contact pressure between a foam mattress and the human body in a supine position

In this paper, we generated finite element (FE) models to predict the contact pressure between a foam mattress and the human body in a supine position. Twenty-year-old males were used for three-dimensional scanning to produce the FE human models, which was composed of skin and muscle tissue. A linear elastic isotropic material model was used for the skin, and the Mooney–Rivlin model was used for the muscle tissue because it can effectively represent the nonlinear behavior of muscle. The contact pressure between the human model and the mattress was predicted by numerical simulation. The human models were validated by comparing the body pressure distribution obtained from the same human subject when he was lying on two different mattress types. The experimental results showed that the slope of the lower part of the mattress caused a decrease in the contact pressure at the heels, and the effect of bone structure was most pronounced in the scapula. After inserting a simple structure to function as the scapula, the contact pressure predicted by the FE human models was consistent with the experimental body pressure distribution for all body parts. These results suggest that the models proposed in this paper will be useful to researchers and designers of products related to the prevention of pressure ulcers.     

A new technique to measure the dynamic contact pressures on the Tibial Plateau

Studies of the load transfer role of the meniscus have been limited to static experimental and analytical approaches. The objective of this study was to develop an experimental technique to allow the contact pressures on the tibial plateau of cadaveric knees to be measured under dynamic physiological loads. Accordingly, we adapted a load-controlled knee joint simulator to accept a cadaveric sheep knee, programmed the simulator with sheep gait kinematics data, and utilized a pressure sensor array to measure the contact pressure distribution on the lateral tibial plateau during gait. The technique was applied to six sheep knees that were tested intact and after meniscectomy. Meniscectomy resulted in a 267% increase in average contact pressure, a 117% increase in peak contact pressure, and an 80% decrease in contact area, all measured at the point of maximum peak contact stress in the gait cycle. It is envisaged that the experimental model herein developed will allow for the screening of candidate materials prior to more expensive and time-consuming animal models.     

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.

     

Development of a statistical shape-function model of the implanted knee for real-time prediction of joint mechanics

Outcomes of total knee arthroplasty (TKA) are dependent on surgical technique, patient variability, and implant design. Non-optimal design or alignment choices may result in undesirable contact mechanics and joint kinematics, including poor joint alignment, instability, and reduced range of motion. Implant design and surgical alignment are modifiable factors with potential to improve patient outcomes, and there is a need for robust implant designs that can accommodate patient variability. Our objective was to develop a statistical shape-function model (SFM) of a posterior stabilized implanted knee to instantaneously predict joint mechanics in an efficient manner. Finite element methods were combined with Latin hypercube sampling and regression analyses to produce modeling equations relating nine implant design and six surgical alignment parameters to tibiofemoral (TF) joint mechanics outcomes during a deep knee bend. A SFM was developed and TF contact mechanics, kinematics, and soft tissue loads were instantaneously predicted from the model. Average normalized root-mean-square error predictions were between 2.79% and 9.42%, depending on the number of parameters included in the model. The statistical shape-function model generated instantaneous joint mechanics predictions using a maximum of 130 training simulations, making it ideally suited for integration into a patient-specific design and alignment optimization pipeline. Such a tool may be used to optimize kinematic function to achieve more natural motion or minimize implant wear, and may aid the engineering and clinical communities in improving patient satisfaction and surgical outcomes.     

Correlation of compartment pressure data from an intraoperative sensing device with postoperative fluoroscopic kinematic results in TKA patients

Fluoroscopy has recently been used to analyze postoperative kinematics in total knee arthroplasty (TKA). These analyses have reported varying results even in patients with similar implant design. In addition, patterns of wear in retrieved tibial polyethylene inserts of similar design have been found to vary substantially. These findings suggest that surgical technique, especially soft tissue balancing, may play a role in postoperative kinematics and implant failure. Accurate soft-tissue balancing is hypothesized to result in similar pressures within the medial and lateral compartments of the knee. However, a method of easily measuring these pressures at TKA has not been developed. In the present study, 32 patients were implanted with a mobile-bearing LCS TKA utilizing the balanced gap technique. An electronic pressure sensor, developed specifically to record pressure magnitude and distribution in the medial and lateral compartments, was incorporated into the implant trials. The knee was then passively taken through a range of motion while pressure data was recorded via computer. Postoperatively, 16 patients underwent active fluoroscopic kinematic analysis to assess for condylar liftoff and femorotibial translation. We found that abnormal compartment pressures and distributions as recorded by the intraoperative pressure sensor were correlated with inappropriate or paradoxical postoperative kinematics. In addition, subjects having similar pressures in both compartments throughout a range of motion did not experience condylar liftoff values greater than 1.0 mm. These data suggest that surgical technique influences the magnitude and distribution of forces at the articulation, postoperative kinematics, and likely, implant longevity.     

The effect of valgus/varus malalignment on load distribution in total knee replacements

Valgus or varus malpositioning of the tibial component of a total knee implant may cause increased propensity for loosening or implant wear and eventually may lead to revision surgery. The aim of this study was to determine the effect of valgus/varus malalignment on tibio-femoral mechanics during surgical trial reduction and simulated gait loading. In seven cadaver legs, posterior cruciate sparing total knee replacements were implanted and tibial inserts representing a neutral alignment and 3 degrees and 5 degrees varus and valgus alignments were sequentially inserted. Each knee with each insert was loaded in a manner representative of a trial reduction performed during knee surgery and loaded in a physiological knee simulator. Simulated gait performed on the simulator demonstrated that internal/external and adduction/abduction rotations showed statistical changes with some of the angled inserts at different points in the walking cycle. Neither medial/lateral nor anterior/posterior translations changed statistically during simulated walking. The pressure distribution and total load in the medial and lateral compartments of the tibial component changed significantly with as little as a 3 degrees variation in angulation when loaded in a manner representative of a trial reduction or with a knee simulator. These results support the need for precise surgical reconstruction of the mechanical axis of the knee and proper alignment of the tibial component. These results further demonstrate that tibial contact pressures measured during a trial reduction method may be predictive of contact mechanics at the higher loading seen in the knee simulator.     

An inverse finite-element model of heel-pad indentation

A numerical-experimental approach has been developed to characterize heel-pad deformation at the material level. Left and right heels of 20 diabetic subjects and 20 nondiabetic subjects matched for age, gender and body mass index were indented using force-controlled ultrasound. Initial tissue thickness and deformation were measured using M-mode ultrasound; indentation forces were recorded simultaneously. An inverse finite-element analysis of the indentation protocol using axisymmetric models adjusted to reflect individual heel thickness was used to extract nonlinear material properties describing the hyperelastic behavior of each heel. Student's t-tests revealed that heel pads of diabetic subjects were not significantly different in initial thickness nor were they stiffer than those from nondiabetic subjects. Another heel-pad model with anatomically realistic surface representations of the calcaneus and soft tissue was developed to estimate peak pressure prediction errors when average rather than individualized material properties were used. Root-mean-square errors of up to 7% were calculated, indicating the importance of subject-specific modeling of the nonlinear elastic behavior of the heel pad. Indentation systems combined with the presented numerical approach can provide this information for further analysis of patient-specific foot pathologies and therapeutic footwear designs.     

In vivo contact kinematics and contact forces of the knee after total knee arthroplasty during dynamic weight-bearing activities

Analysis of polyethylene component wear and implant loosening in total knee arthroplasty (TKA) requires precise knowledge of in vivo articular motion and loading conditions. This study presents a simultaneous in vivo measurement of tibiofemoral articular contact forces and contact kinematics in three TKA patients. These measurements were accomplished via a dual fluoroscopic imaging system and instrumented tibial implants, during dynamic single leg lunge and chair rising-sitting. The measured forces and contact locations were also used to determine mediolateral distribution of axial contact forces. Contact kinematics data showed a medial pivot during flexion of the knee, for all patients in the study. Average axial forces were higher for lunge compared to chair rising-sitting (224% vs. 187% body weight). In this study, we measured peak anteroposterior and mediolateral forces averaging 13.3% BW during lunge and 18.5% BW during chair rising-sitting. Mediolateral distributions of axial contact force were both patient and activity specific. All patients showed equitable medial-lateral loading during lunge but greater loads at the lateral compartment during chair rising-sitting. The results of this study may enable more accurate reproduction of in vivo loads and articular motion patterns in wear simulators and finite element models. This in turn may help advance our understanding of factors limiting longevity of TKA implants, such as aseptic loosening and polyethylene component wear, and enable improved TKA designs.     

Estimation of discretization errors in contact pressure measurements

Contact pressure measurements in total knee replacements are often made using a discrete sensor such as the Tekscan K-Scan sensor. However, no method currently exists for predicting the magnitude of sensor discretization errors in contact force, peak pressure, average pressure, and contact area, making it difficult to evaluate the accuracy of such measurements. This study identifies a non-dimensional area variable, defined as the ratio of the number of perimeter elements to the total number of elements with pressure, which can be used to predict these errors. The variable was evaluated by simulating discrete pressure sensors subjected to Hertzian and uniform pressure distributions with two different calibration procedures. The simulations systematically varied the size of the sensor elements, the contact ellipse aspect ratio, and the ellipse's location on the sensor grid. In addition, contact pressure measurements made with a K-Scan sensor on four different total knee designs were used to evaluate the magnitude of discretization errors under practical conditions. The simulations predicted a strong power law relationship (r(2)>0.89) between worst-case discretization errors and the proposed non-dimensional area variable. In the total knee experiments, predicted discretization errors were on the order of 1-4% for contact force and peak pressure and 3-9% for average pressure and contact area. These errors are comparable to those arising from inserting a sensor into the joint space or truncating pressures with pressure sensitive film. The reported power law regression coefficients provide a simple way to estimate the accuracy of experimental measurements made with discrete pressure sensors when the contact patch is approximately elliptical.     

Development and validation of a computational model of the knee joint for the evaluation of surgical treatments for osteoarthritis

A three-dimensional (3D) knee joint computational model was developed and validated to predict knee joint contact forces and pressures for different degrees of malalignment. A 3D computational knee model was created from high-resolution radiological images to emulate passive sagittal rotation (full-extension to 65°-flexion) and weight acceptance. A cadaveric knee mounted on a six-degree-of-freedom robot was subjected to matching boundary and loading conditions. A ligament-tuning process minimised kinematic differences between the robotically loaded cadaver specimen and the finite element (FE) model. The model was validated by measured intra-articular force and pressure measurements. Percent full scale error between FE-predicted and in vitro-measured values in the medial and lateral compartments were 6.67% and 5.94%, respectively, for normalised peak pressure values, and 7.56% and 4.48%, respectively, for normalised force values. The knee model can accurately predict normalised intra-articular pressure and forces for different loading conditions and could be further developed for subject-specific surgical planning.