Anisotropic properties of human tibial cortical bone as measured by nanoindentation

Fracture behavior of human cortical bone with high glycation content under dynamic loading

The present study simulates the fracture behavior of diabetic cortical bone with high levels of advanced glycation end-products (AGEs) under dynamic loading. We consider that the increased AGEs in diabetic cortical bone degrade the materials heterogeneity of cortical bone through a reduction in critical energy release rates of the microstructural features. To simulate the initiation and propagation of cracks, we implement a phase field fracture framework on 2D models of human tibia cortical microstructure. The simulations show that the mismatch between the fracture properties (e.g., critical energy release rate) of osteons and interstitial tissue due to high AGEs contents can change crack growth trajectories. The results show crack branching in the cortical microstructure under dynamic loading is affected by the mismatches related to AGEs. In addition, we observe cortical features such as osteons and cement lines can prevent multiple cracking under dynamic loading even with changing the mismatches due to high AGEs. Furthermore, under dynamic loading, some toughening mechanisms can be activated and deactivated with different AGEs contents. In conclusion, the current findings present that the combination of the loading type and materials heterogeneity of microstructural features can change the fracture response of diabetic cortical bone and its fragility.

Regional variation of the cortical and trabecular bone material properties in the rabbit skull

The material properties of some bones are known to vary with anatomical location, orientation and position within the bone (e.g., cortical and trabecular bone). Details of the heterogeneity and anisotropy of bone is an important consideration for biomechanical studies that apply techniques such as finite element analysis, as the outcomes will be influenced by the choice of material properties used. Datasets detailing the regional variation of material properties in the bones of the skull are sparse, leaving many finite element analyses of skulls no choice but to employ homogeneous, isotropic material properties, often using data from a different species to the one under investigation. Due to the growing significance of investigating the cranial biomechanics of the rabbit in basic science and clinical research, this study used nanoindentation to measure the elastic modulus of cortical and trabecular bone throughout the skull. The elastic moduli of cortical bone measured in the mediolateral and ventrodorsal direction were found to decrease posteriorly through the skull, while it was evenly distributed when measured in the anteroposterior direction. Furthermore, statistical tests showed that the variation of elastic moduli between separate regions (anterior, middle and posterior) of the skull were significantly different in cortical bone, but was not in trabecular bone. Elastic moduli measured in different orthotropic planes were also significantly different, with the moduli measured in the mediolateral direction consistently lower than that measured in either the anteroposterior or ventrodorsal direction. These findings demonstrate the significance of regional and directional variation in cortical bone elastic modulus, and therefore material properties in finite element models of the skull, particularly those of the rabbit, should consider the heterogeneous and orthotropic properties of skull bone when possible.

Anisotropy, Anatomical Region, and Additional Variables Influence Young's Modulus of Bone: A Systematic Review and Meta-Analysis

The importance of finite element analysis (FEA) is growing in orthopedic research, especially in implant design. However, Young's modulus (E) values, one of the most fundamental parameters, can range across a wide scale. Therefore, our study aimed to identify factors influencing E values in human bone specimens. We report our systematic review and meta-analysis based on the recommendation of the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) 2020 guideline. We conducted the analysis on November 21, 2021. We included studies investigating healthy human bone specimens and reported on E values regarding demographic data, specimen characteristics, and measurement specifics. In addition, we included study types reporting individual specimen measurements. From the acquired data, we created a cohort in which we performed an exploratory data analysis that included the explanatory variables selected by random forest and regression trees methods, and the comparison of groups using independent samples Welch's t test. A total of 756 entries were included from 48 articles. Eleven different bones of the human body were included in these articles. The range of E values is between 0.008 and 33.7 GPa. The E values were most heavily influenced by the cortical or cancellous type of bone tested. Measuring method (compression, tension, bending, and nanoindentation), the anatomical region within a bone, the position of the bone within the skeleton, and the bone specimen size had a decreasing impact on the E values. Bone anisotropy, specimen condition, patient age, and sex were selected as important variables considering the value of E. On the basis of our results, E values of a bone change with bone characteristics, measurement techniques, and demographic variables. Therefore, the evaluation of FEA should be performed after the standardization of in vitro measurement protocol. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Microstructural fatigue fracture behavior of glycated cortical bone

The current study aims to simulate fatigue microdamage accumulation in glycated cortical bone with increased advanced glycation end-products (AGEs) using a phase field fatigue framework. We link the material degradation in the fracture toughness of cortical bone to the high levels of AGEs in this tissue. We simulate fatigue fracture in 2D models of cortical bone microstructure extracted from human tibias. The results present that the mismatch between the critical energy release rate of microstructural features (e.g., osteons and interstitial tissue) can alter crack initiation and propagation patterns. Moreover, the high AGEs content through the increased mismatch ratio can cause the activation or deactivation of bone toughening mechanisms under cyclic loading. The fatigue fracture simulations also show that the lifetime of diabetic cortical bone samples can be dependent on the geometry of microstructural features and the mismatch ratio between the features. Additionally, the results indicate that the trapped cracks in cement lines in the diabetic cortical microstructure can prevent further crack growth under cyclic loading. The present findings show that alterations in the materials heterogeneity of microstructural features can change the fatigue fracture response, lifetime, and fragility of cortical bone with high AGEs contents. Cortical bone models are created from microscopy images taken from the cortical cross-section of human tibias. Increased glycation contents in the cortical bone sample can change the crack growth trajectories.

Multiscale theoretical model shows that aging-related mechanical degradation of cortical bone is driven by microstructural changes in addition to porosity

This study aims to gain mechanistic understanding of how aging-related changes in the microstructure of cortical bone drive mechanical consequences at the macroscale. To that end, cortical bone was modeled as a bundle of elastic-plastic, parallel fibers, which represented osteons and interstitial tissue, loaded in uniaxial tension. Distinct material properties were assigned to each fiber in either the osteon or interstitial fiber "families." Models representative of mature (20-60 yrs.) bone, and elderly (60+) bone were created by modeling aging via the following changes to the input parameters: (i) increasing porosity from 5% to 15%, (ii) increasing the ratio of the number of osteon fibers relative to interstitial fibers from 40% to 50%, and (iii) changing the fiber material properties from representing mature bone samples to representing elderly bone samples (i.e., increased strength and decreased toughness of interstitial fibers together with decreased toughness of osteon fibers). To understand the respective contributions of these changes, additional models isolating one or two of each of these were also created. From the computed stress-strain curve for the fiber bundle, the yield point (ϵy, σy), ultimate point (ϵu, σu), and toughness (UT) for the bundle as a whole were measured. We found that changes to all three input parameters were required for the model to capture the aging-related decline in cortical bone mechanical properties consistent with those previously reported in the literature. In both mature and elderly bundles, rupture of the interstitial fibers drove the initial loss of strength following the ultimate point. Plasticity and more gradual rupture of the osteons drove the remainder of the response. Both the onset and completion of interstitial fiber rupture occurred at lower strains in the elderly vs. mature case. These findings point to the importance of studying microstructural changes beyond porosity, such as the area fraction of osteons and the material properties of osteon and interstitial tissue, in order to further understanding of aging-related changes in bone.

