Rotational stability in screw-fixed scaphoid fractures compared to plate-fixed scaphoid fractures.

BACKGROUND: The literature describes the treatment of scaphoid fractures comparing the volar and dorsal approaches, the advantages and disadvantages of percutaneous screw fixation, as well as the treatment of scaphoid nonunions using different types of cancellous or corticocancellous bone grafts. Yet, to date no studies are available comparing the outcome of rotational stability in screw-fixed scaphoid fractures to angular stable systems. The purpose of this study is to provide reliable data about rotational stability in stabilised scaphoid fractures and to gain information about the rigidity and the stability of the different types of fixation. METHODS: Three groups of different stabilisation methods on standardised scaphoid B2 fractures were tested for rotational stability. Stabilisation was achieved using one or two cannulated compression screws (CCS) or angular stable plating. We performed ten repetitive cycles up to 10°, 20° and 30° rotation, measuring the maximum torque and the average dissipated work at angle level. RESULTS: Our study showed that rotational stability using a two CCS fixation is significantly (p < 0.05) higher than single CCS fixation. Using the angular stable plate system was also superior to the single CCS (p < 0.05). There was, however, no significant difference between two CCS fixation and angular stable plate fixation. CONCLUSION: Even though indications of using screws or plate systems might be different and plate osteosynthesis may be preferable in treatment of dislocated or comminuted fractures as well as for nonunions, our study showed a better rotational stability by choosing more than just one screw for osteosynthesis. Angular stable plating of scaphoid fractures also provides more rotational stability than single CCS fixation. The authors therefore hypothesise higher union rates in scaphoid fractures using more stable fixation systems.

QCT-based finite element models predict human vertebral strength in vitro significantly better than simulated DEXA.

SUMMARY: While dual energy X-ray absorptiometry (DXA) is considered the gold standard to evaluate fracture risk in vivo, in the present study, the quantitative computed tomography (QCT)-based finite element modeling has been found to provide a quantitative and significantly improved prediction of vertebral strength in vitro. This technique might be used in vivo considering however the much larger doses of radiation needed for QCT. INTRODUCTION: Vertebral fracture is a common medical problem in osteoporotic individuals. Bone mineral density (BMD) is the gold standard measure to evaluate fracture risk in vivo. QCT-based finite element (FE) modeling is an engineering method to predict vertebral strength. The aim of this study was to compare the ability of FE and clinical diagnostic tools to predict vertebral strength in vitro using an improved testing protocol. METHODS: Thirty-seven vertebral sections were scanned with QCT and high resolution peripheral QCT (HR-pQCT). Bone mineral content (BMC), total BMD (tBMD), areal BMD from lateral (aBMD-lat), and anterior-posterior (aBMD-ap) projections were evaluated for both resolutions. Wedge-shaped fractures were then induced in each specimen with a novel testing setup. Nonlinear homogenized FE models (hFE) and linear micro-FE (µFE) were generated from QCT and HR-pQCT images, respectively. For experiments and models, both structural properties (stiffness, ultimate load) and material properties (apparent modulus and strength) were computed and compared. RESULTS: Both hFE and µFE models predicted material properties better than structural ones and predicted strength significantly better than aBMD computed from QCT and HR-pQCT (hFE: R² = 0.79, µFE: R² = 0.88, aBMD-ap: R² = 0.48-0.47, aBMD-lat: R² = 0.41-0.43). Moreover, the hFE provided reasonable quantitative estimations of the experimental mechanical properties without fitting the model parameters. CONCLUSIONS: The QCT-based hFE method provides a quantitative and significantly improved prediction of vertebral strength in vitro when compared to simulated DXA. This superior predictive power needs to be verified for loading conditions that simulate even more the in vivo case for human vertebrae.

Selection of animal bone surrogate samples for orthopaedic screw testing based on human radius CT-derived bone morphology.

