An adaptation model for trabecular bone at different mechanical levels

Background  

Bone has the ability to adapt to mechanical usage or other biophysical stimuli in terms of its mass and architecture, indicating that a certain mechanism exists for monitoring mechanical usage and controlling the bone's adaptation behaviors. There are four zones describing different bone adaptation behaviors: the disuse, adaptation, overload, and pathologic overload zones. In different zones, the changes of bone mass, as calculated by the difference between the amount of bone formed and what is resorbed, should be different.     

A theoretical framework for strain-related trabecular bone maintenance and adaptation

It is assumed that density and morphology of trabecular bone is partially controlled by mechanical forces. How these effects are expressed in the local metabolic functions of osteoclast resorption and osteoblast formation is not known. In order to investigate possible mechano-biological pathways for these mechanisms we have proposed a mathematical theory (Nature 405 (2000) 704). This theory is based on hypothetical osteocyte stimulation of osteoblast bone formation, as an effect of elevated strain in the bone matrix, and a role for microcracks and disuse in promoting osteoclast resorption. Applied in a 2-D Finite Element Analysis model, the theory explained the formation of trabecular patterns. In this article we present a 3-D FEA model based on the same theory and investigated its potential morphological predictability of metabolic reactions to mechanical loads. The computations simulated the development of trabecular morphological details during growth, relative to measurements in growing pigs, reasonably realistic. They confirmed that the proposed mechanisms also inherently lead to optimal stress transfer. Alternative loading directions produced new trabecular orientations. Reduction of load reduced trabecular thickness, connectivity and mass in the simulation, as is seen in disuse osteoporosis. Simulating the effects of estrogen deficiency through increased osteoclast resorption frequencies produced osteoporotic morphologies as well, as seen in post-menopausal osteoporosis. We conclude that the theory provides a suitable computational framework to investigate hypothetical relationships between bone loading and metabolic expressions.     

Finite element analysis of a three-dimensional open-celled model for trabecular bone

Based on a regular array of cubic unit cells, each containing a body-centered spherical void, we created an idealized three-dimensional model for both subchondral trabecular bone and a class of porous foams. By considering only face-to-face stacking of unit cells, the inherent symmetry was such that, except at the surface, the displacements and stresses within any one unit cell were representative of the entire porous structure. Using prescribed displacements the model was loaded in both uniaxial compressive strain and uniaxial shear strain. Based on the response to these loads, we found the tensor of elastic constants for an equivalent homogeneous elastic solid with cubic symmetry. We then compared the predicted modulus with our experimental values for bovine trabecular bone and literature values for an open-celled latex rubber foam.     

A morphological model of vertebral trabecular bone

In their micro-structures, typical natural cellular materials such as vertebral trabecular bone have a network of doubly tapered struts, thickening near the strut joints. However, past analytical models for vertebral trabecular bone do not take account of the effect of strut taper on the mechanical properties.This paper presents an analytical cell model comprised of doubly tapered struts to predict the global mechanical properties of vertebral trabecular bone. The predicted results for male, female, and both sexes fit the experimental data well. By considering several strut taper geometries, it is shown that the horizontal Young's modulus and the horizontal uniaxial collapse stress are, in some cases, approximately 1.8- and 2.2-fold higher, respectively, than those of the uniform strut model. This finding illustrates the importance of increased trabecular thickening near the strut joints (i) for improving the accuracy of calculating the mechanical properties and (ii) for the effective treatment of aged bone using drug therapy. It also highlights the need to combine trabecular architecture measurements with information about the morphology near the strut joints.     

Quantification of trabecular bone microdamage using the virtual internal bond model and the individual trabeculae segmentation technique

Trabecular bone microdamage significantly influences the skeletal integrity and bone remodelling process. In this paper a novel constitutive model, called the virtual internal bond model (VIB), was adopted for simulating the damage behaviour of bone tissue. A unique 3D image analysis technique, named individual trabeculae segmentation, was used to analyse the effects of microarchitectures on the damage behaviours of trabecular bone. We demonstrated that the process of initiation and accumulation of microdamage in trabecular bone samples can be captured by the VIB-embedded finite-element method simulation without a separate fracture criterion. Our simulation results showed that the microdamage can occur at as early as about 0.2–0.4% apparent strain, and a large volume of microdamage was accumulated around the apparent yield strain. In addition we found that the plate-like trabeculae, especially the longitudinal ones, take crucial roles in the microdamage behaviours of trabecular bone.     

The canine tail artery as a model for cerebral aneurysm studies

