Quantification of soft tissue balance in total knee arthroplasty using finite element analysis

Unbalanced contact force on the tibial component has been considered a factor leading to loosening of the implant and increased wear of the bearing surface in total knee arthroplasty. Because it has been reported that good alignment cannot guarantee successful clinical outcomes, the soft tissue balance should be checked together with the alignment. Finite element models of patients' lower extremities were developed to analyse the medial and lateral contact force distribution on the tibial insert. The distributions for four out of five patients were not balanced equally, even though the alignment angles were within a clinically acceptable range. Moreover, the distribution was improved by changing soft tissue release and ligament tightening for the specific case. Integration of the biomechanical modelling, image matching and finite element analysis techniques with the patient-specific properties and various dynamic loading would suggest a clinically relevant pre-operative planning for soft tissue balancing.     

Development and validation of a weight-bearing finite element model for total knee replacement

Total knee arthroplasty (TKA) is a successful procedure for osteoarthritis. However, some patients (19%) do have pain after surgery. A finite element model was developed based on boundary conditions of a knee rig. A 3D-model of an anatomical full leg was generated from magnetic resonance image data and a total knee prosthesis was implanted without patella resurfacing. In the finite element model, a restarting procedure was programmed in order to hold the ground reaction force constant with an adapted quadriceps muscle force during a squat from 20° to 105° of flexion. Knee rig experimental data were used to validate the numerical model in the patellofemoral and femorotibial joint. Furthermore, sensitivity analyses of Young’s modulus of the patella cartilage, posterior cruciate ligament (PCL) stiffness, and patella tendon origin were performed. Pearson’s correlations for retropatellar contact area, pressure, patella flexion, and femorotibial ap-movement were near to 1. Lowest root mean square error for retropatellar pressure, patella flexion, and femorotibial ap-movement were found for the baseline model setup with Young’s modulus of 5 MPa for patella cartilage, a downscaled PCL stiffness of 25% compared to the literature given value and an anatomical origin of the patella tendon. The results of the conducted finite element model are comparable with the experimental results. Therefore, the finite element model developed in this study can be used for further clinical investigations and will help to better understand the clinical aspects after TKA with an unresurfaced patella.     

Comparison between mechanical stress and bone mineral density in the femur after total hip arthroplasty by using subject-specific finite element analyses

The mechanism underling bone mineral density (BMD) loss that occurs in the femur after total hip arthroplasty (THA) remains unknown. We compared the equivalent stress and strain energy density (SED) to BMD in the femur after THA using subject-specific finite element analyses. Twenty-four patients who had undergone primary cementless THA were analysed. BMD was measured using dual-energy X-ray absorptiometry (DEXA) at 1 week and 3, 6 and 12 months after THA. Seven regions of interest (ROIs) were defined in accordance with Gruen's system (ROIs 1–7). Computed tomography images of the femurs were acquired pre- and postoperatively, and the images were converted into three-dimensional finite element (FE) models. Equivalent stress and SED were analysed and compared with DEXA data. BMD was maintained 1 year after THA in ROIs 3, 4, 5 and 6, whereas BMD decreased in ROIs 1, 2 and 7. FE analysis revealed that equivalent stress in ROIs 3, 4, 5 and 6 was much higher than that in ROIs 1, 2 and 7. A significant correlation was observed between the rate of changes in BMD and equivalent stress. Reduction of equivalent stress may contribute to decrease in BMD in the femur after THA.     

Computational analysis of knee joint stability following total knee arthroplasty

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

Finite element study of human lumbar disc nucleus replacements

Currently, there are a number of nucleus replacements under development. The important concern is how well these implants duplicate the mechanical function of the native nucleus. This finite element model study aimed to investigate the influence of different nucleus replacements on the mechanical response of the disc. Models included partial, full, over-sized, partially saturated, elastic and poroelastic solid replacements. Over-sized nucleus replacements up to 25% yielded results that were comparable to those in the intact state. Differences were much greater in cases with under-sized nucleus replacements. The effect was most pronounced for the 75% under-sized replacement that resembled the condition with a full nucleotomy. Nucleus implants with elastic properties substantially altered load transmission when 10% under-sized and over-sized replacements were considered. Compared to intact, the under-sized implants should be avoided when using biphasic materials with properties similar to the native nucleus, whereas for elastic replacements both under- and over-sized implants should not be used.     

