Embedding of human vertebral bodies leads to higher ultimate load and altered damage localisation under axial compression

Computer tomography (CT)-based finite element (FE) models of vertebral bodies assess fracture load in vitro better than dual energy X-ray absorptiometry, but boundary conditions affect stress distribution under the endplates that may influence ultimate load and damage localisation under post-yield strains. Therefore, HRpQCT-based homogenised FE models of 12 vertebral bodies were subjected to axial compression with two distinct boundary conditions: embedding in polymethylmethalcrylate (PMMA) and bonding to a healthy intervertebral disc (IVD) with distinct hyperelastic properties for nucleus and annulus. Bone volume fraction and fabric assessed from HRpQCT data were used to determine the elastic, plastic and damage behaviour of bone. Ultimate forces obtained with PMMA were 22% higher than with IVD but correlated highly (R² = 0.99). At ultimate force, distinct fractions of damage were computed in the endplates (PMMA: 6%, IVD: 70%), cortex and trabecular sub-regions, which confirms previous observations that in contrast to PMMA embedding, failure initiated underneath the nuclei in healthy IVDs. In conclusion, axial loading of vertebral bodies via PMMA embedding versus healthy IVD overestimates ultimate load and leads to distinct damage localisation and failure pattern.     

HR-pQCT-based homogenised finite element models provide quantitative predictions of experimental vertebral body stiffness and strength with the same accuracy as μFE models

This study validated two different high-resolution peripheral quantitative computer tomography (HR-pQCT)-based finite element (FE) approaches, enhanced homogenised continuum-level (hFE) and micro-finite element (μFE) models, by comparing them with compression test results of vertebral body sections. Thirty-five vertebral body sections were prepared by removing endplates and posterior elements, scanned with HR-pQCT and tested in compression up to failure. Linear hFE and μFE models were created from segmented and grey-level CT images, and apparent model stiffness values were compared with experimental stiffness as well as strength results. Experimental and numerical apparent elastic properties based on grey-level/segmented CT images (N=35) correlated well for μFE (r2=0.748/0.842) and hFE models (r2=0.741/0.864). Vertebral section stiffness values from the linear μFE/hFE models estimated experimental ultimate apparent strength very well (r2=0.920/0.927). Calibrated hFE models were able to predict quantitatively apparent stiffness with the same accuracy as μFE models. However, hFE models needed no back-calculation of a tissue modulus or any kind of fitting and were computationally much cheaper.     

Removal of the cortical endplates has little effect on ultimate load and damage distribution in QCT-based voxel models of human lumbar vertebrae under axial compression

Every year, 500,000 osteoporotic vertebral compression fractures occur in Europe. Quantitative computed tomography (QCT)-based finite element (FE) voxel models predict ultimate force whether they simulate vertebral bodies embedded in polymethylmethacrylate (PMMA) or vertebral sections with both endplates removed. To assess the effect of endplate removal in those predictions, non-linear FE analyses of QCT-based voxel models of human vertebral bodies were performed. High resolution pQCT images of 11 human lumbar vertebrae without posterior elements were coarsened to clinical resolution and bone volume fraction was used to determine the elastic, plastic and damage behavior of bone tissue. Three model boundary conditions (BCs) were chosen: the endplates were cropped (BC1, BC2) or voxel layers were added on the intact vertebrae to mimic embedding (BC3). For BC1 and BC3, the bottom nodes were fully constrained and the top nodes were constrained transversely while both node sets were freed transversely for BC2. Axial displacement was prescribed to the top nodes. In each model, we compared ultimate force and damage distribution during post-yield loading. The results showed that ultimate forces obtained with BC3 correlated perfectly with those computed with BC1 (R(2)=0.9988) and BC2 (R(2)=0.9987), but were in average 3.4% lower and 6% higher respectively. Moreover, good correlation of damage distribution calculated for BC3 was found with those of BC1 (R(2)=0.92) and BC2 (R(2)=0.73). This study demonstrated that voxel models of vertebral sections provide the same ultimate forces and damage distributions as embedded vertebral bodies, but with less preprocessing and computing time required.     

