Orthotropic properties of cancellous bone modelled as parameterized cellular material

Analysis of stresses and strains in bone tissues and simulation of their adaptive remodelling require exhaustive information about distribution of constitutive properties of cancellous bone and their relationships to microstructural parameters. Homogenization of "equivalent" trabecular microstructures appears to be an advantageous tool for this task. In this study, parameterized orthotropic constitutive models of cancellous bone are derived from finite element analysis of repeatable microstructure cells. The models, based on a space-filling dodecahedron, are fully three-dimensional and are parameterized with four shape parameters. Variation of the parameters allows to imitate most of typical microstructure patterns observed in real bones, along with a variety of intermediate geometries. Finite element models of cells are generated by a special-purpose structured mesh generator for any arbitrary set of shape parameter values. Static numerical tests are performed for an exhaustive number of parameter value sets (microstructure instances). Coefficients of elastic orthotropic stiffness matrix are determined as tabularized functions of elastic constants versus the shape parameters. Additionally, they are correlated to apparent density and principal fabric tensor values. Comparison of the results with micro-FE data obtained for a large set of cancellous bone specimens proves a good agreement.     

Mimetization of the elastic properties of cancellous bone via a parameterized cellular material

Bone tissue mechanical properties and trabecular microarchitecture are the main factors that determine the biomechanical properties of cancellous bone. Artificial cancellous microstructures, typically described by a reduced number of geometrical parameters, can be designed to obtain a mechanical behavior mimicking that of natural bone. In this work, we assess the ability of the parameterized microstructure introduced by Kowalczyk (Comput Methods Biomech Biomed Eng 9:135–147, 2006. doi: 10.1080/10255840600751473) to mimic the elastic response of cancellous bone. Artificial microstructures are compared with actual bone samples in terms of elasticity matrices and their symmetry classes. The capability of the parameterized microstructure to combine the dominant isotropic, hexagonal, tetragonal and orthorhombic symmetry classes in the proportions present in the cancellous bone is shown. Based on this finding, two optimization approaches are devised to find the geometrical parameters of the artificial microstructure that better mimics the elastic response of a target natural bone specimen: a Sequential Quadratic Programming algorithm that minimizes the norm of the difference between the elasticity matrices, and a Pattern Search algorithm that minimizes the difference between the symmetry class decompositions. The pattern search approach is found to produce the best results. The performance of the method is demonstrated via analyses for 146 bone samples.     

Cyclic cryopreservation affects the nanoscale material properties of trabecular bone

Tissues such as bone are often stored via freezing, or cryopreservation. During an experimental protocol, bone may be frozen and thawed a number of times. For whole bone, the mechanical properties (strength and modulus) do not significantly change throughout five freeze-thaw cycles. Material properties at the trabecular and lamellar scales are distinct from whole bone properties, thus the impact of freeze-thaw cycling at this scale is unknown. To address this, the effect of repeated freezing on viscoelastic material properties of trabecular bone was quantified via dynamic nanoindentation. Vertebrae from five cervine spines (1.5-year-old, male) were semi-randomly assigned, three-to-a-cycle, to 0–10 freeze-thaw cycles. After freeze-thaw cycling, the vertebrae were dissected, prepared and tested. ANOVA (factors cycle, frequency, and donor) on storage modulus, loss modulus, and loss tangent, were conducted. Results revealed significant changes between cycles for all material properties for most cycles, no significant difference across most of the dynamic range, and significant differences between some donors. Regression analysis showed a moderate positive correlation between cycles and material property for loss modulus and loss tangent, and weak negative correlation for storage modulus, all correlations were significant. These results indicate that not only is elasticity unpredictably altered, but also that damping and viscoelasticity tend to increase with additional freeze-thaw cycling.     

The contribution of cancellous bone to long bone strength and rigidity

A recent article (Burr and Piotrowski, 1982) suggested that structural analyses of long bone cross-sectional geometry will be inaccurate and should be considered inappropriate when cancellous bone accounts for 10-15% or more of the cross-sectional area. Consideration of material property differences between compact and cancellous bone, however, indicates that even significant proportions of cancellous bone (10-40% of total cross-sectional area) will very likely have negligible effects on bone strength and rigidity, and can be effectively ignored in geometrical analyses of diaphyseal sections. In metaphyseal and epiphyseal regions, however, geometric analyses of section properties such as area moments of inertia are inappropriate, both because of significant trabecular bone effects, and because of the inherent constraints of mechanical beam models.     

