1
|
Investigation on the Differences in the Failure Processes of the Cortical Bone under Different Loading Conditions. Appl Bionics Biomech 2022; 2022:3406984. [DOI: 10.1155/2022/3406984] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2022] [Revised: 09/18/2022] [Accepted: 09/30/2022] [Indexed: 11/18/2022] Open
Abstract
Cortical bone is a transversely isotropic material, and the mechanical properties may be related to the loading direction on the osteon. Therefore, analyzing the differences in the failure processes of cortical bone under different loading conditions is necessary to explore the measures for reducing the incidence of fracture. In this study, to investigate the effects of different loading directions on the fracture performance in the cortical bone, a numerical method that could simultaneously simulate the failure processes in the cortical bone structure under compression and bending loads was established based on continuum damage mechanics theory. The prediction accuracy and feasibility of the numerical method were first verified by comparing with the corresponding experimental results. Then, the differences in the failure process and fracture performance of the same cortical bone structure under compression and bending loads were investigated. The simulation results indicated that for the same structure, the slip-open failure mode appeared under compression load, and the crack propagated along a certain angle to the loading direction; the tension-open failure mode appeared under bending load, and the crack propagated along the direction perpendicular to the loading direction. Meanwhile, the fracture load was greater and the fracture time was later in the compression than in the bending condition. These phenomena stated that discrepant failure processes and fracture patterns occurred in the same cortical bone structure under different loading conditions. The main reason may be related to the tension–compression asymmetry and transversely isotropic characteristics in the cortical bone material. The fracture simulations in the cortical bone under different loading conditions could improve the prediction accuracy in bone biomechanics and provide the prevention method for cortical bone damage and fracture.
Collapse
|
2
|
Sabet FA, Koric S, Idkaidek A, Jasiuk I. High-Performance Computing Comparison of Implicit and Explicit Nonlinear Finite Element Simulations of Trabecular Bone. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2021; 200:105870. [PMID: 33280935 DOI: 10.1016/j.cmpb.2020.105870] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/08/2020] [Accepted: 11/18/2020] [Indexed: 06/12/2023]
Abstract
BACKGROUND AND OBJECTIVE Finite element models built from micro-computed tomography scans have become a powerful tool to investigate the mechanical properties of trabecular bone. There are two types of solving algorithms in the finite element method: implicit and explicit. Both of these methods have been utilized to study the trabecular bone. However, an investigation comparing the results obtained using the implicit and explicit solvers is lacking. Thus, in this paper, we contrast implicit and explicit procedures by analyzing trabecular bone samples as a case study. METHODS Micro-computed tomography-based finite element analysis of trabecular bone under a direct quasi-static compression was done using implicit and explicit methods. The differences in the predictions of mechanical properties and computational time of the two methods were studied using high-performance computing. RESULTS Our findings indicate that the results using implicit and explicit solvers are well comparable, given that similar problem set up is carefully utilized. Also, the parallel scalability of the two methods was similar, while the explicit solver performed about five times faster than the implicit method. Along with faster performance, the explicit method utilized significantly less memory for the analysis, which shows another benefit of using an explicit solver for this case study. CONCLUSIONS The comparison of the implicit and explicit methods for the simulation of trabecular bone samples should be highly valuable to the bone modeling community and researchers studying complex cellular and architectured materials.
Collapse
Affiliation(s)
- Fereshteh A Sabet
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Seid Koric
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA; National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Ashraf Idkaidek
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Iwona Jasiuk
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA.
| |
Collapse
|
3
|
Sadoughi S, Vom Scheidt A, Nawathe S, Zhu S, Moini A, Keaveny TM. Effect of variations in tissue-level ductility on human vertebral strength. Bone 2020; 137:115445. [PMID: 32454256 DOI: 10.1016/j.bone.2020.115445] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/13/2019] [Revised: 05/17/2020] [Accepted: 05/19/2020] [Indexed: 11/28/2022]
Abstract
Although the ductility of bone tissue is a unique element of bone quality and varies with age and across the population, the extent to which and mechanisms by which typical population-variations in tissue-level ductility can alter whole-bone strength remains unclear. To provide insight, we conducted a finite element analysis parameter study of whole-vertebral (monotonic) compressive strength on six human L1 vertebrae. Each model was generated from micro-CT scans, capturing the trabecular micro-architecture in detail, and included a non-linear constitutive model for the bone tissue that allowed for plastic yielding, different strengths in tension and compression, large deformations, and, uniquely, localized damage once a specified limit in tissue-level ultimate strain was exceeded. Those strain limits were based on reported (mean ± SD) values from cadaver experiments (8.8 ± 3.7% strain for trabecular tissue and 2.2 ± 0.9% for cortical tissue). In the parameter study, the strain limits were varied by ±1 SD from their mean values, for a combination of nine analyses per specimen; bounding values of zero and unlimited post-yield strain were also modeled. The main outcomes from the finite element analysis were the vertebral compressive strength and the amount of failed (yielded or damaged) tissue at the overall structure-level failure. Compared to a reference case of using the mean values of ultimate strain, we found that varying both trabecular and cortical tissue ultimate strains by ±1 SD changed the computed vertebral strength by (mean ± SD) ±6.9 ± 1.1% on average. Mechanistically, that modest effect arose because the proportion of yielded tissue (without damage) was 0.9 ± 0.3% of all the bone tissue across the nine cases and the proportion of damaged tissue (i.e. tissue exceeding the prescribed tissue-level ultimate strain) was 0.2 ± 0.1%. If the types of variations in tissue-level ductility investigated here accurately represent real typical variations in the population, the consistency of our results across specimens and the modest effect size together suggest that typical variations in tissue-level ductility only have a modest impact on vertebral compressive strength, in large part because so few trabeculae are damaged at the load capacity of the bone.
Collapse
Affiliation(s)
- Saghi Sadoughi
- Department of Mechanical Engineering, University of California, Berkeley, CA, USA
| | - Annika Vom Scheidt
- Department of Mechanical Engineering, University of California, Berkeley, CA, USA; Department of Osteology and Biomechanics, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Shashank Nawathe
- Department of Mechanical Engineering, University of California, Berkeley, CA, USA
| | - Shan Zhu
- Department of Mechanical Engineering, University of California, Berkeley, CA, USA
| | - Ariana Moini
- Department of Mechanical Engineering, University of California, Berkeley, CA, USA
| | - Tony M Keaveny
- Department of Mechanical Engineering, University of California, Berkeley, CA, USA; Department of Bioengineering, University of California, Berkeley, CA, USA.
| |
Collapse
|
4
|
Salem M, Westover L, Adeeb S, Duke K. Prediction of failure in cancellous bone using extended finite element method. Proc Inst Mech Eng H 2020; 234:988-999. [PMID: 32605523 DOI: 10.1177/0954411920936057] [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: 11/17/2022]
Abstract
The objective of our study is to develop extended finite element method models of cancellous bone specimens that are capable of accurately predicting the onset and propagation of cracks under mechanical loading. In order to do so, previously published three-point bending test results of a single trabecula were replicated using two different extended finite element method approaches, namely, elastic-plastic-fracture and elastic-fracture that considered different configurations of the elasto-plastic properties of bone from which the best approach to fit the experimental data was identified. The behavior of a single trabecula was then used in 2D extended finite element method models to quantify the strength of the trabecular tissue of the forearm along three perpendicular anatomical axes. The results revealed that the elastic-plastic-fracture model better represented the experimental data in the model of a single trabecula. Considering the 2D trabecular specimens, the elastic fracture model predicted higher strength than the elastic-plastic-fracture model and there was no difference in stiffness between the two models. In general, the specimens exhibited higher failure strain and more ductile behavior in compression than in tension. In addition, strength and stiffness were found to be higher in tension than compression on average. It can be concluded that with proper parameters, extended finite element method is capable of simulating the ductile behavior of cancellous bone. The models are able to quantify the tensile strength of trabecular tissue in the various anatomical directions reporting an increased strength in the longitudinal direction of forearm cancellous bone tissue. Extended finite element method of cancellous bone proves to be a valuable tool to predict the mechanical characteristics of cancellous bones as a function of the microstructure.
Collapse
Affiliation(s)
- Mohammad Salem
- Department of Mechanical Engineering, University of Alberta, Edmonton, AB, Canada
| | - Lindsey Westover
- Department of Mechanical Engineering, University of Alberta, Edmonton, AB, Canada
| | - Samer Adeeb
- Department of Civil and Environmental Engineering, University of Alberta, Edmonton, AB, Canada
| | - Kajsa Duke
- Department of Mechanical Engineering, University of Alberta, Edmonton, AB, Canada
| |
Collapse
|
5
|
Efficient materially nonlinear [Formula: see text]FE solver for simulations of trabecular bone failure. Biomech Model Mechanobiol 2019; 19:861-874. [PMID: 31749070 PMCID: PMC7203600 DOI: 10.1007/s10237-019-01254-x] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2019] [Accepted: 11/05/2019] [Indexed: 01/15/2023]
Abstract
An efficient solver for large-scale linear \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mu \hbox {FE}$$\end{document}μFE simulations was extended for nonlinear material behavior. The material model included damage-based tissue degradation and fracture. The new framework was applied to 20 trabecular biopsies with a mesh resolution of \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${36}\,{{\upmu }\hbox {m}}$$\end{document}36μm. Suitable material parameters were identified based on two biopsies by comparison with axial tension and compression experiments. The good parallel performance and low memory footprint of the solver were preserved. Excellent correlation of the maximum apparent stress was found between simulations and experiments (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$R^2 > 0.97$$\end{document}R2>0.97). The development of local damage regions was observable due to the nonlinear nature of the simulations. A novel elasticity limit was proposed based on the local damage information. The elasticity limit was found to be lower than the 0.2% yield point. Systematic differences in the yield behavior of biopsies under apparent compression and tension loading were observed. This indicates that damage distributions could lead to more insight into the failure mechanisms of trabecular bone.
Collapse
|
6
|
Levrero-Florencio F, Pankaj P. Using Non-linear Homogenization to Improve the Performance of Macroscopic Damage Models of Trabecular Bone. Front Physiol 2018; 9:545. [PMID: 29867581 PMCID: PMC5966630 DOI: 10.3389/fphys.2018.00545] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2018] [Accepted: 04/27/2018] [Indexed: 11/13/2022] Open
Abstract
Realistic macro-level finite element simulations of the mechanical behavior of trabecular bone, a cellular anisotropic material, require a suitable constitutive model; a model that incorporates the mechanical response of bone for complex loading scenarios and includes post-elastic phenomena, such as plasticity (permanent deformations) and damage (permanent stiffness reduction), which bone is likely to experience. Some such models have been developed by conducting homogenization-based multiscale finite element simulations on bone micro-structure. While homogenization has been fairly successful in the elastic regime and, to some extent, in modeling the macroscopic plastic response, it has remained a challenge with respect to modeling damage. This study uses a homogenization scheme to upscale the damage behavior from the tissue level (microscale) to the organ level (macroscale) and assesses the suitability of different damage constitutive laws. Ten cubic specimens were each subjected to 21 strain-controlled load cases for a small range of macroscopic post-elastic strains. Isotropic and anisotropic criteria were considered, density and fabric relationships were used in the formulation of the damage law, and a combined isotropic/anisotropic law with tension/compression asymmetry was formulated, based on the homogenized results, as a possible alternative to the currently used single scalar damage criterion. This computational study enhances the current knowledge on the macroscopic damage behavior of trabecular bone. By developing relationships of damage progression with bone's micro-architectural indices (density and fabric) the study also provides an aid for the creation of more precise macroscale continuum models, which are likely to improve clinical predictions.
Collapse
Affiliation(s)
- Francesc Levrero-Florencio
- Computational Cardiovascular Science, Department of Computer Science, University of Oxford, Oxford, United Kingdom.,Institute for Bioengineering, School of Engineering, The University of Edinburgh, Edinburgh, United Kingdom
| | - Pankaj Pankaj
- Institute for Bioengineering, School of Engineering, The University of Edinburgh, Edinburgh, United Kingdom
| |
Collapse
|
7
|
Sabet FA, Jin O, Koric S, Jasiuk I. Nonlinear micro-CT based FE modeling of trabecular bone-Sensitivity of apparent response to tissue constitutive law and bone volume fraction. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2018; 34:e2941. [PMID: 29168345 DOI: 10.1002/cnm.2941] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/10/2017] [Revised: 09/29/2017] [Accepted: 11/06/2017] [Indexed: 06/07/2023]
Abstract
In this study, the sensitivity of the apparent response of trabecular bone to different constitutive models at the tissue level was investigated using finite element (FE) modeling based on micro-computed tomography (micro-CT). Trabecular bone specimens from porcine femurs were loaded under a uniaxial compression experimentally and computationally. The apparent behaviors computed using von Mises, Drucker-Prager, and Cast Iron plasticity models were compared. Secondly, the effect of bone volume fraction was studied by changing the bone volume fraction of a trabecular bone sample while keeping the same basic architecture. Also, constitutive models' parameters of the tissue were calibrated for porcine bone, and the effects of different parameters on resulting apparent response were investigated through a parametric study. The calibrated effective tissue elastic modulus of porcine trabecular bone was 10±1.2 GPa, which is in the lower range of modulus values reported in the literature for human and bovine trabecular bones (4-23.8 GPa). It was also observed that, unlike elastic modulus, yield properties of tissue could not be uniquely calibrated by fitting an apparent response from simulations to experiments under a uniaxial compression. Our results demonstrated that using these 3 tissue constitutive models had only a slight effect on the apparent response. As expected, there was a significant change in the apparent response with varying bone volume fraction. Also, both apparent modulus and maximum stress had a linear relation with bone volume fraction.
Collapse
Affiliation(s)
- F A Sabet
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA
| | - O Jin
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA
| | - S Koric
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA
- National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA
| | - I Jasiuk
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA
| |
Collapse
|
8
|
Hou J, Zhang Y, Sun Y, Xu N, Leng Y. Prediction of Firmness and pH for "Golden Delicious" Apple Based on Elasticity Index from Modal Analysis. J Food Sci 2018; 83:661-669. [PMID: 29437233 DOI: 10.1111/1750-3841.14071] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2017] [Revised: 01/09/2018] [Accepted: 01/10/2018] [Indexed: 11/30/2022]
Abstract
An experimental modal test system was established to extract the natural frequencies of "Golden Delicious" apple, after which the elasticity index was calculated to predict the apple quality parameters based on the orthogonal polynomials method. The elasticity index in every vibration mode changed dramatically (P = 0.01) along time revolution. The multivariate regression methods were used to model the predictive relationship between the elasticity index and the apple quality parameters. The models of the apple juice pH based on support vector regression presented adequate determination coefficients of calibration set (Q2 = 0.68) and prediction set (R2 = 0.55), respectively. The models based on partial least squares regression could be used for predicting the apple firmness parameter offset gradient (Q2 = 0.76 and R2 = 0.72). It helped understanding the fruit dynamic properties of the fruit and spontaneously obtaining the fruit chemical parameters. A nondestructive and portable device was viable for fruit quality estimation by the modal test system during storage, transport, and even growth on the tree. PRACTICAL APPLICATION A nondestructive and portable device was provided for fruit quality detection during storage, transport and even growth based on experimental modal analysis. A systematic statistical analysis method about outlier detection, data set partitioning, parameter optimization, and multiple regression models were provided.
Collapse
Affiliation(s)
- Jumin Hou
- College of Food Science and Engineering, Jilin Univ., No. 5333, Xi'an Road, Changchun, Jilin, China
| | - Yuxia Zhang
- College of Food Science and Engineering, Jilin Univ., No. 5333, Xi'an Road, Changchun, Jilin, China
| | - Yonghai Sun
- College of Food Science and Engineering, Jilin Univ., No. 5333, Xi'an Road, Changchun, Jilin, China
| | - Na Xu
- College of Food Science and Engineering, Jilin Univ., No. 5333, Xi'an Road, Changchun, Jilin, China
| | - Yue Leng
- College of Food Science and Engineering, Jilin Univ., No. 5333, Xi'an Road, Changchun, Jilin, China
| |
Collapse
|
9
|
Effect of including damage at the tissue level in the nonlinear homogenisation of trabecular bone. Biomech Model Mechanobiol 2017; 16:1681-1695. [PMID: 28500359 PMCID: PMC5599493 DOI: 10.1007/s10237-017-0913-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2016] [Accepted: 04/21/2017] [Indexed: 02/05/2023]
Abstract
Being able to predict bone fracture or implant stability needs a proper constitutive model of trabecular bone at the macroscale in multiaxial, non-monotonic loading modes. Its macroscopic damage behaviour has been investigated experimentally in the past, mostly with the restriction of uniaxial cyclic loading experiments for different samples, which does not allow for the investigation of several load cases in the same sample as damage in one direction may affect the behaviour in other directions. Homogenised finite element models of whole bones have the potential to assess complicated scenarios and thus improve clinical predictions. The aim of this study is to use a homogenisation-based multiscale procedure to upscale the damage behaviour of bone from an assumed solid phase constitutive law and investigate its multiaxial behaviour for the first time. Twelve cubic specimens were each submitted to nine proportional strain histories by using a parallel code developed in-house. Evolution of post-elastic properties for trabecular bone was assessed for a small range of macroscopic plastic strains in these nine load cases. Damage evolution was found to be non-isotropic, and both damage and hardening were found to depend on the loading mode (tensile, compression or shear); both were characterised by linear laws with relatively high coefficients of determination. It is expected that the knowledge of the macroscopic behaviour of trabecular bone gained in this study will help in creating more precise continuum FE models of whole bones that improve clinical predictions.
Collapse
|
10
|
Levrero-Florencio F, Manda K, Margetts L, Pankaj P. Nonlinear homogenisation of trabecular bone: Effect of solid phase constitutive model. Proc Inst Mech Eng H 2016; 231:405-414. [PMID: 28427317 DOI: 10.1177/0954411916676220] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Micro-finite element models have been extensively employed to evaluate the elastic properties of trabecular bone and, to a limited extent, its yield behaviour. The macroscopic stiffness tensor and yield surface are of special interest since they are essential in the prediction of bone strength and stability of implants at the whole bone level. While macroscopic elastic properties are now well understood, yield and post-yield properties are not. The aim of this study is to shed some light on what the effect of the solid phase yield criterion is on the macroscopic yield of trabecular bone for samples with different microstructure. Three samples with very different density were subjected to a large set of apparent load cases (which is important since physiological loading is complex and can have multiple components in stress or strain space) with two different solid phase yield criteria: Drucker-Prager and eccentric-ellipsoid. The study found that these two criteria led to small differences in the macroscopic yield strains for most load cases except for those that were compression-dominated; in these load cases, the yield strains for the Drucker-Prager criterion were significantly higher. Higher density samples resulted in higher differences between the two criteria. This work provides a comprehensive assessment of the effect of two different solid phase yield criteria on the macroscopic yield strains of trabecular bone, for a wide range of load cases, and for samples with different morphology.
Collapse
Affiliation(s)
| | - Krishnagoud Manda
- 1 Institute for Bioengineering, School of Engineering, The University of Edinburgh, Edinburgh, UK
| | - Lee Margetts
- 2 School of Mechanical, Aerospace and Civil Engineering The University of Manchester Sackville Street Manchester, UK
| | - Pankaj Pankaj
- 1 Institute for Bioengineering, School of Engineering, The University of Edinburgh, Edinburgh, UK
| |
Collapse
|
11
|
Morphology based anisotropic finite element models of the proximal femur validated with experimental data. Med Eng Phys 2016; 38:1339-1347. [DOI: 10.1016/j.medengphy.2016.08.010] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2016] [Revised: 08/05/2016] [Accepted: 08/30/2016] [Indexed: 11/21/2022]
|
12
|
Fan R, Gong H, Zhang R, Gao J, Jia Z, Hu Y. Quantification of Age-Related Tissue-Level Failure Strains of Rat Femoral Cortical Bones Using an Approach Combining Macrocompressive Test and Microfinite Element Analysis. J Biomech Eng 2016; 138:041006. [DOI: 10.1115/1.4032798] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2015] [Indexed: 12/15/2022]
Abstract
Bone mechanical properties vary with age; meanwhile, a close relationship exists among bone mechanical properties at different levels. Therefore, conducting multilevel analyses for bone structures with different ages are necessary to elucidate the effects of aging on bone mechanical properties at different levels. In this study, an approach that combined microfinite element (micro-FE) analysis and macrocompressive test was established to simulate the failure of male rat femoral cortical bone. Micro-FE analyses were primarily performed for rat cortical bones with different ages to simulate their failure processes under compressive load. Tissue-level failure strains in tension and compression of these cortical bones were then back-calculated by fitting the experimental stress–strain curves. Thus, tissue-level failure strains of rat femoral cortical bones with different ages were quantified. The tissue-level failure strain exhibited a biphasic behavior with age: in the period of skeletal maturity (1–7 months of age), the failure strain gradually increased; when the rat exceeded 7 months of age, the failure strain sharply decreased. In the period of skeletal maturity, both the macro- and tissue-levels mechanical properties showed a large promotion. In the period of skeletal aging (9–15 months of age), the tissue-level mechanical properties sharply deteriorated; however, the macromechanical properties only slightly deteriorated. The age-related changes in tissue-level failure strain were revealed through the analysis of male rat femoral cortical bones with different ages, which provided a theoretical basis to understand the relationship between rat cortical bone mechanical properties at macro- and tissue-levels and decrease of bone strength with age.
Collapse
Affiliation(s)
- Ruoxun Fan
- State Key Laboratory of Automotive Simulation and Control, Jilin University, Changchun 130025, China
- Department of Engineering Mechanics, Jilin University, Nanling Campus, Changchun 130025, China e-mail:
| | - He Gong
- Professor State Key Laboratory of Automotive Simulation and Control, Jilin University, Changchun 130025, China
- Department of Engineering Mechanics, Jilin University, Nanling Campus, Changchun 130025, China e-mail:
| | - Rui Zhang
- Department of Engineering Mechanics, Jilin University, Nanling Campus, Changchun 130025, China
- Key Laboratory for Biomechanics and Mechanobiology of Ministry of Education, School of Biological Science and Medical Engineering, Beihang University, Beijing 10000, China e-mail:
| | - Jiazi Gao
- Department of Engineering Mechanics, Jilin University, Nanling Campus, Changchun 130025, China e-mail:
| | - Zhengbin Jia
- Department of Engineering Mechanics, Jilin University, Nanling Campus, Changchun 130025, China e-mail:
| | - Yanjuan Hu
- School of Mechatronic Engineering, Changchun University of Technology, Changchun 130025, China e-mail:
| |
Collapse
|
13
|
|
14
|
Sanyal A, Scheffelin J, Keaveny TM. The quartic piecewise-linear criterion for the multiaxial yield behavior of human trabecular bone. J Biomech Eng 2015; 137:1937621. [PMID: 25401413 DOI: 10.1115/1.4029109] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2014] [Accepted: 11/17/2014] [Indexed: 11/08/2022]
Abstract
Prior multiaxial strength studies on trabecular bone have either not addressed large variations in bone volume fraction and microarchitecture, or have not addressed the full range of multiaxial stress states. Addressing these limitations, we utilized micro-computed tomography (lCT) based nonlinear finite element analysis to investigate the complete 3D multiaxial failure behavior of ten specimens (5mm cube) of human trabecular bone, taken from three anatomic sites and spanning a wide range of bone volume fraction (0.09–0.36),mechanical anisotropy (range of E3/E1¼3.0–12.0), and microarchitecture. We found that most of the observed variation in multiaxial strength behavior could be accounted for by normalizing the multiaxial strength by specimen-specific values of uniaxial strength (tension,compression in the longitudinal and transverse directions). Scatter between specimens was reduced further when the normalized multiaxial strength was described in strain space.The resulting multiaxial failure envelope in this normalized-strain space had a rectangular boxlike shape for normal–normal loading and either a rhomboidal box like shape or a triangular shape for normal-shear loading, depending on the loading direction. The finite element data were well described by a single quartic yield criterion in the 6D normalized strain space combined with a piecewise linear yield criterion in two planes for normalshear loading (mean error SD: 4.660.8% for the finite element data versus the criterion).This multiaxial yield criterion in normalized-strain space can be used to describe the complete 3D multiaxial failure behavior of human trabecular bone across a wide range of bone volume fraction, mechanical anisotropy, and microarchitecture.
Collapse
|
15
|
Panyasantisuk J, Pahr DH, Zysset PK. Effect of boundary conditions on yield properties of human femoral trabecular bone. Biomech Model Mechanobiol 2015; 15:1043-53. [PMID: 26517986 DOI: 10.1007/s10237-015-0741-6] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2015] [Accepted: 10/20/2015] [Indexed: 11/29/2022]
Abstract
Trabecular bone plays an important mechanical role in bone fractures and implant stability. Homogenized nonlinear finite element (FE) analysis of whole bones can deliver improved fracture risk and implant loosening assessment. Such simulations require the knowledge of mechanical properties such as an appropriate yield behavior and criterion for trabecular bone. Identification of a complete yield surface is extremely difficult experimentally but can be achieved in silico by using micro-FE analysis on cubical trabecular volume elements. Nevertheless, the influence of the boundary conditions (BCs), which are applied to such volume elements, on the obtained yield properties remains unknown. Therefore, this study compared homogenized yield properties along 17 load cases of 126 human femoral trabecular cubic specimens computed with classical kinematic uniform BCs (KUBCs) and a new set of mixed uniform BCs, namely periodicity-compatible mixed uniform BCs (PMUBCs). In stress space, PMUBCs lead to 7-72 % lower yield stresses compared to KUBCs. The yield surfaces obtained with both KUBCs and PMUBCs demonstrate a pressure-sensitive ellipsoidal shape. A volume fraction and fabric-based quadric yield function successfully fitted the yield surfaces of both BCs with a correlation coefficient [Formula: see text]. As expected, yield strains show only a weak dependency on bone volume fraction and fabric. The role of the two BCs in homogenized FE analysis of whole bones will need to be investigated and validated with experimental results at the whole bone level in future studies.
Collapse
Affiliation(s)
- J Panyasantisuk
- Institute for Surgical Technology and Biomechanics, University of Bern, Stauffacherstr. 78, 3014, Bern, Switzerland.
| | - D H Pahr
- Institute of Lightweight Design and Structural Biomechanics, Vienna University of Technology, Getreidemarkt 9, 1060, Vienna, Austria
| | - P K Zysset
- Institute for Surgical Technology and Biomechanics, University of Bern, Stauffacherstr. 78, 3014, Bern, Switzerland
| |
Collapse
|
16
|
Goff MG, Lambers FM, Sorna RM, Keaveny TM, Hernandez CJ. Finite element models predict the location of microdamage in cancellous bone following uniaxial loading. J Biomech 2015; 48:4142-4148. [PMID: 26522622 DOI: 10.1016/j.jbiomech.2015.10.023] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2015] [Revised: 10/10/2015] [Accepted: 10/18/2015] [Indexed: 10/22/2022]
Abstract
High-resolution finite element models derived from micro-computed tomography images are often used to study the effects of trabecular microarchitecture and loading mode on tissue stress, but the degree to which existing finite element methods correctly predict the location of tissue failure is not well characterized. In the current study, we determined the relationship between the location of highly strained tissue, as determined from high-resolution finite element models, and the location of tissue microdamage, as determined from three-dimensional fluoroscopy imaging, which was performed after the microdamage was generated in-vitro by mechanical testing. Fourteen specimens of human vertebral cancellous bone were assessed (8 male donors, 2 female donors, 47-78 years of age). Regions of stained microdamage, were 50-75% more likely to form in highly strained tissue (principal strains exceeding 0.4%) than elsewhere, and generally the locations of the regions of microdamage were significantly correlated (p<0.05) with the locations of highly strained tissue. This spatial correlation was stronger for the largest regions of microdamage (≥1,000,000μm(3) in volume); 87% of large regions of microdamage were located near highly strained tissue. Together, these findings demonstrate that there is a strong correlation between regions of microdamage and regions of high strain in human cancellous bone, particularly for the biomechanically more important large instances of microdamage.
Collapse
Affiliation(s)
- M G Goff
- Department of Biomedical Engineering, Cornell University, Ithaca, NY, USA; Sibley School of Mechanical and Aerospace Engineering, Cornell University, Ithaca, NY, USA
| | - F M Lambers
- Sibley School of Mechanical and Aerospace Engineering, Cornell University, Ithaca, NY, USA
| | - R M Sorna
- Sibley School of Mechanical and Aerospace Engineering, Cornell University, Ithaca, NY, USA
| | - T M Keaveny
- Department of Mechanical Engineering, University of California, Berkeley, CA, USA
| | - C J Hernandez
- Department of Biomedical Engineering, Cornell University, Ithaca, NY, USA; Sibley School of Mechanical and Aerospace Engineering, Cornell University, Ithaca, NY, USA; Hospital for Special Surgery, New York City, NY, USA.
| |
Collapse
|
17
|
Baumann AP, Shi X, Roeder RK, Niebur GL. The sensitivity of nonlinear computational models of trabecular bone to tissue level constitutive model. Comput Methods Biomech Biomed Engin 2015; 19:465-73. [PMID: 25959510 DOI: 10.1080/10255842.2015.1041022] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Abstract
Microarchitectural finite element models have become a key tool in the analysis of trabecular bone. Robust, accurate, and validated constitutive models would enhance confidence in predictive applications of these models and in their usefulness as accurate assays of tissue properties. Human trabecular bone specimens from the femoral neck (n = 3), greater trochanter (n = 6), and lumbar vertebra (n = 1) of eight different donors were scanned by μ-CT and converted to voxel-based finite element models. Unconfined uniaxial compression and shear loading were simulated for each of three different constitutive models: a principal strain-based model, Drucker-Lode, and Drucker-Prager. The latter was applied with both infinitesimal and finite kinematics. Apparent yield strains exhibited minimal dependence on the constitutive model, differing by at most 16.1%, with the kinematic formulation being influential in compression loading. At the tissue level, the quantities and locations of yielded tissue were insensitive to the constitutive model, with the exception of the Drucker-Lode model, suggesting that correlation of microdamage with computational models does not improve the ability to discriminate between constitutive laws. Taken together, it is unlikely that a tissue constitutive model can be fully validated from apparent-level experiments alone, as the calculations are too insensitive to identify differences in the outcomes. Rather, any asymmetric criterion with a valid yield surface will likely be suitable for most trabecular bone models.
Collapse
Affiliation(s)
- Andrew P Baumann
- a Department of Aerospace and Mechanical Engineering, Bioengineering Graduate Program , University of Notre Dame , 147 Multidisciplinary Research Building, Notre Dame , IN 46556 , USA
| | - Xiutao Shi
- a Department of Aerospace and Mechanical Engineering, Bioengineering Graduate Program , University of Notre Dame , 147 Multidisciplinary Research Building, Notre Dame , IN 46556 , USA
| | - Ryan K Roeder
- a Department of Aerospace and Mechanical Engineering, Bioengineering Graduate Program , University of Notre Dame , 147 Multidisciplinary Research Building, Notre Dame , IN 46556 , USA
| | - Glen L Niebur
- a Department of Aerospace and Mechanical Engineering, Bioengineering Graduate Program , University of Notre Dame , 147 Multidisciplinary Research Building, Notre Dame , IN 46556 , USA
| |
Collapse
|
18
|
Theoretical effects of fully ductile versus fully brittle behaviors of bone tissue on the strength of the human proximal femur and vertebral body. J Biomech 2015; 48:1264-9. [PMID: 25828400 DOI: 10.1016/j.jbiomech.2015.02.066] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2014] [Revised: 02/26/2015] [Accepted: 02/28/2015] [Indexed: 11/20/2022]
Abstract
The influence of the ductility of bone tissue on whole-bone strength represents a fundamental issue of multi-scale biomechanics. To gain insight, we performed a computational study of 16 human proximal femurs and 12 T9 vertebral bodies, comparing the whole-bone strength for the two hypothetical bounding cases of fully brittle versus fully ductile tissue-level failure behaviors, all other factors, including tissue-level elastic modulus and yield stress, held fixed. For each bone, a finite element model was generated (60-82 μm element size; up to 120 million elements) and was virtually loaded in habitual (stance for femur, compression for vertebra) and non-habitual (sideways fall, only for femur) loading modes. Using a geometrically and materially non-linear model, the tissue was assumed to be either fully brittle or fully ductile. We found that, under habitual loading, changing the tissue behavior from fully ductile to fully brittle reduced whole-bone strength by 38.3±2.4% (mean±SD) and 39.4±1.9% for the femur and vertebra, respectively (p=0.39 for site difference). These reductions were remarkably uniform across bones, but (for the femur) were greater for non-habitual (57.1±4.7%) than habitual loading (p<0.001). At overall structural failure, there was 5-10-fold less failed tissue for the fully brittle than fully ductile cases. These theoretical results suggest that the whole-bone strength of the proximal femur and vertebra can vary substantially between fully brittle and fully ductile tissue-level behaviors, an effect that is relatively insensitive to bone morphology but greater for non-habitual loading.
Collapse
|
19
|
|
20
|
Sanyal A, Keaveny TM. Biaxial normal strength behavior in the axial-transverse plane for human trabecular bone--effects of bone volume fraction, microarchitecture, and anisotropy. J Biomech Eng 2014; 135:121010. [PMID: 24121715 DOI: 10.1115/1.4025679] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2013] [Indexed: 11/08/2022]
Abstract
The biaxial failure behavior of the human trabecular bone, which has potential relevance both for fall and gait loading conditions, is not well understood, particularly for low-density bone, which can display considerable mechanical anisotropy. Addressing this issue, we investigated the biaxial normal strength behavior and the underlying failure mechanisms for human trabecular bone displaying a wide range of bone volume fraction (0.06-0.34) and elastic anisotropy. Micro-computed tomography (CT)-based nonlinear finite element analysis was used to simulate biaxial failure in 15 specimens (5 mm cubes), spanning the complete biaxial normal stress failure space in the axial-transverse plane. The specimens, treated as approximately transversely isotropic, were loaded in the principal material orientation. We found that the biaxial stress yield surface was well characterized by the superposition of two ellipses--one each for yield failure in the longitudinal and transverse loading directions--and the size, shape, and orientation of which depended on bone volume fraction and elastic anisotropy. However, when normalized by the uniaxial tensile and compressive strengths in the longitudinal and transverse directions, all of which depended on bone volume fraction, microarchitecture, and mechanical anisotropy, the resulting normalized biaxial strength behavior was well described by a single pair of (longitudinal and transverse) ellipses, with little interspecimen variation. Taken together, these results indicate that the role of bone volume fraction, microarchitecture, and mechanical anisotropy is mostly accounted for in determining the uniaxial strength behavior and the effect of these parameters on the axial-transverse biaxial normal strength behavior per se is minor.
Collapse
|
21
|
Nawathe S, Akhlaghpour H, Bouxsein ML, Keaveny TM. Microstructural failure mechanisms in the human proximal femur for sideways fall loading. J Bone Miner Res 2014; 29:507-15. [PMID: 23832419 DOI: 10.1002/jbmr.2033] [Citation(s) in RCA: 65] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/18/2012] [Revised: 05/24/2013] [Accepted: 06/19/2013] [Indexed: 11/09/2022]
Abstract
The etiology of hip fractures remains unclear but might be elucidated by an improved understanding of the microstructural failure mechanisms of the human proximal femur during a sideways fall impact. In this context, we biomechanically tested 12 cadaver proximal femurs (aged 76 ± 10 years; 8 female, 4 male) to directly measure strength for a sideways fall and also performed micro-computed tomography (CT)-based, nonlinear finite element analysis of the same bones (82-micron-sized elements, ∼120 million elements per model) to estimate the amount and location of internal tissue-level failure (by ductile yielding) at initial structural failure of the femur. We found that the correlation between the directly measured yield strength of the femur and the finite element prediction was high (R(2) = 0.94, p < 0.0001), supporting the validity of the finite element simulations of failure. In these simulations, the failure of just a tiny proportion of the bone tissue (1.5% to 6.4% across all bones) led to initial structural failure of the femur. The proportion of failed tissue, estimated by the finite element models, decreased with decreasing measured femoral strength (R(2) = 0.88, p < 0.0001) and was more highly correlated with measured strength than any measure of bone volume, mass, or density. Volume-wise, trabecular failure occurred earlier and was more prominent than cortical failure in all femurs and dominated in the very weakest femurs. Femurs with low measured strength relative to their areal bone mineral density (BMD) (by dual-energy X-ray absorptiometry [DXA]) had a low proportion of trabecular bone compared with cortical bone in the femoral neck (p < 0.001), less failed tissue (p < 0.05), and low structural redundancy (p < 0.005). We conclude that initial failure of the femur during a sideways fall is associated with failure of just a tiny proportion of the bone tissue, failure of the trabecular tissue dominating in the very weakest femurs owing in part to a lack of structural redundancy.
Collapse
Affiliation(s)
- Shashank Nawathe
- Department of Mechanical Engineering, University of California, Berkeley, CA, USA
| | | | | | | |
Collapse
|
22
|
Willems NMBK, Langenbach GEJ, Everts V, Zentner A. The microstructural and biomechanical development of the condylar bone: a review. Eur J Orthod 2013; 36:479-85. [PMID: 24375755 DOI: 10.1093/ejo/cjt093] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
BACKGROUND Bone constantly strives for optimal architecture. Mandibular condyle, which is subjected to various mechanical loads forcing it to be highly adaptive, has a unique structure and a relatively high remodelling rate. Despite the eminent clinical relevance of mandibular condyle, literature on its structural and biomechanical development and on the mechanical role of its mineralized and non-mineralized bone components is scarce. OBJECTIVE The aim of the present review is to provide a brief introduction to basic bone mechanics and a synopsis of the growth and development of human mandibular condyle. Subsequently, the current ideas on the relationship between the structural and biomechanical properties of bone in general and of mandibular condyle in particular are reviewed. Finally, up-to-date knowledge from fundamental bone research will be blended with the current knowledge relevant to clinical dentistry, above all orthodontics. METHODS A comprehensive literature study was performed with an emphasis on recent and innovative work focusing on the interaction between microarchitectural and micromechanical properties of bone. CONCLUSIONS Mandibular condyle is a bone structure with a high bone turnover rate. Mechanical properties of mandibular condyle improve during adolescence and are optimal during adulthood. Local mineralization degree might not be a decisive determinant of the local bone tissue stiffness as was believed hitherto. Bone collagen and its cross links play a role in toughness and tensile strength of bone but not in its compressive properties. Clinical procedures might affect mandibular condyle, which is highly reactive to changes in its mechanical environment.
Collapse
Affiliation(s)
- Nop M B K Willems
- Departments of *Orthodontics and**Oral Cell Biology and Functional Anatomy, MOVE Research Institute, Academic Centre for Dentistry Amsterdam (ACTA), University of Amsterdam and VU University, the Netherlands
| | - Geerling E J Langenbach
- **Oral Cell Biology and Functional Anatomy, MOVE Research Institute, Academic Centre for Dentistry Amsterdam (ACTA), University of Amsterdam and VU University, the Netherlands
| | - Vincent Everts
- **Oral Cell Biology and Functional Anatomy, MOVE Research Institute, Academic Centre for Dentistry Amsterdam (ACTA), University of Amsterdam and VU University, the Netherlands
| | | |
Collapse
|
23
|
Doyle H, Lohfeld S, McHugh P. Predicting the elastic properties of selective laser sintered PCL/β-TCP bone scaffold materials using computational modelling. Ann Biomed Eng 2013; 42:661-77. [PMID: 24057867 DOI: 10.1007/s10439-013-0913-4] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2013] [Accepted: 09/12/2013] [Indexed: 11/26/2022]
Abstract
This study assesses the ability of finite element (FE) models to capture the mechanical behaviour of sintered orthopaedic scaffold materials. Individual scaffold struts were fabricated from a 50:50 wt% poly-ε-caprolactone (PCL)/β-tricalcium phosphate (β-TCP) blend, using selective laser sintering. The tensile elastic modulus of single struts was determined experimentally. High resolution FE models of single struts were generated from micro-CT scans (28.8 μm resolution) and an effective strut elastic modulus was calculated from tensile loading simulations. Three material assignment methods were employed: (1) homogeneous PCL elastic constants, (2) composite PCL/β-TCP elastic constants based on rule of mixtures, and (3) heterogeneous distribution of micromechanically-determined elastic constants. In comparison with experimental results, the use of homogeneous PCL properties gave a good estimate of strut modulus; however it is not sufficiently representative of the real material as it neglects the β-TCP phase. The rule of mixtures method significantly overestimated strut modulus, while there was no significant difference between strut modulus evaluated using the micromechanically-determined elastic constants and experimentally evaluated strut modulus. These results indicate that the multi-scale approach of linking micromechanical modelling of the sintered scaffold material with macroscale modelling gives an accurate prediction of the mechanical behaviour of the sintered structure.
Collapse
Affiliation(s)
- Heather Doyle
- Biomechanics Research Centre (BMEC), Mechanical and Biomedical Engineering, College of Engineering and Informatics, National University of Ireland Galway, Galway, Ireland,
| | | | | |
Collapse
|