Simulation of vertebral trabecular bone loss using voxel finite element analysis

Trabecular bone loss in human vertebral bone is characterised by thinning and eventual perforation of the horizontal trabeculae. Concurrently, vertical trabeculae are completely lost with no histological evidence of significant thinning. Such bone loss results in deterioration in apparent modulus and strength of the trabecular core. In this study, a voxel-based finite element program was used to model bone loss in three specimens of human vertebral trabecular bone. Three sets of analyses were completed. In Set 1, strain adaptive resorption was modelled, whereby elements which were subject to the lowest mechanical stimulus (principal strain) were removed. In Set 2, both strain adaptive and microdamage mechanisms of bone resorption were included. Perforation of vertical trabeculae occurred due to microdamage resorption of elements with strains that exceeded a damage threshold. This resulted in collapse of the trabecular network under compression loading for two of the specimens tested. In Set 3, the damage threshold strain was gradually increased as bone loss progressed, resulting in reduced levels of microdamage resorption. This mechanism resulted in trabecular architectures in which vertical trabeculae had been perforated and which exhibited similar apparent modulus properties compared to experimental values reported in the literature. Our results indicate that strain adaptive remodelling alone does not explain the deterioration in mechanical properties that have been observed experimentally. Our results also support the hypothesis that horizontal trabeculae are lost principally by strain adaptive resorption, while vertical trabeculae may be lost due to perforation from microdamage resorption followed by rapid strain adaptive resorption of the remaining unloaded trabeculae.     

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.     

Stress-concentrating effect of resorption lacunae in trabecular bone

Analyses of the distributions of stress and strain within individual bone trabeculae have not yet been reported. In this study, four trabeculae were imaged and finite elements models were generated in an attempt to quantify the variability of stress/strain in real trabeculae. In three of these trabeculae, cavities were identified with depths comparable to values reported for resorption lacunae ( approximately 50 microm)-although we cannot be certain, it is most probable that they are indeed resorption lacunae. A tensile load was applied to each trabeculum to simulate physiological loading and to ensure that bending was minimized. The force carried by each trabecula was calculated from this value using the average cross sectional area of each trabecula. The analyses predict that very high stresses (>100 MPa) existed within bone trabecular tissue. Stress and strain distributions were highly heterogeneous in all cases, more so in trabeculae with the presumptive resorption lacunae where at least 30% of the tissue had a strain greater than 4000 micoepsilon in all cases. Stresses were elevated at the pit of the lacunae, and peak stress concentrations were located in the longitudinal direction ahead of the lacunae. Given these high strains, we suggest that microdamage is inevitable around resorption lacunae in trabecular bone, and may cause the bone multicellular unit to proceed to resorb a packet of bone in the trabeculum rather than just resorb whatever localized area was initially targeted.     

The effects of tensile-compressive loading mode and microarchitecture on microdamage in human vertebral cancellous bone

The amount of microdamage in bone tissue impairs mechanical performance and may act as a stimulus for bone remodeling. Here we determine how loading mode (tension vs. compression) and microstructure (trabecular microarchitecture, local trabecular thickness, and presence of resorption cavities) influence the number and volume of microdamage sites generated in cancellous bone following a single overload. Twenty paired cylindrical specimens of human vertebral cancellous bone from 10 donors (47–78 years) were mechanically loaded to apparent yield in either compression or tension, and imaged in three dimensions for microarchitecture and microdamage (voxel size 0.7×0.7×5.0 μm³). We found that the overall proportion of damaged tissue was greater (p=0.01) for apparent tension loading (3.9±2.4%, mean±SD) than for apparent compression loading (1.9±1.3%). Individual microdamage sites generated in tension were larger in volume (p<0.001) but not more numerous (p=0.64) than sites in compression. For both loading modes, the proportion of damaged tissue varied more across donors than with bone volume fraction, traditional measures of microarchitecture (trabecular thickness, trabecular separation, etc.), apparent Young?s modulus, or strength. Microdamage tended to occur in regions of greater trabecular thickness but not near observable resorption cavities. Taken together, these findings indicate that, regardless of loading mode, accumulation of microdamage in cancellous bone after monotonic loading to yield is influenced by donor characteristics other than traditional measures of microarchitecture, suggesting a possible role for tissue material properties.     

