1
|
Maroto AB, Navas P, Alfano F. Scaphoid numerical simulation of the critical loading until fracture. J Mech Behav Biomed Mater 2023; 148:106186. [PMID: 37856993 DOI: 10.1016/j.jmbbm.2023.106186] [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: 07/10/2023] [Revised: 10/04/2023] [Accepted: 10/10/2023] [Indexed: 10/21/2023]
Abstract
The numerical study of the scaphoid fracture, although it is relatively unexplored, can be of great clinical interest since it is highly common and can result in temporary or persistent disability. In this manuscript, seven combinations of boundary conditions and contacts between adjacent bones, together with four different loads, simulating real hand movements, are assessed. Three different fracture criteria for bones are employed to study the failure of the scaphoid with the aforementioned combination of interaction conditions. The results offer an interesting view of the accuracy of the possible interaction between adjacent bones. For future calculation, it would be possible to choose a combination of the balance between precision and computational cost savings. This study provides a comprehensive assessment into the modeling of the scaphoid bone and its interactions with adjacent bones. The findings reveal that various choices of interactions can yield similar results, allowing for flexibility in selecting interaction models based on desired accuracy or computational efficiency. Ultimately, this study establishes a foundational understanding for future research on modeling scaphoid motion.
Collapse
Affiliation(s)
- Ana B Maroto
- Continuum Mechanics and Theory of Structures Dep. Civil Eng. School, Universidad Politécnica de Madrid, Prof. Aranguren 3, 28040, Madrid, Spain; Biomedical Image Technologies, Universidad Politécnica de Madrid, Av. Complutense, 30, 28040, Madrid, Spain
| | - Pedro Navas
- Continuum Mechanics and Theory of Structures Dep. Civil Eng. School, Universidad Politécnica de Madrid, Prof. Aranguren 3, 28040, Madrid, Spain.
| | - Felicia Alfano
- Biomedical Image Technologies, Universidad Politécnica de Madrid, Av. Complutense, 30, 28040, Madrid, Spain; CIBER-BBN, ISCIII, Av. Monforte de Lemos, 3-5, 28029, Madrid, Spain
| |
Collapse
|
2
|
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
|
3
|
Abe S, Kouhia R, Nikander R, Narra N, Hyttinen J, Sievänen H. Effect of fall direction on the lower hip fracture risk in athletes with different loading histories: A finite element modeling study in multiple sideways fall configurations. Bone 2022; 158:116351. [PMID: 35131487 DOI: 10.1016/j.bone.2022.116351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/11/2021] [Revised: 02/01/2022] [Accepted: 02/01/2022] [Indexed: 11/24/2022]
Abstract
Physical loading makes bones stronger through structural adaptation. Finding effective modes of exercise to improve proximal femur strength has the potential to decrease hip fracture risk. Previous proximal femur finite element (FE) modeling studies have indicated that the loading history comprising impact exercises is associated with substantially higher fracture load. However, those results were limited only to one specified fall direction. It remains thus unclear whether exercise-induced higher fracture load depends on the fall direction. To address this, using magnetic resonance images of proximal femora from 91 female athletes (mean age 24.7 years with >8 years competitive career) and their 20 non-athletic but physically active controls (mean age 23.7 years), proximal femur FE models were created in 12 different sideways fall configurations. The athletes were divided into five groups by typical loading patterns of their sports: high-impact (H-I: 9 triple- and 10 high-jumpers), odd-impact (O-I: 9 soccer and 10 squash players), high-magnitude (H-M: 17 powerlifters), repetitive-impact (R-I: 18 endurance runners), and repetitive non-impact (R-NI: 18 swimmers). Compared to the controls, the FE models showed that the H-I and R-I groups had significantly (p < 0.05) higher fracture loads, 11-17% and 22-28% respectively, in all fall directions while the O-I group had significantly 10-11% higher fracture loads in four fall directions. The H-M and R-NI groups did not show significant benefit in any direction. Also, the analyses of the minimum fall strength (MFS) among these multiple fall configurations confirmed significantly 15%, 11%, and 14% higher MFSs in these impact groups, respectively, compared to the controls. These results suggest that the lower hip fracture risk indicated by higher fracture loads in athletes engaged in high impact or repetitive impact sports is independent of fall direction whereas the lower fracture risk attributed to odd-impact exercise is more modest and specific to the fall direction. Moreover, in concordance with the literature, the present study also confirmed that the fracture risk increases if the impact is imposed on the more posterolateral aspect of the hip. The present results highlight the importance of engaging in the impact exercises to prevent hip fractures and call for retrospective studies to investigate whether specific impact exercise history in adolescence and young adulthood is also associated with lower incidence of hip fractures in later life.
Collapse
Affiliation(s)
- Shinya Abe
- Structural Mechanics, Faculty of Built Environment, Tampere University, Tampere, Finland.
| | - Reijo Kouhia
- Structural Mechanics, Faculty of Built Environment, Tampere University, Tampere, Finland
| | - Riku Nikander
- Gerontology Research Center, Faculty of Sports Sciences, University of Jyväskylä, Jyväskylä, Finland; Central Hospital of Central Finland, Jyväskylä, Finland
| | - Nathaniel Narra
- BioMediTech Unit, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland
| | - Jari Hyttinen
- BioMediTech Unit, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland
| | - Harri Sievänen
- The UKK Institute for Health Promotion Research, Tampere, Finland
| |
Collapse
|
4
|
Amini M, Reisinger A, Hirtler L, Pahr D. Which experimental procedures influence the apparent proximal femoral stiffness? A parametric study. BMC Musculoskelet Disord 2021; 22:815. [PMID: 34556078 PMCID: PMC8461859 DOI: 10.1186/s12891-021-04656-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/16/2021] [Accepted: 08/26/2021] [Indexed: 11/10/2022] Open
Abstract
Background Experimental validation is the gold standard for the development of FE predictive models of bone. Employing multiple loading directions could improve this process. To capture the correct directional response of a sample, the effect of all influential parameters should be systematically considered. This study aims to determine the impact of common experimental parameters on the proximal femur’s apparent stiffness. Methods To that end, a parametric approach was taken to study the effects of: repetition, pre-loading, re-adjustment, re-fixation, storage, and μCT scanning as random sources of uncertainties, and loading direction as the controlled source of variation in both stand and side-fall configurations. Ten fresh-frozen proximal femoral specimens were prepared and tested with a novel setup in three consecutive sets of experiments. The neutral state and 15-degree abduction and adduction angles in both stance and fall configurations were tested for all samples and parameters. The apparent stiffness of the samples was measured using load-displacement data from the testing machine and validated against marker displacement data tracked by DIC cameras. Results Among the sources of uncertainties, only the storage cycle affected the proximal femoral apparent stiffness significantly. The random effects of setup manipulation and intermittent μCT scanning were negligible. The 15∘ deviation in loading direction had a significant effect comparable in size to that of switching the loading configuration from neutral stance to neutral side-fall. Conclusion According to these results, comparisons between the stiffness of the samples under various loading scenarios can be made if there are no storage intervals between the different load cases on the same samples. These outcomes could be used as guidance in defining a highly repeatable and multi-directional experimental validation study protocol.
Collapse
Affiliation(s)
- Morteza Amini
- Institute of Lightweight Design and Structural Biomechanics, TU Wien, Getreidemarkt 9, Vienna, 1060, Austria.,Division Biomechanics, Karl Landsteiner University of Health Sciences, Dr.-Karl-Dorrek-Straße 30, Krems an der Donau, 3500, Austria
| | - Andreas Reisinger
- Institute of Lightweight Design and Structural Biomechanics, TU Wien, Getreidemarkt 9, Vienna, 1060, Austria.,Division Biomechanics, Karl Landsteiner University of Health Sciences, Dr.-Karl-Dorrek-Straße 30, Krems an der Donau, 3500, Austria
| | - Lena Hirtler
- Center for Anatomy and Cell Biology, Medical University of Vienna, Währinger Straße 13, Vienna, 1090, Austria
| | - Dieter Pahr
- Institute of Lightweight Design and Structural Biomechanics, TU Wien, Getreidemarkt 9, Vienna, 1060, Austria. .,Division Biomechanics, Karl Landsteiner University of Health Sciences, Dr.-Karl-Dorrek-Straße 30, Krems an der Donau, 3500, Austria.
| |
Collapse
|
5
|
Real-time replication of three-dimensional and time-varying physiological loading cycles for bone and implant testing: A novel protocol demonstrated for the proximal human femur while walking. J Mech Behav Biomed Mater 2021; 124:104817. [PMID: 34536802 DOI: 10.1016/j.jmbbm.2021.104817] [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: 03/04/2021] [Revised: 08/26/2021] [Accepted: 09/03/2021] [Indexed: 11/24/2022]
Abstract
In vitro real-time replication of three-dimensional, time-varying load profiles acting on human bones during physical activity can advance bone and implant testing protocols. This study aimed to develop a novel protocol for applying the three-dimensional, time-varying hip contact force while walking to a human femur specimen. The target force profile was obtained from the literature. A proximal femur from an elderly female donor was instrumented using ten rosette strain gages and tested using a custom-made hexapod robot. A load-control algorithm determined the robot position generating the target force at low frequency (0.0004 Hz). Five cycles of the robot position were played back at five intermediate frequencies up to real-time (0.04, 0.08, 0.16, 0.4, and 0.8 Hz). The hip reaction force, the length of the actuators (position), and cortical strains were compared. The error in the load-control force was 0.3 ± 4.2 N (mean ± SD). The last three force, position, and strain cycles varied by less than 1.1% for every frequency analyzed. Across frequencies, the force increased by 28% at 0.8 Hz as a logarithmic function of frequency (R2 = 0.98). The position and strain error linearly increased with frequency up to 0.4 Hz. The median position error and the interquartile range of the strain error reached 15% and 13% at 0.8 Hz. Changes of force and cortical strain at increasing frequencies were linearly related (R2 = 0.99). Therefore, the protocol developed can provide repeatable three-dimensional time-varying load profiles, although the comparison of the specimen deformation obtained across frequencies should be considered with care, particularly in the higher frequency range. This information supports the design of dynamic tests of bone and implants.
Collapse
|
6
|
Bone Density Micro-CT Assessment during Embedding of the Innovative Multi-Spiked Connecting Scaffold in Periarticular Bone to Elaborate a Validated Numerical Model for Designing Biomimetic Fixation of Resurfacing Endoprostheses. MATERIALS 2021; 14:ma14061384. [PMID: 33809176 PMCID: PMC8000753 DOI: 10.3390/ma14061384] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/13/2021] [Revised: 02/26/2021] [Accepted: 03/09/2021] [Indexed: 11/21/2022]
Abstract
Our team has been working for some time on designing a new kind of biomimetic fixation of resurfacing endoprostheses, in which the innovative multi-spiked connecting scaffold (MSC-Scaffold) that mimics the natural interface between articular cartilage and periarticular trabecular bone in human joints is the crucial element. This work aimed to develop a numerical model enabling the design of the considered joint replacement implant that would reflect the mechanics of interacting biomaterials. Thus, quantitative micro-CT analysis of density distribution in bone material during the embedding of MSC-Scaffold in periarticular bone was applied. The performed numerical studies and corresponding mechanical tests revealed, under the embedded MSC-Scaffold, the bone material densification affecting its mechanical properties. On the basis of these findings, the built numerical model was modified by applying a simulated insert of densified bone material. This modification led to a strong correlation between the re-simulation and experimental results (FVU = 0.02). The biomimetism of the MSC-Scaffold prototype that provided physiological load transfer from implant to bone was confirmed based on the Huber–von Mises–Hencky (HMH) stress maps obtained with the validated finite element (FE) model of the problem. The micro-CT bone density assessment performed during the embedding of the MSC-Scaffold prototype in periarticular bone provides insight into the mechanical behaviour of the investigated implant-bone system and validates the numerical model that can be used for the design of material and geometric features of a new kind of resurfacing endoprostheses fixation.
Collapse
|
7
|
Mahmoudi M, Mahbadi H. Numerical investigation of mechanical behavior of human femoral diaphysis in normal and defective geometry: experimental evaluation. Comput Methods Biomech Biomed Engin 2020; 24:637-652. [PMID: 33164564 DOI: 10.1080/10255842.2020.1843639] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Abstract
Failure and major reoperation after internal fixation (IF) in mature femoral bones are common and proper selection of fixation method may reduce the rate of reoperations. Investigating the mechanical behavior of the human femoral diaphysis, this article studies effect of mechanical properties and geometry of the bone on selection of IF method. To this aim, we calculated the bone mineral density in human femurs, and then, using computed tomography scan, we obtained geometry and nonhomogeneous properties of the bone. Finite element (FE) models of osteotomised femurs were reinforced using four types of screws with a locking compression plate (LCP). We performed buckling and 4-point bending simulations, and results of these simulations represent critical buckling loads, maximum von Mises stresses, and strains around the screws and the central defect. To evaluate FE analysis, we employed the compressive experiments and compared load vs. displacement curves with FE results. Results corresponding to intact, osteotomised, and reinforced states are compared together, and the effect of cortical and unicortical screws in LCPs is studied. The FE results showed that application of identical prophylactic IF for two persons with identical injuries in the same conditions bring quite inverse results. As a consequence, evaluation of osteoporosis, elastic modulus, and morphometric data are required before fixation and screw selection. Besides, for short diaphysis, unicortical screws have maximum strengthening factor in bending. While for long samples, these types of screws can be the worst option, application of cortical screws results to maximum strength in comparison with other types.
Collapse
Affiliation(s)
- Moeinoddin Mahmoudi
- Young Researchers and Elite Club, Central Tehran Branch, Islamic Azad University, Tehran, Iran
| | - Hossein Mahbadi
- Department of Mechanical Engineering, Central Tehran Branch, Islamic Azad University, Tehran, Iran
| |
Collapse
|
8
|
Martel DR, Lysy M, Laing AC. Predicting population level hip fracture risk: a novel hierarchical model incorporating probabilistic approaches and factor of risk principles. Comput Methods Biomech Biomed Engin 2020; 23:1201-1214. [PMID: 32687412 DOI: 10.1080/10255842.2020.1793331] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Fall-related hip fractures are a major public health issue. While individual-level risk assessment tools exist, population-level predictive models could catalyze innovation in large-scale interventions. This study presents a hierarchical probabilistic model that predicts population-level hip fracture risk based on Factor of Risk (FOR) principles. Model validation demonstrated that FOR output aligned with a published dataset categorized by sex and hip fracture status. The model predicted normalized FOR for 100000 individuals simulating the Canadian older-adult population. Predicted hip fracture risk was higher for females (by an average of 38%), and increased with age (by15% per decade). Potential applications are discussed.
Collapse
Affiliation(s)
- Daniel R Martel
- Department of Kinesiology, University of Waterloo, Waterloo, Ontario, Canada
| | - Martin Lysy
- Department of Statistics and Actuarial Science, University of Waterloo, Waterloo, Ontario, Canada
| | - Andrew C Laing
- Department of Kinesiology, University of Waterloo, Waterloo, Ontario, Canada
| |
Collapse
|
9
|
Falcinelli C, Whyne C. Image-based finite-element modeling of the human femur. Comput Methods Biomech Biomed Engin 2020; 23:1138-1161. [PMID: 32657148 DOI: 10.1080/10255842.2020.1789863] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Fracture is considered a critical clinical endpoint in skeletal pathologies including osteoporosis and bone metastases. However, current clinical guidelines are limited with respect to identifying cases at high risk of fracture, as they do not account for many mechanical determinants that contribute to bone fracture. Improving fracture risk assessment is an important area of research with clear clinical relevance. Patient-specific numerical musculoskeletal models generated from diagnostic images are widely used in biomechanics research and may provide the foundation for clinical tools used to quantify fracture risk. However, prior to clinical translation, in vitro validation of predictions generated from such numerical models is necessary. Despite adopting radically different models, in vitro validation of image-based finite element (FE) models of the proximal femur (predicting strains and failure loads) have shown very similar, encouraging levels of accuracy. The accuracy of such in vitro models has motivated their application to clinical studies of osteoporotic and metastatic fractures. Such models have demonstrated promising but heterogeneous results, which may be explained by the lack of a uniform strategy with respect to FE modeling of the human femur. This review aims to critically discuss the state of the art of image-based femoral FE modeling strategies, highlighting principal features and differences among current approaches. Quantitative results are also reported with respect to the level of accuracy achieved from in vitro evaluations and clinical applications and are used to motivate the adoption of a standardized approach/workflow for image-based FE modeling of the femur.
Collapse
Affiliation(s)
- Cristina Falcinelli
- Orthopaedic Biomechanics Laboratory, Sunnybrook Research Institute, Toronto, Canada
| | - Cari Whyne
- Orthopaedic Biomechanics Laboratory, Sunnybrook Research Institute, Toronto, Canada
| |
Collapse
|
10
|
Toniolo I, Salmaso C, Bruno G, De Stefani A, Stefanini C, Gracco ALT, Carniel EL. Anisotropic computational modelling of bony structures from CT data: An almost automatic procedure. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2020; 189:105319. [PMID: 31951872 DOI: 10.1016/j.cmpb.2020.105319] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/31/2019] [Revised: 11/27/2019] [Accepted: 01/05/2020] [Indexed: 06/10/2023]
Abstract
BACKGROUND AND OBJECTIVE The use of modelling techniques that combine CT data and bone tissue micromechanics is spreading in computational biomechanics. Finite Element models show great potential in surgical planning of intervention and in prediction of stress and strain fields through a non-invasive method. The main challenge pertains to the reliable characterization of bone mechanical behaviour. An almost automatic procedure is here defined, which provides computational models of bony structures considering the actual anisotropy of bone tissue response. The innovative aspect resides on the automatic detection of the directions of anisotropy as the eigenvectors of a three-dimensional distribution matrix of HU values. METHODS The procedure combines CT data and micromechanics modelling techniques. Regarding a specific location, the procedure reports both the orthotropic elastic constants, by the analysis of the local HU value, and the anisotropic material directions, by the analysis of the HU values distribution around the specific location. RESULTS The procedure returns the distribution of bone tissue orthotropic elasticity tensor. The procedure proves to correctly respect the differentiation between cortical and trabecular bone. Principal directions show to be consistent with experimental data from ultrasound measurements. Regarding the material mapping from voxel to FE model, the developed strategies show to be reliable, leading to marginal errors (lower than 10%) for most of CT voxels (more than 90%). The computational analyses of typical structural loading conditions lead to strain values that are comparable with results from strain gauges experimentations. The development and the exploitation of FE models of different bony structures allow assessing the reliability of the procedure for cortical bone. CONCLUSIONS The results highlight the potentialities of the procedure in providing accurate patient-specific biomechanical models of bony structures starting from CT data. The accuracy and the automatism of the procedure are important factors for the development of real time clinical tools. The main limitations of this work remain the not fully automatism and the reliability assessment, which is based mainly on cortical bone regions only.
Collapse
Affiliation(s)
- Ilaria Toniolo
- The BioRobotics Institute, Scuola Superiore Sant'Anna, Pisa, Italy; Department of Industrial Engineering, University of Padova, Italy.
| | - Claudia Salmaso
- Department of Industrial Engineering, University of Padova, Italy
| | - Giovanni Bruno
- Department of Neurosciences, University of Padova, Italy
| | | | - Cesare Stefanini
- The BioRobotics Institute, Scuola Superiore Sant'Anna, Pisa, Italy
| | | | - Emanuele Luigi Carniel
- Department of Industrial Engineering, University of Padova, Italy; Centre for Mechanics of Biological Materials, University of Padova, Italy
| |
Collapse
|
11
|
Kalajahi SMH, Nazemi SM, Johnston JD. An exclusion approach for addressing partial volume artifacts with quantititive computed tomography-based finite element modeling of the proximal tibia. Med Eng Phys 2019; 76:95-100. [PMID: 31870545 DOI: 10.1016/j.medengphy.2019.10.013] [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/25/2018] [Revised: 10/14/2019] [Accepted: 10/20/2019] [Indexed: 10/25/2022]
Abstract
INTRODUCTION Quantitative computed tomography based finite element modeling (QCT-FE) has potential to clarify the role of subchondral bone stiffness in osteoarthritis. The limited spatial resolution of clinical QCT systems, however, results in partial volume (PV) artifacts and low contrast between cortical and trabecular bone, which adversely affects the accuracy of QCT-FE models. The objective of this research was to evaluate the agreement between stiffness predictions offered by QCT-FE models of proximal tibial subchondral bone (constructed with and without a new voxel-exclusion algorithm) with experimentally-derived local subchondral bone structural stiffness. METHODS Thirteen proximal tibial compartments were obtained and imaged using QCT. Two types of QCT-FE models were developed: (1) standard model, which employed the standard procedure for QCT-FE modeling; and (2) "voxel exclusion (VE)" model, which addressed PV artifacts by excluding low density voxels during the material mapping stage of construction. We assessed agreement between QCT-FE stiffness estimates (using standard and VE approaches) with experimental stiffness by reporting predicted variance from linear regression and mean bias with 95% Limits of Agreement (LOA). RESULTS The standard and VE models explained 81% and 84% of the variance in experimentally measured stiffness, respectively. The standard model showed a mean bias of -268 N/mm (LOA -1210 to 679 N/mm); the VE model showed a mean bias of +59 N/mm (LOA -762 to 910 N/mm). INTERPRETATION The VE model explained more variance in subchondral bone stiffness with less bias. Our findings indicate that the VE method has potential to improve QCT-FE models of bone affected by PV artifacts.
Collapse
Affiliation(s)
| | - S Majid Nazemi
- Biomechanics Research Unit, GIGA In Silico Medicine, University of Liège, Liège, Belgium
| | - James D Johnston
- Department of Mechanical Engineering, University of Saskatchewan, Saskatoon, Canada.
| |
Collapse
|
12
|
Alcântara ACS, Assis I, Prada D, Mehle K, Schwan S, Costa-Paiva L, Skaf MS, Wrobel LC, Sollero P. Patient-Specific Bone Multiscale Modelling, Fracture Simulation and Risk Analysis-A Survey. MATERIALS (BASEL, SWITZERLAND) 2019; 13:E106. [PMID: 31878356 PMCID: PMC6981613 DOI: 10.3390/ma13010106] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/13/2019] [Revised: 12/16/2019] [Accepted: 12/17/2019] [Indexed: 12/26/2022]
Abstract
This paper provides a starting point for researchers and practitioners from biology, medicine, physics and engineering who can benefit from an up-to-date literature survey on patient-specific bone fracture modelling, simulation and risk analysis. This survey hints at a framework for devising realistic patient-specific bone fracture simulations. This paper has 18 sections: Section 1 presents the main interested parties; Section 2 explains the organzation of the text; Section 3 motivates further work on patient-specific bone fracture simulation; Section 4 motivates this survey; Section 5 concerns the collection of bibliographical references; Section 6 motivates the physico-mathematical approach to bone fracture; Section 7 presents the modelling of bone as a continuum; Section 8 categorizes the surveyed literature into a continuum mechanics framework; Section 9 concerns the computational modelling of bone geometry; Section 10 concerns the estimation of bone mechanical properties; Section 11 concerns the selection of boundary conditions representative of bone trauma; Section 12 concerns bone fracture simulation; Section 13 presents the multiscale structure of bone; Section 14 concerns the multiscale mathematical modelling of bone; Section 15 concerns the experimental validation of bone fracture simulations; Section 16 concerns bone fracture risk assessment. Lastly, glossaries for symbols, acronyms, and physico-mathematical terms are provided.
Collapse
Affiliation(s)
- Amadeus C. S. Alcântara
- Department of Computational Mechanics, School of Mechanical Engineering, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-860, Brazil; (A.C.S.A.); (D.P.)
| | - Israel Assis
- Department of Integrated Systems, School of Mechanical Engineering, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-860, Brazil;
| | - Daniel Prada
- Department of Computational Mechanics, School of Mechanical Engineering, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-860, Brazil; (A.C.S.A.); (D.P.)
| | - Konrad Mehle
- Department of Engineering and Natural Sciences, University of Applied Sciences Merseburg, 06217 Merseburg, Germany;
| | - Stefan Schwan
- Fraunhofer Institute for Microstructure of Materials and Systems IMWS, 06120 Halle/Saale, Germany;
| | - Lúcia Costa-Paiva
- Department of Obstetrics and Gynecology, School of Medical Sciences, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-887, Brazil;
| | - Munir S. Skaf
- Institute of Chemistry and Center for Computing in Engineering and Sciences, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-860, Brazil;
| | - Luiz C. Wrobel
- Institute of Materials and Manufacturing, Brunel University London, Uxbridge UB8 3PH, UK;
- Department of Civil and Environmental Engineering, Pontifical Catholic University of Rio de Janeiro, Rio de Janeiro 22451-900, Brazil
| | - Paulo Sollero
- Department of Computational Mechanics, School of Mechanical Engineering, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-860, Brazil; (A.C.S.A.); (D.P.)
| |
Collapse
|
13
|
Hosseini Kalajahi SM, Nazemi SM, Johnston JD. Separate modeling of cortical and trabecular bone offers little improvement in FE predictions of local structural stiffness at the proximal tibia. Comput Methods Biomech Biomed Engin 2019; 22:1258-1268. [DOI: 10.1080/10255842.2019.1661386] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
| | - S. Majid Nazemi
- Department of Mechanical Engineering, University of Saskatchewan, Saskatoon, Canada
| | - James D. Johnston
- Department of Mechanical Engineering, University of Saskatchewan, Saskatoon, Canada
| |
Collapse
|
14
|
Marco M, Giner E, Caeiro-Rey JR, Miguélez MH, Larraínzar-Garijo R. Numerical modelling of hip fracture patterns in human femur. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2019; 173:67-75. [PMID: 31046997 DOI: 10.1016/j.cmpb.2019.03.010] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2019] [Revised: 03/12/2019] [Accepted: 03/13/2019] [Indexed: 06/09/2023]
Abstract
BACKGROUND AND OBJECTIVE Hip fracture morphology is an important factor determining the ulterior surgical repair and treatment, because of the dependence of the treatment on fracture morphology. Although numerical modelling can be a valuable tool for fracture prediction, the simulation of femur fracture is not simple due to the complexity of bone architecture and the numerical techniques required for simulation of crack propagation. Numerical models assuming homogeneous fracture mechanical properties commonly fail in the prediction of fracture patterns. This paper focuses on the prediction of femur fracture based on the development of a finite element model able to simulate the generation of long crack paths. METHODS The finite element model developed in this work demonstrates the capability of predicting fracture patterns under stance loading configuration, allowing the distinction between the main fracture paths: intracapsular and extracapsular fractures. It is worth noting the prediction of different fracture patterns for the same loading conditions, as observed during experimental tests. RESULTS AND CONCLUSIONS The internal distribution of bone mineral density and femur geometry strongly influences the femur fracture morphology and fracture load. Experimental fracture paths have been analysed by means of micro-computed tomography allowing the comparison of predicted and experimental crack surfaces, confirming the good accuracy of the numerical model.
Collapse
Affiliation(s)
- Miguel Marco
- Department of Mechanical Engineering, Universidad Carlos III de Madrid, Avda. de la Universidad 30, 28911 Leganés, Madrid, Spain.
| | - Eugenio Giner
- CIIM-Department of Mechanical and Materials Engineering, Universitat Politècnica de València Camino de Vera, 46022 Valencia, Spain
| | - José Ramón Caeiro-Rey
- Orthopedic Surgery and Traumatology Service, Complejo Hospitalario Universitario de Santiago de Compostela, Rúa de Ramón Baltar, s/n, 15706 Santiago de Compostela, A Coruña, Spain
| | - M Henar Miguélez
- Department of Mechanical Engineering, Universidad Carlos III de Madrid, Avda. de la Universidad 30, 28911 Leganés, Madrid, Spain
| | - Ricardo Larraínzar-Garijo
- Orthopaedics and Trauma Department, Surgery Department, Hospital Universitario Infanta Leonor, Complutense University, Madrid, Spain
| |
Collapse
|
15
|
Huang L, Korhonen RK, Turunen MJ, Finnilä MAJ. Experimental mechanical strain measurement of tissues. PeerJ 2019; 7:e6545. [PMID: 30867989 PMCID: PMC6409087 DOI: 10.7717/peerj.6545] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2018] [Accepted: 01/31/2019] [Indexed: 12/22/2022] Open
Abstract
Strain, an important biomechanical factor, occurs at different scales from molecules and cells to tissues and organs in physiological conditions. Under mechanical strain, the strength of tissues and their micro- and nanocomponents, the structure, proliferation, differentiation and apoptosis of cells and even the cytokines expressed by cells probably shift. Thus, the measurement of mechanical strain (i.e., relative displacement or deformation) is critical to understand functional changes in tissues, and to elucidate basic relationships between mechanical loading and tissue response. In the last decades, a great number of methods have been developed and applied to measure the deformations and mechanical strains in tissues comprising bone, tendon, ligament, muscle and brain as well as blood vessels. In this article, we have reviewed the mechanical strain measurement from six aspects: electro-based, light-based, ultrasound-based, magnetic resonance-based and computed tomography-based techniques, and the texture correlation-based image processing method. The review may help solving the problems of experimental and mechanical strain measurement of tissues under different measurement environments.
Collapse
Affiliation(s)
- Lingwei Huang
- Department of Applied Physics, University of Eastern Finland, Kuopio, Finland
| | - Rami K Korhonen
- Department of Applied Physics, University of Eastern Finland, Kuopio, Finland
| | - Mikael J Turunen
- Department of Applied Physics, University of Eastern Finland, Kuopio, Finland
| | - Mikko A J Finnilä
- Department of Applied Physics, University of Eastern Finland, Kuopio, Finland.,Research Unit of Medical Imaging, Physics and Technology, University of Oulu, Oulu, Finland.,Medical Research Center Oulu, Oulu University Hospital, Oulu, Finland
| |
Collapse
|
16
|
Knowles NK, Langohr GDG, Faieghi M, Nelson AJ, Ferreira LM. A comparison of density-modulus relationships used in finite element modeling of the shoulder. Med Eng Phys 2019; 66:40-46. [PMID: 30833224 DOI: 10.1016/j.medengphy.2019.02.005] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2018] [Revised: 11/29/2018] [Accepted: 02/10/2019] [Indexed: 10/27/2022]
Abstract
Subject- and site-specific modeling techniques greatly improve the accuracy of computational models derived from clinical-resolution quantitative computed tomography (QCT) data. The majority of shoulder finite element (FE) studies use density-modulus relationships developed for alternative anatomical locations. As such, the objectives of this study were to compare the six most commonly used density-modulus relationships in shoulder finite element (FE) studies. To achieve this, ninety-eight (98) virtual trabecular bone cores were extracted from uCT scans of scapulae from 14 cadaveric specimens (7 male; 7 female). Homogeneous tissue moduli of 20 GPa, and heterogeneous tissue moduli scaled by CT-intensity were considered. Micro finite element models (µ-FEMs) of each virtual core were compressively loaded to 0.5% apparent strain and apparent strain energy density (SEDapp) was collected. Each uCT virtual core was then co-registered to clinical QCT images, QCT-FEMs created, and each of the 6 density-modulus relationships applied (6 × 98 = 588 QCT-FEMs). The loading and boundary conditions were replicated and SEDapp was collected and compared to µ-FEM SEDapp. When a homogeneous tissue modulus was considered in the µ-FEMs, SEDapp was best predicted in QCT-FEMs with the density-modulus relationship developed from pooled anatomical locations (QCT-FEM SEDapp = 0.979µ-FEM SEDapp + 0.0066, r2 = 0.933). A different density-modulus relationship best predicted SEDapp (QCT-FEM SEDapp = 1.014µ-FEM SEDapp + 0.0034, r2 = 0.935) when a heterogeneous tissue modulus was considered. This study compared density-modulus relationships used in shoulder FE studies using an independent computational methodology for comparing these relationships.
Collapse
Affiliation(s)
- Nikolas K Knowles
- School of Biomedical Engineering, The University of Western Ontario, London, ON, Canada; Roth
- McFarlane Hand and Upper Limb Centre, Bioengineering Laboratory, Surgical Mechatronics Laboratory, St. Josephs Health Care, 268 Grosvenor St., London, ON, Canada; Collaborative Training Program in MSK Health Research, and Bone and Joint Institute, The University of Western Ontario, London, ON, Canada
| | - G Daniel G Langohr
- Department of Mechanical and Materials Engineering, The University of Western Ontario, London, ON, Canada; Roth
- McFarlane Hand and Upper Limb Centre, Bioengineering Laboratory, Surgical Mechatronics Laboratory, St. Josephs Health Care, 268 Grosvenor St., London, ON, Canada; Collaborative Training Program in MSK Health Research, and Bone and Joint Institute, The University of Western Ontario, London, ON, Canada
| | - Mohammadreza Faieghi
- School of Biomedical Engineering, The University of Western Ontario, London, ON, Canada
| | - Andrew J Nelson
- Department of Anthropology, The University of Western Ontario, London, ON, Canada; Collaborative Training Program in MSK Health Research, and Bone and Joint Institute, The University of Western Ontario, London, ON, Canada; Department of Chemistry, The University of Western Ontario, London, ON, Canada
| | - Louis M Ferreira
- Department of Mechanical and Materials Engineering, The University of Western Ontario, London, ON, Canada; Roth
- McFarlane Hand and Upper Limb Centre, Bioengineering Laboratory, Surgical Mechatronics Laboratory, St. Josephs Health Care, 268 Grosvenor St., London, ON, Canada; Collaborative Training Program in MSK Health Research, and Bone and Joint Institute, The University of Western Ontario, London, ON, Canada.
| |
Collapse
|
17
|
Katz Y, Dahan G, Sosna J, Shelef I, Cherniavsky E, Yosibash Z. Scanner influence on the mechanical response of QCT-based finite element analysis of long bones. J Biomech 2019; 86:149-159. [DOI: 10.1016/j.jbiomech.2019.01.049] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2018] [Revised: 01/03/2019] [Accepted: 01/30/2019] [Indexed: 01/30/2023]
|
18
|
Comparison of HR-pQCT- and microCT-based finite element models for the estimation of the mechanical properties of the calcaneus trabecular bone. Biomech Model Mechanobiol 2018; 17:1715-1730. [DOI: 10.1007/s10237-018-1051-6] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2017] [Accepted: 07/02/2018] [Indexed: 12/13/2022]
|
19
|
Panyasantisuk J, Dall'Ara E, Pretterklieber M, Pahr DH, Zysset PK. Mapping anisotropy improves QCT-based finite element estimation of hip strength in pooled stance and side-fall load configurations. Med Eng Phys 2018; 59:36-42. [PMID: 30131112 DOI: 10.1016/j.medengphy.2018.06.004] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2017] [Revised: 03/26/2018] [Accepted: 06/24/2018] [Indexed: 02/05/2023]
Abstract
Hip fractures are one of the most severe consequences of osteoporosis. Compared to the clinical standard of DXA-based aBMD at the femoral neck, QCT-based FEA delivers a better surrogate of femoral strength and gains acceptance for the calculation of hip fracture risk when a CT reconstruction is available. Isotropic, homogenised voxel-based, finite element (hvFE) models are widely used to estimate femoral strength in cross-sectional and longitudinal clinical studies. However, fabric anisotropy is a classical feature of the architecture of the proximal femur and the second determinant of the homogenised mechanical properties of trabecular bone. Due to the limited resolution, fabric anisotropy cannot be derived from clinical CT reconstructions. Alternatively, fabric anisotropy can be extracted from HR-pQCT images of cadaveric femora. In this study, fabric anisotropy from HR-pQCT images was mapped onto QCT-based hvFE models of 71 human proximal femora for which both HR-pQCT and QCT images were available. Stiffness and ultimate load computed from anisotropic hvFE models were compared with previous biomechanical tests in both stance and side-fall configurations. The influence of using the femur-specific versus a mean fabric distribution on the hvFE predictions was assessed. Femur-specific and mean fabric enhance the prediction of experimental ultimate force for the pooled, i.e. stance and side-fall, (isotropic: r2=0.81, femur-specific fabric: r2=0.88, mean fabric: r2=0.86,p<0.001) but not for the individual configurations. Fabric anisotropy significantly improves bone strength prediction for the pooled configurations, and mapped fabric provides a comparable prediction to true fabric. The mapping of fabric anisotropy is therefore expected to help generate more accurate QCT-based hvFE models of the proximal femur for personalised or multiple load configurations.
Collapse
Affiliation(s)
- J Panyasantisuk
- Institute for Surgical Technology and Biomechanics, University of Bern, Switzerland
| | - E Dall'Ara
- Department of Oncology and Metabolism and INSIGNEO, Institute for in silico Medicine, University of Sheffield, United Kingdom
| | | | - D H Pahr
- Institute for Lightweight Design and Structural Biomechanics, Vienna University of Technology, Austria; Department for Anatomy and Biomechanics, Karl Landsteiner Private University for Health Sciences, Austria
| | - P K Zysset
- Institute for Surgical Technology and Biomechanics, University of Bern, Switzerland.
| |
Collapse
|
20
|
Enns-Bray WS, Ferguson SJ, Helgason B. Strain rate dependency of bovine trabecular bone under impact loading at sideways fall velocity. J Biomech 2018; 75:46-52. [DOI: 10.1016/j.jbiomech.2018.04.042] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2018] [Revised: 04/23/2018] [Accepted: 04/23/2018] [Indexed: 11/16/2022]
|
21
|
Enns-Bray W, Bahaloo H, Fleps I, Ariza O, Gilchrist S, Widmer R, Guy P, Pálsson H, Ferguson S, Cripton P, Helgason B. Material mapping strategy to improve the predicted response of the proximal femur to a sideways fall impact. J Mech Behav Biomed Mater 2018; 78:196-205. [DOI: 10.1016/j.jmbbm.2017.10.033] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2016] [Revised: 09/25/2017] [Accepted: 10/26/2017] [Indexed: 11/29/2022]
|
22
|
On the Failure Initiation in the Proximal Human Femur Under Simulated Sideways Fall. Ann Biomed Eng 2017; 46:270-283. [DOI: 10.1007/s10439-017-1952-z] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2017] [Accepted: 10/30/2017] [Indexed: 11/25/2022]
|