Nanostructure and elastic modulus of single trabecula in bovine cancellous bone

We aimed to investigate the elastic modulus of trabeculae using tensile tests and assess the effects of nanostructure at the hydroxyapatite (HAp) crystal scale on the elastic modulus. In the experiments, 18 trabeculae that were at least 3 mm in length in the proximal epiphysis of three adult bovine femurs were used. Tensile tests were conducted using a small tensile testing device coupled with microscopy under air-dried condition. The c-axis orientation of HAp crystals and the degree of orientation were measured by X-ray diffraction. To observe the deformation behavior of HAp crystals under tensile loading, the same tensile tests were conducted in X-ray diffraction measurements. The mineral content of specimens was evaluated using energy dispersive X-ray spectrometry. The elastic modulus of a single trabecula varied from 4.5 to 23.6 GPa, and the average was 11.5±5.0 GPa. The c-axis of HAp crystals was aligned with the trabecular axis and the crystals were lineally deformed under tensile loading. The ratio of the HAp crystal strain to the tissue strain (strain ratio) had a significant correlation with the elastic modulus (r=0.79; P<0.001). However, the mineral content and the degree of orientation did not vary widely and did not correlate with the elastic modulus in this study. It suggests that the strain ratio may represent the nanostructure of a single trabecula and would determine the elastic modulus as well as mineral content and orientation.     

Influence of boundary conditions on computed apparent elastic properties of cancellous bone

High-resolution finite element models of trabecular bone can be used to study trabecular structure–function relationships, elasticity, multiaxial strength, and tissue remodelling in more detail than experiments. Beside effects of the model size, scan/analysis resolution, segmentation process, etc., the type of the applied boundary conditions (BCs) have a strong influence on the predicted elastic properties. Appropriate BCs have to be applied on hexahedral digital finite element models in order to obtain effective elastic properties. Homogeneous displacement BCs as proposed by Van Rietbergen et al. (J Biomech 29(12):1653–1657, 1996) lead to “apparent” rather than to “effective” elastic properties. This study provides some answers concerning such differences by comparing various BC types (uniform displacement, mixed BCs, periodic BCs), different volume element definitions (original and mirrored models), and several bone volume fractions (BVTV ranging from 6.5 to 37.6%). First, the mixed BCs formulated by Hazanov (Arch Appl Mech 68(6):385–394, 1998) are theoretically extended to shear loading of a porous media. Second, six human bone samples are analyzed, their orthotropic Young’s moduli, shear moduli, and Poisson’s ratios computed and compared. It is found that the proposed mixed BCs give exactly the same effective elastic properties as periodic BCs if a periodic and orthotropic micro-structured material is used and thus denoted as “periodicity compatible” mixed uniform BCs (PMUBCs). As bone samples were shown to be nearly orthotropic for volume element side lengths ≥5 mm the proposed mixed BCs turn out to be the best choice because they give again essentially the same overall elastic properties as periodic BCs. For bone samples of smaller dimensions ( < 5 mm) with a strong anisotropy (beyond orthotropy) uniform displacement BCs remain applicable but they can significantly overestimate the effective stiffness. In Memoriam, Prof. Christian Huet.     

Design and validation of a testing system to assess torsional cancellous bone failure in conjunction with time-lapsed micro-computed tomographic imaging

When compressed axially, cancellous bone often fails at an oblique angle along well-defined bands, highlighting the importance of cancellous bone shear properties. Torsion testing to determine shear properties of cancellous bone has often been conducted under conditions appropriate only for axis-symmetric specimens comprised of homogeneous and isotropic materials. However, most cancellous bone specimens do not meet these stringent test conditions. Therefore, the aim of this study was to design and validate a uniaxial, incremental torsional testing system for non-homogeneous orthotropic or non-axis-symmetric specimens.Precision and accuracy of the newly designed torsion system was validated by using Plexiglas rods and beams, where obtained material properties were compared to those supplied by the manufacturer. Additionally, the incremental step-wise application of angular displacement and simultaneous time-lapsed μCT imaging capability of the system was validated using whale cancellous bone specimens, with step-wise application of angular displacement yielding similar torsional mechanical properties to continuous application of angular displacement in a conventional torsion study.In conclusion, a novel torsion testing system for non-homogeneous, orthotropic materials using the incremental step-wise application of torsion and simultaneous time-lapsed μCT imaging was designed and validated.     