Biomechanical analysis of the Universal 2 implant in total wrist arthroplasty: a finite element study

Little is known about the mechanics of in vivo loading on total wrist prostheses where many studies have looked at the mechanics of other types of arthroplasty such as for the hip and the knee which has contributed to the overall success of these types of procedures. Currently surgeons would prefer to carry out arthrodesis on the wrist rather than consider arthroplasty as clinical data have shown that the outcome of total wrist arthroplasty is poorer than compared to the hip and knee. More research is needed on the loading mechanisms of the implants in order to enhance the design of future generation implants. This study looks at the load transfer characteristics of the Universal 2 implant using a finite element model of a virtually implanted prosthesis during gripping. The results showed that the loading on the implant is higher on the dorsal and ulnar aspect than on the volar and radial aspect of the implant. The whole load is transmitted through the radius and none through the ulna.     

The development,calibration and validation of a numerical total knee replacement kinematics simulator considering laxity and unconstrained flexion motions

Kinematics testing is essential during the development of total knee replacement (TKR) designs. Although computational analysis cannot replace physical testing, it offers repeatability and consistency at a much lower cost and shorter time, making it an excellent complement to experiments. Previous numerical models have been limited by several factors: the validity of the models is usually only considered for a single TKR design, friction models are typically overly simplified and the determination of simulation parameters is often inadequate, or tedious and expensive. The objective of this study is to develop, calibrate and validate a TKR kinematics simulation considering multiple TKR geometries, an accurate friction model and simulation parameters determined using a systematic optimisation method. The calibrated model was able to predict TKR kinematics for different TKR geometries, and is ideal for screening new implant designs, reducing the number of experiments required at the design stage.     

Time-dependent elastohydrodynamic lubrication analysis of total knee replacement under walking conditions

This work is concerned with the lubrication analysis of artificial knee joints, which plays an increasing significant role in clinical performance and longevity of components. Time-dependent elastohydrodynamic lubrication analysis for normal total knee replacement is carried out under the cyclic variation in both load and speed representative of normal walking. An equivalent ellipsoid-on-plane model is adopted to represent an actual artificial knee. A full numerical method is developed to simultaneously solve the Reynolds and elasticity equations using the multigrid technique. The elastic deformation is based on the constrained column model. Results show that, under the combined effect of entraining and squeeze-film actions throughout the walking cycle, the predicted central film thickness tends to decrease in the stance phase but keeps a relatively larger value at the swing phase. Furthermore, the geometry of knee joint implant is verified to play an important role under its lubrication condition, and the length of time period is a key point to influence the lubrication performance of joint components.     

The effect of implant design of linked total elbow arthroplasty on stability and stress: a finite element analysis

Several linked total elbow arthroplasty designs exist, which function similar to a loose hinge joint. Constraint behaviour is an important design consideration, as it affects joint stability, or how much secondary [e.g. varus–valgus (VV)] motion is permitted. Implant durability is also a concern, as bearing failures have been reported. This finite element analysis investigates the constraint characteristics and ultra high molecular weight polyethylene bearing stresses of three linked elbow design concepts [cylindrical (CY), hourglass (HG) and concave cylinder (CC)]. The bearing of the CY design was subjected to elevated Von Mises stresses (2.1–5.4 times higher than the HG and CC designs) due to edge loading. The HG design maintained low stresses, but was unable to provide consistent VV stability. The CC design also maintained low stresses while providing consistent VV stability. These results suggest that CC designs may provide better stability characteristics and durability in vivo, compared to the other two designs.     

Developing simulations to reproduce in vivo fluoroscopy kinematics in total knee replacement patients

