Why incongruous knee replacements do not fail early

Rapid failure of knee prostheses does not usually occur, despite the non-conforming nature of the articulation between femoral and tibial components and the associated large contact pressures. This theoretical study examines the likelihood of fatigue fracture of a layered elastic model loaded by a sliding cylindrical indenter. Cracks (line-defects) were assumed to have nucleated within the layer. The stress intensity factors (SIFs) associated with these cracks were calculated. The values obtained for the SIFs are quite low, with a corresponding low likelihood of crack-growth. When taken in conjunction with the experimentally derived fatigue laws of previous investigators, they suggest that short line-cracks should not grow. It seems that early failures have not simply been due to large shear stresses which occur beneath the prosthesis surface. Other factors, such as the degradation of material through heat-pressing, sterilisation or oxidation, or the deleterious effect of fusion defects, may be required to drive the cracks to delamination.     

Dynamic contact stress patterns on the tibial plateaus during simulated gait: A novel application of normalized cross correlation

The spatial distribution and pattern of local contact stresses within the knee joint during activities of daily living have not been fully investigated. The objective of this study was to determine if common contact stress patterns exist on the tibial plateaus of human knees during simulated gait. To test this hypothesis, we developed a novel normalized cross-correlation (NCC) algorithm and applied it to the contact stresses on the tibial plateaus of 12 human cadaveric knees subjected to multi-directional loads mimicking gait. The contact stress profiles at different locations on the tibial plateaus were compared, where regions with similar contact stress patterns were identified across specimens. Three consistent regional patterns were found, among them two most prominent contact stress patterns were shared by 9–12 of all the knees and the third pattern was shared by 6–8 knees. The first pattern was located at the posterior aspect of the medial tibial plateau and had a single peak stress that occurred during the early stance phase. The second pattern was located at the central-posterior aspects of the lateral plateau and consisted of two peak stresses coincident with the timing of peak axial force at early and late stance. The third pattern was found on the anterior aspect of cartilage-to-cartilage contact region on the medial plateau consisted of double peak stresses. The differences in the location and profile of the contact stress patterns suggest that the medial and lateral menisci function to carry load at different points in the gait cycle: with the posterior aspect of the medial meniscus consistently distributing load only during the early phase of stance, and the posterior aspect of the lateral meniscus consistently distributing load during both the early and late phases of stance. This novel approach can help identify abnormalities in knee contact mechanics and provide a better understanding of the mechanical pathways leading to post-traumatic osteoarthritis.     

A finite element study of the effect of diametral interface gaps on the contact areas and pressures in uncemented cylindrical femoral total hip components

Uncemented femoral total hip components rely entirely on contact with the prepared femur for their initial fixation. The contact areas and stresses between a straight tubular bone and a metal cylindrical prosthesis 12.5 cm long and 13 mm in diameter were calculated in a finite element model which includes uniform diametral gaps varying from 20 to 500 microns, using transverse loads from 100 to 2000 N. Frictionless three-dimensional contact elements were used between the bone and the prosthesis. Contact stresses were high and irregular in all cases, and the contact areas were small. Two regions of contact were apparent for lower loads and larger gaps. A third region of contact occurred near the distal tip of the implant at higher loads. This region of contact markedly increased the contact stresses at the distal tip of the prosthesis. A 20 microns overlap between bone and implant was modelled to assess a slight interference fit. The contact stress distribution in this case was markedly different from the stress distribution with a 20 microns diametral gap. The data collectively indicates that gaps of less than 20 microns between bone and implant can substantially change contact stress distributions.     

Contact forces in several TKA designs during squatting: A numerical sensitivity analysis

Total knee arthroplasty (TKA) is a very successful procedure, but pain or difficulties during activities still persist in patients. Patient outcomes in TKA surgery can be affected by implant design, alignment or patient-related anatomical factors. This paper presents a numerical sensitivity analysis of several TKA types: a fixed bearing, posterior stabilized prosthesis, a high flexion fixed bearing guided motion prosthesis, a mobile bearing prosthesis and a hinge prosthesis. Each prosthesis was virtually implanted on the same cadaver leg model and it underwent a loaded squat, in 10s, between 0° and 120°, similar to several previous experimental tests performed on knee kinematics simulators. The aim of this examination was to investigate the sensitivity of the patello-femoral (PF) and tibio-femoral (TF) contact forces to patient-related anatomical factors, and component position in the different implant types. The following parameters were used for the sensitivity study: the proximo-distal patellar position, the patellar component tilting, the tibial component position and orientation, the locations of the medial and lateral collateral ligaments with respect to femur and tibia and the patellar tendon length. The sensitivity analysis showed that PF contact forces are mostly affected by patella height (increases up to 67% for one TKA type in patella-alta configuration), by an anterior tibial component translation (increases up to 30%), and by patellar component tilting (increases up to 29%); TF contact forces are mostly affected by the anterior displacement of the insertion points of the medial collateral ligament with respect to the reference position (increases up to 48%).     

