Improving the validation of finite element models with quantitative full-field strain comparisons

The techniques used to validate finite element (FE) models against experimental results have changed little during the last decades, even though the traditional approach of using single point measurements from strain gauges has major limitations: the strain distribution across the surface is not captured and the accurate determination of strain gauge positions on the model surface is difficult if the 3D surface topography of the bone surface is not measured. The full-field strain measurement technique of digital speckle pattern interferometry (DSPI) can overcome these problems, but the potential of this technique has not yet been fully exploited in validation studies. Here we explore new ways of quantifying and visualising the variation in strain magnitudes and orientations within and between repeated DSPI measurements as well as between the DSPI measurements and FEA results. We show that our approach provides a much more comprehensive and accurate validation than traditional methods. The measurement repeatability and the correspondence between measured and predicted strains vary to a great degree within and between measurement areas. The two models used in this study predict the measured strain directions and magnitudes surprisingly well considering that homogeneous and isotropic mechanical properties were assigned to the models. However, the full-field comparisons also reveal some discrepancies between measured and predicted strains that are most probably caused by inaccuracies in the models' geometries and the degree of simplification of the modelled material properties.     

Subject-specific finite element model of the pelvis: development, validation and sensitivity studies

A better understanding of the three-dimensional mechanics of the pelvis, at the patient-specific level, may lead to improved treatment modalities. Although finite element (FE) models of the pelvis have been developed, validation by direct comparison with subject-specific strains has not been performed, and previous models used simplifying assumptions regarding geometry and material properties. The objectives of this study were to develop and validate a realistic FE model of the pelvis using subject-specific estimates of bone geometry, location-dependent cortical thickness and trabecular bone elastic modulus, and to assess the sensitivity of FE strain predictions to assumptions regarding cortical bone thickness as well as bone and cartilage material properties. A FE model of a cadaveric pelvis was created using subject-specific computed tomography image data. Acetabular loading was applied to the same pelvis using a prosthetic femoral stem in a fashion that could be easily duplicated in the computational model. Cortical bone strains were monitored with rosette strain gauges in ten locations on the left hemipelvis. FE strain predictions were compared directly with experimental results for validation. Overall, baseline FE predictions were strongly correlated with experimental results (r2=0.824), with a best-fit line that was not statistically different than the line y=x (experimental strains = FE predicted strains). Changes to cortical bone thickness and elastic modulus had the largest effect on cortical bone strains. The FE model was less sensitive to changes in all other parameters. The methods developed and validated in this study will be useful for creating and analyzing patient-specific FE models to better understand the biomechanics of the pelvis.     

Experimental validation of a finite element model of the proximal femur using digital image correlation and a composite bone model

Computational biomechanical models are useful tools for supporting orthopedic implant design and surgical decision making, but because they are a simplification of the clinical scenario they must be carefully validated to ensure that they are still representative. The goal of this study was to assess the validity of the generation process of a structural finite element model of the proximal femur employing the digital image correlation (DIC) strain measurement technique. A finite element analysis model of the proximal femur subjected to gait loading was generated from a CT scan of an analog composite femur, and its predicted mechanical behavior was compared with an experimental model. Whereas previous studies have employed strain gauging to obtain discreet point data for validation, in this study DIC was used for full field quantified comparison of the predicted and experimentally measured strains. The strain predicted by the computational model was in good agreement with experimental measurements, with R(2) correlation values from 0.83 to 0.92 between the simulation and the tests. The sensitivity and repeatability of the strain measurements were comparable to or better than values reported in the literature for other DIC tests on tissue specimens. The experimental-model correlation was in the same range as values obtained from strain gauging, but the DIC technique produced more detailed, full field data and is potentially easier to use. As such, the findings supported the validity of the model generation process, giving greater confidence in the model's predictions, and digital image correlation was demonstrated as a useful tool for the validation of biomechanical models.     

The effect of the density-modulus relationship selected to apply material properties in a finite element model of long bone

