Mechanical properties of trabecular bone. Dependency on strain rate

3D printed trabeculae conditionally reproduce the mechanical properties of the actual trabeculae - A preliminary study

Three-dimensional (3D) printing has been used to fabricate synthetic trabeculae models and to test mechanical behavior that cannot be recognized in the actual sample, but the extent to which 3D printed trabeculae replicate the mechanical behavior of the actual trabeculae remains to be quantified. The aim of this study was to evaluate the accuracy of 3D printed trabeculae in reproducing the mechanical properties of the corresponding actual trabeculae. Twelve human trabecular cubes (5 × 5 × 5 mm) were scanned by micro-CT to form the trabecular 3D model. Each trabecular 3D model was scaled ×2-, ×3-, ×4- and ×5-fold and then printed twice at a layer thickness of 60 μm using poly (lactic acid) (PLA). The actual trabecular cubes and the 3D-printed trabecular cubes were first compressed under a loading rate of 1 mm/min; another replicated stack of 3D-printed trabecular cubes was compressed under a strain rate of 0.2/min. The results showed that the stiffness of the printed cubes tended to increase, while the strength tended to converge when the magnification increased under the two loading conditions. The strain rate effect was found in the printed cubes. The correlation coefficient (R²) of the mechanical properties between the printed and actual trabeculae can reach up to 0.94, especially under ×3-, ×4- and ×5-fold magnification. In conclusion, 3D printing could be a potential tool to evaluate the mechanical behavior of actual trabecular tissue in vitro and may help in the future to predict the risk of fracture and even personalize the treatment evaluation for osteoporosis and other trabecular bone pathologies.

Biomechanical and Microstructural Properties of Subchondral Bone From Three Metacarpophalangeal Joint Sites in Thoroughbred Racehorses

Fatigue-induced subchondral bone (SCB) injury is common in racehorses. Understanding how subchondral microstructure and microdamage influence mechanical properties is important for developing injury prevention strategies. Mechanical properties of the disto-palmar third metacarpal condyle (MCIII) correlate poorly with microstructure, and it is unknown whether the properties of other sites within the metacarpophalangeal (fetlock) joint are similarly complex. We aimed to investigate the mechanical and structural properties of equine SCB from specimens with minimal evidence of macroscopic disease. Three sites within the metacarpophalangeal joint were examined: the disto-palmar MCIII, disto-dorsal MCIII, and proximal sesamoid bone. Two regions of interest within the SCB were compared, a 2 mm superficial and an underlying 2 mm deep layer. Cartilage-bone specimens underwent micro-computed tomography, then cyclic compression for 100 cycles at 2 Hz. Disto-dorsal MCIII specimens were loaded to 30 MPa (n = 10), while disto-palmar MCIII (n = 10) and proximal sesamoid (n = 10) specimens were loaded to 40 MPa. Digital image correlation determined local strains. Specimens were stained with lead-uranyl acetate for volumetric microdamage quantification. The dorsal MCIII SCB had lower bone volume fraction (BVTV), bone mineral density (BMD), and stiffness compared to the palmar MCIII and sesamoid bone (p < 0.05). Superficial SCB had higher BVTV and lower BMD than deeper SCB (p < 0.05), except at the palmar MCIII site where there was no difference in BVTV between depths (p = 0.419). At all sites, the deep bone was stiffer (p < 0.001), although the superficial to deep gradient was smaller in the dorsal MCIII. Hysteresis (energy loss) was greater superficially in palmar MCIII and sesamoid (p < 0.001), but not dorsal MCIII specimens (p = 0.118). The stiffness increased with cyclic loading in total cartilage-bone specimens (p < 0.001), but not in superficial and deep layers of the bone, whereas hysteresis decreased with the cycle for all sites and layers (p < 0.001). Superficial equine SCB is uniformly less stiff than deeper bone despite non-uniform differences in bone density and damage levels. The more compliant superficial layer has an important role in energy dissipation, but whether this is a specific adaptation or a result of microdamage accumulation is not clear.

Three-dimensional finite element analysis of the effect of alveolar cleft bone graft on the maxillofacial biomechanical stabilities of unilateral complete cleft lip and palate

BACKGROUND

The objective is to clarify the effect of alveolar cleft bone graft on maxillofacial biomechanical stabilities, the key areas when bone grafting and in which should be supplemented with bone graft once bone resorption occurred in UCCLP (unilateral complete cleft lip and palate).

METHODS

Maxillofacial CAD (computer aided design) models of non-bone graft and full maxilla cleft, full alveolar cleft bone graft, bone graft in other sites of the alveolar cleft were acquired by processing the UCCLP maxillofacial CT data in three-dimensional modeling software. The maxillofacial bone EQV (equivalent) stresses and bone suture EQV strains under occlusal states were obtained in the finite element analysis software.

RESULTS

Under corresponding occlusal states, the EQV stresses of maxilla, pterygoid process of sphenoid bone on the corresponding side and anterior alveolar arch on the non-cleft side were higher than other maxillofacial bones, the EQV strains of nasomaxillary, zygomaticomaxillary and pterygomaxillary suture on the corresponding side were higher than other maxillofacial bone sutures. The mean EQV strains of nasal raphe, the maximum EQV stresses of posterior alveolar arch on the non-cleft side, the mean and maximum EQV strains of nasomaxillary suture on the non-cleft side in full alveolar cleft bone graft model were all significantly lower than those in non-bone graft model. The mean EQV stresses of bilateral anterior alveolar arches, the maximum EQV stresses of maxilla and its alveolar arch on the cleft side in the model with bone graft in lower 1/3 of the alveolar cleft were significantly higher than those in full alveolar cleft bone graft model.

CONCLUSIONS

