A nonlocal constitutive model for trabecular bone softening in compression

Factual observations of dynamic bone crushing

Dynamic bone-crushing, exemplified by the pig bone rib, is characterized thermo-mechanically in relation to the bone's microstructural characteristics. The cortical bone's dominant role consists of shielding the trabecular component by resisting deformation, sustaining high load levels, and ultimately cracking. Here we present a qualitative factual study to show that this behavior is the absolute opposite of its quasi-static counterpart in which the trabecular bone was found to play the dominant role. Using infrared thermography, we observed for the first time a significant localized temperature rise of up to 11 degrees Celsius in both cortical and trabecular damaging regions. Such observations call for additional clinically oriented research. Such a high contrast between static and dynamic failure mechanisms was not reported previously, and it paves the way for forensic-oriented studies in which the nature of the sustained load must be determined.

QCT-based computational bone strength assessment updated with MRI-derived 'hidden' microporosity

Microdamage accumulated through sustained periods of cyclic loading or single overloading events contributes to bone fragility through a reduction in stiffness and strength. Monitoring microdamage in vivo remains unattainable by clinical imaging modalities. As such, there are no established computational methods for clinical fracture risk assessment that account for microdamage that exists in vivo at any specific timepoint. We propose a method that combines multiple clinical imaging modalities to identify an indicative surrogate, which we term 'hidden porosity', that incorporates pre-existing bone microdamage in vivo. To do so, we use the third metacarpal bone of the equine athlete as an exemplary model for fatigue induced microdamage, which coalesces in the subchondral bone. N = 10 metacarpals were scanned by clinical quantitative computed tomography (QCT) and magnetic resonance imaging (MRI). We used a patch-based similarity method to quantify the signal intensity of a fluid sensitive MRI sequence in bone regions where microdamage coalesces. The method generated MRI-derived pseudoCT images which were then used to determine a pre-existing damage (Dpex) variable to quantify the proposed surrogate and which we incorporate into a nonlinear constitutive model for bone tissue. The minimum, median, and maximum detected Dpex of 0.059, 0.209, and 0.353 reduced material stiffness by 5.9%, 20.9%, and 35.3% as well as yield stress by 5.9%, 20.3%, and 35.3%. Limb-specific voxel-based finite element meshes were equipped with the updated material model. Lateral and medial condyles of each metacarpal were loaded to simulate physiological joint loading during gallop. The degree of detected Dpex correlated with a relative reduction in both condylar stiffness (p = 0.001, R2 > 0.74) and strength (p < 0.001, R2 > 0.80). Our results illustrate the complementary value of looking beyond clinical CT, which neglects the inclusion of microdamage due to partial volume effects. As we use clinically available imaging techniques, our results may aid research beyond the equine model on fracture risk assessment in human diseases such as osteoarthritis, bone cancer, or osteoporosis.

A new mathematical model of compressive stress-strain behaviour of low viscosity and high viscosity bone cement with different strain rates

A new mathematical model of compressive stress-strain behaviour of low viscosity (LV) and high viscosity (HV) bone cement has been proposed to capture large uniaxial deformation under constant applied strain rate by incorporating three-term power law. The modeling capacity of the proposed model has been validated using uniaxial compressive test under eight different low strain rates ranging from 1.39 × 10^-4 s^-1 to 3.53 × 10^-2 s^-1 for low viscosity and high viscosity bone cement. The well agreement between the model and experimental response suggests that the proposed model can successfully predict rate dependent deformation behavior for Poly(methyl methacrylate) (PMMA) bone cement. Additionally, the proposed model was compared with the generalized Maxwell viscoelastic model and found to be in good agreement. The comparison of compressive responses over low strain rates for LV and HV bone cement reveals their rate-dependent compressive yield stress behaviour along with a higher value of compressive yield stress of LV bone cement compared to HV bone cement. For example, at the strain rate of 1.39 × 10^-4 s^-1 the mean value of compressive yield stress of LV bone cement was found to be 64.46 MPa, whereas for HV bone cement it was 54.00 MPa. Moreover, the modeling of experimental compressive yield stress with the Ree-Eyring molecular theory suggests that the variation of yield stress of PMMA bone cement can be predicted using two processes Ree-Eyring theory. The proposed constitutive model might be useful to characterize large deformation behaviour with high accuracy for PMMA bone cement. Finally, both variants of PMMA bone cement also exhibit ductile-like compressive behaviour below the strain rate of 2.1 × 10^-2 s^-1, whereas above this threshold strain rate, brittle-like compressive failure behavior is observed.

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.

