Abstract Additively manufactured (AM) metallic materials often comprise as-printed dislocation cells inside grains. These dislocation cells can give rise to substantial microscale internal stresses in both initial undeformed and plastically deformed samples, thereby affecting the mechanical properties of AM metallic materials. Here we develop models of microscale internal stresses in AM stainless steel by focusing on their back stress components. Three sources of microscale back stresses are considered, including the printing and deformation-induced back stresses associated with as-printed dislocation cells as well as the deformation-induced back stresses associated with grain boundaries. We use a three-dimensional discrete dislocation dynamics model to demonstrate the manifestation of printing-induced back stresses. We adopt a dislocation pile-up model to evaluate the deformation-induced back stresses associated with as-printed dislocation cells. The extracted back stress relation from the pile-up model is incorporated into a crystal plasticity model that accounts for the other two sources of back stresses as well. The crystal plasticity finite element simulation results agree with the experimentally measured tension-compression asymmetry and macroscopic back stress, the latter of which represents the effective resultant of microscale back stresses of different origins. Our results provide an in-depth understanding of the origins and evolution of microscale internal stresses in AM metallic materials.
more »
« less
Phase field theory for fracture at large strains including surface stresses
Phase field theory for fracture is developed at large strains with an emphasis on a correct introduction of surface stresses. This is achieved by multiplying the cohesion and gradient energies by the local ratio of the crack surface areas in the deformed and undeformed configurations and with the gradient energy in terms of the gradient of the order parameter in the reference configuration. This results in an expression for the surface stresses which is consistent with the sharp surface approach. Namely, the structural part of the Cauchy surface stress represents an isotropic biaxial tension, with the magnitude of a force per unit length equal to the surface energy. The surface stresses are a result of the geometric nonlinearities, even when strains are infinitesimal. They make multiple contributions to the Ginzburg-Landau equation for damage evolution, both in the deformed and undeformed configurations. Important connections between material parameters are obtained using an analytical solution for two separating surfaces, as well as an analysis of the stress-strain curves for homogeneous tension for different degradation and interpolation functions. A complete system of equations is presented in the undeformed and deformed configurations. All the phase field parameters are obtained utilizing the existing first principle simulations for the uniaxial tension of Si crystal in the [100] and [111] directions.
more »
« less
- NSF-PAR ID:
- 10220744
- Date Published:
- Journal Name:
- ArXivorg
- ISSN:
- 2331-8422
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
SUMMARY The occurrence of plate tectonics on Earth is rooted in the physics of lithospheric ductile weakening and shear-localization. The pervasiveness of mylonites at lithospheric shear zones is a key piece of evidence that localization correlates with reduction in mineral grain size. Most lithospheric mylonites are polymineralic and the interaction between mineral phases, such as olivine and pyroxene, especially through Zener pinning, impedes normal grain growth while possibly enhancing grain damage, both of which facilitate grain size reduction and weakening, as evident in lab experiments and field observations. The efficacy of pinning, however, relies on the mineral phases being mixed and dispersed at the grain scale, where well-mixed states lead to greater mylonitization. To model grain mixing between different phases at the continuum scale, we previously developed a theory treating grain-scale processes as diffusion between phases, but driven by imposed compressive stresses acting on the boundary between phases. Here we present a new model for shearing rock that combines our theory for diffusive grain mixing, 2-D non-Newtonian flow and two-phase grain damage. The model geometry is designed specifically for comparison to torsional shear-deformation experiments. Deformation is either forced by constant velocity or constant stress boundary conditions. As the layer is deformed, mixing zones between different mineralogical units undergo enhanced grain size reduction and weakening, especially at high strains. For constant velocity boundary experiments, stress drops towards an initial piezometric plateau by a strain of around 4; this is also typical of monophase experiments for which this initial plateau is the final steady state stress. However, polyphase experiments can undergo a second large stress drop at strains of 10–20, and which is associated with enhanced phase mixing and resultant grain size reduction and weakening. Model calculations for polyphase media with grain mixing and damage capture the experimental behaviour when damage to the interface between phases is moderately slower or less efficient than damage to the grain boundaries. Other factors such as distribution and bulk fraction of the secondary phase, as well as grain-mixing diffusivity also influence the timing of the second stress drop. For constant stress boundary conditions, the strain rate increases during weakening and localization. For a monophase medium, there is theoretically one increase in strain rate to a piezometric steady state. But for the polyphase model, the strain rate undergoes a second abrupt increase, the timing for which is again controlled by interface damage and grain mixing. The evolution of heterogeneity through mixing and deformation, and that of grain size distributions also compare well to experimental observations. In total, the comparison of theory to deformation experiments provides a framework for guiding future experiments, scaling microstructural physics to geodynamic applications and demonstrates the importance of grain mixing and damage for the formation of plate tectonic boundaries.more » « less
-
In this work a transformation strain gradient enhancement is introduced into a phenomenological constitutive model for the pseudoelastic behavior of shape memory alloys. The constitutive model is able to capture several unique features of the constitutive response of these materials during the transformation between austenite and martensite during the pseudoelastic response. These features include the asymmetry in the initial transformation stresses in tension versus compression, the asymmetry in the transformation strains in tension and compression, and finally the asymmetry in the hardening behavior in tension and compression. In fact, experiments have shown that untrained NiTi exhibits hardening during its transformation in compression, but softening for tensile loading. It is this softening behavior that motivates the need for the introduction of the transformation strain gradient into the constitutive modeling. Transformation strain gradient effects are introduced via a phase variable that describes the extent of transformation. The free energy of the material then depends on gradients of the phase variable, which introduces a material length scale into the theory. The governing equation for the phase variable is developed from a microforce balance and continuum thermodynamics analysis. The model is implemented in the commercial finite element software Abaqus through user defined subroutines and several numerical simulations are performed to illustrate the model response and lack of numerical mesh-dependency of the results.more » « less
-
A thermodynamically consistent multiphase phase-field approach for stress and temperature-induced martensitic phase transformation at the nanoscale and under large strains is developed. A total of N independent order parameters are considered for materials with N variants, where one of the order parameters describes A ↔ M transformations and the remaining N − 1 independent order parameters describe the transformations between the variants. A non-contradictory gradient energy is used within the free energy of the system to account for the energies of the interfaces. In addition, a non-contradictory kinetic relationships for the rate of the order parameters versus thermodynamic driving forces is suggested. As a result, a system of consistent coupled Ginzburg-Landau equations for the order parameters are derived. The crystallographic solution for twins within twins is presented for the cubic to tetragonal transformations. A 3D complex twins within twins microstructure is simulated using the developed phase-field approach and a large-strain-based nonlinear finite element method. A comparative study between the crystallographic solution and the simulation result is presented.more » « less
-
The nanoscale multiphase phase-field model for stress and temperature-induced multivariant martensitic transformation under large strains developed by the authors in Basak and Levitas (J Mech Phys Solids 113:162–196, 2018) is revisited, the issues related to the gradient energy and coupled kinetic equations for the order parameters are resolved, and a thermodynamically consistent non-contradictory model for the same purpose is developed in this paper. The model considers N+1 order parameters to describe austenite and N martensitic variants. One of the order parameters describes austenite↔martensite transformations, and the remaining N order parameters, whose summation is constrained to the unity, describe the transformations between the variants. A non-contradictory gradient energy is used within the free energy of the system to account for the energies of the interfaces. In addition, a kinetic relationship for the rate of the order parameters versus thermodynamic driving forces is suggested, which leads to a system of consistent coupled Ginzburg–Landau equations for the order parameters. An approximate general crystallographic solution for twins within twins is presented, and the explicit solution for the cubic to tetragonal transformations is derived. A large strain-based finite element method is developed for solving the coupled Ginzburg–Landau and elasticity equations, and it is used to simulate a 3D complex twins within twins microstructure. A comparative study between the crystallographic solution and the simulation results is presented.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- The DOI is not currently available.
