Experimental validation of finite element models of intact and implanted composite hemipelvises using digital image correlation

A detailed understanding of the changes in load transfer due to implantation is necessary to identify potential failure mechanisms of orthopedic implants. Computational finite element (FE) models provide full field data on intact and implanted bone structures, but their validity must be assessed for clinical relevance. The aim of this study was to test the validity of FE predicted strain distributions for the intact and implanted pelvis using the digital image correlation (DIC) strain measurement technique. FE models of an in vitro hemipelvis test setup were produced, both intact and implanted with an acetabular cup. Strain predictions were compared to DIC and strain rosette measurements. Regression analysis indicated a strong linear relationship between the measured and predicted strains, with a high correlation coefficient (R?=?0.956 intact, 0.938 implanted) and a low standard error of the estimate (SE?=?69.53?με, 75.09?με). Moreover, close agreement between the strain rosette and DIC measurements improved confidence in the validity of the DIC technique. The FE model therefore was supported as a valid predictor of the measured strain distribution in the intact and implanted composite pelvis models, confirming its suitability for further computational investigations.     

Determining constitutive behavior of the brain tissue using digital image correlation and finite element modeling

Biomechanics and Modeling in Mechanobiology - Detailed knowledge about the mechanical properties of brain can improve numerical modeling of the brain under various loading conditions. The success...     

Development and validation of subject-specific finite element models for blunt trauma study

This study developed and validated finite element (FE) models of swine and human thoraxes and abdomens that had subject-specific anatomies and could accurately and efficiently predict body responses to blunt impacts. Anatomies of the rib cage, torso walls, thoracic, and abdominal organs were reconstructed from X-ray computed tomography (CT) images and extracted into geometries to build FE meshes. The rib cage was modeled as an inhomogeneous beam structure with geometry and bone material parameters determined directly from CT images. Meshes of soft components were generated by mapping structured mesh templates representative of organ topologies onto the geometries. The swine models were developed from and validated by 30 animal tests in which blunt insults were applied to swine subjects and CT images, chest wall motions, lung pressures, and pathological data were acquired. A comparison of the FE calculations of animal responses and experimental measurements showed a good agreement. The errors in calculated response time traces were within 10% for most tests. Calculated peak responses showed strong correlations with the experimental values. The stress concentration inside the ribs, lungs, and livers produced by FE simulations also compared favorably to the injury locations. A human FE model was developed from CT images from the Visible Human project and was scaled to simulate historical frontal and side post mortem human subject (PMHS) impact tests. The calculated chest deformation also showed a good agreement with the measurements. The models developed in this study can be of great value for studying blunt thoracic and abdominal trauma and for designing injury prevention techniques, equipments, and devices.     

Prediction of the structural response of the femoral shaft under dynamic loading using subject-specific finite element models

The goal of this study was to predict the structural response of the femoral shaft under dynamic loading conditions using subject-specific finite element (SS-FE) models and to evaluate the prediction accuracy of the models in relation to the model complexity. In total, SS-FE models of 31 femur specimens were developed. Using those models, dynamic three-point bending and combined loading tests (bending with four different levels of axial compression) of bare femurs were simulated, and the prediction capabilities of five different levels of model complexity were evaluated based on the impact force time histories: baseline, mass-based scaled, structure-based scaled, geometric SS-FE, and heterogenized SS-FE models. Among the five levels of model complexity, the geometric SS-FE and the heterogenized SS-FE models showed statistically significant improvement on response prediction capability compared to the other model formulations whereas the difference between two SS-FE models was negligible. This result indicated the geometric SS-FE models, containing detailed geometric information from CT images with homogeneous linear isotropic elastic material properties, would be an optimal model complexity for prediction of structural response of the femoral shafts under the dynamic loading conditions. The average and the standard deviation of the RMS errors of the geometric SS-FE models for all the 31 cases was 0.46 kN and 0.66 kN, respectively. This study highlights the contribution of geometric variability on the structural response variation of the femoral shafts subjected to dynamic loading condition and the potential of geometric SS-FE models to capture the structural response variation of the femoral shafts.     

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.     

An accurate estimation of bone density improves the accuracy of subject-specific finite element models