Nonlinear micro finite element models based on digital volume correlation measurements predict early microdamage in newly formed bone

Bone regeneration in critical-sized defects is a clinical challenge, with biomaterials under constant development aiming at enhancing the natural bone healing process. The delivery of bone morphogenetic proteins (BMPs) in appropriate carriers represents a promising strategy for bone defect treatment but optimisation of the spatial-temporal release is still needed for the regeneration of bone with biological, structural, and mechanical properties comparable to the native tissue. Nonlinear micro finite element (μFE) models can address some of these challenges by providing a tool able to predict the biomechanical strength and microdamage onset in newly formed bone when subjected to physiological or supraphysiological loads. Yet, these models need to be validated against experimental data. In this study, experimental local displacements in newly formed bone induced by osteoinductive biomaterials subjected to in situ X-ray computed tomography compression in the apparent elastic regime and measured using digital volume correlation (DVC) were used to validate μFE models. Displacement predictions from homogeneous linear μFE models were highly correlated to DVC-measured local displacements, while tissue heterogeneity capturing mineralisation differences showed negligible effects. Nonlinear μFE models improved the correlation and showed that tissue microdamage occurs at low apparent strains. Microdamage seemed to occur next to large cavities or in biomaterial-induced thin trabeculae, independent of the mineralisation. While localisation of plastic strain accumulation was similar, the amount of damage accumulated in these locations was slightly higher when including material heterogeneity. These results demonstrate the ability of the nonlinear μFE model to capture local microdamage in newly formed bone tissue and can be exploited to improve the current understanding of healing bone and mechanical competence. This will ultimately aid the development of BMPs delivery systems for bone defect treatment able to regenerate bone with optimal biological, mechanical, and structural properties.

Retrospective Evaluation and Framework Development of Bone Anisotropic Material Behavior Compared with Elastic, Elastic-Plastic, and Hyper-Elastic Properties

The main motivation for studying damage in bone tissue is to better understand how damage develops in the bone tissue and how it progresses. Such knowledge may help in the surgical aspects of joint replacement, fracture fixation or establishing the fracture tolerance of bones to prevent injury. Currently, there are no standards that create a realistic bone model with anisotropic material properties, although several protocols have been suggested. This study seeks to retrospectively evaluate the damage of bone tissue with respect to patient demography including age, gender, race, body mass index (BMI), height, and weight, and their role in causing fracture. Investigators believe that properties derived from CT imaging data to estimate the material properties of bone tissue provides more realistic models. Quantifying and associating damage with in vivo conditions will provide the required information to develop mathematical equations and procedures to predict the premature failure and potentially mitigate problems before they begin. Creating a realistic model for bone tissue can predict the premature failure(s), provide preliminary results before getting the surgery, and optimize the design of orthopaedic implants. A comparison was performed between the proposed model and previous efforts, where they used elastic, hyper- elastic, or elastic-plastic properties. Results showed that there was a significant difference between the anisotropic material properties of bone when compared with unrealistic previous methods. The results showed that the density is 50% higher in male subjects than female subjects. Additionally, the results showed that the density is 47.91% higher in Black subjects than Mixed subjects, 53.27% higher than Caucasian subjects and 57.41% higher than Asian. In general, race should be considered during modeling implants or suggesting therapeutic techniques.

Limit analysis of human proximal femur

A limit analysis numerical approach oriented to predict the peak/collapse load of human proximal femur, under two different loading conditions, is presented. A yield criterion of Tsai-Hu-type, expressed in principal stress space, is used to model the orthotropic bone tissues. A simplified human femur 3D model is envisaged to carry on numerical simulation of in-vitro tests borrowed from the relevant literature and to reproduce their findings. A critical discussion, together with possible future developments, is presented.

Length scale parameter of single trabecula in cancellous bone

