Multilevel finite element modeling for the prediction of local cellular deformation in bone

The underlying mechanisms by which bone cells respond to mechanical stimuli or how mechanical loads act on osteocytes housed in lacunae in bone are not well understood. In this study, a multilevel finite element (FE) approach is applied to predict local cell deformations in bone tissue. The local structure of the matrix dictates the local mechanical environment of an osteocyte. Cell deformations are predicted from detailed linear FE analysis of the microstructure, consisting of an arrangement of cells embedded in bone matrix material. This work has related the loads applied to a whole femur during the stance phase of the gait cycle to the strain of a single lacuna and of canaliculi. The predicted bone matrix strains around osteocyte lacunae and canaliculi were nonuniform and differed significantly from the macroscopically measured strains. Peak stresses and strains in the walls of the lacuna were up to six times those in the bulk extracellular matrix. Significant strain concentrations were observed at sites where the process meets the cell body.     

Homogenization theory and digital imaging: A basis for studying the mechanics and design principles of bone tissue

Bone tissue is a complex multilevel composite which has the ability to sense ad respond to its mechanical environment. It is believed that bone cells called osteocytes within the bone matrix sense the mechanical environment and determine whether structural alterations are needed. At present it is not known, however, how loads are transferred from the whole bone level to cells. A computational procedure combining representative volume element (RVE) based homogenization theory with digital imaging is proposed to estimate strains at various levels of bone structure. Bone tissue structural organization and RVE based analysis are briefly reviewed. The digital image based computational procedure was applied to estimate strains in individual trabeculae (first-level microstructure). Homogenization analysis of an idealized model was used to estimate strains at one level of bone structure around osteocyte lacunae (second-level trabecular microstructure). The results showed that strain at one level of bone structure is amplified to a broad range at the next microstructural level. In one case, a zeor-level tensile principal strain of 495 muE engendered strains ranging between -1000 and 7000 muE in individual trabeculae (first-level microstructure). Subsequently, a first-level tensile principal strains of 1325 muE within an inidividual trabecula engendered strains ranging between 782 and 2530 muE around osteocyte lacunae. Lacunar orientation was found to influence strains around osteocyte lacunae much more than lacunar ellipticity. In conclusion, the computational procedure combining homogenization theory with digital imaging can proveide estimates of cell level strains within whole bones. Such results may be used to bridge experimental studies of bone adaptation at the whole bone and cell culture level. (c) 1994 John Wiley & Sons, Inc.     

Observations of microdamage around osteocyte lacunae in bone

Tensile microdamage was examined using laser scanning confocal microscopy in beam specimens of bovine, equine and human long bones loaded in vitro and whole specimens of rat ulnae loaded in vivo. Microcracks were observed to initiate frequently at osteocyte lacunae. The implication is that osteocyte lacunae act as stress concentrating features in bone. This association provides a potential mechanism for the detection of strain and/ or damage by osteocytes in bone.     

Bone remodelling algorithms incorporating both strain and microdamage stimuli

Biomechanical theories to predict bone remodelling have used either mechanical strain or microdamage as the stimulus driving cellular responses. Even though experimental data have implicated both stimuli in bone cell regulation, a mechano-regulatory system incorporating both stimuli has not yet been proposed. In this paper, we test the hypothesis that bone remodelling may be regulated by signals due to both strain and microdamage. Four mechano-regulation algorithms are studied where the stimulus is: strain, damage, combined strain/damage, and either strain or damage with damage-adaptive remodelling prioritised when damage is above a critical level. Each algorithm is implemented with both bone lining cell (surface) sensors and osteocyte cell (internal) sensors. Each algorithm is applied to prediction of a bone multicellular unit (BMU) remodelling on the surface of a bone trabecula. It is predicted that a regulatory system capable of responding to changes in either strain or microdamage but which prioritises removal of damaged bone when damage is above a critical level, is the only one that provides a plausible prediction of BMU behaviour. A mechanism for this may be that, below a certain damage threshold, osteocyte processes can sense changes in strain and fluid flow but above the threshold damage interferes with the signalling mechanism, or causes osteocyte apoptosis so that a remodelling response occurs to remove the dead osteocytes.     

Comparison of strain measurement in the mouse forearm using subject-specific finite element models,strain gaging,and digital image correlation