Theoretical analysis of alendronate and risedronate effects on canine vertebral remodeling and microdamage

Bisphosphonates suppress bone remodeling activity, increase bone volume, and significantly reduce fracture risk in individuals with osteoporosis and other metabolic bone diseases. The objectives of the current study were to develop a mathematical model that simulates control and 1 year experimental results following bisphosphonate treatment (alendronate or risedronate) in the canine fourth lumbar vertebral body, validate the model by comparing simulation predictions to 3 year experimental results, and then use the model to predict potential long term effects of bisphosphonates on remodeling and microdamage accumulation. To investigate the effects of bisphosphonates on bone volume and microdamage, a mechanistic biological model was modified from previous versions to simulate remodeling in a representative volume of vertebral trabecular bone in dogs treated with various doses of alendronate or risedronate, including doses equivalent to those used for treatment of post-menopausal osteoporosis in humans. Bisphosphonates were assumed to affect remodeling by suppressing basic multicellular unit activation and reducing resorption area. Model simulation results for trabecular bone volume fraction, microdamage, and activation frequency following 1 year of bisphosphonate treatment are consistent with experimental measurements. The model predicts that trabecular bone volume initially increases rapidly with 1 year of bisphosphonate treatment, and continues to slowly rise between 1 and 3 years of treatment. The model also predicts that microdamage initially increases rapidly, 0.5–1.5-fold for alendronate or risedronate during the first year of treatment, and reaches its maximum value by 2.5 years before trending downward for all dosages. The model developed in this study suggests that increasing bone volume fraction with long term bisphosphonate treatment may sufficiently reduce strain and damage formation rate so that microdamage does not accumulate above that which is initiated in the first two years of treatment.     

Low bone mass may not be the only cause of skeletal fragility in osteoporosis

Skeletal fragility in postmenopausal osteoporosis is not due solely to reduction in bone mass. This fact explains some of the overlap in bone mineral measurements observed between patients who are fracturing and age- and sex-matched normals who are not. Changes in skeletal architecture and bone remodeling occur with age which can account for some of the fragility. These changes are exaggerated in patients with postmenopausal osteoporosis who are suffering spine fractures. Three abnormalities have been described by histomorphometric methods which can account for skeletal fragility out of proportion to the degree of bone loss. They are: (i) loss of trabecular connectivity such that vertical weight-bearing bars lose their cross-attachments with each other, thus becoming susceptible to buckling; (ii) inefficient and prolonged microdamage repair due to periods of pause in the formation phase of remodeling; and (iii) accumulation of unrepaired microdamage in unremodeled bone tissue in the central part of trabeculae due to reduced osteon wall thickness coupled with maintenance of trabecular thickness. Recognition of these abnormalities should broaden our approach to the study of skeletal fragility in the syndrome of postmenopausal osteoporosis.     

The effect of resorption cavities on bone stiffness is site dependent

Resorption cavities formed during the bone remodelling cycle change the structure and thus the mechanical properties of trabecular bone. We tested the hypotheses that bone stiffness loss due to resorption cavities depends on anatomical location, and that for identical eroded bone volumes, cavities would cause more stiffness loss than homogeneous erosion. For this purpose, we used beam–shell finite element models. This new approach was validated against voxel-based FE models. We found an excellent agreement for the elastic stiffness behaviour of individual trabeculae in axial compression (R² = 1.00) and in bending (R²>0.98), as well as for entire trabecular bone samples to which resorption cavities were digitally added (R² = 0.96, RMSE = 5.2%). After validation, this new method was used to model discrete cavities, with dimensions taken from a statistical distribution, on a dataset of 120 trabecular bone samples from three anatomical sites (4th lumbar vertebra, femoral head, iliac crest). Resorption cavities led to significant reductions in bone stiffness. The largest stiffness loss was found for samples from the 4th lumbar vertebra, the lowest for femoral head samples. For all anatomical sites, resorption cavities caused significantly more stiffness loss than homogeneous erosion did. This novel technique can be used further to evaluate the impact of resorption cavities, which are known to change in several metabolic bone diseases and due to treatment, on bone competence.     

