A novel in vivo mouse model for mechanically stimulated bone adaptation--a combined experimental and computational validation study

To facilitate the investigation of bone formation, in vivo, in response to mechanical loading a caudal vertebra axial compression device (CVAD) has been developed to deliver precise mechanical loads to the fifth caudal vertebra (C5) of the C57BL/6 female mouse. A combined experimental and computational approach was used to quantify the micro-mechanical strain induced in trabecular and cortical components following static and dynamic loading using the CVAD. Cortical bone strains were recorded using micro-strain gages. Finite element (FE) models based on micro-computed tomography were constructed for all C5 vertebrae. Both theoretical and experimental cortical strains correlated extremely well (R(2)>0.96) for a Young's modulus of 14.8 GPa, thus validating the FE model. In this study, we have successfully applied mechanical loads to the C5 murine vertebrae, demonstrating the potential of this model to be used for in vivo loading studies aimed at stimulating both trabecular and cortical bone adaptation.     

Comparison of the elastic and yield properties of human femoral trabecular and cortical bone tissue

The ability to determine trabecular bone tissue elastic and failure properties has biological and clinical importance. To date, trabecular tissue yield strains remain unknown due to experimental difficulties, and elastic moduli studies have reported controversial results. We hypothesized that the elastic and tensile and compressive yield properties of trabecular tissue are similar to those of cortical tissue. Effective tissue modulus and yield strains were calibrated for cadaveric human femoral neck specimens taken from 11 donors, using a combination of apparent-level mechanical testing and specimen-specific, high-resolution, nonlinear finite element modeling. The trabecular tissue properties were then compared to measured elastic modulus and tensile yield strain of human femoral diaphyseal cortical bone specimens obtained from a similar cohort of 34 donors. Cortical tissue properties were obtained by statistically eliminating the effects of vascular porosity. Results indicated that mean elastic modulus was 10% lower (p<0.05) for the trabecular tissue (18.0+/-2.8 GPa) than for the cortical tissue (19.9+/-1.8 GPa), and the 0.2% offset tensile yield strain was 15% lower for the trabecular tissue (0.62+/-0.04% vs. 0.73+/-0.05%, p<0.001). The tensile-compressive yield strength asymmetry for the trabecular tissue, 0.62 on average, was similar to values reported in the literature for cortical bone. We conclude that while the elastic modulus and yield strains for trabecular tissue are just slightly lower than those of cortical tissue, because of the cumulative effect of these differences, tissue strength is about 25% greater for cortical bone.     

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.     

Dependence of yield strain of human trabecular bone on anatomic site

Understanding the dependence of human trabecular bone strength behavior on anatomic site provides insight into structure-function relationships and is essential to the increased success of site-specific finite element models of whole bones. To investigate the hypothesis that the yield strains of human trabecular bone depend on anatomic site, the uniaxial tensile and compressive yield properties were compared for cylindrical specimens from the vertebra (n=61), proximal tibia (n=31), femoral greater trochanter (n=23), and femoral neck (n=27) taken from 61 donors (67+/-15years). Test protocols were used that minimized end artifacts and loaded specimens along the main trabecular orientation. Yield strains by site (mean+/-S.D.) ranged from 0.70+/-0.05% for the trochanter to 0.85+/-0.10% for the femoral neck in compression, from 0.61+/-0.05% for the trochanter to 0.70+/-0.05% for the vertebra in tension, and were always higher in compression than tension (p<0.001). The compressive yield strain was higher for the femoral neck than for all other sites (p<0.001), as was the tensile yield strain for the vertebra (p<0.007). Analysis of covariance, with apparent density as the covariate, indicated that inter-site differences existed in yield stress even after adjusting statistically for density (p<0.035). Coefficients of variation in yield strain within each site ranged from only 5-12%, consistent with the strong linear correlations (r(2)=0.94-0.98) found between yield stress and modulus. These results establish that the yield strains of human trabecular bone can differ across sites, but that yield strain may be considered uniform within a given site despite substantial variation in elastic modulus and yield stress.     

Internal strains and stresses measured in cortical bone via high-energy X-ray diffraction