Noninvasive prediction of vertebral body compressive strength using nonlinear finite element method and an image based technique

Noninvasive prediction of vertebral body strength under compressive loading condition is a valuable tool for the assessment of clinical fractures. This paper presents an effective specimen-specific approach for noninvasive prediction of human vertebral strength using a nonlinear finite element (FE) model and an image based parameter based on the quantitative computed tomography (QCT). Nine thoracolumbar vertebrae excised from three cadavers with an average age of 42 years old were used as the samples. The samples were scanned using the QCT. Then, a segmentation technique was performed on each QCT sectional image. The segmented images were then converted into three-dimensional FE models for linear and nonlinear analyses. A new material model was implemented in our nonlinear model being more compatible with real mechanical behavior of trabecular bone. A new image based MOS (Mechanic of Solids) parameter named minimum sectional strength ((σ_uA)_min) was used for the ultimate compressive strength prediction. Subsequently, the samples were destructively tested under uniaxial compression and their experimental ultimate compressive strengths were obtained. Results indicated that our new implemented FE model can predict ultimate compressive strength of human vertebra with a correlation coefficient (R² = 0.94) better than usual linear and nonlinear FE models (R² = 0.83 and 0.85 respectively). The image based parameter introduced in this study ((σ_uA)_min) was also correlated well with the experimental results (R² = 0.86). Although nonlinear FE method with new implemented material model predicts compressive strength better than the (σ_uA)_min, this parameter is clinically more feasible due to its simplicity and lower computational costs. This can make future applications of the (σ_uA)_min more justified for human vertebral body compressive strength prediction.     

The mechanical behavior of PMMA/bone specimens extracted from augmented vertebrae: a numerical study of interface properties, PMMA shrinkage and trabecular bone damage

Recently published compression tests on PMMA/bone specimens extracted after vertebral bone augmentation indicated that PMMA/bone composites were not reinforced by the trabecular bone at all. In this study, the reasons for this unexpected behavior should be investigated by using non-linear micro-FE models. Six human vertebral bodies were augmented with either standard or low-modulus PMMA cement and scanned with a HR-pQCT system before and after augmentation. Six cylindrical PMMA/bone specimens were extracted from the augmented region, scanned with a micro-CT system and tested in compression. Four different micro-FE models were generated from these images which showed different bone tissue material behavior (with/without damage), interface behavior (perfect bonding, frictionless contact) and PMMA shrinkage due to polymerization. The non-linear stress-strain curves were compared between the different micro-FE models as well as to the compression tests of the PMMA/bone specimens. Micro-FE models with contact between bone and cement were 20% more compliant compared to those with perfect bonding. PMMA shrinkage damaged the trabecular bone already before mechanical loading, which further reduced the initial stiffness by 24%. Progressing bone damage during compression dominated the non-linear part of the stress-strain curves. The micro-FE models including bone damage and PMMA shrinkage were in good agreement with the compression tests. The results were similar with both cements. In conclusion, the PMMA/bone interface properties as well as the initial bone damage due to PMMA polymerization shrinkage clearly affected the stress-strain behavior of the composite and explained why trabecular bone did not contribute to the stiffness and strength of augmented bone.     

A nonlinear finite element model validation study based on a novel experimental technique for inducing anterior wedge-shape fractures in human vertebral bodies in vitro