Bone turnover and trabecular plate survival after artificial menopause.

Considerable bone loss often occurs after menopause, particularly if menopause is induced by surgery. For perhaps two years bone formation fails to keep pace with the rapid acceleration of bone resorption that occurs after sex hormone withdrawal. The threat that this poses to the integrity of the skeleton is not clear. Because ethical constraints limit histological studies in normal women existing normal data and statistical modelling techniques were used to explore the dynamics of iliac trabecular bone after menopause. Trabeculae are breached during remodelling when the osteoclasts resorb to a depth equal to the trabecular thickness. Since holes in trabecular plates cannot normally be bridged such defects are probably permanent. Men lose 7% of their vertebral trabecular bone every 10 years; deeper than average resorption of trabeculae at the thin end of the normal range would account for it. The dramatic losses of trabecular bone that are seen in some postmenopausal women, however, require a period of imbalance between bone formation and bone resorption since this leads rapidly to generalised thinning. The statistical model suggested that an imbalance lasting only two years may account for eventual losses of up to half of the iliac trabecular bone. Further understanding is needed of what determines the amount of bone lost in the immediate postmenopause, which varies considerably among women. A simple mean is needed of identifying women who will lose bone most rapidly at the menopause. This must be suitable for use in general practice because these women should probably be offered long term hormone replacement treatment within a few months of the last menstruation.     

Modeling of dynamic fracture and damage in two-dimensional trabecular bone microstructures using the cohesive finite element method

Trabecular bone fracture is closely related to the trabecular architecture, microdamage accumulation, and bone tissue properties. Micro-finite-element models have been used to investigate the elastic and yield properties of trabecular bone but have only seen limited application in modeling the microstructure dependent fracture of trabecular bone. In this research, dynamic fracture in two-dimensional (2D) micrographs of ovine (sheep) trabecular bone is modeled using the cohesive finite element method. For this purpose, the bone tissue is modeled as an orthotropic material with the cohesive parameters calculated from the experimental fracture properties of the human cortical bone. Crack propagation analyses are carried out in two different 2D orthogonal sections cut from a three-dimensional 8 mm diameter cylindrical trabecular bone sample. The two sections differ in microstructural features such as area fraction (ratio of the 2D space occupied by bone tissue to the total 2D space), mean trabecula thickness, and connectivity. Analyses focus on understanding the effect of the rate of loading as well as on how the rate variation interacts with the microstructural features to cause anisotropy in microdamage accumulation and in the fracture resistance. Results are analyzed in terms of the dependence of fracture energy dissipation on the microstructural features as well as in terms of the changes in damage and stresses associated with the bone architecture variation. Besides the obvious dependence of the fracture behavior on the rate of loading, it is found that the microstructure strongly influences the fracture properties. The orthogonal section with lesser area fraction, low connectivity, and higher mean trabecula thickness is more resistant to fracture than the section with high area fraction, high connectivity, and lower mean trabecula thickness. In addition, it is found that the trabecular architecture leads to inhomogeneous distribution of damage, irrespective of the symmetry in the applied loading with the fracture of the entire bone section rapidly progressing to bone fragmentation once the accumulated damage in any trabeculae reaches a critical limit.     

Age-related changes in human trabecular bone: Relationship between microstructural stress and strain and damage morphology