Elasticity and material anisotropy of lamellar cortical bone in adult bovine tibia characterized via AFM nanoindentation

The research focuses on the evaluation of the mechanical properties of osteonal cortical bone at the lamellar level. Elastic properties of the mid-diaphysis region of the bovine tibia are investigated via cantilever-based nanoindentation at the submicron length scale utilizing Atomic Force Microscopy, where the force-displacement curves are used for the elastic assessment using the Derjaguin-Muller-Toropov model to calculate indentation modulus. Variations of the modulus and the directional mechanical response of the osteonal bone at different distances from the Haversian canal are investigated. Additionally, the effects of demineralization on the indentation modulus are discussed. It was found that in the axial direction, the first and last untreated thick lamella layers show a significant indentation modulus difference compared to all other layers (4.26 ± 0.4 and 4.6 ± 0.3 GPa vs ∼3.5 GPa). On the other hand, the indentation modulus of transverse thick lamella layers shows a periodic variation between ∼3 ± 0.7 GPa and ∼4 ± 0.3 GPa from near the Haversian canal to near the interstitial bone. A periodic variation in the anisotropy ratio was found. Mineral content was quantified via energy-dispersive X-ray microanalysis at different levels of mineralization and shows a positive correlation with the indentation modulus.

Effect of geometrical structure variations on strength and damage onset of cortical bone using multi-scale cohesive zone based finite element method

Three-dimensional multi-scale finite element models were designed to examine the effects of geometrical structure variations on the damage onset in cortical bone at multiple structural scales. A cohesive zone finite element approach, together with anisotropic damage initiation criteria, is used to predict the onset of damage. The finite element models are developed to account for the onset of microdamage from the microscopic length scales consisting of collagen fibres, to the macroscopic level consisting of osteons and the Haversian canals. Numerical results indicated that the yield strain at the initiation of microcracks is independent of variations in the local mineral volume fraction at each structural scale. Further, the yield strain and strength properties of cortical bone are dependent on its structural anisotropy and hierarchical structure. A positive correlation is observed between bone strength and mineral content at each length scale.

Material and nanomechanical properties of bone structural units of cortical and trabecular iliac bone tissues from untreated postmenopausal osteoporotic women

The differences in bone nanomechanical properties between cortical (Ct) and trabecular (Tb) bone remain uncertain, whereas knowing the respective contribution of each compartment is critical to understand the origin of bone strength. Our purpose was to compare bone mechanical and intrinsic properties of Ct and Tb compartments, at the bone structural unit (BSU) level, in iliac bone taken from a homogeneous untreated human population.

Among 60 PMMA-embedded transiliac bone biopsies from untreated postmenopausal osteoporotic women (64 ± 7 year-old), >2000 BSUs were analysed by nanoindentation in physiological wet conditions [indentation modulus (elasticity), hardness, dissipated energy], by Fourier transform infrared (FTIRM) and Raman microspectroscopy (mineral and organic characteristics), and by X-ray microradiography (degree of mineralization of bone, DMB). BSUs were categorized based on tissue age, osteonal (Ost) and interstitial (Int) tissues location and bone compartments (Ct and Tb).

Indentation modulus was higher in Ct than in Tb BSUs, both in Ost and Int. dissipated energy was higher in Ct than Tb, in Int BSUs. Hardness was not different between Ct and Tb BSUs. In Ost or Int BSUs, mineral maturity (conversion of non-apatitic into apatitic phosphates) was higher in Ct than in Tb, as well as for collagen maturity (Ost). Mineral content assessed as mineral/matrix (FTIRM and Raman) or as DMB, was lower in Ct than in Tb. Crystallinity (FTIRM) was similar in BSUs from Ct and Tb, and slightly lower in Ct than in Tb when measured by Raman, indicating that the crystal size/perfection was quite similar between Ct and Tb BSUs. The differences found between Ost and Int tissues were much higher than the difference found between Ct and Tb for all those bone material properties. Multiple regression analysis showed that Indentation modulus and dissipated energy were mainly explained by mineral maturity in Ct and by collagen maturity in Tb, and hardness by mineral content in both Ct and Tb.

In conclusion, in untreated human iliac bone, Ct and Tb BSUs exhibit different characteristics. Ct BSUs have higher indentation modulus, dissipated energy (Int), mineral and organic maturities than Tb BSUs, without difference in hardness. Although those differences are relatively small compared to those found between Ost and Int BSUs, they may influence bone strength at macroscale.

Subchondral bone plate thickness is associated with micromechanical and microstructural changes in the bovine patella osteochondral junction with different levels of cartilage degeneration

The influence of joint degeneration on the biomechanical properties of calcified cartilage and subchondral bone plate at the osteochondral junction is relatively unknown. Common experimental difficulties include accessibility to and visualization of the osteochondral junction, application of mechanical testing at the appropriate length scale, and availability of tissue that provides a consistent range of degenerative changes. This study addresses these challenges. A well-established bovine patella model of early joint degeneration was employed, in which micromechanical testing of fully hydrated osteochondral sections was carried out in conjunction with high-resolution imaging using differential interference contrast (DIC) optical light microscopy. A total of forty-two bovine patellae with different grades of tissue health ranging from healthy to mild, moderate, and severe cartilage degeneration, were selected. From the distal-lateral region of each patella, two adjacent osteochondral sections were obtained for the mechanical testing and the DIC imaging, respectively. Mechanical testing was carried out using a robotic micro-force acquisition system, applying compression tests over an array (area: 200 μm × 1000 μm, step size: 50 μm) across the osteochondral junction to obtain a stiffness map. Morphometric analysis was performed for the DIC images of fully hydrated cryo-sections. The levels of cartilage degeneration, DIC images, and the stiffness maps were used to associate the mechanical properties onto the specific tissue regions of cartilage, calcified cartilage, and subchondral bone plate. The results showed that there were up to 20% and 24% decreases (p < 0.05) in the stiffness of calcified cartilage and subchondral bone plate, respectively, in the severely degenerated group compared to the healthy group. Furthermore, there were increases (p < 0.05) in the number of tidemarks, bone spicules at the cement line, and the mean thickness of the subchondral bone plate with increasing levels of degeneration. The decreasing stiffness in the subchondral bone plate coupled with the presence of bone spicules may be indicative of a subchondral remodeling process involving new bone formation. Moreover, the mean thickness of the subchondral bone plate was found to be the strongest indicator of mechanical and associated structural changes in the osteochondral joint tissues.

Damage analysis of human cortical bone under compressive and tensile loadings

