Stainless steel is a commonly used material in safety-important components of nuclear power plants. In order to study degradation mechanisms in stainless steels, like crack initiation and propagation, it is important to characterize the degree of plastic strain on microstructural level. One way to estimate local plastic strain is by measuring local crystal orientations of the scanned surfaces: the electron backscatter diffraction (EBSD) measurements on stainless steel revealed a strong correlation between the spread of crystal orientations within the individual grains and the imposed macroscopic plastic strain. Similar behavior was also reproduced by finite element simulations where stainless steel was modeled by an anisotropic elasto-plastic constitutive model. In that model the anisotropic Hill’s plasticity function for yield criteria was used and calibrated against the EBSD measurements and macroscopic tensile curve.
In this work the Hill’s phenomenological model is upgraded to a more sophisticated crystal plasticity model where plastic deformation is assumed to be a sum of crystalline slips in all activated slip systems. The hardening laws of Peirce, Asaro and Needleman and of Bassani and Wu are applied in crystal plasticity theory and implemented numerically within the user subroutine in ABAQUS. The corresponding material parameters are taken from literature for 316L stainless steel. Finite element simulations are conducted on the analytical Voronoi tessellation with 100 grains and initial random crystallographic orientations. From the simulations, crystal and modified crystal deformation parameters are calculated, which quantify mean and median spread of crystal orientations within individual grains with respect to central grain orientation. The results are compared to EBSD measurements and previous simulations performed with Hill’s plasticity model.
Studying the effect of grain boundaries in dislocation density based crystal-plasticity finite element simulations