Mechanics of linear microcracking in trabecular bone

Microcracking in trabecular bone is responsible both for the mechanical degradation and remodeling of the trabecular bone tissue. Recent results on trabecular bone mechanics have demonstrated that bone tissue microarchitecture, tissue elastic heterogeneity and tissue-level mechanical anisotropy all should be considered to obtain detailed information on the mechanical stress state. The present study investigated the influence of tissue microarchitecture, tissue heterogeneity in elasticity and material separation properties and tissue-level anisotropy on the microcrack formation process. Microscale bone models were executed with the extended finite element method. It was demonstrated that anisotropy and heterogeneity of the bone tissue contribute significantly to bone tissue toughness and the resistance of trabecular bone to microcrack formation. The compressive strain to microcrack initiation was computed to increase by a factor of four from an assumed homogeneous isotropic tissue to an assumed anisotropic heterogenous tissue.     

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.     

Finite element analysis of peri-implant bone volume affected by stresses around Morse taper implants: effects of implant positioning to the bone crest

Abstract

Objectives: The purpose of the present study was to evaluate the distribution and magnitude of stresses through the bone tissue surrounding Morse taper dental implants at different positioning relative to the bone crest. Materials and Methods: A mandibular bone model was obtained from a computed tomography scan. A three-dimensional (3D) model of Morse taper implant-abutment systems placed at the bone crest (equicrestal) and 2?mm bellow the bone crest (subcrestal) were assessed by finite element analysis (FEA). FEA was carried out on axial and oblique (45°) loading at 150 N relatively to the central axis of the implant. The von Mises stresses were analysed considering magnitude and volume of affected peri-implant bone. Results: On vertical loading, maximum von Mises stresses were recorded at 6-7?MPa for trabecular bone while values ranging from 73 up to 118?MPa were recorded for cortical bone. On oblique loading at the equiquestral or subcrestal positioning, the maximum von Mises stresses ranged from 15 to 21?MPa for trabecular bone while values at 150?MPa were recorded for the cortical bone. On vertical loading, >99.9vol.% cortical bone volume was subjected to a maximum of 2?MPa while von Mises stress values at 15?MPa were recorded for trabecular bone. On oblique loading, >99.9vol.% trabecular bone volume was subjected to maximum stress values at 5?MPa, while von Mises stress values at 35?MPa were recorded for >99.4vol.% cortical bone. Conclusions: Bone volume-based stress analysis revealed that most of the bone volume (>99% by vol) was subjected to significantly lower stress values around Morse taper implants placed at equicrestal or subcrestal positioning. Such analysis is commentary to the ordinary biomechanical assessment of dental implants concerning the stress distribution through peri-implant sites.     

Mineral heterogeneity affects predictions of intratrabecular stress and strain

Knowledge of the influence of mineral variations (i.e., mineral heterogeneity) on biomechanical bone behavior at the trabecular level is limited. The aim of this study is to investigate how this material property affects the intratrabecular distributions of stress and strain in human adult trabecular bone. Two different sets of finite element (FE) models of trabecular samples were constructed; tissue stiffness was either scaled to the local degree of mineralization of bone as measured with microCT (heterogeneous) or tissue stiffness was assumed to be homogeneous. The influence of intratrabecular mineral heterogeneity was analyzed by comparing both models. Interesting effects were seen regarding intratrabecular stress and strain distributions. In the homogeneous model, the highest stresses were found at the surface with a significant decrease towards the core. Higher superficial stresses could indicate a higher predicted fracture risk in the trabeculae. In the heterogeneous model this pattern was different. A significant increase in stress with increasing distance from the trabecular surface was found followed by a significant decrease towards the core. This suggests trabecular bending during a compression. In both models a decrease in strain values from surface to core was predicted, which is consistent with trabecular bending. When mineral heterogeneity was taken into account, the predicted intratrabecular patterns of stress and strain are more consistent with the expected biomechanical behavior as based on mineral variations in trabeculae. Our findings indicate that mineral heterogeneity should not be neglected when performing biomechanical studies on topics such as the (long-term or dose dependent) effects of antiresorptive treatments.     

Mechanisms of uniformity of yield strains for trabecular bone