Developing advanced fracture tools can increase the understanding of crack growth trajectories in human cortical bone. The present study investigates fracture micromechanics of human cortical bone under compressive and tensile loadings utilizing a phase field method. We construct two-dimensional finite element models from cortical microstructure of a human tibia cross section. We apply compression on the cortical bone models to create compressive microcracks. Then, we simulate the fracture of these models under tension to discover influential parameters on microcracks formation and post-yielding behavior. The results show that cement lines are susceptible sites to damage nucleation under compression rather than tension. The findings of this study also indicate a higher accumulation of initial damage (induced by compression) can lead to a lower microscopic stiffness as well as a less resistant material to damage initiation under tension. The simulations further indicate that the post-yielding properties (e.g., toughness) can be dependent on different variables such as morphological information of the osteons, the initial accumulation of microcracks, and the total length of cement lines.

Finite Element Prediction on Fracture Load of Femur with Osteogenesis Imperfecta under Various Loading Conditions

Osteogenesis Imperfecta (OI) is an inherited disorder characterized by extreme bone fragility due to collagen defects. It is an incurable disease. Bone fractures can occur frequently without prior notice, especially among children. Early quantitative prediction of fracture loads due to OI tends to alert patients to avoid unnecessary situations or dangerous conditions. This study is aimed at investigating the fracture loads of femur with OI under various types of loading. Ten finite element models of an OI-affected bone were reconstructed from the normal femur with different bowing angles ranging from 7.5 to 30.0°. The boundary conditions were assigned on an OI-affected femoral head under three types of load: medial-lateral impacts, compression-tension, and internal-external torsions, and various loading direction cases that reflect the stance condition. The fracture load was examined based on the load that can cause bone fracture for each case. The results show that the loads bearable by the femur before fracture were decreased with respect to the increase of OI bowing angles in most of the loading cases. The risk of fracture for the femur with OI was directly proportional to the increase of bowing angles in the frontal plane. This study provides new insights on fracture load prediction in OI-affected bone with respect to various loading types, which could help medical personnel for surgical intervention judgement.

Influence of Molybdenum on the Microstructure, Mechanical Properties and Corrosion Resistance of Ti20Ta20Nb20(ZrHf)20-xMox (Where: x = 0, 5, 10, 15, 20) High Entropy Alloys

The presented work was focused on investigating the influence of the (hafnium and zirconium)/molybdenum ratio on the microstructure and properties of Ti20Ta20Nb20(ZrHf)20-xMox (where: x = 0, 5, 10, 15, 20 at.%) high entropy alloys in an as-cast state. The designed chemical composition was chosen due to possible future biomedical applications. Materials were obtained from elemental powders by vacuum arc melting technique. Phase analysis revealed the presence of dual body-centered cubic phases. X-ray diffraction showed the decrease of lattice parameters of both phases with increasing molybdenum concentration up to 10% of molybdenum and further increase of lattice parameters. The presence of two-phase matrix microstructure and hafnium and zirconium precipitates was proved by scanning and transmission electron microscopy observation. Mechanical property measurements revealed decreased micro- and nanohardness and reduced Young's modulus up to 10% of Mo content, and further increased up to 20% of molybdenum addition. Additionally, corrosion resistance measurements in Ringers' solution confirmed the high biomedical ability of studied alloys due to the presence of stable oxide layers.

Computational study of the mechanical influence of lacunae and perilacunar zones in cortical bone microcracking

The mechanical behavior of cortical bone is influenced by microstructural components such as osteons, Haversian canals, and osteocyte lacunae that arise from biological remodeling processes. This study takes a computational approach to investigate the role of the perilacunar zones formed by the local remodeling processes of lacunar-dwelling osteocytes by utilizing phase-field finite element models based on histological imaging of human bone. The models simulated the microdamage accumulation that occurs in cortical bone under transverse compression in bone without lacunae, with lacunae, and with a perilacunar zone surrounding lacunae in order to investigate the role of these features. The results of the simulations found that while lacunae create stress concentration which initiate further damage, perilacunar regions can delay or prevent the emergence and growth of microcracks.

Experimental study of anisotropic constitutive behavior of β-HMX crystals via nanoindentation and small-scale dynamic impact

For energetic crystals such as HMX, the sensitivity of the material to shock, the possibility of initiation, and the subsequent reaction is known to be controlled by processes occurring at the crystal level. The anisotropic nature of β-HMX can be critical in determining the performance of HMX based polymer bonded explosives, which are widely used across multiple industries as propellant or explosives. In this work, we experimentally obtain constitutive parameters for characterizing the response of multiple crystalline planes of β-HMX crystals to external loading. Nanoindentation and small-scale dynamic impact experiments were performed on multiple planes of β-HMX crystals to comparatively measure the indentation moduli in multiple orientation directions. Anisotropic material behavior, involving constitutive elastic and non-elastic parameters, was measured and studied. Findings regarding material properties for the (100), (010), (001), {110}, and {011} planes of β-HMX are presented and compared with literature data.

Fracture behavior of human cortical bone: Role of advanced glycation end-products and microstructural features

Diabetes is associated with increased fracture risk in human bone, especially in the elderly population. In the present study, we investigate how simulated advanced glycation end-products (AGEs) and materials heterogeneity affect crack growth trajectory in human cortical bone. We used a phase field fracture framework on 2D models of cortical microstructure created from human tibias to analyze crack propagation. The increased AGEs level results in a higher rate of crack formation. The simulations also indicate that the mismatch between the fracture properties (e.g., critical energy release rate) of osteons and interstitial tissue can alter the post-yielding behavior. The results show that if the critical energy release rate of cement lines is lower than that of osteons and the surrounding interstitial matrix, cracks can be arrested by cement lines. Additionally, activation of toughening mechanisms such as crack merging and branching depends on bone microstructural morphology (i.e., osteons geometrical parameters, canals, and lacunae porosities). In conclusion, the present findings suggest that materials heterogeneity of microstructural features and the crack-microstructure interactions can play important roles in bone fragility.

Effects of Fatigue Damage on the Microscopic Modulus of Cortical Bone Using Nanoindentation

Alterations to the bone structure from cycle loadings can undermine its damage resistance at multiple scales. The accumulation of fatigue damage in a bone is commonly characterized by the reduction in the elastic modulus. In this study, nano-indentation was used for investigating microscopic damage evolution of bovine tibia samples subjected to fatigue loading. Indentation tests were conducted in the same 60 μm × 120 μm area with different degrees of damage, including fracture, and the evolution of reduced modulus was observed. The results showed that bone’s reduced modulus decreased significantly during the initial 40% of the life fraction, whereas it proceeded slowly during the remaining period. As the size of the residual indentations was about 4 μm in length, the degradation of bone’s reduced modulus reflected the accumulation of fatigue damage at smaller scales.

Microtensile properties and failure mechanisms of cortical bone at the lamellar level