For UCCLP, bilateral maxillae, pterygoid processes of sphenoid bones and bilateral nasomaxillary, zygomaticomaxillary, pterygomaxillary sutures, anterior alveolar arch on the non-cleft side are the main occlusal load-bearing structures before and after alveolar cleft bone graft. Alveolar cleft bone graft mainly affects biomechanical stabilities of nasal raphe and posterior alveolar arch, nasomaxillary suture on the non-cleft side. The areas near nasal floor and in the middle of the alveolar cleft are the key sites when bone grafting, and should be supplemented with bone graft when the bone resorbed in these areas.

Investigating Torque-Speed Relationship of Self-Tapping Screws

Correctly torquing bone screws is important to prevent fixation failures and ensure positive patient outcomes. It has been proposed that an automatic model-based method may be able to determine the patient-specific material properties of bone, and provide objective and quantitative torquing recommendations. One major part of developing this system is the modelling of the bone-screwing process, and the self-tapping screwing process in general. In this paper, we investigate the relationship between screw insertion torque (Nm) and speed of insertion (RPM). A weak positive correlation was found below approximately 30 RPM. Further research should focus on increasing the precision of the methodology, and this testing must be extended to ex-vivo animal bone testing in addition to the polyurethane foam substitute used here.Clinical relevance: To maximise the accuracy of torque recommendations, the model should account for all important factors. This study investigates and attempts to quantify the relationship between screw insertion speed and torque for later inclusion in modelling if significant.

Development of a crushable foam model for human trabecular bone

Finite element (FE) simulations can be used to evaluate the mechanical behavior of human bone and allow for quantitative prediction of press-fit implant fixation. An adequate material model that captures post-yield behavior is essential for a realistic simulation. The crushable foam (CF) model is a constitutive model that has recently been proposed in this regard. Compression tests under uniaxial and confined loading conditions were performed on 59 human trabecular bone specimens. Three essential material parameters were obtained as a function of bone mineral density (BMD) to develop the isotropic CF model. The related constitutive rule was implemented in FE models and the results were compared to the experimental data. The CF model provided an accurate simulation of uniaxial compression tests and the post-yield behavior of the stress-strain was well-matched with the experimental results. The model was able to reproduce the confined response of the bone up to 15% of strain. This model allows for simulation of the mechanical behavior of the cellular structure of human bone and adequately predicts the post-yield response of trabecular bone, particularly under uniaxial loading conditions. The model can be further improved to simulate bone collapse due to local overload around orthopaedic implants.

Mechanical Characterization of Human Trabecular and Formed Granulate Bone Cylinders Processed by High Hydrostatic Pressure

One main disadvantage of commercially available allogenic bone substitute materials is the altered mechanical behavior due to applied material processing, including sterilization methods like thermal processing or gamma irradiation. The use of high hydrostatic pressure (HHP) might be a gentle alternative to avoid mechanical alteration. Therefore, we compressed ground trabecular human bone to granules and, afterwards, treated them with 250 and 300 MPa for 20 and 30 min respectively. We characterized the formed bone granule cylinders (BGC) with respect to their biomechanical properties by evaluating stiffness and stress at 15% strain. Furthermore, the stiffness and yield strength of HHP-treated and native human trabecular bone cylinders (TBC) as control were evaluated. The mechanical properties of native vs. HHP-treated TBCs as well as HHP-treated vs. untreated BGCs did not differ, independent of the applied HHP magnitude and duration. Our study suggests HHP treatment as a suitable alternative to current processing techniques for allogenic bone substitutes since no negative effects on mechanical properties occurred.

Estimating tissue-level properties of porcine talar subchondral bone

Tissue-level properties of bone play an important role when characterising apparent-level bone biomechanical behaviour and yet little is known about its effect at this hierarchical level. In combination with trabecular morphological data these properties can be used to predict bone strength, which becomes an invaluable tool for clinicians in patient treatment planning. This study developed specimen-specific micro-finite element (μFE) models using validated continuum-level models, containing grayscale-derived material properties, to indirectly establish tissue-level properties of porcine talar subchondral bone. Specimen-specific continuum finite element (hFE) models of subchondral trabecular bone were setup using μCT data of ten cylindrical specimens extracted from juvenile porcine tali. The models were validated using quasi-static uniaxial compression testing. Validated hFE models were used to calibrate the tissue modulus of corresponding μFE models by minimising the difference between the μFE and hFE stiffness values. Key trabecular morphological indices (BV/TV, DA, Conn.D, Tb.Th, EF) were evaluated. Good agreement was observed between hFE models and experiment (CCC = 0.66). Calibrated E_tiss was 504 ± 37.65 MPa. Average BV/TV and DA for μFE specimens were 0.37 ± 0.05 and 0.68 ± 0.11, respectively. BV/TV (r² = 0.667) correlated highly with μFE stiffness. The small intra-specimen variation to tissue-level properties suggests that variations to apparent-level stiffness originate from variations to microarchitecture rather than tissue mechanical properties.

Structure and properties of baleen in the Southern right (Eubalaena australis) and Pygmy right whales (Caperea marginata)

Baleen is a resilient and keratinised filter-feeding structure attached to the maxilla of mysticete whales. It is strong and tough, yet a pliant and resilient material, that withstands extreme pressures in the oral cavity during feeding. We investigated the structure, water content, wettability and mechanical properties of baleen of the Southern right (SRW) and Pygmy right whales (PRW), to understand the effects of hydration on the physical and mechanical properties of baleen. Sixty 25 × 15mm baleen subsamples were prepared from one individual of SRW and PRW. Half were hydrated in circulated natural seawater for 21 days and half were dry. Water content analysis showed that SRW baleen was 21.2% water weight and PRW was 26.1%. Wettability testing indicated that surfaces of both hydrated and dried SRW and PRW baleen were hydrophilic, with hydrated samples of both species having lower contact angle values. For the SRW, the average contact angle of hydrated baleen was 40° ± 13.2 and 73° ± 6 for dried samples. Hydrated PRW baleen had an average contact angle of 44° ± 15.3, which was lower than in dried samples (74° ± 2.9). Three-point bending mechanical tests showed that the average maximum flexural stress of dried SRW (134.1 ± 34.3 MPa) and PRW samples (117.8 ± 22.3 MPa) were significantly higher than those of hydrated SRW (25.7 ± 6.3 MPa) and PRW (19.7 ± 4.8 MPa) baleen. Scanning electron microscope images showed the stratification of the outer cortical layer, with cross-linked keratin fibres observed within and between baleen keratin sheets. Hydrated baleen, as in its natural and functional behaviour, has greater flexibility and strength, attributes necessary for the complex filter feeding mechanism characteristic of whales. Hydration must be considered when addressing the physical and mechanical properties of baleen, especially when using dried museum specimens.

