Contact mechanics and elastohydrodynamic lubrication in a novel metal-on-metal hip implant with an aspherical bearing surface

Diameter and diametral clearance of the bearing surfaces of metal-on-metal hip implants and structural supports have been recognised as key factors to reduce the dry contact and hydrodynamic pressures and improve lubrication performance. On the other hand, application of aspherical bearing surfaces can also significantly affect the contact mechanics and lubrication performance by changing the radius of the curvature of a bearing surface and consequently improving the conformity between the head and the cup. In this study, a novel metal-on-metal hip implant employing a specific aspherical bearing surface, Alpharabola, as the acetabular surface was investigated for both contact mechanics and elastohydrodynamic lubrication under steady-state conditions. When compared with conventional spherical bearing surfaces, a more uniform pressure distribution and a thicker lubricant film thickness within the loaded conjunction were predicted for this novel Alpharabola hip implant. The effects of the geometric parameters of this novel acetabular surface on the pressure distribution and lubricant thickness were investigated. A significant increase in the predicted lubricant film thickness and a significant decrease in the dry contact and hydrodynamic pressures were found with appropriate combinations of these geometric parameters, compared with the spherical bearing surface.     

Contact of dual mobility implants: effects of cup wear and inclination

Cup wear and inclination on the pelvic bone are significant factors, which change the contact of the articulating surfaces, thus, impacting the long-term performance of hip implants. This paper presents a finite element (FE) analysis of the contact of the dual mobility implants under the influence of cup wear and inclination. A 3D FE model of the implant was developed with the application of equivalent physiological loading and boundary conditions. Effects of cup inclination angle ranging from 45° to 60° and the wear depth ranging from 0 to 2.46 mm equivalent to up to 30 years of the implant's life on the contact pressure and von Mises stress were investigated. Simulation results show that the contact pressure and von Mises stress decrease significantly with a modest wear depth and remains quite in-sensitive to the cup inclination angle and wear depth up to 1.64 mm. With wear depth further up to 2.46 mm, the cup thickness (i.e. cup thinning on worn region) may be more predominant than increasing of contact area between the cup and the head. The wear on the inner surface of the cup is found to rule out the overall contact pressure and stress in the implant. Furthermore, individual and combined effects of both important parameters are analysed and discussed with respect to available clinical/laboratory studies.     

Finite Element Analysis of the Contact Mechanics of Ceramic-on-Ceramic Hip Resurfacing Prostheses

<正> Ceramics are good alternative to metal as bearing couple materials because of their better wear resistance. A Finite Element(FE) study was performed to investigate the contact mechanics and stress distribution of Ceramic-on-Ceramic (COC) hip resurfacingprostheses. It was focused in particular on a parametric study to examine the effects of radial clearance, loading,alumina coating on the implants, bone quality, and fixation of cup-bone interface. It was found that a reduction in the radialclearance had the most significant effect on the predicted contact pressure distribution among all of the parameters considered inthis study. It was determined that there was a significant influence of non-metallic materials, such as the bone underneath thebearing components, on the predicted contact mechanics. Stress shielding within the bone tissue was found to be a major concernwhen regarding the use of ceramic as an alternative to metallic resurfacing prostheses. Therefore, using alumina implantswith a metal backing was found to be the best design for ceramic resurfacing prostheses in this study. The loading, bone quality,and acetabular cup fixation conditions were found to have only minor effects on the predicted contact pressure distribution alongthe bearing surfaces.     

Prediction of lubricating film thickness in UHMWPE hip joint replacements

An elastohydrodynamic lubrication model developed for a ball-in-socket configuration in a previous studies by the present authors (Jalali-Vahid et al., Thinning films and tribological interfaces, 26th Leeds-Lyon Symposium on Tribology, 2000, pp. 329-339) was applied to analyse the lubrication problem of a typical artificial hip joint replacement, consisting of an ultra-high molecular weight polyethylene (UHMWPE) acetabular cup against a metallic or ceramic femoral head. The cup was assumed to be stationary whilst the ball was assumed to rotate at a steady angular velocity and under a constant load. A wide range of main design parameters were considered. It has been found that the predicted lubricating film thickness increases with a decrease in the radial clearance, an increase in the femoral head radius, an increase in UHMWPE thickness and a decrease in UHMWPE modulus. However, the predicted lubricating film thicknesses are not found to be sufficiently large in relation to the surface roughness of the cup and head to indicate separation of the two articulating surfaces. It should also be noted that if the design features are unable to secure full fluid film lubrication, it may be preferable to select them for minimum wear rather than maximum film thickness. For example, an increase in head radius will enhance the film thickness, but it will also increase the sliding distance and hence wear in mixed or boundary lubrication conditions. Furthermore, it is pointed out that an increase in the predicted lubricant film thickness is usually associated with an increase in the contact area, and this may cause lubricant starvation and stress concentration at the edge of the cup, and adversely affect the tribological performance of the implant. The effect of running-in process on the lubrication in UHMWPE hip joint replacements is also discussed.     

