Effect of boundary conditions on yield properties of human femoral trabecular bone

Trabecular bone plays an important mechanical role in bone fractures and implant stability. Homogenized nonlinear finite element (FE) analysis of whole bones can deliver improved fracture risk and implant loosening assessment. Such simulations require the knowledge of mechanical properties such as an appropriate yield behavior and criterion for trabecular bone. Identification of a complete yield surface is extremely difficult experimentally but can be achieved in silico by using micro-FE analysis on cubical trabecular volume elements. Nevertheless, the influence of the boundary conditions (BCs), which are applied to such volume elements, on the obtained yield properties remains unknown. Therefore, this study compared homogenized yield properties along 17 load cases of 126 human femoral trabecular cubic specimens computed with classical kinematic uniform BCs (KUBCs) and a new set of mixed uniform BCs, namely periodicity-compatible mixed uniform BCs (PMUBCs). In stress space, PMUBCs lead to 7–72 % lower yield stresses compared to KUBCs. The yield surfaces obtained with both KUBCs and PMUBCs demonstrate a pressure-sensitive ellipsoidal shape. A volume fraction and fabric-based quadric yield function successfully fitted the yield surfaces of both BCs with a correlation coefficient \(R^{2} \ge 0.93\). As expected, yield strains show only a weak dependency on bone volume fraction and fabric. The role of the two BCs in homogenized FE analysis of whole bones will need to be investigated and validated with experimental results at the whole bone level in future studies.     

A three-dimensional elastic plastic damage constitutive law for bone tissue

Motivated by mechanical analysis of bones and bone-implant systems, a 3D constitutive law describing the macroscopic mechanical behaviour of both cortical and trabecular bone in cyclic (not fatigue) overloads is developed. The proposed model which mathematical formulation is established within the framework of generalized standard materials accounts for three distinct material evolution modes where elastic, plastic and damage aspects are closely related. The anisotropic elasticity of bone is described by a morphology-based model and distinct damage behaviour in tension and compression by a halfspacewise generalized Hill criterion. The plastic criterion is based on the intact elastic compliance tensor. The algorithm applies three distinct projections based on the relationship between the internal variables and criteria. Their respective consistent tangent operators are presented. Numerical resolutions of several boundary value problems and a biomechanical application are presented to illustrate the potential of the constitutive model and demonstrate the expected quadratic convergence of the algorithm.     

An experimental and computational investigation of the post-yield behaviour of trabecular bone during vertebral device subsidence

Interbody fusion device subsidence has been reported clinically. An enhanced understanding of the mechanical behaviour of the surrounding bone would allow for accurate predictions of vertebral subsidence. The multiaxial inelastic behaviour of trabecular bone is investigated at a microscale and macroscale level. The post-yield behaviour of trabecular bone under hydrostatic and confined compression is investigated using microcomputed tomography-derived microstructural models, elucidating a mechanism of pressure-dependent yielding at the macroscopic level. Specifically, microstructural trabecular simulations predict a distinctive yield point in the apparent stress–strain curve under uniaxial, confined and hydrostatic compression. Such distinctive apparent stress–strain behaviour results from localised stress concentrations and material yielding in the trabecular microstructure. This phenomenon is shown to be independent of the plasticity formulation employed at a trabecular level. The distinctive response can be accurately captured by a continuum model using a crushable foam plasticity formulation in which pressure-dependent yielding occurs. Vertebral device subsidence experiments are also performed, providing measurements of the trabecular plastic zone. It is demonstrated that a pressure-dependent plasticity formulation must be used for continuum level macroscale models of trabecular bone in order to replicate the experimental observations, further supporting the microscale investigations. Using a crushable foam plasticity formulation in the simulation of vertebral subsidence, it is shown that the predicted subsidence force and plastic zone size correspond closely with the experimental measurements. In contrast, the use of von Mises, Drucker–Prager and Hill plasticity formulations for continuum trabecular bone models lead to over prediction of the subsidence force and plastic zone.     

