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.     

Three-dimensional micro-level computational study of Wolff's law via trabecular bone remodeling in the human proximal femur using design space topology optimization

The law of bone remodeling, commonly referred to as Wolff's Law, asserts that the internal trabecular bone adapts to external loadings, reorienting with the principal stress trajectories to maximize mechanical efficiency creating a naturally optimum structure. The goal of the current study was to utilize an advanced structural optimization algorithm, called design space optimization (DSO), to perform a micro-level three-dimensional finite element bone remodeling simulation on the human proximal femur and analyse the results to determine the validity of Wolff's hypothesis. DSO optimizes the layout of material by iteratively distributing it into the areas of highest loading, while simultaneously changing the design domain to increase computational efficiency. The result is a "fully stressed" structure with minimized compliance and increased stiffness. The large-scale computational simulation utilized a 175 μm mesh resolution and the routine daily loading activities of walking and stair climbing. The resulting anisotropic trabecular architecture was compared to both Wolff's trajectory hypothesis and natural femur samples from literature using a variety of visualization techniques, including radiography and computed tomography (CT). The results qualitatively revealed several anisotropic trabecular regions, that were comparable to the natural human femurs. Quantitatively, the various regional bone volume fractions from the computational results were consistent with quantitative CT analyses. The global strain energy proceeded to become more uniform during optimization; implying increased mechanical efficiency was achieved. The realistic simulated trabecular geometry suggests that the DSO method can accurately predict bone adaptation due to mechanical loading and that the proximal femur is an optimum structure as the Wolff hypothesized.     

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%.     

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.     

Multiscale modeling of bone tissue with surface and permeability control

Natural biological materials usually present a hierarchical arrangement with various structural levels. The biomechanical behavior of the complex hierarchical structure of bone is investigated with models that address the various levels corresponding to different scales. Models that simulate the bone remodeling process concurrently at different scales are in development. We present a multiscale model for bone tissue adaptation that considers the two top levels, whole bone and trabecular architecture. The bone density distribution is calculated at the macroscale (whole bone) level, and the trabecular structure at the microscale level takes into account its mechanical properties as well as surface density and permeability. The bone remodeling process is thus formulated as a material distribution problem at both scales. At the local level, the biologically driven information of surface density and permeability characterizes the trabecular structure. The model is tested by a three-dimensional simulation of bone tissue adaptation for the human femur. The density distribution of the model shows good agreement with the actual bone density distribution. Permeability at the microstructural level assures interconnectivity of pores, which mimics the interconnectivity of trabecular bone essential for vascularization and transport of nutrients. The importance of this multiscale model relays on the flexibility to control the morphometric parameters that characterize the trabecular structure. Therefore, the presented model can be a valuable tool to define bone quality, to assist with diagnosis of osteoporosis, and to support the development of bone substitutes.     

Analysis of microstructural and mechanical alterations of trabecular bone in a simulated three-dimensional remodeling process

Bone remodeling is a complex dynamic process, which modulates both bone mass and bone microstructure. In addition to bone mass, bone microstructure is an important contributor to bone quality in osteoporosis and fragility fractures. However, the quantitative knowledge of evolution of three-dimensional (3D) trabecular microstructure in adaptation to the external forces is currently limited. In this study, a new 3D simulation method of remodeling of human trabecular bone was developed to quantitatively study the dynamic evolution of bone mass and trabecular microstructure in response to different external loading conditions. The morphological features of trabecular plate and rod, such as thickness and number density in different orientations were monitored during the remodeling process using a novel imaging analysis technique, namely Individual Trabecula Segmentation (ITS). We showed that the volume fraction and microstructures of trabecular bone including, trabecular type and orientation, were determined by the applied mechanical load. Particularly, the morphological parameters of trabecular plates were more sensitive to the applied load, indicating that they played the major role in the mechanical properties of the trabecular bone. Reducing the applied load caused severe microstructural deteriorations of trabecular bone, such as trabecular plate perforation, rod breakage, and a conversion from plates to rods.     

Continuum remodeling revisited

