Contact stresses and fatigue life in a knee prosthesis: comparison between in vitro measurements and computational simulations

The evaluation of contact areas and pressures in total knee prosthesis is a key issue to prevent early failure. The first part of this study is based on the hypothesis that the patterns of contact stresses on the tibial insert of a knee prosthesis at different stages of the gait cycle could be an indicator of the wear performances of a knee prosthesis. Contact stresses were calculated for a mobile bearing knee prosthesis by means of finite element method (FEM). Contact areas and stresses were also measured through in vitro tests using Fuji Prescale film in order to support the FEM findings.The second part of this study addresses the long-term structural integrity of metal tibial components in terms of fatigue life by means of experimental tests and FEM simulations. Fatigue experimental evaluations were performed on Cr-Co alloy tibial tray, based on ISO standards. FEM models were used to calculate the stress patterns. The failure risk was estimated with a standard fatigue criterion on the basis of the results obtained from the FEM calculations. Experimental and computational results showed a positive matching.     

Analysis of bone–prosthesis interface micromotion for cementless tibial prosthesis fixation and the influence of loading conditions

A lack of initial stability of the fixation is associated with aseptic loosening of the tibial components of cementless knee prostheses. With sufficient stability after surgery, minimal relative motion between the prosthesis and bone interfaces allows osseointegation to occur thereby providing a strong prosthesis-to-bone biological attachment. Finite element modelling was used to investigate the bone–prosthesis interface micromotion and the relative risk of aseptic loosening. It was anticipated that by prescribing different joint loads representing gait and other activities, and the consideration of varying tibial–femoral contact points during knee flexion, it would influence the computational prediction of the interface micromotion. In this study, three-dimensional finite element models were set up with applied loads representing walking and stair climbing, and the relative micromotions were predicted. These results were correlated to in-vitro measurements and to the results of prior retrieval studies. Two load conditions, (i) a generic vertical joint load of 3×body weight with 70%/30% M/L load share and antero-posterior/medial-lateral shear forces, acted at the centres of the medial and lateral compartments of the tibial tray, and (ii) a peak vertical joint load at 25% of the stair climbing cycle with corresponding antero-posterior shear force applied at the tibial–femoral contact points of the specific knee flexion angle, were found to generate interface micromotion responses which corresponded to in-vivo observations. The study also found that different loads altered the interface micromotion predicted, so caution is needed when comparing the fixation performance of various reported cementless tibial prosthetic designs if each design was evaluated with a different loading condition.     

Simulation of a knee joint replacement during a gait cycle using explicit finite element analysis.

The stress distribution within the polyethylene insert of a total knee joint replacement is dependent on the kinematics, which in turn are dependent on the design of the articulating surfaces, the relative position of the components and the tension of the surrounding soft tissues. Implicit finite element analysis techniques have been used previously to examine the polyethylene stresses. However, these have essentially been static analyses and hence ignored the influence of the kinematics. The aim of this work was to use an explicit finite element approach to simulate both the kinematics and the internal stresses within a single analysis. A simulation of a total knee joint replacement subjected to a single gait cycle within a knee wear simulator was performed and the results were compared with experimental data.The predicted kinematics were in close agreement with the experimental data. Various solution-dependent parameters were found to have little influence on the predicted kinematics. The predicted stresses were found to be dependent on the mesh density. This study has shown that an explicit finite element approach is capable of predicting the kinematics and the stresses within a single analysis at relatively low computational cost.     

Protocol for constructing subject-specific biomechanical models of knee joint

A robust protocol for building subject-specific biomechanical models of the human knee joint is proposed which uses magnetic resonance imaging, motion analysis and force platform data in conjunction with detailed 3D finite element models. The proposed protocol can be used for determining stress and strain distributions and contact kinetics in different knee elements at different body postures during various physical activities. Several examples are provided to highlight the capabilities and potential applications of the proposed protocol. This includes preliminary results on the role of body weight on the stresses and strains induced in the knee articular cartilages and meniscus during single-leg stance and calculations of the induced stresses and ligament forces during the gait cycle.     

Tibiofemoral joint contact force in deep knee flexion and its consideration in knee osteoarthritis and joint replacement