The conflicting requirements of laxity and conformity in total knee replacement

Bearing surfaces of total condylar knees which are designed with a high degree of conformity to produce low stresses in the polyethylene tibial insert may be overconstrained. This study determines femoral and tibial bearing surface geometries which will induce the least destructive fatigue mechanisms in the polyethylene whilst conserving the laxity of the natural knee. Sixteen knee designs were generated by varying four parameters systematically to cover the range of contemporary knee designs. The parameters were the femoral frontal radius (30 or 70 mm), the difference between the femoral and tibial frontal radii (2 or 10 mm), the tibial sagittal radius (56 or 80 mm) and the posterior-distal transition angle (-8 or -20 degrees), which is the angle at which the small posterior arc of the sagittal profile transfers to the larger distal arc. Rigid body analyses determined the anterior-posterior and rotational motions as well as the contact points during the stance phase of gait for the different designs. In addition, a damage function which accumulated the fluctuating maximum shear stresses was used to predict the susceptibility to delamination wear of the polyethylene (damage score). This study predicted that of the 16 designs, the knee with a frontal radius of 70 mm, a difference in femoral and tibial frontal radii of 2 mm, a tibial sagittal radius of 80 mm and a posterior distal transition angle of -20 degrees would satisfy the conflicting needs of both resistance to delamination wear and natural kinematics.     

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.     

An improved method for measuring tibiofemoral contact areas in total knee arthroplasty: a comparison of K-scan sensor and Fuji film.

A computerised, real time, thin-film pressure transducer method is used to measure tibiofemoral contact area in total knee arthroplasty (TKA) devices that is easier and more reliable and reproducible as compared to the Fuji pressure-sensitive film technique. Many authors have suggested that contact areas and pressures within TKA devices can be a predictor of wear and failure of the ultra-high molecular weight polyethylene (UHMWPE) tibial insert. In this study, two contact area measurement techniques (Fuji pressure-sensitive Film and K-scan sensor system) were compared using a custom TKA testing jig designed for freedom of movement so that in any loading configuration the component found and seated in its own "home" position. The K-scan system was used to measure contact areas of one TKA design at several angles from 0 to 110 degrees flexion with loads equating to 4, 4.5, and 5 times body weight. For comparison, four ranges of Fuji film were used to measure areas at the same flexion angles but at 5 times body weight only. Contact areas measured with the Fuji films were 11-36% (p < 0.05) lower than those measured by the K-scan sensor.     

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.     

A mathematical model for the evaluation of the behaviour during flexion of condylar-type knee prostheses

A 3D knee model was developed in order to evaluate the mechanical behaviour during flexion of condylar-type knee prosthesis. Based on the total energy minimization principle, it takes into account the articular surfaces (the tibial surface being deformable), the body weight, and the patello femoral joint. It generates the kinematics of the joint, the motion of the centre of contact, the quadriceps forces, the pressure distribution on the tibial plateau, and ligament lengths and forces between 0 and 120 degrees of flexion. The results for ten digitized knees and the commercially available prostheses are presented. They are in general agreement with experimental results published in the literature. It is concluded that this computer program may be, within its limitations, a useful tool in the preliminary evaluation of new condylar-type knee prosthesis 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.     

Experimental and Numerical Analyses of the Pull-out Response of a Steel Post/Bovine Bone Cementless Fixation

Effect of initial interference fit on pull-out strength in cementless fixation between bovine tibia and smooth stainless steel post was investigated in this study.Compressive behavior of bovine spongious bone was studied using mechanical testing in order to evaluate the elastic-plastic properties in different regions of the proximal tibia.Friction tests were carried out in the aim to evaluate the friction behavior of the contact between bovine spongious bone and stainless steel.A cylindrical stainless steel post inserted in a pre-drilled bovine tibia with an initial interference fit was taken as an in vitro model to assess the contribution of post fixation to the initial stability of the Total Knee Arthroplasty (TKA) tibial component.Pull-out experiments were carried out for different initial interference fits.Finite Element Models (FEM) using local elastic-plastic properties of the bovine bone were developed for the analysis of the experimental ultimate pull-out force results.At the post/bone interface,Coulomb friction was considered in the FEM calculations with pressure-dependent friction coefficient.It was found that the FEM results of the ultimate force are in good agreement with the experimental results.The analysis of the FEM interfacial stresses indicates that the micro-slip initiation depends on the local bone properties.     

Morphological and numerical studies of the stress compatibility of patella implants]