Material property assignment is a critical step in developing subject-specific finite element models of bone. Inhomogeneous material properties are often applied using an equation relating density and elastic modulus, with the density information coming from CT scans of the bone. Very few previous studies have investigated which density-elastic modulus relationships from the literature are most suitable for application in long bone. No such studies have been completed for the ulna. The purpose of this study was to investigate six such density-modulus relationships and compare the results to experimental strains from eight cadaveric ulnae. Subject-specific finite element models were developed for each bone using micro-CT scans. Six density-modulus equations were trialed in each bone, resulting in a total of 48 models. Data from a previously completed experimental study in which each bone was instrumented with twelve strain gauges were used for comparison. Although the relationship that best matched experimental strains was somewhat specimen and location dependent, there were two relations which consistently matched the experimental strains most closely. One of these under-estimated and one over-estimated the experimental strain values, by averages of 15% and 31%, respectively. The results of this study suggest that the ideal relationship for the ulna may lie somewhere in between these two relations.     

Static vs dynamic loads as an influence on bone remodelling

Remodelling activity in the avian ulna was assessed under conditions of disuse alone, disuse with a superimposed continuous compressive load, and disuse interrupted by a short daily period of intermittent loading. The ulnar preparation consisted of the 110 mm section of the bone shaft between two submetaphyseal osteotomies. Each end of the preparation was transfixed by a stainless steel pin and the shaft either protected from normal functional loading with the pins joined by external fixators, loaded continuously in compression by joining the pins with springs, or loaded intermittently in compression for a single 100 s period per day by engaging the pins in an Instron machine. Similar loads (525 N) were used in both static and dynamic cases. The strains engendered were determined by strain gauges, and at their maximum around the bone's midshaft were − 0.002. The intermittent load was applied at a frequency of 1 Hz as a ramped square wave, with a rate of change of strain during the ramp of 0.01 s⁻¹. Peak strain at the midshaft of the ulna during wing flapping in the intact bone was recorded from bone bonded strain gauges in vivo as −0.0033 with a maximum rate of change of strain of 0.056 s⁻¹.

Examination of bone sections from the midpoint of the preparation after an 8 week period indicated that in both non-loaded and statically loaded bones there was an increase in both endosteal diameter and intra cortical porosity. These changes produced a decrease in cross sectional area which was similar in the two groups (− 13%). Intermittently loaded bones however showed a 24% increase in cross sectional area resulting from new bone deposited predominantly, but not exclusively, on the periosteal surface. It appears that in this preparation, a static load sufficient to produce strains in the functional dynamic strain range has no effect on bone remodelling, whereas a similar load applied intermittently for a short daily period may be associated with a substantial increase in bone mass.     

Finite element and experimental cortex strains of the intact and implanted tibia

Finite Element (FE) models for the simulation of intact and implanted bone find their main purpose in accurately reproducing the associated mechanical behavior. FE models can be used for preclinical testing of joint replacement implants, where some biomechanical aspects are difficult, if not possible, to simulate and investigate in vitro. To predict mechanical failure or damage, the models should accurately predict stresses and strains. Commercially available synthetic femur models have been extensively used to validate finite element models, but despite the vast literature available on the characteristics of synthetic tibia, numerical and experimental validation of the intact and implant assemblies of tibia are very limited or lacking. In the current study, four FE models of synthetic tibia, intact and reconstructed, were compared against experimental bone strain data, and an overall agreement within 10% between experimental and FE strains was obtained. Finite element and experimental (strain gauge) models of intact and implanted synthetic tibia were validated based on the comparison of cortex bone strains. The study also includes the analysis carried out on standard tibial components with cemented and noncemented stems of the P.F.C Sigma Modular Knee System. The overall agreement within 10% previously established was achieved, indicating that FE models could be successfully validated. The obtained results include a statistical analysis where the root-mean-square-error values were always <10%. FE models can successfully reproduce bone strains under most relevant acting loads upon the condylar surface of the tibia. Moreover, FE models, once properly validated, can be used for preclinical testing of tibial knee replacement, including misalignment of the implants in the proximal tibia after surgery, simulation of long-term failure according to the damage accumulation failure scenario, and other related biomechanical aspects.     

Subject-specific finite element models implementing a maximum principal strain criterion are able to estimate failure risk and fracture location on human femurs tested in vitro

