3D strain map of axially loaded mouse tibia: a numerical analysis validated by experimental measurements

A combined experimental/numerical study was performed to calculate the 3D octahedral shear strain map in a mouse tibia loaded axially. This study is motivated by the fact that the bone remodelling analysis, in this in vivo mouse model should be performed at the zone of highest mechanical stimulus to maximise the measured effects. Accordingly, it is proposed that quantification of bone remodelling should be performed at the tibial crest and at the distal diaphysis. The numerical model could also be used to furnish a more subtle analysis as a precise correlation between local strain and local biological response can be obtained with the experimentally validated numerical model.     

An analytical and numerical study of the stability of bone remodelling theories: dependence on microstructural stimulus.

The origin of unstable bone remodelling simulations using strain-energy-based remodelling rules was studied mathematically in order to assess whether the unstable behavior was due to the mathematical rules proposed to characterize the processes, or to the numerical approximations used to exercise the mathematical predictions. A condition which is necessary for the stability of a strain-energy-based remodelling theory was derived analytically using the calculus of variation. The analytical result was derived using a simple elastic model which consists of a long beam loaded by an axial force and a bending moment. This loading situation mimics the coupling between local density and global density distributions seen in vivo. A condition necessary for a stable remodelling scheme is arrived at, but the conditions necessary to guarantee a stable remodelling scheme are not. In this remodelling scheme, the elastic modulus is proportional to volumetric density raised to an exponent n, and the microstructural stimulus is taken as the strain energy density divided by volumetric density raised to an exponent m. In order for a remodelling scheme to be stable in this loading situation, m must be greater than n. Finite-difference time-stepping is used to verify the predictions of the analytical study. These numerical studies appear to confirm the analytical studies. Physiologic interpretation of the behavior found with n greater than m indicates that this type of unstable behavior is unlikely to be observed in vivo. Since numerical approximations are not made in deriving this stability condition, we conclude that the mathematical rules proposed to characterize bone remodelling based on strain energy density should meet this condition to be relevant to physiologic bone remodelling.     

Comparison of strain measurement in the mouse forearm using subject-specific finite element models,strain gaging,and digital image correlation

Mechanical loading in bone leads to the activation of bone-forming pathways that are most likely associated with a minimum strain threshold being experienced by the osteocyte. To investigate the correlation between cellular response and mechanical stimuli, researchers must develop accurate ways to measure/compute strain both externally on the bone surface and internally at the osteocyte level. This study investigates the use of finite element (FE) models to compute bone surface strains on the mouse forearm. Strains from three FE models were compared to data collected experimentally through strain gaging and digital image correlation (DIC). Each FE model was assigned subject-specific bone properties and consisted of one-dimensional springs representing the interosseous membrane. After three-point bending was performed on the ulnae and radii, moment of inertia was determined from microCT analysis of the bone region between the supports and then used along with standard beam analyses to calculate the Young’s modulus. Non-contact strain measurements from DIC were determined to be more suitable for validating numerical results than experimental data obtained through conventional strain gaging. When comparing strain responses in the three ulnae, we observed a 3–14% difference between numerical and DIC strains while the strain gage values were 37–56% lower than numerical values. This study demonstrates a computational approach for capturing bone surface strains in the mouse forearm. Ultimately, strains from these macroscale models can be used as inputs for microscale and nanoscale FE models designed to analyze strains directly in the osteocyte lacunae.     

Numerical study of tension and strain distribution around rat molars]

A knowledge of the mechanical processes triggered in the bone and periodontal ligament (PDL) by orthodontic forces applied to a tooth is of decisive importance for an understanding of the subsequent remodelling around the tooth. To investigate these mechanical relationships, three-dimensional finite element (FE) models of the first lower molar in the rat were established. On the basis of digitized serial histological sections, these FE models were generated semi-automatically. Using various simplified geometrical variations, an appropriate FE model for the analysis of the stress and strain distributions was established. The numerical analyses were carried out under a mesially directed force of 0.1 N. Stress distributions in the bone and PDL showed a similar pattern, while strains in the bone were lower than in the PDL by a factor of 10-5. The data confirm the assumption that strain patterns in the PDL may be the key stimulus of bone remodelling.     

