The biomechanical effect of the collar of a femoral stem on total hip arthroplasty

To investigate the biomechanical effect of collars, finite element analyses are carried out through two hip joints that are implanted using collared and collarless stems, respectively, and an intact hip joint model. For the analyses, the sacrum, coxal bone, and the cancellous and cortical bones of a femur are modelled using finite elements based on X-ray computed tomographic images taken from a 27-year-old woman. From the results, it is found that a collar with perfect calcar contact prevents stem subsidence and decreases the proximal–lateral gap and the lateral stem tilting. Therefore, it can impart reasonable biomechanical stability for total hip arthroplasty. However, its low load transmission ability and increased stem tilting effect due to the imperfect contact between the collar and the calcar are found to be serious problems that need to be solved. Results of clinical follow-up are presented for supporting the computational results.     

Stress analysis of the proximo-medial femur after total hip replacement

Stress-induced bone loss in the proximo-medial femur has been identified as a factor leading to loosening in the artificial hip joint. In an effort to develop a quantitative understanding of the stress distribution that causes bone loss, axial and hoop stresses in the medial calcar of the femur have been determined after total hip replacement, using finite element stress analysis. Stress distributions for a high and a low Young's modulus prosthesis material are compared for both collared and uncollared prosthesis designs. The use of a low-modulus material, and of a collar, are predicted to be advantageous, giving rise to proximo-medial stress patterns similar to those of the normal, intact femur.     

Design optimization of a prosthesis stem reinforcing shell in total hip arthroplasty

The use of a perforated, titanium funicular shell to support the proximal femoral cortex in total hip arthroplasty was evaluated with the aid of both analytical and numerical techniques. The principal interactions between the femoral cortex, the metal shell, the implant stem and the acrylic bone cement were modeled using beam on elastic foundations theory and two-dimensional elasticity theory. Subsequent formulation of this model as a nonlinear design optimization problem enabled the determination of the dimensions of the implant and reinforcing shell which minimized an objective function based on a simplified material failure criterion. Two cases were examined, each with two cervico-diaphyseal angles: case A: with a rigid contact between a proximal prosthesis collar and the calcar femorale and case B: no collar contact (a collarless prosthesis or post-operative loosening). Case A achieved an optimal solution at a stem diameter 11-23 percent of the cortex inner diameter, a stem length to diameter ratio of 12-40, shell diameter 22-53 percent and thickness 0.2-7.2 percent of the cortex inner diameter and thickness, respectively. Case B achieved an optimal solution at a stem diameter 67-92 percent of the cortex inner diameter, length to diameter ratio of 4-6, and no shell. In case A the collar support makes the type of internal fixation unimportant, while in the more realistic case B, the shell is not recommended.     

Numerical simulation of the insertion process of an uncemented hip prosthesis in order to evaluate the influence of residual stress and contact distribution on the stem initial stability

The long-term success of a cementless total hip arthroplasty depends on the implant geometry and interface bonding characteristics (fit, coating and ingrowth) and on stem stiffness. This study evaluates the influence of stem geometry and fitting conditions on the evolution and distribution of the bone–stem contact, stress and strain during and after the hip stem insertion, by means of dynamic finite element techniques. Next, the influence of the mechanical state (bone–stem contact, stress and strain) resulted from the insertion process on the stem initial resistance to subsidence is investigated. In addition, a study on the influence of bone–stem interface conditions (friction) on the insertion process and on the initial stem stability under physiological loading is performed. The results indicate that for a stem with tapered shape the contact in the proximal part of the stem was improved, but contact in the calcar region was achieved only when extra press-fit conditions were considered. Changes in stem geometry towards a more tapered shape and extra press fit and variation in the bone–stem interface conditions (contact amount and high friction) led to a raise in the total insertion force. A direct positive relationship was found between the stem resistance to subsidence and stem geometry (tapering and press fit), bone–stem interface conditions (bone–stem contact and friction interface) and the mechanical status at the end of the insertion (residual stress and strain). Therefore, further studies on evaluating the initial performance of different stem types should consider the parameters describing the bone–stem interface conditions and the mechanical state resulted from the insertion process.     