Bone features a remarkable combination of toughness and strength which originates from its complex hierarchical structure and motivates its investigation on multiple length scales. Here, in situ microtensile experiments were performed on dry ovine osteonal bone for the first time at the length scale of a single lamella. The micromechanical response was brittle and revealed larger ultimate tensile strength compared to the macroscale (factor of 2.3). Ultimate tensile strength for axial and transverse specimens was 0.35 ± 0.05 GPa and 0.13 ± 0.02 GPa, respectively. A significantly greater strength anisotropy relative to compression was observed (axial to transverse strength ratio of 2.7:1 for tension, 1.3:1 for compression). Fracture surface and transmission electron microscopic analysis suggested that this may be rationalized by a change in failure mode from fibril-matrix interfacial shearing for axial specimens to fibril-matrix debonding in the transverse direction. An improved version of the classic Hashin's composite failure model was applied to describe lamellar bone strength as a function of fibril orientation. Together with our experimental observations, the model suggests that cortical bone strength at the lamellar level is remarkably tolerant to variations of fibrils orientation of about ±30°. This study highlights the importance of investigating bone's hierarchical organization at several length scales for gaining a deeper understanding of its macroscopic fracture behavior. STATEMENT OF SIGNIFICANCE: Understanding bone deformation and failure behavior at different length scales of its hierarchical structure is fundamental for the improvement of bone fracture prevention, as well as for the development of multifunctional bio-inspired materials combining toughness and strength. The experiments reported in this study shed light on the microtensile properties of dry primary osteonal bone and establish a baseline from which to start further investigations in more physiological conditions. Microtensile specimens were stronger than their macroscopic counterparts by a factor of 2.3. Lamellar bone strength seems remarkably tolerant to variations of the sub-lamellar fibril orientation with respect to the loading direction (±30°). This study underlines the importance of studying bone on all length scales for improving our understanding of bone's macroscopic mechanical response.

Hardness, an Important Indicator of Bone Quality, and the Role of Collagen in Bone Hardness

Bone is a nanocomposite material where the hard inorganic (hydroxyapatite crystallites) and organic (collagen fibrils) components are hierarchically arranged in the nanometer scale. Bone quality is dependent on the spatial distributions in the shape, size and composition of bone constituents (mineral, collagen and water). Bone hardness is an important property of bone, which includes both elastic and plastic deformation. In this study, a microhardness test was performed on a deer bone samples. The deer tibia shaft (diaphysis) was divided into several cross-sections of equal thickness; samples were prepared in untreated, boiled water treatment (100 °C for 30 min) and sodium hypochlorite (NaOCl) treatment conditions. Microhardness tests were performed on various regions of the tibial diaphysis to study the heterogeneous characteristics of bone microhardness and highlight the role of the organic matrix in bone hardness. The results indicated that boiled water treatment has a strong negative correlation with bone hardness. The untreated bone was significantly (+20%) harder than the boiled-water-treated bone. In general, the hardness values near the periosteal surface was significantly (23 to 45%) higher than the ones near the endosteal surface. Samples treated with NaOCl showed a significant reduction in hardness.

Effect of geometrical structure variations on the viscoelastic and anisotropic behaviour of cortical bone using multi-scale finite element modelling

Multi-scale finite element analysis is performed to ascertain the effect of geometrical changes at multiple structural scales on the mechanical properties of cortical bone. Finite element models are developed, with reference to experimental data from existing literature, to account for bone's viscoelastic behaviour and anisotropic structure from the most fundamental level of bone consisting of mineralised collagen fibrils, up to the macroscopic level consisting of osteons and the Haversian canals. A statistical approach is incorporated to perform sensitivity analyses on the effects of different geometrical parameters on the effective material properties of cortical bone at each length scale. Numerical results indicate that there is an exponential correlation between the mineral volume fraction and the effective stiffness constants at each length scale. This contributes to the exponential behaviour of the instantaneous moduli describing cortical bone's two-phase stress relaxation process: a fast and slow response relaxation behaviour. Results indicate that the fast response relaxation time is independent of bone's structural anisotropy, whilst being dependent on variations in the global mineral volume fraction between length scales. However, the slow response relaxation time is independent of the changes in mineral volume fraction. It is also observed that the slow response relaxation time varies with bone's anisotropic structure, and therefore, contributes to the anisotropic properties of bone.

Investigation of bone matrix composition, architecture and mechanical properties reflect structure-function relationship of cortical bone in glucocorticoid induced osteoporosis

Glucocorticoid induced osteoporosis (GIOP) is the most common negative consequence of long-term glucocorticoid treatment, leading to increased fracture risk followed by loss of mobility and high mortality risk. These biologically induced changes in bone quality at molecular level lead to changes both in bone matrix architecture and bone matrix composition. However, the quantitative details of changes in bone quality - and especially their link to reduced macroscale mechanical properties are still largely missing. In this study, a mouse model for glucocorticoid-induced osteoporosis (GIOP) was used to investigate mechanical and material alterations in bone cortex (natural nanocomposite) at different scale. By combining quantitative backscattered electron (qBSE) imaging, nanoindentation and high brilliance synchrotron X-ray nanomechanical imaging on a genetically modified mouse model of GIOP, we were able to quantify the local indentation modulus, mineralization distribution and the alterations of nanoscale structures and deformation mechanisms in the mid-diaphysis of femur, and relate them to the macroscopic mechanical changes. Our results showed clear and significant changes in terms of material quality of bone at nanoscale and microscale, which manifests itself in development of spatial heterogeneities in mineralization and indentation moduli across the bone organ, with potential implications for increased fracture risk.

An indentation-based approach to determine the elastic constants of soft anisotropic tissues

Characterization of the mechanical properties of tissue can help to understand tissue mechanobiology, including disease diagnosis and progression. Indentation is increasingly used to measure the local mechanical properties of tissue, but it has not been fully adapted to capture anisotropic properties. This paper presents an indentation-based method to measure elastic constants of soft anisotropic tissues without additional mechanical tests. The approach uses measurement of the indentation modulus and the aspect ratio of the elliptical contact introduced by anisotropic mechanical properties of tissue to determine the elastic constants from finite element analysis. The imprinted area imparted by a fluorescent bead-coated spherical indenter showed the aspect ratio of the contact area, giving a generalized sense of the level of anisotropy, and instrumented indentation determined the indentation modulus. A parametric study using finite element simulation of the indentation tests established the relationship between the aspect ratio of contact and the non-dimensional ratios, E_x/E_y and G_xy/E_y; here, E_x and E_y are the Young's moduli (E_x > E_y) and G_xy is the shear modulus in the xy plane. For strongly anisotropic materials (E_x/E_y > 150), aspect ratio and indentation modulus are sufficient to determine G_xy and E_y. For weakly anisotropic materials, indentation modulus in the transverse direction, E_y, and the aspect ratio of contact in the anisotropic plane can be used to determine the elastic constants. The proposed approach improves the elastic characterization of soft, anisotropic biological materials from indentation and helps to elucidate the complex mechanical behavior of soft anisotropic tissues.

Functional disuse initiates medullary endosteal micro-architectural impairment in cortical bone characterized by nanoindentation

