Sami P. Väänänen

Experimental Validation Of Finite Element Model For Proximal Composite Femur Using Optical Measurements

Patient-specific finite element models have been used to predict femur strength and fracture risk in individuals. Validation of the adopted finite element modelling procedure against mechanical testing data is a crucial step when aiming for clinical applications. The majority of the works available in literature used data from strain gages to validate the model, thus having up to 15 experimental measurements. Optical techniques, such as digital image correlation, can help to improve the models by providing a continuous field of deformation data over a femoral surface. The main objective of this study was to validate finite element models of six composite femora against strain data from digital image correlation, obtained during fracture tests performed in quasi-axial loading configuration. The finite element models were obtained from CT scans, by means of a semi-automatic segmentation. The principal strains both during the elastic phase and close to the fracture were compared, and showed a correlation coefficient close to 0.9. In the linear region, the slope and intercept were close to zero and unity, while for the case when fracture load was simulated, the slope decreased somewhat. The accuracy of the obtained results is comparable with the state-of-the-art literature, with the significant improvement of having around 50,000 data points for each femur. This large number of measurements allows a more comprehensive validation of the predictions by the finite element models, since thousand of points are tracked along the femoral neck and trochanter region, i.e., the sites that are most critical for femur fracture. Moreover, strain measurement biases due to the strain gage reinforcement effect, were avoided. The combined experimental-numerical approach proved to be ready for application to in-vitro tests of human cadaver femurs, thus helping to develop a suitable mechanistic fracture risk criterion.

Generation of 3D shape, density, cortical thickness and finite element mesh of proximal femur from a DXA image

Areal bone mineral density (aBMD), as measured by dual-energy X-ray absorptiometry (DXA), predicts hip fracture risk only moderately. Simulation of bone mechanics based on DXA imaging of the proximal femur, may help to improve the prediction accuracy. Therefore, we collected three (1-3) image sets, including CT images and DXA images of 34 proximal cadaver femurs (set 1, including 30 males, 4 females), 35 clinical patient CT images of the hip (set 2, including 27 males, 8 females) and both CT and DXA images of clinical patients (set 3, including 12 female patients). All CT images were segmented manually and landmarks were placed on both femurs and pelvises. Two separate statistical appearance models (SAMs) were built using the CT images of the femurs and pelvises in sets 1 and 2, respectively. The 3D shape of the femur was reconstructed from the DXA image by matching the SAMs with the DXA images. The orientation and modes of variation of the SAMs were adjusted to minimize the sum of the absolute differences between the projection of the SAMs and a DXA image. The mesh quality and the location of the SAMs with respect to the manually placed control points on the DXA image were used as additional constraints. Then, finite element (FE) models were built from the reconstructed shapes. Mean point-to-surface distance between the reconstructed shape and CT image was 1.0 mm for cadaver femurs in set 1 (leave-one-out test) and 1.4 mm for clinical subjects in set 3. The reconstructed volumetric BMD showed a mean absolute difference of 140 and 185 mg/cm(3) for set 1 and set 3 respectively. The generation of the SAM and the limitation of using only one 2D image were found to be the most significant sources of errors in the shape reconstruction. The noise in the DXA images had only small effect on the accuracy of the shape reconstruction. DXA-based FE simulation was able to explain 85% of the CT-predicted strength of the femur in stance loading. The present method can be used to accurately reconstruct the 3D shape and internal density of the femur from 2D DXA images. This may help to derive new information from clinical DXA images by producing patient-specific FE models for mechanical simulation of femoral bone mechanics.

Repeatability of digital image correlation for measurement of surface strains in composite long bones

Digital image correlation (DIC) can measure full-field surface strains during mechanical testing of hard and soft tissues. When compared to traditional methods, such as strain gauges, DIC offers larger validation data (∼50,000 points) for, e.g., finite element models. Our main aim was to evaluate the repeatability of surface strain measurements with DIC during compressive testing of composite femurs mimicking human bones. We also studied the similarity of the composite femur samples using CT. Composite femurs were chosen as test material to minimize the uncertainties associated with the use of cadaveric tissues and to understand the variability of the DIC measurement itself. Six medium-sized fourth generation composite human proximal femora (Sawbones) were CT imaged and mechanically tested in stance configuration. The force-displacement curves were recorded and the 3D surface strains were measured with DIC on the anterior surface of the femurs. Five femurs fractured at the neck-trochanter junction and one at the site below the minor trochanter. CT image of this bone showed an air cavity at the initial fracture site. All femurs fractured through a sudden brittle crack. The fracture force for the composite bones was 5751±650N (mean±SD). The maximum von Mises strain during the fractures was 2.4±0.8%. Noise in one experiment was 5-30µε. When applied loads were equalized the variation in strains between the bones was 20-25%, and when the maximum strains were equalized, variation in the other regions was 5-10%. DIC showed that the ability of nominally identical composite bones to bear high strains and loads before fracturing may vary between the samples.