Hip biomechanics in orthopaedic clinical practice

Radiographic and clinical studies, coupled with biomechanical assessment of the hip, are important tools for predicting the development of osteoarthitis of the hip. In order to better understand the treatment of hip dysplasia, it is necessary to determine the contact stress in the hip joint. In this study, a three-dimensional mathematical model was used to determine hip joint contact stress. Because of the discrepancy in the results of analyses of different radiographic indicators of hip dysplasia, the calculation of hip joint contact stress is proposed for a more accurate assessment of the severity of hip dysplasia.     

Shape optimization of a cementless hip stem for a minimum of interface stress and displacement

The primary stem stability is an essential factor for success of cementless hip stems. A correct choice of the stem geometry can improve the stem stability and, consequently, increase the life time of a hip implant. In this work, it is proposed a computational model for shape optimization of cementless hip stems. The optimization problem is formulated by the minimization of relative displacement and stress on bone/stem interface using a multi-criteria objective function. Also multiple loads are considered to incorporate several daily life activities. Design variables are parameters that characterize the geometry of selected cross sections, which are subject to geometric constraints to ensure a clinically admissible shape. The stem/bone set is considered a structure in equilibrium with contact conditions on interface. The contact formulation allows us to analyze different lengths of porous coating. The optimization problem is solved numerically by a steepest descent method. The interface stress and relative displacement are obtained solving the contact problem by the finite element method. Numerical examples are presented for a two-dimensional model of a hip stem, however, the formulation is general and can be applied to the three-dimensional case. The model gives indications about the relation between shape, porous coating and prosthesis stability.     

Shape Optimization of a Cementless Hip Stem for a Minimum of Interface Stress and Displacement

The primary stem stability is an essential factor for success of cementless hip stems. A correct choice of the stem geometry can improve the stem stability and, consequently, increase the life time of a hip implant. In this work, it is proposed a computational model for shape optimization of cementless hip stems. The optimization problem is formulated by the minimization of relative displacement and stress on bone/stem interface using a multi-criteria objective function. Also multiple loads are considered to incorporate several daily life activities. Design variables are parameters that characterize the geometry of selected cross sections, which are subject to geometric constraints to ensure a clinically admissible shape. The stem/bone set is considered a structure in equilibrium with contact conditions on interface. The contact formulation allows us to analyze different lengths of porous coating. The optimization problem is solved numerically by a steepest descent method. The interface stress and relative displacement are obtained solving the contact problem by the finite element method. Numerical examples are presented for a two-dimensional model of a hip stem, however, the formulation is general and can be applied to the three-dimensional case. The model gives indications about the relation between shape, porous coating and prosthesis stability.     

A Three-Dimensional Finite Element Model for Biomechanical Analysis of the Hip

The objective of this study was to construct a three-dimensional (3D) finite element model of the hip. The images of the hip were obtained from Chinese visible human dataset. The hip model includes acetabular bone, cartilage, labrum, and bone. The cartilage of femoral head was constructed using the AutoCAD and Solidworks software. The hip model was imported into ABAQUS analysis system. The contact surface of the hip joint was meshed. To verify the model, the single leg peak force was loaded, and contact area of the cartilage and labrum of the hip and pressure distribution in these structures were observed. The constructed 3D hip model reflected the real hip anatomy. Further, this model reflected biomechanical behavior similar to previous studies. In conclusion, this 3D finite element hip model avoids the disadvantages of other construction methods, such as imprecision of cartilage construction and the absence of labrum. Further, it provides basic data critical for accurately modeling normal and abnormal loads, and the effects of abnormal loads on the hip.     

Subject-specific hip geometry and hip joint centre location affects calculated contact forces at the hip during gait

