Nonlinear micro-CT based FE modeling of trabecular bone-Sensitivity of apparent response to tissue constitutive law and bone volume fraction

Investigation on primary stability of dental implants through considering peri-implant bone damage, caused by small and large deformations: A validated non-linear micro finite element study

Primary stability of a dental implant is defined as its ability to resist the applied load without showing excessive damage in peri-implant bone, which is a prerequisite for secondary stability, and consequently for implantation success. The main goal of this study was to develop a validated micro-finite element (μFE) approach to assess the primary stability of dental implants in terms of stiffness, stiffness reduction, and irreversible displacement of the bone-implant system, subjected to an increasing step-wise quasi-static compressive loading-unloading test. The μFE models were generated based on the μCT images of bone, taken from extracted bovine tibia trabecular bone samples after drilling and implantation. A tissue constitutive model was considered for trabecular bone by describing elasto-plasticity with a modified von Mises yield criterion and element deletion technique to account for trabecular bone damage behavior. Then, the obtained force-displacement curves from the simulation were compared with the in-vitro mechanical test curves to evaluate the validity of the model. The results showed that the proposed μFE model could be properly predict the bone-implant system mechanical response in terms of irreversible displacement (R2 = 0.99), stiffness (R2 = 0.77), and stiffness reduction (R2 = 0.72) of the bone-implant construct for all the applied displacements without a significant difference from the unit slope and zero intercept of the QQ-plot (p-value<0.05). Moreover, a qualitative agreement was seen between the peri-implant bone damage predicted by the μFE model and the observed from μCT images. The adopted methodology used in this study can predict the mechanical failure response of the bone-implant system, which can be employed as a representative tool to study the effects of various dental implant design parameters on the primary stability with the ultimate goal of optimizing dental implants design.

Comparison of linear and nonlinear stepwise μFE displacement predictions to digital volume correlation measurements of trabecular bone biopsies

Digital volume correlation (DVC) enables to evaluate the ability of μFE models in predicting experimental results on the mesoscale. In this study predicted displacement fields of three different linear and materially nonlinear μFE simulation methods were compared to DVC measured displacement fields at specific load steps in the elastic regime (Step_El) and after yield (Step_Ult). Five human trabecular bone biopsies from a previous study were compressed in several displacement steps until failure. At every compression step, μCT images (resolution: 36 μm) were recorded. A global DVC algorithm was applied to compute the displacement fields at all loading steps. The unloaded 3D images were then used to generate homogeneous, isotropic, linear and materially nonlinear μFE models. Three different μFE simulation methods were used: linear (L), nonlinear (NL), and nonlinear stepwise (NLS). Regarding L and NL, the boundary conditions were derived from the interpolated displacement fields at Step_El and Step_Ult, while for the NLS method nonlinear changes of the boundary conditions of the experiments were captured using the DVC displacement field of every available load step until Step_El and Step_Ult. The predicted displacement fields of all μFE simulation methods were in good agreement with the DVC measured displacement fields (individual specimens: R²>0.83 at Step_El and R²>0.59 at Step_Ult; pooled data: R²>0.97 at Step_El and R²>0.92 at Step_Ult). At Step_El, all three simulation methods showed similar intercepts, slopes, and coefficients of determination while the nonlinear μFE models improved the prediction of the displacement fields slightly in all Cartesian directions at Step_Ult (individual specimens: L: R²>0.59 and NL, NLS: R²>0.68; pooled data: L: R²>0.92 and NL, NLS: R²>0.94). Damaged/overstrained elements in L, NL, and NLS occurred at similar locations but the number of overstrained elements was overestimated when using the L simulation method. Considering the increased solving time of the nonlinear μFE models as well as the acceptable performance in displacement prediction of the linear μFE models, one can conclude that for similar use cases linear μFE models represent the best compromise between computational effort and accuracy of the displacement field predictions.

Toward an artificial intelligence-assisted framework for reconstructing the digital twin of vertebra and predicting its fracture response

This article presents an effort toward building an artificial intelligence (AI) assisted framework, coined ReconGAN, for creating a realistic digital twin of the human vertebra and predicting the risk of vertebral fracture (VF). ReconGAN consists of a deep convolutional generative adversarial network (DCGAN), image-processing steps, and finite element (FE) based shape optimization to reconstruct the vertebra model. This DCGAN model is trained using a set of quantitative micro-computed tomography (micro-QCT) images of the trabecular bone obtained from cadaveric samples. The quality of synthetic trabecular models generated using DCGAN are verified by comparing a set of its statistical microstructural descriptors with those of the imaging data. The synthesized trabecular microstructure is then infused into the vertebra cortical shell extracted from the patient's diagnostic CT scans using an FE-based shape optimization approach to achieve a smooth transition between trabecular to cortical regions. The final geometrical model of the vertebra is converted into a high-fidelity FE model to simulate the VF response using a continuum damage model under compression and flexion loading conditions. A feasibility study is presented to demonstrate the applicability of digital twins generated using this AI-assisted framework to predict the risk of VF in a cancer patient with spinal metastasis.

