Benedikt Helgason

Development of a balanced experimental–computational approach to understanding the mechanics of proximal femur fractures

The majority of people who sustain hip fractures after a fall to the side would not have been identified using current screening techniques such as areal bone mineral density. Identifying them, however, is essential so that appropriate pharmacological or lifestyle interventions can be implemented. A protocol, demonstrated on a single specimen, is introduced, comprising the following components; in vitro biofidelic drop tower testing of a proximal femur; high-speed image analysis through digital image correlation; detailed accounting of the energy present during the drop tower test; organ level finite element simulations of the drop tower test; micro level finite element simulations of critical volumes of interest in the trabecular bone. Fracture in the femoral specimen initiated in the superior part of the neck. Measured fracture load was 3760N, compared to 4871N predicted based on the finite element analysis. Digital image correlation showed compressive surface strains as high as 7.1% prior to fracture. Voxel level results were consistent with high-speed video data and helped identify hidden local structural weaknesses. We found using a drop tower test protocol that a femoral neck fracture can be created with a fall velocity and energy representative of a sideways fall from standing. Additionally, we found that the nested explicit finite element method used allowed us to identify local structural weaknesses associated with femur fracture initiation.

Comparison of explicit finite element and mechanical simulation of the proximal femur during dynamic drop-tower testing

Current screening techniques based on areal bone mineral density (aBMD) measurements are unable to identify the majority of people who sustain hip fractures. Biomechanical examination of such events may help determine what predisposes a hip to be susceptible to fracture. Recently, drop-tower simulations of in-vitro sideways falls have allowed the study of the mechanical response of the proximal human femur at realistic impact speeds. This technique has created an opportunity to validate explicit finite element (FE) models against dynamic test data. This study compared the outcomes of 15 human femoral specimens fractured using a drop tower with complementary specimen-specific explicit FE analysis. Correlation coefficient and root mean square error (RMSE) were found to be moderate for whole bone stiffness comparison (R(2)=0.3476 and 22.85% respectively). No correlation was found between experimentally and computationally predicted peak force, however, energy absorption comparison produced moderate correlation and RMSE (R(2)=0.4781 and 29.14% respectively). By comparing predicted strain maps to high speed video data we demonstrated the ability of the FE models to detect vulnerable portions of the bones. Based on our observations, we conclude that there exists a need to extend the current apparent level material models for bone to cover higher strain rates than previously tested experimentally.

Interactive graph-cut segmentation for fast creation of finite element models from clinical ct data for hip fracture prediction

Abstract In this study, we propose interactive graph cut image segmentation for fast creation of femur finite element (FE) models from clinical computed tomography scans for hip fracture prediction. Using a sample of N = 48 bone scans representing normal, osteopenic and osteoporotic subjects, the proximal femur was segmented using manual (gold standard) and graph cut segmentation. Segmentations were subsequently used to generate FE models to calculate overall stiffness and peak force in a sideways fall simulations. Results show that, comparable FE results can be obtained with the graph cut method, with a reduction from 20 to 2–5 min interaction time. Average differences between segmentation methods of 0.22 mm were not significantly correlated with differences in FE derived stiffness (R2 = 0.08, p = 0.05) and weakly correlated to differences in FE derived peak force (R2 = 0.16, p = 0.01). We further found that changes in automatically assigned boundary conditions as a consequence of small segmentation differences were significantly correlated with FE derived results. The proposed interactive graph cut segmentation software MITK-GEM is freely available online at https://simtk.org/home/mitk-gem.

The compressive modulus and strength of saturated calcium sulphate dihydrate cements: Implications for testing standards

Calcium sulphate-based bone cement is a bone filler with proven biological advantages including biodegradability, biocompatibility and osteoconductivity. Mechanical properties of such brittle ceramic cements are frequently determined using the testing standard designed for ductile acrylic cements. The aims of the study were (1) to validate the suitability of this common testing protocol using saturated calcium sulphate dihydrate (CSD), and (2) to compare the strength and effective modulus of non-saturated and saturated CSD, in order to determine the changes in the mechanical behavior of CSD upon saturation. Unconfined compression tests to failure were performed on 190 cylindrical CSD samples. The samples were divided into four groups having different saturation levels (saturated, non-saturated) and end conditions (capped and non-capped). Two effective moduli were calculated per sample, based on the deformations measured using the machine platens and a sample-mounted extensometer. The effective moduli of non-saturated groups were found to be independent of the end conditions. The saturated and capped group showed no difference in the effective moduli derived from different measurement methods, while the saturated and non-capped group showed a significant difference between the machine platen- and extensometer-derived moduli. Strength and modulus values were significantly lower for saturated samples. It was assumed that the existence of water in saturated CSD alters the mechanical response of the material due to the changes in chemical and physical behaviors. These factors are considered to play important roles to decrease the shear strength of CSD. It was proposed that the reduction in CSD shear strength evokes local deformation at the platen-sample boundary, affecting the strength and effective moduli derived from the experiments. The results of this study highlighted the importance of appropriate and consistent testing methods when determining the mechanical properties of saturated ceramic cements.