In this paper, the material length scale parameter of the modified couple stress theory for trabecular bones is studied. For this reason, experimental data for the buckling of single wet and dry trabeculae from the subchondral region of the human medial tibial plateau are used from the literature. A material length scale parameter is extracted using the modified couple stress theory-based buckling relation and the experimental results. This parameter can capture the size-dependent behavior of trabeculae and can be beneficial for micro-mechanical investigation of trabecular bones. In addition, this paper proposes a size-dependent length scale parameter for trabeculae to estimate the buckling behaviors of dry and wet trabeculae, more accurately. Therefore, this paper confirms that the softening (weakening) effects of the modified couple stress theory and size-dependent behavior of the material length scale parameter can be considered for some special cases, and additionally, a correct value for the length scale parameter of the trabecula is revealed.

Bone Mechanical Properties in Healthy and Diseased States

The mechanical properties of bone are fundamental to the ability of our skeletons to support movement and to provide protection to our vital organs. As such, deterioration in mechanical behavior with aging and/or diseases such as osteoporosis and diabetes can have profound consequences for individuals' quality of life. This article reviews current knowledge of the basic mechanical behavior of bone at length scales ranging from hundreds of nanometers to tens of centimeters. We present the basic tenets of bone mechanics and connect them to some of the arcs of research that have brought the field to recent advances. We also discuss cortical bone, trabecular bone, and whole bones, as well as multiple aspects of material behavior, including elasticity, yield, fracture, fatigue, and damage. We describe the roles of bone quantity (e.g., density, porosity) and bone quality (e.g., cross-linking, protein composition), along with several avenues of future research.

Macrodamage Accumulation Model for a Human Femur

The objective of this study was to more fully understand the mechanical behavior of bone tissue that is important to find an alternative material to be used as an implant and to develop an accurate model to predict the fracture of the bone. Predicting and preventing bone failure is an important area in orthopaedics. In this paper, the macrodamage accumulation models in the bone tissue have been investigated. Phenomenological models for bone damage have been discussed in detail. In addition, 3D finite element model of the femur prepared from imaging data with both cortical and trabecular structures is delineated using MIMICS and ANSYS® and simulated as a composite structure. The damage accumulation occurring during cyclic loading was analyzed for fatigue scenario. We found that the damage accumulates sooner in the multiaxial than in the uniaxial loading condition for the same number of cycles, and the failure starts in the cortical bone. The damage accumulation behavior seems to follow a three-stage growth: a primary phase, a secondary phase of damage growth marked by linear damage growth, and a tertiary phase that leads to failure. Finally, the stiffness of the composite bone comprising the cortical and trabecular bone was significantly different as expected.

Experimental validation of finite element analysis of human vertebral collapse under large compressive strains

Osteoporosis-related vertebral fractures represent a major health problem in elderly populations. Such fractures can often only be diagnosed after a substantial deformation history of the vertebral body. Therefore, it remains a challenge for clinicians to distinguish between stable and progressive potentially harmful fractures. Accordingly, novel criteria for selection of the appropriate conservative or surgical treatment are urgently needed. Computer tomography-based finite element analysis is an increasingly accepted method to predict the quasi-static vertebral strength and to follow up this small strain property longitudinally in time. A recent development in constitutive modeling allows us to simulate strain localization and densification in trabecular bone under large compressive strains without mesh dependence. The aim of this work was to validate this recently developed constitutive model of trabecular bone for the prediction of strain localization and densification in the human vertebral body subjected to large compressive deformation. A custom-made stepwise loading device mounted in a high resolution peripheral computer tomography system was used to describe the progressive collapse of 13 human vertebrae under axial compression. Continuum finite element analyses of the 13 compression tests were realized and the zones of high volumetric strain were compared with the experiments. A fair qualitative correspondence of the strain localization zone between the experiment and finite element analysis was achieved in 9 out of 13 tests and significant correlations of the volumetric strains were obtained throughout the range of applied axial compression. Interestingly, the stepwise propagating localization zones in trabecular bone converged to the buckling locations in the cortical shell. While the adopted continuum finite element approach still suffers from several limitations, these encouraging preliminary results towards the prediction of extended vertebral collapse may help in assessing fracture stability in future work.

