1
|
Martel DR, Callaghan JP, Mourtzakis M, Willett TL, Laing AC. Influence of test paradigm on loading dynamics during proximal femur fracture tests simulating sideways falls. J Mech Behav Biomed Mater 2024; 157:106631. [PMID: 38986216 DOI: 10.1016/j.jmbbm.2024.106631] [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: 02/15/2024] [Revised: 06/10/2024] [Accepted: 06/12/2024] [Indexed: 07/12/2024]
Abstract
Fall-related hip fractures are a serious public health issue in older adults. As most mechanistic hip fracture risk prediction models incorporate tissue tolerance, test methods that can accurately characterize the fracture force of the femur (and factors that influence it) are imperative. While bone possesses viscoelastic properties, experimental characterization of rate-dependencies has been inconsistent in the whole-femur literature. The goal of this study was to investigate the influence of experimental paradigm on loading rate and fracture force (both means and variability) during mechanical tests simulating lateral fall loadings on the proximal femur. Six pairs of matched femurs were split randomly between two test paradigms: a 'lower rate' materials testing system (MTS) with a constant displacement rate of 60 mm/s, and a hip impact test system (HIT) comprised of a custom-built vertical drop tower utilizing an impact velocity of 4 m/s. The loading rate was 88-fold higher for the HIT (mean (SD) = 2465.49 (807.38) kN/s) compared to the MTS (27.78 (10.03) kN/s) paradigm. However, no difference in fracture force was observed between test paradigms (mean (SD) = 4096.4 (1272.6) N for HIT, and 3641.3 (1285.8) N for MTS). Within-paradigm variability was not significantly different across paradigms for either loading rate or fracture force (coefficients of variation ranging from 0.311 to 0.361). Within each test paradigm, significant positive relationships were observed between loading rate and fracture force (HIT adjusted R2 = 0.833, p = 0.007; MTS adjusted R2 = 0.983, p < 0.0001). Overall, this study provides evidence that energy-based impact simulators can be a valid method to measure femoral bone strength in the context of fall-related hip fractures. This study motivates future research to characterize potential non-linear relationships between loading rate and fracture threshold at both macro and microscales.
Collapse
Affiliation(s)
- Daniel R Martel
- University of Waterloo, Department of Kinesiology and Health Sciences, Waterloo, ON, Canada
| | - Jack P Callaghan
- University of Waterloo, Department of Kinesiology and Health Sciences, Waterloo, ON, Canada
| | - Marina Mourtzakis
- University of Waterloo, Department of Kinesiology and Health Sciences, Waterloo, ON, Canada
| | - Thomas L Willett
- University of Waterloo, Department of Systems Design Engineering, Waterloo, ON, Canada
| | - Andrew C Laing
- University of Waterloo, Department of Kinesiology and Health Sciences, Waterloo, ON, Canada.
| |
Collapse
|
2
|
Carlsson J, Braesch-Andersen A, Ferguson SJ, Isaksson P. Fracture in porous bone analysed with a numerical phase-field dynamical model. J Mech Behav Biomed Mater 2023; 139:105659. [PMID: 36638634 DOI: 10.1016/j.jmbbm.2023.105659] [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/04/2022] [Revised: 12/13/2022] [Accepted: 01/04/2023] [Indexed: 01/07/2023]
Abstract
A dynamic phase-field fracture finite element model is applied to discretized high-resolution three-dimensional computed tomography images of human trabecular bone to analyse rapid bone fracture. The model is contrasted to quasi-static experimental results and a quasi-static phase-field finite element model. The experiment revealed complex stepwise crack evolution with multiple crack fronts, and crack arrests, as the global tensile displacement load was incrementally increased. The quasi-static phase-field fracture model captures the fractures in the experiment reasonably well, and the dynamic model converges towards the quasi-static model when mechanically loaded at low rates. At higher load rates, i.e., at larger impulses, inertia effects significantly contribute to an increased initial global stiffness, higher peak forces and a larger number of cracks spread over a larger volume. Since the fracture process clearly is different at large impulses compared to small impulses, it is concluded that dynamic fracture models are necessary when simulating rapid bone fracture.
Collapse
Affiliation(s)
- Jenny Carlsson
- Solid Mechanics, Department of Materials Science and Engineering, Uppsala University, Sweden; Now at Cambridge University Engineering Department, Trumpington St., Cambridge, UK
| | - Anna Braesch-Andersen
- Solid Mechanics, Department of Materials Science and Engineering, Uppsala University, Sweden
| | | | - Per Isaksson
- Solid Mechanics, Department of Materials Science and Engineering, Uppsala University, Sweden.
| |
Collapse
|
3
|
The Influence of Static Load and Sideways Impact Fall on Extramedullary Bone Plates Used to Treat Intertrochanteric Femoral Fracture: A Preclinical Strength Assessment. Ann Biomed Eng 2022; 50:1923-1940. [PMID: 35821164 DOI: 10.1007/s10439-022-03013-z] [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: 02/08/2022] [Accepted: 07/06/2022] [Indexed: 12/30/2022]
Abstract
Hip fracture accounts for a large number of hospitalizations, thereby causing substantial economic burden. Majority (> 90%) of all hip fractures are associated to sideways fall. Studies on sideways fall usually involve loading at quasi-static or at constant displacement rate, which neglects the physics of actual fall. Understanding femur resonance frequency and associated mode shapes excited by dynamic loads is also critical. Two commercial extramedullary implants, proximal femoral locking plate (PFLP) and variable angle dynamic hip screw (VA-DHS), were chosen to carry out the preclinical assessments on a simulated Evans-I type intertrochanteric fracture. In this study, we hypothesized that the behavior of the implant depends on the loading types-axial static and transverse impact-and a rigid implanted construct will absorb less impact energy for sideways fall. The in silico models were validated using experimental measurements of full-field strain data obtained from a 2D digital image correlation (DIC) study. Under peak axial load of 3 kN, PFLP construct predicted greater axial stiffness (1.07 kN/mm) as opposed to VA-DHS (0.85 kN/mm), although the former predicted slightly higher proximal stress shielding. Further, with greater mode 2 frequency, PFLP predicted improved performance in resisting bending due to sideways fall as compared to the other implant. Overall, the PFLP implanted femur predicted the least propensity to adverse stress intensities, suggesting better structural rigidity and higher capacity in protecting the fractured femur against fall.
Collapse
|
4
|
The influence of foramina on femoral neck fractures and strains predicted with finite element analysis. J Mech Behav Biomed Mater 2022; 134:105364. [DOI: 10.1016/j.jmbbm.2022.105364] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2022] [Revised: 07/07/2022] [Accepted: 07/09/2022] [Indexed: 11/21/2022]
|
5
|
Revel M, Gardegaront M, Bermond F, Mitton D, Follet H. A credible homogenized finite element model to predict radius fracture in the case of a forward fall. J Mech Behav Biomed Mater 2022; 131:105206. [DOI: 10.1016/j.jmbbm.2022.105206] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2021] [Revised: 03/23/2022] [Accepted: 03/26/2022] [Indexed: 11/15/2022]
|
6
|
Ionic liquid treatment for efficient sample preparation of hydrated bone for scanning electron microscopy. Micron 2021; 153:103192. [PMID: 34896685 DOI: 10.1016/j.micron.2021.103192] [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: 09/08/2021] [Revised: 11/22/2021] [Accepted: 11/23/2021] [Indexed: 11/20/2022]
Abstract
This study presents a new protocol for preparing bone samples for scanning electron microscopy (SEM) using a room temperature ionic liquid (RTIL) treatment method. RTIL-based solutions can be adopted as an alternative to lengthy and laborious traditional means of preparation for SEM due to their unique low-vapour pressure and conductive properties. Applied to biological samples, RTILs can be used quickly and efficiently to observe hydrated, unfixed structures in typical SEM systems. This first-time feasibility study of the optimization of this protocol for bone was explored through various SEM modalities using two distinct ionic liquids, 1-ethyl-3-methylimidazolium tetrafluoroborate ([EMI][BF4]) and 1-butyl-3-methyl imidazolium tetrafluoroborate ([BMI][BF4]), at varying concentrations of 5, 10, and 25 % v/v in aqueous solution through an addition-based method. Based on qualitative observations in the SEM, a 60-second solution addition treatment of 10 % v/v [BMI][BF4] performed the best in imaging hydrated, unfixed bone samples, resulting in minimal charge buildup and no solution pooling on the surface. The treatment was applied effectively to a variety of bone samples, notably flat and polished, as well as highly topographical bone fracture surfaces of both healthy and osteoporotic human bone samples. In comparison to conventionally dehydrated bone, the RTIL treatment better preserved the natural bone structure, resulting in minimal microcracking in observed structures.
Collapse
|
7
|
Bhattacharya P, Li Q, Lacroix D, Kadirkamanathan V, Viceconti M. A systematic approach to the scale separation problem in the development of multiscale models. PLoS One 2021; 16:e0251297. [PMID: 34003842 PMCID: PMC8130972 DOI: 10.1371/journal.pone.0251297] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2020] [Accepted: 04/25/2021] [Indexed: 11/19/2022] Open
Abstract
Throughout engineering there are problems where it is required to predict a quantity based on the measurement of another, but where the two quantities possess characteristic variations over vastly different ranges of time and space. Among the many challenges posed by such 'multiscale' problems, that of defining a 'scale' remains poorly addressed. This fundamental problem has led to much confusion in the field of biomedical engineering in particular. The present study proposes a definition of scale based on measurement limitations of existing instruments, available computational power, and on the ranges of time and space over which quantities of interest vary characteristically. The definition is used to construct a multiscale modelling methodology from start to finish, beginning with a description of the system (portion of reality of interest) and ending with an algorithmic orchestration of mathematical models at different scales within the system. The methodology is illustrated for a specific but well-researched problem. The concept of scale and the multiscale modelling approach introduced are shown to be easily adaptable to other closely related problems. Although out of the scope of this paper, we believe that the proposed methodology can be applied widely throughout engineering.
Collapse
Affiliation(s)
- Pinaki Bhattacharya
- Department of Mechanical Engineering, University of Sheffield, Sheffield, United Kingdom
- INSIGNEO Institute for in silico Medicine, University of Sheffield, Sheffield, United Kingdom
- * E-mail:
| | - Qiao Li
- Department of Mechanical Engineering, University of Sheffield, Sheffield, United Kingdom
- INSIGNEO Institute for in silico Medicine, University of Sheffield, Sheffield, United Kingdom
| | - Damien Lacroix
- Department of Mechanical Engineering, University of Sheffield, Sheffield, United Kingdom
- INSIGNEO Institute for in silico Medicine, University of Sheffield, Sheffield, United Kingdom
| | - Visakan Kadirkamanathan
- INSIGNEO Institute for in silico Medicine, University of Sheffield, Sheffield, United Kingdom
- Department of Automatic Control and Systems Engineering, University of Sheffield, Sheffield, United Kingdom
| | - Marco Viceconti
- Dipartimento di Ingegneria Industriale, Alma Mater Studiorum – University of Bologna, Bologna, Italy
- Laboratorio di Tecnologia Medica, IRCCS Istituto Ortopedico Rizzoli, Bologna, Italy
| |
Collapse
|
8
|
Mohammadi H, Pietruszczak S, Quenneville CE. Numerical analysis of hip fracture due to a sideways fall. J Mech Behav Biomed Mater 2020; 115:104283. [PMID: 33412404 DOI: 10.1016/j.jmbbm.2020.104283] [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] [Received: 08/10/2020] [Revised: 11/19/2020] [Accepted: 12/14/2020] [Indexed: 11/30/2022]
Abstract
The primary purpose of this paper is to outline a methodology for evaluating the likelihood of cortical bone fracture in the proximal femur in the event of a sideways fall. The approach includes conducting finite element (FE) analysis in which the cortical bone is treated as an anisotropic material, and the admissibility of the stress field is validated both in tension and compression regime. In assessing the onset of fracture, two methodologies are used, namely the Critical Plane approach and the Microstructure Tensor approach. The former is employed in the tension regime, while the latter governs the conditions at failure in compression. The propagation of localized damage is modeled using a constitutive law with embedded discontinuity (CLED). In this approach, the localized deformation is described by a homogenization procedure in which the average properties of cortical tissue intercepted by a macrocrack are established. The key material properties governing the conditions at failure are specified from a series of independent material tests conducted on cortical bone samples tested at different orientations relative to the loading direction. The numerical analysis deals with simulations of experiments involving the sideways fall, and the results are compared with the experimental data. This includes both the evolution of fracture pattern and the local load-displacement characteristics. The proposed approach is numerically efficient, and the results do not display a pathological mesh-dependency. Also, in contrast to the XFEM approach, the analysis does not require any extra degrees of freedom.
Collapse
Affiliation(s)
- H Mohammadi
- Department of Civil Engineering, McMaster University, Hamilton, Ontario, Canada
| | - S Pietruszczak
- Department of Civil Engineering, McMaster University, Hamilton, Ontario, Canada.
| | - C E Quenneville
- Department of Mechanical Engineering, McMaster University, Hamilton, Ontario, Canada; School of Biomedical Engineering, McMaster University, Hamilton, Ontario, Canada
| |
Collapse
|
9
|
Jazinizadeh F, Quenneville CE. 3D Analysis of the Proximal Femur Compared to 2D Analysis for Hip Fracture Risk Prediction in a Clinical Population. Ann Biomed Eng 2020; 49:1222-1232. [PMID: 33123827 DOI: 10.1007/s10439-020-02670-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2020] [Accepted: 10/20/2020] [Indexed: 01/22/2023]
Abstract
Due to the adverse impacts of hip fractures on patients' lives, it is crucial to enhance the identification of people at high risk through accessible clinical techniques. Reconstructing the 3D geometry and BMD distribution of the proximal femur could be beneficial in enhancing hip fracture risk predictions; however, it is associated with a high computational burden. It is also not clear whether it provides a better performance than 2D model analysis. Therefore, the purpose of this study was to compare the 2D and 3D model reconstruction's ability to predict hip fracture risk in a clinical population of patients. The DXA scans and CT scans of 16 cadaveric femurs were used to create training sets for the 2D and 3D model reconstruction based on statistical shape and appearance modeling. Subsequently, these methods were used to predict the risk of sustaining a hip fracture in a clinical population of 150 subjects (50 fractured, and 100 non-fractured) that were monitored for five years in the Canadian Multicentre Osteoporosis Study. 3D model reconstruction was able to improve the identification of patients who sustained a hip fracture more accurately than the standard clinical practice (by 40%). Also, the predictions from the 2D statistical model didn't differ significantly from the 3D ones (p > 0.76). These results indicated that to enhance hip fracture risk prediction in clinical practice implementing 2D statistical modeling has comparable performance with lower associated computational load.
Collapse
Affiliation(s)
- Fatemeh Jazinizadeh
- Department of Mechanical Engineering, McMaster University, ABB-C308, 1280 Main St. West, Hamilton, ON, L8S 4L8, Canada
| | - Cheryl E Quenneville
- Department of Mechanical Engineering, McMaster University, ABB-C308, 1280 Main St. West, Hamilton, ON, L8S 4L8, Canada.
- School of Biomedical Engineering, McMaster University, Hamilton, ON, Canada.
| |
Collapse
|
10
|
Jazinizadeh F, Quenneville CE. Enhancing hip fracture risk prediction by statistical modeling and texture analysis on DXA images. Med Eng Phys 2020; 78:14-20. [PMID: 32057626 DOI: 10.1016/j.medengphy.2020.01.015] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2019] [Revised: 01/14/2020] [Accepted: 01/26/2020] [Indexed: 01/09/2023]
Abstract
Each year in the US more than 300,000 older adults suffer from hip fractures. While protective measures exist, identification of those at greatest risk by DXA scanning has proved inadequate. This study proposed a new technique to enhance hip fracture risk prediction by accounting for many contributing factors to the strength of the proximal femur. Twenty-two isolated cadaveric femurs were DXA scanned, 16 of which had been mechanically tested to failure. A function consisting of the calculated modes from the statistical shape and appearance modeling (to consider the shape and BMD distribution), homogeneity index (representing trabecular quality), BMD, age and sex of the donor was created in a training set and used to predict the fracture load in a test group. To classify patients as "high risk" or "low risk", fracture load thresholds were investigated. Hip fracture load estimation was significantly enhanced using the new technique in comparison to using t-score or BMD alone (average R² of 0.68, 0.32, and 0.50, respectively) (P < 0.05). Using a fracture cut-off of 3400 N correctly predicted risk in 94% of specimens, a substantial improvement over t-score classification (38%). Ultimately, by identifying patients at high risk more accurately, devastating hip fractures can be prevented through applying protective measures.
Collapse
Affiliation(s)
- Fatemeh Jazinizadeh
- Department of Mechanical Engineering, McMaster University, 1280 Main St. West, Hamilton, Ontario L8S 4L8, Canada
| | - Cheryl E Quenneville
- Department of Mechanical Engineering, McMaster University, 1280 Main St. West, Hamilton, Ontario L8S 4L8, Canada; School of Biomedical Engineering, McMaster University, Hamilton, Ontario, Canada.
| |
Collapse
|