Analyzing the anisotropic Hooke׳s law for children׳s cortical bone

Immature porcine cortical bone mechanical properties and composition change with maturation and displacement rate

Computational models of mature bone have been used to predict fracture; however, analogous study of immature diaphyseal fracture has not been conducted due to sparse experimental mechanical data. A model of immature bone fracture may be used to aid in the differentiation of accidental and non-accidental trauma fractures in young, newly ambulatory children (0-3 years). The objective of this study was to characterize the evolution of tissue-level mechanical behavior, composition, and microstructure of maturing cortical porcine bone with uniaxial tension, Raman spectroscopy, and light microscopy as a function of maturation. We asked: 1) How do the monotonic uniaxial tensile properties change with maturation and displacement rate; 2) How does the composition and microstructure change with maturation; and 3) Is there a correlation between composition and tensile properties with maturation? Elastic modulus (p < 0.001), fracture stress (p < 0.001), and energy absorption (p < 0.014) increased as a function of maturation at the quasistatic rate by 110%, 86%, and 96%, respectively. Fracture stress also increased by 90% with maturation at the faster rate (p = 0.001). Fracture stress increased as a function of increasing displacement rate by 28% (newborn p = 0.048; 1-month p = 0.004; 3-month p= < 0.001), and fracture strain decreased by 68% with increasing displacement rate (newborn p = 0.002; 1-month p = 0.036; 3-month p < 0.001). Carbonate-to-phosphate ratio was positively linearly related to elastic modulus, and fracture stress was positively related to carbonate-to-phosphate ratio and matrix maturation ratio. The results of this study support that immature bone is strain-rate dependent and becomes more brittle at faster rates, contributing to the foundation upon which a computational model can be built to evaluate immature bone fracture.

Variation in Juvenile Long Bone Properties as a Function of Age: Mechanical and Compositional Characterization

Bone is a heterogeneous, hierarchical biocomposite material made of an organic matrix filled with a mineral component, which plays an important role in bone strength. Although the effect of the mineral/matrix ratio on the mechanical properties of bone during aging has been intensively investigated, the relationship between the mechanical properties and the chemical composition of bone with age requires additional research in juvenile individuals. In this study, bone coupons from bovine and ovine animal species were machined from cortical areas of long bones to quantify whether the variation in mechanical properties at different stages of development is related to the change in the composition of bone tissue. An energy-dispersive X-ray detector (EDX) attached to a scanning electron microscope (SEM) was used to perform a compositional analysis of the tissue. In addition, nanoindentation analyses were carried out to address how the elastic modulus changed with age. Nonparametric statistical analyses found significant differences (p < 0.05) in Ca content and elastic modulus between species, but no differences were found within each species with development. A multiple linear regression model found that the elastic modulus was significantly related to the decrease in P and C in the samples, to the animal species (larger in bovine), and development, although not linearly. This model also found an interaction between Ca and development that could explain the lack of significance of the relationship between the elastic modulus and development in the univariate models.

Towards mechanobiologically optimized mandible reconstruction: CAD/CAM miniplates vs. reconstruction plates for fibula free flap fixation: A finite element study

Due to their advantages in applicability, patient-specific (CAD/CAM) reconstruction plates are increasingly used in fibula free flap mandible reconstruction. In addition, recently, CAD/CAM miniplates, with further advantages in postoperative management, have been introduced. However, biomechanical conditions induced by CAD/CAM systems remain partially unknown. This study aimed to evaluate the primary fixation stability of CAD/CAM fixators. For a patient-specific scenario, the biomechanical conditions induced in a one segmental fibula free flap stabilized using either a CAD/CAM reconstruction plate or CAD/CAM miniplates were determined using finite element analysis. The main output parameters were the strains between intersegmental bone surfaces and stresses in the fixation systems due to different biting scenarios. CAD/CAM miniplates resulted in higher mechanical strains in the mesial interosseous gap, whereas CAD/CAM reconstruction plate fixation resulted in higher strains in the distal interosseous gap. For all investigated fixation systems, stresses in the fixation systems were below the material yield stress and thus material failure would not be expected. While the use of CAD/CAM miniplates resulted in strain values considered adequate to promote bone healing in the mesial interosseous gap, in the distal interosseous gap CAD/CAM reconstruction plate fixation might result in more beneficial tissue straining. A mechanical failure of the fixation systems would not be expected.

