Direct Numerical Simulations in Solid Mechanics for Quantifying the Macroscale Effects of Microstructure and Material Model-Form Error

Two fundamental approximations in macroscale solid-mechanics modeling are (1) the assumption of scale separation in homogenization theory and (2) the use of a macroscopic plasticity material model that represents, in a mean sense, the multitude of inelastic processes occurring at the microscale. With the goal of quantifying the errors induced by these approximations on engineering quantities of interest, we perform a set of direct numerical simulations (DNS) in which polycrystalline microstructures are embedded throughout a macroscale structure. The largest simulations model over 50,000 grains. The microstructure is idealized using a randomly close-packed Voronoi tessellation in which each polyhedral Voronoi cell represents a grain. An face centered cubic crystal-plasticity model is used to model the mechanical response of each grain. The overall grain structure is equiaxed, and each grain is randomly oriented with no overall texture. The detailed results from the DNS simulations are compared to results obtained from conventional macroscale simulations that use homogeneous isotropic plasticity models. The macroscale plasticity models are calibrated using a representative volume element of the idealized microstructure. Ultimately, we envision that DNS modeling will be used to gain new insights into the mechanics of material deformation and failure.

[1]  Nikolas Provatas,et al.  Phase-Field Methods in Materials Science and Engineering , 2010 .

[2]  D. Jeulin,et al.  Determination of the size of the representative volume element for random composites: statistical and numerical approach , 2003 .

[3]  H. Dumontet.  Study of a boundary layer problem in elastic composite materials , 1986 .

[4]  V. G. Kouznetsova,et al.  Multi-scale computational homogenization: Trends and challenges , 2010, J. Comput. Appl. Math..

[5]  Jeffrey S. Vetter.  Contemporary High Performance Computing: From Petascale toward Exascale , 2013 .

[6]  A. Maniatty,et al.  Effect of 3D grain structure representation in polycrystal simulations , 2013 .

[7]  Q. Zheng,et al.  A Quantitative study of minimum sizes of representative volume elements of cubic polycrystals—numerical experiments , 2002 .

[8]  A. Bensoussan,et al.  Asymptotic analysis for periodic structures , 1979 .

[9]  Sia Nemat-Nasser,et al.  Averaging theorems in finite deformation plasticity , 1999 .

[10]  A. A. Gusev.  Representative volume element size for elastic composites: A numerical study , 1997 .

[11]  J. Tinsley Oden,et al.  Multi-scale goal-oriented adaptive modeling of random heterogeneous materials , 2006 .

[12]  Martin Ostoja-Starzewski,et al.  Random field models of heterogeneous materials , 1998 .

[13]  Michael A. Groeber,et al.  Consistent representations of and conversions between 3D rotations , 2015 .

[14]  Jin Qian,et al.  Quality Improvement of Non-manifold Hexahedral Meshes for Critical Feature Determination of Microstructure Materials , 2009, IMR.

[15]  W. Hosford.  A Generalized Isotropic Yield Criterion , 1972 .

[16]  Christian Huet,et al.  Application of variational concepts to size effects in elastic heterogeneous bodies , 1990 .

[17]  Vincent Monchiet,et al.  A micromechanics-based approach for the derivation of constitutive elastic coefficients of strain-gradient media , 2012 .

[18]  P. Wriggers,et al.  An adaptive multiscale resolution strategy for the finite deformation analysis of microheterogeneous structures , 2011 .

[19]  Brad Lee Boyce,et al.  Crystal plasticity simulations of microstructure-induced uncertainty in strain concentration near voids in brass , 2015 .

[20]  Matthew P. Miller,et al.  Influence of slip system hardening assumptions on modeling stress dependence of work hardening , 1997 .

[21]  G. Proust,et al.  Generalised Voronoi tessellation for generating microstructural finite element models with controllable grain‐size distributions and grain aspect ratios , 2015 .

[22]  J. Oden,et al.  Estimation of local modeling error and goal-oriented adaptive modeling of heterogeneous materials: Part II: a computational environment for adaptive modeling of heterogeneous elastic solids , 2001 .

[23]  Wei Chen,et al.  Statistical volume element method for predicting microstructure–constitutive property relations , 2008 .

[24]  Z. Ren,et al.  Effects of grain sizes, shapes, and distribution on minimum sizes of representative volume elements of cubic polycrystals , 2004 .

[25]  G. Cailletaud,et al.  Hardening description for FCC materials under complex loading paths , 2009 .

[26]  M. Diehl,et al.  A spectral method solution to crystal elasto-viscoplasticity at finite strains , 2013 .

[27]  Karel Matouš,et al.  Finite element formulation for modelling large deformations in elasto‐viscoplastic polycrystals , 2004 .

[28]  P. Donato,et al.  An introduction to homogenization , 2000 .

[29]  Steven J. Plimpton,et al.  Crossing the Mesoscale No-Man's Land via Parallel Kinetic Monte Carlo , 2009 .

[30]  G. Cailletaud,et al.  Modeling of latent hardening produced by complex loading paths in FCC alloys , 2013 .

[31]  H. Ledbetter.  Monocrystal-Polycrystal Elastic Constants of a Stainless Steel , 1984, September 16, 1984.

[32]  Dierk Raabe,et al.  Crystal Plasticity Finite Element Methods: in Materials Science and Engineering , 2010 .

[33]  P. Dawson,et al.  A time integration algorithm for elasto‐viscoplastic cubic crystals applied to modelling polycrystalline deformation , 1992 .

[34]  R. Macpherson,et al.  A more accurate two-dimensional grain growth algorithm , 2010 .

[35]  D. Agard,et al.  Microtubule nucleation by γ-tubulin complexes , 2011, Nature Reviews Molecular Cell Biology.

[36]  Fpt Frank Baaijens,et al.  An approach to micro-macro modeling of heterogeneous materials , 2001 .

[37]  William E. Frazier,et al.  Metal Additive Manufacturing: A Review , 2014, Journal of Materials Engineering and Performance.

[38]  W. Drugan,et al.  A micromechanics-based nonlocal constitutive equation and estimates of representative volume element size for elastic composites , 1996 .

[39]  Anthony D. Rollett,et al.  Three-dimensional digital approximations of grain boundary networks in polycrystals , 2014 .

[40]  Surface effects in composite materials: Two simple examples , 1987 .

[41]  Chandrajit L Bajaj,et al.  An Automatic 3D Mesh Generation Method for Domains with Multiple Materials. , 2010, Computer methods in applied mechanics and engineering.

[42]  D. Rixen,et al.  FETI‐DP: a dual–primal unified FETI method—part I: A faster alternative to the two‐level FETI method , 2001 .

[43]  Elizabeth A. Holm,et al.  The computer simulation of microstructural evolution , 2001 .

[44]  V. Zhikov,et al.  Homogenization of Differential Operators and Integral Functionals , 1994 .

[45]  Christopher R. Weinberger,et al.  Direct numerical simulations in solid mechanics for understanding the macroscale effects of microscale material variability , 2015 .

[46]  R. Quey,et al.  Large-scale 3D random polycrystals for the finite element method: Generation, meshing and remeshing , 2011 .

[47]  James Arvo,et al.  Fast Random Rotation matrices , 1992, Graphics Gems III.

[48]  Hojun Lim,et al.  Investigation of grain-scale microstructural variability in tantalum using crystal plasticity-finite element simulations , 2016 .