skip to main content


Title: 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
NSF-PAR ID:
10168983
Author(s) / Creator(s):
;
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
Page Range / eLocation ID:
289 to 298
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Topologically interlocked materials (TIMs) are material systems consisting of one or more repeating unit blocks assembled in a planar configuration such that each block is fully constrained geometrically by its neighbours. The assembly is terminated by a frame that constrains the outermost blocks. The resulting plate-like structure does not use any type of adhesive or fastener between blocks but is capable of carrying transverse loads. These material systems are advantageous due to their potential attractive combination of strength, toughness, and damage tolerance as compared to monolithic plates, especially when using lower strength materials. TIMs are damage tolerant due to the fact that cracks in any single block cannot propagate to neighbouring blocks. Many configurations of TIMs have been conceptualized in the past, particularly in architecture, but less work has been done to understand the mechanics of such varied assembly architectures. This work seeks to expand our knowledge of how TIM architecture is related to TIM mechanics. The present study considers TIMs created from the Archimedean and Laves tessellations. Each tessellation is configured as a TIM by projecting each edge of a tile at alternating angles from the normal to the tiling plane. For each tiling, multiple symmetries exist depending on where the frame is placed relative to the tiling. Six unique tilings and their multiple symmetries and load directions were considered, resulting in 19 unique TIM configurations. All TIM configurations were realized with identical equivalent overall assembly dimensions. The radius of the inscribed circle of the square and hexagon frames were the same, as well as the thickness of the assemblies. The tilings were scaled to possess the similar same number of building blocks within the frame. Finite element models were created for each configuration and subjected to two load types under quasi-static conditions: a prescribed displacement applied at the center of the assembly, and by a gravity load. The force deflection response of all TIM structures was found to be similar to that of a Mises truss, comprised of an initial positive stiffness followed by a period of negative stiffness until failure of the assembly. This response is indeed related to the internal working of load transfer in TIMs. Owing to the granular type character of the TIM assembly, the stress distribution follows a force-network. The key findings of this study are: • The load transfer in TIMs follows from force networks and the geometry of the force network is associated with the dual tessellation of the respective TIM system. • In TIMs based on Laves tessellations (centered around a vertex of the tiling rather than the center of a tile), displayed chirality and exerted a moment normal to the tile plane as they were loaded. • TIMs resulting from tessellations with more than one unique tile, such as squares and octagons, are asymmetric along the normal to the tile plane causing a dependence of the load response to the direction of the transverse load. Work is underway to transform these findings into general rules allowing for a predictive relationship between material architecture and mechanical response of TIM systems. This material is based upon work supported by the National Science Foundation under Grant No. 1662177. 
    more » « less
  2. In this manuscript, we present a multiscale Adaptive Reduced-Order Modeling (AROM) framework to efficiently simulate the response of heterogeneous composite microstructures under interfacial and volumetric damage. This framework builds on the eigendeformation-based reduced-order homogenization model (EHM), which is based on the transformation field analysis (TFA) and operates in the context of computational homogenization with a focus on model order reduction of the microscale problem. EHM pre-computes certain microstructure information by solving a series of linear elastic problems defined over the fully resolved microstructure (i.e., concentration tensors, interaction tensors) and approximates the microscale problem using a much smaller basis spanned over subdomains (also called parts) of the microstructure. Using this reduced basis, and prescribed spatial variation of inelastic response fields over the parts, the microscale problem leads to a set of algebraic equations with part-wise responses as unknowns, instead of node-wise displacements as in finite element analysis. The volumetric and interfacial influence functions are calculated by using the Interface enriched Generalized Finite Element Method (IGFEM) to compute the coefficient tensors, in which the finite element discretization does not need to conform to the material interfaces. AROM takes advantage of pre-computed coefficient tensors associated with the finest ROM and efficiently computes the coefficient tensors of a series of gradually coarsening ROMs. During the multiscale analysis stage, the simulation starts with a coarse ROM which can capture the initial elastic response well. As the loading continues and response in certain parts of the microstructure starts to localize, the analysis adaptively switches to the next level of refined ROM to better capture those local responses. The performance of AROM is evaluated by comparing the results with regular EHM (no adaptive refinement) and IGFEM under different loading conditions and failure modes for various 2D and 3D microstructures. The proposed AROM provides an efficient way to model history-dependent nonlinear responses for composite materials under localized interface failure and phase damage. 
    more » « less
  3. null (Ed.)
    This research presents an experimental program executed to understand the strength and stiffness properties of hollow built-up glass compression members that are intended for use in the modular construction of all glass, compression-dominant, shell-type structures. The proposed compression-dominant geometric form has been developed using the methods of form finding and three-dimensional graphical statics. This research takes the first steps towards a new construction methodology for glass structures where individual hollow glass units (HGU) are assembled using an interlocking system to form large, compression-dominant, shell-type structures, thereby exploiting the high compression strength of glass. In this study, an individual HGU has an elongated hexagonal prism shape and consists of two deck plates, two long side plates, and four short side plates, as is shown in Figure 1. Connections between glass plates are made using a two-sided transparent structural adhesive tape. The test matrix includes four HGUs, two each fabricated with 1 mm and 2 mm thick adhesive tape. All samples are dimensioned 64 cm on the long axis of symmetry, 51 cm on the short axis of symmetry, and are 10 cm in width. Glass plates are all 10 mm thick annealed float glass with geometric fabrication done using 5-axis abrasive water jet cutting. HGU assembly is accomplished using 3D printed truing clips and results in a rigid three-dimensional glass frame. Testing was done with the HGU oriented such that load was introduced on the short side edges of the two deck plates, resulting in an asymmetric load-support condition. A soft interface material was used between the HGU and steel plates of the hydraulic actuator and support for the purpose of avoiding premature cracking from local stress concentrations on the glass edges at the load and support locations. Force was applied in displacement control at 0.25 mm/minute with a full array of displacement and strain sensors. Test results for load vs. center deck plate transverse deflection are shown in Figure 2. All samples failed explosively by flexural buckling with no premature cracking on the load and support edges of the deck plates. Strain and deformation data clearly show the presence of second-order behavior resulting from bending deformation perpendicular to the plane of the deck plates. In general, linear axial behavior transitions to nonlinear second-order behavior, with increasing rates in deflection and strain growth, ultimately ending in glass fracture on the tension surfaces of the buckled deck plates. The failure resulted in near-complete disintegration of the deck plates, but with no observable cracking in any side plates and a secure connection on all adhesive tape. Results of the experimental program clearly demonstrate the feasibility of using HGUs for modular construction of compression dominant all-glass shell-type structures. This method of construction can significantly reduce the self-weight of the structure, and it will inspire the use of sustainable materials in the construction of efficient structures. 
    more » « less
  4. Abstract

    Micromechanics‐based continuum damage criteria have previously been developed to simulate the initiation of ductile fracture in structural steels under conditions with large‐scale plasticity where conventional fracture mechanics indices are invalid. Such models have been combined with methods to simulate ductile crack growth for monotonic loading. In this study, a micromechanics‐based adaptive cohesive zone model for simulating ductile crack propagation under monotonic loading is extended to handle cyclic loading. The proposed model adaptively modifies the cohesive traction–separation relationship for crack opening and closure, as the loading reverses between tension and compression. The approach is implemented into the finite element analysis platform WARP3D, and results of simulations that use the model are compared with data from coupon‐scale tests. The results demonstrate that the proposed model can accurately simulate the effect of crack propagation on specimen response, as well as other key aspects of observed behavior, including crack face closure and crack tunneling.

     
    more » « less
  5. Hydrogels are a class of soft, highly deformable materials formed by swelling a network of polymer chains in water. With mechanical properties that mimic biological materials, hydrogels are often proposed for load bearing biomedical or other applications in which their deformation and failure properties will be important. To study the failure of such materials a means for the measurement of deformation fields beyond simple uniaxial tension tests is required. As a non-contact, full-field deformation measurement method, Digital Image Correlation (DIC) is a good candidate for such studies. The application of DIC to hydrogels is studied here with the goal of establishing the accuracy of DIC when applied to hydrogels in the presence of large strains and large strain gradients. Experimental details such as how to form a durable speckle pattern on a material that is 90% water are discussed. DIC is used to measure the strain field in tension loaded samples containing a central hole, a circular edge notch and a sharp crack. Using a nonlinear, large deformation constitutive model, these experiments are modeled using the finite element method (FEM). Excellent agreement between FEM and DIC results for all three geometries shows that the DIC measurements are accurate up to strains of over 10, even in the presence of very high strain gradients near a crack tip. The method is then applied to verify a theoretical prediction that the deformation field in a cracked sample under relaxation loading, i.e. constant applied boundary displacement, is stationary in time even as the stress relaxes by a factor of three. 
    more » « less