No agreement on the choice of the failure criterion to adopt for the bone tissue can be found in the literature among the finite element studies aiming at predicting fracture risk of bones. The use of stress-based criteria seems to prevail on strain-based ones, while basic bone biomechanics suggest using strain parameters to describe failure. The aim of the present combined experimental-numerical study was to verify, using subject-specific finite element models able to accurately predict strains, if a strain-based failure criterion could identify the failure patterns of bones. Three cadaver femurs were CT-scanned and subsequently fractured in a clinically relevant single-stance loading scenario. Load-displacement curves and high-speed movies were acquired to define the failure load and the location of fracture onset, respectively. Subject-specific finite element models of the three femurs were built from CT data following a validated procedure. A maximum principal strain criterion was implemented in the finite element models, and two stress-based criteria selected for comparison. The failure loads measured were applied to the models, and the computed risks of fracture were compared to the results of the experimental tests. The proposed principal strain criterion managed to correctly identify the level of failure risk and the location of fracture onset in all the modelled specimens, while Von Mises or maximum principal stress criteria did not give significant information. A maximum principal strain criterion can thus be defined a suitable candidate for the in vivo risk factor assessment on long bones.     

Comparison of the linear finite element prediction of deformation and strain of human cancellous bone to 3D digital volume correlation measurements

The mechanical properties of cancellous bone and the biological response of the tissue to mechanical loading are related to deformation and strain in the trabeculae during function. Due to the small size of trabeculae, their motion is difficult to measure. To avoid the need to measure trabecular motions during loading the finite element method has been used to estimate trabecular level mechanical deformation. This analytical approach has been empirically successful in that the analytical models are solvable and their results correlate with the macroscopically measured stiffness and strength of bones. The present work is a direct comparison of finite element predictions to measurements of the deformation and strain at near trabecular level. Using the method of digital volume correlation, we measured the deformation and calculated the strain at a resolution approaching the trabecular level for cancellous bone specimens loaded in uniaxial compression. Smoothed results from linearly elastic finite element models of the same mechanical tests were correlated to the empirical three-dimensional (3D) deformation in the direction of loading with a coefficient of determination as high as 97% and a slope of the prediction near one. However, real deformations in the directions perpendicular to the loading direction were not as well predicted by the analytical models. Our results show, that the finite element modeling of the internal deformation and strain in cancellous bone can be accurate in one direction but that this does not ensure accuracy for all deformations and strains.     

3D characterization of bone strains in the rat tibia loading model

Bone strain is considered one of the factors inducing bone tissue response to loading. Nevertheless, where animal studies can provide detailed data on bone response, they only offer limited information on experimental bone strains. Including micro-CT-based finite element (micro FE) models in the analysis represents a potent methodology for quantifying strains in bone. Therefore, the main objective of this study was to develop and validate specimen-specific micro FE models for the assessment of bone strains in the rat tibia compression model. Eight rat limbs were subjected to axial compression loading; strain at the medio-proximal site of the tibiae was measured by means of strain gauges. Specimen-specific micro FE models were created and analyzed. Repeated measurements on each limb indicated that the effect of limb positioning was small (COV?= 6.45 ± 2.27 %). Instead, the difference in the measured strains between the animals was high (54.2%). The computational strains calculated at the strain gauge site highly correlated to the measured strains (R ²?=?0.95). Maximum peak strains calculated at exactly 25% of the tibia length for all specimens were equal to 435.11 ± 77.88 microstrains (COV?=?17.19%). In conclusion, we showed that strain gauge measurements are very sensitive to the exact strain gauge location on the bone; hence, the use of strain gauge data only is not recommended for studies that address at identifying reliable relationships between tissue response and local strains. Instead, specimen-specific micro FE models of rat tibiae provide accurate estimates of tissue-level strains.     

Computed-tomography-based finite-element models of long bones can accurately capture strain response to bending and torsion