High-energy synchrotron X-ray diffraction was used to study internal stresses in bone under in situ compressive loading. A transverse cross-section of a 12-14 year old beagle fibula was studied with 80.7 keV radiation, and the transmission geometry was used to quantify internal strains and corresponding stresses in the mineral phase, carbonated hydroxyapatite. The diffraction patterns agreed with tabulated patterns, and the distribution of diffracted intensity around 00.2/00.4 and 22.2 diffraction rings was consistent with the imperfect 00.1 fiber texture expected along the axis of a long bone. Residual compressive stress along the bone's longitudinal axis was observed in the specimen prior to testing: for 22.2 this stress equaled -95 MPa and for 00.2/00.4 was between -160 and -240 MPa. Diffraction patterns were collected for applied compressive stresses up to -110 MPa, and, up to about -100 MPa, internal stresses rose proportionally with applied stress but at a higher rate, corresponding to stress concentration in the mineral of 2.8 times the stress applied. The widths of the 00.2 and 00.4 diffraction peaks indicated that crystallite size perpendicular to the 00.1 planes increased from t=41 nm before stress was applied to t=44 nm at -118 MPa applied stress and that rms strain epsilon(rms) rose from 2200 muepsilon before loading to 4600 muepsilon at the maximum applied stress. Small angle X-ray scattering of the unloaded sample, recorded after deformation was complete, showed a collagen D-period of 66.4 nm (along the bone axis).     

Constrained tibial vibration in mice: a method for studying the effects of vibrational loading of bone

Vibrational loading can stimulate the formation of new trabecular bone or maintain bone mass. Studies investigating vibrational loading have often used whole-body vibration (WBV) as their loading method. However, WBV has limitations in small animal studies because transmissibility of vibration is dependent on posture. In this study, we propose constrained tibial vibration (CTV) as an experimental method for vibrational loading of mice under controlled conditions. In CTV, the lower leg of an anesthetized mouse is subjected to vertical vibrational loading while supporting a mass. The setup approximates a one degree-of-freedom vibrational system. Accelerometers were used to measure transmissibility of vibration through the lower leg in CTV at frequencies from 20 Hz to 150 Hz. First, the frequency response of transmissibility was quantified in vivo, and dissections were performed to remove one component of the mouse leg (the knee joint, foot, or soft tissue) to investigate the contribution of each component to the frequency response of the intact leg. Next, a finite element (FE) model of a mouse tibia-fibula was used to estimate the deformation of the bone during CTV. Finally, strain gages were used to determine the dependence of bone strain on loading frequency. The in vivo mouse leg in the CTV system had a resonant frequency of 60 Hz for +/-0.5 G vibration (1.0 G peak to peak). Removing the foot caused the natural frequency of the system to shift from 60 Hz to 70 Hz, removing the soft tissue caused no change in natural frequency, and removing the knee changed the natural frequency from 60 Hz to 90 Hz. By using the FE model, maximum tensile and compressive strains during CTV were estimated to be on the cranial-medial and caudolateral surfaces of the tibia, respectively, and the peak transmissibility and peak cortical strain occurred at the same frequency. Strain gage data confirmed the relationship between peak transmissibility and peak bone strain indicated by the FE model, and showed that the maximum cyclic tibial strain during CTV of the intact leg was 330+/-82microepsilon and occurred at 60-70 Hz. This study presents a comprehensive mechanical analysis of CTV, a loading method for studying vibrational loading under controlled conditions. This model will be used in future in vivo studies and will potentially become an important tool for understanding the response of bone to vibrational loading.     

Mechanical loading of mouse caudal vertebrae increases trabecular and cortical bone mass-dependence on dose and genotype

