A method of quantification of stress shielding in the proximal femur using hierarchical computational modeling

Stress shielding is a biomechanical phenomenon causing adaptive changes in bone strength and stiffness around metallic implants, which potentially lead to implant loosening. Accordingly, there is a need for standard, objective engineering measures of the “stress shielding” performances of an implant that can be employed in the process of computer-aided implant design. To provide and test such measures, we developed hierarchical computational models of adaptation of the trabecular microarchitecture at different sites in the proximal femur, in response to insertion of orthopaedic screws and in response to hypothetical reductions in hip joint and gluteal muscle forces. By identifying similar bone adaptation outcomes from the two scenarios, we were able to quantify the stress shielding caused by screws in terms of analogous hypothetical reductions in hip joint and gluteal muscle forces. Specifically, we developed planar lattice models of trabecular microstructures at five regions of interest (ROI) in the proximal femur. The homeostatic and abnormal loading conditions for the lattices were determined from a finite element model of the femur at the continuum scale and fed to an iterative algorithm simulating the adaptation of each lattice to these loads. When screws were inserted to the femur model, maximal simulated bone loss (17% decrease in apparent density, 10% decrease in thickness of trabeculae) was at the greater trochanter and this effect was equivalent to the effect of 50% reduction in gluteal force and normal hip joint force. We conclude that stress shielding performances can be quantified for different screw designs using model-predicted hypothetical musculoskeletal load fractions that would cause a similar pattern and extent of bone loss to that caused by the implants.     

Mechanical significance of femoral head trabecular bone structure in Loris and Galago evaluated using micromechanical finite element models

Work on the interspecific and intraspecific variation of trabecular bone in the proximal femur of primates demonstrates important architectural variation between animals with different locomotor behaviors. This variation is thought to be related to the processes of bone adaptation whereby bone structure is optimized to the mechanical environment. Micromechanical finite element models were created for the proximal femur of the leaping Galago senegalensis and the climbing and quadrupedal Loris tardigradus by converting bone voxels from high-resolution X-ray computed tomography scans of the femoral head to eight-noded brick elements. The resulting models had approximately 1.8 million elements each. Loading conditions representing takeoff phase of a leap and more generalized load orientations were applied to the models, and the models were solved using the iterative "row-by-row" matrix-vector multiplication algorithm. The principal strain and Von Mises stress results for the leaping model were similar for both species at each load orientation. Similar hip joint reaction forces in the range of 4.9 x to 12 x body weight were calculated for both species under each loading condition, but the hip reaction values estimated for Loris were higher than predicted based on locomotor behavior. These results suggest that functional adaptation to hip joint loading may not fully explain the differences in femoral head trabecular bone structure in Galago and Loris. The finite element method represents a unique and useful tool for analyzing the functional adaptation of trabecular bone in a diversity of animals and for reconstructing locomotor behavior in extinct taxa.     

A Femur-Implant Model for the Prediction of Bone Remodeling Behavior Induced by Cementless Stem

Bone remodeling simulation is an effective tool for the prediction of long-term effect of implant on the bone tissue, as well as the selection of an appropriate implant in terms of architecture and material. In this paper, a finite element model of proximal femur was developed to simulate the structures of internal trabecular and cortical bones by incorporating quantitative bone functional adaptation theory with finite element analysis. Cementless stems made of titanium, two types of Functionally Graded Material (FGM) and flexible ‘iso-elastic’ material as comparison were implanted in the structure of proximal femur respectively to simulate the bone remodeling behaviors of host bone. The distributions of bone density, von Mises stress, and interface shear stress were obtained. All the prosthetic stems had effects on the bone remodeling behaviors of proximal femur, but the degrees of stress shielding were different. The amount of bone loss caused by titanium implant was in agreement with the clinical observation. The FGM stems caused less bone loss than that of the titanium stem, in which FGM I stem (titanium richer at the top to more HAP/Col towards the bottom) could relieve stress shielding effectively, and the interface shear stresses were more evenly distributed in the model with FGM I stem in comparison with those in the models with FGM II (titanium and bioglass) and titanium stems. The numerical simulations in the present study provided theoretical basis for FGM as an appropriate material of femoral implant from a biomechanical point of view. The next steps are to fabricate FGM stem and to conduct animal experiments to investigate the effects of FGM stem on the remodeling behaviors using animal model.     