Vertebral compression fracture is a common medical problem in osteoporotic individuals. The quantitative computed tomography (QCT)-based finite element (FE) method may be used to predict vertebral strength in vivo, but needs to be validated with experimental tests. The aim of this study was to validate a nonlinear anatomy specific QCT-based FE model by using a novel testing setup. Thirty-seven human thoracolumbar vertebral bone slices were prepared by removing cortical endplates and posterior elements. The slices were scanned with QCT and the volumetric bone mineral density (vBMD) was computed with the standard clinical approach. A novel experimental setup was designed to induce a realistic failure in the vertebral slices in vitro. Rotation of the loading plate was allowed by means of a ball joint. To minimize device compliance, the specimen deformation was measured directly on the loading plate with three sensors. A nonlinear FE model was generated from the calibrated QCT images and computed vertebral stiffness and strength were compared to those measured during the experiments. In agreement with clinical observations, most of the vertebrae underwent an anterior wedge-shape fracture. As expected, the FE method predicted both stiffness and strength better than vBMD (R² improved from 0.27 to 0.49 and from 0.34 to 0.79, respectively). Despite the lack of fitting parameters, the linear regression of the FE prediction for strength was close to the 1:1 relation (slope and intercept close to one (0.86 kN) and to zero (0.72 kN), respectively). In conclusion, a nonlinear FE model was successfully validated through a novel experimental technique for generating wedge-shape fractures in human thoracolumbar vertebrae.     

Mapping anisotropy of the proximal femur for enhanced image based finite element analysis

Finite element (FE) models of bone derived from quantitative computed tomography (QCT) rely on realistic material properties to accurately predict bone strength. QCT cannot resolve bone microarchitecture, therefore QCT-based FE models lack the anisotropy apparent within the underlying bone tissue. This study proposes a method for mapping femoral anisotropy using high-resolution peripheral quantitative computed tomography (HR-pQCT) scans of human cadaver specimens. Femur HR-pQCT images were sub-divided into numerous overlapping cubic sub-volumes and the local anisotropy was quantified using a ‘direct-mechanics’ method. The resulting directionality reflected all the major stress lines visible within the trabecular lattice, and provided a realistic estimate of the alignment of Harvesian systems within the cortical compartment. QCT-based FE models of the proximal femur were constructed with isotropic and anisotropic material properties, with directionality interpolated from the map of anisotropy. Models were loaded in a sideways fall configuration and the resulting whole bone stiffness was compared to experimental stiffness and ultimate strength. Anisotropic models were consistently less stiff, but no statistically significant differences in correlation were observed between material models against experimental data. The mean difference in whole bone stiffness between model types was approximately 26%, suggesting that anisotropy can still effect considerable change in the mechanics of proximal femur models. The under prediction of whole bone stiffness in anisotropic models suggests that the orthotropic elastic constants require further investigation. The ability to map mechanical anisotropy from high-resolution images and interpolate information into clinical-resolution models will allow testing of new anisotropic material mapping strategies.     

Preparation of Intact Bovine Tail Intervertebral Discs for Organ Culture

The intervertebral disc (IVD) is the joint of the spine connecting vertebra to vertebra. It functions to transmit loading of the spine and give flexibility to the spine. It composes of three compartments: the innermost nucleus pulposus (NP) encompassing by the annulus fibrosus (AF), and two cartilaginous endplates connecting the NP and AF to the vertebral body on both sides. Discogenic pain possibly caused by degenerative intervertebral disc disease (DDD) and disc herniations has been identified as a major problem in our modern society. To study possible mechanisms of IVD degeneration, in vitro organ culture systems with live disc cells are highly appealing. The in vitro culture of intact bovine coccygeal IVDs has advanced to a relevant model system, which allows the study of mechano-biological aspects in a well-controlled physiological and mechanical environment. Bovine tail IVDs can be obtained relatively easy in higher numbers and are very similar to the human lumbar IVDs with respect to cell density, cell population and dimensions. However, previous bovine caudal IVD harvesting techniques retaining cartilaginous endplates and bony endplates failed after 1-2 days of culture since the nutrition pathways were obviously blocked by clotted blood. IVDs are the biggest avascular organs, thus, the nutrients to the cells in the NP are solely dependent on diffusion via the capillary buds from the adjacent vertebral body. Presence of bone debris and clotted blood on the endplate surfaces can hinder nutrient diffusion into the center of the disc and compromise cell viability. Our group established a relatively quick protocol to "crack"-out the IVDs from the tail with a low risk for contamination. We are able to permeabilize the freshly-cut bony endplate surfaces by using a surgical jet lavage system, which removes the blood clots and cutting debris and very efficiently reopens the nutrition diffusion pathway to the center of the IVD. The presence of growth plates on both sides of the vertebral bone has to be avoided and to be removed prior to culture. In this video, we outline the crucial steps during preparation and demonstrate the key to a successful organ culture maintaining high cell viability for 14 days under free swelling culture. The culture time could be extended when appropriate mechanical environment can be maintained by using mechanical loading bioreactor. The technique demonstrated here can be extended to other animal species such as porcine, ovine and leporine caudal and lumbar IVD isolation.     