The aim of the study was to estimate the tibiofemoral joint force in deep flexion to consider how the mechanical load affects the knee. We hypothesize that the joint force should not become sufficiently large to damage the joint under normal contact area, but should become deleterious to the joint under the limited contact area. Sixteen healthy knees were analyzed using a motion capture system, a force plate, a surface electromyography, and a knee model, and then tibiofemoral joint contact forces were calculated. Also, a contact stress simulation using the contact areas from the literature was performed. The peak joint contact forces (M +/- SD) were 4566 +/- 1932 N at 140 degrees in rising from full squat and 4479 +/- 1478 N at 90 degrees in rising from kneeling. Under normal contact area, the tibiofemoral contact stresses in deep flexion were less than 5 MPa and did not exceed the stress to damage the cartilage. The contact stress simulation suggests that knee prosthesis having the contact area smaller than 200 mm2 may be problematic since the contact stress in deep flexion would become larger than 21 MPa, and it would lead damage or wear of the polyethylene.     

Dynamic contact mechanics on the tibial plateau of the human knee during activities of daily living

Despite significant advances in scaffold design, manufacture, and development, it remains unclear what forces these scaffolds must withstand when implanted into the heavily loaded environment of the knee joint. The objective of this study was to fully quantify the dynamic contact mechanics across the tibial plateau of the human knee joint during gait and stair climbing. Our model consisted of a modified Stanmore knee simulator (to apply multi-directional dynamic forces), a two-camera motion capture system (to record joint kinematics), an electronic sensor (to record contact stresses on the tibial plateau), and a suite of post-processing algorithms. During gait, peak contact stresses on the medial plateau occurred in areas of cartilage–cartilage contact; while during stair climb, peak contact stresses were located in the posterior aspect of the plateau, under the meniscus. On the lateral plateau, during gait and in early stair-climb, peak contact stresses occurred under the meniscus, while in late stair-climb, peak contact stresses were experienced in the zone of cartilage–cartilage contact. At 45% of the gait cycle, and 20% and 48% of the stair-climb cycle, peak stresses were simultaneously experienced on both the medial and lateral compartment, suggesting that these phases of loading warrant particular consideration in any simulation intended to evaluate scaffold performance. Our study suggests that in order to design a scaffold capable of restoring ‘normal’ contact mechanics to the injured knees, the mechanics of the intended site of implantation should be taken into account in any pre-clinical testing regime.     

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.     

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.     

The effects of gastrocnemius–soleus muscle forces on ankle biomechanics during triple arthrodesis

This paper presents a finite element model of the ankle, taking into account the effects of muscle forces, determined by a musculoskeletal analysis, to investigate the contact stress distribution in the tibio-talar joint in patients with triple arthrodesis and in normal subjects. Forces of major ankle muscles were simulated and corresponded well with the trend of their EMG signals. These forces were applied to the finite element model to obtain stress distributions for patients with triple arthrodesis and normal subjects in three stages of the gait cycle, i.e. heel strike, midstance, and heel rise. The results demonstrated that the stress distribution patterns of the tibio-talar joint in patients with triple arthrodesis differ from those of normal subjects in investigated gait cycle stages. The mean and standard deviations for maximum stresses in the tibo-talar joint in the stance phase for patients and normal subjects were 9.398e7 ± 1.75e7 and 7.372e7 ± 4.43e6 Pa, respectively. The maximum von Mises stresses of the tibio-talar joint for all subjects in the stance phase found to be on the lateral side of the inferior surface of the joint. The results also indicate that, in patients with triple arthrodesis, increasing gastrocnemius–soleus muscle force reduces the stress on the medial malleolus compared with normal subjects. Most of stresses in this area are between 45 and 109 kPa, and will decrease to almost 32 kPa in patients after increasing of 40% in gastrocnemius–soleus muscle force.     

Development of an above-knee prosthesis equipped with a microcomputer-controlled knee joint: first test results

The shortcomings of conventional above-knee prostheses are due to their lack of adaptive control. Implementation of a microcomputer controlling the knee joint in a passive way has been suggested to enhance the patient's gait comfort, safety and cosmesis. This approach was used in the design of a new prosthetic system for the above-knee amputee, and tested on one patient. The knee joint of a conventional, modular prosthesis was replaced by a knee joint mechanism, equipped with a controllable brake on the knee joint axis. Sensors and a microcomputer were added, keeping the system self-contained. The modularity of the design permits the use of an alternative, external, PC-based control unit, emulating the self-contained one, and offering extended data monitoring and storage facilities. For both units an operating environment was written, including sensor/actuator interfacing and the implementation of a real-time interrupt, executing the control algorithm. A double finite state approach was used in the design of the control algorithm. On a higher level, the mode identification algorithm reveals the patient's intent. Within a specific mode (lower level), the relevant mode control algorithm looks for the current phase within the gait cycle. Within a particular phase, a specific simple control action with the brake replaces normal knee muscle activity. Test were carried out with one prosthetic patient using a basic control algorithm for level walking, allowing controlled knee flexion during stance phase. The technical feasibility of such a concept is illustrated by the test results, even though only flexion during early stance phase was controlled during the trials. Patient acceptance is not straightforward since knee flexion during stance phase is associated with knee buckling.     