Empirical relationships between bone density and ultimate strength: A literature review

INTRODUCTION

Ultimate strength-density relationships for bone have been reported with widely varying results. Reliable bone strength predictions are crucial for many applications that aim to assess bone failure. Bone density and bone morphology have been proposed to explain most of the variance in measured bone strength. If this holds true, it could lead to the derivation of a single ultimate strength-density-morphology relationship for all anatomical sites.

METHODS

All relevant literature was reviewed. Ultimate strength-density relationships derived from mechanical testing of human bone tissue were included. The reported relationships were translated to ultimate strength-apparent density relationships and normalized with respect to strain rate. Results were grouped based on bone tissue type (cancellous or cortical), anatomical site, and loading mode (tension vs. compression). When possible, the relationships were compared to existing ultimate strength-density-morphology relationships.

RESULTS

Relationships that considered bone density and morphology covered the full spectrum of eight-fold inter-study difference in reported compressive ultimate strength-density relationships for trabecular bone. This was true for studies that tested specimens in different loading direction and tissue from different anatomical sites. Sparse data was found for ultimate strength-density relationships in tension and for cortical bone properties transverse to the main loading axis of the bone.

CONCLUSIONS

Ultimate strength-density-morphology relationships could explain measured strength across anatomical sites and loading directions. We recommend testing of bone specimens in other directions than along the main trabecular alignment and to include bone morphology in studies that investigate bone material properties. The lack of tensile strength data did not allow for drawing conclusions on ultimate strength-density-morphology relationships. Further studies are needed. Ideally, these studies would investigate both tensile and compressive strength-density relationships, including morphology, to close this gap and lead to more accurate evaluation of bone failure.

Effect of loading frequency on deformations at the bone-implant interface

This study considers the time-dependent behaviour of bone in the context of loosening of metal implants, which is one of the typical complications of joint replacement and fracture-fixation surgeries. We employed viscoelastic properties developed from our previous experimental studies for trabecular bone in a representative bone-implant construct, which was subjected to cyclic loading at varying loading frequencies. We found that the separation between the bone and the implant is a function of loading frequency and increases with number of loading cycles applied. Our analysis shows that at the start of cyclic loading, a higher frequency results in a lower displacement response of bone at the bone-implant interface; however, after the bone-implant system has been subjected to a large number of cycles (>500 cycles in this study), higher interfacial displacements are observed at higher loading frequencies. In other words, higher loading frequencies will not result in bone-implant separation if limited number of cycles are applied. In all cases, interfacial displacements increase as bone volume ratio decreases. This simple approach can be used to evaluate the mechanical environment in bone-implant systems due to cyclic loading which commonly used time-independent models that are unable to simulate. The approach can also be used to evaluate implant loosening due to cyclic loading.

High-speed X-ray visualization of dynamic crack initiation and propagation in bone

The initiation and propagation of physiological cracks in porcine cortical and cancellous bone under high rate loading were visualized using high-speed synchrotron X-ray phase-contrast imaging (PCI) to characterize their fracture behaviors under dynamic loading conditions. A modified Kolsky compression bar was used to apply dynamic three-point flexural loadings on notched specimens and images of the fracture processes were recorded using a synchronized high-speed synchrotron X-ray imaging set-up. Three-dimensional synchrotron X-ray tomography was conducted to examine the initial microstructure of the bone before high-rate experiments. The experimental results showed that the locations of fracture initiations were not significantly different between the two types of bone. However, the crack velocities in cortical bone were higher than in cancellous bone. Crack deflections at osteonal cement lines, a prime toughening mechanism in bone at low rates, were observed in the cortical bone under dynamic loading in this study. Fracture toughening mechanisms, such as uncracked ligament bridging and bridging in crack wake were also observed for the two types of bone. The results also revealed that the fracture toughness of cortical bone was higher than cancellous bone. The crack was deflected to some extent at osteon cement line in cortical bone instead of comparatively penetrating straight through the microstructures in cancellous bone. STATEMENT OF SIGNIFICANCE: Fracture toughness is with great importance to study for crack risk prediction in bone. For those cracks in bone, most of them are associated with impact events, such as sport accidents. Consequently, we visualized, in real-time, the entire processes of dynamic fractures in notched cortical bone and cancellous bone specimens using synchrotron X-ray phase contrast imaging. The onset location of crack initiation was found independent on the bone type. We also found that, although the extent was diminished, crack deflections at osteon cement lines, a major toughening mechanism in transversely orientated cortical bone at quasi-static rate, were still played a role in resisting cracking in dynamically loaded specimen. These finding help researchers to understand the dynamic fracture behaviors in bone.

Effects of loading rate on the of mechanical behavior of the femur in falling condition