Micromechanics of the human vertebral body for forward flexion

To provide mechanistic insight into the etiology of osteoporotic wedge fractures, we investigated the spatial distribution of tissue at the highest risk of initial failure within the human vertebral body for both forward flexion and uniform compression loading conditions. Micro-CT-based linear elastic finite element analysis was used to virtually load 22 human T9 vertebral bodies in either 5° of forward flexion or uniform compression; we also ran analyses replacing the simulated compliant disc (E=8 MPa) with stiff polymethylmethacrylate (PMMA, E=2500 MPa). As expected, we found that, compared to uniform compression, forward flexion increased the overall endplate axial load on the anterior half of the vertebra and shifted the spatial distribution of high-risk tissue within the vertebra towards the anterior aspect of the vertebral body. However, despite that shift, the high-risk tissue remained primarily within the central regions of the trabecular bone and endplates, and forward flexion only slightly altered the ratio of cortical-to-trabecular load sharing at the mid-vertebral level (mean±SD for n=22: 41.3±7.4% compression; 44.1±8.2% forward flexion). When the compliant disc was replaced with PMMA, the anterior shift of high-risk tissue was much more severe. We conclude that, for a compliant disc, a moderate degree of forward flexion does not appreciably alter the spatial distribution of stress within the vertebral body.     

Development and validation of a bioreactor system for dynamic loading and mechanical characterization of whole human intervertebral discs in organ culture

Intervertebral disc (IVD) degeneration is a common cause of back pain, and attempts to develop therapies are frustrated by lack of model systems that mimic the human condition. Human IVD organ culture models can address this gap, yet current models are limited since vertebral endplates are removed to maintain cell viability, physiological loading is not applied, and mechanical behaviors are not measured. This study aimed to (i) establish a method for isolating human IVDs from autopsy with intact vertebral endplates, and (ii) develop and validate an organ culture loading system for human or bovine IVDs. Human IVDs with intact endplates were isolated from cadavers within 48 h of death and cultured for up to 21 days. IVDs remained viable with ~80% cell viability in nucleus and annulus regions. A dynamic loading system was designed and built with the capacity to culture 9 bovine or 6 human IVDs simultaneously while applying simulated physiologic loads (maximum force: 4 kN) and measuring IVD mechanical behaviors. The loading system accurately applied dynamic loading regimes (RMS error <2.5 N and total harmonic distortion <2.45%), and precisely evaluated mechanical behavior of rubber and bovine IVDs. Bovine IVDs maintained their mechanical behavior and retained >85% viable cells throughout the 3 week culture period. This organ culture loading system can closely mimic physiological conditions and be used to investigate response of living human and bovine IVDs to mechanical and chemical challenges and to screen therapeutic repair techniques.     

A patient-specific finite element methodology to predict damage accumulation in vertebral bodies under axial compression, sagittal flexion and combined loads