Influence of interface condition and implant design on bone remodelling and failure risk for the resurfaced femoral head

Resurfacing of the femur has experienced a revival, particularly in younger and more active patients. The implant is generally cemented onto the reamed trabecular bone and theoretical remodelling for this configuration, as well as uncemented variations, has been studied with relation to component positioning for the most common designs. The purpose of this study was to investigate the influence of different interface conditions, for alternative interior implant geometries, on bone strains in comparison to the native femur, and its consequent remodelling. A cylindrical interior geometry, two conical geometries and a spherical cortex-preserving design were compared with a standard implant (ASR, DePuy International, Ltd., UK), which has a 3° cone. Cemented as well as uncemented line to line and press-fit conditions were modelled for each geometry. A patient-specific finite element model of the proximal femur was used with simulated walking loads. Strain energy density was compared between the reference and resurfaced femur, and input into a remodelling algorithm to predict density changes post-operatively. The common cemented designs (cylindrical, slightly conical) had strain shielding in the superior femoral head (>35% reduction) as well as strain concentrations (strain>5%) in the neck regions near the implant rim. The cortex-preserving (spherical) and strongly conical designs showed less strain shielding. In contrast to the cemented implants, line to line implants showed a density decrease at the centre of the femoral head, while all press-fit versions showed a density increase (>100%) relative to the native femur, which suggests that uncemented press-fit implants could limit bone resorption.     

Epiphyseal-based designs for tibial plateau components--I. Stress analysis in the frontal plane

Two-dimensional, finite element studies were conducted of the proximal tibia before and after joint arthroplasty. Equivalent-thickness models projected onto the mid-frontal plane were created for the natural, proximal tibia and for the proximal tibia with four different types of tibial plateau components. All components simulated bony ingrowth fixation, i.e. no cement layer existed between component and bone. In addition, the interface between component and bone was assumed to be intimately connected, representing complete bony ingrowth and a rigid state of fixation. Loads consisted of bi-condylar and uni-condylar forces. Results indicated that conventional plateau designs with central posts or multiple pegs led to higher stress magnitudes in the trabecular bone near the distal ends of the post/pegs and stress shielding at more proximal locations. A design without posts or pegs whose interface geometry mimics the epiphyseal plate minimizes bone stress shielding. An implant consisting of separate components covering each condyle was found effective in limiting component tilting and the consequent tensile stresses caused by non-symmetrical, uni-condylar loading.     

Analysis of a femoral hip prosthesis designed to reduce stress shielding

The natural stress distribution in the femur is significantly altered after total hip arthroplasty (THA). When an implant is introduced, it will carry a portion of the load, causing a reduction of stress in some regions of the remaining bone. This phenomenon is commonly known as stress shielding. In response to the changed mechanical environment the shielded bone will remodel according to Wolff's law, resulting in a loss of bone mass through the biological process called resorption. Resorption can, in turn, cause or contribute to loosening of the prosthesis. The problem is particularly common among younger THA recipients. This study explores the hypothesis that through redesign, a total hip prosthesis can be developed to substantially reduce stress shielding. First, we describe the development of a new femoral hip prosthesis designed to alleviate this problem through a new geometry and system of proximal fixation. A numerical comparison with a conventional intramedullary prosthesis as well as another proximally fixed prosthesis, recently developed by Munting and Verhelpen (1995. Journal of Biomechanics 28(8), 949–961) is presented. The results show that the new design produces a more physiological stress state in the proximal femur.     

A finite element analysis of hollow stemmed hip prostheses as a means of reducing stress shielding of the femur.

