Modeling of bending and torsional stiffnesses of bone at sub-microscale: Effect of curved mineral lamellae.

Recent transmission electron microscopy images of transverse sections of human cortical bone showed that mineral lamellae (polycrystalline sheets of apatite crystals) form arcuate multi-radius patterns around collagen fibrils. The 3-6 nm thick mineral lamellae are arranged in stacks of 3-20 layers and curve around individual fibrils, few fibrils, and higher numbers of collagen fibrils. We evaluate the effect of these stacked mineral lamellae with various radius of curvature patterns on the elastic bending and torsional responses of bone at the sub-microscale using a finite element method. We find that the curved multi-radius stack patterns increased the bending and torsional stiffnesses by 7% and 23%, respectively, compared to when the stacks of mineral lamellae only encircle individual fibrils for the idealized geometric models considered. This study provides new insights into the structure-property relations for the bone ultrastructure.

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.

Modeling of Osteoprobe indentation on bone.

Osteoprobe (ActiveLife, Santa Barbara, CA) is a novel handheld microindentation instrument designed to test bone in vivo by measuring a Bone Material Strength index (BMSi). In this paper, the Osteoprobe indentation on a cortical bone is modeled computationally to gain insights into the physical interpretation of the BMSi output. The analysis is conducted using an axisymmetric finite element model with an isotropic viscoelastic-plastic constitutive law with continuum damage. The computational model is validated by comparing it with experimental data from the literature. Experimental factors (indenter tip radius and friction coefficient between the indenter and the bone) and four mechanical properties of bone (Young's modulus, compressive yield stress, and damage and viscosity constants) are varied to study their influence on the BMSi. We find that varying the friction coefficient can proportionally change the BMSi up to 3%. The indenter tip radius is proportional to the BMSi, with a more pronounced proportional relation when it is greater than 30 µm. Young's modulus has a proportional relation with the BMSi, where decreasing it by 73% reduces the BMSi by 41%. The damage constant has an inversely proportional relation to the BMSi, where increasing it from 0.5 to 0.96 reduces the BMSi by 29%. The compressive yield stress and the viscosity constant have a close proportional relation to the BMSi, where increasing the compressive yield stress from 50 MPa to 200 MPa increases the Osteoprobe BMSi by 21%. In summary, the friction coefficient and the indenter tip radius (when smaller than 30 µm) have a small effect on BMSi. Young's modulus and damage have stronger relations with the BMSi than compressive yield stress and viscosity constant. This fundamental study provides new insights into the BMSi measurement and serves as a basis for further computational and experimental investigations on the Osteoprobe technique. Such research is needed to facilitate the embrace of this technique by the clinical community.

Cortical bone fracture analysis using XFEM - case study.

We aim to achieve an accurate simulation of human cortical bone fracture using the extended finite element method within a commercial finite element software abaqus. A two-dimensional unit cell model of cortical bone is built based on a microscopy image of the mid-diaphysis of tibia of a 70-year-old human male donor. Each phase of this model, an interstitial bone, a cement line, and an osteon, are considered linear elastic and isotropic with material properties obtained by nanoindentation, taken from literature. The effect of using fracture analysis methods (cohesive segment approach versus linear elastic fracture mechanics approach), finite element type, and boundary conditions (traction, displacement, and mixed) on cortical bone crack initiation and propagation are studied. In this study cohesive segment damage evolution for a traction separation law based on energy and displacement is used. In addition, effects of the increment size and mesh density on analysis results are investigated. We find that both cohesive segment and linear elastic fracture mechanics approaches within the extended finite element method can effectively simulate cortical bone fracture. Mesh density and simulation increment size can influence analysis results when employing either approach, and using finer mesh and/or smaller increment size does not always provide more accurate results. Both approaches provide close but not identical results, and crack propagation speed is found to be slower when using the cohesive segment approach. Also, using reduced integration elements along with the cohesive segment approach decreases crack propagation speed compared with using full integration elements. Copyright © 2016 John Wiley & Sons, Ltd.

Finite element simulation of Reference Point Indentation on bone.

Reference Point Indentation (RPI) is a novel technique aimed to assess bone quality. Measurements are recorded by the BioDent instrument that applies multiple indents to the same location of cortical bone. Ten RPI parameters are obtained from the resulting force-displacement curves. Using the commercial finite element analysis software Abaqus, we assess the significance of the RPI parameters. We create an axisymmetric model and employ an isotropic viscoelastic-plastic constitutive relation with damage to simulate indentations on a human cortical bone. Fracture of bone tissue is not simulated for simplicity. The RPI outputs are computed for different simulated test cases and then compared with experimental results, measured using the BioDent, found in literature. The number of cycles, maximum indentation load, indenter tip radius, and the mechanical properties of bone: Young׳s modulus, compressive yield stress, and viscosity and damage constants, are varied. The trends in the RPI parameters are then investigated. We find that the RPI parameters are sensitive to the mechanical properties of bone. An increase in Young׳s modulus of bone causes the force-displacement loading and unloading slopes to increase and the total indentation distance (TID) to decrease. The compressive yield stress is inversely proportional to a creep indentation distance (CID1) and the TID. The viscosity constant is proportional to the CID1 and an average of the energy dissipated (AvED). The maximum indentation load is proportional to the TID, CID1, loading and unloading slopes, and AvED. The damage parameter is proportional to the TID, but it is inversely proportional to both the loading and unloading slopes and the AvED. The value of an indenter tip radius is proportional to the CID1 and inversely proportional to the TID. The number of load cycles is inversely proportional to an average of a creep indentation depth (AvCID) and the AvED. The indentation distance increase (IDI) is strongly inversely proportional to the compressive yield stress, and strongly proportional to the viscosity constant and maximum applied load, but has weak relation with the damage parameter, indenter tip radius, and elastic modulus. This computational study advances our understanding of the RPI outputs and provides a starting point for more comprehensive computational studies of the RPI technique.

Toward high-speed 3D nonlinear soft tissue deformation simulations using Abaqus software.

We aim to achieve a fast and accurate three-dimensional (3D) simulation of a porcine liver deformation under a surgical tool pressure using the commercial finite element software Abaqus. The liver geometry is obtained using magnetic resonance imaging, and a nonlinear constitutive law is employed to capture large deformations of the tissue. Effects of implicit versus explicit analysis schemes, element type, and mesh density on computation time are studied. We find that Abaqus explicit and implicit solvers are capable of simulating nonlinear soft tissue deformations accurately using first-order tetrahedral elements in a relatively short time by optimizing the element size. This study provides new insights and guidance on accurate and relatively fast nonlinear soft tissue simulations. Such simulations can provide force feedback during robotic surgery and allow visualization of tissue deformations for surgery planning and training of surgical residents.