Image based weighted center of proximity versus directly measured knee contact location during simulated gait

To understand the mechanical consequences of knee injury requires a detailed analysis of the effect of that injury on joint contact mechanics during activities of daily living. Three-dimensional (3D) knee joint geometric models have been combined with knee joint kinematics to dynamically estimate the location of joint contact during physiological activities—using a weighted center of proximity (WCoP) method. However, the relationship between the estimated WCoP and the actual location of contact has not been defined. The objective of this study was to assess the relationship between knee joint contact location as estimated using the image-based WCoP method, and a directly measured weighted center of contact (WCoC) method during simulated walking. To achieve this goal, we created knee specific models of six human cadaveric knees from magnetic resonance imaging. All knees were then subjected to physiological loads on a knee simulator intended to mimic gait. Knee joint motion was captured using a motion capture system. Knee joint contact stresses were synchronously recorded using a thin electronic sensor throughout gait, and used to compute WCoC for the medial and lateral plateaus of each knee. WCoP was calculated by combining knee kinematics with the MRI-based knee specific model. Both metrics were compared throughout gait using linear regression. The anteroposterior (AP) location of WCoP was significantly correlated with that of WCoC on both tibial plateaus in all specimens (p<0.01, 95% confidence interval of Pearson?s coefficient r>0), but the correlation was not significant in the mediolateral (ML) direction for 4/6 knees (p>0.05). Our study demonstrates that while the location of joint contact obtained from 3D knee joint contact model, using the WCoP method, is significantly correlated with the location of actual contact stresses in the AP direction, that relationship is less certain in the ML direction.     

Thigh–calf contact: Does it affect the loading of the knee in the high-flexion range?

Recently, high-flexion knee implants have been developed to provide for a large range of motion (ROM>120°) after total knee arthroplasty (TKA). Since knee forces typically increase with larger flexion angles, it is commonly assumed that high-flexion knee implants are subjected to larger loads than conventional knee implants. However, most high-flexion studies do not consider thigh–calf contact which occurs during high-flexion activities such as squatting and kneeling. In this study, we hypothesized that thigh–calf contact reduces the knee forces during deep knee flexion as the tibio-femoral load shifts from occurring inside the knee towards the thigh–calf contact interface. Hence, the effect of thigh–calf contact on the knee loading was evaluated using a free body diagram and a finite element model and both the knee forces and polyethylene stresses were analyzed. Thigh–calf contact force characteristics from an earlier study were included and a squatting movement was simulated. In general, we found thigh–calf contact considerably reduced both the knee forces and polyethylene stresses during deep knee flexion. At maximal flexion (155°), the compressive knee force decreased from 4.89 to 2.90 times the bodyweight (BW) in case thigh–calf contact was included and the polyethylene contact stress at the tibial post decreased from 49.3 to 28.1 MPa. Additionally, there was a clear correlation between a subject's thigh and calf circumference and the force reduction at maximal flexion due to thigh–calf contact (R=0.89). The findings presented in this study can be used to optimize the mechanical behavior of high-flexion total knee arthroplasty designs.     

The load-bearing area in the knee joint

Measurements were made of the location and size of the contact areas in cadaver knee joints, for a load of 150 Kgf applied for 5 sec down the long axis of the tibia. Results were obtained from a total of 4 knees, considering flexion angles from 0 to 120°. The methods used were to measure directly from castings of the joint cavity; and to calculate from measurements of radii of curvature and joint deflection. Average contact areas for lateral and medial condyles were 1·4 and 1·8 cm² respectively. Areas for the medial condyle were greater than for the lateral condyle and also the areas diminished as flexion angle increased. The implications of the results to contact stresses, joint lubrication and ‘condylar replacement’ knee prosthesis design were discussed.     

On the estimation of joint kinematics during gait.

