SummaryThe brittle-ductile behavior of cracks has long been an area of intensive study. In particular, body centered cubic (bcc) iron has been studied often in the recent past in the framework of continuum (e.g. ) and atomistic models (e.g. , and as well in experiments (e.g. Marsh et al. [9], Šmida and Bošanský [10]) due to the structural steels applications. The models consider usually plane strain conditions (2D) along the crack front and temperature of 0 K. However, free sample surface (where plane stress conditions are expected in the normal direction) may influence the ductile-brittle behavior, as well as the thermal atomic motion. Such studies require 3D atomistic simulations.This contribution is devoted to crack simulations by molecular dynamic (MD) technique in 3D bcc iron crystals. We use our new MD code for parallel processing in MPI. Interatomic interactions in bcc iron are described using N-body potentials of Finnis-Sinclair type (Ackland et al. [11,12]). The 3D codes have been tested in perfect samples under simple uni-axial tension and in thermal simulations. The simulated thermal expansion in bcc iron agrees well with experimental data (see Machová [13]), as well as the phonon frequency spectra [11,12].In crack simulations we consider a central pre-existing Griffith (through) crack loaded in tension mode I. The relatively long crack is embedded in a thin bcc iron crystal having the basic cubic {100} orientation. The crack is introduced by removing part of atoms from the central plane, i.e. its initial blunting corresponds to the lattice parameter. Crack surface lies on a (001) plane, crack front is oriented along the [010] direction, and the direction of the potential crack extension is [100]. The interatomic interactions across the free crack faces are not allowed. Surface relaxation has been performed before a loading to avoid its influence on crack tip processes.We performed the 3D crack simulations at temperatures of 0 K and 300 K. At these temperatures, the samples were loaded symmetrically in the <001> directions by prescribing external forces F ext distributed homogeneously at individual atoms lying in several surface layers, similar to [11]. The samples were loaded slowly, gradually in time. While at the temperature of 0 K brittle crack initiation has been observed, at 300 K dislocation emission and slip processes on {110} and {112} planes have been detected. The slip processes start at the free sample surfaces, which is in agreement with our stress analysis and MD simulations by Zhou et al. [14]. The process begins on {110} planes via emission of a curved dislocation (~ quarter circle loop) from the corners, where the crack penetrates the free sample surface. We further show that the slip processes on the inclined {110} planes cause crack tip blunting and hinders crack growth. The slip processes on the oblique {112} planes make jogs in the crack front and enable a slow plastic crack growth.
ECF16Our results are in agreement with continuum predictions that the microscopic processes at the ...