The occluded canine tail artery, which comes off in the same plane as the aortoiliac junction, has been used as a flow model for cerebral aneurysms. These experiments were designed to determine if it is a realistic distensible model of human intracranial aneurysms. Distensibility studies were done on the aorta, and the iliac and tail arteries of four dogs. From these pressure-volume studies, tension-strain curves, elastances, and collagen slack were obtained. The tail artery is stiffer longitudinally and more distensible circumferentially than the other vessels. The iliac arteries and the aorta are not significantly different. The elastance of elastin and collagen is lower in the tail artery, and the collagen is more wavy circumferentially. Longitudinally, the collagen slack is least for the tail artery, and the elastance of elastin is not different in all three vessels. The number of elastin layers in the iliac and tail arteries seen in cross section is not significantly different, but the aorta is different from both these vessels. In another four dogs the aorta proximal to the trifurcation was cannulated and infused with saline to increase pressure. India ink marks were put on the surface to measure changes in length. Photographs were taken at intervals of 10 mmHg(1 mmHg = 133.3 Pa). This was done with the vessels tethered and untethered in the body and then taken out and studied with the same method in vitro. Arteries tethered in the body expanded circumferentially more than longitudinally. The tail artery becomes less distensible if untethered in the body and therefore acts more like an aneurysm.(ABSTRACT TRUNCATED AT 250 WORDS)     

A structural model for the mechanical behavior of trabecular bone

A quantitative model is developed for trabecular bone by approximating the trabecular geometry with a hypothetical network of compact bone. For the region immediately beneath the articular cartilage in the distal end of the femur, finite element analyses were performed with a high speed computer, assuming a physiological static load. The results indicate that bending and buckling of trabeculae are considerable in any elastic deformation of the bone; that fatigue fracture in some fraction of suitably oriented trabeculae is inevitable in normal ambulation; and that the stiffness varies considerably with lateral position across the subchondral plate. The latter depends totally on trabecular arrangement and may play a role in joint function and degeneration. The adjustments necessary to bring the gross stiffness into agreement with experiment imply that the intertrabecular soft tissues are of no consequence to the mechanical properties and that the compact bone of which trabeculae are made is probably not as stiff as cortical bone.     

A coupled mechano-biochemical model for bone adaptation

Bone remodelling is a fundamental biological process that controls bone microrepair, adaptation to environmental loads and calcium regulation among other important processes. It is not surprising that bone remodelling has been subject of intensive both experimental and theoretical research. In particular, many mathematical models have been developed in the last decades focusing in particular aspects of this complicated phenomenon where mechanics, biochemistry and cell processes strongly interact. In this paper, we present a new model that combines most of these essential aspects in bone remodelling with especial focus on the effect of the mechanical environment into the biochemical control of bone adaptation mainly associated to the well known RANKL-RANK-OPG pathway. The predicted results show a good correspondence with experimental and clinical findings. For example, our results indicate that trabecular bone is more severely affected both in disuse and disease than cortical bone what has been observed in osteoporotic bones. In future, the methodology proposed would help to new therapeutic strategies following the evolution of bone tissue distribution in osteoporotic patients.     

A generic 3-dimensional system to mimic trabecular bone surface adaptation

The paper presents the structure optimizing system based on surface remodelling. The grounds for algorithm formulation are given by the bone remodelling phenomenon leading to optimization of trabecular network in the bone. The assumptions, algorithms and limitations of the own mesh generator Cosmoprojector are described. Unlike other approaches, the system is able to mimic real bone evolution including tissue consolidation and separation. The article presents a closed system consisting of finite element mesh generation, decision criteria for structure adaptation and finite element analysis in parallel environment. It also provides some computation results obtained by using specially designed software.     

A generic 3-dimensional system to mimic trabecular bone surface adaptation

The paper presents the structure optimizing system based on surface remodelling. The grounds for algorithm formulation are given by the bone remodelling phenomenon leading to optimization of trabecular network in the bone. The assumptions, algorithms and limitations of the own mesh generator Cosmoprojector are described. Unlike other approaches, the system is able to mimic real bone evolution including tissue consolidation and separation. The article presents a closed system consisting of finite element mesh generation, decision criteria for structure adaptation and finite element analysis in parallel environment. It also provides some computation results obtained by using specially designed software.     

A nonlocal constitutive model for trabecular bone softening in compression