Most in vivo studies addressing the skeletal responses of mice to mechanical loading have targeted cortical bone. To investigate trabecular bone responses also we have developed a caudal vertebral axial compression device (CVAD) that transmits mechanical loads to compress the fifth caudal vertebra via stainless steel pins inserted into the forth and sixth caudal vertebral bodies. Here, we used the CVAD in C57BL/6 (B6) and C3H/Hej (C3H) female mice (15 weeks of age) to investigate whether the effect of regular bouts of mechanical stimulation on bone modeling and bone mass was dependent on dose and genotype. A combined micro-computed tomographic and dynamic histomorphometric analysis was carried out at the end of a 4-week loading regimen (3,000 cycles, 10 Hz, 3 x week) for load amplitudes of 0N, 2N, 4N and 8N. Significant increases in trabecular bone mass of 9 and 21% for loads of 4N and 8N, respectively, were observed in B6 mice. A significant increase of 10% in trabecular bone mass occurred for a load of 8N in the C3H strain. For other loads, no significant increases were detected. Both mouse strains exhibited substantial increases in trabecular bone formation rates for all loads, B6: 111% (2N), 86% (4N), 164% (8N), C3H: 41% (2N), 38% (4N), 141% (8N). Significant decreases in osteoclast number of 146 and 93% for a load of 8N were detected in B6 and C3H mice, respectively. These findings demonstrate that the effect of loading on the structural and functional parameters of bone is dose and genotype dependent. The caudal vertebral loading model established here is proposed for further studies addressing the molecular processes involved in the skeletal responses to mechanical stimuli.     

Subject-specific finite element model of the pelvis: development, validation and sensitivity studies

A better understanding of the three-dimensional mechanics of the pelvis, at the patient-specific level, may lead to improved treatment modalities. Although finite element (FE) models of the pelvis have been developed, validation by direct comparison with subject-specific strains has not been performed, and previous models used simplifying assumptions regarding geometry and material properties. The objectives of this study were to develop and validate a realistic FE model of the pelvis using subject-specific estimates of bone geometry, location-dependent cortical thickness and trabecular bone elastic modulus, and to assess the sensitivity of FE strain predictions to assumptions regarding cortical bone thickness as well as bone and cartilage material properties. A FE model of a cadaveric pelvis was created using subject-specific computed tomography image data. Acetabular loading was applied to the same pelvis using a prosthetic femoral stem in a fashion that could be easily duplicated in the computational model. Cortical bone strains were monitored with rosette strain gauges in ten locations on the left hemipelvis. FE strain predictions were compared directly with experimental results for validation. Overall, baseline FE predictions were strongly correlated with experimental results (r2=0.824), with a best-fit line that was not statistically different than the line y=x (experimental strains = FE predicted strains). Changes to cortical bone thickness and elastic modulus had the largest effect on cortical bone strains. The FE model was less sensitive to changes in all other parameters. The methods developed and validated in this study will be useful for creating and analyzing patient-specific FE models to better understand the biomechanics of the pelvis.     

Load transfer analysis of the distal radius from in-vivo high-resolution CT-imaging.

Prevention of osteoporotic bone fractures requires accurate diagnostic methods to detect the increase in bone fragility at an early stage of osteoporosis. However, today's bone fracture risk prediction, primarily based on bone density measurement, is not sufficiently precise. There is increasing evidence that, in addition to bone density, also the bone microarchitecture and its mechanical loading conditions are important factors determining the fracture risk. Recently, it has been shown that new high-resolution imaging techniques in combination with new computer modeling techniques based on the finite-element (FE) method can account for these additional factors. These techniques might provide information that is more relevant for the prediction of bone fracture risk. So far, however, these new imaged-based FE techniques have not been feasible in-vivo. The objectives of this study were to quantify the load transfer through the trabecular network in a distal radius using a computer model based on in-vivo high-resolution images and to determine if common regions of fractures can be explained as a result of high tissue loading in these regions. The left distal radius and the two adjacent carpal bones of a healthy volunteer were imaged using a high-resolution three-dimensional CT system providing an isotropic resolution of 165 microm. The bone representation was converted into a FE-model that was used to calculate stresses and strains in the trabecular network. The two carpal bones were loaded using different load ratios (for each load case 1000 N in total) representing impact forces on the hand either in near-neutral position or ulnar/radial deviation. The load transfer through the trabecular network of the radius was characterized by the tissue strain energy density (SED) distribution for all load cases. It was found that the distribution of the tissue loading depends on the ratio of the forces acting on the carpal bones. For all load cases the higher SED values (on average: 0.02 +/- 0.08 (S.D.) N mm(-2)) are found in a 10 mm region adjacent to the articular surface which corresponds well with the region where Colles- or Chauffeur-fractures occur. We expect that, eventually, this new approach can lead to a better prediction of the fracture risk than methods based on bone density alone since it accounts for the bone microstructure as well as its loading conditions.     