Study of micromotion in modular acetabular components during gait and subluxation: a finite element investigation

Polyethylene wear after total hip arthroplasty may occur as a result of normal gait and as a result of subluxation and relocation with impact. Relocation of a subluxed hip may impart a moment to the cup creating sliding as well as compression at the cup liner interface. The purpose of the current study is to quantify, by a validated finite element model, the forces generated in a hip arthroplasty as a result of subluxation relocation and compare them to the forces generated during normal gait. The micromotion between the liner and acetabular shell was quantified by computing the sliding track and the deformation at several points of the interface. A finite element analysis of polyethylene liner stress and liner/cup micromotion in total hip arthroplasty was performed under two dynamic profiles. The first profile was a gait loading profile simulating the force vectors developed in the hip arthroplasty during normal gait. The second profile is generated during subluxation and subsequent relocation of the femoral head. The forces generated by subluxation relocation of a total hip arthroplasty can exceed those forces generated during normal gait. The induced micromotion at the cup polyethylene interface as a result of subluxation can exceed micromotion as a result of the normal gait cycle. This may play a significant role in the generation of backsided wear. Minimizing joint subluxation by restoring balance to the hip joint after arthroplasty should be explored as a strategy to minimize backsided wear.     

Finite element analysis of a total ankle replacement during the stance phase of gait

Total ankle replacement (TAR) designs have still several important issues to be addressed before the treatment becomes fully acceptable clinically. Very little is known about the performance, in terms of the contact pressures and kinematics of TAR when subjected to daily activities such as level gait. For this purpose, an explicit finite element model of a novel 3-component TAR was developed, which incorporated a previously validated mechanical model of the ankle ligament apparatus. The intermediate mobile polyethylene meniscal bearing was modelled as an elastic-plastic continuum while the articulating surfaces of the tibial and talar metal components as rigid bodies. Overall kinematics, contact pressures and ligament forces were analysed during passive, i.e. virtually unloaded, and active, i.e. stance phase of gait, conditions. Simulation of passive motion predicted similar kinematics as reported previously in an analytical four-bar linkage model. The meniscal bearing was observed to move 5.6 mm posteriorly during the simulated stance and the corresponding antero-posterior displacement of the talar component was 8.3 mm. The predicted pattern and the amount (10.6 degrees ) of internal-external rotation of the ankle complex were found to be in good agreement with corresponding in vivo measurements on normal ankles. A peak contact pressure of 16.8 MPa was observed, with majority of contact pressures below 10 MPa. For most ligaments, reaction forces remain within corresponding physiological ranges. A first realistic representation of the biomechanical behaviour of the human ankle when replaced by prosthetic joints is provided. The applied methodology can potentially be applied to other TAR designs.     

Elastohydrodynamic lubrication analysis of metal-on-metal hip-resurfacing prostheses

The elastohydrodynamic lubrication analysis was carried out in this study for a typical metal-on-metal hip-resurfacing prosthesis under a simple steady-state rotation. Both the Reynolds equation and the elasticity equation were coupled and solved numerically by the finite difference method. The finite element method was used to determine the elastic deformation of both the femoral and the acetabular components required for the lubrication analysis. The effect of the radial clearance between the femoral head and the acetabular cup on the predicted film thickness and pressure distribution was investigated. The predicted minimum lubricating film thickness was found to compare favourably with the prediction using the Hamrock and Dowson [J. Lubrication Technol. 100 (1978) 236] formula based on the assumption of ball-on-plane semi-infinite solids. This implies that the non-metallic materials such as bone and cement underlying the metallic components have a small effect on the predicted lubrication performance for the particular metal-on-metal hip-resurfacing prosthesis considered in this study. Under realistic physiological walking conditions, a decrease in the radial clearance from 150 to 50 microm resulted in a 137% increase in the predicted minimum film thickness from 19 to 45 nm. Consequently, given a surface roughness of 0.01 microm for both the metallic femoral and acetabular bearing surfaces, the predicted mixed lubrication regime for the larger clearance was changed to a full fluid film lubrication regime for the smaller clearance. This clearly highlighted the importance of the design and manufacturing parameters on the tribological performance of these hard-on-hard hip prostheses.     

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.     