Failure of total knee arthroplasty is relatively often caused by problems of the patellofemoral replacement. The purpose of this study was to analyze the distribution of stresses within an anatomical patella and the changes in stress distribution after patellar resurfacing with a Miller-Galante I patellar implant using two- and three dimensional finite element models (FEM). To assess validity, FEM results were compared with morphological findings from contact radiographs and densitographs. Internal orientation of bone trabeculae is in good agreement with the orientation of theoretically calculated principal stresses. Almost unchanged principal tensile stresses after implantation, together with the lack of extreme stress peaks within the cancellous bone ensure stress compatibility of the implant. In the case of a firmly seated implant with good bone ingrowth, increased von Mises stresses are found near the fixation peg/plate junction. Their relevance for improved bone ingrowth near this part of the interface is emphasized. At the same time, material failure at the peg/plate junction can be better understood. An analysis of the early postoperative period assuming nonlinear interface conditions failed to demonstrate an uniform distribution of normal and tangential interface forces.     

Patient-Specific Design and Biomechanical Evaluation of a Novel Bipolar Femoral Hemi-Knee Prosthesis

While total knee replacement is successful, hemiarthroplasty is necessary for some young, obese and active patients who are especially not suitable for unicompartmental or total knee prostheses. Hemiarthroplasty also provides an opportunity for children with bone tumors. The design ofhemiarthroplasty should be patient-specific to reduce contact stress and friction as well as instability, compared to conventional hemi-knee prosthesis. A novel bipolar hemi-knee prosthesis with two flexion stages was developed according to a healthy male＇s knee morphological profile. The motion mode of the bipolar hemi-knee prosthesis was observed through roentgenoscopy in vitro experiment. The biomechanical properties in one gait cycle were evaluated though finite element simulation. The bipolar hemi-knee prosthesis was found to produce knee flexion at two stages through X-ray images. The first stage is the motion from upright posture to a specified 60~ flexion, followed by the second stage of motion subsequently to deep flexion. The finite element simulation results also show that the designed hemi-knee prosthesis has the ability to reduce stresses on the joint contact surfaces. Therefore, it is possible for the bipolar hemi-knee prosthesis to provide better biotribological performances because it can reduce stresses and potentially wear on the opposing contacting surface during a gait cycle, orovidin~ a t~romisin~ treatment strate~v in future Joint renair znd renlneement     

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.     

Effect of superficial collagen patterns and fibrillation of femoral articular cartilage on knee joint mechanics-a 3D finite element analysis

Collagen fibrils of articular cartilage have specific depth-dependent orientations and the fibrils bend in the cartilage surface to exhibit split-lines. Fibrillation of superficial collagen takes place in osteoarthritis. We aimed to investigate the effect of superficial collagen fibril patterns and collagen fibrillation of cartilage on stresses and strains within a knee joint. A 3D finite element model of a knee joint with cartilage and menisci was constructed based on magnetic resonance imaging. The fibril-reinforced poroviscoelastic material properties with depth-dependent collagen orientations and split-line patterns were included in the model. The effects of joint loading on stresses and strains in cartilage with various split-line patterns and medial collagen fibrillation were simulated under axial impact loading of 1000 N. In the model, the collagen fibrils resisted strains along the split-line directions. This increased also stresses along the split-lines. On the contrary, contact and pore pressures were not affected by split-line patterns. Simulated medial osteoarthritis increased tissue strains in both medial and lateral femoral condyles, and contact and pore pressures in the lateral femoral condyle. This study highlights the importance of the collagen fibril organization, especially that indicated by split-line patterns, for the weight-bearing properties of articular cartilage. Osteoarthritic changes of cartilage in the medial femoral condyle created a possible failure point in the lateral femoral condyle. This study provides further evidence on the importance of the collagen fibril organization for the optimal function of articular cartilage.     

A dynamic finite element analysis of impulsive loading of the extension-splinted rabbit knee

A dynamic nonlinear finite element model was developed to study juxtarticular stresses in the splinted rabbit knee, an established laboratory model for creating osteoarthrosis due to impulsive loading. Plane strain finite element results were validated by comparison with corresponding experimental data. Parametric effects studied included the input tibial displacement speed, the local bone density distribution, and the modulus of cartilage and subchondral bone. While the computed resultant contact force magnitude was sensitive to a number of model parameters, the stress patterns, when normalized to a given resultant force magnitude, were not. Despite comparable force peaks, the finite element results showed approximately six-fold higher effective strain rate levels for a severely impulsive loading protocol known to induce rapid osteoarthrosis, versus those for a mildly impulsive loading protocol not usually associated with cartilage damage. A propensity for elevated shear in the deep cartilage layer near the contact periphery, observed in nearly all computed stress distributions, is consistent with previous experimental findings of fissuring at that level in the impulsively loaded rabbit knee.     

The effects of tibial component inclination on bone stress after unicompartmental knee arthroplasty