Damage Function of a Quasi-Brittle Material, Damage Rate, Acceleration and Jerk during Uniaxial Compression: Model and Application to Analysis of Trabecular Bone Tissue Destruction

A diversity of quasi-brittle materials can be observed in various engineering structures and natural objects (rocks, frozen soil, concrete, ceramics, bones, etc.). In order to predict the condition and safety of these objects, a large number of studies aimed at analyzing the strength of quasi-brittle materials has been conducted and presented in publications. However, at the modeling level, the problem of estimating the rate and acceleration of destruction of a quasi-brittle material under loading remains relevant. The purpose of the study was to substantiate the function of damage to a quasi-brittle material under uniaxial compression, determine the rate, acceleration and jerk of the damage process, and also to apply the results obtained to predicting the destruction of trabecular bone tissue. In accordance with the purpose of the study, the basic concepts of fracture mechanics and standard methods of mathematical modeling were used. The proposed model is based on the application of the previously obtained differentiable damage function without parameters. The results of the study are presented in the form of plots and analytical relations for computing the rate, acceleration and jerk of the damage process. Examples are given. The predicted peak of the combined effect of rate, acceleration and jerk of the damage process are found to be of practical interest as an additional criterion for destruction. The simulation results agree with the experimental data known from the available literature.

A Damage Model to Trabecular Bone and Similar Materials: Residual Resource, Effective Elasticity Modulus, and Effective Stress under Uniaxial Compression

Experimental research of bone strength remains costly and limited for ethical and technical reasons. Therefore, to predict the mechanical state of bone tissue, as well as similar materials, it is desirable to use computer technology and mathematical modeling. Yet, bone tissue as a bio-mechanical object with a hierarchical structure is difficult to analyze for strength and rigidity; therefore, empirical models are often used, the disadvantage of which is their limited application scope. The use of new analytical solutions overcomes the limitations of empirical models and significantly improves the way engineering problems are solved. Aim of the paper: the development of analytical solutions for computer models of the mechanical state of bone and similar materials. Object of research: a model of trabecular bone tissue as a quasi-brittle material under uniaxial compression (or tension). The new ideas of the fracture mechanics, as well as the methods of mathematical modeling and the biomechanics of bone tissues were used in the work. Compression and tension are considered as asymmetric mechanical states of the material. Results: a new nonlinear function that simulates both tension and compression is justified, analytical solutions for determining the effective and apparent elastic modulus are developed, the residual resource function and the damage function are justified, and the dependences of the initial and effective stresses on strain are obtained. Using the energy criterion, it is proven that the effective stress continuously increases both before and after the extremum point on the load-displacement plot. It is noted that the destruction of bone material is more likely at the inflection point of the load-displacement curve. The model adequacy is explained by the use of the energy criterion of material degradation. The results are consistent with the experimental data available in the literature.

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

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.

A Review on Recent Advances in the Constitutive Modeling of Bone Tissue

PURPOSE OF REVIEW

Image-based finite element analysis (FEA) to predict and understand the biomechanical response has become an essential methodology in musculoskeletal research. An important part of such simulation models is the constitutive material model of which recent advances are summarized in this review.

RECENT FINDINGS

The review shows that existing models from other fields were introduced, such as cohesion zone (cortical bone) or phase-field models (trabecular bone). Some progress has been made in describing cortical bone involving physical mechanisms such as microcracks. Problems with validations at different length scales remain a problem. The improvement of recent constitutive models is partially obscured by uncertainties that affect overall predictions, such as image quality and calibration or boundary conditions. Nevertheless, in vivo CT-based FEA simulations based on a sophisticated constitutive behavior are a very valuable tool for clinical-related osteoporosis research.

Image-based finite-element modeling of the human femur

Fracture is considered a critical clinical endpoint in skeletal pathologies including osteoporosis and bone metastases. However, current clinical guidelines are limited with respect to identifying cases at high risk of fracture, as they do not account for many mechanical determinants that contribute to bone fracture. Improving fracture risk assessment is an important area of research with clear clinical relevance. Patient-specific numerical musculoskeletal models generated from diagnostic images are widely used in biomechanics research and may provide the foundation for clinical tools used to quantify fracture risk. However, prior to clinical translation, in vitro validation of predictions generated from such numerical models is necessary. Despite adopting radically different models, in vitro validation of image-based finite element (FE) models of the proximal femur (predicting strains and failure loads) have shown very similar, encouraging levels of accuracy. The accuracy of such in vitro models has motivated their application to clinical studies of osteoporotic and metastatic fractures. Such models have demonstrated promising but heterogeneous results, which may be explained by the lack of a uniform strategy with respect to FE modeling of the human femur. This review aims to critically discuss the state of the art of image-based femoral FE modeling strategies, highlighting principal features and differences among current approaches. Quantitative results are also reported with respect to the level of accuracy achieved from in vitro evaluations and clinical applications and are used to motivate the adoption of a standardized approach/workflow for image-based FE modeling of the femur.

Patient-Specific Bone Multiscale Modelling, Fracture Simulation and Risk Analysis-A Survey

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.