Locally measured microstructural parameters are better associated with vertebral strength than whole bone density

Whole vertebrae areal and volumetric bone mineral density (BMD) measurements are not ideal predictors of vertebral fractures. We introduce a technique which enables quantification of bone microstructural parameters at precisely defined anatomical locations. Results show that local assessment of bone volume fraction at the optimal location can substantially improve the prediction of vertebral strength.

INTRODUCTION

Whole vertebrae areal and volumetric BMD measurements are not ideal predictors of vertebral osteoporotic fractures. Recent studies have shown that sampling bone microstructural parameters in smaller regions may permit better predictions. In such studies, however, the sampling location is described only in general anatomical terms. Here, we introduce a technique that enables the quantification of bone volume fraction and microstructural parameters at precisely defined anatomical locations. Specific goals of this study were to investigate at what anatomical location within the vertebrae local bone volume fraction best predicts vertebral-body strength, whether this prediction can be improved by adding microstructural parameters and to explore if this approach could better predict vertebral-body strength than whole bone volume fraction and finite element (FE) analyses.

METHODS

Eighteen T12 vertebrae were scanned in a micro-computed tomography (CT) system and FE meshes were made using a mesh-morphing tool. For each element, bone microstructural parameters were measured and correlated with vertebral compressive strength as measured experimentally. Whole bone volume fraction and FE-predicted vertebral strength were also compared to the experimental measurements.

RESULTS

A significant association between local bone volume fraction measured at a specific central region and vertebral-body strength was found that could explain up to 90% of the variation. When including all microstructural parameters in the regression, the predictive value of local measurements could be increased to 98%. Whole bone volume fraction could explain only 64% and FE analyses 76% of the variation in bone strength.

CONCLUSIONS

A local assessment of volume fraction at the optimal location can substantially improve the prediction of bone strength. Local assessment of other microstructural parameters may further improve this prediction but is not clinically feasible using current technology.

Assessment of trabecular bone yield and post-yield behavior from high-resolution MRI-based nonlinear finite element analysis at the distal radius of premenopausal and postmenopausal women susceptible to osteoporosis

RATIONALE AND OBJECTIVES

To assess the performance of a nonlinear microfinite element model on predicting trabecular bone yield and post-yield behavior based on high-resolution in vivo magnetic resonance images via the serial reproducibility.

MATERIALS AND METHODS

The nonlinear model captures material nonlinearity by iteratively adjusting tissue-level modulus based on tissue-level effective strain. It enables simulations of trabecular bone yield and post-yield behavior from micro magnetic resonance images at in vivo resolution by solving a series of nonlinear systems via an iterative algorithm on a desktop computer. Measures of mechanical competence (yield strain/strength, ultimate strain/strength, modulus of resilience, and toughness) were estimated at the distal radius of premenopausal and postmenopausal women (N = 20, age range 50-75) in whom osteoporotic fractures typically occur. Each subject underwent three scans (20.2 ± 14.5 days). Serial reproducibility was evaluated via coefficient of variation (CV) and intraclass correlation coefficient (ICC).

RESULTS

Nonlinear simulations were completed in an average of 14 minutes per three-dimensional image data set involving analysis of 61 strain levels. The predicted yield strain/strength, ultimate strain/strength, modulus of resilience, and toughness had a mean value of 0.78%, 3.09 MPa, 1.35%, 3.48 MPa, 14.30 kPa, and 32.66 kPa, respectively, covering a substantial range by a factor of up to 4. Intraclass correlation coefficient ranged from 0.986 to 0.994 (average 0.991); CV ranged from 1.01% to 5.62% (average 3.6%), with yield strain and toughness having the lowest and highest CV values, respectively.

CONCLUSIONS

The data suggest that the yield and post-yield parameters have adequate reproducibility to evaluate treatment effects in interventional studies within short follow-up periods.

A novel approach to estimate trabecular bone anisotropy using a database approach