Using the three-dimensional morphological data provided by computed tomography, finite element (FE) models can be generated and used to compute the stiffness and strength of whole bones. Three-dimensional constitutive laws capturing the main features of bone mechanical behavior can be developed and implemented into FE software to enable simulations on complex bone structures. For this purpose, a constitutive law is proposed, which captures the compressive behavior of trabecular bone as a porous material with accumulation of irreversible strain and loss of stiffness beyond its yield point and softening beyond its ultimate point. To account for these features, a constitutive law based on damage coupled with hardening anisotropic elastoplasticity is formulated using density and fabric-based tensors. To prevent mesh dependence of the solution, a nonlocal averaging technique is adopted. The law has been implemented into a FE software and some simple simulations are first presented to illustrate its behavior. Finally, examples dealing with compression of vertebral bodies clearly show the impact of softening on the localization of the inelastic process.     

The limitations of canine trabecular bone as a model for human: a biomechanical study

Distal canine femurs were sectioned into 8 mm cubic specimens. Orthogonal compression tests were performed to preyield in two or three directions and to failure in a third. Apparent density and ash weight density were measured for a subset of specimens. The results were compared to the human distal femur results of Ciarelli et al. (Transactions of the 32nd Annual Meeting of the Orthopaedic Research Society, Vol. 11, p. 42, 1986). Quantitative similarities existed in the fraction of components comprising the trabecular tissue of the two species. Qualitative similarities were seen in the positional and anisotropic variation of the mechanical properties, and also in the form and strength of the relationships between the mean modulus and bone density, ultimate stress and density, and ultimate stress and modulus. However, significantly different regression equations resulted for the mean modulus-density, and ultimate stress modulus relationships, indicating that for the same density, canine trabecular bone displays a lower modulus than human, and may achieve greater compressive strains before failure.     

Artificial composite bone as a model of human trabecular bone: the implant-bone interface

The use of artificial bones in implant testing has become popular due to their low variability and ready availability. However, friction coefficients, which are critical to load transfer in uncemented implants, have rarely been compared between human and artificial bone, particularly for wet and dry conditions. In this study, the static and dynamic friction coefficients for four commercially used titanium surfaces (polished, Al(2)O(3) blasted, plasma sprayed, beaded) acting on the trabecular component of artificial bones (Sawbones) were compared to those for human trabecular bone. Artificial bones were tested in dry and wet conditions and normal interface stress was varied (0.25, 0.5, 1.0MPa). Friction coefficients were mostly lower for artificial bones than real bone. In particular, static friction coefficients for the dry polished surface were 20% of those for real bone and 42-61% for the dry beaded surface, with statistical significance (alpha<0.05). Less marked differences were observed for dynamic friction coefficients. Significant but non-systematic effects of normal stress or wet/dry condition on friction coefficients were observed within each surface type. These results indicate that the use of artificial bone models for pre-clinical implant testing that rely on interface load transfer with trabecular bone for mechanical integrity can be particularly sensitive to surface finish and lubrication conditions.     

A phenomenological model for predicting fatigue life in bovine trabecular bone

Cyclic loading of bone during daily activities can lead to fatigue degradation and increased risk of fracture in both the young and elderly population. Damage processes under cyclic loading in trabecular bone result in the reduction of the elastic modulus and accumulation of residual strain. These effects increase with increasing stress levels, leading to a progressive reduction in fatigue life. The present work analyzes the effect of stress and strain variation on the above damage processes in bovine trabecular bone, and develops a phenomenological model relating fatigue life to the imposed stress level. The elastic modulus reduction of the bone specimens was observed to depend on the maximum compressive strain, while the rate of residual strain accumulation was a function of the stress level. A model was developed for the upper and lower bounds of bone elastic modulus reduction with increasing number of cycles, at each stress range. The experimental observations were described well by the model. The model predicted the bounds of the fatigue life with change in fatigue stress. The decrease in the fatigue life with increasing stress was related to corresponding increases in the residual strain accumulation rates at the elevated stress levels. The model shows the validity of fatigue predictions from relatively few cyclic experiments, by combining trends observed in the monotonic and the cyclic tests. The model also presents a relatively simple procedure for predicting the endurance limit for bovine trabecular bone specimens.     

A novel mathematical model of bone remodelling cycles for trabecular bone at the cellular level