The Surgeon General estimates that by 2020, half of all Americans could have weak bones due to bone loss. Osteoporosis causes more than 1.5 million fractures every year. Identifying effective interventions based on individual patient characteristics remains a major challenge. Proximal femur fractures are common and devastating events for individuals with osteoporosis. Since fracture is primarily a mechanical event, noninvasive predictions of fracture strength and location would be useful both for identifying at-risk individuals and evaluating treatment effects. However, bone fracture prediction is complicated due to the complex microstructure and nanostructure of bone. Bone is a highly heterogeneous material with rate-dependent mechanical behavior and large inter-individual variation. In this study, we designed two mechanical test procedures to understand the mechanical response of bone under impact and quasi-static load tests. The boundary conditions of the tests were designed in a way to simulate a fall to the side. The present study consists of three main parts: cadaver testing, quantitative image analysis, and finite element (FE) modeling. We obtained ten human femur bones and used high-resolution CT to quantify the microstructure and density of each sample. Specimen-specific FE models were created to evaluate the ability of various failure criteria to predict experimental fracture. Afterward, the samples were tested and their failure patterns were recorded. The fractured samples were rescanned to analyze the fractured surfaces. Our experimental results show that the loading necessary to fracture the femur samples is much higher in the impact tests. However, the toughening mechanisms are more pronounced in quasi-static tests. We found that FE model formulations were able to accurately predict femur stiffness and strength for quasi-static and impact conditions separately, but that no single formulation could account for the rate-dependent outcomes.

Standardizing compression testing for measuring the stiffness of human bone

Objectives

The ability to determine human bone stiffness is of clinical relevance in many fields, including bone quality assessment and orthopaedic prosthesis design. Stiffness can be measured using compression testing, an experimental technique commonly used to test bone specimens in vitro. This systematic review aims to determine how best to perform compression testing of human bone.

Methods

A keyword search of all English language articles up until December 2017 of compression testing of bone was undertaken in Medline, Embase, PubMed, and Scopus databases. Studies using bulk tissue, animal tissue, whole bone, or testing techniques other than compression testing were excluded.

Results

A total of 4712 abstracts were retrieved, with 177 papers included in the analysis; 20 studies directly analyzed the compression testing technique to improve the accuracy of testing. Several influencing factors should be considered when testing bone samples in compression. These include the method of data analysis, specimen storage, specimen preparation, testing configuration, and loading protocol.

Conclusion

Compression testing is a widely used technique for measuring the stiffness of bone but there is a great deal of inter-study variation in experimental techniques across the literature. Based on best evidence from the literature, suggestions for bone compression testing are made in this review, although further studies are needed to establish standardized bone testing techniques in order to increase the comparability and reliability of bone stiffness studies.

Cite this article: S. Zhao, M. Arnold, S. Ma, R. L. Abel, J. P. Cobb, U. Hansen, O. Boughton. Standardizing compression testing for measuring the stiffness of human bone. Bone Joint Res 2018;7:524–538. DOI: 10.1302/2046-3758.78.BJR-2018-0025.R1.

Bone Mechanical Properties in Healthy and Diseased States

The mechanical properties of bone are fundamental to the ability of our skeletons to support movement and to provide protection to our vital organs. As such, deterioration in mechanical behavior with aging and/or diseases such as osteoporosis and diabetes can have profound consequences for individuals' quality of life. This article reviews current knowledge of the basic mechanical behavior of bone at length scales ranging from hundreds of nanometers to tens of centimeters. We present the basic tenets of bone mechanics and connect them to some of the arcs of research that have brought the field to recent advances. We also discuss cortical bone, trabecular bone, and whole bones, as well as multiple aspects of material behavior, including elasticity, yield, fracture, fatigue, and damage. We describe the roles of bone quantity (e.g., density, porosity) and bone quality (e.g., cross-linking, protein composition), along with several avenues of future research.

Finite-Element Analysis of Bone Stresses on Primary Impact in a Large-Animal Model: The Distal End of the Equine Third Metacarpal

OBJECTIVE

To assess whether the transient stresses of foot impact with the ground are similar to those found during midstance loading and if the location of high stress correlate with the sites most commonly associated with mechanically induced osteoarthritis (OA). We compared impact stresses in subchondral bone between two subject-specific, three-dimensional, finite-element models of the equine metacarpophalangeal (MCP) joint-one with advanced OA and one healthy, and with similar published data on the stresses that occur at midstance.

METHODS

Two right MCP joints (third metacarpal and proximal phalanx) were scanned using micro-computed tomography (μCT). Images were segmented, and meshed using modified 10-node quadratic tetrahedral elements. Bone material properties were assigned based on the bone density. An impact velocity of 3.55 m/s was applied to each model and contact pressures and stress distribution were calculated for each. In a separate iteration, the third metacarpal was loaded statically. A sampling grid of 160 equidistant points was superimposed over selected slices, and average peak stresses were calculated for 6 anatomical regions. Within-region maximal peak and average von Mises stresses were compared between healthy and OA bones in both midstance and impact loading.

RESULTS

Average impact stresses across all regions, in both locations (palmar and dorsal) were greater in the OA model. Highest impact stresses were located in the dorsal medial condyle in the healthy (12.8 MPa) and OA (14.1MPa) models, and were lowest in the palmar medial and lateral parasagittal grooves in the healthy (5.94 MPa) and OA (7.07 MPa) models. The healthy static model had higher peak (up to 49.7% greater) and average (up to 38.6% greater) stresses in both locations and across all regions compared to the OA static model.

CONCLUSIONS

Under simulated footfall a trot, loading on the dorsal aspect of the third metacarpal at impact created stresses similar to those found during midstance. The high accelerations that occur under impact loading are likely responsible for creating the high stresses, as opposed to midstance loading where the high stresses are the result of high mass loading. Although the stress magnitudes were found to be similar among the two loading conditions, the location of the high stress loading occurred in sites that are not typically associated with osteoarthritic changes.

Individual and combined effects of OA-related subchondral bone alterations on proximal tibial surface stiffness: a parametric finite element modeling study