Correlation between pre-operative periprosthetic bone density and post-operative bone loss in THA can be explained by strain-adaptive remodelling.

Periprosthetic adaptive bone remodelling after total hip arthroplasty can be simulated in computer models, combining bone remodelling theory with finite element analysis. Patient specific three-dimensional finite element models of retrieved bone specimens from an earlier bone densitometry (DEXA) study were constructed and bone remodelling simulations performed. Results of the simulations were analysed both qualitatively and quantitatively. Patterns of predicted bone loss corresponded very well with the DEXA measurements on the retrievals. The amount of predicted bone loss, measured quantitatively by simulating DEXA on finite element models, was found to be inversely correlated with the initial bone mineral content. It was concluded that the same clinically observed correlation can therefore be explained by mechanically induced remodelling. This finding extends the applicability of numerical pre-clinical testing to the analysis of interaction between implant design and initial state of the bone.     

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.     

Interstitial fluid flow in the osteon with spatial gradients of mechanical properties: a finite element study

Bone remodelling is the process that maintains bone structure and strength through adaptation of bone tissue mechanical properties to applied loads. Bone can be modelled as a porous deformable material whose pores are filled with cells, organic material and interstitial fluid. Fluid flow is believed to play a role in the mechanotransduction of signals for bone remodelling. In this work, an osteon, the elementary unit of cortical bone, is idealized as a hollow cylinder made of a deformable porous matrix saturated with an interstitial fluid. We use Biot’s poroelasticity theory to model the mechanical behaviour of bone tissue taking into account transverse isotropic mechanical properties. A finite element poroelastic model is developed in the COMSOL Multiphysics software. Elasticity equations and Darcy’s law are implemented in this software; they are coupled through the introduction of an interaction term to obtain poroelasticity equations. Using numerical simulations, the investigation of the effect of spatial gradients of permeability or Poisson’s ratio is performed. Results are discussed for their implication on fluid flow in osteons: (i) a permeability gradient affects more the fluid pressure than the velocity profile; (ii) focusing on the fluid flow, the key element of loading is the strain rate; (iii) a Poisson’s ratio gradient affects both fluid pressure and fluid velocity. The influence of textural and mechanical properties of bone on mechanotransduction signals for bone remodelling is also discussed.     

Multiscale methodology for bone remodelling simulation using coupled finite element and neural network computation

The aim of this paper is to develop a multiscale hierarchical hybrid model based on finite element analysis and neural network computation to link mesoscopic scale (trabecular network level) and macroscopic (whole bone level) to simulate the process of bone remodelling. As whole bone simulation, including the 3D reconstruction of trabecular level bone, is time consuming, finite element calculation is only performed at the macroscopic level, whilst trained neural networks are employed as numerical substitutes for the finite element code needed for the mesoscale prediction. The bone mechanical properties are updated at the macroscopic scale depending on the morphological and mechanical adaptation at the mesoscopic scale computed by the trained neural network. The digital image-based modelling technique using μ-CT and voxel finite element analysis is used to capture volume elements representativeof 2 mm³ at the mesoscale level of the femoral head. The input data for the artificial neural network are a set of bone material parameters, boundary conditions and the applied stress. The output data are the updated bone properties and some trabecular bone factors. The current approach is the first model, to our knowledge, that incorporates both finite element analysis and neural network computation to rapidly simulate multilevel bone adaptation.     

Trabecular bone fracture healing simulation with finite element analysis and fuzzy logic

Trabecular bone fractures heal through intramembraneous ossification. This process differs from diaphyseal fracture healing in that the trabecular marrow provides a rich vascular supply to the healing bone, there is very little callus formation, woven bone forms directly without a cartilage intermediary, and the woven bone is remodelled to form trabecular bone. Previous studies have used numerical methods to simulate diaphyseal fracture healing or bone remodelling, however not trabecular fracture healing, which involves both tissue differentiation and trabecular formation. The objective of this study was to determine if intramembraneous bone formation and remodelling during trabecular bone fracture healing could be simulated using the same mechanobiological principles as those proposed for diaphyseal fracture healing. Using finite element analysis and the fuzzy logic for diaphyseal healing, the model simulated formation of woven bone in the fracture gap and subsequent remodelling of the bone to form trabecular bone. We also demonstrated that the trabecular structure is dependent on the applied loading conditions. A single model that can simulate bone healing and remodelling may prove to be a useful tool in predicting musculoskeletal tissue differentiation in different vascular and mechanical environments.     