In this study, we evaluated the effect of functional disuse-induced bone remodeling on its mechanical properties, individually at periosteum and medullary endosteum regions of the cortical bone. Left middle tibiae were obtained from 5-month-old female Sprague-Dawley rats for the baseline control as well as hindlimb suspended (disuse) groups. Micro-nano-mechanical elastic moduli (at lateral region) was evaluated along axial (Z), circumferential (C) and radial (R) orientations using nanoindentation. Results indicated an anisotropic microstructure with axial orientation having the highest and radial orientation with the lowest moduli at periosteum and medullary endosteum for both baseline control as well as disuse groups. Between the groups: at periosteum, an insignificant difference was evaluated for each of the orientations (p > 0.05) and at endosteum, a significant decrease of elastic moduli in the radial (p < 0.0001), circumferential (p < 0.001) and statistically insignificant difference in axial (p > 0.05) orientation. For the moduli ratios between groups: at periosteum, only significant difference in the Z/R (p < 0.05) anisotropy ratio, whereas at endosteum, a statistically significant difference in Z/C (p < 0.001), and Z/R (p < 0.001), as well as C/R (p < 0.05) anisotropy ratios, was evaluated. The results suggested initial bone remodeling impaired bone micro-architecture predominantly at the medullary endosteum with possible alterations in the geometric orientations of collagen and mineral phases inside the bone. The findings could be significant for studying the mechanotransduction pathways involved in maintaining the bone micro-architecture and possibly have high clinical significance for drug use against impairment from functional disuse.

Quantitative ultrasound (QUS) axial transmission method reflects anisotropy in micro-arrangement of apatite crystallites in human long bones: A study with 3-MHz-frequency ultrasound

Anisotropic arrangement of apatite crystallites, i.e., preferential orientation of the apatite c-axis, is known to be an important bone quality parameter that governs the mechanical properties. However, noninvasive evaluation of apatite orientation has not been achieved to date. The present paper reports the potential of quantitative ultrasound (QUS) for noninvasive evaluation of the degree of apatite orientation in human bone for the first time. A novel QUS instrument for implementation of the axial transmission (AT) method is developed, so as to achieve precise measurement of the speed of sound (SOS) in the cortex (cSOS) of human long bone. The advantages of our QUS instrument are the following: (i) it is equipped with a cortical bone surface-morphology detection system to correct the ultrasound transmission distance, which should be necessary for AT measurement of long bone covered by soft tissue of non-uniform thickness; and (ii) ultrasound with a relatively high frequency of 3 MHz is employed, enabling thickness-independent cSOS measurement even for the thin cortex by preventing guide wave generation. The reliability of the proposed AT measurement system is confirmed through comparison with the well-established direct transmission (DT) method. The cSOS in human long bone is found to exhibit considerable direction-dependent anisotropy; the axial cSOS (3870 ± 66 m/s) is the highest, followed by the tangential (3411 ± 94 m/s) and radial (3320 ± 85 m/s) cSOSs. The degree of apatite orientation exhibits the same order, despite the unchanged bone mineral density. Multiple regression analysis reveals that the cSOS of human long bone strongly reflects the apatite orientation. The cSOS determined by the AT method is positively correlated with that determined by the DT method and sensitively reflects the apatite orientation variation, indicating the validity of the AT instrument developed in this study. Our instrument will be beneficial for noninvasive evaluation of the material integrity of the human long-bone cortex, as determined by apatite c-axis orientation along the axial direction.

Application of the Johnson-Cook plasticity model in the finite element simulations of the nanoindentation of the cortical bone

The mechanical behavior of the cortical bone in nanoindentation is a complicated mechanical problem. The finite element analysis has commonly been assumed to be the most appropriate approach to this issue. One significant problem in nanoindentation modeling of the elastic-plastic materials is pile-up deformation, which is not observed in cortical bone nanoindentation testing. This phenomenon depends on the work-hardening of materials; it doesn't occur for work-hardening materials, which suggests that the cortical bone could be considered as a work-hardening material. Furthermore, in a recent study [59], a plastic hardening until failure was observed on the micro-scale of a dry ovine osteonal bone samples subjected to micropillar compression. The purpose of the current study was to apply an isotropic hardening model in the finite element simulations of the nanoindentation of the cortical bone to predict its mechanical behavior. The Johnson-Cook (JC) model was chosen as the constitutive model. The finite element modeling in combination with numerical optimization was used to identify the unknown material constants and then the finite element solutions were compared to the experimental results. A good agreement of the numerical curves with the target loading curves was found and no pile-up was predicted. A Design Of Experiments (DOE) approach was performed to evaluate the linear effects of the material constants on the mechanical response of the material. The strain hardening modulus and the strain hardening exponent were the most influential parameters. While a positive effect was noticed with the Young's modulus, the initial yield stress and the strain hardening modulus, an opposite effect was found with the Poisson's ratio and the strain hardening exponent. Finally, the JC model showed a good capability to describe the elastoplastic behavior of the cortical bone.

Analysis of Mechanical Properties and Mechanical Anisotropy in Canine Bone Tissues of Various Ages

The effect of age on mechanical behavior and microstructure anisotropy of bone is often ignored by researchers engaged in the study of biomechanics. The objective of our study was to determine the variations in mechanical properties of canine femoral cortical bone with age and the mechanical anisotropy between the longitudinal and transverse directions. Twelve beagles divided into three age groups (6, 12, and 36 months) were sacrificed and all femurs were extracted. The longitudinal and transverse samples of cortical bone were harvested from three regions of diaphysis (proximal, central, and distal). A nanoindentation technique was used for simultaneously measuring force and displacement of a diamond tip pressed 2000nm into the hydrated bone tissue. An elastic modulus was calculated from the unloading curve with an assumed Poisson ratio of 0.3, while hardness was defined as the maximal force divided by the corresponding contact area. The mechanical properties of cortical bone were determined from 852 indents on two orthogonal cross-sectional surfaces. Mean elastic modulus ranged from 7.56±0.32 GPa up to 21.56±2.35 GPa, while mean hardness ranged from 0.28±0.057 GPa up to 0.84±0.072 GPa. Mechanical properties of canine femoral cortical bone tended to increase with age, but the magnitudes of these increase for each region might be different. The longitudinal mechanical properties were significantly higher than that of transverse direction (P<0.01). A significant anisotropy was found in the mechanical properties while there was no significant correlation between the two orthogonal directions in each age group (r²<0.3). Beyond that, the longitudinal mechanical properties of the distal region in each age group were lower than the proximal and central regions. Hence, mechanical properties in nanostructure of bone tissue must differ mainly among age, sample direction, anatomical sites, and individuals. These results may help a number of researchers develop more accurate constitutive micromechanics models of bone tissue in future studies.

Bone fracture healing under Ilizarov fixator: Influence of fixator configuration, fracture geometry, and loading