Due to the inherent limitations of DXA, assessment of the biomechanical properties of vertebral bodies relies increasingly on CT-based finite element (FE) models, but these often use simplistic material behaviour and/or single loading cases. In this study, we applied a novel constitutive law for bone elasticity, plasticity and damage to FE models created from coarsened pQCT images of human vertebrae, and compared vertebral stiffness, strength and damage accumulation for axial compression, anterior flexion and a combination of these two cases. FE axial stiffness and strength correlated with experiments and were linearly related to flexion properties. In all loading modes, damage localised preferentially in the trabecular compartment. Damage for the combined loading was higher than cumulated damage produced by individual compression and flexion. In conclusion, this FE method predicts stiffness and strength of vertebral bodies from CT images with clinical resolution and provides insight into damage accumulation in various loading modes.     

Harmonizing finite element modelling for non-invasive strength estimation by high-resolution peripheral quantitative computed tomography

The finite element (FE) method based on high-resolution peripheral quantitative computed tomography (HR-pQCT) use a variety of tissue constitutive properties and boundary conditions at different laboratories making comparison of mechanical properties difficult. Furthermore, the advent of a second-generation HR-pQCT poses challenges due to improved resolution and a larger region of interest (ROI). This study addresses the need to harmonize results across FE models. The aims are to establish the relationship between FE results as a function of boundary conditions and a range of tissue properties for the first-generation HR-pQCT system, and to determine appropriate model parameters for the second-generation HR-pQCT system. We implemented common boundary conditions and tissue properties on a large cohort (N = 1371), and showed the relationships were highly linear (R² > 0.99) for yield strength and reaction force between FE models. Cadaver radii measured on both generation HR-pQCT with matched ROIs were used to back-calculate a tissue modulus that accounts for the increased resolution (61 µm versus 82 µm), resulting in a modulus of 8748 MPa for second-generation HR-pQCT to produce bone yield strength and reaction force equivalent to using 6829 MPa for first-generation HR-pQCT. Finally, in vivo scans (N = 61) conducted on both generations demonstrated that the larger ROI in the second-generation system results in stronger bone outcome measures, suggesting it is not advisable to convert FE results across HR-pQCT generations without matching ROIs. Together, these findings harmonize FE results by providing a means to compare findings with different boundary conditions and tissue properties, and across scanner generations.     

The effects of bone and pore volume fraction on the mechanical properties of PMMA/bone biopsies extracted from augmented vertebrae

Vertebroplasty forms a porous PMMA/bone composite which was shown to be weaker and less stiff than pure PMMA. It is not known what determines the mechanical properties of such composites in detail. This study investigated the effects of bone volume fraction (BV/TV), cement porosity (PV/(TV-BV), PV…pore volume) and cement stiffness. Nine human vertebral bodies were augmented with either standard or low-modulus PMMA cement and scanned with a HR-pQCT system before and after augmentation. Fourteen cylindrical PMMA/bone biopsies were extracted from the augmented region, scanned with a micro-CT system and tested in compression until failure. Micro-finite element (FE) models of the complete biopsies, of the trabecular bone alone as well as of the porous cement alone were generated from CT images to gain more insight into the role of bone and pores. PV/(TV-BV) and experimental moduli of standard/low-modulus cement (R(2)=0.91/0.98) as well as PV/(TV-BV) and yield stresses (R(2)=0.92/0.83) were highly correlated. No correlation between BV/TV (ranging from 0.057 to 0.138) and elastic moduli was observed (R(2)< 0.05). Interestingly, the micro-FE models of the porous cement alone reproduced the experimental elastic moduli of the standard/low-modulus cement biopsies (R(2)=0.75/0.76) more accurately than the models with bone (R(2)=0.58/0.31). In conclusion, the mechanical properties of the biopsies were mainly determined by the cement porosity and the cement material properties. The study showed that bone tissue inside the biopsies was mechanically "switched off" such that load was carried essentially by the porous PMMA.     

Evaluation of filler materials used for uniform load distribution at boundaries during structural biomechanical testing of whole vertebrae