An experimental-numerical study was performed to investigate the relationships between computed tomography (CT)-density and ash density, and between ash density and apparent density for bone tissue, to evaluate their influence on the accuracy of subject-specific FE models of human bones. Sixty cylindrical bone specimens were examined. CT-densities were computed from CT images while apparent and ash densities were measured experimentally. The CT/ash-density and ash/apparent-density relationships were calculated. Finite element models of eight human femurs were generated considering these relationships to assess their effect on strain prediction accuracy. CT and ash density were linearly correlated (R(2)=0.997) over the whole density range but not equivalent (intercep t <0, slope >1). A constant ash/apparent-density ratio (0.598+/-0.004) was found for cortical bone. A lower ratio, with a larger dispersion, was found for trabecular bone (0.459+/-0.100), but it became less dispersed, and equal to that of cortical tissue, when testing smaller trabecular specimens (0.598+/-0.036). This suggests that an experimental error occurred in apparent-density measurements for large trabecular specimens and a constant ratio can be assumed valid for the whole density range. Introducing the obtained relationships in the FE modelling procedure improved strain prediction accuracy (R(2)=0.95, RMSE=7%). The results suggest that: (i) a correction of the densitometric calibration should be used when evaluating bone ash-density from clinical CT scans, to avoid ash-density underestimation and overestimation for low- and high-density bone tissue, respectively; (ii) the ash/apparent-density ratio can be assumed constant in human femurs and (iii) the correction improves significantly the model accuracy and should be considered in subject-specific bone modelling.     

Sensitivity of periprosthetic stress-shielding to load and the bone density-modulus relationship in subject-specific finite element models

Subject-specific finite element (FE) computer models of the proximal femur in hip replacement could potentially predict stress-shielding and subsequent bone loss in individual patients. Before such predictions can be made, it is important first to determine if between subject differences in stress-shielding are sensitive to poorly defined parameters such as the load and the bone material properties. In this study we investigate if subject-specific FE models provide consistent stress-shielding patterns in the bone, independent of the choice of the loading conditions and the bone density-modulus relationship used in the computer model. FE models of two right canine femurs with and without implants were constructed based on contiguous computed tomography (CT) scans so that subject-specific estimates of stress-shielding could be calculated. Four different loading conditions and two bone density-modulus relationships were tested. Stress-shielding was defined as the decrease of strain energy per gram bone mass in the femur with the implant in place relative to the intact femur.The analyses showed that for the four loading conditions and two bone density-modulus relationships the difference in stress-shielding between the two subjects was essentially constant (1% variation) when the same loading condition and density-modulus relationship was used for both subjects. The severity of stress-shielding within a subject was sensitive to these input parameters, varying up to 20% in specific regions with a change in loading conditions and up to 10% for a change in the assumed density-modulus relationship. We conclude that although the choice of input parameters can substantially affect stress-shielding in an individual, this choice had virtually no effect on the relative differences in femoral periprosthetic stress-shielding between individuals. Thus, while care should be taken in the interpretation of the absolute value of stress-shielding calculated with these type of models, subject-specific FE models may be useful for explaining the variation in bone adaptation responsiveness between different subjects in experimental or clinical studies.     

Toward subject-specific evaluation: methods of evaluating finite element brain models using experimental high-rate rotational brain motion

Biomechanics and Modeling in Mechanobiology - Computational models of the brain have become the gold standard in biomechanics to understand, predict, and mitigate traumatic brain injuries. Many...     

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.     

Effect of round curvature of anterior implant-supported zirconia frameworks: finite element analysis and in vitro study using digital image correlation

Two groups of 4-unit zirconia frameworks were produced by CAD/CAM to simulate the restoration of an anterior edentulous gap supported by 2 implant-abutment assemblies. Group 1 comprised straight configuration frameworks and group 2 consisted of arched frameworks. Specimens were made with the same connector cross-section area and were cemented and submitted to static loads. Displacements were captured with two high-speed photographic cameras and analysed with video correlation system. Frameworks and the implant-abutment assembly were scanned and converted to 3DCAD objects by reverse engineering process. A specimen of each group was veneered and the corresponding 3D geometry was similarly obtained after scanning. Numerical models were created from the CAD objects and the FE analysis was performed on the zirconia frameworks and on the FPDs bi-layered with porcelain (veneered frameworks). Displacements were higher for the curved frameworks group, under any load. The predicted displacements correlated well with the experimental values of the two framework groups, but on the straight framework the experimental vertical displacements were superior to those predicted by the FEA. The results showed that the round curvature of zirconia anterior implant-supported FPDs plays a significant role on the deformation/stress of FPDs that cannot be neglected neither in testing nor in simulation and should be considered in the clinical setting.     

Automatic generation of accurate subject-specific bone finite element models to be used in clinical studies