Mechanical loading in bone leads to the activation of bone-forming pathways that are most likely associated with a minimum strain threshold being experienced by the osteocyte. To investigate the correlation between cellular response and mechanical stimuli, researchers must develop accurate ways to measure/compute strain both externally on the bone surface and internally at the osteocyte level. This study investigates the use of finite element (FE) models to compute bone surface strains on the mouse forearm. Strains from three FE models were compared to data collected experimentally through strain gaging and digital image correlation (DIC). Each FE model was assigned subject-specific bone properties and consisted of one-dimensional springs representing the interosseous membrane. After three-point bending was performed on the ulnae and radii, moment of inertia was determined from microCT analysis of the bone region between the supports and then used along with standard beam analyses to calculate the Young’s modulus. Non-contact strain measurements from DIC were determined to be more suitable for validating numerical results than experimental data obtained through conventional strain gaging. When comparing strain responses in the three ulnae, we observed a 3–14% difference between numerical and DIC strains while the strain gage values were 37–56% lower than numerical values. This study demonstrates a computational approach for capturing bone surface strains in the mouse forearm. Ultimately, strains from these macroscale models can be used as inputs for microscale and nanoscale FE models designed to analyze strains directly in the osteocyte lacunae.     

Decrease in the osteocyte lacunar density accompanied by hypermineralized lacunar occlusion reveals failure and delay of remodeling in aged human bone

Aging decreases the human femur’s fatigue resistance, impact energy absorption, and the ability to withstand load. Changes in the osteocyte distribution and in their elemental composition might be involved in age‐related bone impairment. To address this question, we carried out a histomorphometric assessment of the osteocyte lacunar distribution in the periosteal and endosteal human femoral cortexes of 16 female and 16 male donors with regard to age‐ and sex‐related bone remodeling. Measurements of the bone mineral density distribution by quantitative backscattered electron imaging and energy dispersive X‐ray analysis were taken to evaluate the osteocyte lacunar mineral composition and characteristics. Age‐dependent decreases in the total osteocyte lacunar number were measured in all of the cases. This change signifies a risk for the bone’s safety. Cortical subdivision into periosteal and endosteal regions of interest emphasized that, in both sexes, primarily the endosteal cortex is affected by age‐dependent reduction in number of osteocyte lacunae, whereas the periosteal compartment showed a less pronounced osteocyte lacunar deficiency. In aged bone, osteocyte lacunae showed an increased amount of hypermineralized calcium phosphate occlusions in comparison with younger cases. With respect to Frost’s early delineation of micropetrosis, our microanalyses revealed that the osteocyte lacunae are subject to hypermineralization. Intralacunar hypermineralization accompanied by a decrease in total osteocyte lacunar density may contribute to failure or delayed bone repair in aging bone. A decreased osteocyte lacunar density may cause deteriorations in the canalicular fluid flow and reduce the detection of microdamage, which counteracts the bone’s structural integrity, while hypermineralized osteocyte lacunae may increase bone brittleness and render the bone fragile.     

Modeling the onset and propagation of trabecular bone microdamage during low-cycle fatigue

Relatively small amounts of microdamage have been suggested to have a major effect on the mechanical properties of bone. A significant reduction in mechanical properties (e.g. modulus) can occur even before the appearance of microcracks. This study uses a novel non-linear microdamaging finite-element (FE) algorithm to simulate the low-cycle fatigue behavior of high-density trabecular bone. We aimed to investigate if diffuse microdamage accumulation and concomitant modulus reduction, without the need for complete trabecular strut fracture, may be an underlining mechanism for low-cycle fatigue failure (defined as a 30% reduction in apparent modulus). A microCT constructed FE model was subjected to a single cycle monotonic compression test, and constant and variable amplitude loading scenarios to study the initiation and accumulation of low-cycle fatigue microdamage. Microcrack initiation was simulated using four damage criteria: 30%, 40%, 50% and 60% reduction in bone element modulus (el-MR). Evaluation of structural (apparent) damage using the four different tissue level damage criteria resulted in specimen fatigue failure at 72, 316, 969 and 1518 cycles for the 30%, 40%, 50% and 60% el-MR models, respectively. Simulations based on the 50% el-MR model were consistent with previously published experimental findings. A strong, significant non-linear, power law relationship was found between cycles to failure (N) and effective strain (Deltasigma/E(0)): N=1.394x10(-25)(Deltasigma/E(0))(-12.17), r(2)=0.97, p<0.0001. The results suggest that microdamage and microcrack propagation, without the need for complete trabecular strut fracture, are mechanisms for high-density trabecular bone failure. Furthermore, the model is consistent with previous numerical fatigue simulations indicating that microdamage to a small number of trabeculae results in relatively large specimen modulus reductions and rapid failure.     