Finite element (FE) models of long bones constructed from computed-tomography (CT) data are emerging as an invaluable tool in the field of bone biomechanics. However, the performance of such FE models is highly dependent on the accurate capture of geometry and appropriate assignment of material properties. In this study, a combined numerical-experimental study is performed comparing FE-predicted surface strains with strain-gauge measurements. Thirty-six major, cadaveric, long bones (humerus, radius, femur and tibia), which cover a wide range of bone sizes, were tested under three-point bending and torsion. The FE models were constructed from trans-axial volumetric CT scans, and the segmented bone images were corrected for partial-volume effects. The material properties (Young's modulus for cortex, density-modulus relationship for trabecular bone and Poisson's ratio) were calibrated by minimizing the error between experiments and simulations among all bones. The R(2) values of the measured strains versus load under three-point bending and torsion were 0.96-0.99 and 0.61-0.99, respectively, for all bones in our dataset. The errors of the calculated FE strains in comparison to those measured using strain gauges in the mechanical tests ranged from -6% to 7% under bending and from -37% to 19% under torsion. The observation of comparatively low errors and high correlations between the FE-predicted strains and the experimental strains, across the various types of bones and loading conditions (bending and torsion), validates our approach to bone segmentation and our choice of material properties.     

Using digital image correlation to determine bone surface strains during loading and after adaptation of the mouse tibia

Previous models of cortical bone adaptation, in which loading is imposed on the bone, have estimated the strains in the tissue using strain gauges, analytical beam theory, or finite element analysis. We used digital image correlation (DIC), tracing a speckle pattern on the surface of the bone during loading, to determine surface strains in a murine tibia during compressive loading through the knee joint. We examined whether these surface strains in the mouse tibia are modified following two weeks of load-induced adaptation by comparison with contralateral controls. Results indicated non-uniform strain patterns with isolated areas of high strain (0.5%), particularly on the medial side. Strain measurements were reproducible (standard deviation of the error 0.03%), similar between specimens, and in agreement with strain gauge measurements (between 0.1 and 0.2% strain). After structural adaptation, strains were more uniform across the tibial surface, particularly on the medial side where peak strains were reduced from 0.5% to 0.3%. Because DIC determines local strains over the entire surface, it will provide a better understanding of how strain stimulus influences the bone response during adaptation.     

Thermographic stress analysis in cortical bone.

Under adiabatic (or near adiabatic) conditions a volumetric change in an elastic material will produce a corresponding change in temperature. Based upon this principle, thermographic stress analysis (TSA) measures changes in surface heat flux (which are related to changes in surface temperature) and relates them to a coupled form of strains or stresses. To demonstrate the feasibility of using this technique for biomechanical applications, we thermographically measured heat flux from loaded specimens of cortical bone and correlated the results with strain gage data. Regular parallelepipeds were cut from the cortex of bovine femora and loaded sinusoidally at 20 Hz. At this rate of loading, mechanically induced changes in surface temperature could be sampled (via heat flux) prior to a measureable attenuation of the thermoelastic effect. Correlation coefficients demonstrated a significant linear relationship between TSA and measured and computed mechanical parameters (stress, strain, first strain invariant, and strain energy density). TSA therefore appears to be a promising technology for experimental stress analysis in cortical bone.     

A comparison of bone strain measurements at anatomically relevant sites using surface gauges versus strain gauged bone staples

Instrumented bone staples were first introduced as an alternative to surface-mounted strain gauges for use in human in vivo bone strain measurements because their fixation to bone is secure and requires not only minimally invasive surgery. Bench-top bone bending models have shown that the output from strain gauged bone staples compares favorably to that of traditional mounted gauges. However their within- and across-subject performance at sites typically instrumented in vivo has never been examined. This study used seven human cadaver lower extremities with an age range of 23-81 years old and a dynamic gait simulator to examine and compare axial strains in the mid tibial diaphysis and on the dorsal surface of the second metatarsal as measured simultaneously with strain gauged bone staples and with traditional surface-mounted gauges. Rosette configurations were used at the tibial site for deriving principal compression and tension, and shear strains. Axial outputs from the two gauge types demonstrated strong linear relationships for the tibia (r(2)=0.78-0.94) and the second metatarsal (r(2)=0.96-0.99), but coefficients (slopes) for the relationship were variable (range 7-20), across subjects and across sites. The apparent low reliability of strain gauged staples may be explained by the fact that both strain gauged staples and surface strain gauges are inexact to some degree, do not measure strains from exactly the same areas and strain gauged staples reflect surface strains as well as deformations within the cortex. There were no relationships for the principal tibia compression, tension or shear strain measurements derived from the two rosette gauge types, reflecting the very different anatomical areas measured by each of the constructs in this study. Strain gauged bone staples may be most useful in comparing relative axial intra-subject differences between activities, but inter-subject variability may require larger sample sizes to detect differences between populations.     