Animal bones are commonly used to test the mechanical competence of bone screws since they are easier to obtain compared to human bones. Nevertheless, selecting an appropriate animal sample that correctly represents the human bone architecture where the screw is implanted is frequently overlooked. This study presents a protocol for bone sample selection for screw mechanical testing based on a characterization of the local CT-derived bone morphology. For this, 36 human radii were used to quantify the local peri-implant bone morphology of 360 osteosynthesis screws, 10 per bone, whose implantation site and depth were fully known. A cylindrical volume of interest was created along the screw path and used to measure the local morphology. With this, 10 average peri-implant bone morphologies were defined. Additionally, two animal models, pig, and sheep, were selected and used as potential sample sources. From each model, six bones were selected for analysis. Based on a surface mesh of each bone a computational algorithm was created to automatically extract cylindrical probes in several locations from which the local bone morphometry was calculated. A multi-parametric bone similarity score was developed and used to compare the local morphology of each animal bone to that of the human average peri-implant bone morphology. The score was then mapped to the surface of the bone thus allowing to visually identify regions on the animal bone with human-like bone morphology. By using this methodology, the use of human bones can be avoided since samples with human-like bone morphologies can be found on animal bones. This is not only useful in cases where strict ethical constrains must be fulfilled, but also in studies where the relationship between morphology and screw competence is to be studied, something that is hard to replicate with commercially available synthetic alternatives.

A calibration methodology of QCT BMD for human vertebral body with registered micro-CT images.

