Trabecular surface remodeling simulation for cancellous bone using microstructural voxel finite element models

A computational simulation method for three-dimensional trabecular surface remodeling was proposed, using voxel finite element models of cancellous bone, and was applied to the experimental data. In the simulation, the trabecular microstructure was modeled based on digital images, and its morphological changes due to surface movement at the trabecular level were directly expressed by removing/adding the voxel elements from/to the trabecular surface. A remodeling simulation at the single trabecular level under uniaxial compressive loading demonstrated smooth morphological changes even though the trabeculae were modeled with discrete voxel elements. Moreover, the trabecular axis rotated toward the loading direction with increasing stiffness, simulating functional adaptation to the applied load. In the remodeling simulation at the trabecular structural level, a cancellous bone cube was modeled using a digital image obtained by microcomputed tomography (microCT), and was uniaxially compressed. As a result, the apparent stiffness against the applied load increased by remodeling, in which the trabeculae reoriented to the loading direction. In addition, changes in the structural indices of the trabecular architecture coincided qualitatively with previously published experimental observations. Through these studies, it was demonstrated that the newly proposed voxel simulation technique enables us to simulate the trabecular surface remodeling and to compare the results obtained using this technique with the in vivo experimental data in the investigation of the adaptive bone remodeling phenomenon.     

A model of bone adaptation as an optimization process

The internal bone adaptation of the proximal femur is considered. A three-dimensional finite element model of the proximal femur is used. The bone remodeling in this work is numerically described by an evolutionary remodeling scheme with anisotropic material parameters and time-dependent loading. The memory of past loading is included in the model to account for the delay in the bone response from the load changes. The remodeling rate equation is derived from the structural optimization task of maximizing the stiffness in each time step. Additional information can be extracted from the optimization process; the remodeling equilibrium parameter where no apposition or resorption takes place, is defined as the element optimality conditions and the optimal design is used as an initial design for the onset of the remodeling simulation. Two examples of bone adaptation resulting from load changes are given, and the irreversible nature of the model is illustrated.     

Computer simulation of orthodontic tooth movement using CT image-based voxel finite element models with the level set method

Orthodontic tooth movement (OTM) is an adaptive biomechanical response of dentoalveolar components to orthodontic forces, in which remodeling of the alveolar bone occurs in response to changes in the surrounding mechanical environment. In this study, we developed a framework for OTM simulation by combining an image-based voxel finite element method, with a surface-tracking level set method using three-dimensional computer models. For a case study to demonstrate its capability of expressing clinical tooth movement, we observed displacement and rotation of the tooth under three types of force conditions. The simulation results demonstrate that the proposed simulation method has the potential to predict clinical OTM.     

Osteocyte-viability-based simulations of trabecular bone loss and recovery in disuse and reloading

Osteocyte apoptosis is known to trigger targeted bone resorption. In the present study, we developed an osteocyte-viability-based trabecular bone remodeling (OVBR) model. This novel remodeling model, combined with recent advanced simulation methods and analysis techniques, such as the element-by-element 3D finite element method and the ITS technique, was used to quantitatively study the dynamic evolution of bone mass and trabecular microstructure in response to various loading and unloading conditions. Different levels of unloading simulated the disuse condition of bed rest or microgravity in space. The amount of bone loss and microstructural deterioration correlated with the magnitude of unloading. The restoration of bone mass upon the reloading condition was achieved by thickening the remaining trabecular architecture, while the lost trabecular plates and rods could not be recovered by reloading. Compared to previous models, the predictions of bone resorption of the OVBR model are more consistent with physiological values reported from previous experiments. Whereas osteocytes suffer a lack of loading during disuse, they may suffer overloading during the reloading phase, which hampers recovery. The OVBR model is promising for quantitative studies of trabecular bone loss and microstructural deterioration of patients or astronauts during long-term bed rest or space flight and thereafter bone recovery.     

Capturing microscopic features of bone remodeling into a macroscopic model based on biological rationales of bone adaptation