Accumulation of microdamage in aging and disease can cause skeletal fragility and is one of several factors contributing to osteoporotic fractures. To better understand the role of microdamage in fragility fracture, the mechanisms of bone failure must be elucidated on a tissue-level scale where interactions between bone matrix properties, the local biomechanical environment, and bone architecture are concurrently examined for their contributions to microdamage formation. A technique combining histological damage assessment of individual trabeculae with linear finite element solutions of trabecular von Mises and principal stress and strain was used to compare the damage initiation threshold between pre-menopausal (32-37 years, n=3 donors) and post-menopausal (71-80 years, n=3 donors) femoral cadaveric bone. Strong associations between damage morphology and stress and strain parameters were observed in both groups, and an age-related decrease in undamaged trabecular von Mises stress was detected. In trabeculae from younger donors, the 95% CI for von Mises stress on undamaged regions ranged from 50.7-67.9MPa, whereas in trabeculae from older donors, stresses were significantly lower (38.7-50.2, p<0.01). Local microarchitectural analysis indicated that thinner, rod-like trabeculae oriented along the loading axis are more susceptible to severe microdamage formation in older individuals, while only rod-like architecture was associated with severe damage in younger individuals. This study therefore provides insight into how damage initiation and morphology relate to local trabecular microstructure and the associated stresses and strains under loading. Furthermore, by comparison of samples from pre- and post-menopausal women, the results suggest that trabeculae from younger individuals can sustain higher stresses prior to microdamage initiation.     

Micro-CT analysis of cancellous bone fragments from the distal radius fracture zone in osteoporosis]

The incidence and prevalence of osteoporosis must be considered to continue to increase significantly due to the expected demographic development and environmental changes. In the diagnosis and staging of osteoporosis the three-dimensional bone structure should be as important as the bone mass or the mineral content of the bone. In this study, microfragments were taken from distal radius fracture zones and investigated in Micro-CT scans. Patients in which osteodensitometry of the lumbal spine had revealed osteoporosis in were found to have significantly reduced bone mass, bone density and trabecular thickness. Trabecular fractures which were found in non-osteoporotic patients even in robust trabeculae were detected by the two-dimensional analysis in thin locations and arborizations. Despite some trabeculae turned out to be very small the differences in the histomorphometry and the quality of trabecular fractures in osteoporotic as well as non-osteoporotic patients could be visualized very good in the Micro-CT analysis.     

Mechanical properties of single bovine trabeculae are unaffected by strain rate

For a better understanding of traumatic bone fractures, knowledge of the mechanical behavior of bone at high strain rates is required. Importantly, it needs to be clarified how quasistatic mechanical testing experiments relate to real bone fracture. This merits investigating the mechanical behavior of bone with an increase in strain rate. Various studies examined how cortical and trabecular bone behave at varying strain rates, but no one has yet looked at this question using individual trabeculae. In this study, three-point bending tests were carried out on bovine single trabeculae excised from a proximal femur to test the trabecular material's strain rate sensitivity. An experimental setup was designed, capable of measuring local strains at the surface of such small specimens, using digital image correlation. Microdamage was detected using the bone whitening effect. Samples were tested through two orders of magnitude, at strain rates varying between 0.01 and 3.39 s(-1). No linear relationship was observed between the strain rate and the Young's modulus (1.13-16.46 GPa), the amount of microdamage, the maximum tensile strain at failure (14.22-61.65%) and at microdamage initiation (1.95-12.29%). The results obtained in this study conflict with previous studies reporting various trends for macroscopic cortical and trabecular bone samples with an increase in strain rate. This discrepancy might be explained by the bone type, the small sample geometry, the limited range of strain rates tested here, the type of loading and the method of microdamage detection. Based on the results of this study, the strain rate can be ignored when modeling trabecular bone.     

Adaptive diffuse domain approach for calculating mechanically induced deformation of trabecular bone

Remodelling of trabecular bone is essentially affected by the mechanical load of the trabeculae. Mathematical modelling and simulation of the remodelling process have to include time-consuming calculations of the displacement field within the complex trabecular structure under loading. We present an adaptive diffuse domain approach for calculating the elastic bone deformation based on micro computer tomogram data of real trabecular bone structures and compared it with a conventional voxel-based finite element method. In addition to allowing for higher computational efficiency, the adaptive approach is characterised by a very smooth representation of the bone surface, which suggests that this approach would be suitable as a basis for future simulations of bone resorption and formation processes within the trabecular structure.     