The role of subchondral bone in OA pathogenesis is unclear. While some OA-related changes to morphology and material properties in different bone regions have been described, the effect of these alterations on subchondral bone surface stiffness has not been investigated. The objectives of this study were to characterize the individual (Objective 1) and combined (Objective 2) effects of OA-related morphological and mechanical alterations to subchondral and epiphyseal bone on surface stiffness of the proximal tibia. We developed and validated a parametric FE model of the proximal tibia using quantitative CT images of 10 fresh-frozen cadaveric specimens and in situ macro-indentation testing. Using this validated FE model, we estimated the individual and combined roles of OA-related alterations in subchondral cortical thickness and elastic modulus, and subchondral trabecular and epiphyseal trabecular elastic moduli on local surface stiffness. A 20% increase in subchondral cortical or subchondral trabecular elastic moduli resulted in little change in stiffness (1% increase). A 20% reduction in epiphyseal trabecular elastic modulus, however, resulted in an 11% reduction in stiffness. Our parametric analysis suggests that subchondral bone stiffness is affected primarily by epiphyseal trabecular bone elastic modulus rather than subchondral cortical and trabecular morphology or mechanical properties. Our results suggest that observed OA-related alterations to epiphyseal trabecular bone (e.g., lower mineralization, bone volume fraction, density and elastic modulus) may contribute to OA proximal tibiae being less stiff than normal.

Comparison of proximal femur and vertebral body strength improvements in the FREEDOM trial using an alternative finite element methodology

Denosumab reduced the incidence of new fractures in postmenopausal women with osteoporosis by 68% at the spine and 40% at the hip over 36 months compared with placebo in the FREEDOM study. This efficacy was supported by improvements from baseline in vertebral (18.2%) strength in axial compression and femoral (8.6%) strength in sideways fall configuration at 36 months, estimated in Newtons by an established voxel-based finite element (FE) methodology. Since FE analyses rely on the choice of meshes, material properties, and boundary conditions, the aim of this study was to independently confirm and compare the effects of denosumab on vertebral and femoral strength during the FREEDOM trial using an alternative smooth FE methodology. Unlike the previous FE study, effects on femoral strength in physiological stance configuration were also examined. QCT data for the proximal femur and two lumbar vertebrae were analyzed by smooth FE methodology at baseline, 12, 24, and 36 months for 51 treated (denosumab) and 47 control (placebo) subjects. QCT images were segmented and converted into smooth FE models to compute bone strength. L1 and L2 vertebral bodies were virtually loaded in axial compression and the proximal femora in both fall and stance configurations. Denosumab increased vertebral body strength by 10.8%, 14.0%, and 17.4% from baseline at 12, 24, and 36 months, respectively (p<0.0001). Denosumab also increased femoral strength in the fall configuration by 4.3%, 5.1%, and 7.2% from baseline at 12, 24, and 36 months, respectively (p<0.0001). Similar improvements were observed in the stance configuration with increases of 4.2%, 5.2%, and 5.2% from baseline (p≤0.0007). Differences between the increasing strengths with denosumab and the decreasing strengths with placebo were significant starting at 12 months (vertebral and femoral fall) or 24 months (femoral stance). Using an alternative smooth FE methodology, we confirmed the significant improvements in vertebral body and proximal femur strength previously observed with denosumab. Estimated increases in strength with denosumab and decreases with placebo were highly consistent between both FE techniques.

Prediction of local proximal tibial subchondral bone structural stiffness using subject-specific finite element modeling: Effect of selected density-modulus relationship

BACKGROUND

Quantitative computed tomography based subject-specific finite element modeling has potential to clarify the role of subchondral bone alterations in knee osteoarthritis initiation, progression, and pain initiation. Calculation of bone elastic moduli from image data is a basic step when constructing finite element models. However, different relationships between elastic moduli and imaged density (known as density-modulus relationships) have been reported in the literature. The objective of this study was to apply seven different trabecular-specific and two cortical-specific density-modulus relationships from the literature to finite element models of proximal tibia subchondral bone, and identify the relationship(s) that best predicted experimentally measured local subchondral structural stiffness with highest explained variance and least error.

METHODS

Thirteen proximal tibial compartments were imaged via quantitative computed tomography. Imaged bone mineral density was converted to elastic moduli using published density-modulus relationships and mapped to corresponding finite element models. Proximal tibial structural stiffness values were compared to experimentally measured stiffness values from in-situ macro-indentation testing directly on the subchondral bone surface (47 indentation points).

FINDINGS

Regression lines between experimentally measured and finite element calculated stiffness had R(2) values ranging from 0.56 to 0.77. Normalized root mean squared error varied from 16.6% to 337.6%.

INTERPRETATION

Of the 21 evaluated density-modulus relationships in this study, Goulet combined with Snyder and Schneider or Rho appeared most appropriate for finite element modeling of local subchondral bone structural stiffness. Though, further studies are needed to optimize density-modulus relationships and improve finite element estimates of local subchondral bone structural stiffness.

Bone properties of the humeral head and resistance to screw cutout

Surgical treatment of fractures involving the proximal humeral head is hampered by complications. Screw cutout is the major pitfall seen in connection with rigid plating. We have exploited a bony explanation for this phenomenon.

Development of an Inertia-Driven Model of Sideways Fall for Detailed Study of Femur Fracture Mechanics

A new method for laboratory testing of human proximal femora in conditions simulating a sideways fall was developed. Additionally, in order to analyze the strain state in future cadaveric tests, digital image correlation (DIC) was validated as a tool for strain field measurement on the bone of the femoral neck. A fall simulator which included models for the body mass, combined lateral femur and pelvis mass, pelvis stiffness, and trochanteric soft tissue was designed. The characteristics of each element were derived and developed based on human data from the literature. The simulator was verified by loading a state-of-the-art surrogate femur and comparing the resulting force-time trace to published, human volunteer experiments. To validate the DIC, 20 human proximal femora were prepared with a strain rosette and speckle paint pattern, and loaded to 50% of their predicted failure load at a low compression rate. Strain rosettes were taken as the gold standard, and minimum principal strains from the DIC and the rosettes were compared using descriptive statistics. The initial slope of the force-time curve obtained in the fall simulator matched published human volunteer data, with local peaks superimposed in the model due to internal vibrations of the spring used to model the pelvis stiffness. Global force magnitude and temporal characteristics were within 2% of published volunteer experiments. The DIC minimum principal strains were found to be accurate to 127±239μɛ. These tools will allow more biofidelic laboratory simulation of falls to the side, and more detailed analysis of proximal femur failure mechanisms using human cadaver specimens.