Variations in yield strains for trabecular bone within a specific anatomic site are only a small fraction of the substantial variations that exist for elastic modulus and strength, and yet the source of this uniformity is not known. Our goal was to investigate the underlying mechanisms by using high-resolution, materially nonlinear finite element models of 12 human femoral neck trabecular bone specimens. The finite element models, used to obtain apparent yield strains in both tension and compression, assumed that the tissue-level yield strains were the same across all specimens. Comparison of the model predictions with the experimental data therefore enabled us to isolate the combined roles of volume fraction and architecture from the role of tissue material properties. Results indicated that, for both tensile and compressive loading, natural variations in volume fraction and architecture produced a negligible coefficient of variation (less than 3%) in apparent yield strains. Analysis of tissue-level strains showed that while bending of individual trabeculae played only a minor role in the apparent elastic behavior, the combined effects of this bending and tissue-level strength asymmetry produced apparent-level failure strains in compression that were 14% lower than those at the tissue level. By contrast, tissue and apparent-level yield strains were equivalent for tensile loading. We conclude that the uniformity of apparent yield strains is primarily the result of the highly oriented architecture that minimizes bending. Most of the variation that does occur is the result of the non-uniformity of the tissue-level yield strains.     

Trabecular shear stress in human vertebral cancellous bone: intra- and inter-individual variations

Correlation of the mean and standard deviation of trabecular stresses has been proposed as a mechanism by which a strong relationship between the apparent strength and stiffness of cancellous bone can be achieved. The current study examined whether the relationship between the mean and standard deviation of trabecular von Mises stresses can be generalized for any group of cancellous bone. Cylindrical human vertebral cancellous bone specimens were cut in the infero-superior direction from T12 of 23 individuals (inter-individual group). Thirty nine additional specimens were prepared similarly from the T4-T12 and L2-L5 vertebrae of a 63 year old male (intra-individual group). The specimens were scanned by micro-computed tomography (microCT) and trabecular von Mises stresses were calculated using finite element modeling. The expected value, standard deviation and coefficient of variation of the von Mises stress were calculated form a three-parameter Weibull function fitted to von Mises stress data from each specimen. It was found that the average and standard deviation of trabecular von Mises shear stress were: (i) correlated with each other, supporting the idea that high correlation between the apparent strength and stiffness of cancellous bone can be achieved through controlling the trabecular level shear stress variations, (ii) dependent on anatomical site and sample group, suggesting that the variation of stresses are correlated to the mean stress to different degrees between vertebrae and individuals, and (iii) dependent on bone volume fraction, consistent with the idea that shear stress is less well controlled in bones with low BV/TV. The conversion of infero-superior loading into trabecular von Mises stresses was maximum for the tissue at the junction of the thoracic and lumbar spine (T12-L1) consistent with this junction being a common site of vertebral fracture.     

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.     

Anisotropic post-yield response of cancellous bone simulated by stress–strain curves of bulk equivalent structures

During the last decade, finite element (FE) modelling has become ubiquitous in understanding complex mechanobiological phenomena, e.g. bone–implant interactions. The extensive computational effort required to achieve biorealistic results when modelling the post-yield behaviour of microstructures like cancellous bone is a major limitation of these techniques. This study describes the anisotropic biomechanical response of cancellous bone through stress–strain curves of equivalent bulk geometries. A cancellous bone segment, reverse engineered by micro computed tomography, was subjected to uniaxial compression. The material's constitutive law, obtained by nano-indentations, was considered during the simulation of the experimental process. A homodimensionally bulk geometry was employed to determine equivalent properties, resulting in a similar anisotropic response to the trabecular structure. The experimental verification of our model sustained that the obtained stress–strain curves can adequately reflect the post-yield behaviour of the sample. The introduced approach facilitates the consideration of nonlinearity and anisotropy of the tissue, while reducing the geometrical complexity of the model to a minimum.     

Biomechanical effect of mineral heterogeneity in trabecular bone