Trabecular bone remodelling under pathological conditions based on biochemical and mechanical processes involved in BMU activity

In adulthood, bone tissue is continuously renewed by processes governed by basic multicellular units composed of osteocytes, osteoclasts and osteoblasts, which are subjected to local mechanical loads. Osteocytes are known to be integrated mechanosensors that regulate the activation of the osteoclasts and osteoblasts involved in bone resorption and apposition processes, respectively. After collagen tissue apposition, a process of collagen mineralisation takes place, gradually increasing the effective stiffness of bone. This study presents a new model based on physicochemical parameters involved in spongy bone remodelling under pathological conditions. Our model simulates the transient evolution of both geometry and effective Young's modulus of the trabeculae, also taking turnover into account. Various loads were applied on a trabecula in order to determine the evolution of bone volume fraction under pathological conditions. A parametric study performed on the model showed that one key parameter here is the kinetic constant of hydroxyapatite crystallisation. We subsequently tested our model on a pathological case approaching osteoporosis, involving a decrease in the number of viable osteocytes present in bone. The model converges to a lower value ( ? 5%) for bone volume fraction than with a normal quantity of osteocytes. This useful tool offers new perspectives for predicting bone remodelling deficits on a local scale in patients with pathological conditions such as osteoporosis and in bedridden patients, as well as for astronauts subjected to weightlessness in space.     

Maxillary expansion treatment using bone anchors: development and validation of a 3D finite element model

Objective: Develop a finite element (FE) model of a skull to perform biomechanical studies of maxillary expansion using bone anchors (BA).

Materials and methods: A skull model was developed and assigned material properties based on Hounsfield unit (HU) values of cone-beam computerized tomography (CBCT) images. A 3 mm diameter cylindrical BA was modelled and inserted in the palatal bone. A 4 mm transverse displacement was applied on the anchor. An evaluation on the effect on local stresses of BA implantation inclination angle was performed.

Results: Proper displacement results and strain–stress trends for the expansion process were present. Stress distribution patterns were similar as reported in the literature. No significant difference between BA inclination angles was found.

Conclusion: This work leads to a better understanding and prediction of craniofacial and maxillary bone remodelling during ME with BA treatments and is a first step towards the development of patient specific treatments.     

A meshless microscale bone tissue trabecular remodelling analysis considering a new anisotropic bone tissue material law

In this work, a novel anisotropic material law for the mechanical behaviour of the bone tissue is proposed. This new law, based on experimental data, permits to correlate the bone apparent density with the obtained level of stress. Combined with the proposed material law, a biomechanical model for predicting bone density distribution was developed, based on the assumption that the bone structure is a gradually self-optimising anisotropic biological material that maximises its own structural stiffness. The strain and the stress field required in the iterative remodelling process are obtained by means of an accurate meshless method, the Natural Neighbour Radial Point Interpolation Method (NNRPIM). Comparing with other numerical approaches, the inclusion of the NNRPIM presents numerous advantages such as the high accuracy and the smoother stress and strain field distribution. The natural neighbour concept permits to impose organically the nodal connectivity and facilitates the analysis of convex boundaries and extremely irregular meshes. The viability and efficiency of the model were tested on several trabecular benchmark patch examples. The results show that the pattern of the local bone apparent density distribution and the anisotropic bone behaviour predicted by the model for the microscale analysis are in good agreement with the expected structural architecture and bone apparent density distribution.     

Comparative analysis of numerical integration schemes of density equation for a computational model of bone remodelling