This study was designed to compare the compressive mechanical properties of filler materials, Wood's metal, dental stone, and polymethylmethacrylate (PMMA), which are widely used for performing structural testing of whole vertebrae. The effect of strain rate and specimen size on the mechanical properties of the filler materials was examined using standardized specimens and mechanical testing. Because Wood's metal can be reused after remelting, the effect of remelting on the mechanical properties was tested by comparing them before and after remelting. Finite element (FE) models were built to simulate the effect of filler material size and properties on the stiffness of vertebral body construct in compression. Modulus, yield strain, and yield strength were not different between batches (melt-remelt) of Wood's metal. Strain rate had no effect on the modulus of Wood's metal, however, Young's modulus decreased with increasing strain rate in dental stone whereas increased in PMMA. Both Wood's metal and dental stone were significantly stiffer than PMMA (12.7 +/- 1.8 GPa, 10.4 +/- 3.4 GPa, and 2.9 +/- 0.4 GPa, respectively). PMMA had greater yield strength than Wood's metal (62.9 +/- 8.7 MPa and 26.2 +/- 2.6 MPa). All materials exhibited size-dependent modulus values. The FE results indicated that filler materials, if not accounted for, could cause more than 9% variation in vertebral body stiffness. We conclude that Wood's metal is a superior moldable bonding material for biomechanical testing of whole bones, especially whole vertebrae, compared to the other candidate materials.     

Confocal microscopy demonstrates association of LTBP-2 in fibrillin-1 microfibrils and colocalisation with perlecan in the disc cell pericellular matrix

Comparative immunolocalisations of latent transforming growth factor-beta-1 binding protein (LTBP)-2, fibrillin-1, versican and perlecan were undertaken in foetal human and wild type C57BL/6 mouse and Hspg2 exon 3 null HS deficient mouse intervertebral discs (IVDs). LTBP-2 was a prominent pericellular component of annular fibrochondrocytes in the posterior annulus fibrosus (AF), interstitial matrix adjacent to nucleus pulposus (NP) cells and to fibrillar and cell associated material in the anterior AF of the human foetal IVD and also displayed a pericellular localisation pattern in murine IVDs. Perlecan and LTBP-2 displayed strong pericellular colocalisation patterns in the posterior AF and to fibrillar material in the outer anterior AF in the foetal human IVD. Versican was a prominent fibril-associated component in the posterior and anterior AF, localised in close proximity to fibrillin-1 in fibrillar arrangements in the cartilaginous vertebral rudiments around paraspinal blood vessels, to major collagen fibre bundles in the anterior and posterior AF and shorter fibres in the NP. Fibrillin-1 was prominent in the outer anterior AF of the human foetal IVD and in fibres extending from the AF into the cartilaginous vertebral rudiments. LTBP-2 was prominently associated with annular fibrils containing fibrillin-1, versican was localised in close proximity to these but not specifically with LTBP-2. The similar deposition levels of LTBP-2 observed in the AF of the Hspg2 exon 3 null and wild type murine IVDs indicated that perlecan HS was not essential for LTBP-2 deposition but colocalisation of LTBP-2 with perlecan in the foetal human IVD was consistent with HS mediated interactions which have already been demonstrated in-vitro.     

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.     

Shh signaling from the nucleus pulposus is required for the postnatal growth and differentiation of the mouse intervertebral disc

Intervertebral discs (IVD) are essential components of the vertebral column. They maintain separation, and provide shock absorbing buffers, between adjacent vertebrae, while also allowing movements between them. Each IVD consists of a central semi-liquid nucleus pulposus (NP) surrounded by a multi-layered fibrocartilagenous annulus fibrosus (AF). Although the IVDs grow and differentiate after birth along with the vertebral column, little is known about the mechanism of this. Understanding the signals that control normal IVD growth and differentiation would also provide potential therapies for degenerative disc disease, which is the major cause of lower back pain and affects a large proportion of the population. In this work, we show that during postnatal growth of the mouse, Sonic hedgehog (Shh) signaling from the NP cells controls many aspects of growth and differentiation of both the NP cells themselves and of the surrounding AF, and that it acts, at least partly, by regulating other signaling pathways in the NP and AF. Recent studies have shown that the NP cells arise from the embryonic notochord, which acts as a major signaling center in the embryo. This work shows that this notochord-derived tissue continues to carry out a major signaling function in the postnatal body and that the IVDs are signaling centers, in addition to their already known functions in the mechanics of vertebral column function.     