In gait analysis, the concepts of Euler and helical (screw) angles are used to define the three-dimensional relative joint angular motion of lower extremities. Reliable estimation of joint angular motion depends on the accurate definition and construction of embedded axes within each body segment. In this paper, using sensitivity analysis, we quantify the effects of uncertainties in the definition and construction of embedded axes on the estimation of joint angular motion during gait. Using representative hip and knee motion data from normal subjects and cerebral palsy patients, the flexion-extension axis is analytically perturbed +/- 15 degrees in 5 degrees steps from a reference position, and the joint angles are recomputed for both Euler and helical angle definitions. For the Euler model, hip and knee flexion angles are relatively unaffected while the ab/adduction and rotation angles are significantly affected throughout the gait cycle. An error of 15 degrees in the definition of flexion-extension axis gives rise to maximum errors of 8 and 12 degrees for the ab/adduction angle, and 10-15 degrees for the rotation angles at the hip and knee, respectively. Furthermore, the magnitude of errors in ab/adduction and rotation angles are a function of the flexion angle. The errors for the ab/adduction angles increase with increasing flexion angle and for the rotation angle, decrease with increasing flexion angle. In cerebral palsy patients with flexed knee pattern of gait, this will result in distorted estimation of ab/adduction and rotation. For the helical model, similar results are obtained for the helical angle and associated direction cosines.(ABSTRACT TRUNCATED AT 250 WORDS)     

A subject-specific finite element musculoskeletal framework for mechanics analysis of a total knee replacement

Concurrent use of finite element (FE) and musculoskeletal (MS) modeling techniques is capable of considering the interactions between prosthetic mechanics and subject dynamics after a total knee replacement (TKR) surgery is performed. However, it still has not been performed in terms of favorable prediction accuracy and systematic experimental validation. In this study, we presented a methodology to develop a subject-specific FE-MS model of a human right lower extremity including the interactions among the subject-specific MS model, the knee joint model with ligament bundles, and the deformable FE prosthesis model. In order to evaluate its accuracy, the FE-MS model was compared with a traditional hinge-constraint MS model and experimentally verified over a gait cycle. Both models achieved good temporal agreement between the predicted muscle force and the electromyography results, though the magnitude on models is different. A higher predicted accuracy, quantified by the root-mean-square error (RMSE) and the squared Pearson correlation coefficient (r²), was found in the FE-MS model (RMSE = 177.2 N, r² = 0.90) when compared with the MS model (RMSE = 224.1 N, r² = 0.81) on the total tibiofemoral contact force. The contact mechanics, including the contact area, pressure, and stress were synchronously simulated, and the maximum contact pressure, 22.06 MPa, occurred on the medial side of the tibial insert without exceeding the yield strength of the ultra-high-molecular-weight polyethylene, 24.79 MPa. The approach outlines an accurate knee joint biomechanics analysis and provides an effective method of applying individualized prosthesis design and verification in TKR.     

Stress and micromotion in the taper lock joint of a modular segmental bone replacement prosthesis

The stress distribution within the components and the micromotion of the interface significantly influence the long-term function of the taper lock joint in a modular segmental bone replacement prosthesis. Bending-induced gap opening between the cone and the sleeve can lead to an inflow of biological fluids, and thus accelerate implant corrosion. Local areas of high stress can also accelerate the corrosive processes and initiate local yielding, which may lead to a fracture in one of the components. In this study, a 3-D finite element (FE) model of a modular segmental bone replacement prosthesis was developed to study the interface micromotion and component stress distribution under the maximum loads applied during gait for a taper lock joint with multiple material combinations. Bending was the main cause of the local high stresses and interface separation within the taper joint. For Ti6A14V components, cortical bone bridging and ingrowth across the taper lock gap reduced the peak stress by 45% and reduced the contact interface separation by 55%. Such tissue formation around the taper lock joint could also form a closed capsule to restrict the migration of potential wear particles and thus prevent the biologic process of bone resorption induced by metal debris.     

The effect of meniscal tears and resultant partial meniscectomies on the knee contact stresses: a finite element analysis