Biomechanics of immature human cortical bone: A systematic review

The whole bone geometry, microstructure, and mechanical properties of mature human bone are widely reported; however, immature bone (0-18 years) has not been similarly robustly characterized. There is an interest in analyzing and predicting the mechanical loading conditions associated with long bone diaphyseal fractures attributed to trauma in children. Thus, understanding the mechanical properties of immature bone in a temporal reference frame is an essential first step to understand diaphyseal fractures of pediatric long bones. The purpose of this systematic review was to ask, what is the state of knowledge regarding the 1) evolution of whole bone geometry and microstructure of immature pediatric bone as a function of maturation and 2) cortical bone density and experimental quasi-static mechanical properties at the tissue level in the diaphyseal region of immature pediatric long bones? The systematic search yielded 36 studies of the whole bone geometry, microstructure, and mechanical properties of immature pediatric long bones. The elastic modulus, yield stress, and ultimate stress were shown to generally increase with maturation, whereas the yield strain was approximately invariant; however, the specific year-to-year progression of these properties could not be characterized from the limited studies available. The results of this systematic search indicate there is a dearth of knowledge associated with the biomechanics of cortical bone from immature pediatric long bones; it also provides a basis for computational studies of immature human long bones. Additional biomechanical studies of immature human bone are necessary to develop a robust catalogue, which can be used in broad applications to understand fracture mechanics, bone pathologies, and athletic injury in the pediatric setting.

Ultrasounds could be considered as a future tool for probing growing bone properties