A theoretical model to predict distribution of the fabric tensor and apparent density in cancellous bone

 The adaptation of cancellous bone to mechanical forces is well recognized. Theoretical models for predicting cancellous bone architecture have been developed and have mainly focused on the distribution of trabecular mass or the apparent density. The purpose of this study was to develop a theoretical model which can simultaneously predict the distribution of trabecular orthotropy/orientation, as represented by the fabric tensor, along with apparent density. Two sets of equations were derived under the assumption that cancellous bone is a biological self-optimizing material which tends to minimize strain energy. The first set of equations provide the relationship between the fabric tensor and stress tensor, and have been verified to be consistent with Wolff’s law of trabecular architecture, that is, the principal directions of the fabric tensor coincide with the principal stress trajectories. The second set of equations yield the apparent density from the stress tensor, which was shown to be identical to those obtained based on local optimization with strain energy density of true bone tissue as the objective function. These two sets of equations, together with elasticity field equations, provide a complete mathematical formulation for the adaptation of cancellous bone. Received: 25 February 1997/Revised version: 23 September 1997     

A Model of Bone Adaptation Using a Global Optimisation Criterion Based on the Trajectorial Theory of Wolff

Julius Wolff originally proposed that trabecular bone was influenced by mechanical stresses during the formative processes of growth and repair such that trabeculae were required to intersect at right angles. In this work, we have developed an analytical parametric microstructural model, which captures this restriction. Using homogenisation theory, a global material model was obtained. An optimal structure constructed of the homogenised material could then be found by optimising a cost function accounting for both the structural stiffness and the biological cost associated with metabolic maintenance of the bone tissue. The formulation was applied to an example problem of the proximal femur. Optimal densities and orientations were obtained for single load cases. The situation of multiple loads was also considered. In this case, we observe that the alignment of principal strains with the material orthotropy direction is, in general, not possible for all load cases. Thus less restrictive microstructures (nonorthotropic) will yield higher structural stiffnesses than strictly orthotropic microstructures.     

Numerical modeling of bone tissue adaptation—A hierarchical approach for bone apparent density and trabecular structure

In this work, a three-dimensional model for bone remodeling is presented, taking into account the hierarchical structure of bone. The process of bone tissue adaptation is mathematically described with respect to functional demands, both mechanical and biological, to obtain the bone apparent density distribution (at the macroscale) and the trabecular structure (at the microscale). At global scale bone is assumed as a continuum material characterized by equivalent (homogenized) mechanical properties. At local scale a periodic cellular material model approaches bone trabecular anisotropy as well as bone surface area density. For each scale there is a material distribution problem governed by density-based design variables which at the global level can be identified with bone relative density. In order to show the potential of the model, a three-dimensional example of the proximal femur illustrates the distribution of bone apparent density as well as microstructural designs characterizing both anisotropy and bone surface area density. The bone apparent density numerical results show a good agreement with Dual-energy X-ray Absorptiometry (DXA) exams. The material symmetry distributions obtained are comparable to real bone microstructures depending on the local stress field. Furthermore, the compact bone porosity is modeled giving a transversal isotropic behavior close to the experimental data. Since, some computed microstructures have no permeability one concludes that bone tissue arrangement is not a simple stiffness maximization issue but biological factors also play an important role.     

A Novel 3D Microstructural Model for Trabecular Bone: II. The Relationship Between Fabric and the Yield Surface

A novel 3D microstructural model was proposed and validated in part I of this publication. In part II, the model was used to identify the yield surface of a representative volume element of human trabecular bone as a function of volume fraction and degree of anisotropy. Finite element models of open and closed cells geometries were used to calculate effective yield stresses for a variety of loading cases with periodic boundary conditions. The postyield behaviour of the trabecular tissue was assumed from data available for cortical tissue. The yield stresses defined by a 0.2% offset in the global stress-strain curve were fit to an orthotropic Hill criterion and the parameters of the surface calculated. Similarly to the previous elastic analysis, distinct but strong relationships were obtained between volume fraction, fabric and the yield surface parameters for both the open and closed cell geometries. This finding suggests that volume fraction and fabric may be used to predict the initiation of mechanical damage in human trabecular bone at the continuum level.     

