Incorporating physically-based microstructures in materials modeling: Bridging phase field and crystal plasticity frameworks

Abstract

The mechanical properties of materials systems are highly influenced by various features at the microstructural level. The ability to capture these heterogeneities and incorporate them into continuum-scale frameworks of the deformation behavior is considered a key step in the development of complex non-local models of failure. In this study, we present a modeling framework that incorporates physically-based realizations of polycrystalline aggregates from a phase field (PF) model into a crystal plasticity finite element (CP-FE) framework. Simulated annealing via the PF model yields ensembles of materials microstructures with various grain sizes and shapes. With the aid of a novel FE meshing technique, FE discretizations of these microstructures are generated, where several key features, such as conformity to interfaces, and triple junction angles, are preserved. The discretizations are then used in the CP-FE framework to simulate the mechanical response of polycrystalline α-iron. It is shown that the conformal discretization across interfaces reduces artificial stress localization commonly observed in non-conformal FE discretizations. The work presented herein is a first step towards incorporating physically-based microstructures in lieu of the overly simplified representations that are commonly used. In broader terms, the proposed framework provides future avenues to explore bridging models of materials processes, e.g. additive manufacturing and microstructure evolution of multi-phase multi-component systems, into continuum-scale frameworks of the mechanical properties.

Reduced partitioning of plastic strain for strong and yet ductile precipitate-strengthened alloys

“…8a . The cuboidal precipitate was reproduced by 163,400 hexahedral finite elements, and the smooth interfaces between two phases were constructed while maintaining precipitate volume fraction of 50.6% 59 . Taking an orthogonal ( x , y , z ) coordinate system with origin at a top corner of the cube and axes normal to its faces, as shown in the same figure, the compression was applied with a rate of 1 × 10 −3 s −1 by displacing the y = 0 face along - y direction while imposing u x = 0 on the x = 0 face, u z = 0 on the z = 0 face and u y = 0 on the - y face.…”

Progress in modeling solidification in molten salt coolants

“…The initial concentration is determined by a Voronoi tessellation method [69] following the eutectic composition for each of the phases. With the generated concentration field a 2D phase-field simulation is performed in the base plane (z * = 0) [70] in order homogenize the phase distribution. This allows to reduce the sensitivity of the system to the initial conditions.…”

Lower length scale informed improvements to Bison U-Pu-Zr fuel swelling model

“…Each grain within the polycrystalline structure is a unique mesh block and the interface between the grains is sharp. The conformal mesh of the grain structure is generated via several steps using the Sculpt and Cubit software tools from Sandia National Laboratory [17]. First, a phase field simulation of grain growth generates a realistic diffuse-interface equiaxed structure.…”