In this study, a computational model of bone remodelling problem as proposed by Weinans et al. (1992) is described and solved by other temporal integration techniques different from the Euler scheme. This model considers three types of numerical integration schemes of the evolution of the material density during the remodelling: Euler, Heun and Runge-Kutta methods. Also the strain and the density field are obtained inside each element, at Gauss points or at the nodes of the mesh. A square plate with 1.00?m of side subjected to non-uniform pressure is simulated with two meshes of quadrilateral element with size [Formula: see text] and [Formula: see text]?m. Two increments time size: [Formula: see text] and [Formula: see text] days are used. The results show that Euler, Heun and Runge-Kutta's methods correctly approached the problem of bone remodelling and that there were no appreciable differences in the patterns obtained by the mesh and time step used. In contrast, using an element-based approach and node-based approach, substantial differences were produced in bone remodelling density pattern. 'Chess board' type discontinuities were found in the element approach near the applied pressure area, as were well-defined columns away from this. The node-based approach showed continuity in density distribution. These patterns were well represented by the methods for resolving the density equation. This study concluded that any method of time integration could be used for these meshes and time steps size.     

Comparative analysis of bone remodelling models with respect to computerised tomography-based finite element models of bone

Subject-specific finite element models are an extensively used tool for the numerical analysis of the biomechanical behaviour of human bones. However, bone modelling is not an easy task due to the complex behaviour of bone tissue, involving non-homogeneous and anisotropic mechanical properties. Moreover, bone is a living tissue and therefore its microstructure and mechanical properties evolve with time in a known process called bone remodelling. This phenomenon has been widely studied, many being the numerical models that have been formulated to predict density distribution and its evolution in several bones. The aim of the present study is to assess the capability of a bone remodelling model to predict the bone density distribution of different types of human bone (femur, tibia and mandible) comparing the obtained results with the bone density estimated by means of computerised tomography. Good accuracy was observed for the bone remodelling predictions including the thickness of the cortical layer.     

A full body musculoskeletal model based on flexible multibody simulation approach utilised in bone strain analysis during human locomotion

Load-induced strains applied to bone can stimulate its development and adaptation. In order to quantify the incident strains within the skeleton, in vivo implementation of strain gauges on the surfaces of bone is typically used. However, in vivo strain measurements require invasive methodology that is challenging and limited to certain regions of superficial bones only such as the anterior surface of the tibia. Based on our previous study [Al Nazer et al. (2008) J Biomech. 41:1036–1043], an alternative numerical approach to analyse in vivo strains based on the flexible multibody simulation approach was proposed. The purpose of this study was to extend the idea of using the flexible multibody approach in the analysis of bone strains during physical activity through integrating the magnetic resonance imaging (MRI) technique within the framework. In order to investigate the reliability and validity of the proposed approach, a three-dimensional full body musculoskeletal model with a flexible tibia was used as a demonstration example. The model was used in a forward dynamics simulation in order to predict the tibial strains during walking on a level exercise. The flexible tibial model was developed using the actual geometry of human tibia, which was obtained from three-dimensional reconstruction of MRI. Motion capture data obtained from walking at constant velocity were used to drive the model during the inverse dynamics simulation in order to teach the muscles to reproduce the motion in the forward dynamics simulation. Based on the agreement between the literature-based in vivo strain measurements and the simulated strain results, it can be concluded that the flexible multibody approach enables reasonable predictions of bone strain in response to dynamic loading. The information obtained from the present approach can be useful in clinical applications including devising exercises to prevent bone fragility or to accelerate fracture healing.     

Revisiting the links between bone remodelling and osteocytes: insights from across phyla

We question two major tenets of bone biology: that the primary role of remodelling is to remove damage in the bone (so‐called damage‐driven remodelling) and that osteocytes are the only strain‐sensing orchestrators of this process. These concepts are distilled largely from research on model mammal species, but in fact, there are a number of features of various bones, from mammalian and non‐mammalian species, that do not accord with these ‘rules’. Here, we assemble a variety of examples, ranging from species that lack osteocytes but that still seem capable of remodelling their bones, to species with osteocytic bones that do not remodel, and to instances of inter‐species, inter‐bone and/or intra‐bone variation in bone remodelling that show that this purported repair process is not always where the ‘rules’ tell us it should be. This collection of points argues that our understanding of the advantages, roles and primary drivers of remodelling are inadequate and biased to quite a small phylogenetic cross section of the species that possess bone. We suggest a variety of new directions for bone research that would provide us with a better understanding of bone remodelling, tying together the interests of comparative biologists, palaeontologists and medical researchers.     