Hip loading affects the development of hip osteoarthritis, bone remodelling and osseointegration of implants. In this study, we analyzed the effect of subject-specific modelling of hip geometry and hip joint centre (HJC) location on the quantification of hip joint moments, muscle moments and hip contact forces during gait, using musculoskeletal modelling, inverse dynamic analysis and static optimization. For 10 subjects, hip joint moments, muscle moments and hip loading in terms of magnitude and orientation were quantified using three different model types, each including a different amount of subject-specific detail: (1) a generic scaled musculoskeletal model, (2) a generic scaled musculoskeletal model with subject-specific hip geometry (femoral anteversion, neck-length and neck-shaft angle) and (3) a generic scaled musculoskeletal model with subject-specific hip geometry including HJC location. Subject-specific geometry and HJC location were derived from CT. Significant differences were found between the three model types in HJC location, hip flexion–extension moment and inclination angle of the total contact force in the frontal plane. No model agreement was found between the three model types for the calculation of contact forces in terms of magnitude and orientations, and muscle moments. Therefore, we suggest that personalized models with individualized hip joint geometry and HJC location should be used for the quantification of hip loading. For biomechanical analyses aiming to understand modified hip joint loading, and planning hip surgery in patients with osteoarthritis, the amount of subject-specific detail, related to bone geometry and joint centre location in the musculoskeletal models used, needs to be considered.     

Normal hip joint contact pressure distribution in single-leg standing--effect of gender and anatomic parameters.

A practical and easy-to-use analysis technique that can study the patient's hip joint contact force/pressure distribution would be useful to assess the effect of abnormal biomechanical conditions and anatomical deformities on joint contact stress for treatment planning purpose. This technique can also help to establish the normative database on hip joint contact pressure distribution in men and women in different age groups. Twelve anatomic parameters and seven biomechanical parameters of the hip joint in a normal population (41 females, 15 males) were calculated. The inter-parameter correlations were investigated. The pressure distribution in the hip joint was calculated using a three-dimensional discrete element analysis (DEA) technique. The 3D contact geometry of the hip joint was estimated from a 2D radiograph by assuming that the femoral head and the acetabular surface were spherical in shape. The head-trochanter ratio (HT), femoral head radius, pelvic height, the joint contact area, the normalized peak contact pressure, abductor force, and the joint contact force were significantly different between men and women. The normalized peak contact pressure was correlated both with acetabular coverage and head-trochanter ratio. Change of abductor force direction within normal variation did not affect the joint peak contact pressure. However, in simulated dysplastic conditions when the CE angle is small or negative, abductor muscle direction becomes very sensitive in joint contact pressure estimation. The models and the results presented can be used as the reference base in computer simulation for preoperative planning in pelvic or femoral osteotomy.     

Probabilistic finite element analysis of the uncemented hip replacement—effect of femur characteristics and implant design geometry

In the present study, a probabilistic finite element tool was assessed using an uncemented total hip replacement model. Fully bonded and frictional interfaces were investigated for combinations of three proximal femurs and two implant designs, the Proxima short stem and the IPS hip stem prostheses. The Monte Carlo method was used with two performance indicators: the percentage of bone volume that exceeded specified strain limits and the maximum nodal micromotion. The six degrees of freedom of bone-implant relative position, magnitude of the hip contact force (L), and spatial direction of L were the random variables. The distal portion of the proximal femurs was completely constrained and some of the main muscle forces acting in the hip were applied. The coefficients of the linear approximation between the random variables and the output were used as the sensitivity values. In all cases, bone-implant position related parameters were the most sensitive parameters. The results varied depending on the femur, the implant design and the interface conditions. Values of maximum nodal micromotion agreed with results from previous studies, confirming the robustness of the implemented computational tool. It was demonstrated that results from a single model study should not be generalised to the entire population of femurs and that bone variability is an important factor that should be investigated in such analyses.     

Hip contact stress during normal and staircase walking: the influence of acetabular anteversion angle and lateral coverage of the acetabulum

Hip contact stress is considered to be an important biomechanical factor related to development of coxarthrosis. The effect of the lateral coverage of the acetabulum on the hip contact stress has been demonstrated in several studies of hip dysplasia, whereas the effect of the anterior anteversion remains unclear. Therefore, the joint hip contact stress during normal level walking and staircase walking, in normal and dysplastic hips, for small and large acetabular anteversion angle was computed. For small acetabular anteversion angle, the hip contact stress is slightly increased (less than 15%) in staircase walking when compared with normal walking. In hips with large angle of acetabular anteversion, walking downstairs significantly increases the maximal peak contact stress (70% in normal hips and 115% in dysplastic hips) whereas walking upstairs decreases the peak contact stress (4% in normal hips and 34% in dysplastic hips) in comparison to normal walking. Based on the presented results, we suggest that the acetabular anteversion should be considered in biomechanical evaluation of the hips, especially when the lateral coverage of the acetabulum is small.     