An algorithm for bone mechanoresponsiveness: implementation to study the effect of patient-specific cell mechanosensitivity on trabecular bone loss

The rate of bone loss is subject to considerable variation between individuals. With the 'mechanostat' model of Frost, genetic variations in bone mechanoresponsiveness are modelled by different mechanostat 'setpoints'--which may also change with age or disease. In this paper, the following setpoints are used: epsilonmin (strain below which resorption is triggered); epsilonmax (strain above which deposition occurs); omegacrit (microdamage-level above which damage-stimulated resorption occurs). To simulate decreased mechanosensitivity, epsilonmax is increased. Analyses carried out on a simplified model of a trabecula show that epsilonmax is a critical parameter: if it is higher in an individual (genetics) or increases (with age) the mass deficit each remodelling cycle increases. Furthermore, there is a value of epsilonmax above which trabecular perforation occurs, leading to rapid loss of bone mass. Maintaining bone cell mechanosensitivity could therefore be a therapeutic target for the prevention of osteoporosis.     

Heterogeneous glycation of cancellous bone and its association with bone quality and fragility

Non-enzymatic glycation (NEG) and enzymatic biochemical processes create crosslinks that modify the extracellular matrix (ECM) and affect the turnover of bone tissue. Because NEG affects turnover and turnover at the local level affects microarchitecture and formation and removal of microdamage, we hypothesized that NEG in cancellous bone is heterogeneous and accounts partly for the contribution of microarchitecture and microdamage on bone fragility. Human trabecular bone cores from 23 donors were subjected to compression tests. Mechanically tested cores as well as an additional 19 cores were stained with lead-uranyl acetate and imaged to determine microarchitecture and measure microdamage. Post-yield mechanical properties were measured and damaged trabeculae were extracted from a subset of specimens and characterized for the morphology of induced microdamage. Tested specimens and extracted trabeculae were quantified for enzymatic and non-enzymatic crosslink content using a colorimetric assay and Ultra-high Performance Liquid Chromatography (UPLC). Results show that an increase in enzymatic crosslinks was beneficial for bone where they were associated with increased toughness and decreased microdamage. Conversely, bone with increased NEG required less strain to reach failure and were less tough. NEG heterogeneously modified trabecular microarchitecture where high amounts of NEG crosslinks were found in trabecular rods and with the mechanically deleterious form of microdamage (linear microcracks). The extent of NEG in tibial cancellous bone was the dominant predictor of bone fragility and was associated with changes in microarchitecture and microdamage.     

Modeling modulus reduction in bovine trabecular bone damaged in compression.

Loading bone beyond its yield point creates microdamage, leading to reduction in stiffness. Previously, we related microdamage accumulation to changes in mechanical properties. Here, we develop a model that predicts stiffness loss based on the presence of microdamage. Modeling is done at three levels: (1) a single trabecula, (2) a cellular solid consisting of intact, damaged, and fractured trabeculae, and (3) a specimen with a localized damage band. Predictions of a reduced modulus agree well with experimental measured modulus reductions of post-yield compression of bovine trabecular bone. The predicted reduced modulus is relatively insensitive to changes in the input parameters.     

Age-related differences in the morphology of microdamage propagation in trabecular bone

Microdamage density has been shown to increase with age in trabecular bone and is associated with decreased fracture toughness. Numerous studies of crack propagation in cortical bone have been conducted, but data in trabecular bone is lacking. In this study, propagation of severe, linear, and diffuse damage was examined in trabecular bone cores from the femoral head of younger (61.3±3.1 years) and older (75.0±3.9 years) men and women. Using a two-step mechanical testing protocol, damage was first initiated with static uniaxial compression to 0.8% strain then propagated at a normalized stress level of 0.005 to a strain endpoint of 0.8%. Coupling mechanical testing with a dual-fluorescent staining technique, the number and length/area of propagating cracks were quantified. It was found that the number of cycles to the test endpoint was substantially decreased in older compared to younger samples (younger: 77,372±15,984 cycles; older: 34,944±11,964 cycles, p=0.06). This corresponded with a greater number of severely damaged trabeculae expanding in area during the fatigue test in the older group. In the younger group, diffusely damaged trabeculae had a greater damage area, which illustrates an efficient energy dissipation mechanism. These results suggest that age-related differences in fatigue life of human trabecular bone may be due to differences in propagated microdamage morphology.     

