High-Performance Computing Comparison of Implicit and Explicit Nonlinear Finite Element Simulations of Trabecular Bone

Comput Methods Programs Biomed. 2021 Mar:200:105870. doi: 10.1016/j.cmpb.2020.105870. Epub 2020 Nov 23.

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.

Keywords: Explicit solver; High performance computing; Implicit solver; Nonlinear finite element analysis; Trabecular bone.

MeSH terms

  • Cancellous Bone* / diagnostic imaging
  • Computer Simulation
  • Finite Element Analysis
  • Stress, Mechanical
  • X-Ray Microtomography