An official website of the United States government
Mechanical behavior of materials with granular microstructures is confounded by unique features of their grain-scale mechano-morphology, such as the tension–compression asymmetry of grain interactions and irregular grain structure. Continuum models, necessary for the macro-scale description of these materials, must link to the grain-scale behavior to describe the consequences of this mechano-morphology. Here, we consider the damage behavior of these materials based upon purely mechanical concepts utilizing energy and variational approach. Granular micromechanics is accounted for through Piola’s ansatz and objective kinematic descriptors obtained for grain-pair relative displacement in granular materials undergoing finite deformations. Karush–Kuhn–Tucker (KKT)-type conditions that provide the evolution equations for grain-pair damage and Euler–Lagrange equations for evolution of grain-pair relative displacement are derived based upon a non-standard (hemivariational) variational approach. The model applicability is illustrated for particular form of grain-pair elastic energy and dissipation functionals through numerical examples. Results show interesting damage-induced anisotropy evolution including the emergence of a type of chiral behavior and formation of finite localization zones.
more »
« less
Computational analysis of tensile damage and failure of mineralized tissue assisted with experimental observations
In this study, deformation and failure mechanisms of mineralized tissue (bone) were investigated both experimentally and computationally by performing diametral compression tests on millimetric disk specimens and conducting finite element analysis in which a granular micromechanics-based nonlinear user-defined material model is implemented. The force–displacement relationship obtained in the simulation agreed well with the experimental results. The simulation was also able to capture location of the failure initiation observed in the experiment, which is inside out from the hole along the loading axis. Furthermore, propagation of micro-sized cracks into failure was observed both in the experiment using simultaneous slow-motion microscopy imaging and in the simulation analyzing the local distortion and local volume change within the specimen. The anisotropy evolution was found to be significant around the hole along the loading axis by evaluating the anisotropy index computed using finite element results. In conclusion, this work revealed that the prediction capability of granular micromechanics-based user-defined nonlinear material model (UMAT) is promising considering the match between the results and observations from the physical experiment and finite element analysis such as force–displacement relationship and failure initiation/pattern. This work has also shown that the tensile damage and failure of mineralized tissues can be characterized using diametral compression (split tension) test.
more »
« less
- Award ID(s):
- 1727433
- PAR ID:
- 10549494
- Publisher / Repository:
- SAGE Publications
- Date Published:
- Journal Name:
- Proceedings of the Institution of Mechanical Engineers, Part H: Journal of Engineering in Medicine
- Volume:
- 234
- Issue:
- 3
- ISSN:
- 0954-4119
- Format(s):
- Medium: X Size: p. 289-298
- Size(s):
- p. 289-298
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
At the nanoscale bone is composed of aligned mineralized collagen fibrils organized into packets along the surface of trabecular bone creating an anisotropic tissue microstructure. Newer packets at the trabecular surfaces are usually less mineralized than older bone in the interior of the trabeculae, which along with irregular mineral deposition within packets, forms a heterogeneous material across the span of a trabeculae. However, finite element (FE) models of bone typically use homogenous isotropic material properties, because it is challenging to build anisotropy and heterogeneity into a model in a way that is applicable to the complex geometries of trabecular bone. Both the material anisotropy and heterogeneity may influence the stress state of trabecular bone, and it is important to understand the implications of such differences for determining bone biomechanical failure. It was hypothesized that taking into consideration both the tissue anisotropy and heterogeneity of bone’s biomechanical properties would alter the expected failure locations by reducing tensile stress on near surface elements of an FE model of canine trabecular bone. The objective of this study was to test this hypothesis and to develop a method to apply anisotropic and heterogeneous material properties to a model automatically from micro-computed tomography (μCT) data.more » « less
-
Trabecular bone is composed of organized mineralized collagen fibrils, which results in heterogeneous and anisotropic mechanical properties at the tissue level. Recently, biomechanical models computing stresses and strains in trabecular bone have indicated a significant effect of tissue heterogeneity on predicted stresses and strains. How-ever, the effect of the tissue-level mechanical anisotropy on the trabecular bone biomechanical response is unknown. Here, a computational method was established to automatically impose physiologically relevant orientation inherent in trabecular bone tissue on a trabecular bone microscale finite element model. Spatially varying tissue-level anisotropic elastic properties were then applied according to the bone mineral density and the local tissue orientation. The model was used to test the hypothesis that anisotropy in both homogeneous and heterogeneous models alters the predicted distribution of stress invariants. Linear elastic finite element computations were performed on a 3 mm cube model isolated from a microcomputed tomography scan of human trabecular bone from the distal femur. Hydrostatic stress and von Mises equivalent stress were recorded at every element, and the distributions of these values were analyzed. Anisotropy reduced the range of hydrostatic stress in both tension and compression more strongly than the associated increase in von Mises equivalent stress. The effect of anisotropy was independent of the spatial redistribution high compressive stresses due to tissue elastic heterogeneity. Tissue anisotropy and heterogeneity are likely important mechanisms to protect bone from failure and should be included for stress analyses in trabecular bone.more » « less
-
The work reported in ``Granular micromechanics-based identification of isotropic strain gradient parameters for elastic geometrically nonlinear deformations" misidentified key terms in the grain-pair objective relative displacement when accounting for the second gradient of placement. In this paper, we correct that oversight by deriving a revised expression for the grain-pair objective relative displacement within the granular micromechanics framework. The amended terms, which resemble Christoffel symbols expressed in terms of strain gradients, modify the contributions of both the normal and tangential components to the strain energy and, consequently, alter the identified strain gradient elastic parameters. Importantly, the identification of the standard (first gradient) elastic tensor remains unchanged. This brief paper presents the corrected derivation, the resulting stiffness tensors for anisotropic strain gradient elasticity, and updated analytical expressions for the material parameters in both 2D and 3D isotropic settings.more » « less
-
Abstract A preceding 2023 study argued that the resistance of a heterogeneous material to the curvature of the displacement field is the most physically realistic localization limiter for softening damage. The curvature was characterized by the second gradient of the displacement vector field, which includes the material rotation gradient, and was named the “sprain” tensor, while the term “spress” is here proposed as the force variable work-conjugate to “sprain.” The partial derivatives of the associated sprain energy density yielded in the preceeding study, sets of curvature resisting self-equilibrated nodal sprain forces. However, the fact that the sprain forces had to be applied on the adjacent nodes of a finite element greatly complicated the programming and extended the simulation time in a commercial code such as abaqus by almost two orders of magnitude. In the present model, Smooth Lagrangian Crack Band Model (slCBM), these computational obstacles are here overcome by using finite elements with linear shape functions for both the displacement vector and for an approximate displacement gradient tensor. A crucial feature is that the nodal values of the approximate gradient tensor are shared by adjacent finite elements. The actual displacement gradient tensor calculated from the nodal displacement vectors is constrained to the approximate displacement gradient tensor by means of a Lagrange multiplier tensor, either one for each element or one for each node. The gradient tensor of the approximate gradient tensor then represents the approximate third-order displacement curvature tensor, or Hessian of the displacement field. Importantly, the Lagrange multiplier behaves as an externally applied generalized moment density that, similar to gravity, does not affect the total strain-plus-sprain energy density of material. The Helmholtz free energy of the finite element and its associated stiffness matrix are formulated and implemented in a user’s element of abaqus. The conditions of stationary values of the total free energy of the structure with respect to the nodal degrees-of-freedom yield the set of equilibrium equations of the structure for each loading step. One- and two-dimensional examples of crack growth in fracture specimens are given. It is demonstrated that the simulation results of the three-point bend test are independent of the orientation of a regular square mesh, capture the width variation of the crack band, the damage strain profile across the band, and converge as the finite element mesh is refined.more » « less
