We present a phase-field (PF) model to simulate the microstructure evolution occurring in polycrystalline materials with a variation in the intra-granular dislocation density. The model accounts for two mechanisms that lead to the grain boundary migration: the driving force due to capillarity and that due to the stored energy arising from a spatially varying dislocation density. In addition to the order parameters that distinguish regions occupied by different grains, we introduce dislocation density fields that describe spatial variation of the dislocation density. We assume that the dislocation density decays as a function of the distance the grain boundary has migrated. To demonstrate and parameterize the model, we simulate microstructure evolution in two dimensions, for which the initial microstructure is based on real-time experimental data. Additionally, we applied the model to study the effect of a cyclic heat treatment (CHT) on the microstructure evolution. Specifically, we simulated stored-energy-driven grain growth during three thermal cycles, as well as grain growth without stored energy that serves as a baseline for comparison. We showed that the microstructure evolution proceeded much faster when the stored energy was considered. A non-self-similar evolution was observed in this case, while a nearly self-similar evolution was found when the microstructure evolution is driven solely by capillarity. These results suggest a possible mechanism for the initiation of abnormal grain growth during CHT. Finally, we demonstrate an integrated experimental-computational workflow that utilizes the experimental measurements to inform the PF model and its parameterization, which provides a foundation for the development of future simulation tools capable of quantitative prediction of microstructure evolution during non-isothermal heat treatment.
- Award ID(s):
- 1729336
- NSF-PAR ID:
- 10104817
- Date Published:
- Journal Name:
- Proceedings of NUMIFORM 2019: The 13th International Conference on Numerical Methods in Industrial Forming Processes
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
Scale-free phase-field approach and corresponding finite element method simulations for multivariant martensitic phase transformation from cubic Si I to tetragonal Si II in a polycrystalline aggregate are presented. Important features of the model are large and very anisotropic transformation strain tensor 𝜺𝑡 = {0.1753; 0.1753; −0.447} and stress-tensor dependent athermal dissipative threshold for transformation, which produce essential challenges for computations. 3D polycrystals with stochastically oriented grains are subjected to uniaxial strain- and stress-controlled loadings under periodic boundary conditions and zero averaged lateral strains. Coupled evolution of discrete martensitic microstructure, volume fractions of martensitic variants and Si II, stress and transformation strain tensors, and texture are presented and analyzed. Macroscopic variables effectively representing multivariant transformational behavior are introduced. Macroscopic stress–strain and transformational behavior for 55 and 910 grains are close. Large transformation strains and grain boundaries lead to huge internal stresses of tens GPa, which affect microstructure evolution and macroscopic behavior. In contrast to a single crystal, the local mechanical instabilities due to phase transformation and negative local tangent modulus are stabilized at the macroscale by arresting/slowing the growth of Si II regions by the grain boundaries. This leads to increasing stress during transformation. The developed methodology can be used for studying similar phase transformations with large transformation strains and for further development by including plastic strain and strain-induced transformations.more » « less
-
Abstract The phase-field (PF) method is a physics-based computational approach for simulating interfacial morphology. It has been used to model powder melting, rapid solidification, and grain structure evolution in metal additive manufacturing (AM). However, traditional direct numerical simulation (DNS) of the PF method is computationally expensive due to sufficiently small mesh size. Here, a physics-embedded graph network (PEGN) is proposed to leverage an elegant graph representation of the grain structure and embed the classic PF theory into the graph network. By reformulating the classic PF problem as an unsupervised machine learning task on a graph network, PEGN efficiently solves temperature field, liquid/solid phase fraction, and grain orientation variables to minimize a physics-based loss/energy function. The approach is at least 50 times faster than DNS in both CPU and GPU implementation while still capturing key physical features. Hence, PEGN allows to simulate large-scale multi-layer and multi-track AM build effectively.
-
Abstract The present work utilizes Orientation Imaging Microscopy and Finite Element Modelling to analyse microstructure evolution in grains near defects during plane strain indentation of direct metal laser sintered Inconel 718. Defects are inevitably produced during printing of metals and they degrade the mechanical behaviour of parent components. Understanding microstructure evolution of grains present near defects can help create better predictive models of mechanical behaviour of components resulting from additive manufacturing. In this work, an ex-situ study of microstructure evolution during plane strain indentation of DMLS Inconel 718 specimens is performed. Regions that lie near volumetric porosity defects were studied. Grain Orientation Spread was utilized as a metric to quantify intra-granular deformation. It was seen that microstructure evolution of grains near defects is enhanced due to strain concentrations whereby they exhibit larger orientation spread after plastic deformation. Finite Element Analysis was used to simulate the plane strain indentation test on the specimen in which, porosity defects and roughness textures similar to those seen in the as-received specimen were programmed using the python scripting interface of Abaqus. Results from finite element analysis were compared with insights from microstructure analysis to describe evolution of microstructure during deformation near defects.
-
Abstract The influence of grain size distribution on ductile intergranular crack growth resistance is investigated using full-field microstructure-based finite element calculations and a simpler model based on discrete unit events and graph search. The finite element calculations are carried out for a plane strain slice with planar grains subjected to mode I small-scale yielding conditions. The finite element formulation accounts for finite deformations, and the constitutive relation models the loss of stress carrying capacity due to progressive void nucleation, growth, and coalescence. The discrete unit events are characterized by a set of finite element calculations for crack growth at a single-grain boundary junction. A directed graph of the connectivity of grain boundary junctions and the distances between them is used to create a directed graph in J-resistance space. For a specified grain boundary distribution, this enables crack growth resistance curves to be calculated for all possible crack paths. Crack growth resistance curves are calculated based on various path choice criteria and compared with the results of full-field finite element calculations of the initial boundary value problem. The effect of unimodal and bimodal grain size distributions on intergranular crack growth is considered. It is found that a significant increase in crack growth resistance is obtained if the difference in grain sizes in the bimodal grain size distribution is sufficiently large.more » « less