Most of the finite element models of bones used in orthopaedic biomechanics research are based on generic anatomies. However, in many cases it would be useful to generate from CT data a separate finite element model for each subject of a study group. In a recent study a hexahedral mesh generator based on a grid projection algorithm was found very effective in terms of accuracy and automation. However, so far the use of this method has been documented only on data collected in vitro and only for long bones. The present study was aimed at verifying if this method represents a procedure for the generation of finite element models of human bones from data collected in vivo, robust, accurate, automatic and general enough to be used in clinical studies. Robustness, automation and numerical accuracy of the proposed method were assessed on five femoral CT data sets of patients affected by various pathologies. The generality of the method was verified by processing a femur, an ileum, a phalanx, a proximal femur reconstruction, and the micro-CT of a small sample of spongy bone. The method was found robust enough to cope with the variability of the five femurs, producing meshes with a numerical accuracy and a computational weight comparable to those found in vitro. Even when the method was used to process the other bones the levels of mesh conditioning remained within acceptable limits. Thus, it may be concluded that the method presents a generality sufficient to cope with almost any orthopaedic application.     

High resolution three-dimensional strain mapping of bioprosthetic heart valves using digital image correlation

Transcatheter aortic valve replacement (TAVR) is a safe and effective treatment option for patients deemed at high and intermediate risk for surgical aortic valve replacement. Similar to surgical aortic valves (SAVs), transcatheter aortic valves (TAVs) undergo calcification and mechanical wear over time. However, to date, there have been limited publications on the long-term durability of TAV devices. To assess longevity and mechanical strength of TAVs in comparison to surgical bioprosthetic valves, three-dimensional deformation analysis and strain measurement of the leaflets become an inevitable part of the evaluation. The goal of this study was to measure and compare leaflet displacement and strain of two commonly used TAVs in a side-by-side comparison with a commonly used SAV using a high-resolution digital image correlation (DIC) system. 26-mm Edwards SAPIEN 3, 26-mm Medtronic CoreValve, and 25-mm Carpentier-Edwards PERIMOUNT Magna surgical bioprosthesis were examined in a custom-made valve testing apparatus. A time-varying, spatially uniform pressure was applied to the leaflets at different loading rates. GOM ARAMIS® software was used to map leaflet displacement and strain fields during loading and unloading. High displacement regions were found to be at the leaflet belly region of the three bioprosthetic valves. In addition, the frame of the surgical bioprosthesis was found to be remarkably flexible, in contrary to CoreValve and SAPIEN 3 in which the stent was nearly rigid under a similar loading condition. The experimental DIC measurements can be used to characterize the anisotropic materiel behavior of the bioprosthetic heart valve leaflets and validate heart valve computational simulations.     

The accuracy of digital image-based finite element models.

Digital image-based finite element meshing is an alternative approach to time-consuming conventional meshing techniques for generating realistic three-dimensional (3D) models of complex structures. Although not limited to biological applications, digital image-based modeling has been used to generate structure-specific (i.e., non-generic) models of whole bones and trabecular bone microstructures. However, questions remain regarding the solution accuracy provided by the digital meshing approach, particularly at model or material boundaries. The purpose of this study was to compare the accuracy of digital and conventional smooth boundary models based on theoretical solutions for a two-dimensional (2D) compression plate and a 3D circular cantilever beam. For both the plate and beam analyses, the predicted solution at digital model boundaries was characterized by local oscillations, which produced potentially high errors within individual boundary elements. Significantly, however, the digital model boundary solution oscillated approximately about the theoretical solution. A marked improvement in solution accuracy was therefore achieved by considering average results within a region composed of several elements. Absolute errors for Von Mises stress averaged over the beam cross section, for example, converged to less than 4 percent, and the predicted free-end displacement of the cantilever beam was within 1 percent of the theoretical solution. Analyses at several beam orientations and mesh resolutions suggested a minimum discretization of three to four digital finite elements through the beam cross section to avoid high numerical stiffening errors under bending.     

Comparison between mechanical stress and bone mineral density in the femur after total hip arthroplasty by using subject-specific finite element analyses

