1
|
Gardegaront M, Sas A, Brizard D, Levillain A, Bermond F, Confavreux CB, Pialat JB, van Lenthe GH, Follet H, Mitton D. Inter-laboratory reproduction and sensitivity study of a finite element model to quantify human femur failure load: Case of metastases. J Mech Behav Biomed Mater 2024; 158:106676. [PMID: 39121530 DOI: 10.1016/j.jmbbm.2024.106676] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2023] [Revised: 04/19/2024] [Accepted: 07/24/2024] [Indexed: 08/12/2024]
Abstract
INTRODUCTION Metastases increase the risk of fracture when affecting the femur. Consequently, clinicians need to know if the patient's femur can withstand the stress of daily activities. The current tools used in clinics are not sufficiently precise. A new method, the CT-scan-based finite element analysis, gives good predictive results. However, none of the existing models were tested for reproducibility. This is a critical issue to address in order to apply the technique on a large cohort around the world to help evaluate bone metastatic fracture risk in patients. The aim of this study is then to evaluate 1) the reproducibility 2) the transposition of the reproduced model to another dataset and 3) the global sensitivity of one of the most promising models of the literature (original model). METHODS The model was reproduced based on the paper describing it and discussion with authors to avoid reproduction errors. The reproducibility was evaluated by comparing the results given in the original model by the original first team (Leuven, Belgium) and the reproduced model made by another team (Lyon, France) on the same dataset of CT-scans of ex vivo femurs. The transposition of the model was evaluated by comparing the results of the reproduced model on two different datasets. The global sensitivity analysis was done by using the Morris method and evaluates the influence of the density calibration coefficient, the segmentation, the orientations and the length of the femur. RESULTS The original and reproduced models are highly correlated (r2 = 0.95), even though the reproduced model gives systematically higher failure loads. When using the reproduced model on another dataset, predictions are less accurate (r2 with the experimental failure load decreases, errors increase). The global sensitivity analysis showed high influence of the density calibration coefficient (mean variation of failure load of 84 %) and non-negligible influence of the segmentation, orientation and length of the femur (mean variation of failure load between 7 and 10 %). CONCLUSION This study showed that, although being validated, the reproduced model underperformed when using another dataset. The difference in performance depending on the dataset is commonly the cause of overfitting when creating the model. However, the dataset used in the original paper (Sas et al., 2020a) and the Leuven's dataset gave similar performance, which indicates a lesser probability for the overfitting cause. Also, the model is highly sensitive to density parameters and automation of measurement may minimize the uncertainty on failure load. An uncertainty propagation analysis would give the actual precision of such model and improve our understanding of its behavior and is part of future work.
Collapse
Affiliation(s)
- Marc Gardegaront
- Univ Lyon, Univ Claude Bernard Lyon 1, INSERM, LYOS UMR 1033, 69008, Lyon, France; Univ Lyon, Univ Eiffel, Univ Claude Bernard Lyon 1, LBMC UMR_T9406, 69622, Lyon, France.
| | - Amelie Sas
- Biomechanics Section, Dept. Mechanical Engineering, KU Leuven, Leuven, Belgium.
| | - Denis Brizard
- Univ Lyon, Univ Eiffel, Univ Claude Bernard Lyon 1, LBMC UMR_T9406, 69622, Lyon, France.
| | - Aurélie Levillain
- Univ Lyon, Univ Eiffel, Univ Claude Bernard Lyon 1, LBMC UMR_T9406, 69622, Lyon, France.
| | - François Bermond
- Univ Lyon, Univ Eiffel, Univ Claude Bernard Lyon 1, LBMC UMR_T9406, 69622, Lyon, France.
| | - Cyrille B Confavreux
- Univ Lyon, Univ Claude Bernard Lyon 1, INSERM, LYOS UMR 1033, 69008, Lyon, France; Centre Expert des Métastases Osseuses (CEMOS), Hôpital Lyon Sud, Hospices Civils de Lyon, France.
| | - Jean-Baptiste Pialat
- Centre Expert des Métastases Osseuses (CEMOS), Hôpital Lyon Sud, Hospices Civils de Lyon, France; Creatis CNRS UMR 5220, INSERM, U1294, Université Lyon 1, Villeurbanne, France.
| | - G Harry van Lenthe
- Biomechanics Section, Dept. Mechanical Engineering, KU Leuven, Leuven, Belgium.
| | - Hélène Follet
- Univ Lyon, Univ Claude Bernard Lyon 1, INSERM, LYOS UMR 1033, 69008, Lyon, France.
| | - David Mitton
- Univ Lyon, Univ Eiffel, Univ Claude Bernard Lyon 1, LBMC UMR_T9406, 69622, Lyon, France.
| |
Collapse
|
2
|
Dudle A, Ith M, Egli R, Heverhagen J, Gugler Y, Wapp C, Frauchiger DA, Lippuner K, Jackowski C, Zysset P. Asynchronous calibration of a CT scanner for bone mineral density estimation: sources of error and correction. JBMR Plus 2024; 8:ziae096. [PMID: 39183821 PMCID: PMC11344033 DOI: 10.1093/jbmrpl/ziae096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/13/2024] [Revised: 06/20/2024] [Accepted: 07/12/2024] [Indexed: 08/27/2024] Open
Abstract
The estimation of BMD with CT scans requires a calibration method, usually based on a phantom. In asynchronous calibration, the phantom is scanned separately from the patient. A standardized acquisition protocol must be used to avoid variations between patient and phantom. However, variations can still be induced, for example, by temporal fluctuations or patient characteristics. Based on the further use of 739 forensic and 111 clinical CT scans, this study uses the proximal femur BMD value ("total hip") to assess asynchronous calibration accuracy, using in-scan calibration as ground truth. It identifies the parameters affecting the calibration accuracy and quantifies their impact. For time interval and table height, the impact was measured by calibrating the CT scan twice (once using the phantom scan with closest acquisition parameters and once using a phantom scan with standard values) and comparing the calibration accuracy. For other parameters such as body weight, the impact was measured by computing a linear regression between parameter values and calibration accuracy. Finally, this study proposes correction methods to reduce the effect of these parameters and improve the calibration accuracy. The BMD error of the asynchronous calibration, using the phantom scan with the closest acquisition parameters, was -1.2 ± 1.7% for the forensic and - 1.6 ± 3.5% for the clinical dataset. Among the parameters studied, time interval and body weight were identified as the main sources of error for asynchronous calibration, followed by table height and reconstruction kernel. Based on these results, a correction method was proposed to improve the calibration accuracy.
Collapse
Affiliation(s)
- Alice Dudle
- ARTORG Center for Biomedical Engineering Research, University of Bern, sitem-insel, Freiburgstrasse 3, Bern 3010, Switzerland
| | - Michael Ith
- Department of Diagnostic, Interventional and Pediatric Radiology, Inselspital, University of Bern, Freiburgstrasse 10, Bern 3010, Switzerland
| | - Rainer Egli
- Department of Diagnostic, Interventional and Pediatric Radiology, Inselspital, University of Bern, Freiburgstrasse 10, Bern 3010, Switzerland
| | - Johannes Heverhagen
- Department of Diagnostic, Interventional and Pediatric Radiology, Inselspital, University of Bern, Freiburgstrasse 10, Bern 3010, Switzerland
| | - Yvan Gugler
- ARTORG Center for Biomedical Engineering Research, University of Bern, sitem-insel, Freiburgstrasse 3, Bern 3010, Switzerland
| | - Christina Wapp
- ARTORG Center for Biomedical Engineering Research, University of Bern, sitem-insel, Freiburgstrasse 3, Bern 3010, Switzerland
| | - Daniela A Frauchiger
- ARTORG Center for Biomedical Engineering Research, University of Bern, sitem-insel, Freiburgstrasse 3, Bern 3010, Switzerland
- Department of Osteoporosis, Inselspital, Bern University Hospital, University of Bern, Freiburgstrsasse 4, Bern 3010, Switzerland
| | - Kurt Lippuner
- Department of Osteoporosis, Inselspital, Bern University Hospital, University of Bern, Freiburgstrsasse 4, Bern 3010, Switzerland
| | - Christian Jackowski
- Institute of Forensic Medicine, University of Bern, Murtenstrasse 28, Bern 3008, Switzerland
| | - Philippe Zysset
- ARTORG Center for Biomedical Engineering Research, University of Bern, sitem-insel, Freiburgstrasse 3, Bern 3010, Switzerland
| |
Collapse
|
3
|
Amini M, Reisinger A, Synek A, Hirtler L, Pahr D. The predictive ability of a QCT-FE model of the proximal femoral stiffness under multiple load cases is strongly influenced by experimental uncertainties. J Mech Behav Biomed Mater 2023; 139:105664. [PMID: 36657193 DOI: 10.1016/j.jmbbm.2023.105664] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2022] [Revised: 01/03/2023] [Accepted: 01/05/2023] [Indexed: 01/11/2023]
Abstract
Despite significant improvements in terms of the predictive ability of Quantitative Computed Tomography based Finite Element (QCT-FE) models in estimating femoral strength (fracture load and stiffness), no substantial clinical adoption of this method has taken place to date. Narrowing the wide variability of FE results by standardizing the methodology and validation protocols, as well as reducing the uncertainties in the FEA process have been proposed as routes towards improved reliability. The aim of this study was to: First, validate a QCT-FE model of proximal femoral stiffness in multiple stance load cases, and second, using a parametric approach, determine the influence of select experimental and modeling parameters on the predictive ability of our model. Ten fresh frozen human femoral samples were tested in neutral stance, 15° adducted and 15° abducted load cases. Voxel-based linear-elastic QCT-FE models of the samples were generated to predict the models' stiffness values in all load cases. The base FE models were validated against the experimental results using linear regression. Thirty six deviated models were created using the minimum and maximum values of experiment-based "plausible range" for 18 parameters in 4 categories of embedding, loading, material, and segmentation. The predictive ability of the models were compared in terms of the coefficient of determination (R2) of the linear regression between the measured and predicted stiffness values in all load cases. Our model was capable of capturing 90% of the variation in the experimental stiffness of the samples in neutral stance position (R2 = 0.9, concordance correlation coefficient (CCC) = 0.93, percent root mean squared error (RMSE%) = 8.4%, slope and intercept not significantly different from unity and zero, respectively). Embedding and loading categories strongly affected the predictive ability of the models with an average percent difference in R2 of 4.36% ± 2.77 and 2.96% ± 1.69 for the stance-neutral load case, respectively. The performance of the models were significantly different in adducted and abducted load cases with their R2 dropping to 71% and 70%, respectively. Similarly, off-axes load cases were affected by the parameters differently compared to the neutral load case, with the loading parameter category imposing more than 10% difference on their R2, larger than all other categories. We also showed that automatically selecting the best performing plausible value for each parameter and each sample would result in a perfectly linear correlation (R2> 0.99) between the "tuned" model's predicted stiffness and experimental results. Based on our results, high sensitivity of the model performance to experimental parameters requires extra diligence in modeling the embedding geometry and the loading angles since these sources of uncertainty could dwarf the effects of material modeling and image processing parameters. The results of this study could help in improving the robustness of the QCT-FE models of proximal femur by limiting the uncertainties in the experimental and modeling steps.
Collapse
Affiliation(s)
- Morteza Amini
- Institute of Lightweight Design and Structural Biomechanics, TU Wien, Getreidemarkt 9, 1060, Vienna, Austria.
| | - Andreas Reisinger
- Division Biomechanics, Karl Landsteiner University of Health Sciences, Dr.-Karl-Dorrek-Straße 30, 3500 Krems an der Donau, Austria.
| | - Alexander Synek
- Institute of Lightweight Design and Structural Biomechanics, TU Wien, Getreidemarkt 9, 1060, Vienna, Austria.
| | - Lena Hirtler
- Center for Anatomy and Cell Biology, Medical University of Vienna, Währinger Straße 13, 1090, Vienna, Austria.
| | - Dieter Pahr
- Institute of Lightweight Design and Structural Biomechanics, TU Wien, Getreidemarkt 9, 1060, Vienna, Austria; Division Biomechanics, Karl Landsteiner University of Health Sciences, Dr.-Karl-Dorrek-Straße 30, 3500 Krems an der Donau, Austria.
| |
Collapse
|
4
|
Eliyahu L, Yosibash Z, Avivi I, Cohen YC, Ariel G, Sadovnic O, Sternheim A. On the influence of computed tomography's slice thickness on computer tomography based finite element analyses results. Clin Biomech (Bristol, Avon) 2023; 102:105889. [PMID: 36774735 DOI: 10.1016/j.clinbiomech.2023.105889] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/03/2022] [Revised: 12/06/2022] [Accepted: 01/09/2023] [Indexed: 01/19/2023]
Abstract
BACKGROUND Patient-specific autonomous finite element analyses of femurs, based on clinical computed tomography scans may be used to monitor the progression of bone-related diseases. Some CT scan protocols provide lower resolution (slice thickness of 3 mm) that affects the accuracy. To investigate the impact of low-resolution scans on the CT-based finite element analyses results, identical CT raw data were reconstructed twice to generate a 1 mm ("gold standard") and a 3 mm slice thickness scans. METHODS CT-based finite element analyses of twenty-four femurs (twelve patients) under stance and sideways fall loads were performed based on 1 and 3 mm slice thickness scans. Bone volume, load direction, and strains were extracted at different locations along the femurs and differences were evaluated. FINDINGS Average differences in bone volume were 1.0 ± 1.5%. The largest average difference in strains in stance position was in the neck region (11.0 ± 13.4%), whereas in other regions these were much smaller. For sidewise fall loading, the average differences were at most 9.2 ± 16.0%. INTERPRETATION Whole-body low dose CT scans (3 mm-slice thickness) are suboptimal for monitoring strain changes in patient's femurs but may allow longitudinal studies if larger than 5% in all areas and larger than 12% in the upper neck. CT-based finite element analyses with slice thickness of 3 mm may be used in clinical practice for patients with smoldering myeloma to associate changes in strains with progression to active myeloma if above ∼10%.
Collapse
Affiliation(s)
- Leetal Eliyahu
- Computational Mechanics and Experimental Biomechanics Lab, School of Mechanical Engineering, Tel-Aviv University, Israel
| | - Zohar Yosibash
- Computational Mechanics and Experimental Biomechanics Lab, School of Mechanical Engineering, Tel-Aviv University, Israel.
| | - Irit Avivi
- Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel; Hematology Division, Tel Aviv Sourasky Medical Center, Tel Aviv, Israel
| | - Yael C Cohen
- Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel; Hematology Division, Tel Aviv Sourasky Medical Center, Tel Aviv, Israel
| | - Gal Ariel
- Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel; National Unit of Orthopaedic Oncology, Tel Aviv Sourasky Medical Center, Tel Aviv, Israel
| | - Ofer Sadovnic
- Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel; Radiology Division, Tel Aviv Sourasky Medical Center, Tel Aviv, Israel
| | - Amir Sternheim
- Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel; National Unit of Orthopaedic Oncology, Tel Aviv Sourasky Medical Center, Tel Aviv, Israel
| |
Collapse
|
5
|
Finite Element Method for the Evaluation of the Human Spine: A Literature Overview. J Funct Biomater 2021; 12:jfb12030043. [PMID: 34449646 PMCID: PMC8395922 DOI: 10.3390/jfb12030043] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2021] [Revised: 07/23/2021] [Accepted: 07/29/2021] [Indexed: 02/07/2023] Open
Abstract
The finite element method (FEM) represents a computer simulation method, originally used in civil engineering, which dates back to the early 1940s. Applications of FEM have also been used in numerous medical areas and in orthopedic surgery. Computing technology has improved over the years and as a result, more complex problems, such as those involving the spine, can be analyzed. The spine is a complex anatomical structure that maintains the erect posture and supports considerable loads. Applications of FEM in the spine have contributed to the understanding of bone biomechanics, both in healthy and abnormal conditions, such as scoliosis, fractures (trauma), degenerative disc disease and osteoporosis. However, since FEM is only a digital simulation of the real condition, it will never exactly simulate in vivo results. In particular, when it concerns biomechanics, there are many features that are difficult to represent in a FEM. More FEM studies and spine research are required in order to examine interpersonal spine stiffness, young spine biomechanics and model accuracy. In the future, patient-specific models will be used for better patient evaluations as well as for better pre- and inter-operative planning.
Collapse
|
6
|
Prado M, Khosla S, Chaput C, Giambini H. Opportunistic application of phantom-less calibration methods for fracture risk prediction using QCT/FEA. Eur Radiol 2021; 31:9428-9435. [PMID: 34047849 DOI: 10.1007/s00330-021-08071-w] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2021] [Revised: 04/30/2021] [Accepted: 05/11/2021] [Indexed: 11/25/2022]
Abstract
OBJECTIVES Quantitative computed tomography (QCT)-based finite element analysis (FEA) implements a calibration phantom to estimate bone mineral density (BMD) and assign material properties to the models. The objectives of this study were to (1) propose robust phantom-less calibration methods, using subject-specific tissues, to obtain vertebral fracture properties estimations using QCT/FEA; and (2) correlate QCT/FEA predictions to DXA values of areal BMD. METHODS Eighty of a cohort of 111 clinical QCT scans were used to obtain subject-specific parameters using a phantom calibration approach and for the development of the phantom-less calibration equations. Equations were developed based on the HU measured from various soft tissues and regions, and using multiple linear regression analyses. Thirty-one additional QCT scans were used for cross-validation of QCT/FEA estimated fracture loads from the L3 vertebrae based on the phantom and phantom-less equations. Finally, QCT/FEA-predicted fracture loads were correlated with aBMD obtained from DXA. RESULTS Overall, 217 QCT/FEA models from 31 subjects (20 females, 11 men) with mean ages of 69.6 (13.1) and 67.3 (14) were used to cross-validate the phantom-less equations and assess bone strength. The proposed phantom-less equations showed high correlations with phantom-based estimates of BMD (99%). Cross-validation of QCT/FEA-predicted fracture loads from phantom-less equations and phantom-specific outcomes resulted in high correlations for all proposed methods (0.94-0.99). QCT/FEA correlation outcomes from the phantom-less equations and DXA-aBMD were moderately high (0.64-0.68). CONCLUSIONS The proposed QCT/FEA subject-specific phantom-less calibration methods demonstrated the potential to be applied to both prospective and retrospective applications in the clinical setting. KEY POINTS • QCT/FEA overcomes the disadvantages of DXA and improves fracture properties predictions of vertebrae. • QCT/FEA fracture estimates using the phantom-less approach highly correlated to values obtained using a calibration phantom. • QCT/FEA prediction using a phantom-less approach is an accurate alternative over phantom-based methods.
Collapse
Affiliation(s)
- Maria Prado
- Department of Biomedical Engineering and Chemical Engineering, University of Texas at San Antonio, One UTSA Circle, San Antonio, TX, 78249, USA
| | - Sundeep Khosla
- Kogod Center on Aging and Division of Endocrinology, Mayo Clinic College of Medicine, Mayo Clinic, Rochester, MN, USA
| | - Christopher Chaput
- Department of Orthopedics, The University of Texas Health Science Center at San Antonio, San Antonio, TX, USA
| | - Hugo Giambini
- Department of Biomedical Engineering and Chemical Engineering, University of Texas at San Antonio, One UTSA Circle, San Antonio, TX, 78249, USA.
| |
Collapse
|
7
|
Gardegaront M, Allard V, Confavreux C, Bermond F, Mitton D, Follet H. Variabilities in µQCT-based FEA of a tumoral bone mice model. J Biomech 2021; 118:110265. [PMID: 33545571 DOI: 10.1016/j.jbiomech.2021.110265] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2020] [Revised: 01/04/2021] [Accepted: 01/16/2021] [Indexed: 01/13/2023]
Abstract
A finite element analysis based on Micro-Quantitative Computed Tomography (µQCT) is a method with high potential to improve fracture risk prediction. However, the segmentation process and model generation are generally not automatized in their entirety. Even with a rigorous protocol, the operator might add uncertainties during the creation of the model. The aim of this study was to evaluate a µQCT-based model of mice tumoral and sham tibias in terms of the variabilities induced by the operator and sensitivity to operator-dependent variables (such as model orientation or length). Two different operators generated finite element (FE) models from µCT images of 8 female Balb/c nude mice tibias aged 10 weeks old with bone tumors induced in the right tibia and with sham injection in the left. From these models, predicted failure load was determined for two different boundary conditions: fixed support and spherical joints. The difference between the predicted and experimental failure load of both operators was large (-122% to 93%). The difference in the predicted failure load between operators was less for the spherical joints boundary conditions (9.8%) than for the fixed support (58.3%), p < 0.001, whereas varying the orientation of bone tibia caused more variability for the fixed support boundary condition (44.7%) than for the spherical joints (9.1%), p < 0.002. Varying tibia length had no significant effect, regardless of boundary conditions (<4%). When using the same mesh and same orientation, the difference between operators is non-significant (<6%) for each model. This study showed that the operator influences the failure load assessed by a µQCT-based finite element model of the tumoral and sham mice tibias. The results suggest that automation is needed for better reproducibility.
Collapse
Affiliation(s)
- M Gardegaront
- Univ Lyon, Université Claude Bernard Lyon 1, INSERM, LYOS UMR 1033, 69008 Lyon, France
| | - V Allard
- Univ Lyon, Université Claude Bernard Lyon 1, INSERM, LYOS UMR 1033, 69008 Lyon, France
| | - C Confavreux
- Univ Lyon, Université Claude Bernard Lyon 1, INSERM, LYOS UMR 1033, 69008 Lyon, France; Centre Expert des Métastases et d'Oncologie Osseuses (CEMOS), Service de Rhumatologie Sud, Centre Hospitalier Lyon Sud, Hospices Civils de Lyon, Lyon, France
| | - F Bermond
- Univ Lyon, Université Claude Bernard Lyon 1, Univ Gustave Eiffel, IFSTTAR, LBMC UMR_T9406, 69622 Lyon, France
| | - D Mitton
- Univ Lyon, Université Claude Bernard Lyon 1, Univ Gustave Eiffel, IFSTTAR, LBMC UMR_T9406, 69622 Lyon, France
| | - H Follet
- Univ Lyon, Université Claude Bernard Lyon 1, INSERM, LYOS UMR 1033, 69008 Lyon, France.
| |
Collapse
|