Improving stress shielding following total hip arthroplasty by using a femoral stem made of β type Ti-33.6Nb-4Sn with a Young’s modulus gradation

Stress shielding-related bone loss occurs after total hip arthroplasty because the stiffness of metallic implants differs from that of the host femur. Although reducing stem stiffness can ameliorate the bone resorption, it increases stress at the bone–implant interface and can inhibit fixation. To overcome this complication, a novel cementless stem with a gradient in Young’s modulus was developed using Ti-33.6Nb-4Sn (TNS) alloy. Local heat treatment applied at the neck region for increasing its strength resulted in a gradual decrease in Young’s modulus from the proximal to the distal end, from 82.1 to 51.0 GPa as calculated by a heat transfer simulation. The Young’s modulus gradient did not induce the excessive interface stress which may cause the surface debonding. The main purpose of this study was to evaluate bone remodeling with the TNS stem using a strain-adaptive bone remodeling simulation based on finite element analysis. Our predictions showed that, for the TNS stem, bone reduction in the calcar region (Gruen zone 7) would be 13.6% at 2 years, 29.0% at 5 years, and 45.8% at 10 years postoperatively. At 10 years, the bone mineral density for the TNS stem would be 42.6% higher than that for the similar Ti-6Al-4V alloy stem. The stress–strength ratio would be lower for the TNS stem than for the Ti-6Al-4V stem. These results suggest that although proximal bone loss cannot be eliminated completely, the TNS stem with a Young’s modulus gradient may have bone-preserving effects and sufficient stem strength, without the excessive interface stress.     

Hip-joint and abductor-muscle forces adequately represent in vivo loading of a cemented total hip reconstruction.

Using finite element analyses, we investigated which muscle groups acting around the hip-joint most prominently affected the load distributions in cemented total hip reconstructions with a bonded and debonded femoral stem. The purpose was to determine which muscle groups should be included in pre-clinical tests, predicting bone adaptation and mechanical failure of cemented reconstructions, ensuring an adequate representation of in vivo loading of the reconstruction. Loads were applied as occurring during heel-strike, mid-stance and push-off phases of gait. The stress/strain distributions within the reconstruction, produced by the hip-joint contact force, were compared to ones produced after sequentially including the abductors, the iliotibial tract and the adductors and vastii. Inclusion of the abductors had the most pronounced effect. They neutralized lateral bending of the reconstruction at heel-strike and increased medial bending at mid-stance and push-off. Bone strains and stem stresses were changed accordingly. Peak tensile cement stresses were reduced during all gait phases by amounts up to 50% around a bonded stem and 11% around a debonded one. Additional inclusion of the iliotibial tract, the adductors and the vastii produced relatively small effects during all gait phases. Their most prominent effect was a slight reduction of bone strains at the level of the stem tip during heel-strike. These results suggest that a loading configuration including the hip-joint contact force and the abductor forces can adequately reproduce in vivo loading of cemented total hip reconstructions in pre-clinical tests.     

Primary stability of cementless stem in THA improved with reduced interfacial gaps

Large interfacial gaps between the stem and the bone in cementless total hip arthroplasty may prevent successful bone ingrowth at the sites, and can also be a passage for wear particles. Furthermore, interfacial gaps between the stem and the bone are believed to compromise the primary stability of the implant. Thus, a broaching method that serves to reduce gaps is expected to give clinically preferable results. A modified broach system with a canal guide is introduced to enhance the accuracy of femoral canal shaping in comparison with the conventional broach system for a Versys fibermetal taper stem. The primary stability of the hip systems and the ratios of the stem surface in contact with the femur were measured in a composite femur model. With the conventional method, an average of 67% of the stem surface was shown to be in contact with the bone, and an average stem micromotion/migration of 35 microm 290 microm was observed under 1000 cycles of stair climbing loads. With the modified method, the stem-bone contact ratio significantly increased to 82% (p<0.05), and the average micromotion/migration reduced to 29 microm 49 microm, respectively (p<0.05 for migration). Our finite element models of the hip systems supported that the difference in micromotion could be attributed to the difference in interfacial contact. Interfacial gaps occurring with the conventional broach system were effectively reduced by the proposed method, resulting in improved primary stability.     