Recent research effort in bone remodeling has been directed toward describing interstitial fluid flow in the lacuno-canalicular system and its potential as a cellular stimulus. Regardless of the precise contents of the mechanotransduction “black box”, it seems clear that the fluid flow on which the remodeling is predicated cannot occur under static loading conditions. In an attempt to help continuum remodeling simulations catch up with cellular and subcellular research, this paper presents a simple, strain rate driven remodeling algorithm for density allocation and principal material direction rotations. An explicit finite element code was written and deployed on a supercomputer which discretizes the remodeling process and uses an objective hypoelastic constitutive law to simulate trabecular realignment. Results indicate that a target strain rate for this dynamic approach is |D _I ^*| = 1.7% per second which seems reasonable when compared to observed strain rates. Simulations indicate that a morpho-mechanically realistic three-dimensional bone can be synthesized by applying a few dynamic loads at the envelope of common daily physiological rates, even with no static loading component.     

Comparison of different hip prosthesis shapes considering micro-level bone remodeling and stress-shielding criteria using three-dimensional design space topology optimization

Since the late 1980s, computational analysis of total hip arthroplasty (THA) prosthesis components has been completed using macro-level bone remodeling algorithms. The utilization of macro-sized elements requires apparent bone densities to predict cancellous bone strength, thereby, preventing visualization and analysis of realistic trabecular architecture. In this study, we utilized a recently developed structural optimization algorithm, design space optimization (DSO), to perform a micro-level three-dimensional finite element bone remodeling simulation on the human proximal femur pre- and post-THA. The computational simulation facilitated direct performance comparison between two commercially available prosthetic implant stems from Zimmer Inc.: the Alloclassic and the Mayo conservative. The novel micro-level approach allowed the unique ability to visualize the trabecular bone adaption post-operation and to quantify the changes in bone mineral content by region. Stress-shielding and strain energy distribution were also quantified for the immediate post-operation and the stably fixated, post-remodeling conditions. Stress-shielding was highest in the proximal region and remained unchanged post-remodeling; conversely, the mid and distal portions show large increases in stress, suggesting a distal shift in the loadpath. The Mayo design conserves bone mass, while simultaneously reducing the incidence of stress-shielding compared to the Alloclassic, revealing a key benefit of the distinctive geometry. Several important factors for stable fixation, determined in clinical evaluations from the literature, were evident in both designs: high levels of proximal bone loss and distal bone densification. The results suggest this novel computational framework can be utilized for comparative hip prosthesis shape, uniquely considering the post-operation bone remodeling as a design criterion.     

A theoretical model to predict distribution of the fabric tensor and apparent density in cancellous bone

 The adaptation of cancellous bone to mechanical forces is well recognized. Theoretical models for predicting cancellous bone architecture have been developed and have mainly focused on the distribution of trabecular mass or the apparent density. The purpose of this study was to develop a theoretical model which can simultaneously predict the distribution of trabecular orthotropy/orientation, as represented by the fabric tensor, along with apparent density. Two sets of equations were derived under the assumption that cancellous bone is a biological self-optimizing material which tends to minimize strain energy. The first set of equations provide the relationship between the fabric tensor and stress tensor, and have been verified to be consistent with Wolff’s law of trabecular architecture, that is, the principal directions of the fabric tensor coincide with the principal stress trajectories. The second set of equations yield the apparent density from the stress tensor, which was shown to be identical to those obtained based on local optimization with strain energy density of true bone tissue as the objective function. These two sets of equations, together with elasticity field equations, provide a complete mathematical formulation for the adaptation of cancellous bone. Received: 25 February 1997/Revised version: 23 September 1997     

Computer simulation of trabecular remodeling in human proximal femur using large-scale voxel FE models: Approach to understanding Wolff's law