Homogenized stiffness matrices for mineralized collagen fibrils and lamellar bone using unit cell finite element models

Mineralized collagen fibrils have been usually analyzed like a two-phase composite material where crystals are considered as platelets that constitute the reinforcement phase. Different models have been used to describe the elastic behavior of the material. In this work, it is shown that when Halpin–Tsai equations are applied to estimate elastic constants from typical constituent properties, not all crystal dimensions yield a model that satisfy thermodynamic restrictions. We provide the ranges of platelet dimensions that lead to positive definite stiffness matrices. On the other hand, a finite element model of a mineralized collagen fibril unit cell under periodic boundary conditions is analyzed. By applying six canonical load cases, homogenized stiffness matrices are numerically calculated. Results show a monoclinic behavior of the mineralized collagen fibril. In addition, a 5-layer lamellar structure is also considered where crystals rotate in adjacent layers of a lamella. The stiffness matrix of each layer is calculated applying Lekhnitskii transformations, and a new finite element model under periodic boundary conditions is analyzed to calculate the homogenized 3D anisotropic stiffness matrix of a unit cell of lamellar bone. Results are compared with the rule-of-mixtures showing in general good agreement.     

Multi-axial mechanical properties of human trabecular bone

In the context of osteoporosis, evaluation of bone fracture risk and improved design of epiphyseal bone implants rely on accurate knowledge of the mechanical properties of trabecular bone. A multi-axial loading chamber was designed, built and applied to explore the compressive multi-axial yield and strength properties of human trabecular bone from different anatomical locations. A thorough experimental protocol was elaborated for extraction of cylindrical bone samples, assessment of their morphology by micro-computed tomography and application of different mechanical tests: torsion, uni-axial traction, uni-axial compression and multi-axial compression. A total of 128 bone samples were processed through the protocol and subjected to one of the mechanical tests up to yield and failure. The elastic data were analyzed using a tensorial fabric–elasticity relationship, while the yield and strength data were analyzed with fabric-based, conewise generalized Hill criteria. For each loading mode and more importantly for the combined results, strong relationships were demonstrated between volume fraction, fabric and the elastic, yield and strength properties of human trabecular bone. Despite the reviewed limitations, the obtained results will help improve the simulation of the damage behavior of human bones and bone-implant systems using the finite element method.     

The sensitivity of nonlinear computational models of trabecular bone to tissue level constitutive model

Microarchitectural finite element models have become a key tool in the analysis of trabecular bone. Robust, accurate, and validated constitutive models would enhance confidence in predictive applications of these models and in their usefulness as accurate assays of tissue properties. Human trabecular bone specimens from the femoral neck (n = 3), greater trochanter (n = 6), and lumbar vertebra (n = 1) of eight different donors were scanned by μ-CT and converted to voxel-based finite element models. Unconfined uniaxial compression and shear loading were simulated for each of three different constitutive models: a principal strain-based model, Drucker–Lode, and Drucker–Prager. The latter was applied with both infinitesimal and finite kinematics. Apparent yield strains exhibited minimal dependence on the constitutive model, differing by at most 16.1%, with the kinematic formulation being influential in compression loading. At the tissue level, the quantities and locations of yielded tissue were insensitive to the constitutive model, with the exception of the Drucker–Lode model, suggesting that correlation of microdamage with computational models does not improve the ability to discriminate between constitutive laws. Taken together, it is unlikely that a tissue constitutive model can be fully validated from apparent-level experiments alone, as the calculations are too insensitive to identify differences in the outcomes. Rather, any asymmetric criterion with a valid yield surface will likely be suitable for most trabecular bone models.     

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.     

Biaxial failure behavior of bovine tibial trabecular bone