Mathematical modelling of stress in the hip during gait.

A mathematical model is developed for calculating the contact stress distribution in the hip for a known resultant hip force and characteristic geometrical parameters. Using a relatively simple single nonlinear algebraic equation, the model can be readily applied in clinical practice to estimate the stress distribution in the most frequent body positions of everyday activities. This is demonstrated by analyzing the data on the resultant hip force obtained from laboratory observations where a stance period of gait is considered.     

An elasto-plastic finite element model for polyethylene wear in total hip arthroplasty

A new finite element model (FEM) based on an elasto-plastic behavior of ultra high molecular weight polyethylene (UHMWPE) was used to study the wear behavior of UHMWPE acetabular cup, which has a 32 mm diameter femoral head. The model imposed a plastic yield stress of 8 MPa on the UHMWPE so that any stresses beyond this would automatically be redistributed to its neighbor. The FEM model adopted a unique mesh design based on an open cube concept which eliminated the problems of singularities. Wear prediction combined the influences of contact stress, sliding distance and a surface wear coefficient. The new model predicted significantly higher volumetric wear rate (57 mm(3)/yr) well within the average reported clinical values. The model was also used to study the effect of friction and clearance between the acetabular cup and the femoral head. Increase in friction increased the volumetric wear rate but did not appear to affect the linear wear rate, which remained at 0.12 +/- 0.02 mm/yr. The predicted wear was sensitive to clearance. It was found that when the clearance was close to 0 and >0.5mm, severe wear occurred. The best clearance range was between 0.1 and 0.15 mm where the average linear wear rate was 0.1mm/yr and the volumetric wear was 55 mm(3)/yr. The present work indicates the importance of avoiding too tight or too loose a diametrical clearance.     

Transient elastohydrodynamic lubrication analysis of metal-on-metal hip implant under simulated walking conditions

The transient elastohydrodynamic lubrication (EHL) analysis was performed in this study for a typical metal-on-metal bearing employing a polyethylene backing underneath a metallic cup inlay under dynamic operating conditions of load and speed representative of normal walking. A ball-in-socket configuration was adopted to represent the articulation between the femoral head and the acetabular cup. The governing Reynolds and elasticity equations were solved simultaneously by using both finite difference and finite element methods. The predicted transient film thickness from the present study was compared with the estimation based on the quasi-static analysis. It was found that the polyethylene backing employed in the typical metal-on-metal hip bearing, combined with dynamic squeeze-film action, significantly improved the transient lubricant film thickness under cyclic walking and consequently a fluid film lubrication regime was possible for smooth bearing surfaces with an average roughness less than 0.005 microm.     

Ceramic cups for hip endoprostheses. 6: Cup design, inclination and antetorsion angle modify range of motion and impingement]

The present investigation focuses on total hip replacement using ceramic acetabular components. The relationship between the position of the cup and the range of motion (ROM) was investigated. A limited range of motion may cause impingement, which is defined as contact between the femoral neck and the rim of the acetabular cup. Impingement may result in wear, chipping, fracture or dislocation of the femoral head. Joint movements were simulated in a three-dimensional CAD program. The results obtained underscore the importance of correct positioning and design of the cup for achieving a ROM as close to the physiological situation as possible. With ceramic cups, the inclination angle should not be more than 45 degrees, and the antetorsion angle between 10 and 15 degrees. If the cup is too vertical, the risk of dislocation and fracture of the ceramic increases. If, on the other hand, the angle of inclination is too small, flexion and abduction will be greatly limited. The study shows that acetabular components with non-recessed ceramic inserts should not be used. Slight recession of the insert helps to avoid impingement. The ROM is reduced and the risk of impingement appreciably increased when mushroom-shaped femoral heads (XL heads) or ceramic inserts protected by a polyethylene ring are used.     

Effect of cup inclination on predicted contact stress-induced volumetric wear in total hip replacement

