Splitting fracture in bovine bone using a porosity-based spring network model

Elastic response of trabecular bone under compression calculated using the firm and floppy boundary lattice element method

Micro-Finite Element analysis (μFEA) has become widely used in biomechanical research as a reliable tool for the prediction of bone mechanical properties within its microstructure such as apparent elastic modulus and strength. However, this method requires substantial computational resources and processing time. Here, we propose a computationally efficient alternative to FEA that can provide an accurate estimation of bone trabecular mechanical properties in a fast and quantitative way. A lattice element method (LEM) framework based on the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) open-source software package is employed to calculate the elastic response of trabecular bone cores. A novel procedure to handle pore-material boundaries is presented, referred to as the Firm and Floppy Boundary LEM (FFB-LEM). Our FFB-LEM calculations are compared to voxel- and geometry-based FEA benchmarks incorporating bovine and human trabecular bone cores imaged by micro Computed Tomography (μCT). Using 14 computer cores, the apparent elastic modulus calculation of a trabecular bone core from a μCT-based input with FFB-LEM required about 15 min, including conversion of the μCT data into a LAMMPS input file. In contrast, the FEA calculations on the same system including the mesh generation, required approximately 30 and 50 min for voxel- and geometry-based FEA, respectively. There were no statistically significant differences between FFB-LEM and voxel- or geometry-based FEA apparent elastic moduli (+24.3% or +7.41%, and +0.630% or -5.29% differences for bovine and human samples, respectively).

Role of composition in fracture behavior of two-phase solids

In a two-phase solid, we examine the growth of a preexisting macroscopic crack based on simulations of a random spring network model. We find that the enhancement in toughness, as well as strength, is strongly dependent on the ratio of elastic moduli as well as on the relative proportion of the phases. We find that the mechanism that leads to enhancement in toughness is not the same as that for enhancement in strength; however, the overall enhancement is similar in mode I and mixed-mode loading. Based on the crack paths, and the spread of the fracture process zone, we identify the type of fracture to transition from nucleation type, for close to single-phase compositions, whether hard or soft, to avalanche type for more mixed compositions. We also show that the associated avalanche distributions exhibit power-law statistics with different exponents for each phase. The significance of variations in the avalanche exponents with the relative proportion of phases and possible connections to the fracture types are discussed in detail.

Interplay between disorder and hardening during tensile fracture of a quasi-brittle solid

We examine the specific role of the interplay between hardening and disorder characteristics of a representative quasi-brittle material on its failure mechanisms using a random spring network model. The model incorporates quasi-brittleness in its spring constants and disorder in the failure strain threshold and is shown to be effective in simulating the experimentally observed tensile and fracture behaviour of a quasi-brittle epoxy resin-based polymer. It is shown that rapid localization of deformation and associated damage growth occurs for a weakly hardening solid while for a linear elastic material, damage nucleates at multiple independent sites, and there is significant growth of damage, independent of other nucleating sites, prior to maximum load. The failure mechanism is shown to crossover from an avalanche-dominated fracture for a linear elastic material to nucleation-dominated fracture for a weakly hardening material.

Image-Based Polygonal Lattices for Mechanical Modeling of Biological Materials: 2D Demonstrations

Understanding the structure-property relationship of biological materials, such as bones, teeth, cells, and biofilms, is critical for diagnosing diseases and developing bioinspired materials and structures. The intrinsic multiphase heterogeneity with interfaces places great challenges for mechanical modeling. Here, we develop an image-based polygonal lattice model for simulating the mechanical deformation of biological materials with complicated shapes and interfaces. The proposed lattice model maintains the uniform meshes inside the homogeneous phases and restricts the irregular polygonal meshes near the boundaries or interfaces. This approach significantly simplifies the mesh generation from images of biological structures with complicated geometries. The conventional finite element simulations validate this polygonal lattice model. We further demonstrate that the image-based polygonal lattices generate meshes from images of composite structures with multiple inclusions and capture the nonlinear mechanical deformation. We conclude the paper by highlighting a few future research directions that will benefit from the functionalities of polygonal lattice modeling.

Role of spatial patterns in fracture of disordered multiphase materials

Multiphase materials, such as composite materials, exhibit multiple competing failure mechanisms during the growth of a macroscopic defect. For the simulation of the overall fracture process in such materials, we develop a two-phase spring network model that accounts for the architecture between the different components as well as the respective disorders in their failure characteristics. In the specific case of a plain weave architecture, we show that any offset between the layers reduces the delocalization of the stresses at the crack tip and thereby substantially lowers the strength and fracture toughness of the overall laminate. The avalanche statistics of the broken springs do not show a distinguishable dependence on the offsets between layers. The power-law exponents are found to be much smaller than that of disordered spring network models in the absence of a crack. A discussion is developed on the possibility of the avalanche statistics being those near breakdown.

[Research progress on mechanical performance evaluation of artificial intervertebral disc]

The mechanical properties of artificial intervertebral disc (AID) are related to long-term reliability of prosthesis. There are three testing methods involved in the mechanical performance evaluation of AID based on different tools: the testing method using mechanical simulator, in vitro specimen testing method and finite element analysis method. In this study, the testing standard, testing equipment and materials of AID were firstly introduced. Then, the present status of AID static mechanical properties test (static axial compression, static axial compression-shear), dynamic mechanical properties test (dynamic axial compression, dynamic axial compression-shear), creep and stress relaxation test, device pushout test, core pushout test, subsidence test, etc. were focused on. The experimental techniques using in vitro specimen testing method and testing results of available artificial discs were summarized. The experimental methods and research status of finite element analysis were also summarized. Finally, the research trends of AID mechanical performance evaluation were forecasted. The simulator, load, dynamic cycle, motion mode, specimen and test standard would be important research fields in the future.

Using Non-linear Homogenization to Improve the Performance of Macroscopic Damage Models of Trabecular Bone

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.

Role of porosity and matrix behavior on compressive fracture of Haversian bone using random spring network model

Haversian remodeling is known to result in improved resistance to compressive fracture in healthy cortical bone. Here, we examine the individual roles of the mean porosity, structure of the network of pores and remodeled bone matrix properties in the fracture behavior of Haversian bone. The detailed structure of porosity network is obtained both pre- and post-testing of dry cubical bone samples using micro-Computed Tomography. Based on the periodicity in the features of porosity along tangential direction, we develop a two dimensional porosity-based random spring network model for Haversian bone. The model is shown to capture well the macroscopic response and reproduce the avalanche statistics similar to recently reported experiments on porcine bone. The predictions suggest that at the millimeter scale, the remodeled bone matrix of Haversian bone is less stiff but tougher than that of plexiform/primary bone.

Role of matrix behavior in compressive fracture of bovine cortical bone

In compressive fracture of dry plexiform bone, we examine the individual roles of overall mean porosity, the connectivity of the porosity network, and the elastic as well as the failure properties of the nonporous matrix, using a random spring network model (RSNM). Porosity network structure is shown to reduce the compressive strength by up to 30%. However, the load-bearing capacity increases with an increase in either of the matrix properties-the elastic modulus or the failure strain threshold. To validate the porosity-based RSNM model with available experimental data, bone-specific failure strain thresholds for the ideal matrix of similar elastic properties were estimated to be within 60% of each other. Further, we observe the avalanche size exponents to be independent of the bone-dependent parameters as well as the structure of the porosity network.