To understand Wolff’s law, bone adaptation by remodeling at the cellular and tissue levels has been discussed extensively through experimental and simulation studies. For the clinical application of a bone remodeling simulation, it is significant to establish a macroscopic model that incorporates clarified microscopic mechanisms. In this study, we proposed novel macroscopic models based on the microscopic mechanism of osteocytic mechanosensing, in which the flow of fluid in the lacuno-canalicular porosity generated by fluid pressure gradients plays an important role, and theoretically evaluated the proposed models, taking biological rationales of bone adaptation into account. The proposed models were categorized into two groups according to whether the remodeling equilibrium state was defined globally or locally, i.e., the global or local uniformity models. Each remodeling stimulus in the proposed models was quantitatively evaluated through image-based finite element analyses of a swine cancellous bone, according to two introduced criteria associated with the trabecular volume and orientation at remodeling equilibrium based on biological rationales. The evaluation suggested that nonuniformity of the mean stress gradient in the local uniformity model, one of the proposed stimuli, has high validity. Furthermore, the adaptive potential of each stimulus was discussed based on spatial distribution of a remodeling stimulus on the trabecular surface. The theoretical consideration of a remodeling stimulus based on biological rationales of bone adaptation would contribute to the establishment of a clinically applicable and reliable simulation model of bone remodeling.     

Evaluation of bone remodeling around single dental implants of different lengths: a mechanobiological numerical simulation and validation using clinical data

Algorithmic models have been proposed to explain adaptive behavior of bone to loading; however, these models have not been applied to explain the biomechanics of short dental implants. Purpose of present study was to simulate bone remodeling around single implants of different lengths using mechanoregulatory tissue differentiation model derived from the Stanford theory, using finite elements analysis (FEA) and to validate the theoretical prediction with the clinical findings of crestal bone loss. Loading cycles were applied on 7-, 10-, or 13-mm-long dental implants to simulate daily mastication and bone remodeling was assessed by changes in the strain energy density of bone after a 3, 6, and 12 months of function. Moreover, clinical findings of marginal bone loss in 45 patients rehabilitated with same implant designs used in the simulation (n = 15) were computed to validate the theoretical results. FEA analysis showed that although the bone density values reduced over time in the cortical bone for all groups, bone remodeling was independent of implant length. Clinical data showed a similar pattern of bone resorption compared with the data generated from mathematical analyses, independent of implant length. The results of this study showed that the mechanoregulatory tissue model could be employed in monitoring the morphological changes in bone that is subjected to biomechanical loads. In addition, the implant length did not influence the bone remodeling around single dental implants during the first year of loading.     

Computational study on the effect of loading alteration caused by disc degeneration on the trabecular architecture in human lumbar spine

A thorough understanding of age-related phenomena on the trabecular architecture in the human lumbar spine can help the diagnosis and prognosis of age-related architectural changes, and provide an insight into the corresponding clinical assessments. In this paper we considered the different loading conditions of the young and old lumbar spines mainly caused by disc degeneration and studied the effect of loading alteration on trabecular architecture in lumbar spines. A two-dimensional μFE models with a 40 μm pixel resolution were built to represent the full trabecular architecture in the human lumbar spine, and a topology optimization with the aid of finite element method was conducted to numerically investigate the trabecular morphological changes. Topology optimization iteratively distributes material in a design domain producing optimal layout or configuration, and it has been widely and successfully used for the study of bone remodeling. As a result of adaptive response of bone remodeling due to different loading conditions, we obtained two distinctively different trabecular architectures for the young and old lumbar spines, and we observed a strong correlation between our numerical results and the actual trabecular architecture in the literature. The proposed numerical framework and results demonstrated the potential use of the topology optimization-based numerical tool for putative treatments in advance of actual clinical procedures for the patients.     

Adaptive glenoid bone remodeling simulation

Glenoid prosthesis loosening is the most common cause for revision total shoulder arthroplasty. Stress-induced bone remodeling may compromise long-term prosthesis fixation and significantly contribute to loosening. Realistic, robust analysis of bone-prosthesis constructs need to look beyond initial post-implantation mechanics provided by static finite element (FE) simulation. Adaptive bone remodeling simulations based on Wolff's law are needed for evaluating long-term glenoid prostheses fixation. The purpose of this study was to take a first step towards this goal and create and validate two-dimensional FE simulations, using the intact glenoid, for computing subject-specific adaptive glenoid remodeling. Two-dimensional glenoid FE models were created from scapulae computed tomography images. Two distinct processes, “element” and “node” simulations, used the forward-Euler method to compute bone remodeling. Initial bone density was homogeneous. Center and offset load combinations were iteratively applied. To validate the simulations we performed location-specific statistical comparisons between predicted and actual bone density, load combinations, and “element” and “node” processes. Visually and quantitatively “element” simulations produced better results (p>0.22), and correlation coefficients ranged 0.51–0.69 (p<0.001). Having met this initial work's goals, we expect subject-specific FE glenoid bone remodeling simulations together with static FE stress analyses to be effective tools for designing and evaluating glenoid prostheses.     

The behavior of adaptive bone-remodeling simulation models