Due to daily loading, trabecular bone is subjected to deformations (i.e., strain), which lead to stress in the bone tissue. When stress and/or strain deviate from the normal range, the remodeling process leads to adaptation of the bone architecture and its degree of mineralization to effectively withstand the sustained altered loading. As the apparent mechanical properties of bone are assumed to depend on the degree and distribution of mineralization, the goal of the present study was examine the influences of mineral heterogeneity on the biomechanical properties of trabecular bone in the human mandibular condyle. For this purpose nine right condyles from human dentate mandibles were scanned and evaluated with a microCT system. Cubic regional volumes of interest were defined, and each was transformed into two different types of finite element (FE) models, one homogeneous and one heterogeneous. In the heterogeneous models the element tissue moduli were scaled to the local degree of mineralization, which was determined using microCT. Compression and shear tests were simulated to determine the apparent elastic moduli in both model types. The incorporation of mineralization variation decreased the apparent Young's and shear moduli by maximally 21% in comparison to the homogeneous models. The heterogeneous model apparent moduli correlated significantly with bone volume fraction and degree of mineralization. It was concluded that disregarding mineral heterogeneity may lead to considerable overestimation of apparent elastic moduli in FE models.     

The impact of the parameters of the constitutive model on the distribution of strain in the femoral head

The rapid spread of the finite element method has caused that it has become, among other methods, the standard tool for pre-clinical estimates of bone properties. This paper presents an application of this method for the calculation and prediction of strain and stress fields in the femoral head. The aim of the work is to study the influence of the considered anisotropy and heterogeneity of the modeled bone on the mechanical fields during a typical gait cycle. Three material models were tested with different properties of porous bone carried out in literature: a homogeneous isotropic model, a heterogeneous isotropic model, and a heterogeneous anisotropic model. In three cases studied, the elastic properties of the bone were determined basing on the Zysset-Curnier approach. The tensor of elastic constants defining the local properties of porous bone is correlated with a local porosity and a second order fabric tensor describing the bone microstructure. In the calculations, a model of the femoral head generated from high-resolution tomographic scans was used. Experimental data were drawn from publicly available database “Osteoporotic Virtual Physiological Human Project.” To realistically reflect the load on the femoral head, main muscles were considered, and their contraction forces were determined based on inverse kinematics. For this purpose, the results from OpenSim packet were used. The simulations demonstrated that differences between the results predicted by these material models are significant. Only the anisotropic model allowed for the plausible distribution of stresses along the main trabecular groups. The outcomes also showed that the precise evaluation of the mechanical fields is critical in the context of bone tissue remodeling under mechanical stimulations.

     

Trabecular bone remodeling around smooth and porous implants in an equine patellar model

The objective of this investigation was to examine the stress-morphology relationships for trabecular bone around implants with different surface characteristics. Stainless steel spheres with either a polished surface or a sintered-bead porous coating were implanted unilaterally into equine patellae and maintained for a 6 month period. Stereological methods were used to quantify the trabecular bone morphology and finite element analyses were performed to predict the trabecular bone stresses. In general, the remodeling response around the smooth implants was greater than that around those porous implants that exhibited bone ingrowth. In accordance with these differences, the finite element models predicted greater changes in the stresses adjacent to the smooth implants due to the nonlinear boundary conditions. However, it did not appear that the trajectorial theory, in its simplest form, was applicable to the remodeling induced by the implants. A linear relationship between the change in bone areal density and the change in von Mises effective stress provides support for the hypothesis that the architecture of trabecular bone corresponds to an optimal structure. The results also demonstrated that, under certain circumstances, small changes in the stress state may result in large changes in the principal material orientation.     

Influence of anisotropic bone properties on the biomechanical behavior of the acetabular cup implant: a multiscale finite element study

Although the biomechanical behavior of the acetabular cup (AC) implant is determinant for the surgical success, it remains difficult to be assessed due to the multiscale and anisotropic nature of bone tissue. The aim of the present study was to investigate the influence of the anisotropic properties of peri-implant trabecular bone tissue on the biomechanical behavior of the AC implant at the macroscopic scale. Thirteen bovine trabecular bone samples were imaged using micro-computed tomography (μCT) with a resolution of 18 μm. The anisotropic biomechanical properties of each sample were determined at the scale of the centimeter based on a dedicated method using asymptotic homogenization. The material properties obtained with this multiscale approach were used as input data in a 3D finite element model to simulate the macroscopic mechanical behavior of the AC implant under different loading conditions. The largest stress and strain magnitudes were found around the equatorial rim and in the polar area of the AC implant. All macroscopic stiffness quantities were significantly correlated (R² > 0.85, p < 6.5 e-6) with BV/TV (bone volume/total volume). Moreover, the maximum value of the von Mises stress field was significantly correlated with BV/TV (R² > 0.61, p < 1.6 e-3) and was always found at the bone-implant interface. However, the mean value of the microscopic stress (at the scale of the trabeculae) decrease as a function of BV/TV for vertical and torsional loading and do not depend on BV/TV for horizontal loading. These results highlight the importance of the anisotropic properties of bone tissue.     