Juvenile bone growth is well described (physiological and anatomical) but there are still lacks of knowledge on intrinsic material properties. Our group has already published, on different samples, several studies on the assessment of intrinsic material properties of juvenile bone compared to material properties of adult bone. The purpose of this study was finally to combine different experimental modalities available (ultrasonic measurement, micro-Computed Tomography analysis, mechanical compression tests and biochemical measurements) applied on small cubic bone samples in order to gain insight into the multiparametric evaluation of bone quality. Differences were found between juvenile and adult groups in term of architectural parameters (Porosity Separation), Tissue Mineral Density (TMD), diagonal stiffness coefficients (C₃₃, C_44, C_55, C₆₆) and ratio between immature and mature cross-links (CX). Diagonal stiffness coefficients are more representative of the microstructural and biochemical parameters of child bone than of adult bone. We also found that compression modulus E was highly correlated with several microstructure parameters and CX in children group while it was not at all correlated in the adult group. Similar results were found for the CX which was linked to several microstructure parameters (TMD and E) only in the juvenile group. To our knowledge, this is the first time that, on a same sample, ultrasonic measurements have been combined with the assessment of mechanical and biochemical properties. It appears that ultrasonic measurements can provide relevant indicators of child bone quality (microstructural and biochemical parameters) which is promising for clinical application since, B-mode ultrasound is the preferred first-line modality over other more constraining imaging modalities (radiation, parent–child accessibility and access to the patient's bed) for pediatric patients.

Compositional and mechanical properties of growing cortical bone tissue: a study of the human fibula

Human cortical bone contains two types of tissue: osteonal and interstitial tissue. Growing bone is not well-known in terms of its intrinsic material properties. To date, distinctions between the mechanical properties of osteonal and interstitial regions have not been investigated in juvenile bone and compared to adult bone in a combined dataset. In this work, cortical bone samples obtained from fibulae of 13 juveniles patients (4 to 18 years old) during corrective surgery and from 17 adult donors (50 to 95 years old) were analyzed. Microindentation was used to assess the mechanical properties of the extracellular matrix, quantitative microradiography was used to measure the degree of bone mineralization (DMB), and Fourier transform infrared microspectroscopy was used to evaluate the physicochemical modifications of bone composition (organic versus mineral matrix). Juvenile and adult osteonal and interstitial regions were analyzed for DMB, crystallinity, mineral to organic matrix ratio, mineral maturity, collagen maturity, carbonation, indentation modulus, indicators of yield strain and tissue ductility using a mixed model. We found that the intrinsic properties of the juvenile bone were not all inferior to those of the adult bone. Mechanical properties were also differently explained in juvenile and adult groups. The study shows that different intrinsic properties should be used in case of juvenile bone investigation.

Estimation of the elastic modulus of child cortical bone specimens via microindentation

Purpose: Non-pathological child cortical bone (NPCCB) studies can provide clinicians with vital information and insights. However, assessing the anisotropic elastic properties of NPCCB remains a challenge for the biomechanical engineering community. For the first time, this paper provides elastic moduli values for NPCCB specimens in two perpendicular directions (longitudinal and transverse) and for two different structural components of bone tissue (osteon and interstitial lamellae). Materials and Methods: Microindentation is one of the reference methods used to measure bone stiffness. Here, 8 adult femurs (mean age 82 ± 8.9 years), 3 child femurs (mean age 13.3 ± 2.1 years), and 16 child fibulae (mean age 10.2 ± 3.9 years) were used to assess the elastic moduli of adult and child bones by microindentation. Results: For adult specimens, the mean moduli measured in this study are 18.1 (2.6) GPa for osteons, 21.3 (2.3) GPa for interstitial lamellae, and 13.8 (1.7) GPa in the transverse direction. For child femur specimens, the mean modulus is 14.1 (0.8) GPa for osteons, lower than that for interstitial lamellae: 15.5 (1.5) GPa. The mean modulus is 11.8 (0.7) GPa in the transverse direction. Child fibula specimens show a higher elastic modulus for interstitial lamellae 15.8 (1.5) than for osteons 13.5 (1.6), with 10.2 (1) GPa in the transverse direction. Conclusion: For the first time, NPCCB elastic modulus values are provided in longitudinal and transverse directions at the microscale level.

Anisotropic elastic properties of human femoral cortical bone and relationships with composition and microstructure in elderly

The strong dependence between cortical bone elasticity at the millimetre-scale (mesoscale) and cortical porosity has been evidenced by previous studies. However, bone is an anisotropic composite material made by mineral, proteins and water assembled in a hierarchical structure. Whether the variations of structural and compositional properties of bone affect the different elastic coefficients at the mesoscale is not clear. Aiming to understand the relationships between bone elastic properties and compositions and microstructure, we applied state-of-the-art experimental modalities to assess these aspects of bone characteristics. All elastic coefficients (stiffness tensor of the transverse isotropic bone material), structure of the vascular pore network, collagen and mineral properties were measured in 52 specimens from the femoral diaphysis of 26 elderly donors. Statistical analyses and micromechanical modeling showed that vascular pore volume fraction and the degree of mineralization of bone are the most important determinants of cortical bone anisotropic mesoscopic elasticity. Though significant correlations were observed between collagen properties and elasticity, their effects in bone mesoscopic elasticity were minor in our data. This work also provides a unique set of data exhibiting a range of variations of compositional and microstructural cortical bone properties in the elderly and gives strong experimental evidence and basis for further development of biomechanical models for human cortical bone. STATEMENT OF SIGNIFICANCE: This study reports the relationships between microstructure, composition and the mesoscale anisotropic elastic properties of human femoral cortical bone in elderly. For the first time, we provide data covering the complete anisotropic elastic tensor, the microstructure of cortical vascular porosity, mineral and collagen characteristics obtained from the same or adjacent samples in each donor. The results revealed that cortical vascular porosity and degree of mineralization of bone are the most important determinants of bone anisotropic stiffness at the mesoscale. The presented data gives strong experimental evidence and basis for further development of biomechanical models for human cortical bone.

Optimization of Taylor spatial frame half-pins diameter for bone deformity correction: Application to femur

Using external fixtures for bone deformity correction takes advantages of less soft tissue injury, better bone alignment and enhances strain development for bone formation on cutting section, which cause shorter healing time. Among these fixtures, Taylor spatial frame is widely used and includes two rings and six adjustable struts developing 6 degrees of freedom, making them very flexible for this type of application. The current study describes a method to optimize Taylor spatial frame pin-sizes currently chosen from the surgeon's experiences. A three-dimensional model of femur was created from computed tomography images; segmentation of the medical images was made based on the Hounsfield unit (gray scale) in order to allocate adequate mechanical properties into cortical and trabecular bone sections. Both the cortical and trabecular sections were assumed to be isotropic and homogeneous. The diameter optimization of Taylor spatial frame's half-pins was carried out by coupling genetic algorithm and finite element analysis. The finite element analysis was based on a static mechanical load corresponding to a standing person's body weight. Finite element analysis results were validated with experimentally measured strains obtained from bone compression tests. A cost function, based on the developed bone stresses, was defined close to the Taylor spatial frame's half-pins. The calculated cost function showed a decrease of over 33% from the initial half-pin selection by the surgeon and the genetic algorithm optimization. Consequently, the maximum stresses experienced by the bone in the connected location of the half-pins decreased from 121.4 MPa in the surgeon's selection to 73.07 MPa as a result of the optimization process.

Assessment of elastic coefficients of child cortical bone using resonant ultrasound spectroscopy

The assessment of the anisotropic elastic properties of non-pathological child cortical bone remains a challenge for the biomechanical engineering community and an important clinical issue. Resonant ultrasound spectroscopy (RUS) can be used to determine bone stiffness coefficients from the mechanical resonances of bone specimens. Here, a RUS protocol was used on 7 fibula specimens from children (mean age 14 ± 3 years) to estimate the whole elastic stiffness tensor of non-pathological child cortical bone considered as orthotropic. Despite a small number of sample, results are consistent with this hypothesis, even if a trend towards transverse isotropy is discussed. Indeed, the average values of the 9 independent stiffness coefficients obtained in this study for child bone are: C₁₁ = 16.73 ± 0.19 GPa, C₂₂ = 16.19 ± 0.12 GPa, C₃₃ = 24.47 ± 0.30 GPa, C₄₄ = 4.14 ± 0.08 GPa, C₅₅ = 4.16 ± 0.07 GPa, C₆₆ = 3.13 ± 0.05 GPa, C₁₂ = 10.14 ± 0.20 GPa, C₁₃ = 10.67 ± 0.27 GPa, C₂₃ = 10.25 ± 0.14 GPa.

A Review of Pediatric Lower Extremity Data for Pedestrian Numerical Modeling: Injury Epidemiology, Anatomy, Anthropometry, Structural, and Mechanical Properties

Pedestrian injuries are the fourth leading cause of unintentional injury-related death among children aged 1 to 19. The lower extremity represents the most frequently injured body region in car-to-pedestrian accidents. The goal of this study was to perform a systematic review of the data related to pedestrian lower extremity injuries, anatomy, anthropometry, structural, and mechanical properties, which can be used in the development of new pediatric computational models. The study began with a review of epidemiologic data related to pediatric pedestrian accidents. Anatomy of the child lower extremity and age-related anthropometry data were presented as well. Then, both the mechanical and structural properties of the lower extremity main components (e.g., bones, cartilages, knee ligaments, muscles, tendons, and growth plates) available in literature were summarized. The study concluded with a brief description of current child pedestrian models, which included a discussion about their limitations. We believe that data included in this review study can help in improving the biofidelity of current child models and support the development and validation of new child models used by safety researchers for protection of pediatric population.

A Novel Multiscale Mathematical Model for Building Bone Substitute Materials for Children

Bone is an engineering marvel that achieves a unique combination of stiffness and toughness exceeding that of synthesized materials. In orthopedics, we are currently challenged for the child population that needs a less stiff but a tougher bone substitute than adults. Recent evidence suggests that the relationship between inter-molecular connections that involve the two main bone building blocks, TropoCollagen molecules (TC) and carbonated Hydroxyapatite (cAp), and bone macroscopic mechanical properties, stiffness and toughness, are key to building bone substitute materials for children. The goal of our study is to establish how inter-molecular connections that occur during bone mineralization are related to macroscopic mechanical properties in child bones. Our aim is to link the biological alterations of the TC-cAp self assembly process happening during bone mineralization to the bone macroscopic mechanical properties' alterations during aging. To do so, we have developed a multiscale mathematical model that includes collagen cross links (TC⁻TC interface) from experimental studies of bone samples to forecast bone macroscopic mechanical properties. Our results support that the Young's modulus cannot be a linear parameter if we want to solve our system. In relation to bone substitute material with innovative properties for children, our results propose values of several biological parameters, such as the number of crystals and their size, and collagen crosslink maturity for the desired bone mechanical competence. Our novel mathematical model combines mineralization and macroscopic mechanical behavior of bone and is a step forward in building mechanically customized biomimetic bone grafts that would fit children's orthopedic needs.

The different distribution of enzymatic collagen cross-links found in adult and children bone result in different mechanical behavior of collagen

Enzymatic collagen cross-linking has been shown to play an important role in the macroscopic elastic and plastic deformation of bone across ages. However, its direct contribution to collagen fibril deformation is unknown. The aim of this study is to determine how covalent intermolecular connections from enzymatic collagen cross-links contribute to collagen fibril elastic and plastic deformation of adults and children's bone matrix. We used ex vivo data previously obtained from biochemical analysis of children and adults bone samples (n = 14; n = 8, respectively) to create 22 sample-specific computational models of cross-linked collagen fibrils. By simulating a tensile test for each fibril, we computed the modulus of elasticity (E), ultimate tensile and yield stress (σ_u and σ_y), and elastic, plastic and total work (W_e, W_p and W_tot) for each collagen fibril. We present a novel difference between children and adult bone in the deformation of the collagen phase and suggest a link between collagen fibril scale and macroscale for elastic behavior in children bone under the influence of immature enzymatic cross-links. We show a parametric linear correlation between W_e and immature enzymatic collagen cross-links at the collagen fibril scale in the children population that is similar to the one we found at the macroscale in our previous study. Finally, we suggest the key role of covalent intermolecular connections to stiffness parameters (e.g. elastic modulus and W_e) in children's collagen fibril and to toughness parameters in adult's collagen fibril, respectively.

Quantification of stiffness measurement errors in resonant ultrasound spectroscopy of human cortical bone

Resonant ultrasound spectroscopy (RUS) is the state-of-the-art method used to investigate the elastic properties of anisotropic solids. Recently, RUS was applied to measure human cortical bone, an anisotropic material with low Q-factor (20), which is challenging due to the difficulty in retrieving resonant frequencies. Determining the precision of the estimated stiffness constants is not straightforward because RUS is an indirect method involving minimizing the distance between measured and calculated resonant frequencies using a model. This work was motivated by the need to quantify the errors on stiffness constants due to different error sources in RUS, including uncertainties on the resonant frequencies and specimen dimensions and imperfect rectangular parallelepiped (RP) specimen geometry. The errors were first investigated using Monte Carlo simulations with typical uncertainty values of experimentally measured resonant frequencies and dimensions assuming a perfect RP geometry. Second, the exact specimen geometry of a set of bone specimens were recorded by synchrotron radiation micro-computed tomography. Then, a "virtual" RUS experiment is proposed to quantify the errors induced by imperfect geometry. Results show that for a bone specimen of ∼1° perpendicularity and parallelism errors, an accuracy of a few percent ( <6.2%) for all the stiffness constants and engineering moduli is achievable.

A critical assessment of the in-vitro measurement of cortical bone stiffness with ultrasound

Elasticity assessment based on bulk wave velocity (BWV) measurements is the most popular technique to characterize the anisotropic stiffness tensor in cortical bone. Typically, a cuboid bone specimen is cut with its sides along the different anatomical directions. Then, the velocity of shear and longitudinal waves propagating along different directions are assessed, from which stiffness coefficients are calculated. Despite the importance of obtaining accurate elasticity values for bone research, there is no generally accepted protocol to measure BWV and the precision of the technique has been seldom investigated. The purpose of this work is to critically assess the method to measure BWV on cuboid specimens in terms of ultrasound frequency, specimen size and signal processing technique. In this study, we measured polycarbonate specimens of different dimensions and 55 human bone specimens with different transducers using frequencies ranging from 2.25 to 10MHz and 1-5MHz for longitudinal and shear waves, respectively. We compared four signal processing methods to detect the wave arrival time. The main results are that, (1) the measurement of shear waves is more complex than that of longitudinal wave, being less precise and more sensitive to sample size; (2) the estimated stiffness depends on the signal processing technique used (up to 10% variation for shear coefficients of bone); and (3) bone stiffness assessed from BWV using the first arrival of the signal to determine the time-of-flight is not different from stiffness assessed using resonant ultrasound spectroscopy (RUS). These results evidence that the measurement method can have an effect on the stiffness values estimates and hence, a well-defined protocol is needed to accurately measure bone stiffness coefficients based on BWV.