Simplified boundary conditions alter cortical-trabecular load sharing at the distal radius; A multiscale finite element analysis

High-resolution peripheral quantitative computed tomography (HR-pQCT) derived micro-finite element (FE) modeling is used to evaluate mechanical behavior at the distal radius microstructure. However, these analyses typically simulate non-physiologic simplified platen-compression boundary conditions on a small section of the distal radius. Cortical and trabecular regions contribute uniquely to distal radius mechanical behavior, and various factors affect these regions distinctly. Generalized strength predictions from standardized platen-compression analyses may not adequately capture region specific responses in bone. Our goal was to compare load sharing within the cortical-trabecular compartments between the standardized platen-compression BC simulations, and physiologic BC simulations using a validated multiscale approach. Clinical- and high-resolution images were acquired from nine cadaveric forearm specimens using an HR-pQCT scanner. Multiscale FE models simulating physiologic BCs, and micro-FE only models simulating platen-compression BCs were created for each specimen. Cortical and trabecular loads (N) along the length of the distal radius micro-FE section were compared between BCs using correlations. Principal strain distributions were also compared quantitatively. Cortical and trabecular loads from the platen-compression BC simulations were strongly correlated to the physiologic BC simulations. However, a 30% difference in cortical loads distally, and a 53% difference in trabecular loads proximally was observed under platen BC simulations. Also, distribution of principal strains was clearly different. Our data indicated that platen-compression BC simulations alter cortical-trabecular load sharing. Therefore, results from these analyses should be interpreted in the appropriate mechanical context for clinical evaluations of normal and pathologic mechanical behavior at the distal radius.     

High-resolution finite element models with tissue strength asymmetry accurately predict failure of trabecular bone

The ability to predict trabecular failure using microstructure-based computational models would greatly facilitate study of trabecular structure–function relations, multiaxial strength, and tissue remodeling. We hypothesized that high-resolution finite element models of trabecular bone that include cortical-like strength asymmetry at the tissue level, could predict apparent level failure of trabecular bone for multiple loading modes. A bilinear constitutive model with asymmetric tissue yield strains in tension and compression was applied to simulate failure in high-resolution finite element models of seven bovine tibial specimens. Tissue modulus was reduced by 95% when tissue principal strains exceeded the tissue yield strains. Linear models were first calibrated for effective tissue modulus against specimen-specific experimental measures of apparent modulus, producing effective tissue moduli of (mean±S.D.) 18.7±3.4 GPa. Next, a parameter study was performed on a single specimen to estimate the tissue level tensile and compressive yield strains. These values, 0.60% strain in tension and 1.01% strain in compression, were then used in non-linear analyses of all seven specimens to predict failure for apparent tensile, compressive, and shear loading. When compared to apparent yield properties previously measured for the same type of bone, the model predictions of both the stresses and strains at failure were not statistically different for any loading case (p>0.15). Use of symmetric tissue strengths could not match the experimental data. These findings establish that, once effective tissue modulus is calibrated and uniform but asymmetric tissue failure strains are used, the resulting models can capture the apparent strength behavior to an outstanding level of accuracy. As such, these computational models have reached a level of fidelity that qualifies them as surrogates for destructive mechanical testing of real specimens.     

An accurate estimation of bone density improves the accuracy of subject-specific finite element models

An experimental-numerical study was performed to investigate the relationships between computed tomography (CT)-density and ash density, and between ash density and apparent density for bone tissue, to evaluate their influence on the accuracy of subject-specific FE models of human bones. Sixty cylindrical bone specimens were examined. CT-densities were computed from CT images while apparent and ash densities were measured experimentally. The CT/ash-density and ash/apparent-density relationships were calculated. Finite element models of eight human femurs were generated considering these relationships to assess their effect on strain prediction accuracy. CT and ash density were linearly correlated (R(2)=0.997) over the whole density range but not equivalent (intercep t <0, slope >1). A constant ash/apparent-density ratio (0.598+/-0.004) was found for cortical bone. A lower ratio, with a larger dispersion, was found for trabecular bone (0.459+/-0.100), but it became less dispersed, and equal to that of cortical tissue, when testing smaller trabecular specimens (0.598+/-0.036). This suggests that an experimental error occurred in apparent-density measurements for large trabecular specimens and a constant ratio can be assumed valid for the whole density range. Introducing the obtained relationships in the FE modelling procedure improved strain prediction accuracy (R(2)=0.95, RMSE=7%). The results suggest that: (i) a correction of the densitometric calibration should be used when evaluating bone ash-density from clinical CT scans, to avoid ash-density underestimation and overestimation for low- and high-density bone tissue, respectively; (ii) the ash/apparent-density ratio can be assumed constant in human femurs and (iii) the correction improves significantly the model accuracy and should be considered in subject-specific bone modelling.     