A Femur-Implant Model for the Prediction of Bone Remodeling Behavior Induced by Cementless Stem

Bone remodeling simulation is an effective tool for the prediction of long-term effect of implant on the bone tissue, as well as the selection of an appropriate implant in terms of architecture and material. In this paper, a finite element model of proximal femur was developed to simulate the structures of internal trabecular and cortical bones by incorporating quantitative bone functional adaptation theory with finite element analysis. Cementless stems made of titanium, two types of Functionally Graded Material (FGM) and flexible ‘iso-elastic’ material as comparison were implanted in the structure of proximal femur respectively to simulate the bone remodeling behaviors of host bone. The distributions of bone density, von Mises stress, and interface shear stress were obtained. All the prosthetic stems had effects on the bone remodeling behaviors of proximal femur, but the degrees of stress shielding were different. The amount of bone loss caused by titanium implant was in agreement with the clinical observation. The FGM stems caused less bone loss than that of the titanium stem, in which FGM I stem (titanium richer at the top to more HAP/Col towards the bottom) could relieve stress shielding effectively, and the interface shear stresses were more evenly distributed in the model with FGM I stem in comparison with those in the models with FGM II (titanium and bioglass) and titanium stems. The numerical simulations in the present study provided theoretical basis for FGM as an appropriate material of femoral implant from a biomechanical point of view. The next steps are to fabricate FGM stem and to conduct animal experiments to investigate the effects of FGM stem on the remodeling behaviors using animal model.     

Stress and strain adaptation in load-dependent remodeling of the embryonic left ventricle

Altered pressure in the developing left ventricle (LV) results in altered morphology and tissue material properties. Mechanical stress and strain may play a role in the regulating process. This study showed that confocal microscopy, three-dimensional reconstruction, and finite element analysis can provide a detailed model of stress and strain in the trabeculated embryonic heart. The method was used to test the hypothesis that end-diastolic strains are normalized after altered loading of the LV during the stages of trabecular compaction and chamber formation. Stage-29 chick LVs subjected to pressure overload and underload at stage 21 were reconstructed with full trabecular morphology from confocal images and analyzed with finite element techniques. Measured material properties and intraventricular pressures were specified in the models. The results show volume-weighted end-diastolic von Mises stress and strain averaging 50–82 % higher in the trabecular tissue than in the compact wall. The volume-weighted-average stresses for the entire LV were 115, 64, and 147 Pa in control, underloaded, and overloaded models, while strains were 11, 7, and 4 %; thus, neither was normalized in a volume-weighted sense. Localized epicardial strains at mid-longitudinal level were similar among the three groups and to strains measured from high-resolution ultrasound images. Sensitivity analysis showed changes in material properties are more significant than changes in geometry in the overloaded strain adaptation, although resulting stress was similar in both types of adaptation. These results emphasize the importance of appropriate metrics and the role of trabecular tissue in evaluating the evolution of stress and strain in relation to pressure-induced adaptation.     

Flat medial-lateral conformity in total knee replacements does not minimize contact stresses.

The potential for wear in UHMWPE components for total knee replacements can be reduced by decreasing the stresses and strains arising from tibial-femoral contact. The conformity of the articular surfaces has a large effect on the resultant stresses, and components that achieve flat medial-lateral contact have been assumed to produce the lowest stresses due to their perfect conformity. We computed the stresses arising from curved and flat contact on a half-space using two-dimensional, plane strain elasticity solutions and finite element analyses to compare the performance of curved and flat indenters. These indenters were represented by a polynomial so the profiles could be continuously varied from curved to flat. Curved contact resulted in maximum stresses at the center of contact, while flat contact produced maximum stresses at the edge of contact. In addition, three contemporary tibial configurations (flat-on-flat, curved-on-flat, and curved-on-curved geometries) were analyzed using the finite element method with nonlinear material properties. The maximum contact stress, von Mises stress, and von Mises strain were lowest for the curved-on-curved model. The other configurations resulted in higher contact stresses, von Mises stresses, and von Mises strains. The perfect conformity arising from flat contact did not reduce the contact stresses in the UHMWPE component. The tensile stresses, however, were lowest for the flat-on-flat geometry compared with the other two configurations. Relating these distinct differences could prove useful in interpretation of data from simulator and retrieval studies.     

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

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