The process of adaptive bone remodeling can be described mathematically and simulated in a computer model, integrated with the finite element method. In the model discussed here, cortical and trabecular bone are described as continuous materials with variable density. The remodeling rule applied to simulate the remodeling process in each element individually is, in fact, an objective function for an optimization process, relative to the external load. Its purpose is to obtain a constant, preset value for the strain energy per unit bone mass, by adapting the density. If an element in the structure cannot achieve that, it either turns to its maximal density (cortical bone) or resorbs completely.

It is found that the solution obtained in generally a discontinuous patchwork. For a two-dimensional proximal femur model this patchwork shows a good resemblance with the density distribution of a real proximal femur.

It is shown that the discontinuous end configuration is dictated by the nature of the differential equations describing the remodeling process. This process can be considered as a nonlinear dynamical system with many degrees of freedom, which behaves divergent relative to the objective, leading to many possible solutions. The precise solution is dependent on the parameters in the remodeling rule, the load and the initial conditions. The feedback mechanism in the process is self-enhancing; denser bone attracts more strain energy, whereby the bone becomes even more dense. It is suggested that this positive feedback of the attractor state (the strain energy field) creates order in the end configuration. In addition, the process ensures that the discontinuous end configuration is a structure with a relatively low mass, perhaps a minimal-mass structure, although this is no explicit objective in the optimization process.

It is hypothesized that trabecular bone is a chaotically ordered structure which can be considered as a fractal with characteristics of optimal mechanical resistance and minimal mass, of which the actual morphology depends on the local (internal) loading characteristics, the sensor-cell density and the degree of mineralization.     

Bone morphology allows estimation of loading history in a murine model of bone adaptation

Bone adapts its morphology (density/micro- architecture) in response to the local loading conditions in such a way that a uniform tissue loading is achieved (‘Wolff’s law’). This paradigm has been used as a basis for bone remodeling simulations to predict the formation and adaptation of trabecular bone. However, in order to predict bone architectural changes in patients, the physiological external loading conditions, to which the bone was adapted, need to be determined. In the present study, we developed a novel bone loading estimation method to predict such external loading conditions by calculating the loading history that produces the most uniform bone tissue loading. We applied this method to murine caudal vertebrae of two groups that were in vivo loaded by either 0 or 8 N, respectively. Plausible load cases were sequentially applied to micro-finite element models of the mice vertebrae, and scaling factors were calculated for each load case to derive the most uniform tissue strain-energy density when all scaled load cases are applied simultaneously. The bone loading estimation method was able to predict the difference in loading history of the two groups and the correct load magnitude for the loaded group. This result suggests that the bone loading history can be estimated from its morphology and that such a method could be useful for predicting the loading history for bone remodeling studies or at sites where measurements are difficult, as in bone in vivo or fossil bones.     

Interstitial fluid flow in canaliculi as a mechanical stimulus for cancellous bone remodeling: in silico validation

Cancellous bone has a dynamic 3-dimensional architecture of trabeculae, the arrangement of which is continually reorganized via bone remodeling to adapt to the mechanical environment. Osteocytes are currently believed to be the major mechanosensory cells and to regulate osteoclastic bone resorption and osteoblastic bone formation in response to mechanical stimuli. We previously developed a mathematical model of trabecular bone remodeling incorporating the possible mechanisms of cellular mechanosensing and intercellular communication in which we assumed that interstitial fluid flow activates the osteocytes to regulate bone remodeling. While the proposed model has been validated by the simulation of remodeling of a single trabecula, it remains unclear whether it can successfully represent in silico the functional adaptation of cancellous bone with its multiple trabeculae. In the present study, we demonstrated the response of cancellous bone morphology to uniaxial or bending loads using a combination of our remodeling model with the voxel finite element method. In this simulation, cancellous bone with randomly arranged trabeculae remodeled to form a well-organized architecture oriented parallel to the direction of loading, in agreement with the previous simulation results and experimental findings. These results suggested that our mathematical model for trabecular bone remodeling enables us to predict the reorganization of cancellous bone architecture from cellular activities. Furthermore, our remodeling model can represent the phenomenological law of bone transformation toward a locally uniform state of stress or strain at the trabecular level.     

The effect of muscle loading on the simulation of bone remodelling in the proximal femur

