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: Multiscale Modeling of Composite Materials under Volumetric and Interfacial Damage: Achieving Adaptive Model Order Reduction
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
Award ID(s):
2114822
PAR ID:
10396899
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
AIAA SCITECH 2023 Forum
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Metal alloys frequently contain distributions of second-phase particles that deleteriously affect the material behavior by acting as sites for void nucleation. These distributions are often extremely complex and processing can induce high levels of anisotropy. The particle length-scale precludes high-fidelity microstructure modeling in macroscale simulations, so computational homogenization methods are often employed. These, however, involve simplifying assumptions to make the problem tractable and many rely on periodic microstructures. Here we propose a methodology to bridge the gap between realistic microstructures composed of anisotropic, spatially varying second-phase void morphologies and idealized periodic microstructures with roughly equivalent mechanical responses. We create a high-throughput, parametric study to investigate 96 unique bridging methods. We apply our proposed solution to a rolled AZ31B magnesium alloy, for which we have a rich dataset of microstructure morphology and mechanical behavior. Our methodology converts aµ-CT scan of the realistic microstructure to idealized periodic unit cell microstructures that are specific to the loading orientation. We recreate the unit cells for each parameter set in a commercial finite element software, subject them to macroscopic uniaxial loading conditions, and compare our results to the datasets for the various loading orientations. We find that certain combinations of our parameters capture the overall stress–strain response, including anisotropy effects, with some degree of success. The effect of different parameter options are explored in detail and we find that excluding certain particle populations from the analysis can give improved results. 
    more » « less
  2. 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
  3. 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
  4. 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
  5. 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