Effect of the strain rate on the twisting of trabecular bone from women with hip fracture

As one of the major functions of bone is to provide structural support for the musculoskeletal system, it is important to evaluate its mechanical strength. Bones may be subjected to multiaxial stresses due to bone pathologies, accidental loads which may lead to hip, wrist fracture, or to a prosthetic joint replacement. Twist loading may lead to fractures, especially involving long bones from lower limbs. The aim of this work was to study the effect of the strain rate on the shear properties of trabecular bone samples from women with hip fracture (from 65 to 100 years). Cylindrical samples were core drilled from human femoral heads along the primary trabecular direction. The cylinder's ends were polished and embedded in blocks of polymeric material which fit the grips of the testing device. Deformation rates of 0.005, 0.01, 0.015, and 0.05 s⁻¹ were applied. Twisting tests were conducted with or without an applied axial load of 500 N. From the torque-angular displacement curves, the shear stress-strain curves were obtained. The maximum shear strength and the shear modulus (i.e. the slope of the linear region) were determined. A large scatter of the results of the shear strength and the shear modulus was found, which is probably related to the heterogeneity of nonhealthy human bone samples. There is no significant effect of the strain rate on the maximum shear stress and the shear modulus, either in tests undertaken with or without the application of an axial load. The effect of strain rate on nonhealthy bone trabecular twisting properties did not follow the trend observed on the effect of strain rate in healthy bone, where an increase is detected.

Low elastic modulus titanium-nickel scaffolds for bone implants

The superelastic nature of repeating the human bones is crucial to the ideal artificial biomedical implants to ensure smooth load transfer and foster the ingrowth of new bone tissues. Three dimensional interconnected porous TiNi scaffolds, which have the tailorable porous structures with micro-hole, were fabricated by slurry immersing with polymer sponge and sintering method. The crystallinity and phase composition of scaffolds were studied by X-ray diffraction. The pore morphology, size and distribution in the scaffolds were characterized by scanning electron microscopy. The porosity ranged from 65 to 72%, pore size was 250-500μm. Compressive strength and elastic modulus of the scaffolds were ~73MPa and ~3GPa respectively. The above pore structural and mechanical properties are similar to those of cancellous bone. In the initial cell culture test, osteoblasts adhered well to the scaffold surface during a short time, and then grew smoothly into the interconnected pore channels. These results indicate that the porous TiNi scaffolds fabricated by this method could be bone substitute materials.

Nanomechanical Characterization of Canine Femur Bone for Strain Rate Sensitivity in the Quasistatic Range under Dry versus Wet Conditions

As a strain rate-dependent material, bone has a different mechanical response to various loads. Our aim was to evaluate the effect of water and different loading/unloading rates on the nanomechanical properties of canine femur cortical bone. Six cross-sections were cut from the diaphysis of six dog femurs and were nanoindented in their cortical area. Both dry and wet conditions were taken into account for three quasistatic trapezoid profiles with a maximum force of 2000 μN (holding time = 30 s) at loading/unloading rates of 10, 100, and 1000 μN/s, respectively. For each specimen, 254 ± 9 (mean ± SD) indentations were performed under different loading conditions. Significant differences were found for the elastic modulus and hardness between wet and dry conditions (P < 0.001). No influence of the loading/unloading rates was observed between groups except for the elastic modulus measured at 1000 μN/s rate under dry conditions (P < 0.001) and for the hardness measured at a rate of 10 μN/s under wet conditions (P < 0.001). Therefore, for a quasistatic test with peak load of 2000 μN held for 30 s, it is recommended to nanoindent under wet conditions at a loading/unloading rate of 100–1000 μN/s, so the reduced creep effect allows for a more accurate computation of mechanical properties.

High speed insertion of bone fracture fixation pins: a finite element penetration model with experimental comparisons

A new method of bone fracture fixation is considered in which small pins/darts are dynamically inserted into bone to prevent translation and rotation at the fracture site. An ABAQUS model was developed to analyze dart penetration in cortical and cancellous bone for varying dart diameter, material, and velocity, and cortical thickness. The method is advocated for bioresorbable darts, so polylactide (PLA) and magnesium are the materials examined in this study. Numerical results showed that magnesium darts can achieve full penetration in bone while suffering little damage. The PLA darts penetrated thin bone well, but substantial deformation was seen as the cortical thickness increased, especially for small diameter darts. As partial validation, prototype PLA fixation darts were fired into cadaveric bone with a custom nailer. As in the model, the PLA darts could penetrate thin cortices but saw gross deformation when impacted against thicker bone.

Sustained-release rhBMP-2 increased bone mass and bone strength in an ovine model of postmenopausal osteoporosis

OBJECTIVE

The purpose of this study was to analyze the local treatment effects of rhBMP-2 combined with fibrin sealant (FS) on bone mineral density, microarchitectural and mechanical properties in osteoporotic ovine spine.

MATERIALS AND METHODS

Postmenopausal osteoporosis was induced in eight sheep through ovariectomy (OVX) and a low-calcium diet for a period of 12 months. According to the Latin square design, L3-L6 vertebrae were randomly assigned to four treatment groups: A (rhBMP-2/FS), B (rhBMP-2), C (FS) and D (blank control). All materials were injected into the assigned vertebra transpedicularly. All animals were euthanized 3 months after treatment. Bone mineral density (BMD), microarchitectural and mechanical properties were assessed. ANOVA analysis of variance was used to determine effects of rhBMP-2/FS (α = 0.05).