Mechanisms of initial endplate failure in the human vertebral body

Endplate failure occurs frequently in osteoporotic vertebral fractures and may be related to the development of high tensile strain. To determine whether the highest tensile strains in the vertebra occur in the endplates, and whether such high tensile strains are associated with the material behavior of the intervertebral disc, we used micro-CT-based finite element analysis to assess tissue-level strains in 22 elderly human vertebrae (81.5±9.6 years) that were compressed through simulated intervertebral discs. In each vertebra, we compared the highest tensile and compressive strains across the different compartments: endplates, cortical shell, and trabecular bone. The influence of Poisson-type expansion of the disc on the results was determined by compressing the vertebrae a second time in which we suppressed the Poisson expansion. We found that the highest tensile strains occurred within the endplates whereas the highest compressive strains occurred within the trabecular bone. The ratio of strain to assumed tissue-level yield strain was the highest for the endplates, indicating that the endplates had the greatest risk of initial failure. Suppressing the Poisson expansion of the disc decreased the amount of highly tensile-strained tissue in the endplates by 79.4±11.3%. These results indicate that the endplates are at the greatest risk of initial failure due to the development of high tensile strains, and that such high tensile strains are associated with the Poisson expansion of the disc. We conclude that initial failure of the vertebra is associated with high tensile strains in the endplates, which in turn are influenced by the material behavior of the disc.     

The modified super-ellipsoid yield criterion for human trabecular bone

Despite the importance of multiaxial failure of trabecular bone in many biomechanical applications, to date no complete multiaxial failure criterion for human trabecular bone has been developed. By using experimentally validated nonlinear high-resolution, micromechanical finite-element models as a surrogate for multiaxial loading experiments, we determined the three-dimensional normal strain yield surface and all combinations of the two-dimensional normal-shear strain yield envelope. High-resolution finite-element models of three human femoral neck trabecular bone specimens obtained through microcomputed tomography were used. In total, 889 multiaxial-loading cases were analyzed, requiring over 41,000 CPU hours on parallel supercomputers. Our results indicated that the multiaxial yield behavior of trabecular bone in strain space was homogeneous across the specimens and nearly isotropic. Analysis of stress-strain curves along each axis in the 3-D normal strain space indicated uncoupled yield behavior whereas substantial coupling was seen for normal-shear loading. A modified super-ellipsoid surface with only four parameters fit the normal strain yield data very well with an arithmetic error +/-SD less than -0.04 +/- 5.1%. Furthermore, the principal strains associated with normal-shear loading showed excellent agreement with the yield surface obtained for normal strain loading (arithmetic error +/- SD < 2.5 +/- 6.5%). We conclude that the four-parameter "Modified Super-Ellipsoid" yield surface presented here describes the multiaxial failure behavior of human femoral neck trabecular bone very well.     

Damage in trabecular bone at small strains

Evidence that damage decreases bone quality, increases fracture susceptibility, and serves as a remodeling stimulus motivates further study of what loading magnitudes induce damage in trabecular bone. In particular, whether damage occurs at the smaller strains characteristic of habitual, as opposed to traumatic, loading is not known. The overall goal of this study was to characterize damage accumulation in trabecular bone at small strains (0.20 - 0.45% strain). A continuum damage mechanics approach was taken whereby damage was quantified by changes in modulus and residual strain. Human vertebral specimens (n = 7) were tested in compression using a multi-cycle load - unload protocol in which the maximum applied strain for each cycle, epsilonmax, was increased incrementally from epsilonmax = 0.20% on the first loading cycle to epsilonmax = 0.45% on the last cycle. Modulus and residual strain were measured for each cycle. Both changes in modulus and residual strains commenced at small strains, beginning as early as 0.24 and 0.20% strain, respectively. Strong correlations between changes in modulus and residual strains were observed (r = 0.51 - 0.98). Fully nonlinear, high-resolution finite element analyses indicated that even at small apparent strains, tissue-level strains were sufficiently high to cause local yielding. These results demonstrate that damage in trabecular bone occurs at apparent strains less than half the apparent compressive yield strain reported previously for human vertebral trabecular bone. Further, these findings imply that, as a consequence of the highly porous trabecular structure, tissue yielding can initiate at very low apparent strains and that this local failure has detectable and negative consequences on the apparent mechanical properties of trabecular bone.     