In order to increase the lifetime of the total hip endoprosthesis, it is necessary to understand mechanisms leading to its failure. In this work, we address volumetric wear of the artificial cup, in particular the effect of its inclination with respect to the vertical. Volumetric wear was calculated by using mathematical models for resultant hip force, contact stress and penetration of the prosthesis head into the cup. Relevance of the dependence of volumetric wear on inclination of the cup (its abduction angle ?_A) was assessed by the results of 95 hips with implanted endoprosthesis. Geometrical parameters obtained from standard antero-posterior radiographs were taken as input data. Volumetric wear decreases with increasing cup abduction angle ?_A. The correlation within the population of 95 hips was statistically significant (P = 0.006). Large cup abduction angle minimises predicted volumetric wear but may increase the risk for dislocation of the artificial head from the cup in the one-legged stance. Cup abduction angle and direction of the resultant hip force may compensate each other to achieve optimal position of the cup with respect to wear and dislocation in the one-legged stance for a particular patient.     

Development of computational wear simulation of metal-on-metal hip resurfacing replacements

As one of the alternatives to traditional metal-on-polyethylene total hip replacements, metal-on-metal hip resurfacing prostheses demonstrating lower wear have been introduced for younger and more active patients during the past decade. However, in vitro hip simulator testing for the predicted increased lifetime of these surface replacements is time-consuming and costly. Computational wear modelling based on the Archard wear equation and finite element contact analysis was developed in this study for artificial hip joints and particularly applied to metal-on-metal resurfacing bearings under simulator testing conditions to address this issue. Wear factors associated with the Archard wear equation were experimentally determined and based on the short-term hip simulator wear results. The computational wear simulation was further extended to a long-term evaluation up to 50 million cycles assuming that the wear rate stays constant. The prediction from the computational model shows good agreement with the corresponding simulator study in terms of volumetric wear and the wear geometry. The simulation shows the progression of linear wear penetrations, and the complexity of contact stress distribution on the worn bearing surfaces. After 50 million cycles, the maximum linear wear was predicted to be approximately 6 and 8 microm for the cup and head, respectively, and no edge contact was found.     

Computer determination of contact stress distribution and size of weight bearing area in the human hip joint

The mathematical models and the corresponding computer program for determination of the hip joint contact force, the contact stress distribution, and the size of the weight bearing area from a standard anteroposterior radiograph are described. The described method can be applied in clinical practice to predict an optimal stress distribution after different operative interventions in the hip joint and to analyze the short and long term outcome of the treatment of various pathological conditions in the hip. A group of dysplastic hips and a group of normal hips were examined, with respect to the peak contact stress normalized by the body weight, and with respect to the functional angle of the weight bearing area. It is shown that both these parameters can be used in the assessment of hip dysplasia.     

The transmission of load through the human hip joint

This paper describes the results of loading experiments carried out on human hip joints. The unloaded surfaces of the femoral head and the acetabulum are slightly incongruous. The location and magnitude of the contact areas between the surfaces therefore depend on the magnitude and direction of the applied load. The contact areas were determined experimentally for a variety of loads typical of normal walking. Two distinct contact areas were found on the anterior and posterior aspects of the acetabulum at light loads, gradually merging with increasing load until, at a certain transition load, the dome of the acetabulum comes into contact and contact is then complete. The value of the transition load depends on the rate of loading, due to creep of the cartilage, and was found to vary from 50 per cent of body weight in young specimens to 25 per cent of body weight for elderly specimens for rates of loading typical of normal walking. Thus, the dome of the acetabulum is out of contact for a substantial portion of the swing phase of normal walking.

The analysis of a much simplified model of the hip joint is presented. The dependence of contact area on load is demonstrated, but also a method of determining the transition load for complete contact from the load/deflection relation for the hip is suggested. The values of the transition load quoted above were obtained by this method. The analysis further indicates that the distribution of pressure between the articular surfaces depends critically on the distribution of cartilage thickness throughout the joint. It is suggested that the distribution of cartilage thickness is such as to lead to a state of uniform pressure at the upper end of the physiological load range. Some experimental evidence is presented in support of this suggestion.

It is concluded that the function of joint incongruity is to allow the articular surfaces to come out of contact at light loads so that the cartilage may be exposed to synovial fluid for the purposes of nutrition and lubrication. At large loads, the distribution of cartilage thickness ensures that a state of hydrostatic pressure is achieved in order that cartilage, with a large fluid content, may transmit large pressures without flow and consequent loss of its integrity.     

Computer Determination of Contact Stress Distribution and Size of Weight Bearing Area in the Human Hip Joint