Tensile behavior of cortical bone: dependence of organic matrix material properties on bone mineral content

A porous composite model is developed to analyze the tensile mechanical properties of cortical bone. The effects of microporosity (volksman's canals, osteocyte lacunae) on the mechanical properties of bone tissue are taken into account. A simple shear lag theory, wherein tensile loads are transferred between overlapped mineral platelets by shearing of the organic matrix, is used to model the reinforcement provided by mineral platelets. It is assumed that the organic matrix is elastic in tension and elastic-perfectly plastic in shear until it fails. When organic matrix shear stresses at the ends of mineral platelets reach their yield values, the stress-strain curve of bone tissue starts to deviate from linear behavior. This is referred as the microscopic yield point. At the point where the stress-strain behavior of bone shows a sharp curvature, the organic phase reaches its shear yield stress value over the entire platelet. This is referred as the macroscopic yield point. It is assumed that after macroscopic yield, mineral platelets cannot contribute to the load bearing capacity of bone and that the mechanical behavior of cortical bone tissue is determined by the organic phase only. Bone fails when the principal stress of the organic matrix is reached. By assuming that mechanical properties of the organic matrix are dependent on bone mineral content below the macroscopic yield point, the model is used to predict the entire tensile mechanical behavior of cortical bone for different mineral contents. It is found that decreased shear yield stresses and organic matrix elastic moduli are required to explain the mechanical behavior of bones with lowered mineral contents. Under these conditions, the predicted values (elastic modulus, 0.002 yield stress and strain, and ultimate stress and strain) are within 15% of experimental data.     

Three‐dimensional characterization of osteocyte volumes at multiple scales,and its relationship with bone biology and genome evolution in ray‐finned fishes

Osteocytes, cells embedded within the bone mineral matrix, inform on key aspects of vertebrate biology. In particular, a relationship between volumes of the osteocytes and bone growth and/or genome size has been proposed for several tetrapod lineages. However, the variation in osteocyte volume across different scales is poorly characterized and mostly relies on incomplete, two‐dimensional information. In this study, we characterize the variation of osteocyte volumes in ray‐finned fishes (Actinopterygii), a clade including more than half of modern vertebrate species in which osteocyte biology is poorly known. We use X‐ray synchrotron micro‐computed tomography (SRµCT) to achieve a three‐dimensional visualization of osteocyte lacunae and direct measurement of their size (volumes). Our specimen sample is designed to characterize variation in osteocyte lacuna morphology at three scales: within a bone, among the bones of one individual and among species. At the intra‐bone scale, we find that osteocyte lacunae vary noticeably in size between zones of organized and woven bone (being up to six times larger in woven bone), and across cyclical bone deposition. This is probably explained by differences in bone deposition rate, with larger osteocyte lacunae contained in bone that deposits faster. Osteocyte lacuna volumes vary 3.5‐fold among the bones of an individual, and this cannot readily be explained by variation in bone growth rate or other currently observable factors. Finally, we find that genome size provides the best explanation of variation in osteocyte lacuna volume among species: actinopterygian taxa with larger genomes (polyploid taxa in particular) have larger osteocyte lacunae (with a ninefold variation in median osteocyte volume being measured). Our findings corroborate previous two‐dimensional studies in tetrapods that also observed similar patterns of intra‐individual variation and found a correlation with genome size. This opens new perspectives for further studies on bone evolution, physiology and palaeogenomics in actinopterygians, and vertebrates as a whole.     

Measurement of microstructural strain in cortical bone

It is well known that mechanical factors affect bone remodeling such that increased mechanical demand results in net bone formation, whereas decreased demand results in net bone resorption. Current theories suggest that bone modeling and remodeling is controlled at the cellular level through signals mediated by osteocytes. The objective of this study was to investigate how macroscopically applied bone strains similar in magnitude to those that occur in vivo are manifest at the microscopic level in the bone matrix. Using a digital image correlation strain measurement technique, experimentally determined bone matrix strains around osteocyte lacuna resulting from macroscopic strains of approximately 2,000 microstrain (0.2%) reach levels of over 30,000 microstrain (3%) over fifteen times greater than the applied macroscopic strain. Strain patterns were highly heterogeneous and in some locations similar to observed microdamage around osteocyte lacuna indicating the resulting strains may represent the precursors to microdamage. This information may lead to a better understanding of how bone cells are affected by whole bone functional loading.     

Osteocyte lacuna size and shape in women with and without osteoporotic fracture

Osteocytes have been hypothesized to control the amount and location of bone tissue which is resorbed or formed, based on the strain magnitude they perceive, and therefore may play a role in the bone loss of osteoporosis. The shape of osteocyte lacunae influences the mechanical strain applied to the osteocyte; thus, it is important to quantify their shape to further understand the mechanical environment of this cell. Previous studies of the size and shape of lacunae have been contradictory and limited to two-dimensional measurements on iliac crest biopsies. This investigation measured the size and shape of osteocyte lacunae in trabecular bone near a typical fracture site from three-dimensional image sets obtained by confocal microscopy. Bone tissue specimens were obtained from individuals undergoing hip replacement subsequent to fracture, and matched cadaveric specimens without fracture. After extensive image processing to differentiate the lacunae from the matrix, the volume and anisotropy of the lacuna were determined. No significant difference was found in the size (volume) or shape (anisotropy) of the lacunae between women with and without osteoporotic fracture, although there was a large range of sizes and shapes in both groups. These results suggest that the size or shape of the lacunae, which influences the strain in osteocytes, does not play a role in osteoporotic fracture. In addition, this study provides geometric measures of lacunae that are important in computational modeling of the mechanical environment of osteocytes.     

Age-related properties at the microscale affect crack propagation in cortical bone

The increased risk for fracture with age is associated not only with reduced bone mass but also with impaired bone quality. At the microscale, bone quality is related to porosity, microstructural organization, accumulated microdamage and intrinsic material properties. However, the link between these characteristics and fracture behavior is still missing. Bone tissue has a complex structure and as age-related compositional and structural changes occur at all hierarchical length scales it is difficult to experimentally identify and discriminate the effect of each mechanism. The aim of this study was therefore to use computational models to analyze how microscale characteristics in terms of porosity, intrinsic toughness properties and microstructural organization affect the mechanical behavior of cortical bone. Tensile tests were simulated using realistic microstructural geometries based on microscopy images of human cortical bone. Crack propagation was modelled using the extended finite element method where cement lines surrounding osteons were modelled with an interface damage law to capture crack deflections along osteon boundaries. Both increased porosity and impaired material integrity resulted in straighter crack paths with cracks penetrating osteons, similar to what is seen experimentally for old cortical bone. However, only the latter predicted a more brittle failure behavior. Furthermore, the local porosity influenced the crack path more than the macroscopic porosity. In conclusion, age-related changes in cortical bone affect the crack path and the mechanical response. However, increased porosity alone was not driving damage in old bone, but instead impaired tissue integrity was required to capture brittle failure in aging bone.     

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.     

Postfailure modulus strongly affects microcracking and mechanical property change in human iliac cancellous bone: a study using a 2D nonlinear finite element method

A two-dimensional (2D) finite element (FE) method was used to estimate the ability of bone tissue to sustain damage as a function of postfailure modulus. Briefly, 2D nonlinear compact-tension FE models were created from quantitative back-scattered electron images taken of human iliac crest bone specimens. The effects of different postfailure moduli on predicted microcrack propagation were examined. The 2D FE models were used as surrogates for real bone tissues. The crack number was larger in models with higher postfailure modulus, while mean crack length and area were smaller in these models. The rate of stiffness reduction was greater in the models with lower postfailure modulus. Hence, the current results supported the hypothesis that hard tissue postfailure properties have strong effects on bone microdamage morphology and the rate of change in apparent mechanical properties.     

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.     

Finite element and experimental cortex strains of the intact and implanted tibia

Finite Element (FE) models for the simulation of intact and implanted bone find their main purpose in accurately reproducing the associated mechanical behavior. FE models can be used for preclinical testing of joint replacement implants, where some biomechanical aspects are difficult, if not possible, to simulate and investigate in vitro. To predict mechanical failure or damage, the models should accurately predict stresses and strains. Commercially available synthetic femur models have been extensively used to validate finite element models, but despite the vast literature available on the characteristics of synthetic tibia, numerical and experimental validation of the intact and implant assemblies of tibia are very limited or lacking. In the current study, four FE models of synthetic tibia, intact and reconstructed, were compared against experimental bone strain data, and an overall agreement within 10% between experimental and FE strains was obtained. Finite element and experimental (strain gauge) models of intact and implanted synthetic tibia were validated based on the comparison of cortex bone strains. The study also includes the analysis carried out on standard tibial components with cemented and noncemented stems of the P.F.C Sigma Modular Knee System. The overall agreement within 10% previously established was achieved, indicating that FE models could be successfully validated. The obtained results include a statistical analysis where the root-mean-square-error values were always <10%. FE models can successfully reproduce bone strains under most relevant acting loads upon the condylar surface of the tibia. Moreover, FE models, once properly validated, can be used for preclinical testing of tibial knee replacement, including misalignment of the implants in the proximal tibia after surgery, simulation of long-term failure according to the damage accumulation failure scenario, and other related biomechanical aspects.     

Osteocyte lacunae tissue strain in cortical bone

Current theories suggest that bone modeling and remodeling are controlled at the cellular level through signals mediated by osteocytes. However, the specific signals to which bone cells respond are still unknown. Two primary theories are: (1) osteocytes are stimulated via the mechanical deformation of the perilacunar bone matrix and (2) osteocytes are stimulated via fluid flow generated shear stresses acting on osteocyte cell processes within canaliculi. Recently, much focus has been placed on fluid flow theories since in vitro experiments have shown that bone cells are more responsive to analytically estimated levels of fluid shear stress than to direct mechanical stretching using macroscopic strain levels measured on bone in vivo. However, due to the complex microstructural organization of bone, local perilacunar bone tissue strains potentially acting on osteocytes cannot be reliably estimated from macroscopic bone strain measurements. Thus, the objective of this study was to quantify local perilacunar bone matrix strains due to macroscopically applied bone strains similar in magnitude to those that occur in vivo. Using a digital image correlation strain measurement technique, experimentally measured bone matrix strains around osteocyte lacunae resulting from macroscopic strains of approximately 2000 microstrain are significantly greater than macroscopic strain on average and can reach peak levels of over 30,000 microstrain locally. Average strain concentration factors ranged from 1.1 to 3.8, which is consistent with analytical and numerical estimates. This information should lead to a better understanding of how bone cells are affected by whole bone functional loading.     

High-resolution image-based simulation reveals membrane strain concentration on osteocyte processes caused by tethering elements

Osteocytes are vital for regulating bone remodeling by sensing the flow-induced mechanical stimuli applied to their cell processes. In this mechanosensing mechanism, tethering elements (TEs) connecting the osteocyte process with the canalicular wall potentially amplify the strain on the osteocyte processes. The ultrastructure of the osteocyte processes and canaliculi can be visualized at a nanometer scale using high-resolution imaging via ultra-high voltage electron microscopy (UHVEM). Moreover, the irregular shapes of the osteocyte processes and the canaliculi, including the TEs in the canalicular space, should considerably influence the mechanical stimuli applied to the osteocytes. This study aims to characterize the roles of the ultrastructure of osteocyte processes and canaliculi in the mechanism of osteocyte mechanosensing. Thus, we constructed a high-resolution image-based model of an osteocyte process and a canaliculus using UHVEM tomography and investigated the distribution and magnitude of flow-induced local strain on the osteocyte process by performing fluid–structure interaction simulation. The analysis results reveal that local strain concentration in the osteocyte process was induced by a small number of TEs with high tension, which were inclined depending on the irregular shapes of osteocyte processes and canaliculi. Therefore, this study could provide meaningful insights into the effect of ultrastructure of osteocyte processes and canaliculi on the osteocyte mechanosensing mechanism.

     

A nonlocal constitutive model for trabecular bone softening in compression

Using the three-dimensional morphological data provided by computed tomography, finite element (FE) models can be generated and used to compute the stiffness and strength of whole bones. Three-dimensional constitutive laws capturing the main features of bone mechanical behavior can be developed and implemented into FE software to enable simulations on complex bone structures. For this purpose, a constitutive law is proposed, which captures the compressive behavior of trabecular bone as a porous material with accumulation of irreversible strain and loss of stiffness beyond its yield point and softening beyond its ultimate point. To account for these features, a constitutive law based on damage coupled with hardening anisotropic elastoplasticity is formulated using density and fabric-based tensors. To prevent mesh dependence of the solution, a nonlocal averaging technique is adopted. The law has been implemented into a FE software and some simple simulations are first presented to illustrate its behavior. Finally, examples dealing with compression of vertebral bodies clearly show the impact of softening on the localization of the inelastic process.