Secondary hip dysplasia increases risk for early coxarthritis after Legg-Calve-Perthes disease. A study of 255 hips

Abstract

The biomechanical parameters of the hip joint articular surface were analysed in 141 adult hips after Legg-Calve Perthes Disease, and 114 contralateral unaffected hips (controls), by using HIPSTRESS mathematical models. Geometrical parameters, assessed from anteroposterior and axial radiograms, were used as input to models for resultant hip force and contact hip stress. Results confirm previous indications that head enlargement after the Legg-Calve-Perthes Disease compensates the values of hip stress. Furthermore, it was found that an increased risk for coxarthritis development after the disease is secondary to concomitant hip dysplasia, with considerable and statistically significantly lower centre-edge angle and unfavourable distribution of stress.     

Fatigue failure of the sliding hip screw – clinical and biomechanical analysis

The study was aimed at the clinical and biomechanical analyses of the sliding hip screw breakage with the use of finite element method. We have identified two patients with the sliding screw breakage. In the first patient, the biomechanical analysis revealed the reduced stress values σ_HMH not exceeding the yield limit or strength limit of the implant. The yield limit was exceeded in second one. Clinical and biomechanical analyses have demonstrated that adherence to technical requirements of the appropriate osteosynthesis implementation is the principal condition of timely healing since it prevents the material failure.     

Finite element algorithm for frictionless contact of porous permeable media under finite deformation and sliding

This study formulates and implements a finite element contact algorithm for solid-fluid (biphasic) mixtures, accommodating both finite deformation and sliding. The finite element source code is made available to the general public. The algorithm uses a penalty method regularized with an augmented Lagrangian method to enforce the continuity of contact traction and normal component of fluid flux across the contact interface. The formulation addresses the need to automatically enforce free-draining conditions outside of the contact interface. The accuracy of the implementation is verified using contact problems, for which exact solutions are obtained by alternative analyses. Illustrations are also provided that demonstrate large deformations and sliding under configurations relevant to biomechanical applications such as articular contact. This study addresses an important computational need in the biomechanics of porous-permeable soft tissues. Placing the source code in the public domain provides a useful resource to the biomechanics community.     

Long-term study of bone remodelling after femoral stem: a comparison between dexa and finite element simulation

A hip replacement with a cemented or cementless femoral stem produces an effect on the bone called adaptive remodelling, attributable to mechanical and biological factors. The objective of all of cementless prostheses designs has been to achieve a perfect transfer of loads in order to avoid stress-shielding, which produces an osteopenia. In order to quantify this, the long term and mass-produced study with dual energy X-ray absorptiometry (DEXA) is necessary. Finite element (FE) simulation makes possible the explanation of the biomechanical changes which are produced in the femur after stem implantation. The good correlation obtained between the results of the FE simulation and the densitometric study allow, on one hand, to explain from the point of view of biomechanical performance the changes observed in bone density in the long-term, where it is clear that these are due to a different transfer of load in the implanted model compared to the healthy femur; on the other hand, it validates the simulation model, in a way that it can be used in different conditions and at different time periods, to carry out a sufficiently precise prediction of the evolution of the bone density from the biomechanical behaviour in the interaction between the prosthesis and femur.     

New design of hip prosthesis using carbon fibre reinforced composite

We present a new design for a hip prosthesis using polyethylene-hydroxyapatite, a composite material with layered stiffness and good bone compatibility. In order to overcome the low stiffness of the matrix, carbon fibre reinforcement is suggested. This new collarless prosthesis was designed to give maximum stiffness at the core and lower stiffness at the bone-implant interface. For this design a 2D finite element analysis was undertaken; its results were compared with those for a titanium alloy prosthesis. The effect of the variation in the collar stiffness was also analysed.