The mathematical models and the corresponding computer program for determination of the hip joint contact force, the contact stress distribution, and the size of the weight bearing area from a standard anteroposterior radiograph are described. The described method can be applied in clinical practice to predict an optimal stress distribution after different operative interventions in the hip joint and to analyze the short and long term outcome of the treatment of various pathological conditions in the hip. A group of dysplastic hips and a group of normal hips were examined, with respect to the peak contact stress normalized by the body weight, and with respect to the functional angle of the weight bearing area. It is shown that both these parameters can be used in the assessment of hip dysplasia.     

Effect of bearing geometry and structure support on transient elastohydrodynamic lubrication of metal-on-metal hip implants

An effective lubrication can significantly reduce wear of metal-on-metal artificial hip joints. The improvement of the lubrication can be achieved through the optimisation of the bearing geometry in terms of a small clearance and/or the structural support such as a polyethylene backing underneath a metallic bearing in a sandwich acetabular cup form. The separate effects of these two factors on fluid film lubrication of 28 mm diameter metal-on-metal total hip joints under walking conditions were numerically investigated in this paper. The results show that a larger lubricant film due to the polyethylene backing can be significantly enhanced by the transient squeeze-film action, particularly during the stance phase, and a similar lubricant film can be developed for both the monolithic cup relying on the smaller clearance and the sandwich cup benefiting from the polyethylene backing. Both cup systems can function in a wide range of lubrication regimes, covering both mixed and fluid film, under the current design and manufacture conditions.     

Effects of implant design parameters on fluid convection, potentiating third-body debris ingress into the bearing surface during THA impingement/subluxation

Aseptic loosening from polyethylene wear debris is the leading cause of failure for metal-on-polyethylene total hip implants. Third-body debris ingress to the bearing space results in femoral head roughening and acceleration of polyethylene wear. How third-body particles manage to enter the bearing space between the closely conforming articulating surfaces of the joint is not well understood. We hypothesize that one such mechanism is from convective fluid transport during subluxation of the total hip joint. To test this hypothesis, a three-dimensional (3D) computational fluid dynamics (CFD) model was developed and validated, to quantify fluid ingress into the bearing space during a leg-cross subluxation event. The results indicated that extra-articular joint fluid could be drawn nearly to the pole of the cup with even very small separations of the femoral head (<0.60mm). Debris suspended near the equator of the cup at the site of maximum fluid velocity just before the subluxation began could be transported to within 11 degrees from the cup pole. Larger head diameters resulted in increased fluid velocity at all sites around the entrance to the gap compared to smaller head sizes, with fluid velocity being greatest along the anterosuperolateral cup edge, for all head sizes. Fluid pathlines indicated that suspended debris would reach similar angular positions in the bearing space regardless of head size. Increased inset of the femoral head into the acetabular cup resulted both in higher fluid velocity and in transport of third-body debris further into the bearing space.     

The safe-zones for combined cup and neck anteversions that fulfill the essential range of motion and their optimum combination in total hip replacements

Reduction of the range of motion (ROM) until prosthetic impingement of a total hip replacement may lead to frequent impingement, subluxation and dislocation especially for patients with good hip movement. The ROM until prosthetic impingement can be calculated using the technical ROM (theta) and the cup and neck positions by a previously created mathematical formula. A larger (theta) with proper cup and neck positions results in a larger ROM. However there was only one paper written in English, which revealed the optimum theoretical combination of cup and neck anteversions. ROM of more than 110 degrees flexion, 30 degrees internal-rotation at 90 degrees flexion, 30 degrees extension and 40 degrees external-rotation were defined as the criteria for essential ROM for ADL. The safe-zones for combined cup anteversion (betaanat) and neck anteversion (b) were defined as the areas that fulfill all the criteria of ROM without prosthetic impingement. The safe-zones were created for 35 degrees , 45 degrees and 55 degrees cup abductions (alpha) and for 120 degrees and 135 degrees (theta). The safe-zones for combined (betaanat) and (b) were much larger for a 135 degrees (theta) than a 120 degrees (theta). Their safe-zones showed that (b) should be reduced if (betaanat) is increased and choosing a lower (alpha) requires that the sum of (betaanat) and (b) should be higher and vice versa. A (theta) of more than 135 degrees is recommended as it further increases the size of the safe-zone and provides a larger ROM, and the optimum values of combined cup and neck anteversions can be estimated by the formula: (alpha) + (betaanat) + 0.77(b) = 84.3.