Stress shielding of the femur is known to be a principal factor in aseptic loosening of hip replacements. This paper considers the use of a hollow stemmed hip implant for reducing the effects of stress shielding, while maintaining acceptably low levels of stress in the cement. Using finite element modelling, the stresses in the proximal femur using different shapes of hollow stem were compared with those produced using comparable sizes of solid stem with different values of elastic modulus. A reduction in stress shielding could be achieved with a hollow stem. A cylindrical hollow stem design was then optimised in order to control the maximum allowable stress in the cement, the minimum allowable stresses in the bone, and a combination of the two. The resulting stems achieved an increase in proximal bone stress of about 15% for the first case and 32% for a model using high strength cement, compared with solid stems of the same nominal outside diameter. The gains of these theoretically optimised designs dropped off rapidly further down the stem. Linearly tapered hollow stems reached a 22% gain, which could be a good compromise between acceptable cement stresses and ease of manufacture.     

Topological optimization in hip prosthesis design

With particular interest on total hip arthroplasty (THA), optimization of orthopedic prostheses is employed in this work to minimize the probability of implant failure or maximize prosthesis reliability. This goal is often identified with the reduction of stress concentrations at the interface between bone and these devices. However, aseptic loosening of the implant is mainly influenced by bone resorption phenomena revealed in some regions of the femur when a prosthesis is introduced. As a consequence, bone resorption appears due to stress shielding, that is to say the decrease of the stress level in the implanted femur caused by the significant load carrying of the prosthesis due to its higher stiffness. A maximum stiffness topological optimization-based (TO) strategy is utilized for non-linear static finite element (FE) analyses of the femur–implant assembly, with the goal of reducing stress shielding in the femur and to furnish guidelines for re-designing hip prostheses. This is accomplished by employing an extreme accuracy for both the three-dimensional reconstruction of the femur geometry and the material properties maps assigned as explicit functions of the local densities.     

The discrete nature of trabecular bone microarchitecture affects implant stability

Small endosseous implants, such as screws, are important components of modern orthopedics and dentistry. Hence they have to reliably fulfill a variety of requirements, which makes the development of such implants challenging. Finite element analysis is a widely used computational tool used to analyze and optimize implant stability in bone. For these purposes, bone is generally modeled as a continuum material. However, bone failure and bone adaptation processes are occurring at the discrete level of individual trabeculae; hence the assessment of stresses and strains at this level is relevant. Therefore, the aim of the present study was to investigate how peri-implant strain distribution and load transfer between implant and bone are affected by the continuum assumption. We performed a computational study in which cancellous screws were inserted in continuum and discrete models of trabecular bone; axial loading was simulated. We found strong differences in bone-implant stiffness between the discrete and continuum bone model. They depended on bone density and applied boundary conditions. Furthermore, load transfer from the screw to the surrounding bone differed strongly between the continuum and discrete models, especially for low-density bone. Based on our findings we conclude that continuum bone models are of limited use for finite element analysis of peri-implant mechanical loading in trabecular bone when a precise quantification of peri-implant stresses and strains is required. Therefore, for the assessment and improvement of trabecular bone implants, finite element models which accurately represent trabecular microarchitecture should be used.     

Computational simulation of trabecular adaptation progress in human proximal femur during growth

There are a large number of clinical and experimental studies that analyzed trabecular architecture as a result of bone adaptation. However, only a limited amount of quantitative data is currently available on the progress of trabecular adaptation during growth. In this paper, we proposed a two-step numerical simulation method that predicts trabecular adaptation progress during growth using a recently developed topology optimization algorithm, design space optimization (DSO), under the hypothesis that the mechanisms of DSO are functionally equivalent to those of bone adaptation. We applied the proposed scheme to trabecular adaptation simulation in human proximal femur. For the simulation, the full trabecular architecture in human proximal femur was represented by a two-dimensional μFE model with 50 μm resolution. In Step 1, we determined a reference value that regulates trabecular adaptation in human proximal femur. In Step 2, we simulated trabecular adaptation in human proximal femur during growth with the reference value derived in Step 1. We analyzed the architectural and mechanical properties of trabecular patterns through iterations. From the comparison with experimental data in the literature, we showed that in the early growth stage trabecular adaptation was achieved mainly by increasing bone volume fraction (or trabecular thickness), while in the later stage of the development the trabecular architecture gained higher structural efficiency by increasing structural anisotropy with a relatively low level of bone volume fraction (or trabecular thickness). We demonstrated that the proposed numerical framework predicted the growing progress of trabecular bone that has a close correlation with experimental data.     