Spatial correlations of trabecular bone microdamage with local stresses and strains using rigid image registration

Although microdamage is known to accumulate in trabecular bone with overloading and aging, the tissue-level stresses and strains associated with local bone failure are not well known. Local correlation of microdamage with microstructural stresses and strains requires methods to accurately register histological sections with micro-computed tomography (micro-CT) based finite element models. In addition, the resolution of correlation (i.e., grid size) selected for analysis may affect the observed results. Therefore, an automated, repeatable, and accurate image registration algorithm was developed to determine the range of local stresses and strains associated with microdamage initiation. Using a two-dimensional rigid registration algorithm, bone structures from histology and micro-CT imaging were aligned. Once aligned, microdamaged regions were spatially correlated with local stresses and strains obtained from micro-CT based finite element analysis. Using this more sophisticated registration technique, we were able to analyze the effects of varying spatial grid resolution on local stresses and strains initiating microdamage. The results indicated that grid refinement to the individual pixel level (pixel-by-pixel method) more precisely defined the range of microdamage initiation compared to manually selected individual damaged and undamaged trabeculae. Using the pixel-by-pixel method, we confirmed that trabecular bone from younger cows sustained higher local strains prior to microdamage initiation compared to older bone.     

Determination of dynamically adapting anisotropic material properties of bone under cyclic loading

Because bone tissue adapts to loading conditions, finite element simulations of remodelling bone require a precise prediction of dynamically changing anisotropic elastic parameters. We present a phenomenological theory that refers to the tissue in terms of the tendency of the structure to align with principal stress directions. We describe the material parameters of remodelling bone. This work follows findings by the same research group and independently by Danilov (1971) in the field of plasticity, where the dependencies of the components of the stiffness tensor in terms of time are based on Hill's anisotropy. We modify such an approach in this novel theory that addresses bone tissue that can regenerate. The computational assumption of the theory is that bone trabeculae have the tendency to orient along one of the principal stress directions but during remodelling the principal stresses change continuously and the resulting orientation of the trabeculae can differ from the principal stress direction at any given time. The novelty of this work consists in the limited number of parameters needed to compute the twenty-one anisotropic material parameters at any given location in the bone tissue. In addition to the theory, we present here two cases of simplified geometry, loading and boundary conditions to show the effect of (1) time on the material properties; and (2) change of loading conditions on the anisotropic parameters. The long term goal is to experimentally verify that the predictions generated by theory provide a reliable simulation of cancellous bone properties.     

Comparison of HR-pQCT- and microCT-based finite element models for the estimation of the mechanical properties of the calcaneus trabecular bone

The calcaneus bone is formed of extensive trabecular bone and is therefore well suited to be used as an example of loaded bone to establish the ability of combining microfinite element (microFE) technique with high-resolution peripheral quantitative computed tomography (HR-pQCT) in determining its mechanical properties. HR-pQCT is increasingly used as a tool for in vivo bone clinical research, but its use has been limited to the distal radius and tibia. The goal of this study was to determine the applicability of HR-pQCT-derived microFE models of the calcaneus trabecular bone with 82 μm voxel size with reference to higher-resolution microCT-based models taken as gold standard. By comparing the outputs of microFE models generated from both HR-pQCT and microCT images of the trabecular bone of five calcaneus cadaveric specimens, it was found that the HR-pQCT-based models predicted mechanical properties for fracture load, total reaction force and von Mises stress are considerably different from microCT-based counterparts by 33, 64 and 70%, respectively. Also, the morphological analysis showed a comprehensive geometrical difference between HR-pQCT-based microFE models and their microCT-based equivalents. The results of the HR-pQCT-based models were found to have strong dependency on the threshold value chosen to binarise the images prior to finite element modelling. In addition, it was found that the voxel size has a strong impact on accuracy of imaged-based microFE models compared to other factors such as the presence of soft tissue and image scanning integration time. Therefore, although HR-pQCT has shown to be useful to predict overall structural and biomechanical changes, it is limited in providing local accurate biomechanical properties of trabecular bone and therefore should be used with caution when assessing bone remodelling through local changes of trabecular bone apposition and resorption in disease treatment monitoring.