The effect of charge density on the velocity and attenuation of ultrasound waves in human cancellous bone

Cancellous bone is a highly porous material, and two types of waves, fast and slow, are observed when ultrasound is used for detecting bone diseases. There are several possible stimuli for bone remodelling processes, including bone fluid flow, streaming potential, and piezoelectricity. Poroelasticity has been widely used for elucidating the bone fluid flow phenomenon, but the combination of poroelasticity with charge density has not been introduced. Theoretically, general poroelasticity with a varying charge density is employed for determining the relationship between wave velocity and attenuation with charge density. Fast wave velocity and attenuation are affected by porosity as well as charge density; however, for a slow wave, both slow wave velocity and attenuation are not as sensitive to the effect of charge density as they are for a fast wave. Thus, employing human femoral data, we conclude that charged ions gather on trabecular struts, and the fast wave, which moves along the trabecular struts, is significantly affected by charge density.     

A method for the determination of the residual chloroform in defatted cancellous bone transplants

The removal of fat from cancellous bone tissue promotes the clinical healing of the transplant and improves the penetration of chemical sterilisation media into the tissue. Treatment using chloroform/methanol (2:1, 2 h) is frequently used as a defatting procedure. Eight rinses with methanol followed by two rinses with aqua ad iniectabilia (20 min each, with ultrasonic effect) ensure depletion in the level of chloroform from defatted cancellous bone to a concentration below 25 ppm (limit value). For the necessary routine quality checks on the production process, a gas chromatography method has been developed that determines the level of chloroform in cancellous bone, for which the detection limit is 0.003 ppm.     

Mechanical properties of compacted morselized cancellous bone graft using one-dimensional consolidation testing

Failures of orthopaedic procedures that use morselized cancellous bone (MCB) graft for load bearing are often due to gross displacement within the graft material. For this reason the mechanical behavior of MCB must be better understood. Our purpose is to present a detailed testing methodology for the mechanical characterization of MCB, and to illustrate how this methodology can be used to study the influence of water and fat content. Complete one-dimensional consolidation testing was performed on bovine cancellous bone processed to represent MCB typically used in surgery (52% water, 31% fat). The one-dimensional consolidation strain under a stress of 1.09 MPa was 30.9% and the confined modulus was 8.0 MPa. The coefficient of consolidation (rate of consolidation) was 2.2×10⁻⁵ cm²/s and the coefficient of secondary strain (steady-state creep rate) was 1.9%. While reducing the water content alone had some influence on properties, reducing the fat content improved both the static and dynamic behavior. A sample of MCB which had fat intentionally minimized and a lower overall moisture content (56% water, 5% fat) demonstrated 23.1% strain, a confined modulus of 9.6 MPa, a coefficient of consolidation of 3.4×10⁻³ cm²/s, and a coefficient of secondary strain of 0.9%. The test methods described in this technical note can be used to evaluate the influence of fluid content on the mechanical behavior of MCB.     

The effects of tensile-compressive loading mode and microarchitecture on microdamage in human vertebral cancellous bone

The amount of microdamage in bone tissue impairs mechanical performance and may act as a stimulus for bone remodeling. Here we determine how loading mode (tension vs. compression) and microstructure (trabecular microarchitecture, local trabecular thickness, and presence of resorption cavities) influence the number and volume of microdamage sites generated in cancellous bone following a single overload. Twenty paired cylindrical specimens of human vertebral cancellous bone from 10 donors (47–78 years) were mechanically loaded to apparent yield in either compression or tension, and imaged in three dimensions for microarchitecture and microdamage (voxel size 0.7×0.7×5.0 μm³). We found that the overall proportion of damaged tissue was greater (p=0.01) for apparent tension loading (3.9±2.4%, mean±SD) than for apparent compression loading (1.9±1.3%). Individual microdamage sites generated in tension were larger in volume (p<0.001) but not more numerous (p=0.64) than sites in compression. For both loading modes, the proportion of damaged tissue varied more across donors than with bone volume fraction, traditional measures of microarchitecture (trabecular thickness, trabecular separation, etc.), apparent Young?s modulus, or strength. Microdamage tended to occur in regions of greater trabecular thickness but not near observable resorption cavities. Taken together, these findings indicate that, regardless of loading mode, accumulation of microdamage in cancellous bone after monotonic loading to yield is influenced by donor characteristics other than traditional measures of microarchitecture, suggesting a possible role for tissue material properties.     