The mechanism underling bone mineral density (BMD) loss that occurs in the femur after total hip arthroplasty (THA) remains unknown. We compared the equivalent stress and strain energy density (SED) to BMD in the femur after THA using subject-specific finite element analyses. Twenty-four patients who had undergone primary cementless THA were analysed. BMD was measured using dual-energy X-ray absorptiometry (DEXA) at 1 week and 3, 6 and 12 months after THA. Seven regions of interest (ROIs) were defined in accordance with Gruen's system (ROIs 1–7). Computed tomography images of the femurs were acquired pre- and postoperatively, and the images were converted into three-dimensional finite element (FE) models. Equivalent stress and SED were analysed and compared with DEXA data. BMD was maintained 1 year after THA in ROIs 3, 4, 5 and 6, whereas BMD decreased in ROIs 1, 2 and 7. FE analysis revealed that equivalent stress in ROIs 3, 4, 5 and 6 was much higher than that in ROIs 1, 2 and 7. A significant correlation was observed between the rate of changes in BMD and equivalent stress. Reduction of equivalent stress may contribute to decrease in BMD in the femur after THA.     

A three-dimensional digital image correlation technique for strain measurements in microstructures

A three-dimensional digital image correlation technique is presented for strain measurements in open-cell structures such as trabecular bone. The technique uses high-resolution computed tomography images for displacement measurements in the solid structure. In order to determine the local strain-state within single trabeculae, a tetrahedronization method is used to fill the solid structure with tetrahedrae. Displacements are calculated at the nodes of the tetrahedrae. The displacement data is subsequently converted to a deformation tensor in each of the tetrahedral element centers with a least-squares estimation method. Because the trabeculae are represented by a mesh, it is possible to deform this mesh according to the deformation tensor and, at the same time, visualize the calculated local strain in the deformed mesh with a finite element post-processing tool. In this way, the deformation of a single trabecula from an aluminum foam sample was determined and validated with rendered images of the three-dimensional sample. A precision analysis showed that a rigid translation or rotation does not affect the accuracy. Typical values for the standard deviation in the displacement and strain components are 2.0 microm and 0.01, respectively. Presently, the precision limits the technique to strain measurements beyond the yield strain.     

A subject-specific anisotropic visco-hyperelastic finite element model of female pelvic floor stress and strain during the second stage of labor

ObjectivesTo develop an improved model representation of the biomechanics of the levator muscles during the second stage of labor and to use a sensitivity analysis to explore the pathomechanics of levator muscle injury.MethodsA subject-specific finite element model of human pelvic floor and fetal head was developed based on in vivo MRI data of a fetal head and maternal pelvis. An anisotropic visco-hyperelastic constitutive model employed material parameters estimated from biaxial tests on pelvic floor tissues. Boundary conditions reflected both anatomic constraints and the curve of Carus. A short second stage of labor, scaled to 10 min, was then simulated using a single expulsive push made in the absence of levator co-contraction.ResultsLarge levator stresses occurred near the levator hiatus reaching 9 MPa at the pubovisceral muscle enthesis. The dominant principal stresses were located at, and aligned with, the edge of the hiatus. Muscle stretch bordering the levator hiatus was inhomogeneous: the average levator stretch was 3.55 with a high of 4.64 at the pubovisceral muscle enthesis. Decreasing perineal body stiffness by 40%, 50%, and 60% led to reductions in the maximum principal stretch ratio at the pubovisceral muscle enthesis of 8%, 13%, and 18%, respectively.ConclusionsThe pubovisceral muscle enthesis and the muscle near the perineal body are the regions of greatest strain thereby placing them at highest risk for stretch-related injury. Decreasing perineal body tissue stiffness significantly reduced tissue stress and strain, and therefore injury risk, in those regions.     

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.     

Analysis of the compressive strain below the removable and fixed prosthesis in the posterior mandible using a digital image correlation method

It was the purpose of this study to determine and analyse strains in the bone of posterior mandible below the removable and fixed partial dentures using Digital Image Correlation Method. Dried mandible with shortened dental arch (Kennedy Class 1) was used in the experiment. The mandible model was prepared for accepting two types of restorations for bilaterally missing molars conventional therapy, and removable and fixed partial dentures were manufactured following standard prosthetic protocol. The models with prosthetic restorations placed in situ were subjected to loading of 50-300 N, and photographed using two digital cameras as part of the digital image correlation method equipment. Afterwards, the obtained data for strains within restored mandible during loading ciclus were analysed in the software Aramis and graphically presented. Percentage size of the mandible strains within the line section are from 0.14 to 0.19% for the removable partial denture experimental model and 0-0.14% for the fixed partial denture experimental model. The study has demonstrated that Digital Image Correlation method can be used to measure strain on the mandible surface and replacements during loading and that from biomechanical perspective both therapeutic modalities can be equally useful for the restoration of the mandible with bilaterally distal edentulous spaces.