A large number of finite element analyses of the proximal femur rely on a simplified set of muscle and joint contact loads to represent the boundary conditions of the model. In the context of bone remodelling analysis around hip implants, muscle loading affects directly the spatial distribution of the remodelling signal. In the present study we performed a sensitivity analysis on the effect of different muscle loading configurations on the outcome of the bone remodelling simulation. An anatomical model of the femur with the implanted stem in place was constructed using the CT data of the Visible Human Project dataset of the National Institute of Health. The model was loaded with three muscle force configurations with increasing level of complexity. A strain adaptive remodelling rule was employed to simulate the post-operative bone changes around the implant stem and the results of the simulation were assessed quantitatively in terms of the bone mineral content changes in 18 periprosthetic regions of interest. The results showed considerable differences in the amount of bone loss predicted between the three cases. The simplified models generally predicted more pronounced bone loss. Although the overall remodelling patterns observed were similar, the bone conserving effect of additional muscle forces in the vicinity of their areas of attachment was clear. The results of this study suggest that the loading configuration of the FE model does play an important role in the outcome of the remodelling simulation.     

Changes in the Fabric and Compliance Tensors of Cancellous Bone due to Trabecular Surface Remodeling,Predicted by a Digital Image-based Model

Fabric and compliance tensors of a cube of cancellous bone with a complicated three-dimensional trabecular structure were obtained for trabecular surface remodeling by using a digital image-based model combined with a large-scale finite element method. Using mean intercept length and a homogenization method, the fabric and compliance tensors were determined for the trabecular structure obtained in the computer remodeling simulation. The tensorial quantities obtained indicated that anisotropic structural changes occur in cancellous bone adapting to the compressive loading condition. There were good correlations between the fabric tensor, bone volume fraction, and compliance tensor in the remodeling process. The result demonstrates that changes in the structural and mechanical properties of cancellous bone are essentially anisotropic and should be expressed by tensorial quantities.     

Changes in the fabric and compliance tensors of cancellous bone due to trabecular surface remodeling, predicted by a digital image-based model

Fabric and compliance tensors of a cube of cancellous bone with a complicated three-dimensional trabecular structure were obtained for trabecular surface remodeling by using a digital image-based model combined with a large-scale finite element method. Using mean intercept length and a homogenization method, the fabric and compliance tensors were determined for the trabecular structure obtained in the computer remodeling simulation. The tensorial quantities obtained indicated that anisotropic structural changes occur in cancellous bone adapting to the compressive loading condition. There were good correlations between the fabric tensor, bone volume fraction, and compliance tensor in the remodeling process. The result demonstrates that changes in the structural and mechanical properties of cancellous bone are essentially anisotropic and should be expressed by tensorial quantities.     

Periosteal and endosteal control of bone remodeling under torsional loading.

The shape changes that occur in the mid-diaphysis of a long bone due to adaptive remodeling induced by increasing or decreasing the axial and/or torsional loading of the bone are investigated using a simple model. In this model the mid-diaphysis of a long bone is represented as a hollow thick-walled right-circular cylinder, and different optimal strategies for bone remodeling are considered. It is shown that if such a thick-walled right-circular cylinder capable of surface remodeling is subjected to an axial compressive load and a twisting torque, then the remodeling patterns depend on whether the periosteal surface or the endosteal surface controls the limits of the remodeling process. It is shown that the effect of increasing the torque is always opposite to the effect of increasing the compressive load. Thus, similar remodeling patterns are obtained by increasing one type of loading and decreasing the other. Aside from the restriction of idealized cylindrical geometry, the only assumptions made are that the bone tissue is linearly elastic and that there exists a finite range of remodeling equilibrium stresses. Only those loading situations which maintain the bone in remodeling equilibrium are considered in this work. It follows that the results presented are independent of the specific type of rule governing the temporal evolution of the bone shape, since any such rule applies only in situations where there is active remodeling and, hence, no remodeling equilibrium.     

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.     

Activity intensity,assistive devices and joint replacement influence predicted remodelling in the proximal femur

Bone morphology and density changes are commonly observed following joint replacement, may contribute to the risks of implant loosening and periprosthetic fracture and reduce the available bone stock for revision surgery. This study was presented in the ‘Bone and Cartilage Mechanobiology across the scales’ WCCM symposium to review the development of remodelling prediction methods and to demonstrate simulation of adaptive bone remodelling around hip replacement femoral components, incorporating intrinsic (prosthesis) and extrinsic (activity and loading) factors. An iterative bone remodelling process was applied to finite element models of a femur implanted with a cementless total hip replacement (THR) and a hip resurfacing implant. Previously developed for a cemented THR implant, this modified process enabled the influence of pre- to post-operative changes in patient activity and joint loading to be evaluated. A control algorithm used identical pre- and post-operative conditions, and the predicted extents and temporal trends of remodelling were measured by generating virtual X-rays and DXA scans. The modified process improved qualitative and quantitative remodelling predictions for both the cementless THR and resurfacing implants, but demonstrated the sensitivity to DXA scan region definition and appropriate implant–bone position and sizing. Predicted remodelling in the intact femur in response to changed activity and loading demonstrated that in this simplified model, although the influence of the extrinsic effects were important, the mechanics of implantation were dominant. This study supports the application of predictive bone remodelling as one element in the range of physical and computational studies, which should be conducted in the preclinical evaluation of new prostheses.     