PURPOSE: The accuracy of QCT-based homogenized finite element (FE) models is strongly related to the accuracy of the prediction of bone volume fraction (BV/TV) from bone mineral density (BMD). The goal of this study was to establish a calibration methodology to relate the BMD computed with QCT with the BV/TV computed with micro-CT (microCT) over a wide range of bone mineral densities and to investigate the effect of region size in which BMD and BV/TV are computed. METHODS: Six human vertebral bodies were dissected from the spine of six donors and scanned submerged in water with QCT (voxel size: 0.391 x 0.391 x 0.450 mm3) and microCT (isotropic voxel size: 0.018(3) mm3). The microCT images were segmented with a single level threshold. Afterward, QCT-grayscale, microCT-grayscale, and microCT-segmented images were registered. Two isotropic grids of 1.230 mm (small) and 4.920 mm (large) were superimposed on every image, and QCT(BMD) was compared both with microCT(BMD) and microCT(BV/TV) for each grid cell. RESULTS: The ranges of QCT(BMD) for large and small regions were 9-559 mg/cm3 and -90 to 1006 mg/cm3, respectively. QCT(BMD) was found to overestimate microCT(BMD). No significant differences were found between the QCT(BMD)-microCT(BV/TV) regression parameters of the two grid sizes. However, the R2 was higher, and the standard error of the estimate (SEE) was lower for large regions when compared to small regions. For the pooled data, an extrapolated QCTBMD value equal to 1062 mg/ cm3 was found to correspond to 100% microCT(BV/TV). CONCLUSIONS: A calibration method was defined to evaluate BV/TV from QCTBMD values for cortical and trabecular bone in vitro. The QCT(BMD-microCT(BV/TV) calibration was found to be dependent on the scanned vertebral section but not on the size of the regions. However, the higher SEE computed for small regions suggests that the deleterious effect of QCT image noise on FE modelling increases with decreasing voxel size.

Mapping anisotropy improves QCT-based finite element estimation of hip strength in pooled stance and side-fall load configurations.

Hip fractures are one of the most severe consequences of osteoporosis. Compared to the clinical standard of DXA-based aBMD at the femoral neck, QCT-based FEA delivers a better surrogate of femoral strength and gains acceptance for the calculation of hip fracture risk when a CT reconstruction is available. Isotropic, homogenised voxel-based, finite element (hvFE) models are widely used to estimate femoral strength in cross-sectional and longitudinal clinical studies. However, fabric anisotropy is a classical feature of the architecture of the proximal femur and the second determinant of the homogenised mechanical properties of trabecular bone. Due to the limited resolution, fabric anisotropy cannot be derived from clinical CT reconstructions. Alternatively, fabric anisotropy can be extracted from HR-pQCT images of cadaveric femora. In this study, fabric anisotropy from HR-pQCT images was mapped onto QCT-based hvFE models of 71 human proximal femora for which both HR-pQCT and QCT images were available. Stiffness and ultimate load computed from anisotropic hvFE models were compared with previous biomechanical tests in both stance and side-fall configurations. The influence of using the femur-specific versus a mean fabric distribution on the hvFE predictions was assessed. Femur-specific and mean fabric enhance the prediction of experimental ultimate force for the pooled, i.e. stance and side-fall, (isotropic: r2=0.81, femur-specific fabric: r2=0.88, mean fabric: r2=0.86,p<0.001) but not for the individual configurations. Fabric anisotropy significantly improves bone strength prediction for the pooled configurations, and mapped fabric provides a comparable prediction to true fabric. The mapping of fabric anisotropy is therefore expected to help generate more accurate QCT-based hvFE models of the proximal femur for personalised or multiple load configurations.

The effect of the extensor mechanism on maximum isometric fingertip forces: A numerical study on the index finger.

The extensor mechanism is a tendinous network connecting intrinsic and extrinsic muscles of the finger and its function has not yet been fully understood. The goal of this study was to assess the effect of the extensor mechanism on the maximum isometric fingertip forces - a parameter which is essential for grasping. For this purpose, maximum fingertip forces in all directions (i.e. feasible force sets) of two musculoskeletal models of the index finger were compared: the wEM model included a full representation of the extensor mechanism, whereas in the noEM model the extensor mechanism was replaced by a single extensor tendon without connectivity to intrinsic muscles. The feasible force sets were computed in the flexion-extension plane for nine postures. Forces in four predefined directions (palmar, proximal, dorsal, and distal), and the peak resultant forces were evaluated. Averaged forces in all four predefined directions were considerably larger in the wEM model (+187.6%). However, peak resultant forces were slightly lower in the wEM model (-4.3% on average). The general advantage of the wEM model could be explained by co-contraction of intrinsic and extrinsic extensor muscles which allowed reaching larger activation levels of the extrinsic flexors. Only within a narrow range of force directions the co-contraction of intrinsic muscles limited the fingertip forces and lead to lower peak resultant forces in the wEM model. Rather than maximizing peak resultant forces, it appears that the extensor mechanism is a sophisticated tool for increasing maximum fingertip forces over a broad range of postures and force directions - making the finger more versatile during grasping.

Effect of boundary conditions on yield properties of human femoral trabecular bone.

Trabecular bone plays an important mechanical role in bone fractures and implant stability. Homogenized nonlinear finite element (FE) analysis of whole bones can deliver improved fracture risk and implant loosening assessment. Such simulations require the knowledge of mechanical properties such as an appropriate yield behavior and criterion for trabecular bone. Identification of a complete yield surface is extremely difficult experimentally but can be achieved in silico by using micro-FE analysis on cubical trabecular volume elements. Nevertheless, the influence of the boundary conditions (BCs), which are applied to such volume elements, on the obtained yield properties remains unknown. Therefore, this study compared homogenized yield properties along 17 load cases of 126 human femoral trabecular cubic specimens computed with classical kinematic uniform BCs (KUBCs) and a new set of mixed uniform BCs, namely periodicity-compatible mixed uniform BCs (PMUBCs). In stress space, PMUBCs lead to 7-72 % lower yield stresses compared to KUBCs. The yield surfaces obtained with both KUBCs and PMUBCs demonstrate a pressure-sensitive ellipsoidal shape. A volume fraction and fabric-based quadric yield function successfully fitted the yield surfaces of both BCs with a correlation coefficient [Formula: see text]. As expected, yield strains show only a weak dependency on bone volume fraction and fabric. The role of the two BCs in homogenized FE analysis of whole bones will need to be investigated and validated with experimental results at the whole bone level in future studies.

Experimental validation of DXA-based finite element models for prediction of femoral strength.

Osteoporotic fractures are a major clinical problem and current diagnostic tools have an accuracy of only 50%. The aim of this study was to validate dual energy X-rays absorptiometry (DXA)-based finite element (FE) models to predict femoral strength in two loading configurations. Thirty-six pairs of fresh frozen human proximal femora were scanned with DXA and quantitative computed tomography (QCT). For each pair one femur was tested until failure in a one-legged standing configuration (STANCE) and one by replicating the position of the femur in a fall onto the greater trochanter (SIDE). Subject-specific 2D DXA-based linear FE models and 3D QCT-based nonlinear FE models were generated for each specimen and used to predict the measured femoral strength. The outcomes of the models were compared to standard DXA-based areal bone mineral density (aBMD) measurements. For the STANCE configuration the DXA-based FE models (R(2)=0.74, SEE=1473N) outperformed the best densitometric predictor (Neck_aBMD, R(2)=0.66, SEE=1687N) but not the QCT-based FE models (R(2)=0.80, SEE=1314N). For the SIDE configuration both QCT-based FE models (R(2)=0.85, SEE=455N) and DXA neck aBMD (R(2)=0.80, SEE=502N) outperformed DXA-based FE models (R(2)=0.77, SEE=529N). In both configurations the DXA-based FE model provided a good 1:1 agreement with the experimental data (CC=0.87 for SIDE and CC=0.86 for STANCE), with proper optimization of the failure criteria. In conclusion we found that the DXA-based FE models are a good predictor of femoral strength as compared with experimental data ex vivo. However, it remains to be investigated whether this novel approach can provide good predictions of the risk of fracture in vivo.

Experimental validation of a nonlinear µFE model based on cohesive-frictional plasticity for trabecular bone.

Trabecular bone is a porous mineralized tissue playing a major load bearing role in the human body. Prediction of age-related and disease-related fractures and the behavior of bone implant systems needs a thorough understanding of its structure-mechanical property relationships, which can be obtained using microcomputed tomography-based finite element modeling. In this study, a nonlinear model for trabecular bone as a cohesive-frictional material was implemented in a large-scale computational framework and validated by comparison of µFE simulations with experimental tests in uniaxial tension and compression. A good correspondence of stiffness and yield points between simulations and experiments was found for a wide range of bone volume fraction and degree of anisotropy in both tension and compression using a non-calibrated, average set of material parameters. These results demonstrate the ability of the model to capture the effects leading to failure of bone for three anatomical sites and several donors, which may be used to determine the apparent behavior of trabecular bone and its evolution with age, disease, and treatment in the future.

The influence of bone density and anisotropy in finite element models of distal radius fracture osteosynthesis: Evaluations and comparison to experiments.

Continuum-level finite element (FE) models can be used to analyze and improve osteosynthesis procedures for distal radius fractures (DRF) from a biomechanical point of view. However, previous models oversimplified the bone material and lacked thorough experimental validation. The goal of this study was to assess the influence of local bone density and anisotropy in FE models of DRF osteosynthesis for predictions of axial stiffness, implant plate stresses, and screw loads. Experiments and FE analysis were conducted in 25 fresh frozen cadaveric radii with DRFs treated by volar locking plate osteosynthesis. Specimen specific geometries were captured using clinical quantitative CT (QCT) scans of the prepared samples. Local bone material properties were computed based on high resolution CT (HR-pQCT) scans of the intact radii. The axial stiffness and individual screw loads were evaluated in FE models, with (1) orthotropic inhomogeneous (OrthoInhom), (2) isotropic inhomogeneous (IsoInhom), and (3) isotropic homogeneous (IsoHom) bone material and compared to the experimental axial stiffness and screw-plate interface failures. FE simulated and experimental axial stiffness correlated significantly (p<0.0001) for all three model types. The coefficient of determination was similar for OrthoInhom (R(2)=0.807) and IsoInhom (R(2)=0.816) models but considerably lower for IsoHom models (R(2)=0.500). The peak screw loads were in qualitative agreement with experimental screw-plate interface failure. Individual loads and implant plate stresses of IsoHom models differed significantly (p<0.05) from OrthoInhom and IsoInhom models. In conclusion, including local bone density in FE models of DRF osteosynthesis is essential whereas local bone anisotropy hardly effects the models׳ predictive abilities.

Orthotropic HR-pQCT-based FE models improve strength predictions for stance but not for side-way fall loading compared to isotropic QCT-based FE models of human femurs.

Quantitative computed tomography (QCT) based nonlinear homogenized finite element (hFE) models of the human femur do not take bone׳s microstructure into account due to the low resolution of the QCT images. Models based on high-resolution peripheral quantitative computed tomography (HR-pQCT) are able to include trabecular orientation and allow the modeling of a cortical shell. Such a model showed improvements compared to QCT-based models when studying human vertebral bodies. The goal of this study was to compare the femoral strength prediction ability of subject specific nonlinear homogenized FE (hFE) models based on HR-pQCT and QCT images. Thirty-six pairs of femurs were scanned with QCT as well as HR-pQCT, and tested in one-legged stance (STANCE) and side-ways fall (SIDE) configurations up to failure. Non-linear hFE models were generated from HR-pQCT images (smooth meshes) and compared to recently published QCT based models (voxel meshes) as well as experiments with respect to ultimate force. HR-pQCT-based hFE models improved ultimate force (R(2)=0.87 vs 0.80, p=0.02) predictions only in STANCE configuration but not in SIDE (R(2)=0.86 vs 0.84, p=0.6). Damage locations were similar for both types of models. In conclusion, it was shown for the first time on a large femur dataset that a more accurate representation of trabecular orientation and cortex only improve FE predictions in STANCE configuration, where the main trabecular orientation is aligned with the load direction. In the clinically more relevant SIDE configuration, the improvements were not significant.

Identification of a crushable foam material model and application to strength and damage prediction of human femur and vertebral body.

Finite element (FE) models allow quantitative predictions of bone strength and fracture location and, thus, became increasingly popular for assessing fracture risk or effectiveness of osteoporosis therapies. However, predictions crucially depend on the used material models, which are usually complex and rely on a large number of parameters. Therefore, the goal of this study was to propose a simple crushable foam (CF) material model and to perform an extensive comparison with data from the literature. Material parameters of the CF plasticity model were identified based on previously published yield stress data. Voxel-based FE models of thirty-six femora pairs and thirty-eight vertebral bodies were generated from QCT images. The femora models were analyzed with boundary conditions simulating one-legged stance and fall on the greater trochanter. The vertebral body models were subjected to uniaxial compression. Load-displacement curves, ultimate forces and damage distributions computed with the CF material model were compared to a reference material model as well as to in vitro experiments. The result showed that the FE models with CF material provided reasonable quantitative predictions of the ultimate forces measured in the experiments (R(2)>0.80). Comparison of the FE results obtained with CF and reference material model showed very similar outcomes regarding ultimate force, load-displacement behavior and damage patterns for all investigated anatomical sites and loading conditions. In conclusion, the identified CF material model provided good strength and damage predictions, required only few material parameters and is already implemented in many commercial FE solvers. Thus, it can be easily used in other studies.

A nonlinear QCT-based finite element model validation study for the human femur tested in two configurations in vitro.

PURPOSE: Femoral fracture is a common medical problem in osteoporotic individuals. Bone mineral density (BMD) is the gold standard measure to evaluate fracture risk in vivo. Quantitative computed tomography (QCT)-based homogenized voxel finite element (hvFE) models have been proved to be more accurate predictors of femoral strength than BMD by adding geometrical and material properties. The aim of this study was to evaluate the ability of hvFE models in predicting femoral stiffness, strength and failure location for a large number of pairs of human femora tested in two different loading scenarios. METHODS: Thirty-six pairs of femora were scanned with QCT and total proximal BMD and BMC were evaluated. For each pair, one femur was positioned in one-legged stance configuration (STANCE) and the other in a sideways configuration (SIDE). Nonlinear hvFE models were generated from QCT images by reproducing the same loading configurations imposed in the experiments. For experiments and models, the structural properties (stiffness and ultimate load), the failure location and the motion of the femoral head were computed and compared. RESULTS: In both configurations, hvFE models predicted both stiffness (R(2)=0.82 for STANCE and R(2)=0.74 for SIDE) and femoral ultimate load (R(2)=0.80 for STANCE and R(2)=0.85 for SIDE) better than BMD and BMC. Moreover, the models predicted qualitatively well the failure location (66% of cases) and the motion of the femoral head. CONCLUSIONS: The subject specific QCT-based nonlinear hvFE model cannot only predict femoral apparent mechanical properties better than densitometric measures, but can additionally provide useful qualitative information about failure location.

DXA predictions of human femoral mechanical properties depend on the load configuration.

The aim of this study was to evaluate the ability of dual energy X-rays absorptiometry (DXA) areal bone mineral density (aBMD) measured in different regions of the proximal part of the human femur for predicting the mechanical properties of matched proximal femora tested in two different loading configurations. 36 pairs of fresh frozen femora were DXA scanned and tested until failure in two loading configurations: a fall on the side or a one-legged standing. The ability of the DXA output from four different regions of the proximal femur in predicting the femoral mechanical properties was measured and compared for the two loading scenarios. The femoral neck DXA BMD was best correlated to the femoral ultimate force for both configurations and predicted significantly better femoral failure load (R(2)=0.80 vs. R(2)=0.66, P<0.05) when simulating a side than when simulating a standing configuration. Conversely, the work to failure was predicted similarly for both loading configurations (R(2)=0.54 vs. R(2)=0.53, P>0.05). Therefore, neck BMD should be considered as one of the key factors for discriminating femoral fracture risk in vivo. Moreover, the better predictive ability of neck BMD for femoral strength if tested in a fall compared to a one-legged stance configuration suggests that DXA's clinical relevance may not be as high for spontaneous femoral fractures than for fractures associated to a fall.

Tissue properties of the human vertebral body sub-structures evaluated by means of microindentation.

PURPOSE: The better understanding of vertebral mechanical properties can help to improve the diagnosis of vertebral fractures. As the bone mechanical competence depends not only from bone mineral density (BMD) but also from bone quality, the goal of the present study was to investigate the anisotropic indentation moduli of the different sub-structures of the healthy human vertebral body and spondylophytes by means of microindentation. METHODS: Six human vertebral bodies and five osteophytes (spondylophytes) were collected and prepared for microindentation test. In particular, indentations were performed on bone structural units of the cortical shell (along axial, circumferential and radial directions), of the endplates (along the anterio-posterior and lateral directions), of the trabecular bone (along the axial and transverse directions) and of the spondylophytes (along the axial direction). A total of 3164 indentations down to a maximum depth of 2.5 µm were performed and the indentation modulus was computed for each measurement. RESULTS: The cortical shell showed an orthotropic behavior (indentation modulus, Ei, higher if measured along the axial direction, 14.6±2.8 GPa, compared to the circumferential one, 12.3±3.5 GPa, and radial one, 8.3±3.1 GPa). Moreover, the cortical endplates (similar Ei along the antero-posterior, 13.0±2.9 GPa, and along the lateral, 12.0±3.0 GPa, directions) and the trabecular bone (Ei= 13.7±3.4 GPa along the axial direction versus Ei=10.9±3.7 GPa along the transverse one) showed transversal isotropy behavior. Furthermore, the spondylophytes showed the lower mechanical properties measured along the axial direction (Ei=10.5±3.3 GPa). CONCLUSIONS: The original results presented in this study improve our understanding of vertebral biomechanics and can be helpful to define the material properties of the vertebral substructures in computational models such as FE analysis.

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.

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.

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(2) 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.

Validation of an anatomy specific finite element model of Colles' fracture.

Osteoporotic fractures are harmful injuries and their number is on the rise. Distal radius fractures are precursors of other osteoporotic fractures. The wrist's bony geometry and trabecular architecture can be assessed in vivo using the recently introduced HR-pQCT. The goal of this study was the validation of a newly developed HR-pQCT based anatomy specific FE technique including separation of cortical and trabecular bone regions using an experimental model for producing Colles' fractures. Mechanical compression tests of 21 embalmed human radii were conducted. Continuum level FE models were built using HR-pQCT images of the bones and nonlinear analyses were performed using boundary conditions highly similar to the mechanical tests. Density and fabric based material properties were taken from previous tests on biopsies and no adjustments were made. Numerical results provided good prediction of the experimental stiffness (R(2)=0.793) and even better for strength (R(2)=0.874). High damage zones of the FE models coincided with the actual failure patterns of the specimens. These encouraging results allow to conclude that the developed method represents an attractive and efficient tool for simulation of Colles' fracture.

Relationship between ultrasonic parameters and apparent trabecular bone elastic modulus: a numerical approach.

The physical principles underlying quantitative ultrasound (QUS) measurements in trabecular bone are not fully understood. The translation of QUS results into bone strength remains elusive. However, ultrasound being mechanical waves, it is likely to assess apparent bone elasticity. The aim of this study is to derive the sensitivity of QUS parameters to variations of apparent bone elasticity, a surrogate for strength. The geometry of 34 human trabecular bone samples cut in the great trochanter was reconstructed using 3-D synchrotron micro-computed tomography. Finite-difference time-domain simulations coupled to 3-D micro-structural models were performed in the three perpendicular directions for each sample and each direction. A voxel-based micro-finite element linear analysis was employed to compute the apparent Young's modulus (E) of each sample for each direction. For the antero-posterior direction, the predictive power of speed of sound and normalized broadband ultrasonic attenuation to assess E was equal to 0.9 and 0.87, respectively, which is better than what is obtained using bone density alone or coupled with micro-architectural parameters and of the same order of what can be achieved with the fabric tensor approach. When the direction of testing is parallel to the main trabecular orientation, the predictive power of QUS parameters decreases and the fabric tensor approach always gives the best results. This decrease can be explained by the presence of two longitudinal wave modes. Our results, which were obtained using two distinct simulation tools applied on the same set of samples, highlight the potential of QUS techniques to assess bone strength.