RESULTS

The BMD in group A (rhBMP-2/FS) was 18.8, 30.4 and 27.9% higher than that in group B, C and D, respectively. Analysis of bone structure by micro-CT revealed higher trabecular bone volume (BV/TV), trabecular thickness (Tb.Th) and trabecular number (Tb.N) in the rhBMP-2/FS group (P < 0.01). In addition, vertebrae treated with rhBMP-2/FS exhibited higher yield stress, ultimate stress, energy absorption and bone modulus compared to the control groups.

CONCLUSIONS

Local administration of rhBMP-2/FS showed a positive trend in improving BMD, microarchitectural parameters and mechanical strength of osteoporotic vertebra. Slow release of rhBMP-2 using FS appeared to be an effective method of protein delivery. The local treatment of osteoporosis in the spine can increase bone strength and reduce fracture risk quickly.

A new method to evaluate the elastic modulus of cortical bone by using a combined computed tomography and finite element approach

This study attempted to estimate the elastic modulus of cortical specimens directly from the computed tomography (CT) number (CT). Cubic specimens of fresh bovine femora were tested under compressive loading. The corresponding three-dimensional mesh models of these specimens were established with the aid of CT images. The elastic modulus of each element was determined using the following transfer formula: E=Ax(CT)(B). The A and B were determined by matching the simulation results with the experimental force-displacement curves. An optimization approach was used to match the entire specimen sets between simulations and experiments. The results indicated that the elastic modulus of cortical bone can be effectively estimated using a power relationship (E=1.26 x 10(-3) x (CT)(1.93)) directly from the CT number. This procedure eliminates the need to determine the bone density, and therefore reduces the computational time required to define mechanical properties in finite element methods.

Microarchitectural adaptations in aging and osteoarthrotic subchondral bone issues

The human skeleton optimizes its microarchitecture by elaborate adaptations to mechanical loading during development and growth. The mechanisms for adaptation involve a multistep process of cellular mechanotransduction stimulating bone modelling, and remodeling resulting in either bone formation or resorption. This process causes appropriate microarchitectural changes tending to adjust and improve the bone structure to its prevailing mechanical environment. Normal individual reaches peak bone mass at age between 25 and 30 years, and thereafter bone mass declines with age in both genders. The bone loss is accompanied by microarchitectural deterioration resulting in reduced mechanical strength likely leading to fragility fractures. With aging, inevitable bone loss occurs, which is frequently the cause of osteoporosis; and inevitable bone and joint degeneration happens, which often results in osteoarthrosis. These diseases are among the major health care problems in terms of socio-economic costs. The overall goals of the current series of studies were to investigate the age-related and osteoarthrosis (OA) related changes in the 3-D microarchitectural properties, mechanical properties, collagen and mineral quality of subchondral cancellous and cortical bone tissues. The studies included mainly two parts. For human subjects: aging- (I–IV) and early OArelated (V–VI) changes in cancellous bone properties were assessed. For OA guinea pig models (VII–IX), three topics were studied: firstly, the spontaneous, age-related development of guinea pig OA; secondly, the potential effects of hyaluronan on OA subchondral bone tissues; and thirdly, the effects on OA progression of an increase in subchondral bone density by inhibition of bone remodeling with a bisphosphonate. These investigations aimed to obtain more insight into the age-related and OA-related subchondral bone adaptations.

Crack revision improves fixation of uncemented HA-coated implants compared with reaming: an experiment in dogs

The crack procedure is a surgical technique for preparing the implant cavity at revision of loose joint replacement components. It disrupts the neocortical bone shell that typically forms around the cavity. Using an animal model, we compared the crack technique with reaming. Twenty micromotion implants were inserted bilaterally into the knees of 10 dogs according to our revision protocol, allowing formation of a standardized revision cavity (loose implant, fibrous tissue, and sclerotic bone rim). Eight weeks later we performed revision surgery. On the control side, in which the neocortex was removed, the cavity was reamed. On the intervention side, in which the neocortex was perforated but left in situ, the cavity was cracked. For revision we used non-motioning hydroxyapatite (HA)-coated, plasma-sprayed titanium implants. Observation after revision was 4 weeks. The implants revised by the crack technique had better mechanical fixation in all mechanical parameters by the push-out test. The crack revisions also provided more new bone formation around the implants compared with the reamed revisions but had no effect on new bone ongrowth. The data suggest using this bone-sparing technique may be superior to reaming in terms of achieving improved early implant fixation of uncemented HA-coated revision implants.

Role of trabecular microarchitecture in whole-vertebral body biomechanical behavior

The role of trabecular microarchitecture in whole-vertebral biomechanical behavior remains unclear, and its influence may be obscured by such factors as overall bone mass, bone geometry, and the presence of the cortical shell. To address this issue, 22 human T(9) vertebral bodies (11 female; 11 male; age range: 53-97 yr, 81.5 +/- 9.6 yr) were scanned with microCT and analyzed for measures of trabecular microarchitecture, BMC, cross-sectional area, and cortical thickness. Sixteen of the vertebrae were biomechanically tested to measure compressive strength. To estimate vertebral compressive stiffness with and without the cortical shell for all 22 vertebrae, two high-resolution finite element models per specimen-one intact model and one with the shell removed-were created from the microCT scans and virtually compressed. Results indicated that BMC and the structural model index (SMI) were the individual parameters most highly associated with strength (R(2) = 0.57 each). Adding microarchitecture variables to BMC in a stepwise multiple regression model improved this association (R(2) = 0.85). However, the microarchitecture variables in that regression model (degree of anisotropy, bone volume fraction) differed from those when BMC was not included in the model (SMI, mean trabecular thickness), and the association was slightly weaker for the latter (R(2) = 0.76). The finite element results indicated that the physical presence of the cortical shell did not alter the relationships between microarchitecture and vertebral stiffness. We conclude that trabecular microarchitecture is associated with whole-vertebral biomechanical behavior and that the role of microarchitecture is mediated by BMC but not by the cortical shell.