For clinically predictive testing and design-phase evaluation of prospective total knee replacement (TKR) implants, devices should ideally be evaluated under physiological loading conditions which incorporate population-level variability. A challenge exists for experimental and computational researchers in determining appropriate loading conditions for wear and kinematic knee simulators which reflect in vivo joint loading conditions. There is a great deal of kinematic data available from fluoroscopy studies. The purpose of this work was to develop computational methods to derive anterior–posterior (A–P) and internal–external (I–E) tibiofemoral (TF) joint loading conditions from in vivo kinematic data. Two computational models were developed, a simple TF model, and a more complex lower limb model. These models were driven through external loads applied to the tibia and femur in the TF model, and applied to the hip, ankle and muscles in the lower limb model. A custom feedback controller was integrated with the finite element environment and used to determine the external loads required to reproduce target kinematics at the TF joint. The computational platform was evaluated using in vivo kinematic data from four fluoroscopy patients, and reproduced in vivo A–P and I–E motions and compressive force with a root-mean-square (RMS) accuracy of less than 1 mm, 0.1°, and 40 N in the TF model and in vivo A–P and I–E motions, TF flexion, and compressive loads with a RMS accuracy of less than 1 mm, 0.1°, 1.4°, and 48 N in the lower limb model. The external loading conditions derived from these models can ultimately be used to establish population variability in loading conditions, for eventual use in computational as well as experimental activity simulations.     

A novel formulation for scratch-based wear modelling in total hip arthroplasty

Damage to the femoral head in total hip arthroplasty often takes the form of discrete scratches, which can lead to dramatic wear acceleration of the polyethylene (PE) liner. Here, a novel formulation is reported for finite element (FE) analysis of wear acceleration due to scratch damage. A diffused-light photography technique was used to globally locate areas of damage, providing guidance for usage of high-magnification optical profilometry to determine individual scratch morphology. This multiscale image combination allowed comprehensive input of scratch-based damage patterns to an FE Archard wear model, to determine the wear acceleration associated with specific retrieval femoral heads. The wear algorithm imposed correspondingly elevated wear factors on areas of PE incrementally overpassed by individual scratches. Physical validation was provided by agreement with experimental data for custom-ruled scratch patterns. Illustrative wear acceleration results are presented for four retrieval femoral heads.     

Translation of biomechanical concepts in bone tissue engineering: from animal study to revision knee arthroplasty

Bone defects in revision knee arthroplasty are often located in load-bearing regions. The goal of this study was to determine whether a physiologic load could be used as an in situ osteogenic signal to the scaffolds filling the bone defects. In order to answer this question, we proposed a novel translation procedure having four steps: (1) determining the mechanical stimulus using finite element method, (2) designing an animal study to measure bone formation spatially and temporally using micro-CT imaging in the scaffold subjected to the estimated mechanical stimulus, (3) identifying bone formation parameters for the loaded and non-loaded cases appearing in a recently developed mathematical model for bone formation in the scaffold and (4) estimating the stiffness and the bone formation in the bone-scaffold construct. With this procedure, we estimated that after 3 years mechanical stimulation increases the bone volume fraction and the stiffness of scaffold by 1.5- and 2.7-fold, respectively, compared to a non-loaded situation.     

In-vivo analysis of sliding distance and cross-shear in Bi-cruciate retaining total knee arthroplasty