Knee osteoarthritis (OA) is believed to result from high levels of contact stresses on the articular cartilage and meniscus after meniscal damage. This study investigated the effect of meniscal tears and partial meniscectomies on the peak compressive and shear stresses in the human knee joint. An elaborate three-dimensional finite element model of knee joint including bones, articular cartilages, menisci and main ligaments was developed from computed tomography and magnetic resonance imaging images. This model was used to model four types of meniscal tears and their resultant partial meniscectomies and analysed under an axial 1150 N load at 0° flexion. Three different conditions were compared: a healthy knee joint, a knee joint with medial meniscal tears and a knee joint following partial meniscectomies. The numerical results showed that each meniscal tear and its resultant partial meniscectomy led to an increase in the peak compressive and shear stresses on the articular cartilages and meniscus in the medial knee compartment, especially for partial meniscectomy. Among the four types of meniscal tears, the oblique tear resulted in the highest values of the peak compressive and shear stresses. For the four partial meniscectomies, longitudinal meniscectomy led to the largest increase in these two stresses. The lateral compartment was minimally affected by all the simulations. The results of this study demonstrate meniscal tear and its resultant partial meniscectomy has a positive impact on the maintenance of high levels of contact stresses, which may improve the progression of knee OA, especially for partial meniscectomy. Surgeons should adopt a prudent strategy to preserve the greatest amount of meniscus possible.     

Tibiofemoral kinematics and condylar motion during the stance phase of gait

Accurate knowledge of the dynamic knee motion in-vivo is instrumental for understanding normal and pathological function of the knee joint. However, interpreting motion of the knee joint during gait in other than the sagittal plane remains controversial. In this study, we utilized the dual fluoroscopic imaging technique to investigate the six-degree-of-freedom kinematics and condylar motion of the knee during the stance phase of treadmill gait in eight healthy volunteers at a speed of 0.67 m/s. We hypothesized that the 6DOF knee kinematics measured during gait will be different from those reported for non-weightbearing activities, especially with regards to the phenomenon of femoral rollback. In addition, we hypothesized that motion of the medial femoral condyle in the transverse plane is greater than that of the lateral femoral condyle during the stance phase of treadmill gait. The rotational motion and the anterior–posterior translation of the femur with respect to the tibia showed a clear relationship with the flexion–extension path of the knee during the stance phase. Additionally, we observed that the phenomenon of femoral rollback was reversed, with the femur noted to move posteriorly with extension and anteriorly with flexion. Furthermore, we noted that motion of the medial femoral condyle in the transverse plane was greater than that of the lateral femoral condyle during the stance phase of gait (17.4±2.0 mm vs. 7.4±6.1 mm, respectively; p<0.01). The trend was opposite to what has been observed during non-weightbearing flexion or single-leg lunge in previous studies. These data provide baseline knowledge for the understanding of normal physiology and for the analysis of pathological function of the knee joint during walking. These findings further demonstrate that knee kinematics is activity-dependent and motion patterns of one activity (non-weightbearing flexion or lunge) cannot be generalized to interpret a different one (gait).     

Knee joint biomechanics in closed-kinetic-chain exercises

Effective management of knee joint disorders demands appropriate rehabilitation programs to restore function while strengthening muscles. Excessive stresses in cartilage/menisci and forces in ligaments should be avoided to not exacerbate joint condition after an injury or reconstruction. Using a validated 3D nonlinear finite element model, detailed biomechanics of the entire joint in closed-kinetic-chain squat exercises are investigated at different flexion angles, weights in hands, femur-tibia orientations and coactivity in hamstrings. Predictions are in agreement with results of earlier studies. Estimation of small forces in cruciate ligaments advocates the use of squat exercises at all joint angles and external loads. In contrast, large contact stresses, especially at the patellofemoral joint, that approach cartilage failure threshold in compression suggest avoiding squatting at greater flexion angles, joint moments and weights in hands. Current results are helpful in comprehensive evaluation and design of effective exercise therapies and trainings with minimal risk to various components.     

Prediction of total knee motion using a three-dimensional computer-graphics model

Twenty-three knees were sectioned, digitized, and standardized to determine the 'average' three-dimensional bony geometry and ligamentous attachments. Data on normal knee motion were obtained from a cadaveric study. An algorithm was written to simulate three-dimensional patella motion. Verification of the knee model was achieved by determining femoro-tibial and patello-femoral contact locations, as well as ligament length patterns, and comparing the results with published data. The criterion for maximum predicted knee motion with a prosthesis in place was the length of the posterior cruciate ligament. Three total knee replacement surfaces were mathematically generated: flat, laxity and conforming. A greater flexion angle was obtained with a flat tibial surface than for the laxity or conforming. Posterior tibial component displacement increased the range of motion, but only slightly. For all tibial surfaces, increased range of motion was achieved with a 10 degrees posterior tilt of the tibial tray. Anterior femoral component displacement increased motion due to reduction in posterior cruciate tension during flexion. The results are applicable to the design and surgical technique of total knee replacement.