skip to main content

Title: Numerical modeling of Ti-6Al-4V microstructure evolution for thermomechanical process control
Microstructure evolution modeling using finite element crystal plasticity (FECP), Monte- Carlo (MC), and phase field (PF) methods are being used to simulate microstructure evolution in Ti-6Al-4V under thermomechanical loading conditions. FECP is used to simulate deformation induced evolution of the microstructure and compute heterogeneous stored energy providing additional source of energy to MC and PF models. The MC grain growth model, calibrated using literature and experimental data, is used to simulate α+𝛽 grain growth. A multi-phase field, augmented with crystallographic symmetry and orientation relationship between α-𝛽, is employed to model simultaneous evolution and growth of all twelve α-variants in 3D. The influence of transformation and coherency strain energy on α-variant selection is studied by coupling the model with the Khachaturyan-Shatalov formalism for elastic strain calculation. This FECP/MC/PF suite will be able to simulate evolution of grains in the microstructure and within individual 𝛽- grains during typical thermomechanical processing conditions.
; ; ;
Award ID(s):
Publication Date:
Journal Name:
Proceedings of NUMIFORM 2019: The 13th International Conference on Numerical Methods in Industrial Forming Processes
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 amore »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.« less
  2. Among metal additive manufacturing technologies, additive friction stir deposition stands out for its ability to create freeform and fully-dense structures without melting and solidification. Here, we employ a comparative approach to investigate the process-microstructure linkages in additive friction stir deposition, utilizing two materials with distinct thermomechanical behavior—an Al-Mg-Si alloy and Cu—both of which are challenging to print using beam-based additive processes. The deposited Al-Mg-Si is shown to exhibit a relatively homogeneous microstructure with extensive subgrain formation and a strong shear texture, whereas the deposited Cu is characterized by a wide distribution of grain sizes and a weaker shear texture. Wemore »show evidence that the microstructure in Al-Mg-Si primarily evolves by continuous dynamic recrystallization, including geometric dynamic recrystallization and progressive lattice rotation, while the heterogeneous microstructure of Cu results from discontinuous recrystallization during both deposition and cooling. In Al-Mg-Si, the continuous recrystallization progresses with an increase of the applied strain, which correlates with the ratio between the tool rotation rate and travel velocity. Conversely, the microstructure evolution in Cu is found to be less dependent on , instead varying more with changes to . This difference originates from the absence of Cu rotation in the deposition zone, which reduces the influence of tool rotation on strain development. We attribute the distinct process-microstructure linkages and the underlying mechanisms between Al-Mg-Si and Cu to their differences in intrinsic thermomechanical properties and interactions with the tool head.« less
  3. Solid solutions of Mg 2 Si and Mg 2 Sn are promising thermoelectric materials owing to their high thermoelectric figures-of-merit and non-toxicity, but they may undergo phase separation under thermal cycling due to the presence of miscibility gaps, implying that the thermoelectric properties could be significantly degraded during thermoelectric device operation. Herein, this study investigates the strain-induced suppression of the miscibility gap in solid solutions of Mg 2 Si and Mg 2 Sn. Separately prepared Mg 2 Si and Mg 2 Sn powders were made into (Mg 2 Si) 0.7 (Mg 2 Sn) 0.3 mixtures using a high energy ball-millingmore »method followed by spark plasma sintering. Afterwards, the phase evolution of the mixtures, depending on thermal annealing and mixing conditions, was studied experimentally and theoretically. Transmission electron microscopy and X-ray diffraction results show that, despite the presence of a miscibility gap in the pseudo-binary phase diagram, the initial mixture of Mg 2 Si and Mg 2 Sn evolved towards a solid solution state after annealing for 3 hours at 720 °C. Thermodynamic analysis as well as phase-field microstructure simulations show that the strain energy due to the coherent spinodal effect suppresses the chemical spinodal entirely and prevents phase separation. This strategy to suppress the miscibility gap induced by lattice strain through non-equilibrium processing can benefit the thermoelectric figure-of-merit by maximizing phonon alloy scattering. Furthermore, stable solid solutions by engineering phase diagrams have the potential to facilitate the reliable long term operation of thermoelectric generators under continuous thermal loads.« less
  4. Abstract. Snow in the environment acts as a host to rich chemistry and provides a matrix for physical exchange of contaminants within the ecosystem. The goal of this review is to summarise the current state of knowledge of physical processes and chemical reactivity in surface snow with relevance to polar regions. It focuses on a description of impurities in distinct compartments present in surface snow, such as snow crystals, grain boundaries, crystal surfaces, and liquid parts. It emphasises the microscopic description of the ice surface and its link with the environment. Distinct differences between the disordered air–ice interface, often termedmore »quasi-liquid layer, and a liquid phase are highlighted. The reactivity in these different compartments of surface snow is discussed using many experimental studies, simulations, and selected snow models from the molecular to the macro-scale.

    Although new experimental techniques have extended our knowledge of the surface properties of ice and their impact on some single reactions and processes, others occurring on, at or within snow grains remain unquantified. The presence of liquid or liquid-like compartments either due to the formation of brine or disorder at surfaces of snow crystals below the freezing point may strongly modify reaction rates. Therefore, future experiments should include a detailed characterisation of the surface properties of the ice matrices. A further point that remains largely unresolved is the distribution of impurities between the different domains of the condensed phase inside the snowpack, i.e. in the bulk solid, in liquid at the surface or trapped in confined pockets within or between grains, or at the surface. While surface-sensitive laboratory techniques may in the future help to resolve this point for equilibrium conditions, additional uncertainty for the environmental snowpack may be caused by the highly dynamic nature of the snowpack due to the fast metamorphism occurring under certain environmental conditions.

    Due to these gaps in knowledge the first snow chemistry models have attempted to reproduce certain processes like the long-term incorporation of volatile compounds in snow and firn or the release of reactive species from the snowpack. Although so far none of the models offers a coupled approach of physical and chemical processes or a detailed representation of the different compartments, they have successfully been used to reproduce some field experiments. A fully coupled snow chemistry and physics model remains to be developed.

    « less
  5. Abstract The rheology of the upper mantle impacts a variety of geodynamic processes, including postseismic deformation following great earthquakes and post-glacial rebound. The deformation of upper mantle rocks is controlled by the rheology of olivine, the most abundant upper mantle mineral. The mechanical properties of olivine at steady state are well constrained. However, the physical mechanism underlying transient creep, an evolutionary, hardening phase converging to steady state asymptotically, is still poorly understood. Here, we constrain a constitutive framework that captures transient creep and steady state creep consistently using the mechanical data from laboratory experiments on natural dunites containing at leastmore »94% olivine under both hydrous and anhydrous conditions. The constitutive framework represents a Burgers assembly with a thermally activated nonlinear stress-versus-strain-rate relationship for the dashpots. Work hardening is obtained by the evolution of a state variable that represents internal stress. We determine the flow law parameters for dunites using a Markov chain Monte Carlo method. We find the activation energy $$430\pm 20$$ 430 ± 20   and $$250\pm 10$$ 250 ± 10  kJ/mol for dry and wet conditions, respectively, and the stress exponent $$2.0\pm 0.1$$ 2.0 ± 0.1 for both the dry and wet cases for transient creep, consistently lower than those of steady-state creep, suggesting a separate physical mechanism. For wet dunites in the grain-boundary sliding regime, the grain-size dependence is similar for transient creep and steady-state creep. The lower activation energy of transient creep could be due to a higher jog density of the corresponding soft-slip system. More experimental data are required to estimate the activation volume and water content exponent of transient creep. The constitutive relation used and its associated flow law parameters provide useful constraints for geodynamics applications. Graphical Abstract« less