Polyethylene remains the most popular bearing material for total knee arthroplasty (TKA). Despite its widespread use, wear continue to be one of major factors implicated in revision surgery. Sliding distance, cross-shear, and contact stress are the major factors influencing polyethylene wear. As previous studies have either relied on wear simulations, computational modeling, or in vitro measurements to quantify sliding distance and cross-shear, in vivo subject-specific sliding distance and cross-shear after bi-cruciate retaining (BCR) TKA has not been previously reported. The objective of this study was to quantify the 6°-of-freedom (6DOF) in vivo kinematics, sliding distance, and cross-shear in BCR TKA patients during gait. Twenty-nine unilateral BCR TKA patients performed level walking on a treadmill under dual fluoroscopic imaging system (DFIS) surveillance. Cumulative normalized sliding distances between the lateral and medial compartments did not change significantly (p > 0.05) during the gait cycle. Although the total normalized sliding distance was similar between the lateral and medial compartments, the cross-shear at the lateral compartment differed significantly from that at the medial compartment (p < 0.001). Significant differences in the relative length positions of the peak sliding distance and cross-shear were found between the lateral and medial bearing components. The flexion-extension motion of the reconstructed knee was more associated with the linear displacements (anterior-posterior, R² = 0.6; lateral-medial, R² = 0.8, proximal-distal, R² = 0.7) than the angular displacement (varus-valgus, R² = 0.18; internal-external rotation, R² = 0.28). Despite some differences in peak sliding distance and cross-shear positons, our results suggest similar articular contact patterns between the lateral and medial compartments in BCR TKA patients during gait. The data could provide insights into understanding the potential wear patterns in BCR TKAs.     

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.     

Surgeon variability in total knee arthroplasty component alignment: a Monte Carlo analysis

Component mal-alignment in total knee arthroplasty has been associated with increased revision rates and poor clinical outcomes. A significant source of variability in traditional, jig-based total knee arthroplasty is the performance of the surgeon. The purpose of this study was to determine the most sensitive steps in the femoral and tibia arthroplasty procedures. A computational model of the total knee arthroplasty procedure was created, and Monte Carlo simulations were performed that included surgeon variability in each step of the procedure. The proportion of well-aligned components from the model agrees with clinical literature in most planes. When components must be aligned within ±3° in all planes, component alignment was most sensitive to the accuracy of identifying the lateral epicondyle for the femoral component, and to the precision of the transverse plane alignment of the extramedullary guide for the tibial component. This model can be used as a tool for evaluating different procedural approaches or sources of variability to improve the quality of the total knee arthroplasty procedure.     

Biomechanical behaviour of cancellous bone on patellofemoral arthroplasty with Journey prosthesis: a finite element study

Isolated patellofemoral (PF) arthritis of the knee is a common cause of anterior knee pain and disability. Patellofemoral arthroplasty (PFA) is a bone conserving solution for patients with PF degeneration. Failure mechanisms of PFA include growing tibiofemoral arthritis and loosening of components. The implant loosening can be associated with bone resorption or fatigue-failure of bone by overload. This research work aims at determining the structural effects of the implantation of PF prosthesis Journey PFJ (Smith & Nephew, Inc., Memphis, TN, USA) on femoral cancellous bone. For this purpose, the finite element method is considered to perform computational simulations for different conditions, such as well-fixed and loosening scenarios. From the global results obtained, in the well-fixed scenario, a decrease in strain on cancellous bone was noticed, which can be related to bone resorption. In the loosening scenario, when the cement layer becomes inefficient, a significant increase in cancellous bone strain was observed, which can be associated with bone fatigue-failure.These strain changes suggest a weakness of the femur after PFA.     

Finding and defining the ideal patellar resection plane in total knee arthroplasty

Asymmetric resection of the patella during total knee arthroplasty (TKA) correlates with anterior knee pain, bony impingement and patellar maltracking. Despite this, there is no consensus regarding the desired landmarks; the cut is often done freehand; and there has been no quantitative comparison of proposed resection planes. The objectives of this study were to: determine the intra- and inter-surgeon repeatability of two radiographic resection definitions (medial–divot, MD, and medial–lateral extents, MLE); calculate two additional definitions from the radiographic patellar circumferences (parallel to the anterior surface, ANT, and perpendicular to the anteroposterior tangent points, PERP); compare the clinical resection line to the previous four definitions before and after introducing the MD method clinically; and identify distinguishing features of patellae with better vs. worse resection angles. We hypothesized that the MD method would improve repeatability both radiographically and clinically, that the different radiographic definitions would produce comparable angles, and that we could identify distinguishing features. For the radiographic study, three surgeons drew lines on 40 preoperative X-rays plus 9 interspersed repetitions of 3 of these X-rays. For the clinical study, we compared the patellar resection angle for 20 patients immediately before and after implementing the new method. Given that the clinical goal is to have equal distances from the resection surface to the anterior surface, we compared all results to the ANT definition as the theoretically ideal definition. Confirming the first hypothesis, intra-surgeon repeatability (10 repetitions of 3 X-rays) and inter-surgeon repeatability (3 surgeons×40 X-rays) were both significantly better using the new MD method compared to the MLE method (p<0.001). Contrary to the second hypothesis, clinical use of the MD method did not improve resection symmetry. Contrary to the third hypothesis, the PERP definition was significantly different from the other three definitions. In agreement with the fourth hypothesis, female patellae and more deformed patella had significantly greater asymmetry (p<0.001). Given the inherent variability shown in drawing the ‘patellar horizon’, we encourage researchers to draw the line several times and average the results when comparing tilt or the resection angle to this horizon. Based on the distinguishing characteristics of asymmetrically resurfaced patellae in our series, we recommend that clinicians be particularly careful when resecting laterally deformed patellae and the patellae of female patients.     