How accurately can subject-specific finite element models predict strains and strength of human femora? Investigation using full-field measurements.

Subject-specific finite element models have been proposed as a tool to improve fracture risk assessment in individuals. A thorough laboratory validation against experimental data is required before introducing such models in clinical practice. Results from digital image correlation can provide full-field strain distribution over the specimen surface during in vitro test, instead of at a few pre-defined locations as with strain gauges. The aim of this study was to validate finite element models of human femora against experimental data from three cadaver femora, both in terms of femoral strength and of the full-field strain distribution collected with digital image correlation. The results showed a high accuracy between predicted and measured principal strains (R(2)=0.93, RMSE=10%, 1600 validated data points per specimen). Femoral strength was predicted using a rate dependent material model with specific strain limit values for yield and failure. This provided an accurate prediction (<2% error) for two out of three specimens. In the third specimen, an accidental change in the boundary conditions occurred during the experiment, which compromised the femoral strength validation. The achieved strain accuracy was comparable to that obtained in state-of-the-art studies which validated their prediction accuracy against 10-16 strain gauge measurements. Fracture force was accurately predicted, with the predicted failure location being very close to the experimental fracture rim. Despite the low sample size and the single loading condition tested, the present combined numerical-experimental method showed that finite element models can predict femoral strength by providing a thorough description of the local bone mechanical response.

Full-Field Strain Measurement During Mechanical Testing of the Human Femur at Physiologically Relevant Strain Rates

Understanding the mechanical properties of human femora is of great importance for the development of a reliable fracture criterion aimed at assessing fracture risk. Earlier ex vivo studies have been conducted by measuring strains on a limited set of locations using strain gauges (SGs). Digital image correlation (DIC) could instead be used to reconstruct the full-field strain pattern over the surface of the femur. The objective of this study was to measure the full-field strain response of cadaver femora tested at a physiological strain rate up to fracture in a configuration resembling single stance. The three cadaver femora were cleaned from soft tissues, and a white background paint was applied with a random black speckle pattern over the anterior surface. The mechanical tests were conducted up to fracture at a constant displacement rate of 15 mm/s, and two cameras recorded the event at 3000 frames per second. DIC was performed to retrieve the full-field displacement map, from which strains were derived. A low-pass filter was applied over the measured displacements before the crack opened in order to reduce the noise level. The noise levels were assessed using a dedicated control plate. Conversely, no filtering was applied at the frames close to fracture to get the maximum resolution. The specimens showed a linear behavior of the principal strains with respect to the applied force up to fracture. The strain rate was comparable to the values available in literature from in vivo measurements during daily activities. The cracks opened and fully propagated in less than 1 ms, and small regions with high values of the major principal strains could be spotted just a few frames before the crack opened. This corroborates the hypothesis of a strain-driven fracture mechanism in human bone. The data represent a comprehensive collection of full-field strains, both at physiological load levels and up to fracture. About 10,000 points were tracked on each bone, providing superior spatial resolution compared to ∼15 measurements typically collected using SGs. These experimental data collection can be further used for validation of numerical models, and for experimental verification of bone constitutive laws and fracture criteria.

Effect of bone inhomogeneity on tibiofemoral contact mechanics during physiological loading

It is not known how inhomogeneous mechanical properties of bone affect contact mechanics and cartilage response during physiological loading of the knee joint. In this study, a finite element model of a cadaver knee joint was constructed based on quantitative computed tomography (QCT). The mechanical properties of bone were altered and their effect on tibiofemoral contact mechanics and cartilage stresses, strains and pore pressures were evaluated during the first 20% of stance. For this purpose, models with rigid, homogeneous and inhomogeneous bones were created. When bone was modeled to be rigid, the resulting contact pressures were substantially higher in the medial side of the joint, as compared to the non-rigid bones. Similar changes were revealed also in stresses, strains and pore pressures throughout the cartilage depth at the cartilage-cartilage contact area. Furthermore, the mechanical response of medial tibial cartilage was found to be highly dependent on the bone properties. When Young׳s modulus in the model with homogeneous bone was 5GPa, cartilage mechanical response approached to that of the model with inhomogeneous bone. Finally, when the apparent bone mineral densities were decreased globally in the inhomogeneous bone, stresses, strains and pore pressures were decreased at all layers of medial tibial cartilage. Similar changes were observed also in cartilage-cartilage contact area of the lateral compartment but with a lesser extent. These results indicate that during physiological loading Young׳s modulus of bone has a substantial influence on cartilage stresses and strains, especially in the medial compartment.

Estimation of 3D rotation of femur in 2D hip radiographs