The Effectiveness of Percutaneous Vertebroplasty Is Determined by the Patient-Specific Bone Condition and the Treatment Strategy

Purpose Vertebral fragility fractures are often treated by injecting bone cement into the collapsed vertebral bodies (vertebroplasty). The mechanisms by which vertebroplasty induces pain relief are not completely understood yet and recent debates cast doubt over the outcome of the procedure. The controversy is intensified by inconsistent results of randomized clinical trials and biomechanical studies that have investigated the effectiveness or the change in biomechanical response due to the reinforcement. The purpose of this study was to evaluate the effectiveness of vertebroplasty, by varying the relevant treatment parameters and (a) computationally predicting the improvement of the fracture risk depending on the chosen treatment strategy, and (b) identifying the determinants of a successful treatment. Methods A Finite Element model with a patient-specific failure criterion and direct simulation of PMMA infiltration in four lumbar vertebrae was used to assess the condition of the bone under compressive load before and after the virtual treatment, simulating in a total of 12000 virtual treatments. Results The results showed that vertebroplasty is capable of reducing the fracture risk by magnitudes, but can also have a detrimental effect. Effectiveness was strongly influenced by interactions between local bone quality, cement volume and injection location. However, only a moderate number of the investigated treatment strategies were able to achieve the necessary improvement for preventing a fracture. Conclusions We conclude that the effectiveness of vertebroplasty is sensitive to the patient’s condition and the treatment strategy.

Compressive mechanical properties and cytocompatibility of bone-compliant, linoleic acid-modified bone cement in a bovine model

Adjacent vertebral fractures are a common complication experienced by osteoporosis patients shortly after vertebroplasty. Whether these fractures are due to the bone cement properties, the cement filling characteristics or to the natural course of the disease is still unclear. However, some data suggests that such fractures might occur because of an imbalance in the load distribution due to a mismatch between the elastic modulus (E) of the bone-cement composite, and that of the vertebral cancellous bone. In this study, the properties of bone-compliant linoleic acid-modified bone cements were assessed using a bovine vertebroplasty model. Two groups of specimens (cement-only and bone-cement composites), and four subgroups comprising bone cements with elastic moduli in the range of 870-3500MPa were tested to failure in uniaxial compression. In addition, monomer release as well as time and concentration-dependent cytocompatibility was assessed through the cement extracts using a Saos-2 cell model. Composites augmented with bone-compliant cements exhibited a reduction in E despite their relatively high bone volume fraction (BVF). Moreover, a significant positive correlation between the BVF and the E for the composites augmented with 870MPa modulus cements was found. This was attributed to the increased relative contribution of the bone to the mechanical properties of the composites with a decrease in E of the bone cement. The use of linoleic acid reduced monomer conversion resulting in six times more monomer released after 24h. However, the cytocompatibility of the bone-compliant cements was comparable to that of the unmodified cements after the extracts were diluted four times. This study represents an important step towards introducing viable bone-compliant bone cements into vertebroplasty practice.

The influence of the modulus–density relationship and the material mapping method on the simulated mechanical response of the proximal femur in side-ways fall loading configuration

Contributing to slow advance of finite element (FE) simulations for hip fracture risk prediction, into clinical practice, could be a lack of consensus in the biomechanics community on how to map properties to the models. Thus, the aim of the present study was first, to systematically quantify the influence of the modulus-density relationship (E-ρ) and the material mapping method (MMM) on the predicted mechanical response of the proximal femur in a side-ways fall (SWF) loading configuration and second, to perform a model-to-model comparison of the predicted mechanical response within the femoral neck for all the specimens tested in the present study, using three different modelling techniques that have yielded good validation outcome in terms of surface strain prediction and whole bone response according to the literature. We found the outcome to be highly dependent on both the E-ρ relationship and the MMM. In addition, we found that the three modelling techniques that have resulted in good validation outcome in the literature yielded different principal strain prediction both on the surface as well as internally in the femoral neck region of the specimens modelled in the present study. We conclude that there exists a need to carry out a more comprehensive validation study for the SWF loading mode to identify which combination of MMMs and E-ρ relationship leads to the best match for whole bone and local mechanical response. The MMMs tested in the present study have been made publicly available at https://simtk.org/home/mitk-gem.