Relation between subject-specific hip joint loading, stress distribution in the proximal femur and bone mineral density changes after total hip replacement

In the prediction of bone remodelling processes after total hip replacement (THR), modelling of the subject-specific geometry is now state-of-the-art. In this study, we demonstrate that inclusion of subject-specific loading conditions drastically influences the calculated stress distribution, and hence influences the correlation between calculated stress distributions and changes in bone mineral density (BMD) after THR.For two patients who received cementless THR, personalized finite element (FE) models of the proximal femur were generated representing the pre- and post-operative geometry. FE analyses were performed by imposing subject-specific three-dimensional hip joint contact forces as well as muscle forces calculated based on gait analysis data. Average values of the von Mises stress were calculated for relevant zones of the proximal femur. Subsequently, the load cases were interchanged and the effect on the stress distribution was evaluated. Finally, the subject-specific stress distribution was correlated to the changes in BMD at 3 and 6 months after THR.We found subject-specific differences in the stress distribution induced by specific loading conditions, as interchanging of the loading also interchanged the patterns of the stress distribution. The correlation between the calculated stress distribution and the changes in BMD were affected by the two-dimensional nature of the BMD measurement.Our results confirm the hypothesis that inclusion of subject-specific hip contact forces and muscle forces drastically influences the stress distribution in the proximal femur. In addition to patient-specific geometry, inclusion of patient-specific loading is, therefore, essential to obtain accurate input for the analysis of stress distribution after THR.     

Hip contact forces and gait patterns from routine activities.

In vivo loads acting at the hip joint have so far only been measured in few patients and without detailed documentation of gait data. Such information is required to test and improve wear, strength and fixation stability of hip implants. Measurements of hip contact forces with instrumented implants and synchronous analyses of gait patterns and ground reaction forces were performed in four patients during the most frequent activities of daily living. From the individual data sets an average was calculated. The paper focuses on the loading of the femoral implant component but complete data are additionally stored on an associated compact disc. It contains complete gait and hip contact force data as well as calculated muscle activities during walking and stair climbing and the frequencies of daily activities observed in hip patients. The mechanical loading and function of the hip joint and proximal femur is thereby completely documented. The average patient loaded his hip joint with 238% BW (percent of body weight) when walking at about 4 km/h and with slightly less when standing on one leg. This is below the levels previously reported for two other patients (Bergmann et al., Clinical Biomechanics 26 (1993) 969-990). When climbing upstairs the joint contact force is 251% BW which is less than 260% BW when going downstairs. Inwards torsion of the implant is probably critical for the stem fixation. On average it is 23% larger when going upstairs than during normal level walking. The inter- and intra-individual variations during stair climbing are large and the highest torque values are 83% larger than during normal walking. Because the hip joint loading during all other common activities of most hip patients are comparably small (except during stumbling), implants should mainly be tested with loading conditions that mimic walking and stair climbing.     

Modulation of bone ingrowth of rabbit femur titanium implants by in vivo axial micromechanical loading.