A round-robin finite element analysis of human femur mechanics between seven participating laboratories with experimental validation

Finite element analysis is a common tool that has been used for the past few decades to predict the mechanical behavior of bone. However, to our knowledge, there are no round-robin finite element analyses of long human bones with more than two participating biomechanics laboratories published yet, where the results of the experimental tests were not known in advance. We prepared a fresh-frozen human femur for a compression test in a universal testing machine measuring the strains at 10 bone locations as well as the deformation of the bone in terms of the displacement of the loading point at a load of 2?kN. The computed tomography data of the bone with a calibration phantom as well as the orientation of the bone in the testing machine with the according boundary conditions were delivered to seven participating laboratories. These were asked to perform a finite element analysis simulating the experimental setup and deliver their results to the coordinator without knowing the experimental results. Resultantly, four laboratories had deviations from the experimentally measured strains of less than 40%, and three laboratories had deviations of their numerically determined values compared to the experimental data of more than 120%. These deviations are thought to be based on different material laws and material data, as well as different material mapping methods. Investigations will be conducted to clarify and assess the reasons for the large deviations in the numerical data. It was shown that the precision of finite element models of the human femur is not yet as developed as desired by the biomechanics community.     

A full body musculoskeletal model based on flexible multibody simulation approach utilised in bone strain analysis during human locomotion

Load-induced strains applied to bone can stimulate its development and adaptation. In order to quantify the incident strains within the skeleton, in vivo implementation of strain gauges on the surfaces of bone is typically used. However, in vivo strain measurements require invasive methodology that is challenging and limited to certain regions of superficial bones only such as the anterior surface of the tibia. Based on our previous study [Al Nazer et al. (2008) J Biomech. 41:1036-1043], an alternative numerical approach to analyse in vivo strains based on the flexible multibody simulation approach was proposed. The purpose of this study was to extend the idea of using the flexible multibody approach in the analysis of bone strains during physical activity through integrating the magnetic resonance imaging (MRI) technique within the framework. In order to investigate the reliability and validity of the proposed approach, a three-dimensional full body musculoskeletal model with a flexible tibia was used as a demonstration example. The model was used in a forward dynamics simulation in order to predict the tibial strains during walking on a level exercise. The flexible tibial model was developed using the actual geometry of human tibia, which was obtained from three-dimensional reconstruction of MRI. Motion capture data obtained from walking at constant velocity were used to drive the model during the inverse dynamics simulation in order to teach the muscles to reproduce the motion in the forward dynamics simulation. Based on the agreement between the literature-based in vivo strain measurements and the simulated strain results, it can be concluded that the flexible multibody approach enables reasonable predictions of bone strain in response to dynamic loading. The information obtained from the present approach can be useful in clinical applications including devising exercises to prevent bone fragility or to accelerate fracture healing.     

Experimental modal analysis on fresh-frozen human hemipelvic bones employing a 3D laser vibrometer for the purpose of modal parameter identification

To provide a close-to-reality simulation model, such as for improved surgery planning, this model has to be experimentally verified. The present article describes the use of a 3D laser vibrometer for determining modal parameters of human pelvic bones that can be used for verifying a finite elements model. Compared to previously used sensors, such as acceleration sensors or strain gauges, the laser vibrometric procedure used here is a non-contact and non-interacting measuring method that allows a high density of measuring points and measurement in a global coordinate system. Relevant modal parameters were extracted from the measured data and provided for verifying the model. The use of the 3D laser vibrometer allowed the establishment of a process chain for experimental examination of the pelvic bones that was optimized with respect to time and effort involved. The transfer functions determined feature good signal quality. Furthermore, a comparison of the results obtained from pairs of pelvic bones showed that repeatable measurements can be obtained with the method used.     

Application of the digital volume correlation technique for the measurement of displacement and strain fields in bone: A literature review