This study aims to enhance the understanding of the relationship between Ilizarov fixator configuration and its effects on bone fracture healing. Using Taylor spatial frame (TSF) as an example, the roles of critical parameters (ie, TSF ring diameter, wire pre-tension, fracture gap size, and axial load) that govern fracture healing during the early stages were investigated by using computational modelling in conjunction with mechanical testing involving an advanced 3D optical measurement system. The computational model was first validated using the mechanical test results and then used to simulate mesenchymal stem cell (MSC) differentiations within different regions of the fracture site under various combinations of TSF ring diameter, wire pre-tension, fracture gap size, and axial load values. Predicted spatially dependent MSC differentiation patterns and the influence of each parameter on differentiations were compared with in vivo results, and good agreement was seen between the two. Gap size was identified as the most influential parameter in MSC differentiation, and the influence of axial loading and TSF configuration (ie, ring diameter and wire pre-tension) on cell differentiation was seen to be gap size dependent. Most changes in cell differentiation were predicted in the external callus (periosteal), which is the crucial region of the callus in the early stages. However, for small gap sizes (eg, 1 mm), significant changes were predicted in the endosteal callus as well. The study exhibits the potential of computational models in assessing the performance of Ilizarov fixators as well as assisting surgeons in patient-specific clinical treatment planning.

Mechanical Stimuli in the Local In Vivo Environment in Bone: Computational Approaches Linking Organ-Scale Loads to Cellular Signals

PURPOSE OF REVIEW

Connecting organ-scale loads to cellular signals in their local in vivo environment is a current challenge in the field of bone (re)modelling. Understanding this critical missing link would greatly improve our ability to anticipate mechanotransduction during different modes of stimuli and the resultant cellular responses. This review characterises computational approaches that could enable coupling links across the multiple scales of bone.

RECENT FINDINGS

Current approaches using strain and fluid shear stress concepts have begun to link organ-scale loads to cellular signals; however, these approaches fail to capture localised micro-structural heterogeneities. Furthermore, models that incorporate downstream communication from osteocytes to osteoclasts, bone-lining cells and osteoblasts, will help improve the understanding of (re)modelling activities. Incorporating this potentially key information in the local in vivo environment will aid in developing multiscale models of mechanotransduction that can predict or help describe resultant biological events related to bone (re)modelling. Progress towards multiscale determination of the cell mechanical environment from organ-scale loads remains elusive. Construction of organ-, tissue- and cell-scale computational models that include localised environmental variation, strain amplification and intercellular communication mechanisms will ultimately help couple the hierarchal levels of bone.

Mechanical Characterization of Bone: State of the Art in Experimental Approaches-What Types of Experiments Do People Do and How Does One Interpret the Results?

PURPOSE OF REVIEW

The mechanical integrity of bone is determined by the direct measurement of bone mechanical properties. This article presents an overview of the current, most common, and new and upcoming experimental approaches for the mechanical characterization of bone. The key outcome variables of mechanical testing, as well as interpretations of the results in the context of bone structure and biology are also discussed.

RECENT FINDINGS

Quasi-static tests are the most commonly used for determining the resistance to structural failure by a single load at the organ (whole bone) level. The resistance to crack initiation or growth by fracture toughness testing and fatigue loading offers additional and more direct characterization of tissue material properties. Non-traditional indentation techniques and in situ testing are being increasingly used to probe the material properties of bone ultrastructure. Destructive ex vivo testing or clinical surrogate measures are considered to be the gold standard for estimating fracture risk. The type of mechanical test used for a particular investigation depends on the length scale of interest, where the outcome variables are influenced by the interrelationship between bone structure and composition. Advancement in the sensitivity of mechanical characterization techniques to detect changes in bone at the levels subjected to modifications by aging, disease, and/or pharmaceutical treatment is required. As such, a number of techniques are now available to aid our understanding of the factors that contribute to fracture risk.

Modal analysis of nanoindentation data, confirming that reduced bone turnover may cause increased tissue mineralization/elasticity

It is widely believed that the activities of bone cells at the tissue scale not only govern the size of the vascular pore spaces (and hence, the amount of bone tissue available for actually carrying the loads), but also the characteristics of the extracellular bone matrix itself. In this context, increased mechanical stimulation (in mediolateral regions of human femora, as compared to anteroposterior regions) may lead to increased bone turnover, lower bone matrix mineralization, and therefore lower tissue modulus. On the other hand, resorption-only processes (in endosteal versus periosteal regions) may have the opposite effect. A modal analysis of nanoindentation data obtained on femurs from the Melbourne Femur Research Collection (MFRC) indeed confirms that bone is stiffer in endosteal regions compared to periosteal regions (E̅_endost = 29.34 ± 0.75 GPa >E̅_periost = 24.67 ± 1.63 GPa), most likely due to the aging-related increase in resorption modeling on endosteal surfaces resulting in trabecularization of cortical bone. The results also show that bone is stiffer along the anteroposterior direction compared the mediolateral direction (E̅_anteropost = 28.89 ± 1.08 GPa >E̅_mediolat = 26.03 ± 2.31 GPa), the former being aligned with the neutral bending axis of the femur and, thus, undergoing more resorption modeling and consequently being more mineralized.

Determination of a tissue-level failure evaluation standard for rat femoral cortical bone utilizing a hybrid computational-experimental method

Macro-level failure in bone structure could be diagnosed by pain or physical examination. However, diagnosing tissue-level failure in a timely manner is challenging due to the difficulty in observing the interior mechanical environment of bone tissue. Because most fractures begin with tissue-level failure in bone tissue caused by continually applied loading, people attempt to monitor the tissue-level failure of bone and provide corresponding measures to prevent fracture. Many tissue-level mechanical parameters of bone could be predicted or measured; however, the value of the parameter may vary among different specimens belonging to a kind of bone structure even at the same age and anatomical site. These variations cause difficulty in representing tissue-level bone failure. Therefore, determining an appropriate tissue-level failure evaluation standard is necessary to represent tissue-level bone failure. In this study, the yield and failure processes of rat femoral cortical bones were primarily simulated through a hybrid computational-experimental method. Subsequently, the tissue-level strains and the ratio between tissue-level failure and yield strains in cortical bones were predicted. The results indicated that certain differences existed in tissue-level strains; however, slight variations in the ratio were observed among different cortical bones. Therefore, the ratio between tissue-level failure and yield strains for a kind of bone structure could be determined. This ratio may then be regarded as an appropriate tissue-level failure evaluation standard to represent the mechanical status of bone tissue.

Selective effect of hydroxyapatite nanoparticles on osteoporotic and healthy bone formation correlates with intracellular calcium homeostasis regulation