Multiaxial failure properties of trabecular bone are important for modeling of whole bone fracture and can provide insight into structure-function relationships. There is currently no consensus on the most appropriate form of multiaxial yield criterion for trabecular bone. Using experimentally validated, high-resolution, non-linear finite element models, biaxial plain strain boundary conditions were applied to seven bovine tibial specimens. The dependence of multiaxial yield properties on volume fraction was investigated to quantify the interspecimen heterogeneity in yield stresses and strains. Two specimens were further analyzed to determine the yield properties for a wide range of biaxial strain loading conditions. The locations and quantities of tissue level yielding were compared for on-axis, transverse, and biaxial apparent level yielding to elucidate the micromechanical failure mechanisms. As reported for uniaxial loading of trabecular bone, the yield strains in multiaxial loading did not depend on volume fraction, whereas the yield stresses did. Micromechanical analysis indicated that the failure mechanisms in the on-axis and transverse loading directions were mostly independent. Consistent with this, the biaxial yield properties were best described by independent curves for on-axis and transverse loading. These findings establish that the multiaxial failure of trabecular bone is predominantly governed by the strain along the loading direction, requiring separate analytical expressions for each orthotropic axis to capture the apparent level yield behavior.     

Yield behavior of bovine cancellous bone

The compressive yield strain was measured for 61 specimens of bovine cancellous bone from three distal femora. There was no significant relationship (p = 0.08, R2 = 0.051) between yield strain and the degree of trabecular orientation. There was a significant positive correlation (p less than 0.00001, R2 = 0.319) between yield strain and structural (apparent) density and significant negative correlation (p less than 0.0025, R2 = 0.145) between yield strain and bone density. Yield strain correlated best with bone solid volume fraction Vv (epsilon y = 0.592 +2- 1.446vv, R2 = 0.337). The quantity, yield strain, is highly dependent on specific definitions of the yield point and the point of zero strain. For this study the yield point was defined by a 0.0003 offset criterion, and the point of zero strain was defined as the point where the tangent at 15 percent of yield crosses zero. The results using these definitions were compared with results using yield strain values determined by other definitions of the yield point and zero strain. The correlations between yield strain and trabecular orientation, structural density and bone density changed very little for differing definitions of yield. The results suggest that yield strain in cancellous bone is isotropic or independent of textural anisotropy, so the yield behaviour may be characterized by a maximum strain yield criterion. The results also suggest that the primary mode of yield in cancellous bone is buckling of the trabeculae.     

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.     

A global relationship between trabecular bone morphology and homogenized elastic properties.

An alternative concept of the relationship between morphological and elastic properties of trabecular bone is presented and applied to human tissue from several anatomical locations using a digital approach. The three-dimensional morphology of trabecular bone was assessed with a microcomputed tomography system and the method of directed secants as well as the star volume procedure were used to compute mean intercept length (MIL) and average bone length (ABL) of 4 mm cubic specimens. Assuming isotropic elastic properties for the trabecular tissue, the general elastic tensors of the bone specimens were determined using the homogenization method and the closest orthotropic tensors were calculated with an optimization algorithm. The assumption of orthotropy for trabecular bone was found to improve with specimen size and hold within 6.1 percent for a 4 mm cube size. A strong global relationship (r2 = 0.95) was obtained between fabric and the orthotropic elastic tensor with a minimal set of five constants. Mean intercept length and average bone length provided an equivalent power of prediction. These results support the hypothesis that the elastic properties of human trabecular bone from an arbitrary anatomical location can be estimated from an approximation of the anisotropic morphology and a prior knowledge of tissue properties.     

Micro-finite element simulation of trabecular-bone post-yield behaviour – effects of material model,element size and type

Micro-finite element (micro-FE) analysis became a standard tool for the evaluation of trabecular bone mechanical properties. The accuracy of micro-FE models for linear analyses is well established. However, the accuracy of recently developed nonlinear micro-FE models for simulations of trabecular bone failure is not known. In this study, a trabecular bone specimen was compressed beyond the apparent yield point. The experiment was simulated using different micro-FE meshes with different element sizes and types, and material models based on cortical bone. The results from the simulations were compared with experimental results to study the effects of the different element and material models. It was found that a decrease in element size from 80 to 40 μm had little effect on predicted post-yield behaviour. Element type and material model had significant effects. Nevertheless, none of the established material models for cortical bone were able to predict the typical descent in the load-displacement curve seen during compression of trabecular bone.     