Rib fractures under anterior–posterior dynamic loads: Experimental and finite-element study

The purpose of this study was to investigate whether using a finite-element (FE) mesh composed entirely of hexahedral elements to model cortical and trabecular bone (all-hex model) would provide more accurate simulations than those with variable thickness shell elements for cortical bone and hexahedral elements for trabecular bone (hex–shell model) in the modeling human ribs. First, quasi-static non-injurious and dynamic injurious experiments were performed using the second, fourth, and tenth human thoracic ribs to record the structural behavior and fracture tolerance of individual ribs under anterior–posterior bending loads. Then, all-hex and hex–shell FE models for the three ribs were developed using an octree-based and multi-block hex meshing approach, respectively. Material properties of cortical bone were optimized using dynamic experimental data and the hex–shell model of the fourth rib and trabecular bone properties were taken from the literature. Overall, the reaction force–displacement relationship predicted by both all-hex and hex–shell models with nodes in the offset middle-cortical surfaces compared well with those measured experimentally for all the three ribs. With the exception of fracture locations, the predictions from all-hex and offset hex–shell models of the second and fourth ribs agreed better with experimental data than those from the tenth rib models in terms of reaction force at fracture (difference <15.4%), ultimate failure displacement and time (difference <7.3%), and cortical bone strains. The hex–shell models with shell nodes in outer cortical surfaces increased static reaction forces up to 16.6%, compared to offset hex–shell models. These results indicated that both all-hex and hex–shell modeling strategies were applicable for simulating rib responses and bone fractures for the loading conditions considered, but coarse hex–shell models with constant or variable shell thickness were more computationally efficient and therefore preferred.     

Fracture prediction for the proximal femur using finite element models: Part II--Nonlinear analysis.

In Part I we reported the results of linear finite element models of the proximal femur generated using geometric and constitutive data collected with quantitative computed tomography. These models demonstrated excellent agreement with in vitro studies when used to predict ultimate failure loads. In Part II, we report our extension of those finite element models to include nonlinear behavior of the trabecular and cortical bone. A highly nonlinear material law, originally designed for representing concrete, was used for trabecular bone, while a bilinear material law was used for cortical bone. We found excellent agreement between the model predictions and in vitro fracture data for both the onset of bone yielding and bone fracture. For bone yielding, the model predictions were within 2 percent for a load which simulated one-legged stance and 1 percent for a load which simulated a fall. For bone fracture, the model predictions were within 1 percent and 17 percent, respectively. The models also demonstrated different fracture mechanisms for the two different loading configurations. For one-legged stance, failure within the primary compressive trabeculae at the subcapital region occurred first, leading to load transfer and, ultimately, failure of the surrounding cortical shell. However, for a fall, failure of the cortical and trabecular bone occurred simultaneously within the intertrochanteric region. These results support our previous findings that the strength of the subcapital region is primarily due to trabecular bone whereas the strength of the intertrochanteric region is primarily due to cortical bone.     

3D characterization of bone strains in the rat tibia loading model

