skip to main content


This content will become publicly available on December 16, 2024

Title: Microstructurally-informed stochastic inhomogeneity of material properties and material symmetries in 3D-printed 316 L stainless steel
Stochastic mesoscale inhomogeneity of material properties and material symmetries are investigated in a 3D-printed material. The analysis involves a spatially-dependent characterization of the microstructure in 316 L stainless steel, obtained through electron backscatter diffraction imaging. These data are subsequently fed into a Voigt–Reuss–Hill homogenization approxima- tion to produce maps of elasticity tensor coefficients along the path of experimental probing. Information-theoretic stochastic models corresponding to this stiffness random field are then introduced. The case of orthotropic fields is first defined as a high-fidelity model, the realizations of which are consistent with the elasticity maps. To investigate the role of material symmetries, an isotropic approximation is next introduced through ad-hoc projections (using various metrics). Both stochastic representations are identified using the dataset. In particular, the correlation length along the characterization path is identified using a maximum likelihood estimator. Uncertainty propagation is finally performed on a complex geometry, using a Monte Carlo analysis. It is shown that mechanical predictions in the linear elastic regime are mostly sensitive to material symmetry but weakly depend on the spatial correlation length in the considered propagation scenario.  more » « less
Award ID(s):
1942928
NSF-PAR ID:
10480298
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
Springer
Date Published:
Journal Name:
Computational Mechanics
ISSN:
0178-7675
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Current thoroughly described biodegradable and cross‐linkable polymers mainly rely on acrylate cross‐linking. However, despite the swift cross‐linking kinetics of acrylates, the concomitant brittleness of the resulting materials limits their applicability. Here, photo‐cross‐linkable poly(ε‐caprolactone) networks through orthogonal thiol‐ene chemistry are introduced. The step‐growth polymerized networks are tunable, predictable by means of the rubber elasticity theory and it is shown that their mechanical properties are significantly improved over their acrylate cross‐linked counterparts. Tunability is introduced to the materials, by alteringMc(or the molar mass between cross‐links), and its effect on the thermal properties, mechanical strength and degradability of the materials is evaluated. Moreover, excellent volumetric printability is illustrated and the smallest features obtained via volumetric 3D‐printing to date are reported, for thiol‐ene systems. Finally, by means of in vitro and in vivo characterization of 3D‐printed constructs, it is illustrated that the volumetrically 3D‐printed materials are biocompatible. This combination of mechanical stability, tunability, biocompatibility, and rapid fabrication by volumetric 3D‐printing charts a new path toward bedside manufacturing of biodegradable patient‐specific implants.

     
    more » « less
  2. Accurate characterization of the mechanical properties of the human brain at both microscopic and macroscopic length scales is a critical requirement for modeling of traumatic brain injury and brain folding. To date, most experimental studies that employ classical tension/compression/shear tests report the mechanical properties of the brain averaged over both the gray and white matter within the macroscopic regions of interest. As a result, there is a missing correlation between the independent mechanical properties of the microscopic constituent elements and the composite bulk macroscopic mechanical properties of the tissue. This microstructural computational study aims to inversely predict the hyperelastic mechanical properties of the axonal fibers and their surrounding extracellular matrix (ECM) from the bulk tissue's mechanical properties. We develop a representative volume element (RVE) model of the bulk tissue consisting of axonal fibers and ECM with the embedded element technique. A multiobjective optimization technique is implemented to calibrate the model and establish the independent mechanical properties of axonal fibers and ECM based on seven previously reported experimental mechanical tests for bulk white matter tissue from the corpus callosum. The result of the study shows that the discrepancy between the reported values for the elastic behavior of white matter in literature stems from the anisotropy of the tissue at the microscale. The shear modulus of the axonal fiber is seven times larger than the ECM, with axonal fibers that also show greater nonlinearity, contrary to the common assumption that both components exhibit identical nonlinear characteristics. Statement of significance The reported mechanical properties of white matter microstructure used in traumatic brain injury or brain mechanics studies vary widely, in some cases by up to two orders of magnitude. Currently, the material parameters of the white matter microstructure are identified by a single loading mode or ultimately two modes of the bulk tissue. The presented material models only define the response of the bulk and homogenized white matter at a macroscopic scale and cannot explicitly capture the connection between the material properties of microstructure and bulk structure. To fill this knowledge gap, our study characterizes the hyperelastic material properties of axonal fibers and ECM using microscale computational modeling and multiobjective optimization. The hyperelastic material properties for axonal fibers and ECM presented in this study are more accurate than previously proposed because they have been optimized using seven or six loading modes of the bulk tissue, which were previously limited to only two of the seven possible loading modes. As such, the predicted values with high accuracy could be used in various computational modeling studies. The systematic characterization of the material properties of the human brain tissue at both macro- and microscales will lead to more accurate computational predictions, which will enable a better understanding of injury criteria, and has a positive impact on the improved development of smart protection systems, and more accurate prediction of brain development and disease progression. 
    more » « less
  3. Abstract

    Materials based on minimal surface geometries have shown superior strength and stiffness at low densities, which makes them promising continuous‐based material platforms for a variety of engineering applications. In this work, it is demonstrated how these mechanical properties can be complemented by dynamic functionalities resulting from robust topological guiding of elastic waves at interfaces that are incorporated into the considered material platforms. Starting from the definition of Schwarz P minimal surface, geometric parametrizations are introduced that break spatial symmetry by forming 1D dimerized and 2D hexagonal minimal surface‐based materials. Breaking of spatial symmetries produces topologically non‐trivial interfaces that support the localization of vibrational modes and the robust propagation of elastic waves along pre‐defined paths. These dynamic properties are predicted through numerical simulations and are illustrated by performing vibration and wave propagation experiments on additively manufactured samples. The introduction of symmetry‐breaking topological interfaces through parametrizations that modify the geometry of periodic minimal surfaces suggests a new strategy to supplement the load‐bearing properties of this class of materials with novel dynamic functionalities.

     
    more » « less
  4. null (Ed.)
    Recent theoretical and computational progress has led to unprecedented understanding of symmetry-breaking instabilities in 2D dynamic fracture. At the heart of this progress resides the identification of two intrinsic, near crack tip length scales — a nonlinear elastic length scale ℓ and a dissipation length scale ξ — that do not exist in Linear Elastic Fracture Mechanics (LEFM), the classical theory of cracks. In particular, it has been shown that at a propagation velocity v of about 90% of the shear wave-speed, cracks in 2D brittle materials undergo an oscillatory instability whose wavelength varies linearly with ℓ, and at larger loading levels (corresponding to yet higher propagation velocities), a tip-splitting instability emerges, both in agreements with experiments. In this paper, using phase-field models of brittle fracture, we demonstrate the following properties of the oscillatory instability: (i) It exists also in the absence of near-tip elastic nonlinearity, i.e. in the limit ℓ→0, with a wavelength determined by the dissipation length scale ξ. This result shows that the instability crucially depends on the existence of an intrinsic length scale associated with the breakdown of linear elasticity near crack tips, independently of whether the latter is related to nonlinear elasticity or to dissipation. (ii) It is a supercritical Hopf bifurcation, featuring a vanishing oscillations amplitude at onset. (iii) It is largely independent of the phenomenological forms of the degradation functions assumed in the phase-field framework to describe the cohesive zone, and of the velocity-dependence of the fracture energy Γ(v) that is controlled by the dissipation time scale in the Ginzburg-Landau-type evolution equation for the phase-field. These results substantiate the universal nature of the oscillatory instability in 2D. In addition, we provide evidence indicating that the tip-splitting instability is controlled by the limiting rate of elastic energy transport inside the crack tip region. The latter is sensitive to the wave-speed inside the dissipation zone, which can be systematically varied within the phase-field approach. Finally, we describe in detail the numerical implementation scheme of the employed phase-field fracture approach, allowing its application in a broad range of materials failure problems. 
    more » « less
  5. A comparative study is presented to solve the inverse problem in elasticity for the shear modulus (stiffness) distribution utilizing two constitutive equations: (1) linear elasticity assuming small strain theory, and (2) finite elasticity with a hyperelastic neo-Hookean material model. Assuming that a material undergoes large deformations and material nonlinearity is assumed negligible, the inverse solution using (2) is anticipated to yield better results than (1). Given the fact that solving a linear elastic model is significantly faster than a nonlinear model and more robust numerically, we posed the following question: How accurately could we map the shear modulus distribution with a linear elastic model using small strain theory for a specimen undergoing large deformations? To this end, experimental displacement data of a silicone composite sample containing two stiff inclusions of different sizes under uniaxial displacement controlled extension were acquired using a digital image correlation system. The silicone based composite was modeled both as a linear elastic solid under infinitesimal strains and as a neo-Hookean hyperelastic solid that takes into account geometrically nonlinear finite deformations. We observed that the mapped shear modulus contrast, determined by solving an inverse problem, between inclusion and background was higher for the linear elastic model as compared to that of the hyperelastic one. A similar trend was observed for simulated experiments, where synthetically computed displacement data were produced and the inverse problem solved using both, the linear elastic model and the neo-Hookean material model. In addition, it was observed that the inverse problem solution was inclusion size-sensitive. Consequently, an 1-D model was introduced to broaden our understanding of this issue. This 1-D analysis revealed that by using a linear elastic approach, the overestimation of the shear modulus contrast between inclusion and background increases with the increase of external loads and target shear modulus contrast. Finally, this investigation provides valuable information on the validity of the assumption for utilizing linear elasticity in solving inverse problems for the spatial distribution of shear modulus associated with soft solids undergoing large deformations. Thus, this work could be of importance to characterize mechanical property variations of polymer based materials such as rubbers or in elasticity imaging of tissues for pathology. 
    more » « less