In this paper, we present concurrent atomistic-continuum (CAC) simulations of the hydrogen (H) diffusion along a grain boundary (GB), nearby which a large population of dislocations are piled up, in a plastically deformed bi-crystalline bcc iron sample. With the microscale dislocation slip and the atomic structure evolution at the GB being simultaneously retained, our main findings are: (i) the accumulation of tens of dislocations near the H-charged GB can induce a local internal stress as high as 3 GPa; (ii) the more dislocations piled up at the GB, the slower the H diffusion ahead of the slip–GB intersection; and (iii) H atoms diffuse fast behind the pileup tip, get trapped within the GB, and diffuse slowly ahead of the pileup tip. The CAC simulation-predicted local H diffusivity, Dpileup−tip, and local stresses, σ, are correlated with each other. We then consolidate such correlations into a mechanics model by considering the dislocation pileup as an Eshelby inclusion. These findings will provide researchers with opportunities to: (a) characterize the interplay between plasticity, H diffusion, and crack initiation underlying H-induced cracking (HIC); (b) develop mechanism-based constitutive rules to be used in diffusion–plasticity coupling models for understanding the interplay between mechanical and mass transport in materials at the continuum level; and (c) connect the atomistic deformation physics of polycrystalline materials with their performance in aqueous environments, which is currently difficult to achieve in experiments.
- NSF-PAR ID:
- 10349055
- Date Published:
- Journal Name:
- Modelling and Simulation in Materials Science and Engineering
- ISSN:
- 0965-0393
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
Abstract This work studies the use of laser shock peening (LSP) to improve back stress in additively manufactured (AM) 316L parts. Unusual hardening behavior in AM metal due to tortuous microstructure and strong texture poses additional design challenges. Anisotropic mechanical behavior complicates application for mechanical design because 3D printed parts will behave differently than traditionally manufactured parts under the same loading conditions. The prevalence of back-stress hardening or the Bauschinger effect causes reduced fatigue life under random loading and dissipates beneficial compressive residual stresses that prevent crack propagation. LSP is known to improve fatigue life by inducing compressive residual stress and has been applied with promising results to AM metal parts. It is here demonstrated that LSP may also be used as a tool for mitigating tensile back-stress hardening in AM parts, thereby reducing anisotropic hardening behavior and improving design use. It is also shown that the method of application of LSP to additively manufactured parts is key for achieving effective back-stress reduction. Back stress is extracted from additively manufactured dog bone samples built in both XY and XZ directions using hysteresis tensile. Both LSPed and as-built conditions are tested and compared, showing that LSPed samples exhibit a significant reduction to back stress when the laser processing is applied to the sample along the build direction. Electron backscatter diffraction (EBSD) performed under these conditions elucidates how grain morphologies and texture contribute to the observed improvement. Crystal plasticity finite element (CPFE) modeling develops insights as to the mechanisms by which this reduction is achieved in comparison with EBSD results. In particular, the difference in plastic behavior across build orientations of identified crystal planes and grain families are shown to impact the degree of LSP-induced back-stress reduction that is sustained through tensile loading.more » « less
-
more » « less
Abstract For some polycrystalline materials such as austenitic stainless steel, magnesium, TATB, and HMX, twinning is a crucial deformation mechanism when the dislocation slip alone is not enough to accommodate the applied strain. To predict this coupling effect between crystal plasticity and deformation twinning, we introduce a mathematical model and the corresponding monolithic and operator splitting solvers that couple the crystal plasticity material model with a phase field twining model such that the twinning nucleation and propagation can be captured via an implicit function. While a phase field order parameter is introduced to quantify the twinning induced shear strain and corresponding crystal reorientation, the evolution of the order parameter is driven by the resolved shear stress on the twinning system. To avoid introducing an additional set of slip systems for dislocation slip within the twinning region, we introduce a Lie algebra averaging technique to determine the Schmid tensor throughout the twinning transformation. Three different numerical schemes are proposed to solve the coupled problem, including a monolithic scheme, an alternating minimization scheme, and an operator splitting scheme. Three numerical examples are utilized to demonstrate the capability of the proposed model, as well as the accuracy and computational cost of the solvers.
-
In this paper, we perform concurrent atomistic-continuum (CAC) simulations to (i) characterize the internal stress induced by the microscale dislocation pileup at an atomically structured interface; (ii) decompose this stress into two parts, one of which is from the dislocations behind the pileup tip according to the Eshelby model and the other is from the dislocations at the pileup tip according to a super-dislocation model; and (iii) assess how such internal stresses contribute to the atomic-scale phase transformations (PTs), reverse PTs, and twinning. The main novelty of this work is to unify the atomistic description of the interface and the coarse-grained (CG) description of the lagging dislocations away from the interface within one single framework. Our major findings are: (a) the interface dynamically responds to a pileup by forming steps/ledges, the height of which is proportional to the number of dislocations arriving at the interface; (b) when the pre-sheared sample is compressed, a direct square-to-hexagonal PT occurs ahead of the pileup tip and eventually grows into a wedge shape; (c) upon a further increase of the loading, part of the newly formed hexagonal phase transforms back to the square phase. The square product phase resulting from this reverse PT forms a twin with respect to the initial square phase. All phase boundaries (PBs) and twin boundaries (TBs) are stationary and correspond to zero thermodynamic Eshelby driving forces; and (d) the stress intensity induced by a pileup consisting of 16 dislocations reduces the stress required for initiating a PT by a factor of 5.5, comparing with that in the sample containing no dislocations. This work is the first characterization of the behavior of PTs/twinning resulting from the reaction between a microscale dislocation slip and an atomically structured interface.more » « less
-
more » « less
Abstract Plastic deformation of olivine at relatively low temperatures (i.e., low‐temperature plasticity) likely controls the strength of the lithospheric mantle in a variety of geodynamic contexts. Unfortunately, laboratory estimates of the strength of olivine deforming by low‐temperature plasticity vary considerably from study to study, limiting confidence in extrapolation to geological conditions. Here we present the results of deformation experiments on olivine single crystals and aggregates conducted in a deformation‐DIA at confining pressures of 5 to 9 GPa and temperatures of 298 to 1473 K. These results demonstrate that, under conditions in which low‐temperature plasticity is the dominant deformation mechanism, fine‐grained samples are stronger at yield than coarse‐grained samples, and the yield stress decreases with increasing temperature. All samples exhibited significant strain hardening until an approximately constant flow stress was reached. The magnitude of the increase in stress from the yield stress to the flow stress was independent of grain size and temperature. Cyclical loading experiments revealed a Bauschinger effect, wherein the initial yield strength is higher than the yield strength during subsequent cycles. Both strain hardening and the Bauschinger effect are interpreted to result from the development of back stresses associated with long‐range dislocation interactions. We calibrated a constitutive model based on these observations, and extrapolation of the model to geological conditions predicts that the strength of the lithosphere at yield is low compared to previous experimental predictions but increases significantly with increasing strain. Our results resolve apparent discrepancies in recent observational estimates of the strength of the oceanic lithosphere.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1088/1361-651X/ac8698