Micromechanical modeling of elastic properties of cortical bone accounting for anisotropy of dense tissue

The paper analyzes the connection between microstructure of the osteonal cortical bone and its overall elastic properties. The existing models either neglect anisotropy of the dense tissue or simplify cortical bone microstructure (accounting for Haversian canals only). These simplifications (related mostly to insufficient mathematical apparatus) complicate quantitative analysis of the effect of microstructural changes – produced by age, microgravity, or some diseases – on the overall mechanical performance of cortical bone. The present analysis fills this gap; it accounts for anisotropy of the dense tissue and uses realistic model of the porous microstructure. The approach is based on recent results of Sevostianov et al. (2005) and Saadat et al. (2012) on inhomogeneities in a transversely-isotropic material. Bone?s microstructure is modeled according to books of Martin and Burr (1989), Currey (2002), and Fung (1993) and includes four main families of pores. The calculated elastic constants for porous cortical bone are in agreement with available experimental data. The influence of each of the pore types on the overall moduli is examined.     

Functional adaptation of cancellous bone in human proximal femur predicted by trabecular surface remodeling simulation toward uniform stress state

Two-dimensional simulation of trabecular surface remodeling was conducted for a human proximal femur to investigate the structural change of cancellous bone toward a uniform stress state. Considering that a local mechanical stimulus plays an important role in cellular activities in bone remodeling, local stress nonuniformity was assumed to drive trabecular structural change to seek a uniform stress state. A large-scale pixel-based finite element model was used to simulate structural changes of individual trabeculae over the entire bone. As a result, the initial structure of trabeculae changed from isotropic to anisotropic due to trabecular microstructural changes caused by surface remodeling according to the mechanical environment in the proximal femur. Under a single-loading condition, it was shown that the apparent structural property evaluated by fabric ellipses corresponded to the apparent stress state in cancellous bone. As is observed in the actual bone, a distributed trabecular structure was obtained under a multiple-loading condition. Through these studies, it was concluded that trabecular surface remodeling toward a local uniform stress state at the trabecular level could naturally bring about functional adaptation phenomenon at the apparent tissue level. The proposed simulation model would be capable of providing insight into the hierarchical mechanism of trabecular surface remodeling at the microstructural level up to the apparent tissue level.     

Assessment of vertebral wedge strength using cancellous textural properties derived from digital tomosynthesis and density properties from dual energy X-ray absorptiometry and high resolution computed tomography

The purpose of this study was to examine the potential of digital tomosynthesis (DTS) derived cancellous bone textural measures to predict vertebral strength under conditions simulating a wedge fracture. 40 vertebral bodies (T6, T8, T11, and L3 levels) from 5 male and 5 female cadaveric donors were utilized. The specimens were scanned using dual energy X-ray absorptiometry (DXA) and high resolution computed tomography (HRCT) to obtain measures of bone mineral density (BMD) and content (BMC), and DTS to obtain measures of bone texture. Using a custom loading apparatus designed to deliver a nonuniform displacement resulting in a wedge deformity similar to those observed clinically, the specimens were loaded to fracture and their fracture strength was recorded. Mixed model regressions were used to determine the associations between wedge strength and DTS derived textural variables, alone and in the presence of BMD or BMC information. DTS derived fractal, lacunarity and mean intercept length variables correlated with wedge strength, and individually explained up to 53% variability. DTS derived textural variables, notably fractal dimension and lacunarity, contributed to multiple regression models of wedge strength independently from BMC and BMD. The model from a scan orientation transverse to the spine axis and in the anterior-posterior view resulted in highest explanatory capability (R²_adj = 0.91), with a scan orientation parallel to the spine axis and in the lateral view offering an alternative (R²_adj = 0.88). In conclusion, DTS can be used to examine cancellous texture relevant to vertebral wedge strength, and potentially complement BMD in assessment of vertebral fracture risk.     

Finite element calculated uniaxial apparent stiffness is a consistent predictor of uniaxial apparent strength in human vertebral cancellous bone tested with different boundary conditions