The effects of misalignment during in vivo loading of bone: Techniques to detect the proximity of objects in three-dimensional models

Theories of mechanical adaptation of bone suggest that mechanical loading causes bone formation at discrete locations within bone microstructure experiencing the greatest mechanical stress/strain. Experimental testing of such theories requires in vivo loading experiments and high-resolution finite element models to determine the distribution of mechanical stresses. Finite element models of in vivo loading experiments typically assume idealized boundary conditions with applied load perfectly oriented on the bone, however small misalignments in load orientation during an in vivo experiment are unavoidable, and potentially confound the ability of finite element models to predict locations of bone formation at the scale of micrometers. Here we demonstrate two different three-dimensional spatial correlation methods to determine the effects of misalignment in load orientation on the locations of high mechanical stress/strain in the rodent tail loading model. We find that, in cancellous bone, the locations of tissue with high stress are maintained under reasonable misalignments in load orientation (p<0.01). In cortical bone, however, angular misalignments in the dorsal direction can alter the locations of high mechanical stress, but the locations of tissue with high stress are maintained under other misalignments (p<0.01). We conclude that, when using finite element models of the rodent tail loading model, small misalignments in loading orientation do not affect the predicted locations of high mechanical stress within cancellous bone.     

Effects of Materials of Cementless Femoral Stem on the Functional Adaptation of Bone

The objective of this paper is to identify the effects of materials of cementless femoral stem on the functional adaptive behaviors of bone.The remodeling behaviors of a two-dimensional simplified model of cementless hip prosthesis with stiff stem,flexible 'iso-elastic' stem,one-dimensional Functionally Graded Material (FGM) stem and two-dimensional FGM stem for the period of four years after prosthesis replacement were quantified by incorporating the bone remodeling algorithm with finite element analysis.The distributions of bone density,von Mises stress,and interface shear stress were obtained.The results show that two-dimensional FGM stem may produce more mechanical stimuli and more uniform interface shear stress compared with the stems made of other materials,thus the host bone is well preserved.Accordingly,the two-dimensional FGM stem is an appropriate femoral implant from a biomechanical point of view.The numerical simulation in this paper can provide a quantitative computational paradigm for the changes of bone morphology caused by implants,which can help to improve the design of implant to reduce stress shielding and the risk of bone-prosthesis interface failure.     

Sensitivity of damage predictions to tissue level yield properties and apparent loading conditions

High-resolution finite element models of trabecular bone failure could be used to augment current techniques for measuring damage in trabecular bone. However, the sensitivity of such models to the assumed tissue yield properties and apparent loading conditions is unknown. The goal of this study was to assess the sensitivity of the amount and mode (tension vs. compression) of tissue level yielding in trabecular bone to these factors. Linear elastic, high-resolution finite element models of nine bovine tibial trabecular bone specimens were used to calculate the fraction of the total tissue volume that exceeded each criterion for apparent level loading to the reported elastic limit in both on-axis and transverse compression and tension, and in shear. Four candidate yield criteria were studied, based on values suggested in the literature. Both the amount and the failure mode of yielded tissue were sensitive to the magnitudes of the tissue yield strains, the degree of tension-compression asymmetry of the yield criterion, and the applied apparent loads. The amount of yielded tissue was most sensitive to the orientation of the applied apparent loading, with the most tissue yielding for loading along the principal trabecular orientation and the least for loading perpendicular to it, regardless of the assumed tissue level yield criterion. Small changes in the magnitudes and the degree of asymmetry of the tissue yield criterion resulted in much larger changes in the amount of yielded tissue in the model. The results indicate that damage predictions based on high-resolution finite element models are highly sensitive to the assumed tissue yield properties. As such, good estimates of these values are needed before high-resolution finite element models can be applied to the study of trabecular bone damage. Regardless of the assumed tissue yield properties, the amount and type of damage that occurs in trabecular bone depends on the relative orientations of the applied apparent loads to the trabecular architecture, and this parameter should be controlled for both experimental and computational damage studies.