Numerical simulation of load-induced bone structural remodelling using stress-limit criterion

A simple and efficient numerical method for predicting the remodelling of adaptive materials and structures under applied loading was presented and implemented within a finite element framework. The model uses the trajectorial architecture theory of optimisation to predict the remodelling of material microstructure and structural organisation under mechanical loading. We used the proposed model to calculate the density distribution of proximal femur in the frontal plane. The loading considered was the hip joint contact forces and muscular forces at the attachment sites of the muscles to the bone. These forces were estimated from a separate finite element calculation using a heterogeneous three-dimensional model of the proximal femur. The density distributions obtained by this procedure has a qualitative similarity with in vivo observations. Solutions displayed the characteristic high-density channels that are evident in the Dual X-ray Absorptiometry scan. There is also evidence of the intramedullary canal, as well as low-density regions in the femoral neck. Several parametric studies were carried out to highlight the advantages of the proposed method, which includes fast convergence and low-computational cost. The potential applications of the proposed method in predicting bone structural remodelling in cancer are also briefly discussed.     

Allosteric control model of bone remodelling containing periodical modes

To help to understand the modelling process that occurs when a scaffold is implanted it is vital to understand the rather complex bone remodelling process prevalent in native bone. We have formulated a mathematical model that predicts osteoactivity both in scaffolds, as well as in bone in vivo and could set a basis for the more detailed allosteric models. The model is extended towards a bio-cybernetic vision of basic multicellular unit (BMU) action, when some of the regulation loops have been modified to reflect the allosteric control mechanisms, developed by Michaels-Menten, Hill, Koshland-Nemethy-Filmer, Monod-Wyman-Changeux. By implementation of this approach a four-dimensional system was obtained that shows steady cyclic behaviour using a wide range of constants with clear biological meaning. We have observed that a local steady state appears as a limiting cycle in multi-dimensional phase space and this is discussed in this paper. Physiological interpretation of this limiting four-dimension cycle possibly related to a conservative-like value has been proposed. Analysis and simulation of the model has shown an analogy between this conservative value, as a kind of substrate-energy regenerative potential of the bone remodelling system with a molecular nature, and to the classical physical value--energy. This dynamic recovery potential is directed against both mechanical and biomechanical damage to the bone. Furthermore, the current model has credibility when compared to the normal bone remodelling process. In the framework of widely recognised Hill mechanisms of allosteric regulation the cyclic attractor, described formerly for a pure cellular model, prevails for different forms of feedback control. This result indicates the viability of the proposed existence of a conservative value (analogous to energy) that characterises the recovery potential of the bone remodelling cycle. Linear stability analysis has been performed in order to determine the robustness of the basic state, however, additional work is required to study a wider range of constants.     

Bone remodelling algorithms incorporating both strain and microdamage stimuli

Biomechanical theories to predict bone remodelling have used either mechanical strain or microdamage as the stimulus driving cellular responses. Even though experimental data have implicated both stimuli in bone cell regulation, a mechano-regulatory system incorporating both stimuli has not yet been proposed. In this paper, we test the hypothesis that bone remodelling may be regulated by signals due to both strain and microdamage. Four mechano-regulation algorithms are studied where the stimulus is: strain, damage, combined strain/damage, and either strain or damage with damage-adaptive remodelling prioritised when damage is above a critical level. Each algorithm is implemented with both bone lining cell (surface) sensors and osteocyte cell (internal) sensors. Each algorithm is applied to prediction of a bone multicellular unit (BMU) remodelling on the surface of a bone trabecula. It is predicted that a regulatory system capable of responding to changes in either strain or microdamage but which prioritises removal of damaged bone when damage is above a critical level, is the only one that provides a plausible prediction of BMU behaviour. A mechanism for this may be that, below a certain damage threshold, osteocyte processes can sense changes in strain and fluid flow but above the threshold damage interferes with the signalling mechanism, or causes osteocyte apoptosis so that a remodelling response occurs to remove the dead osteocytes.