After an initial phase of growth and development, bone undergoes a continuous cycle of repair, renewal and optimisation by a process called remodelling. This paper describes a novel mathematical model of the trabecular bone remodelling cycle. It is essentially formulated to simulate a remodelling event at a fixed position in the bone, integrating bone removal by osteoclasts and formation by osteoblasts. The model is developed to construct the variation in bone thickness at a particular point during the remodelling event, derived from standard bone histomorphometric analyses. The novelties of the approach are the adoption of a predator-prey model to describe the dynamic interaction between osteoclasts and osteoblasts, using a genetic algorithm-based solution; quantitative reconstruction of the bone remodelling cycle; and the introduction of a feedback mechanism in the bone formation activity to co-regulate bone thickness. The application of the model is first demonstrated by using experimental data recorded for normal (healthy) bone remodelling to predict the temporal variation in the number of osteoblasts and osteoclasts. The simulated histomorphometric data and remodelling cycle characteristics compare well with the specified input data. Sensitivity studies then reveal how variations in the model's parameters affect its output; it is hoped that these parameters can be linked to specific biochemical factors in the future. Two sample pathological conditions, hypothyroidism and primary hyperparathyroidism, are examined to demonstrate how the model could be applied more broadly, and, for the first time, the osteoblast and osteoclast populations are predicted for these conditions. Further data are required to fully validate the model's predictive capacity, but this work shows it has potential, especially in the modelling of pathological conditions and the optimisation of the treatment of those conditions.     

Experimental and finite element analysis of the mouse caudal vertebrae loading model: prediction of cortical and trabecular bone adaptation

In this study, we attempt to predict cortical and trabecular bone adaptation in the mouse caudal vertebrae loading model using knowledge of bone’s local mechanical environment at the onset of loading. In a previous study, we demonstrated appreciable 25.9 and 11% increases in both trabecular and cortical bone volume density, respectively, when subjecting the fifth caudal vertebrae (C5) of C57BL/6 (B6) mice to an acute loading regime (amplitude of 8N, 3000 cycles, 10 Hz, 3 times a week for 4 weeks). We have also established a validated finite element (FE) model of the C5 vertebra using micro-computed tomography (micro-CT), which characterizes, in 3D, the micro-mechanical strains present in both cortical and trabecular compartments due to the applied loads. To investigate the relationship between load-induced bone adaptation and mechanical strains in-vivo and in-silico data sets were compared. Using data from the previous cross-sectional study, we divided cortical and trabecular compartments into 15 subregions and determined, for each region, a bone formation parameter ΔBV/BS (a cross-sectional measure of the bone volume added to cortical and trabecular surfaces following the described loading regime). Linear regression was then used to correlate mean regional values of ΔBV/BS with mean values of mechanical strains derived from the FE models which were similarly regionalized. The mechanical parameters investigated were strain energy density (SED), the orthogonal strains (e _x , e _y , e _z ) and the three shear strains (e _xy , e _yz , e _zx ). For cortical regions, regression analysis showed SED to correlate extremely well with ΔBV/BS (R ² =?0.82) and e _z (R ²?=?0.89). Furthermore, SED was found to predict expansion of the cortical shell correlating significantly with the regional percentage increases in cortical tissue volume (R ² = 0.92), cortical marrow volume (R ² =?0.91) and cortical thickness (R ² = 0.56). For trabecular regions, FE parameters were found not to correlate with load-induced trabecular bone morphology. These results indicate that load-induced cortical morphology can be predicted from population data, whereas the prediction of trabecular morphology requires subject-specific micro- architecture.     

A 3D damage model for trabecular bone based on fabric tensors

Motivated by the role of damage in normal and pathological conditions of trabecular bone, a novel 3D constitutive law was developed that describes anisotropic elasticity and the rate-independent degradation in mechanical properties resulting from the growth of cracks or voids in the trabecular tissue. The theoretical model was formulated within the framework of continuum damage mechanics and based on two fabric tensors characterizing the local trabecular morphology. Experimental validation of the model was achieved by uniaxial and torsional testing of waisted bovine trabecular bone specimens. Strong correlations were found between cumulated permanent strain, reduction in elastic moduli and nonlinear postyield stress which support the hypothesis that these variables reflect the same underlying damage process.     

Describing force-induced bone growth and adaptation by a mathematical model

This work proposes a mathematical model that qualitative describes the process of mechanically force-induced bone growth and adaptation. The mathematical model includes osteocytes as the key interfacing layer connecting tissue, cellular and molecular signaling levels. Specifically, in the presence of an increase in the mechanical stimuli, osteocytes respond by mechano-transduction releasing the local factors nitric oxide (NO) and prostaglandin E(2) (PGE(2)). These local factors act as the signaling recruitment signals for bone cells progenitors and influence the coupling activity among osteoblasts and osteoclasts during the process of bone remodeling. The model is in agreement with qualitative observations found in the literature concerning the process of bone adaptation and the cellular interactions during a local bone remodeling cycle induced by mechanical stimulation.