Adequate bone substitutes osseointegration has been difficult to achieve in osteoporosis. Hydroxyapatite of the osteoporotic bone, secreted by pathologic osteoblasts, had a smaller crystal size and lower crystallinity than that of the normal. To date, little is known regarding the interaction of synthetic hydroxyapatite nanoparticles (HANPs) with osteoblasts born in bone rarefaction. The present study investigated the biological effects of HANPs on osteoblastic cells derived from osteoporotic rat bone (OVX-OB), in comparison with the healthy ones (SHM-OB). A selective effect of different concentrations of HANPs on the two cell lines was observed that the osteoporotic osteoblasts had a higher tolerance. Reductions in cell proliferation, ALP activity, collagen secretion and osteoblastic gene expressions were found in the SHM-OB when administered with HANPs concentration higher than 25µg/ml. In contrast, those of the OVX-OB suffered no depression but benefited from 25 to 250µg/ml HANPs in a dose-dependent manner. We demonstrated that the different effects of HANPs on osteoblasts were associated with the intracellular calcium influx into the endoplasmic reticulum. The in vivo bone defect model further confirmed that, with a critical HANPs concentration administration, the osteoporotic rats had more and mechanically matured new bone formation than the non-treated ones, whilst the sham rats healed no better than the natural healing control. Collectively, the observed epigenetic regulation of osteoblastic cell function by HANPs has significant implication on defining design parameters for a potential therapeutic use of nanomaterials.

STATEMENT OF SIGNIFICANCE

In this study, we investigated the biological effects of hydroxyapatite nanoparticles (HANPs) on osteoporotic rat bone and the derived osteoblast. Our findings revealed a previously unrecognized phenomenon that the osteoporotic individuals could benefit from higher concentrations of HANPs, as compared with the healthy individuals. The in vivo bone defect model confirmed that, with a critical HANPs concentration administration, the osteoporotic rats had more mechanically matured new bone formation than the non-treated ones, whilst the sham rats healed no better than the natural healing control. The selective effect of HANPs might be associated with the intracellular calcium influx into the endoplasmic reticulum. Collectively, the observed epigenetic regulation by HANPs has significant implication on defining design parameters for a potential therapeutic use of nanomaterials in a pathological condition.

Micromechanical anisotropy and heterogeneity of the meniscus extracellular matrix

To understand how the complex biomechanical functions of the meniscus are endowed by the nanostructure of its extracellular matrix (ECM), we studied the anisotropy and heterogeneity in the micromechanical properties of the meniscus ECM. We used atomic force microscopy (AFM) to quantify the time-dependent mechanical properties of juvenile bovine meniscus at deformation length scales corresponding to the diameters of collagen fibrils. At this scale, anisotropy in the elastic modulus of the circumferential fibers, the major ECM structural unit, can be attributed to differences in fibril deformation modes: uncrimping when normal to the fiber axis, and laterally constrained compression when parallel to the fiber axis. Heterogeneity among different structural units is mainly associated with their variations in microscale fiber orientation, while heterogeneity across anatomical zones is due to alterations in collagen fibril diameter and alignment at the nanoscale. Unlike the elastic modulus, the time-dependent properties are more homogeneous and isotropic throughout the ECM. These results enable a detailed understanding of the meniscus structure-mechanics at the nanoscale, and can serve as a benchmark for understanding meniscus biomechanical functions, documenting disease progression and designing tissue repair strategies.

STATEMENT OF SIGNIFICANCE

Meniscal damage is a common cause of joint injury, which can lead to the development of post-traumatic osteoarthritis among young adults. Restoration of meniscus function requires repairing its highly heterogeneous and complex extracellular matrix. Employing AFM, this study quantifies the anisotropic and heterogeneous features of the meniscus ECM structure and mechanics. The micromechanical properties are interpreted within the context of the collagen fibril nanostructure and its variation with tissue anatomical locations. These results provide a fundamental structure-mechanics knowledge benchmark, against which, repair and regeneration strategies can be developed and evaluated with respect to the specialized structural and functional complexity of the native tissue.

Nanoindentation analysis of the micromechanical anisotropy in mouse cortical bone

Studies investigating micromechanical properties in mouse cortical bone often solely focus on the mechanical behaviour along the long axis of the bone. Therefore, data on the anisotropy of mouse cortical bone is scarce. The aim of this study is the first-time evaluation of the anisotropy ratio between the longitudinal and transverse directions of reduced modulus and hardness in mouse femurs by using the nanoindentation technique. For this purpose, nine 22-week-old mice (C57BL/6) were sacrificed and all femurs extracted. A total of 648 indentations were performed with a Berkovich tip in the proximal (P), central (C) and distal (D) regions of the femoral shaft in the longitudinal and transverse directions. Higher values for reduced modulus are obtained for indentations in the longitudinal direction, with anisotropy ratios of 1.72 ± 0.40 (P), 1.75 ± 0.69 (C) and 1.34 ± 0.30 (D). Hardness is also higher in the longitudinal direction, with anisotropic ratios of 1.35 ± 0.27 (P), 1.35 ± 0.47 (C) and 1.17 ± 0.19 (D). We observed a significant anisotropy in the micromechanical properties of the mouse femur, but the correlation for reduced modulus and hardness between the two directions is low (r² < 0.3) and not significant. Therefore, we highly recommend performing independent indentation testing in both the longitudinal and transverse directions when knowledge of the tissue mechanical behaviour along multiple directions is required.

Nanoindentation analysis of the micromechanical anisotropy in mouse cortical bone

Studies investigating micromechanical properties in mouse cortical bone often solely focus on the mechanical behaviour along the long axis of the bone. Therefore, data on the anisotropy of mouse cortical bone is scarce. The aim of this study is the first-time evaluation of the anisotropy ratio between the longitudinal and transverse directions of reduced modulus and hardness in mouse femurs by using the nanoindentation technique. For this purpose, nine 22-week-old mice (C57BL/6) were sacrificed and all femurs extracted. A total of 648 indentations were performed with a Berkovich tip in the proximal (P), central (C) and distal (D) regions of the femoral shaft in the longitudinal and transverse directions. Higher values for reduced modulus are obtained for indentations in the longitudinal direction, with anisotropy ratios of 1.72 ± 0.40 (P), 1.75 ± 0.69 (C) and 1.34 ± 0.30 (D). Hardness is also higher in the longitudinal direction, with anisotropic ratios of 1.35 ± 0.27 (P), 1.35 ± 0.47 (C) and 1.17 ± 0.19 (D). We observed a significant anisotropy in the micromechanical properties of the mouse femur, but the correlation for reduced modulus and hardness between the two directions is low (r² < 0.3) and not significant. Therefore, we highly recommend performing independent indentation testing in both the longitudinal and transverse directions when knowledge of the tissue mechanical behaviour along multiple directions is required.

Microscopic assessment of bone toughness using scratch tests

Bone is a composite material with five distinct structural levels: collagen molecules, mineralized collagen fibrils, lamellae, osteon and whole bone. However, most fracture testing methods have been limited to the macroscopic scale and there is a need for advanced characterization methods to assess toughness at the osteon level and below. The goal of this investigation is to present a novel framework to measure the fracture properties of bone at the microscopic scale using scratch testing. A rigorous experimental protocol is articulated and applied to examine cortical bone specimens from porcine femurs. The observed fracture behavior is very complex: we observe a strong anisotropy of the response with toughening mechanisms and a competition between plastic flow and brittle fracture. The challenge consists then in applying nonlinear fracture mechanics methods such as the J-integral or the energetic Size Effect Law to quantify the fracture toughness in a rigorous fashion. Our result suggests that mixed-mode fracture is instrumental in determining the fracture resistance. There is also a pronounced coupling between fracture and elasticity. Our methodology opens the door to fracture assessment at multiple structural levels, microscopic and potentially nanometer length scale, due to the scalability of scratch tests.