Micro-finite element simulation of trabecular-bone post-yield behaviour--effects of material model, element size and type

Micro-finite element (micro-FE) analysis became a standard tool for the evaluation of trabecular bone mechanical properties. The accuracy of micro-FE models for linear analyses is well established. However, the accuracy of recently developed nonlinear micro-FE models for simulations of trabecular bone failure is not known. In this study, a trabecular bone specimen was compressed beyond the apparent yield point. The experiment was simulated using different micro-FE meshes with different element sizes and types, and material models based on cortical bone. The results from the simulations were compared with experimental results to study the effects of the different element and material models. It was found that a decrease in element size from 80 to 40 mum had little effect on predicted post-yield behaviour. Element type and material model had significant effects. Nevertheless, none of the established material models for cortical bone were able to predict the typical descent in the load-displacement curve seen during compression of trabecular bone.     

Application of the Tsai-Wu quadratic multiaxial failure criterion to bovine trabecular bone

As a first step toward development of a multiaxial failure criterion for human trabecular bone, the Tsai-Wu quadratic failure criterion was modified as a function of apparent density and applied to bovine tibial trabecular bone. Previous data from uniaxial compressive, tensile, and torsion tests (n = 139 total) were combined with those from new triaxial tests (n = 17) to calibrate and then verify the criterion. Combinations of axial compression and radial pressure were used to produce the triaxial compressive stress states. All tests were performed with minimal end artifacts in the principal material coordinate system of the trabecular network. Results indicated that the stress interaction term F12 exhibited a strong nonlinear dependence on apparent density (r2 > 0.99), ranging from -0.126 MPa-2 at low densities (0.29 g/cm3) to 0.005 MPa-2 at high densities (0.63 g/cm3). After calibration and when used to predict behavior of new-specimens without any curve-fitting, the Tsai-Wu criterion had a mean (+/- SD) error of -32.6 +/- 10.6 percent. Except for the highest density triaxial specimens, most (15/17 specimens) failed at axial stresses close to their predicted uniaxial values, and some reinforcement for transverse loading was observed. We conclude that the Tsai-Wu quadratic criterion, as formulated here, is at best only a reasonable predictor of the multiaxial failure behavior of trabecular bone, and further work is required before it can be confidently applied to human bone.     

Diagnostic uses of the Hill (Logit and Nernst) plots.

Analysis of the Hill plot for ligand-binding studies shows that it can adopt a variety of shapes other than a straight line. The shape of the curve can yield valuable information about the details of the binding process. In some cases it is possible to discriminate between models of co-operativity without the need for curve-fitting by computer.     

Computational simulation of simultaneous cortical and trabecular bone change in human proximal femur during bone remodeling

In this study, we developed a numerical framework that computationally determines simultaneous and interactive structural changes of cortical and trabecular bone types during bone remodeling, and we investigated the structural correlation between the two bone types in human proximal femur. We implemented a surface remodeling technique that performs bone remodeling in the exterior layer of the cortical bone while keeping its interior area unchanged. A micro-finite element (μFE) model was constructed that represents the entire cortical bone and full trabecular architecture in human proximal femur. This study simulated and compared the bone adaptation processes of two different structures: (1) femoral bone that has normal cortical bone shape and (2) perturbed femoral bone that has an artificial bone lump in the inferomedial cortex. Using the proposed numerical method in conjunction with design space optimization, we successfully obtained numerical results that resemble actual human proximal femur. The results revealed that actual cortical bone, as well as the trabecular bone, in human proximal femur has structurally optimal shapes, and it was also shown that a bone abnormality that has little contribution to bone structural integrity tends to disappear. This study also quantitatively determined the structural contribution of each bone: when the trabecular adaptation was complete, the trabecular bone supported 54% of the total load in the human proximal femur while the cortical bone carried 46%.