1
|
Asbai-Ghoudan R, Nasello G, Pérez MÁ, Verbruggen SW, Ruiz de Galarreta S, Rodriguez-Florez N. In silico assessment of the bone regeneration potential of complex porous scaffolds. Comput Biol Med 2023; 165:107381. [PMID: 37611419 DOI: 10.1016/j.compbiomed.2023.107381] [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: 04/28/2023] [Revised: 07/21/2023] [Accepted: 08/14/2023] [Indexed: 08/25/2023]
Abstract
Mechanical environment plays a crucial role in regulating bone regeneration in bone defects. Assessing the mechanobiological behavior of patient-specific orthopedic scaffolds in-silico could help guide optimal scaffold designs, as well as intra- and post-operative strategies to enhance bone regeneration and improve implant longevity. Additively manufactured porous scaffolds, and specifically triply periodic minimal surfaces (TPMS), have shown promising structural properties to act as bone substitutes, yet their ability to induce mechanobiologially-driven bone regeneration has not been elucidated. The aim of this study is to i) explore the bone regeneration potential of TPMS scaffolds made of different stiffness biocompatible materials, to ii) analyze the influence of pre-seeding the scaffolds and increasing the post-operative resting period, and to iii) assess the influence of patient-specific parameters, such as age and mechanosensitivity, on outcomes. To perform this study, an in silico model of a goat tibia is used. The bone ingrowth within the scaffold pores was simulated with a mechano-driven model of bone regeneration. Results showed that the scaffold's architectural properties affect cellular diffusion and strain distribution, resulting in variations in the regenerated bone volume and distribution. The softer material improved the bone ingrowth. An initial resting period improved the bone ingrowth but not enough to reach the scaffold's core. However, this was achieved with the implantation of a pre-seeded scaffold. Physiological parameters like age and health of the patient also influence the bone regeneration outcome, though to a lesser extent than the scaffold design. This analysis demonstrates the importance of the scaffold's geometry and its material, and highlights the potential of using mechanobiological patient-specific models in the design process for bone substitutes.
Collapse
Affiliation(s)
- Reduan Asbai-Ghoudan
- Department of Mechanical Engineering and Materials, Universidad de Navarra, TECNUN Escuela de Ingenieros, Paseo Manuel de Lardizabal, 13, 20018, San Sebastian, Spain.
| | - Gabriele Nasello
- Prometheus Division of Skeletal Tissue Engineering, KU Leuven, O&N1, Herestraat 49, PB 813, 3000, Leuven, Belgium
| | - María Ángeles Pérez
- Multiscale in Mechanical and Biological Engineering, Instituto de Investigación en Ingeniería de Aragón (I3A), Instituto de Investigación Sanitaria Aragón (IIS Aragón), University of Zaragoza, 50018, Zaragoza, Spain
| | - Stefaan W Verbruggen
- Centre for Bioengineering, School of Engineering and Materials Science, Queen Mary University of London, London, E1 4NS, UK; Department of Mechanical Engineering and INSIGNEO Institute for in Silico Medicine, University of Sheffield, Sheffield, S1 3JD, UK
| | - Sergio Ruiz de Galarreta
- Department of Mechanical Engineering and Materials, Universidad de Navarra, TECNUN Escuela de Ingenieros, Paseo Manuel de Lardizabal, 13, 20018, San Sebastian, Spain
| | - Naiara Rodriguez-Florez
- Department of Mechanical Engineering and Materials, Universidad de Navarra, TECNUN Escuela de Ingenieros, Paseo Manuel de Lardizabal, 13, 20018, San Sebastian, Spain; IKERBASQUE, Basque Foundation for Science, Plaza Euskadi 5, 48009, Bilbao, Spain
| |
Collapse
|
2
|
Henyš P, Vořechovský M, Stebel J, Kuchař M, Exner P. From computed tomography to finite element space: A unified bone material mapping strategy. Clin Biomech (Bristol, Avon) 2022; 97:105704. [PMID: 35849946 DOI: 10.1016/j.clinbiomech.2022.105704] [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: 11/19/2021] [Revised: 03/23/2022] [Accepted: 06/07/2022] [Indexed: 02/07/2023]
Abstract
BACKGROUND The spatially varying mechanical properties in finite element models of bone are most often derived from bone density data obtained via quantitative computed tomography. The key step is to accurately and efficiently map the density given in voxels to the finite element mesh. METHODS The density projection is first formulated in least-squares terms and then discretized using a continuous and discontinuous variant of the finite element method. Both discretization variants are compared with the nodal and element approaches known from the literature. FINDINGS In terms of accuracy in the L2 norm, energy distance and efficiency, the discontinuous zero-order variant appears to be the most advantageous. The proposed variant sufficiently preserves the spectrum of density at the edges, while keeping computational cost low. INTERPRETATION The continuous finite element method is analogous to the nodal formulation in the literature, while the discontinuous finite element method is analogous to the element formulation. The two variants differ in terms of implementation, computational cost and ability to preserve the density spectrum. These differences cannot be described and measured by known indirect methods from the literature.
Collapse
Affiliation(s)
- Petr Henyš
- Institute of New Technologies and Applied Informatics, Faculty of Mechatronics, Informatics and Interdisciplinary Studies, Technical University of Liberec, Studentská 1402/2, 46117 Liberec, Czech Republic
| | - Miroslav Vořechovský
- Institute of Structural Mechanics, Faculty of Civil Engineering, Brno University of Technology, Veveří 331/95, 60200 Brno, Czech Republic
| | - Jan Stebel
- Institute of New Technologies and Applied Informatics, Faculty of Mechatronics, Informatics and Interdisciplinary Studies, Technical University of Liberec, Studentská 1402/2, 46117 Liberec, Czech Republic
| | - Michal Kuchař
- Department of Anatomy, Faculty of Medicine in Hradec Králové, Charles University, Šimkova 870, 50003, Hradec Králové, Czech Republic
| | - Pavel Exner
- Institute of New Technologies and Applied Informatics, Faculty of Mechatronics, Informatics and Interdisciplinary Studies, Technical University of Liberec, Studentská 1402/2, 46117 Liberec, Czech Republic.
| |
Collapse
|
3
|
Baines AJ, Babazadeh-Naseri A, Dunbar NJ, Lewis VO, Fregly BJ. Bilateral asymmetry of bone density adjacent to pelvic sarcomas: A retrospective study using computed tomography. J Orthop Res 2022; 40:644-653. [PMID: 33914952 DOI: 10.1002/jor.25067] [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: 12/23/2020] [Revised: 04/12/2021] [Accepted: 04/26/2021] [Indexed: 02/04/2023]
Abstract
Limb-salvaging hemipelvectomy surgeries involving allograft or custom prosthesis reconstruction require high quality remaining pelvic bone for adequate device fixation. Modeling studies of custom pelvis prosthesis designs typically mirror contralateral pelvic bone material properties to the ipsilateral pelvis. However, the extent of bone material property and geometric symmetry, and thus the appropriateness of mirroring, remains unknown and should be considered when designing or analyzing the performance of pelvic prostheses. This study investigates preoperative differences between ipsilateral and contralateral pelvic bone for patients with a pelvic sarcoma. Computed tomography (CT) data were obtained retrospectively from eight patients with a pelvic sarcoma. Subject-specific computational models of the pelvic bones were constructed from the CT data. Bilateral asymmetry of bone material properties and cross-sectional areas between the ipsilateral and contralateral hemipelvis were quantified at points adjacent to the pelvic sarcoma. Large bilateral asymmetry (>20%) in trabecular but not cortical bone density was observed within 20 mm of the tumor location. Differences in trabecular bone density typically declined with increased distance from the tumor. The greatest bilateral difference in cross-sectional area occurred within 10 mm of the tumor boundary for three patients and within 40 mm from the tumor site for four patients. Our results suggest that pelvic sarcomas can cause significant bilateral asymmetries in trabecular bone density for patients with a pelvic sarcoma. These differences should be taken into account when designing custom implants for this patient population.
Collapse
Affiliation(s)
- Andrew J Baines
- Department of Mechanical Engineering, Rice University, Houston, Texas, USA
| | | | - Nicholas J Dunbar
- Department of Mechanical Engineering, Rice University, Houston, Texas, USA
| | - Valerae O Lewis
- Department of Orthopaedic Oncology, MD Anderson Cancer Center, Houston, Texas, USA
| | - Benjamin J Fregly
- Department of Mechanical Engineering, Rice University, Houston, Texas, USA
| |
Collapse
|
4
|
Babazadeh Naseri A, Dunbar NJ, Baines AJ, Akin JE, Higgs Iii CF, Fregly BJ. Heterogeneous material mapping methods for patient-specific finite element models of pelvic trabecular bone: A convergence study. Med Eng Phys 2021; 96:1-12. [PMID: 34565547 DOI: 10.1016/j.medengphy.2021.07.012] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2020] [Revised: 07/18/2021] [Accepted: 07/20/2021] [Indexed: 12/17/2022]
Abstract
Patient-specific finite element (FE) models of bone require the assignment of heterogeneous material properties extracted from the subject's computed tomography (CT) images. Though node-based (NB) and element-based (EB) material mapping methods (MMMs) have been proposed, the sensitivity and convergence of FE models to MMM for varying mesh sizes are not well understood. In this work, CT-derived and synthetic bone material data were used to evaluate the effect of MMM on results from FE analyses. Pelvic trabecular bone data was extracted from CT images of six subjects, while synthetic data were created to resemble trabecular bone properties. The numerical convergence of FE bone models using different MMMs were evaluated for strain energy, von-Mises stress, and strain. NB and EB MMMs both demonstrated good convergence regarding total strain energy, with the EB method having a slight edge over the NB. However, at the local level (e.g., maximum stress and strain), FE results were sensitive to the field type, MMM, and the FE mesh size. The EB method exhibited superior performance in finer meshes relative to the voxel size. The NB method converged better than did the EB method for coarser meshes. These findings may lead to higher-fidelity patient-specific FE bone models.
Collapse
Affiliation(s)
| | - Nicholas J Dunbar
- Department of Mechanical Engineering, Rice University, Houston, TX 77005, USA
| | - Andrew J Baines
- Department of Mechanical Engineering, Rice University, Houston, TX 77005, USA
| | - John E Akin
- Department of Mechanical Engineering, Rice University, Houston, TX 77005, USA
| | - C Fred Higgs Iii
- Department of Mechanical Engineering, Rice University, Houston, TX 77005, USA
| | - Benjamin J Fregly
- Department of Mechanical Engineering, Rice University, Houston, TX 77005, USA.
| |
Collapse
|
5
|
CHEN WENMING, CAI YIHENG, YU YUE, GENG XIANG, MA XIN. OPTIMAL MESH CRITERIA IN FINITE ELEMENT MODELING OF HUMAN FOOT: THE DEPENDENCE FOR MULTIPLE MODEL OUTPUTS ON MESH DENSITY AND LOADING BOUNDARY CONDITIONS. J MECH MED BIOL 2021. [DOI: 10.1142/s0219519421400340] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The use of finite element models has gained popularity in the field of foot and footwear biomechanics to predict the stress–strain distribution and the treatment effectiveness of therapeutic insoles for pathological foot conditions. However, a comprehensive evaluation of mesh quality is often ignored, meanwhile no golden standard exists for the mesh density and selection of element size at an acceptable accuracy. Here, we make a convergence test and established anatomically-realistic foot models at different mesh densities. The study compared the discrepancy in output variables to the changes of element type and mesh density under barefoot and footwear conditions with compressive and shear loads, which are commonly encountered in foot and footwear biomechanics simulations. For a range of loading conditions simulated in 125 finite element models, the peak plantar pressure consistently converged with optimal mesh size determined at 2.5[Formula: see text]mm. The convergence variable of principal strains and stress tensors, however, varies significantly. The max von-Mises stress showed strong sensitive behavior to the changes of the mesh density. The pattern for contact pressure distribution became less accurate when the element sizes increase to 6.0[Formula: see text]mm; in particular, the locations of the pressure peak do not show remarkable changes, but the size of the area of contact still changes. The current study could offer a general guideline when generating a reasonable accurate finite element models for the analysis of plantar pressure distributions and stress/strain states employed for foot and footwear biomechanics evaluations.
Collapse
Affiliation(s)
- WEN-MING CHEN
- Academy for Engineering & Technology, Fudan University, 220 Handan Road, Shanghai, P. R. China
| | - YI-HENG CAI
- Department of Biomedical Engineering, University of Shanghai for Science and Technology, Shanghai, P. R. China
| | - YUE YU
- Department of medical imaging, Suzhou Institute of Biomedical Engineering and Technology, Chinese Academy of Sciences, Suzhou, P. R. China
| | - XIANG GENG
- Department of Orthopaedics, Huashan Hospital, Fudan University, 12 Middle Wulumuqi Road, Shanghai, P. R. China
| | - XIN MA
- Department of Orthopaedics, Huashan Hospital, Fudan University, 12 Middle Wulumuqi Road, Shanghai, P. R. China
| |
Collapse
|
6
|
Nasello G, Vautrin A, Pitocchi J, Wesseling M, Kuiper JH, Pérez MÁ, García-Aznar JM. Mechano-driven regeneration predicts response variations in large animal model based on scaffold implantation site and individual mechano-sensitivity. Bone 2021; 144:115769. [PMID: 33276152 DOI: 10.1016/j.bone.2020.115769] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/24/2020] [Revised: 11/02/2020] [Accepted: 11/20/2020] [Indexed: 02/02/2023]
Abstract
It is well founded that the mechanical environment may regulate bone regeneration in orthopedic applications. The purpose of this study is to investigate the mechanical contributions of the scaffold and the host to bone regeneration, in terms of subject specificity, implantation site and sensitivity to the mechanical environment. Using a computational approach to model mechano-driven regeneration, bone ingrowth in porous titanium scaffolds was simulated in the distal femur and proximal tibia of three goats and compared to experimental results. The results showed that bone ingrowth shifted from a homogeneous distribution pattern, when scaffolds were in contact with trabecular bone (max local ingrowth 12.47%), to a localized bone ingrowth when scaffolds were implanted in a diaphyseal location (max local ingrowth 20.64%). The bone formation dynamics revealed an apposition rate of 0.37±0.28%/day in the first three weeks after implantation, followed by limited increase in bone ingrowth until the end of the experiment (12 weeks). According to in vivo data, we identified one animal whose sensitivity to mechanical stimulation was higher than the other two. Moreover, we found that the stimulus initiating bone formation was consistently higher in the femur than in the tibia for all the individuals. Overall, the dependence of the osteogenic response on the host biomechanics means that, from a mechanical perspective, the regenerative potential depends on both the scaffold and the host environment. Therefore, this work provides insights on how the mechanical conditions of both the recipient and the scaffold contribute to meet patient and location-specific characteristics.
Collapse
Affiliation(s)
- Gabriele Nasello
- Multiscale in Mechanical and Biological Engineering, Instituto de Investigación en Ingeniería de Aragón (I3A), Instituto de Investigación Sanitaria Aragón (IIS Aragón), University of Zaragoza, Zaragoza, Spain; Biomechanics Section, KU Leuven, Leuven, Belgium
| | - Antoine Vautrin
- Ecole Nationale d'Ingénieurs de Metz, University of Lorraine, Metz, France
| | - Jonathan Pitocchi
- Multiscale in Mechanical and Biological Engineering, Instituto de Investigación en Ingeniería de Aragón (I3A), Instituto de Investigación Sanitaria Aragón (IIS Aragón), University of Zaragoza, Zaragoza, Spain; Biomechanics Section, KU Leuven, Leuven, Belgium; Materialise NV, Leuven, Belgium
| | | | - Jan Herman Kuiper
- Institute for Science and Technology in Medicine, Keele University, Keele, UK; The Robert Jones and Agnes Hunt Orthopaedic Hospital, NHS Foundation Trust, Oswestry, UK
| | - María Ángeles Pérez
- Multiscale in Mechanical and Biological Engineering, Instituto de Investigación en Ingeniería de Aragón (I3A), Instituto de Investigación Sanitaria Aragón (IIS Aragón), University of Zaragoza, Zaragoza, Spain
| | - José Manuel García-Aznar
- Multiscale in Mechanical and Biological Engineering, Instituto de Investigación en Ingeniería de Aragón (I3A), Instituto de Investigación Sanitaria Aragón (IIS Aragón), University of Zaragoza, Zaragoza, Spain.
| |
Collapse
|
7
|
QCT-FE modeling of the proximal tibia: Effect of mapping strategy on convergence time and model accuracy. Med Eng Phys 2021; 88:41-46. [PMID: 33485512 DOI: 10.1016/j.medengphy.2020.12.006] [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] [Received: 12/10/2019] [Revised: 12/04/2020] [Accepted: 12/22/2020] [Indexed: 11/21/2022]
Abstract
Quantitative computed tomography (QCT) based finite element (FE) modeling, referred to as QCT-FE, has seen rapid growth and application for modeling bone mechanics. With this approach, varying bone material properties are set via experimentally-derived density-modulus equations. One challenge though associated with QCT-FE is to identify the appropriate mapping strategy for assigning elastic moduli to elements. The goal of this study was to evaluate different QCT-FE mapping strategies to identify the optimum approach with fastest convergence rate and highest accuracy. Four proximal tibial medial compartments were imaged using QCT and experimentally tested to characterize proximal tibial subchondral bone stiffness at four surface points, resulting in a total of 16 indentation measures. Three material mapping methods were analyzed: (1) constant-E where an average elastic modulus was assigned to each element; (2) node-based where the material properties were first mapped on nodes then interpolated to Gaussian integration points; and (3) element-based in which the material properties were directly assigned to Gaussian integration points. Different element sizes were assessed with edge-lengths ranging from 0.9 to 3 mm. Results indicated that all converged models showed similar coefficient-of-determination (R2) and normalized root-mean-square errors (RMSE%). Though, the constant-E and node-based methods converged with the element edge-length of 1.5 mm (prediction error of 4.8% and 2.5%, respectively) whereas the element-based method converged with a larger element having an edge-length 2.5 mm (error = 4.9%). In conclusion, the element-based method, with a larger element size, resulted in similar predictive accuracy, faster convergence and shorter run-times relative to the constant-E and node-based approaches. As such, we recommend the element-based method for future subject-specific QCT-FE modeling.
Collapse
|
8
|
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: 9] [Impact Index Per Article: 1.8] [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
|
9
|
Knowles NK, Ip K, Ferreira LM. The Effect of Material Heterogeneity, Element Type, and Down-Sampling on Trabecular Stiffness in Micro Finite Element Models. Ann Biomed Eng 2018; 47:615-623. [PMID: 30362084 DOI: 10.1007/s10439-018-02152-6] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2018] [Accepted: 10/05/2018] [Indexed: 11/30/2022]
Abstract
Preclinical and clinical bone strength predictions can be elucidated by understanding bone mechanics at a variety of hierarchical levels. As such, down-sampled micro-CT images are often used to make comparisons across image resolutions or used to reduce computational resources in micro finite element models (µFEMs). Therefore, the objectives of this study were to compare trabecular apparent modulus among (i) hexahedral and tetrahedral µFEMs, (ii) µFEMs generated from 32, 64, and 64 µm down-sampled from 32 µm µCT scans, and (iii) µFEMs with homogeneous and heterogeneous tissue moduli. Trabecular µFEMs were generated from scans at the three spatial resolutions taken from the glenoid vault of 14 cadaveric specimens. Simulated unconstrained compression was performed and used to calculate and compare the apparent modulus of each µFEM. It was found that models derived from high-resolution images that account for material heterogeneity are nearly equivalent whether hexahedral or tetrahedral elements are used. However, translation of stiffness from down-sampled scans are not equivalent to scans performed at the down-sampled resolution, or that account for trabecular material heterogeneity. Material heterogeneity is most representative of in vivo trabecular bone and to accurately model trabecular mechanical properties, material heterogeneity should be considered in future µFEM development.
Collapse
Affiliation(s)
- Nikolas K Knowles
- School of Biomedical Engineering, The University of Western Ontario, London, ON, Canada.,Surgical Mechatronics Laboratory, Roth
- McFarlane Hand and Upper Limb Centre, 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
| | - Kenneth Ip
- School of Biomedical Engineering, The University of Western Ontario, London, ON, Canada.,Surgical Mechatronics Laboratory, Roth
- McFarlane Hand and Upper Limb Centre, St. Josephs Health Care, 268 Grosvenor St., London, ON, Canada
| | - Louis M Ferreira
- Department of Mechanical and Materials Engineering, The University of Western Ontario, London, ON, Canada. .,Surgical Mechatronics Laboratory, Roth
- McFarlane Hand and Upper Limb Centre, 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
|