Unlike the case with total knee arthroplasty, the femorotibial angle (FTA) after unicompartmental knee arthroplasty (UKA) does not directly depend on the inclination of the tibial component when the height of the joint line is maintained. This study analyzed the effects of the inclination of the tibial component in the coronal plane on the contact pressure of the implant-bone surface and the stresses on the proximal tibia. A two-dimensional, coronal plane model of the proximal tibia was subjected to finite-element analysis. Sixteen patterns of finite-element models of equal FTA were developed in which the inclination of tibial components ranged from 5 degrees valgus to 10 degrees varus in increments of 1 degrees. Stress concentration at the proximal medial diaphyseal cortex gradually increased as the inclination changed from valgus to varus. Maximum contact pressure on the metal-bone interface similarly changed and shifted from the lateral edge to the medial edge of the implant as the inclination changed to varus. It was found that even without changing FTA, the inclination of the tibial component might affect stress concentration and contact pressure in the proximal tibia after UKA. The results suggested that slight valgus inclination of the tibial component might be preferable to varus and even to 0 degrees (square) inclination so far as the stress distribution is concerned.     

The effect of the interface on the bone stresses beneath tibial components

It was proposed that the stresses in the layer of bone immediately beneath a tibial component are an important determinant of fixation durability. Using finite element analysis, (ANSYS), the stresses were determined as a function of the amount of bone resection, the localization or completeness of implant-bone contact, and the interface material. The model was of two-dimensional sagittal slices consisting of quadrilateral elements (1 mm) with a range of seventeen material properties determined by CT scans. Typical prosthesis designs shifted the center of pressure more centrally rather than posteriorly, and thus increased anterior bone stresses. Resection up to 10 mm could actually decrease bone stresses due to an increase in bone surface area as long as complete coverage was obtained. A cement interface or direct metal on bone produced identical stresses. However a 1 mm complian: interface significantly reduced stresses in regions of high elastic modulus gradient. For rigid interfaces, the contact can be irregular, which leads to areas of over and understressing of bone. These conclusions have implications related to implant design.     

Discrete element analysis in musculoskeletal biomechanics

This paper is written to honor Professor Y. C. Fung, the applied mechanician who has made seminal contributions in biomechanics. His work has generated great spin-off utility in the field of musculoskeletal biomechanics. Following the concept of the Rigid Body-Spring Model theory by T. Kawai (1978) for non-linear analysis of beam, plate, and shell structures and the soil-gravel mixture foundation, we have derived a generalized Discrete Element Analysis (DEA) method to determine human articular joint contact pressure, constraining ligament tension and bone-implant interface stresses. The basic formulation of DEA to solve linear problems is reviewed. The derivation of non-linear springs for the cartilage in normal diarthrodial joint contact problem was briefly summarized. Numerical implementation of the DEA method for both linear and non-linear springs is presented. This method was able to generate comparable results to the classic contact stress problem (the Hertzian solution) and the use of Finite Element Modeling (FEM) technique on selected models. Selected applications in human knee and hip joints are demonstrated. In addition, the femoral joint prosthesis stem/bone interface stresses in a non-cemented fixation were analyzed using a 2D plane-strain approach. The DEA method has the advantages of ease in creating the model and reducing computational time for joints of irregular geometry. However, for the analysis of joint tissue stresses, the FEA technique remains the method of choice.     

Computational wear prediction of artificial knee joints based on a new wear law and formulation

Laboratory joint wear simulator testing has become the standard means for preclinical evaluation of wear resistance of artificial knee joints. Recent simulator designs have been advanced and become successful at reproducing the wear patterns observed in clinical retrievals. However, a single simulator test can be very expensive and take a long time to run. On the other hand computational wear modelling is an alternative attractive solution to these limitations. Computational models have been used extensively for wear prediction and optimisation of artificial knee designs. However, all these models have adopted the classical Archard's wear law, which was developed for metallic materials, and have selected wear factors arbitrarily. It is known that such an approach is not generally true for polymeric bearing materials and is difficult to implement due to the high dependence of the wear factor on the contact pressure. Therefore, these studies are generally not independent and lack general predictability. The objective of the present study was to develop a new computational wear model for the knee implants, based on the contact area and an independent experimentally determined non-dimensional wear coefficient. The effects of cross-shear and creep on wear predictions were also considered. The predicted wear volume was compared with the laboratory simulation measurements. The model was run under two different kinematic inputs and two different insert designs with curved and custom designed flat bearing surfaces. The new wear model was shown to be capable of predicting the difference of the wear volume and wear pattern between the two kinematic inputs and the two tibial insert designs. Conversely, the wear factor based approach did not predict such differences. The good agreement found between the computational and experimental results, on both the wear scar areas and volumetric wear rates, suggests that the computational wear modelling based on the new wear law and the experimentally calculated non-dimensional wear coefficient should be more reliable and therefore provide a more robust virtual modelling platform.