Mechanics of bone/PMMA composite structures: an in vitro study of human vertebrae

The goal of this study was to provide material property data for the cement/bone composite resulting from the introduction of PMMA bone cement into human vertebral bodies. A series of quasistatic tensile and compressive mechanical tests were conducted using cement/bone composite structures machined from cement-infiltrated vertebral bodies. Experiments were performed both at room temperature and at body temperature. We found that the modulus of the composite structures was lower than bulk cement (p<0.0001). For compression at 37( composite function)C: composite =2.3+/-0.5GPa, cement =3.1+/-0.2GPa; at 23( composite function)C: composite =3.0+/-0.3GPa, cement =3.4+/-0.2GPa. Specimens tested at room temperature were stiffer than those tested at body temperature (p=0.0004). Yield and ultimate strength factors for the composite were all diminished (55-87%) when compared to cement properties. In general, computational models have assumed that cement/bone composite had the same modulus as cement. The results of this study suggest that computational models of cement infiltrated vertebrae and cemented arthroplasties could be improved by specifying different material properties for cement and cement/bone composite.     

Femoral strength is better predicted by finite element models than QCT and DXA.

Clinicians and patients would benefit if accurate methods of predicting and monitoring bone strength in-vivo were available. A group of 51 human femurs (age range 21-93; 23 females, 28 males) were evaluated for bone density and geometry using quantitative computed tomography (QCT) and dual energy X-ray absorptiometry (DXA). Regional bone density and dimensions obtained from QCT and DXA were used to develop statistical models to predict femoral strength ex vivo. The QCT data also formed the basis of a three-dimensional finite element (FE) models to predict structural stiffness. The femurs were separated into two groups; a model training set (n = 25) was used to develop statistical models to predict ultimate load, and a test set (n = 26) was used to validate these models. The main goal of this study was to test the ability of DXA, QCT and FE techniques to predict fracture load non-invasively, in a simple load configuration which produces predominantly femoral neck fractures. The load configuration simulated the single stance phase portion of normal gait; in 87% of the specimens, clinical appearing sub-capital fractures were produced. The training/test study design provided a tool to validate that the predictive models were reliable when used on specimens with "unknown" strength characteristics. The FE method explained at least 20% more of the variance in strength than the DXA models. Planned refinements of the FE technique are expected to further improve these results. Three-dimensional FE models are a promising method for predicting fracture load, and may be useful in monitoring strength changes in vivo.     

The influence of material property and morphological parameters on specimen-specific finite element models of porcine vertebral bodies

The use of finite element (FE) methods in spinal research is increasing, but there is only limited information available on the influence of different input parameters on the model predictions. The aim of this study was to investigate the role of these parameters in FE models of the vertebral body. Experimental tests were undertaken on porcine lumbar vertebral bodies and scans of the specimens were used to create specimen-specific FE models. Three models were created for each specimen with combinations of generic and specimen-specific parameters. Stiffness and strength predictions were also made directly from the specimen trabecular bone volume fraction (BVF) and cross-sectional area (CSA). The agreement between the experimental results and the FE models with generic morphology was poorer (concordance coefficients = 0.058, 0.125 for stiffness, strength) than those made from the BVF and CSA (concordance coefficients = 0.638, 0.609). The greatest levels of agreement were found with the morphologically specific models including element-specific material properties (concordance coefficients = 0.881, 0.752). This indicates that highly specific models, both in terms of morphology and bone quality, are necessary if the FE tool is to be used effectively for spinal research and clinical practice.