Estimating total knee replacement joint load ratios from kinematics

Accurate prediction of loads acting at the joint in total knee replacement (TKR) patients is key to developing experimental or computational simulations which evaluate implant designs under physiological loading conditions. In vivo joint loads have been measured for a small number of telemetric TKR patients, but in order to assess device performance across the entire patient population, a larger patient cohort is necessary. This study investigates the accuracy of predicting joint loads from joint kinematics. Specifically, the objective of the study was to assess the accuracy of internal–external (I–E) and anterior–posterior (A–P) joint load predictions from I–E and A–P motions under a given compressive load, and to evaluate the repeatability of joint load ratios (I–E torque to compressive force (I–E:C), and A–P force to compressive force (A–P:C)) for a range of compressive loading profiles. A tibiofemoral finite element model was developed and used to simulate deep knee bend, chair-rise and step-up activities for five patients. Root-mean-square (RMS) differences in I–E:C and A–P:C load ratios between telemetric measurements and model predictions were less than 1.10e–3 Nm/N and 0.035 N/N for all activities. I–E:C and A–P:C load ratios were consistently reproduced regardless of the compressive force profile applied (RMS differences less than 0.53e–3 Nm/N and 0.010 N/N, respectively). When error in kinematic measurement was introduced to the model, joint load predictions were forgiving to kinematic measurement error when conformity between femoral and tibial components was low. The prevalence of kinematic data, in conjunction with the analysis presented here, facilitates determining the scope of A–P and I–E joint loading ratios experienced by the TKR population.     

Influence of anisotropic bone properties on the biomechanical behavior of the acetabular cup implant: a multiscale finite element study

Although the biomechanical behavior of the acetabular cup (AC) implant is determinant for the surgical success, it remains difficult to be assessed due to the multiscale and anisotropic nature of bone tissue. The aim of the present study was to investigate the influence of the anisotropic properties of peri-implant trabecular bone tissue on the biomechanical behavior of the AC implant at the macroscopic scale. Thirteen bovine trabecular bone samples were imaged using micro-computed tomography (μCT) with a resolution of 18 μm. The anisotropic biomechanical properties of each sample were determined at the scale of the centimeter based on a dedicated method using asymptotic homogenization. The material properties obtained with this multiscale approach were used as input data in a 3D finite element model to simulate the macroscopic mechanical behavior of the AC implant under different loading conditions. The largest stress and strain magnitudes were found around the equatorial rim and in the polar area of the AC implant. All macroscopic stiffness quantities were significantly correlated (R² > 0.85, p < 6.5 e-6) with BV/TV (bone volume/total volume). Moreover, the maximum value of the von Mises stress field was significantly correlated with BV/TV (R² > 0.61, p < 1.6 e-3) and was always found at the bone-implant interface. However, the mean value of the microscopic stress (at the scale of the trabeculae) decrease as a function of BV/TV for vertical and torsional loading and do not depend on BV/TV for horizontal loading. These results highlight the importance of the anisotropic properties of bone tissue.