Bone strain is considered one of the factors inducing bone tissue response to loading. Nevertheless, where animal studies can provide detailed data on bone response, they only offer limited information on experimental bone strains. Including micro-CT-based finite element (micro FE) models in the analysis represents a potent methodology for quantifying strains in bone. Therefore, the main objective of this study was to develop and validate specimen-specific micro FE models for the assessment of bone strains in the rat tibia compression model. Eight rat limbs were subjected to axial compression loading; strain at the medio-proximal site of the tibiae was measured by means of strain gauges. Specimen-specific micro FE models were created and analyzed. Repeated measurements on each limb indicated that the effect of limb positioning was small (COV?= 6.45 ± 2.27 %). Instead, the difference in the measured strains between the animals was high (54.2%). The computational strains calculated at the strain gauge site highly correlated to the measured strains (R ²?=?0.95). Maximum peak strains calculated at exactly 25% of the tibia length for all specimens were equal to 435.11 ± 77.88 microstrains (COV?=?17.19%). In conclusion, we showed that strain gauge measurements are very sensitive to the exact strain gauge location on the bone; hence, the use of strain gauge data only is not recommended for studies that address at identifying reliable relationships between tissue response and local strains. Instead, specimen-specific micro FE models of rat tibiae provide accurate estimates of tissue-level strains.     

Structural and Mechanical Improvements to Bone Are Strain Dependent with Axial Compression of the Tibia in Female C57BL/6 Mice

Strain-induced adaption of bone has been well-studied in an axial loading model of the mouse tibia. However, most outcomes of these studies are restricted to changes in bone architecture and do not explore the mechanical implications of those changes. Herein, we studied both the mechanical and morphological adaptions of bone to three strain levels using a targeted tibial loading mouse model. We hypothesized that loading would increase bone architecture and improve cortical mechanical properties in a dose-dependent fashion. The right tibiae of female C57BL/6 mice (8 week old) were compressively loaded for 2 weeks to a maximum compressive force of 8.8N, 10.6N, or 12.4N (generating periosteal strains on the anteromedial region of the mid-diaphysis of 1700 με, 2050 με, or 2400 με as determined by a strain calibration), while the left limb served as an non-loaded control. Following loading, ex vivo analyses of bone architecture and cortical mechanical integrity were assessed by micro-computed tomography and 4-point bending. Results indicated that loading improved bone architecture in a dose-dependent manner and improved mechanical outcomes at 2050 με. Loading to 2050 με resulted in a strong and compelling formation response in both cortical and cancellous regions. In addition, both structural and tissue level strength and energy dissipation were positively impacted in the diaphysis. Loading to the highest strain level also resulted in rapid and robust formation of bone in both cortical and cancellous regions. However, these improvements came at the cost of a woven bone response in half of the animals. Loading to the lowest strain level had little effect on bone architecture and failed to impact structural- or tissue-level mechanical properties. Potential systemic effects were identified for trabecular bone volume fraction, and in the pre-yield region of the force-displacement and stress-strain curves. Future studies will focus on a moderate load level which was largely beneficial in terms of cortical/cancellous structure and cortical mechanical function.     

Effects of suppression of bone turnover on cortical and trabecular load sharing in the canine vertebral body

The relative biomechanical effects of antiresorptive treatment on cortical thickness vs. trabecular bone microarchitecture in the spine are not well understood. To address this, T-10 vertebral bodies were analyzed from skeletally mature female beagle dogs that had been treated with oral saline (n=8 control) or a high dose of oral risedronate (0.5 mg/kg/day, n=9 RIS-suppressed) for 1 year. Two linearly elastic finite element models (36-μm voxel size) were generated for each vertebral body—a whole-vertebra model and a trabecular-compartment model—and subjected to uniform compressive loading. Tissue-level material properties were kept constant to isolate the effects of changes in microstructure alone. Suppression of bone turnover resulted in increased stiffness of the whole vertebra (20.9%, p=0.02) and the trabecular compartment (26.0%, p=0.01), while the computed stiffness of the cortical shell (difference between whole-vertebra and trabecular-compartment stiffnesses, 11.7%, p=0.15) was statistically unaltered. Regression analyses indicated subtle but significant changes in the relative structural roles of the cortical shell and the trabecular compartment. Despite higher average cortical shell thickness in RIS-suppressed vertebrae (23.1%, p=0.002), the maximum load taken by the shell for a given value of shell mass fraction was lower (p=0.005) for the RIS-suppressed group. Taken together, our results suggest that—in this canine model—the overall changes in the compressive stiffness of the vertebral body due to suppression of bone turnover were attributable more to the changes in the trabecular compartment than in the cortical shell. Such biomechanical studies provide an unique insight into higher-scale effects such as the biomechanical responses of the whole vertebra.