Titanium implants commonly used in orthopedics and dentistry integrate into host bone by a complex and coordinated process. Despite increasingly well illustrated molecular healing processes, mechanical modulation of implant bone ingrowth is poorly understood. The objective of the present study was to determine whether micromechanical forces applied axially to titanium implants modulate bone ingrowth surrounding intraosseous titanium implants. We hypothesized that small doses of micromechanical forces delivered daily to the bone-implant interface enhance implant bone ingrowth. Small titanium implants were placed transcortically in the lateral aspect of the proximal femur in 15 New Zealand White rabbits under general anesthesia and allowed to integrate with the surrounding bone for 6 wk. Micromechanical forces at 200 mN and 1 Hz were delivered axially to the right femur implants for 10 min/day over 12 consecutive days, whereas the left femur implants served as controls. The average bone volume 1 mm from mechanically loaded implants (n = 15) was 73 +/- 12%, which was significantly greater than the average bone volume (52 +/- 21%) of the contralateral controls (n = 15) (P < 0.01). The average number of osteoblast-like cells per endocortical bone surface was 55 +/- 8 cells/mm(2) for mechanically loaded implants, which was significantly greater than the contralateral controls (35 +/- 6 cells/mm(2)) (P < 0.01). Dynamic histomorphometry showed a significant increase in mineral apposition rate and bone-formation rate of mechanically stressed implants (3.8 +/- 1.2 microm/day and 2.4 +/- 1.0 microm(3).microm(-2).day(-1), respectively) than contralateral controls (2.2 +/- 0.92 microm/day and 1.2 +/- 0.60 microm(3).microm(-2).day(-1), respectively; P < 0.01). Collectively, these data suggest that micromechanical forces delivered axially on intraosseous titanium implants may have anabolic effects on implant bone ingrowth.     

Functional adaptation of cancellous bone in human proximal femur predicted by trabecular surface remodeling simulation toward uniform stress state

Two-dimensional simulation of trabecular surface remodeling was conducted for a human proximal femur to investigate the structural change of cancellous bone toward a uniform stress state. Considering that a local mechanical stimulus plays an important role in cellular activities in bone remodeling, local stress nonuniformity was assumed to drive trabecular structural change to seek a uniform stress state. A large-scale pixel-based finite element model was used to simulate structural changes of individual trabeculae over the entire bone. As a result, the initial structure of trabeculae changed from isotropic to anisotropic due to trabecular microstructural changes caused by surface remodeling according to the mechanical environment in the proximal femur. Under a single-loading condition, it was shown that the apparent structural property evaluated by fabric ellipses corresponded to the apparent stress state in cancellous bone. As is observed in the actual bone, a distributed trabecular structure was obtained under a multiple-loading condition. Through these studies, it was concluded that trabecular surface remodeling toward a local uniform stress state at the trabecular level could naturally bring about functional adaptation phenomenon at the apparent tissue level. The proposed simulation model would be capable of providing insight into the hierarchical mechanism of trabecular surface remodeling at the microstructural level up to the apparent tissue level.     

Computational simulation of simultaneous cortical and trabecular bone change in human proximal femur during bone remodeling

In this study, we developed a numerical framework that computationally determines simultaneous and interactive structural changes of cortical and trabecular bone types during bone remodeling, and we investigated the structural correlation between the two bone types in human proximal femur. We implemented a surface remodeling technique that performs bone remodeling in the exterior layer of the cortical bone while keeping its interior area unchanged. A micro-finite element (μFE) model was constructed that represents the entire cortical bone and full trabecular architecture in human proximal femur. This study simulated and compared the bone adaptation processes of two different structures: (1) femoral bone that has normal cortical bone shape and (2) perturbed femoral bone that has an artificial bone lump in the inferomedial cortex. Using the proposed numerical method in conjunction with design space optimization, we successfully obtained numerical results that resemble actual human proximal femur. The results revealed that actual cortical bone, as well as the trabecular bone, in human proximal femur has structurally optimal shapes, and it was also shown that a bone abnormality that has little contribution to bone structural integrity tends to disappear. This study also quantitatively determined the structural contribution of each bone: when the trabecular adaptation was complete, the trabecular bone supported 54% of the total load in the human proximal femur while the cortical bone carried 46%.     

Densitometric, morphometric and mechanical distributions in the human proximal femur