Influence of cement stiffness and bone morphology on the compressive properties of bone–cement composites in simulated vertebroplasty†

Vertebroplasty is widely used to treat vertebral compression fractures. Little is known about the influence of morphological parameters of the bone on the properties of the bone-cement composite. Furthermore, although generic finite element (FE) models have been suggested as a way to compute the values of these properties, their accuracy has not been established. In the experimental part of this study, we tested bovine cancellous bone and three different polymethylmethacrylate bone cements and determined six quasi-static uniaxial compressive properties of bone-cement composite specimens and 10 morphological parameters of the bone. For the FE work, we used two simulations, one being μFE and the other unit cell FE. In conclusion, we found that (1) for composite specimens, that relative contribution of the cement to the overall response of the composite increases with increasing cement stiffness; (2) the anisotropy ratio is the bone morphological property that exerts the most significant influence on the experimentally obtained compressive properties of the bone-cement composites determined; (3) the accuracy of the computed compressive properties of the composites ranged from low to high, depending on simulation method used. The largest errors, however, can partially be explained by difference in boundary conditions between the experimental testing and the simulation techniques. The lattermost finding points to the potential for simplified FE models being incorporated into automatic material mapping schemes in whole bone vertebra FE simulations.

Morphology based anisotropic finite element models of the proximal femur validated with experimental data

Finite element analysis (FEA) of bones scanned with Quantitative Computed Tomography (QCT) can improve early detection of osteoporosis. The accuracy of these models partially depends on the assigned material properties, but anisotropy of the trabecular bone cannot be fully captured due to insufficient resolution of QCT. The inclusion of anisotropy measured from high resolution peripheral QCT (HR-pQCT) could potentially improve QCT-based FEA of the femur, although no improvements have yet been demonstrated in previous experimental studies. This study analyzed the effects of adding anisotropy to clinical resolution femur models by constructing six sets of FE models (two isotropic and four anisotropic) for each specimen from a set of sixteen femurs that were experimentally tested in sideways fall loading with a strain gauge on the superior femoral neck. Two different modulus-density relationships were tested, both with and without anisotropy derived from mean intercept length analysis of HR-pQCT scans. Comparing iso- and anisotropic models to the experimental data resulted in nearly identical correlation and highly similar linear regressions for both whole bone stiffness and strain gauge measurements. Anisotropic models contained consistently greater principal compressive strains, approximately 14% in magnitude, in certain internal elements located in the femoral neck, greater trochanter, and femoral head. In summary, anisotropy had minimal impact on macroscopic measurements, but did alter internal strain behavior. This suggests that organ level QCT-based FE models measuring femoral stiffness have little to gain from the addition of anisotropy, but studies considering failure of internal structures should consider including anisotropy to their models.

Material mapping strategy to improve the predicted response of the proximal femur to a sideways fall impact

Sideways falls are largely responsible for the highly prevalent osteoporotic hip fractures in todays society. These injuries are dynamic events, therefore dynamic FE models validated with dynamic ex vivo experiments provide a more realistic simulation than simple quasi-static analysis. Drop tower experiments using cadaveric specimens were used to identify the material mapping strategy that provided the most realistic mechanical response under impact loading. The present study tested the addition of compression-tension asymmetry, tensile bone damage, and cortical-specific strain rate dependency to the material mapping strategy of fifteen dynamic FE models of the proximal femur, and found improved correlations and reduced error for whole bone stiffness (R2 = 0.54, RSME = 0.87kN/mm) and absolute maximum force (R2 = 0.56, RSME =0.57kN), and a high correlation in impulse response (R2 = 0.82, RSME =12.38kg/s). Simulations using fully bonded nodes between the rigid bottom plate and PMMA cap supporting the femoral head had higher correlations and less error than simulations using a frictionless sliding at this contact surface. Strain rates over 100/s were observed in certain elements in the femoral neck and trochanter, indicating that additional research is required to better quantify the strain rate dependencies of both trabecular and cortical bone at these strain rates. These results represent the current benchmark in dynamic FE modeling of the proximal femur in sideways falls. Future work should also investigate improvements in experimental validation techniques by developing better displacement measurements and by enhancing the biofidelity of the impact loading wherever possible.