Femoral radiographs are affected by the degree of rotation of the femur with respect to the plane of projection. We aimed to determine the 3D rotation of the proximal femur in 2D radiographs. A 3D Statistical Appearance Model (SAM), which was built from CT images of cadaver proximal femurs (n=33) was randomly sampled to form a training set of 500 bones. Nineteen clinical CT images were collected for testing. All CT images were rotated to ±20° in 2° division around the shaft axis, ±10° around medial-lateral axis, and by simultaneous rotation of both axes (±16° and ±8° around shaft and medial-lateral axes). In each orientation, a 2D projection was recorded for generating a 2D SAM. The outcome parameters of the 2D SAM were used as input for a linear regression model and an artificial neural network to predict the rotation. The artificial neural network estimated the rotation more accurately than the linear regression. For artificial neural networks the mean errors were 4.0° and 2.0° around the shaft and medial-lateral axes, respectively. For an individual radiograph, the confidence interval of estimation was still relatively large. However, this method has high potential to differentiate the amount of rotations in two image sets.

Semi-automated International Cartilage Repair Society scoring of equine articular cartilage lesions in optical coherence tomography images

BACKGROUND Arthroscopic optical coherence tomography (OCT) is a promising tool for the detailed evaluation of articular cartilage injuries. However, OCT-based articular cartilage scoring still relies on the operators visual estimation. OBJECTIVES To test the hypothesis that semi-automated International Cartilage Repair Society (ICRS) scoring of chondral lesions seen in OCT images could enhance intra- and interobserver agreement of scoring and its accuracy. STUDY DESIGN Validation study using equine cadaver tissue. METHODS Osteochondral samples (n = 99) were prepared from 18 equine metacarpophalangeal joints and imaged using OCT. Custom-made software was developed for semi-automated ICRS scoring of cartilage lesions on OCT images. Scoring was performed visually and semi-automatically by five observers, and levels of inter- and intraobserver agreement were calculated. Subsequently, OCT-based scores were compared with ICRS scores based on light microscopy images of the histological sections of matching locations (n = 82). RESULTS When semi-automated scoring of the OCT images was performed by multiple observers, mean levels of intraobserver and interobserver agreement were higher than those achieved with visual OCT scoring (83% vs. 77% and 74% vs. 33%, respectively). Histology-based scores from matching regions of interest agreed better with visual OCT-based scoring than with semi-automated OCT scoring; however, the accuracy of the software was improved by optimising the threshold combinations used to determine the ICRS score. MAIN LIMITATIONS Images were obtained from cadavers. CONCLUSIONS Semi-automated scoring software improved the reproducibility of ICRS scoring of chondral lesions in OCT images and made scoring less observer-dependent. The image analysis and segmentation techniques adopted in this study warrant further optimisation to achieve better accuracy with semi-automated ICRS scoring. In addition, studies on in vivo applications are required.

Tissue viscoelasticity is related to tissue composition but may not fully predict the apparent-level viscoelasticity in human trabecular bone – An experimental and finite element study

Trabecular bone is viscoelastic under dynamic loading. However, it is unclear how tissue viscoelasticity controls viscoelasticity at the apparent-level. In this study, viscoelasticity of cylindrical human trabecular bone samples (n=11, male, age 18-78 years) from 11 proximal femurs were characterized using dynamic and stress-relaxation testing at the apparent-level and with creep nanoindentation at the tissue-level. In addition, bone tissue elasticity was determined using scanning acoustic microscope (SAM). Tissue composition and collagen crosslinks were assessed using Raman micro-spectroscopy and high performance liquid chromatography (HPLC), respectively. Values of material parameters were obtained from finite element (FE) models by optimizing tissue-level creep and apparent-level stress-relaxation to experimental nanoindentation and unconfined compression testing values, respectively, utilizing the second order Prony series to depict viscoelasticity. FE simulations showed that tissue-level equilibrium elastic modulus (Eeq) increased with increasing crystallinity (r=0.730, p=.011) while at the apparent-level it increased with increasing hydroxylysyl pyridinoline content (r=0.718, p=.019). In addition, the normalized shear modulus g1 (r=-0.780, p=.005) decreased with increasing collagen ratio (amide III/CH2) at the tissue-level, but increased (r=0.696, p=.025) with increasing collagen ratio at the apparent-level. No significant relations were found between the measured or simulated viscoelastic parameters at the tissue- and apparent-levels nor were the parameters related to tissue elasticity determined with SAM. However, only Eeq, g2 and relaxation time τ1 from simulated viscoelastic values were statistically different between tissue- and apparent-levels (p<.01). These findings indicate that bone tissue viscoelasticity is affected by tissue composition but may not fully predict the macroscale viscoelasticity in human trabecular bone.

FUNCTIONAL IMAGING OF THE PROXIMAL FEMUR USING SHAPE TEMPLATE AND A BONE MINERAL DENSITY IMAGE

Measurement of bone mineral density (BMD) by DXA (dual-energy X-ray absorptiometry) is generally considered to be the clinical gold standard to diagnose osteoporosis. However, BMD alone is only a moderate predictor of fracture risk. Finite element analyses (FEA) of bone mechanics can contribute to a more accurate prediction of fracture risk (Cody et al. 1999). However, CT imaging is relatively expensive and inflicts larger radiation doses on the patient.Copyright