Digital volume correlation (DVC) provides experimental measurements of displacements and strains throughout the interior of porous materials such as trabecular bone. It can provide full-field continuum- and tissue-level measurements, desirable for validation of finite element models, by comparing image volumes from subsequent µCT scans of a sample in unloaded and loaded states.     

Subject-specific finite element models can accurately predict strain levels in long bones

The prediction of the stress-state and fracture risk induced in bones by various loading conditions in individual patients using subject-specific finite element models still represents a challenge in orthopaedic biomechanics. The accuracy of the strain predictions reported in the literature is variable and generally not satisfactory. The aim of the present study was to evaluate if a proper choice of the density-elasticity relationship can lead to accurate strain predictions in the frame of an automatic subject-specific model generation strategy. To this aim, a combined numerical-experimental study was performed comparing finite element predicted strains with strain-gauges measurements obtained on eight cadaver proximal femurs, each instrumented with 15 rosettes mostly concentrated in the bone metaphyses, tested non-destructively in vitro under six different loading scenarios. Three different density-elasticity power relationships were selected from the literature and implemented in the finite element models derived from computed tomography data. The results of the present study confirm the great influence of the density-elasticity relationship used on the accuracy of numerical predictions. One of the tested constitutive laws provided a very good agreement (R(2)=0.91, RMSE lower than 10% of the maximum measured value) between numerical calculations and experimental measurements. The presented results show, in addition, that the adoption of a single density-elasticity relationship over the whole bone density range is adequate to obtain an accuracy that is already suitable for many applications.     

Comparison of strain measurement in the mouse forearm using subject-specific finite element models,strain gaging,and digital image correlation

Mechanical loading in bone leads to the activation of bone-forming pathways that are most likely associated with a minimum strain threshold being experienced by the osteocyte. To investigate the correlation between cellular response and mechanical stimuli, researchers must develop accurate ways to measure/compute strain both externally on the bone surface and internally at the osteocyte level. This study investigates the use of finite element (FE) models to compute bone surface strains on the mouse forearm. Strains from three FE models were compared to data collected experimentally through strain gaging and digital image correlation (DIC). Each FE model was assigned subject-specific bone properties and consisted of one-dimensional springs representing the interosseous membrane. After three-point bending was performed on the ulnae and radii, moment of inertia was determined from microCT analysis of the bone region between the supports and then used along with standard beam analyses to calculate the Young’s modulus. Non-contact strain measurements from DIC were determined to be more suitable for validating numerical results than experimental data obtained through conventional strain gaging. When comparing strain responses in the three ulnae, we observed a 3–14% difference between numerical and DIC strains while the strain gage values were 37–56% lower than numerical values. This study demonstrates a computational approach for capturing bone surface strains in the mouse forearm. Ultimately, strains from these macroscale models can be used as inputs for microscale and nanoscale FE models designed to analyze strains directly in the osteocyte lacunae.     

Patient-specific finite element analysis of the human femur--a double-blinded biomechanical validation

Patient-specific finite element (PSFE) models based on quantitative computer tomography (qCT) are generally used to "predict" the biomechanical response of human bones with the future goal to be applied in clinical decision-making. However, clinical applications require a well validated tool that is free of numerical errors and furthermore match closely experimental findings. In previous studies, not all measurable data (strains and displacements) were considered for validation. Furthermore, the same research group performed both the experiments and PSFE analyses; thus, the validation may have been biased. The aim of the present study was therefore to validate PSFE models with biomechanical experiments, and to address the above-mentioned issues of measurable data and validation bias. A PSFE model (p-method) of each cadaver femur (n = 12) was generated based on qCT scans of the specimens. The models were validated by biomechanical in-vitro experiments, which determined strains and local displacements on the bone surface and the axial stiffness of the specimens. The validation was performed in a double-blinded manner by two different research institutes to avoid any bias. Inspecting all measurements (155 values), the numerical results correlated well with the experimental results (R(2) = 0.93, slope 1.0093, mean of absolute deviations 22%). In conclusion, a method to generate PSFE models from qCT scans was used in this study on a sample size not yet considered in the past, and compared to experiments in a douple-blinded manner. The results demonstrate that the presented method is in an advanced stage, and can be used in clinical computer-aided decision-making.