•

We measure bone toughness at the osteon level using microscopic scratch tests.

•

Scratch tests induce a ductile-to-brittle transition in cortical bone.

•

Toughening mechanisms are observed at the microscopic scale.

•

The J-integral and the size effect law are employed to predict the fracture energy.

Spatial distribution and remodeling of elastic modulus of bone in micro-regime as prediction of early stage osteoporosis

We assessed the local distribution of bone mechanical properties on a micro-nano-scale and its correlation to strain distribution. Left tibia samples were obtained from 5-month old female Sprague Dawley rats, including baseline control (n=9) and hindlimb suspended (n=9) groups. Elastic modulus was measured by nanoindentation at the dedicated locations. Three additional tibias from control rats were loaded axially to measure bone strain, with 6-10N at 1Hz on a Bose machine for strain measurements. In the control group, the difference of the elastic modulus between periosteum and endosteum was much higher at the anterior and posterior regions (2.6GPa), where higher strain differences were observed (45μɛ). Minimal elastic modulus difference between periosteum and endosteum was observed at the medial region (0.2GPa), where neutral axis of the strain distribution was oriented with lower strain difference (5μɛ). In the disuse group, however, the elastic modulus differences in the anterior posterior regions reduced to 1.2GPa from 2.6GPa in the control group, and increased in the medial region to 2.7GPa from 0.2GPa. It is suggested that the remodeling rate in a region of bone is possibly influenced by the strain gradient from periosteum to endosteum. Such pattern of moduli gradients was compromised in disuse osteopenia, suggesting that the remodeling in distribution of micro-nano-elastic moduli among different regions may serve as a predictor for early stage of osteoporosis.

Synthesis and microstructural properties of the scaffold based on a 3-(trimethoxysilyl)propyl methacrylate–POSS hybrid towards potential tissue engineering applications

The biomimetic organic–inorganic scaffold with the chemical composition, structural dimensions, topography, and microstructural properties that fulfills the requirements for hard-tissue engineering was developed.

A Comparative Study of the Mechanical Properties of Hierarchical Trabecular Bone with other Approaches and Existing Experimental Data

The bone is a hierarchically structured material with mechanical properties depending on its architecture at all scales. Water plays an important role in the bio-mineralization process and serves as a plasticizer, enhancing the toughness of bone. In this paper, a trabecular bone multiscale model based on finite element analysis was developed to link scales from sub-nanoscopic scale (Microfibril) to sub-microscopic (Lamella) in order to predict the orthotropic properties of bone at different structural level. To identify the orthotropic properties, an inverse identification algorithm is used. Furthermore, the effect of water is incorporated. Good agreement is found between theoretical and experimental results.

Shear deformation and fracture of human cortical bone

Bone can be viewed as a nano-fibrous composite with complex hierarchical structures. Its deformation and fracture behaviors depend on both the local structure and the type of stress applied. In contrast to the extensive studies on bone fracture under compression and tension, there is a lack of knowledge on the fracture process under shear, a stress state often exists in hip fracture. This study investigated the mechanical behavior of human cortical bone under shear, with the focus on the relation between the fracture pattern and the microstructure. Iosipescu shear tests were performed on notched rectangular bar specimens made from human cortical bone. They were prepared at different angles (i.e. 0°, 30°, 60° and 90°) with respect to the long axis of the femoral shaft. The results showed that human cortical bone behaved as an anisotropic material under shear with the highest shear strength (~50MPa) obtained when shearing perpendicular to the Haversian systems or secondary osteons. Digital image correlation (DIC) analysis found that shear strain concentration bands had a close association with long bone axis with an average deviation of 11.8° to 18.5°. The fracture pattern was also greatly affected by the structure with the crack path generally following the direction of the long axes of osteons. More importantly, we observed unique peripheral arc-shaped microcracks within osteons, using laser scanning confocal microscopy (LSCM). They were generally long cracks that developed within a lamella without crossing the boundaries. This microcracking pattern clearly differed from that created under either compressive or tensile stress: these arc-shaped microcracks tended to be located away from the Haversian canals in early-stage damaged osteons, with ~70% developing in the outer third osteonal wall. Further study by second harmonic generation (SHG) and two-photon excitation fluorescence (TPEF) microscopy revealed a strong influence of the organization of collagen fibrils on shear microcracking. This study concluded that shear-induced microcracking of human cortical bone follows a unique pattern that is governed by the lamellar structure of the osteons.

Biomechanics and mechanobiology of trabecular bone: a review

Trabecular bone is a highly porous, heterogeneous, and anisotropic material which can be found at the epiphyses of long bones and in the vertebral bodies. Studying the mechanical properties of trabecular bone is important, since trabecular bone is the main load bearing bone in vertebral bodies and also transfers the load from joints to the compact bone of the cortex of long bones. This review article highlights the high dependency of the mechanical properties of trabecular bone on species, age, anatomic site, loading direction, and size of the sample under consideration. In recent years, high resolution micro finite element methods have been extensively used to specifically address the mechanical properties of the trabecular bone and provide unique tools to interpret and model the mechanical testing experiments. The aims of the current work are to first review the mechanobiology of trabecular bone and then present classical and new approaches for modeling and analyzing the trabecular bone microstructure and macrostructure and corresponding mechanical properties such as elastic properties and strength.

To what extent can cortical bone millimeter-scale elasticity be predicted by a two-phase composite model with variable porosity?

An evidence gap exists in fully understanding and reliably modeling the variations in elastic anisotropy that are observed at the millimeter scale in human cortical bone. The porosity (pore volume fraction) is known to account for a large part, but not all, of the elasticity variations. This effect may be modeled by a two-phase micromechanical model consisting of a homogeneous matrix pervaded by cylindrical pores. Although this model has been widely used, it lacks experimental validation. The aim of the present work is to revisit experimental data (elastic coefficients, porosity) previously obtained from 21 cortical bone specimens from the femoral mid-diaphysis of 10 donors and test the validity of the model by proposing a detailed discussion of its hypotheses. This includes investigating to what extent the experimental uncertainties, pore network modeling, and matrix elastic properties influence the model's predictions. The results support the validity of the two-phase model of cortical bone which assumes that the essential source of variations of elastic properties at the millimeter-scale is the volume fraction of vascular porosity. We propose that the bulk of the remaining discrepancies between predicted stiffness coefficients and experimental data (RMSE between 6% and 9%) is in part due to experimental errors and part due to small variations of the extravascular matrix properties. More significantly, although most of the models that have been proposed for cortical bone were based on several homogenization steps and a large number of variable parameters, we show that a model with a single parameter, namely the volume fraction of vascular porosity, is a suitable representation for cortical bone. The results could provide a guide to build specimen-specific cortical bone models. This will be of interest to analyze the structure-function relationship in bone and to design bone-mimicking materials.