Stéroides sexuels et ostéoporose chez l'homme

Several epidemiological studies have shown that about 25 per cent of hip fractures and 20 per cent of symptomatic vertebral fractures occur in men. The lifetime risk of hip fracture was estimated to be about 6 to 8 per cent and the risk of any osteoporotic fracture was estimated to be about 18 per cent in 50-year-old white men. In about 60% of cases in men, bone loss is secondary to several pathological conditions, such as long-term steroid therapy, severe hypogonadism, smoking or alcohol abuse or gastrointestinal disorders. In 40% of cases, osteoporosis is primary or idiopathic in men between the ages of 40 and 60 years. Genetic factors, a defect of boneforming cells or abnormal serum levels of bioavailable sex steroids could explain bone loss and fragility fractures in these men. It has been shown that hypogonadism is associated with a marked increase in bone remodelling and particularly in bone resorption with a dramatic loss in trabecular bone. It is now known that testosterone is partly transformed into estradiol by aromatase. Testosterone may therefore act on bone in two ways: it directly stimulates bone formation and estradiol regulates bone remodelling and inhibits bone resorption. Finally, in men over the age of 60 without hypogonadism, it has been shown that bone mineral density and fracture risk were better correlated with serum levels of bioavailable estradiol and Sex Hormone Binding Globulin than with serum testosterone levels.     

A comparison of the fatigue behavior of human trabecular and cortical bone tissue

The fatigue properties of trabecular bone tissue (single trabeculae) and similarly sized cortical bone specimens from human tibia were experimentally determined on a microstructural level using four-point bending cyclic tests, and they were compared based on modulus, mineral density, and microstructural characteristics. The results showed that trabecular specimens had significantly lower moduli and lower fatigue strength than cortical specimens, despite their higher mineral density values. Fracture surface and microdamage analyses illustrated different fracture and damage patterns between trabecular and cortical bone tissue, depending upon their microstructural characteristics. Based on the results from mechanical tests and qualitative observations, a possible mechanical role of the cement lines in trabecular tissue microfracture was suggested.     

Trabecular bone microdamage and microstructural stresses under uniaxial compression

The balance between local remodeling and accumulation of trabecular bone microdamage is believed to play an important role in the maintenance of skeletal integrity. However, the local mechanical parameters associated with microdamage initiation are not well understood. Using histological damage labeling, micro-CT imaging, and image-based finite element analysis, regions of trabecular bone microdamage were detected and registered to estimated microstructural von Mises effective stresses and strains, maximum principal stresses and strains, and strain energy density (SED). Bovine tibial trabecular bone cores underwent a stepwise uniaxial compression routine in which specimens were micro-CT imaged following each compression step. The results indicate that the mode of trabecular failure observed by micro-CT imaging agreed well with the polarity and distribution of stresses within an individual trabecula. Analysis of on-axis subsections within specimens provided significant positive relationships between microdamage and each estimated tissue stress, strain and SED parameter. In a more localized analysis, individual microdamaged and undamaged trabeculae were extracted from specimens loaded within the elastic region and to the apparent yield point. As expected, damaged trabeculae in both groups possessed significantly higher local stresses and strains than undamaged trabeculae. The results also indicated that microdamage initiation occurred prior to apparent yield at local principal stresses in the range of 88-121 MPa for compression and 35-43 MPa for tension and local principal strains of 0.46-0.63% in compression and 0.18-0.24% in tension. These data provide an important step towards understanding factors contributing to microdamage initiation and establishing local failure criteria for normal and diseased trabecular bone.