Development of the fetal ilium--challenging concepts of bipedality

Macroradiographs of 30 human fetal and neonatal ilia were analysed to investigate the early pattern of trabecular bone organization prior to the influences of direct weight-bearing locomotion. Consistent and well-defined patterns of internal organization were identified within the fetal and neonatal ilium, which correspond with previously recognized regions that have been attributed directly to forces associated with bipedal locomotion. This study proposes that patterns previously attributed to weight-bearing locomotive responses are present in the earliest stages of the development of this bone. It is suggested that the rudimentary scaffold seen in the fetal and neonatal ilium could indicate a predetermined template upon which locomotive influences may be superimposed and perhaps reinforced at a later age. Alternatively, this early pattern may mimic the adult form due to the effects of in-utero limb movement activity even though it is not weight bearing. This is a preliminary study that will be supported in a further communication with three-dimensional micro-computed trabecular analysis.

Freeze casting of porous hydroxyapatite scaffolds. II. Sintering, microstructure, and mechanical behavior

In Part I, the influence of processing parameters on the general microstructure of freeze-cast hydroxyapatite (HA) constructs was explored. This work is an extension of Part I to investigate the effect of sintering conditions on the microstructure and mechanical behavior of freeze-cast HA. For constructs prepared from aqueous suspensions (5-20 vol % HA), sintering for 3 h at temperatures from 1250 degrees C to 1375 degrees C produced a decrease in porosity of <5% but an increase in strength of nearly 50%. Constructs with a porosity of 52% had compressive strengths of 12 +/- 1 MPa and 5 +/- 1 MPa in the directions parallel and perpendicular to the freezing direction, respectively. The mechanical response showed high strain tolerance (5-10% at the maximum stress), high strain to failure (>20%), and high strain rate sensitivity. Manipulation of the freeze-cast microstructure, achieved by additions of glycerol and 1,4-dioxane to the aqueous suspensions, produced changes in the magnitude of the mechanical response, but little change in the general nature of the response. The favorable mechanical behavior of the porous constructs, coupled with the ability to modify their microstructure, indicates the potential of the present freeze-casting route for the production of porous scaffolds for bone tissue engineering.

Staged-injection procedure to prevent cement leakage during vertebroplasty: an in vitro study

STUDY DESIGN

Fibrin sealant (FS) combined with or without growth factor was used to improve the micro-architectural and biomechanical properties of vertebral body in osteoporotic ovine spine.

OBJECTIVE

To analyze the treatment effects of bovine bone morphogenetic protein (bBMP) combined with FS on osteopenic ovine vertebral architecture, bone mineral density, and biomechanics in vivo.

SUMMARY OF BACKGROUND DATA

Vertebroplasty and kyphoplasty were used to treat spinal osteoporosis. They can increase strength of vertebrae physically. However, each has specific disadvantages. bBMP could rapidly increasing bone formation and suppressing bone resorption, but little is known about its effect on ovariectomized-induced osteoporosis.

METHODS

Six sheep underwent ovariectomy and were placed on a low-calcium diet. Twelve months later, according to Ladin square design, L4-L6 vertebrae in all sheep were randomly assigned to 3 treatment groups: A (30 mg bBMP/1.5 mL FS), B (30 mg bBMP) and C (1.5 mL FS). All materials were injected into the assigned vertebra transpedicularly. Animals were killed 3 months after injection, and bone mineral density (BMD), biomechanics, and histomorphometry were assessed. Analysis of variance was used to determine effects of bBMP/FS (alpha = 0.05).

RESULTS

The BMD in Group 1 was 17.1% and 14.7% higher than that in Group 2 and Group 3, respectively. The micro-CT reconstruction analysis showed that the density and connectivity of trabecular bone in bBMP/FS treated vertebrae were higher than those in control groups. The mechanical properties (yield stress, ultimate stress, energy absorption, bone modulus) of the vertebrae were also significantly higher. In this study, local bBMP/FS treatment showed a positive trend in improving BMD, histomorphometric parameters, and mechanical strength of osteoporotic vertebra. Slow release of bBMP using FS appeared to be an effective method of protein delivery.

CONCLUSION

The use of bBMP/FS in the treatment of vertebral osteoporosis in an attempt to enhance bone strength merits further study.

Computational Determination of the Critical Microcrack Size That Causes a Remodeling Response in a Trabecula: A Feasibility Study

Bone is a living tissue, which undergoes continuous renewal to repair local defects. Two separate processes, adaptation and remodeling, are involved when a defect appears. The defect produces stress concentrations that provoke regional adaptation, and is gradually repaired, first by resorption and then by deposition of new bone. Using a mathematical formulation of the adaptation mechanism in trabeculae of cancellous bone, we hypothesize that in some cases, where a microcrack is small enough relative to the dimensions of the trabecula, the adaptation response of the whole trabecula may be sufficient to regain homeostatic mechanical conditions (with no need for a remodeling process). The simulation results showed that for trabeculae with nominal length of 900 µm and nominal thickness of 80–800 µm, a microcrack with minimal length of 48 µm and minimal depth of 13% of the trabecula’s thickness was required to initiate a remodeling process. A longer (100 µm) but shallower (depth of 7% of the trabecula’s thickness) crack also triggered remodeling. These computational results support our hypothesis that when a microcrack small enough relative to the dimensions of the trabecula occurs, adaptation of the whole trabecula may be sufficient to regain homeostatic mechanical conditions with no need for a local remodeling process.