Ever since Julius Wolff proposed the law of bone transformation in the 19th century, it has been widely known that the trabecular structure of cancellous bone adapts functionally to the loading environment. To understand the mechanism of Wolff's law, a three-dimensional (3D) computer simulation of trabecular structural changes due to surface remodeling was performed for a human proximal femur. A large-scale voxel finite element model was constructed to simulate the structural changes of individual trabeculae over the entire cancellous region. As a simple remodeling model that considers bone cellular activities regulated by the local mechanical environment, nonuniformity of local stress was assumed to drive the trabecular surface remodeling to seek a uniform stress state. Simulation results demonstrated that cell-scale (～10 μm) remodeling in response to mechanical stimulation created complex 3D trabecular structures of the entire bone-scale (～10 cm), as illustrated in the reference of Wolff. The bone remodeling reproduced the characteristic anisotropic structure in the coronal cross section and the isotropic structures in other cross sections. The principal values and axes of a structure characterized by fabric ellipsoids corresponded to those of the apparent stress of the structure. The proposed large-scale computer simulation indicates that in a complex mechanical environment of a hierarchical bone structure of over 10⁴ length scale (from ～10 μm to ～10 cm), a simple remodeling at the cellular/trabecular levels creates a highly complex and functional trabecular structure, as characterized by bone density and orientation.     

Application of homogenization theory to the study of trabecular bone mechanics.

It is generally accepted that the strength and stiffness of trabecular bone is strongly affected by trabecular microstructure. It has also been hypothesized that stress induced adaptation of trabecular bone is affected by trabecular tissue level stress and/or strain. At this time, however, there is no generally accepted (or easily accomplished) technique for predicting the effect of microstructure on trabecular bone apparent stiffness and strength or estimating tissue level stress or strain. In this paper, a recently developed mechanics theory specifically designed to analyze microstructured materials, called the homogenization theory, is presented and applied to analyze trabecular bone mechanics. Using the homogenization theory it is possible to perform microstructural and continuum analyses separately and then combine them in a systematic manner. Stiffness predictions from two different microstructural models of trabecular bone show reasonable agreement with experimental results, depending on metaphyseal region, (R2 greater than 0.5 for proximal humerus specimens, R2 less than 0.5 for distal femur and proximal tibia specimens). Estimates of both microstructural strain energy density (SED) and apparent SED show that there are large differences (up to 30 times) between apparent SED (as calculated by standard continuum finite element analyses) and the maximum microstructural or tissue SED. Furthermore, a strut and spherical void microstructure gave very different estimates of maximum tissue SED for the same bone volume fraction (BV/TV). The estimates from the spherical void microstructure are between 2 and 20 times greater than the strut microstructure at 10-20% BV/TV.     

Trabecular bone remodeling around smooth and porous implants in an equine patellar model

The objective of this investigation was to examine the stress-morphology relationships for trabecular bone around implants with different surface characteristics. Stainless steel spheres with either a polished surface or a sintered-bead porous coating were implanted unilaterally into equine patellae and maintained for a 6 month period. Stereological methods were used to quantify the trabecular bone morphology and finite element analyses were performed to predict the trabecular bone stresses. In general, the remodeling response around the smooth implants was greater than that around those porous implants that exhibited bone ingrowth. In accordance with these differences, the finite element models predicted greater changes in the stresses adjacent to the smooth implants due to the nonlinear boundary conditions. However, it did not appear that the trajectorial theory, in its simplest form, was applicable to the remodeling induced by the implants. A linear relationship between the change in bone areal density and the change in von Mises effective stress provides support for the hypothesis that the architecture of trabecular bone corresponds to an optimal structure. The results also demonstrated that, under certain circumstances, small changes in the stress state may result in large changes in the principal material orientation.     

Computational study of Wolff's law with trabecular architecture in the human proximal femur using topology optimization

In the field of bone adaptation, it is believed that the morphology of bone is affected by its mechanical loads, and bone has self-optimizing capability; this phenomenon is well known as Wolff's law of the transformation of bone. In this paper, we simulated trabecular bone adaptation in the human proximal femur using topology optimization and quantitatively investigated the validity of Wolff's law. Topology optimization iteratively distributes material in a design domain producing optimal layout or configuration, and it has been widely and successfully used in many engineering fields. We used a two-dimensional micro-FE model with 50 microm pixel resolution to represent the full trabecular architecture in the proximal femur, and performed topology optimization to study the trabecular morphological changes under three loading cases in daily activities. The simulation results were compared to the actual trabecular architecture in previous experimental studies. We discovered that there are strong similarities in trabecular patterns between the computational results and observed data in the literature. The results showed that the strain energy distribution of the trabecular architecture became more uniform during the optimization; from the viewpoint of structural topology optimization, this bone morphology may be considered as an optimal structure. We also showed that the non-orthogonal intersections were constructed to support daily activity loadings in the sense of optimization, as opposed to Wolff's drawing.     