With the prevalent use of DXA-measured BMD to assess pathologic hip fractures and its recently reported lack of reliability to predict fracture or account for efficacy of anti-resorptive therapy, it is reasonable to assess whether variations in the primary and secondary tensile and compressive trabecular microstructure can account for variations in proximal femur strength in comparison to DXA-measured BMD. To that end, microstructural and densitometric measures of trabecular bone specimens, from discrete sites within the proximal femur, were correlated with their mechanical properties. We hypothesize that accounting for regional variations in trabecular microstructure will improve predictions of proximal femur strength and stiffness compared to bone density measured by DXA. Forty-seven samples (seven donors) from seven distinct sites of human proximal femur underwent DXA and muCT imaging and mechanical testing. The results revealed significant variations in BMC, morphometric indices and mechanical properties within the proximal femur. This work has demonstrated that the mechanical performance of each sub-region is highly dependent on the corresponding trabecular microstructure. BMD measured by DXA at standard regions of interest cannot resolve the variations in trabecular density and microstructure that govern the mechanical behavior of the proximal femur. This work suggests that a quantitative Singh index that uses high resolution QCT to monitor the trabecular microstructure at specific sub-regions of the proximal femur may allow better predictions of hip fracture risk in individual patients and an improved assessment of changing bone structure in response to pharmacological interventions.     

Trabecular bone adaptation with an orthotropic material model.

Most bone adaptation algorithms, that attempt to explain the connection between bone morphology and loads, assume that bone is effectively isotropic. An isotropic material model can explain the bone density distribution, but not the structure and pattern of trabecular bone, which clearly has a mechanical significance. In this paper, an orthotropic material model is utilized to predict the proximal femur trabecular structure. Two hypotheses are combined to determine the local orientation and material properties of each element in the model. First, it is suggested that trabecular directions, which correspond to the orthotropic material axes, are determined locally by the maximal principal stress directions due to the multiple load cases (MLC) the femur is subject to. The second hypothesis is that material properties in each material direction can be determined using directional stimuli, thus extending existing adaptation algorithms to include directionality. An algorithm is utilized, where each iteration comprises of two stages. First, material axes are rotated to the direction of the largest principal stress that occurs from a multiple load scheme applied to the proximal femur. Next, material properties are modified in each material direction, according to a directional stimulus. Results show that local material directions correspond with known trabecular patterns, reproducing all main groups of trabeculae very well. The local directional stiffnesses, degree of anisotropy and density distribution are shown to conform to real femur morphology.     

On the optimal shape of hip implants

The success of a total hip arthroplasty is strongly related to the initial stability of the femoral component and to the stress shielding effect. In fact, for cementless stems, initial stability is essential to promote bone ingrowth into the stem coating. An inefficient primary stability is also a cause of thigh pain. In addition, the bone adaptation after the surgery can lead to an excessive bone loss and, consequently, can compromise the success of the implant. These factors depend on prosthesis design, namely on material, interface conditions and shape. Although, surgeons use stems with very different geometries, new computational tools using structural optimization methods have been used to achieve a better design in order to improve initial stability and therefore, the implant durability. In this work, a multi-criteria shape optimization process is developed to study the relationship between implants performance and geometry. The multi-criteria objective function takes into account the initial stability of the femoral stem and the effect of stress shielding on bone adaptation after the surgery. Then, the optimized stems are tested using a concurrent model for bone remodeling and osseointegration to evaluate long-term performance. Additionally, the sensitivity to misalignments is analyzed, since femoral stems are often placed in varus or valgus position. Results show that the different criteria are contradictory resulting in different characteristics for the hip stem. However, the multi-criteria formulation leads to compromise solutions, with a combination of the geometric characteristics obtained for each criterion separately.     

Implant–bone interface healing and adaptation in resurfacing hip replacement

Hip resurfacing demonstrates good survivorship as a treatment for young patients with osteoarthritis, but occasional implant loosening failures occur. On the femoral side there is radiographic evidence suggesting that the implant stem bears load, which is thought to lead to proximal stress shielding and adaptive bone remodelling. Previous attempts aimed at reproducing clinically observed bone adaptations in response to the implant have not recreated the full set of common radiographic changes, so a modified bone adaptation algorithm was developed in an attempt to replicate more closely the effects of the prosthesis on the host bone. The algorithm features combined implant–bone interface healing and continuum bone remodelling. It was observed that remodelling simulations that accounted for progressive gap filling at the implant–bone interface predicted the closest periprosthetic bone density changes to clinical X-rays and DEXA data. This model may contribute to improved understanding of clinical failure mechanisms with traditional hip resurfacing designs and enable more detailed pre-clinical analysis of new designs.