Strong correspondence between the uniaxial apparent strength and stiffness of cancellous bone allows the use of stiffness as a predictor of bone strength. Measured values of mechanical properties in cancellous bone can be different between experiments due to different experimental conditions. In the current study, bone volume fraction, experimentally determined and finite element (FE) predicted stiffness were examined as predictors of cancellous bone ultimate strength in two different groups each of which was tested using a different end constraint. It is demonstrated that, although always significant, the relationships of strength with bone volume fraction and experimentally determined stiffness are different between test groups. Apparent stiffness, estimated by FE modeling, predicts the ultimate strength of human cancellous bone consistently for all examined experimental protocols.     

Effects of damage on the orthotropic material symmetry of bovine tibial trabecular bone

The macroscopic mechanical properties of trabecular bone can be predicted by its architecture using theoretical relationships between the elastic and architectural properties. Microdamage caused by overloading or fatigue decreases the apparent elastic moduli of trabecular bone requiring these relationships to be modified to predict the damaged elastic properties. In the case of isotropic damage, the apparent level elastic properties could be determined by multiplying all of the elastic constants by a single scalar factor. If the damage is anisotropic, the elastic constants may change by differing factors and the material coordinate system could become misaligned with the fabric coordinate system. High-resolution finite element models were used to simulate damage overloading on seven trabecular bone specimens subjected to pure shear strain in two planes. Comparison of the apparent elastic moduli of the specimens before and after damage showed that the reduction of the elastic moduli was anisotropic. This suggests that the microdamage within the specimens was inhomogeneous. However, after damage the specimens exhibited nearly orthotropic material symmetry as they did before damage. Changes in the orientation of the orthotropic material coordinate system were also small and occurred primarily in the transverse plane. Thus, while damage in trabecular bone is anisotropic, the material coordinate system remains aligned with the fabric tensor.     

Macroscopic anisotropic bone material properties in children with severe osteogenesis imperfecta

Children with severe osteogenesis imperfecta (OI) typically experience numerous fractures and progressive skeletal deformities over their lifetime. Recent studies proposed finite element models to assess fracture risk and guide clinicians in determining appropriate intervention in children with OI, but lack of appropriate material property inputs remains a challenge. This study aimed to characterize macroscopic anisotropic cortical bone material properties and investigate relationships with bone density measures in children with severe OI. Specimens were obtained from tibial or femoral shafts of nine children with severe OI and five controls. The specimens were cut into beams, characterized in bending, and imaged by synchrotron radiation X-ray micro-computed tomography. Longitudinal modulus of elasticity, yield strength, and bending strength were 32–65% lower in the OI group (p < 0.001). Yield strain did not differ between groups (p ≥ 0.197). In both groups, modulus and strength were lower in the transverse direction (p ≤ 0.009), but anisotropy was less pronounced in the OI group. Intracortical vascular porosity was almost six times higher in the OI group (p < 0.001), but no differences were observed in osteocyte lacunar porosity between the groups (p = 0.086). Volumetric bone mineral density was lower in the OI group (p < 0.001), but volumetric tissue mineral density was not (p = 0.770). Longitudinal OI bone modulus and strength were correlated with volumetric bone mineral density (p ≤ 0.024) but not volumetric tissue mineral density (p ≥ 0.099). Results indicate that cortical bone in children with severe OI yields at the same strain as normal bone, and that their decreased bone material strength is associated with reduced volumetric bone mineral density. These results will enable the advancement of fracture risk assessment capability in children with severe OI.     

Multi-axial mechanical properties of human trabecular bone

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

Age-dependent material properties of the porcine cerebrum: effect on pediatric inertial head injury criteria

During growth and development, the immature central nervous system undergoes rapid alterations in constituents and structure. We hypothesize that these alterations are accompanied by changes in the mechanical properties of brain tissue which, in turn, influence the response of the brain to traumatic inertial loads. Samples of frontal cerebrum from neonatal (2–3 days) and adult pigs were harvested and tested within 3 h post-mortem. The complex shear modulus of the samples was measured in a custom-designed oscillatory shear testing device at engineering shear strain amplitudes of 2.5% or 5% from 20–200 Hz, at 25°C and 100% humidity. In this range, the elastic and viscous components of the complex shear modulus increased significantly with the development of the cerebral region of the brain. Using an idealized model of the developing head, the age-dependent material properties of brain tissue were shown to affect the mechanical response of the brain to inertial loading. This study is a first step toward developing head injury tolerance criteria specifically for the pediatric population.