A microstructural finite element simulation of mechanically induced bone formation.

A finite element method to simulate the formation of an interconnected trabectular bone microstructure oriented with respect to applied in vivo mechanical forces is introduced and quantitatively compared to experimental data from a hydraulic bone chamber implant model. Randomly located 45 microm mineralized nodules were used as the initial condition for the model simulations to represent an early stage of intramembranous bone formation. Boundary conditions were applied consistent with the mechanical environment provided by the in vivo bone chamber model. A two-dimensional repair simulation algorithim that incorporated strain energy density (SED), SED gradient, principal strain, or principal strain gradient as the local objective criterion was utilized to simulate the formation of an oriented trabecular bone microstructure. The simulation solutions were convergent, unique, and relatively insensitive to the assumed initial distribution of mineralized nodules. Model predictions of trabecular bone morphology and anisotropy were quantitatively compared to experimental results. All simulations produced structures that qualitatively resembled oriented trabecular bone. However only simulations utilizing a gradient objective criterion yielded results quantitatively similar to in vivo observations. This simulation approach coupled with an experimental model that delivers controlled in vivo mechanical stimuli can be utilized to study the relationship between physical factors and microstructural adaptation during bone repair.     

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.     

Age and gender effects on bone mass density variation: finite elements simulation

Bone remodeling is a physiological process by which bone constantly adapts its structure to changes in long-term loading manifested by interactions between osteoclasts and osteoblasts. This process can be influenced by many local factors, via effects on bone cells differentiation and proliferation, which are produced by bone cells and act in a paracrine or autocrine way. The aim of the current work is to provide mechanobiological finite elements modeling coupling both cellular activities and mechanical behavior in order to investigate age and gender effects on bone remodeling evolution. A series of computational simulations have been performed on a 2D and 3D human proximal femur. An age- and gender-related impacts on bulk density alteration of trabecular bone have been noticed, and the major actors responsible of this phenomenon have been then discussed.     

Localized damage in vertebral bone is most detrimental in regions of high strain energy density

It was hypothesized that damage to bone tissue would be most detrimental to the structural integrity of the vertebral body if it occurred in regions with high strain energy density, and not necessarily in regions of high or low trabecular bone apparent density, or in a particular anatomic location. The reduction in stiffness due to localized damage was computed in 16 finite element models of 10-mm-thick human vertebral sections. Statistical analyses were performed to determine which characteristic at the damage location--strain energy density, apparent density, or anatomic location--best predicted the corresponding stiffness reduction. There was a strong positive correlation between regional strain energy density and structural stiffness reduction in all 16 vertebral sections for damage in the trabecular centrum (p < 0.05, r2 = 0.43-0.93). By contrast, regional apparent density showed a significant negative correlation to stiffness reduction in only four of the sixteen bones (p < 0.05, r2 = 0.47-0.58). While damage in different anatomic locations did lead to different reductions in stiffness (p < 0.0001, ANOVA), no single location was consistently the most critical location for damage. Thus, knowledge of the characteristics of bone that determine strain energy density distributions can provide an understanding of how damage reduces whole bone mechanical properties. A patient-specific finite element model displaying a map of strain energy density can help optimize surgical planning and reinforcement of bone in individuals with high fracture risk.     

Mechanics of linear microcracking in trabecular bone

Microcracking in trabecular bone is responsible both for the mechanical degradation and remodeling of the trabecular bone tissue. Recent results on trabecular bone mechanics have demonstrated that bone tissue microarchitecture, tissue elastic heterogeneity and tissue-level mechanical anisotropy all should be considered to obtain detailed information on the mechanical stress state. The present study investigated the influence of tissue microarchitecture, tissue heterogeneity in elasticity and material separation properties and tissue-level anisotropy on the microcrack formation process. Microscale bone models were executed with the extended finite element method. It was demonstrated that anisotropy and heterogeneity of the bone tissue contribute significantly to bone tissue toughness and the resistance of trabecular bone to microcrack formation. The compressive strain to microcrack initiation was computed to increase by a factor of four from an assumed homogeneous isotropic tissue to an assumed anisotropic heterogenous tissue.     

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.