Continuum finite element (FE) models of bones have become a standard pre-clinical tool to estimate bone strength. These models are usually based on clinical CT scans and material properties assigned are chosen as isotropic based only on the density distribution. It has been shown, however, that trabecular bone elastic behavior is best described as orthotropic. Unfortunately, the use of orthotropic models in FE analysis derived from CT scans is hampered by the fact that the measurement of a trabecular orientation (fabric) is not possible from clinical CT images due to the low resolution of such images. In this study, we explore the concept of using a database (DB) of high-resolution bone models to derive the fabric information that is missing in clinical images. The goal of this study was to investigate if models with fabric derived from a relatively small database can already produce more accurate results than isotropic models. A DB of 33 human proximal femurs was generated from micro-CT scans with a nominal isotropic resolution of 82 µm. Continuum FE models were generated from the images using a pre-defined mesh template in combination with an iso-anatomic mesh morphing tool. Each element within the mesh template is at a specific anatomical location. For each element within the cancellous bone, a spherical region around the element centroid with a radius of 2mm was defined. Bone volume fraction and the mean-intercept-length fabric tensor were analyzed for that region. Ten femurs were used as test cases. For each test femur, four different models were generated: (1) an orthotropic model based on micro-CT fabric measurements (gold standard), (2) an orthotropic model based on the fabric derived from the best-matched database model, (3) an isotropic-I model in which the fabric tensor was set to the identity tensor, and (4) a second isotropic-II model with its total bone stiffness fitted to the gold standard. An elastic-plastic damage model was used to simulate failure and post failure behavior during a fall to the side. The results show that all models produce a similar stress distribution. However, compared to the gold standard, both isotropic-I and II models underestimated the stress/damage distributions significantly. We found no significant difference between DB-derived and gold standard models. Compared to the gold standard, the isotropic-I models further underestimated whole bone stiffness by 26.3% and ultimate load by 14.5%, while these differences for the DB-derived orthotropic models were only 4.9% and 3.1% respectively. The results indicate that the concept of using a DB to estimate patient-specific anisotropic material properties can considerably improve the results. We expect that this approach can lead to more accurate results in particular for cases where bone anisotropy plays an important role, such as in osteoporotic patients and around implants.

Towards patient-specific material modeling of trabecular bone post-yield behavior

Bone diseases such as osteoporosis are one of the main causes of bone fracture and often result in hospitalization and long recovery periods. Researchers are aiming to develop new tools that consider the multiple determinants acting at the different scales of bone, and which can be used to clinically estimate patient-specific fracture risk and also assess the efficacy of new therapies. The main step towards this goal is a deep understanding of the bone organ, and is achieved by modeling the complexity of the structure and the high variability of the mechanical outcome. This review uses a hierarchical approach to evaluate bone mechanics at the macroscale, microscale, and nanoscale levels and the interactions between scales. The first section analyzes the experimental evidence of bone mechanics in the elastic and inelastic regions, microdamage generation, and post-yield toughening mechanisms from the organ level to the ultrastructural level. On the basis of these observations, the second section provides an overview of the constitutive models available to describe bone mechanics and predict patient-specific outcomes. Overall, the role of the hierarchical structure of bone and the interplay between each level is highlighted, and their effect is evaluated in terms of modeling biological variability and patient specificity.

A non-linear homogeneous model for bone-like materials under compressive load

Finite element (FE) models accurately compute the mechanical response of bone and bone-like materials when the models include their detailed microstructure. In order to simulate non-linear behavior, which currently is only feasible at the expense of extremely high computational costs, coarser models can be used if the local morphology has been linked to the apparent mechanical behavior. The aim of this paper is to implement and validate such a constitutive law. This law is able to capture the non-linear structural behavior of bone-like materials through the use of fabric tensors. It also allows for irreversible strains using an elastoplastic material model incorporating hardening. These features are expressed in a constitutive law based on the anisotropic continuum damage theory coupled with isotropic elastoplasticity in a finite strain framework. This material model was implemented into metafor (LTAS-MNNL, University of Liège, Belgium), a non-linear FE software. The implementation was validated against experimental data of cylindrical samples subjected to compression. Three materials with bone-like microstructure were tested: aluminum foams of variable density (ERG, Oakland, CA, USA), polylactic acid foam (CERM, University of Liège, Liège, Belgium), and cancellous bone tissue of a deer antler (Faculty of Veterinary Medicine, University of Liège, Liège, Belgium).