Microstructure-sensitive uncertainty quantification for crystal plasticity finite element constitutive models using stochastic collocation methods

TL;DR

This paper applies stochastic collocation methods to quantify uncertainty in crystal plasticity finite element models, analyzing the robustness of different constitutive models across various crystal structures and sensitivities.

Contribution

It introduces a stochastic collocation approach to rigorously quantify uncertainty and sensitivity in CPFEM constitutive models for different crystal structures.

Findings

Quantified uncertainty in stress-strain predictions for various models.

Identified key parameters influencing initial yield behavior.

Provided insights for robust model calibration.

Abstract

Uncertainty quantification (UQ) plays a major role in verification and validation of computational engineering models and simulations, and establishes trust in the predictive capability of computational models. In the materials science and engineering context, where the process-structure-property-performance linkage is well known to be the only road mapping from manufacturing to engineering performance, numerous integrated computational materials engineering (ICME) models have been developed across a wide spectrum of length-scales and time-scales to relieve the burden of resource-intensive experiments. Within the structure-property linkage, crystal plasticity finite element method (CPFEM) models have been widely used since they are one of a few ICME toolboxes that allow numerical predictions, providing the bridge from microstructure to properties and performances. Several constitutive models have been proposed to capture the mechanics and plasticity behavior of materials. While some UQ studies have been performed, the robustness and uncertainty of these constitutive models have not been rigorously established. In this work, we apply a stochastic collocation (SC) method to quantify the uncertainty of the three most commonly used constitutive models in CPFEM, namely phenomenological models (with and without twinning), and dislocation-density-based constitutive models, for three different types of crystal structures, namely face-centered cubic (fcc) copper (Cu), body-centered cubic (bcc) tungsten (W), and hexagonal close packing (hcp) magnesium (Mg). Our numerical results not only quantify the uncertainty of these constitutive models in stress-strain curves, but also analyze the global sensitivity of the underlying constitutive parameters with respect to the initial yield behavior, which may be helpful for robust constitutive model calibration works in the future.