skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Data-Driven Calibration of Multifidelity Multiscale Fracture Models Via Latent Map Gaussian Process
Abstract Fracture modeling of metallic alloys with microscopic pores relies on multiscale damage simulations which typically ignore the manufacturing-induced spatial variabilities in porosity. This simplification is made because of the prohibitive computational expenses of explicitly modeling spatially varying microstructures in a macroscopic part. To address this challenge and open the doors for the fracture-aware design of multiscale materials, we propose a data-driven framework that integrates a mechanistic reduced-order model (ROM) with a calibration scheme based on random processes. Our ROM drastically accelerates direct numerical simulations (DNS) by using a stabilized damage algorithm and systematically reducing the degrees of freedom via clustering. Since clustering affects local strain fields and hence the fracture response, we calibrate the ROM by constructing a multifidelity random process based on latent map Gaussian processes (LMGPs). In particular, we use LMGPs to calibrate the damage parameters of an ROM as a function of microstructure and clustering (i.e., fidelity) level such that the ROM faithfully surrogates DNS. We demonstrate the application of our framework in predicting the damage behavior of a multiscale metallic component with spatially varying porosity. Our results indicate that microstructural porosity can significantly affect the performance of macro-components and hence must be considered in the design process.  more » « less
Award ID(s):
2103708
PAR ID:
10442718
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Mechanical Design
Volume:
145
Issue:
1
ISSN:
1050-0472
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Predicting the fracture behavior of macroscale components containing microscopic porosity relies on multiscale damage models which typically ignore the manufacturing-induced spatial variabilities in porosity. This simplification is made due to the prohibitive computational costs associated with explicitly modeling spatially varying microstructures in a macroscopic component. To address this challenge, we propose a data-driven framework that integrates a mechanistic reduced-order model (ROM) with a calibration scheme based on latent map Gaussian processes (LMGPs). Our ROM drastically accelerates direct numerical simulations (DNS) by using a stabilized damage algorithm and systematically reducing the degrees of freedom via clustering. Since clustering affects local strain fields and hence the fracture response, we construct a multi-fidelity LMGP to inversely estimate the damage parameters of an ROM as a function of microstructure and clustering level such that the ROM faithfully surrogates DNS. We demonstrate the application of our framework in predicting the damage behavior of a multiscale metallic component with spatially varying porosity. 
    more » « less
  2. Abstract Aluminum alloys are increasingly utilized as lightweight materials in the automobile industry due to their superior capability in withstanding high mechanical loads. A significant challenge impeding the large-scale use of these alloys in high-performance applications is the presence of manufacturing-induced, spatially varying porosity defects. In order to understand the impacts of these defects on the macro-mechanical properties of cast alloys, multiscale simulations are often required. In this paper, we introduce a computationally efficient reduced-order multiscale framework to simulate the behavior of metallic components containing process-induced porosity under irreversible nonlinear deformations. In our approach, we start with a data compression scheme that significantly reduces the number of unknown macroscale and microscale variables by agglomerating close-by finite element nodes into a limited number of clusters. Then, we use deflation methods to project these variables into a lower-dimensional space where the material’s elastoplastic behaviors are approximated. Finally, we solve for the unknown variables and map them back to the original, high-dimensional space. We call our method deflated clustering analysis and by comparing it to direct numerical simulations we demonstrate that it accurately captures macroscale deformations and microscopic effective responses. To illustrate the effect of microscale pores on the macroscopic response of a cast component, we conduct multi-scale simulations with spatially varying local heterogeneities that are modeled with a microstructure characterization and reconstruction algorithm. 
    more » « less
  3. Abstract The vagina undergoes large finite deformations and has complex geometry and microstructure, resulting in material and geometric nonlinearities, complicated boundary conditions, and nonhomogeneities within finite element (FE) simulations. These nonlinearities pose a significant challenge for numerical solvers, increasing the computational time by several orders of magnitude. Simplifying assumptions can reduce the computational time significantly, but this usually comes at the expense of simulation accuracy. This study proposed the use of reduced order modeling (ROM) techniques to capture experimentally measured displacement fields of rat vaginal tissue during inflation testing in order to attain both the accuracy of higher‐fidelity models and the speed of simpler simulations. The proper orthogonal decomposition (POD) method was used to extract the significant information from FE simulations generated by varying the luminal pressure and the parameters that introduce the anisotropy in the selected constitutive model. A new data‐driven (DD) variational multiscale (VMS) ROM framework was extended to obtain the displacement fields of rat vaginal tissue under pressure. For comparison purposes, we also investigated the classical Galerkin ROM (G‐ROM). In our numerical study, both the G‐ROM and the DD‐VMS‐ROM decreased the FE computational cost by orders of magnitude without a significant decrease in numerical accuracy. Furthermore, the DD‐VMS‐ROM improved the G‐ROM accuracy at a modest computational overhead. Our numerical investigation showed that ROM has the potential to provide efficient and accurate computational tools to describe vaginal deformations, with the ultimate goal of improving maternal health. 
    more » « less
  4. Abstract Fault‐damage zones comprise multiscale fracture networks that may slip dynamically and interact with the main fault during earthquake rupture. Using 3D dynamic rupture simulations and scale‐dependent fracture energy, we examine dynamic interactions of more than 800 intersecting multiscale fractures surrounding a listric fault, emulating a major listric fault and its damage zone. We investigate 10 distinct orientations of maximum horizontal stress, probing the conditions necessary for sustained slip within the fracture network or activating the main fault. Additionally, we assess the feasibility of nucleating dynamic rupture earthquake cascades from a distant fracture and investigate the sensitivity of fracture network cascading rupture to the effective normal stress level. We model either pure cascades or main fault rupture with limited off‐fault slip. We find that cascading ruptures within the fracture network are dynamically feasible under certain conditions, including: (a) the fracture energy scales with fracture and fault size, (b) favorable relative pre‐stress of fractures within the ambient stress field, and (c) close proximity of fractures. We find that cascading rupture within the fracture network discourages rupture on the main fault. Our simulations suggest that fractures with favorable relative pre‐stress, embedded within a fault damage zone, may lead to cascading earthquake rupture that shadows main fault slip. We find that such off‐fault events may reach moment magnitudes up toMw ≈ 5.5, comparable to magnitudes that can be otherwise hosted by the main fault. Our findings offer insights into physical processes governing cascading earthquake dynamic rupture within multiscale fracture networks. 
    more » « less
  5. Abstract Double-network (DN) hydrogels, consisting of two contrasting and interpenetrating polymer networks, are considered as perhaps the toughest soft-wet materials. Current knowledge of DN gels from synthesis methods to toughening mechanisms almost exclusively comes from chemically-linked DN hydrogels by experiments. Molecular modeling and simulations of inhomogeneous DN structure in hydrogels have proved to be extremely challenging. Herein, we developed a new multiscale simulation platform to computationally investigate the early fracture of physically-chemically linked agar/polyacrylamide (agar/PAM) DN hydrogels at a long timescale. A “random walk reactive polymerization” (RWRP) was developed to mimic a radical polymerization process, which enables to construct a physically-chemically linked agar/PAM DN hydrogel from monomers, while conventional and steered MD simulations were conducted to examine the structural-dependent energy dissipation and fracture behaviors at the relax and deformation states. Collective simulation results revealed that energy dissipation of agar/PAM hydrogels was attributed to a combination of the pulling out of agar chains from the DNs, the disruption of massive hydrogen bonds between and within DN structures, and the strong association of water molecules with both networks, thus explaining a different mechanical enhancement of agar/PAM hydrogels. This computational work provided atomic details of network structure, dynamics, solvation, and interactions of a hybrid DN hydrogel, and a different structural-dependent energy dissipation mode and fracture behavior of a hybrid DN hydrogel, which help to design tough hydrogels with new network structures and efficient energy dissipation modes. Additionally, the RWRP algorithm can be generally applied to construct the radical polymerization-produced hydrogels, elastomers, and polymers. 
    more » « less