metals Article Increasing Low-Temperature Toughness of 09Mn2Si Steel through Lamellar Structuring by Helical Rolling Sergey Panin 1,2,* , Ilya Vlasov 1 , Dmitry Moiseenko 1,3, Pavel Maksimov 1, Pavlo Maruschak 4 , Alexander Yakovlev 2,5, Julia Gomorova 1, Ivan Mishin 1 and Siegfried Schmauder 3 ���������� ������� Citation: Panin, S.; Vlasov, I.; Moiseenko, D.; Maksimov, P.; Maruschak, P.; Yakovlev, A.; Gomorova, J.; Mishin, I.; Schmauder, S. Increasing Low-Temperature Toughness of 09Mn2Si Steel through Lamellar Structuring by Helical Rolling. Metals 2021, 11, 352. https://doi.org/10.3390/met11020352 Academic Editor: Marcello Cabibbo Received: 27 January 2021 Accepted: 16 February 2021 Published: 19 February 2021 Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affil- iations. Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/). 1 Institute of Strength Physics and Materials Sciences SB RAS, 634055 Tomsk, Russia; vlasov.ilya.viktorovich@gmail.com (I.V.); moisey@rocketmail.com (D.M.); mpv@ispms.ru (P.M.); julia.gomorova@gmail.com (J.G.); mishinv1@yandex.ru (I.M.) 2 School of Advanced Manufacturing Technologies, Tomsk Polytechnic University, 634055 Tomsk, Russia; alexandryakovl@gmail.com 3 Institute for Materials Testing, Materials Science and Strength of Materials (IMWF), University of Stuttgart, 70569 Stuttgart, Germany; siegfried.schmauder@imwf.uni-stuttgart.de 4 Department of Industrial Automation, Ternopil Ivan Puluj National Technical University, 46001 Ternopil, Ukraine; maruschak.tu.edu@gmail.com 5 TomskNIPIneft JSC, 634027 Tomsk, Russia * Correspondence: svp@ispms.ru; Tel.: +7-3822-286-904 Abstract: The aim of the paper was to investigate the helical rolling parameters (a number of passes) for the microstructural modification and the low-temperature impact toughness improvement of the 09Mn2Si High Strength Low-Alloyed (HSLA) steel. In order to achieve this purpose, work spent to crack initiation and propagation was analyzed and compared with patterns of fracture surfaces. The microstructure and impact toughness values were presented in the temperature range from +20 to –70◦C. Also, the fracture mechanisms in individual regions on the fracture surfaces were discussed. In addition, a methodology for computer simulation of the process was developed and implemented within the framework of the excitable cellular automata method and its integration with the kinetic theory of fracture. Finally, a theoretical analysis of the effect of grain shapes and orientations on the strain response patterns of a certain meso-volume simulating the material after the helical rolling was carried out. Keywords: high strength low alloyed steels; helical rolling; low-temperature impact toughness; lamellar microstructure; computer simulation; excitable cellular automata method 1. Introduction The active development of the Arctic region justifies enhanced requirements for the mechanical properties of materials used for application in the cold climate. To solve this issue, in particular upon the construction of pipelines, steels of new generations (the X80 and X100 grades, or higher) are being designed, which are characterized by a continuously increasing ratio between ductility and strength [1,2]. Achieving high fracture resistance is possible due to the microstructure refinement by manufacturing processes that causes smaller bainitic grain sizes, reduced inclusion contents, etc. [3]. However, high strength low alloyed steels (HSLA) of the X52 grade, which have been used in the pipeline construction for a long time, continue to be of interest. Their properties, especially impact toughness at low climatic temperatures, have to be improved. It is known that the X65 and X70 grade steels exhibit a high ductility-to-strength ratio, whereas the X100 ones have its lower level and their fracture is, typically, rather brittle. The above metallurgical (technological) aspects of increasing impact toughness at low climatic temperatures have also been discussed in various papers [4–6]. A special role in improving this property is assigned to the bainite formation [7]. Moreover, upper bainite, formed due Metals 2021, 11, 352. https://doi.org/10.3390/met11020352 https://www.mdpi.com/journal/metals https://www.mdpi.com/journal/metals https://www.mdpi.com https://orcid.org/0000-0001-7623-7360 https://orcid.org/0000-0001-9110-8313 https://orcid.org/0000-0002-0880-2898 https://orcid.org/0000-0001-6022-3853 https://orcid.org/0000-0002-3001-0512 https://doi.org/10.3390/met11020352 https://doi.org/10.3390/met11020352 https://creativecommons.org/ https://creativecommons.org/licenses/by/4.0/ https://creativecommons.org/licenses/by/4.0/ https://doi.org/10.3390/met11020352 https://www.mdpi.com/journal/metals https://www.mdpi.com/2075-4701/11/2/352?type=check_update&version=2 Metals 2021, 11, 352 2 of 28 to the presence of martensite–austenite (M–A) and an increased amount of low-angle lath boundaries, significantly decreases the ductility in these conditions. To quantify crack resistance of a material, the ductile-brittle transition (DBT) tempera- ture (DBTT) under impact bending test is applied [8]. This parameter actually characterizes the permissible operating temperature range. In [9], the microstructural aspects of em- brittlement and toughening have been discussed in terms of the DBT. It has been shown that micro-voids hardly nucleate at fine M–A regions or carbides (they factually grow, but through bainite laths with high misorientation between them). Also, it is known from [10] that the DBTT is not intrinsic to a material, but to a significant extent depends on the specimen type and test conditions, which actually determine the plastic constraint. The DBTT aspect is typically discussed in relation to materials used in nuclear power plants [11], since its development is due to aging through thermal and/or irradiation effects during long-term operation. Also, it is actively being studied for high-strength stainless pipeline steels. The key parameters that affect the DBT value are grain sizes, as well as the crack propagation direction relative to the boundaries of microstructural elements (determined by the rolling direction) [12]. If strength of specimens cut at different angles to the rolling axis has not changed significantly, the DBTT is extremely sensitive to the predominant grain orientation. Thus, the targeted formation of the boundaries of microstructural elements (especially, high-angle ones) can be an effective way to control the DBTT. Finally, many authors have discussed the DBTT aspect for relatively inexpensive stainless steels with low nickel contents [13]. It has been shown that the DBT (as expected) is strongly affected by the crack propagation direction, depending on crystallographic texture and microstructural morphology. Moreover, the DBTT becomes extremely low in the case of fracture accompanying delamination. Nowadays, “strain nano-technologies” [14] are being actively developed that based on the formation of micro-scale structural elements in metals and alloys, according to the ‘top-bottom’ principle. It has been shown that (a) the Hall–Petch law on the relationship between grain sizes and yield stresses ceases to be fulfilled in the dimension range below 10 nm; (b) it is not necessary to reduce the sizes of all microstructural elements below the “generally accepted” threshold of 100 nm to ensure high strength properties and maintain acceptable plasticity [15], but it is important to form a hierarchically organized microstruc- ture in this case [16]. In addition, one of the critical factors is the strong interaction of the microstructural elements with each other. This causes an increase in strength and excludes rapid strain localization processes. Within this formulation of the problem, a promising direction is associated with the formation of “heterogeneous materials”, which is widely discussed in recent papers in the field of materials science [17]. It should be noted that all structural materials are inherently heterogonous in one way or another. According to [18], the principles of the heterogonous microstructure formation are determined both by the available technologies for processing metals and alloys, and by the possibility of its hierarchical alignment. Enhancing strength of these materials is achieved by heterogeneity of their microstructures, crystal structures, or compositions. The domain sizes can be in the range from micro- to millimeters, and their shapes vary to form very diverse material systems. Additional factors to ensure extraordinary strength and ductility are “back-stress strengthening and back-stress work hardening”. From the standpoint of the mechanisms providing the relaxation-ability of such materials, “statistically stored dislocations and geometrically necessary dislocations” are typically discussed. At a higher (meso-scale) level, grain shapes and sizes, as well as texture and grain size variations, are crucial. The “lamellar microstructure with elongated soft coarse-grained domains embedded in the ultrafine grain matrix” is considered to be the “optimal” grain one [19]. Thus, the principle of the hierarchical microstructure formation in the hetero- geneous materials is determined by the back-stress produced due to dislocation pile-up at the domain boundaries, when their interface density is high. On the other hand, the Metals 2021, 11, 352 3 of 28 spacing between them should be sufficient to ensure the possibility of dislocation strains. The authors have called this principle “strain partitioning among heterogeneous domains”, when the strain gradient increase the back-stress work hardening. Undoubtedly, one of the factors limiting the widespread industrial use of the hetero- geneous materials is the lack of high-performance procedures for the formation of such a microstructure. The above-described principle of the lamellar microstructure formation at various scale levels is clearly illustrated in [20]. It is shown that following the natural principles of the hierarchical microstructure alignment (like a seashell) enables to achieve high strength properties (with satisfactory ductility) for structural steels. This implies the formation of hierarchical sandwich microstructures, including gradient grain sizes, as well as both martensite and austenite phases. In this case, the lamellar microstructure has been observed both at the highest structural level (over the cross-sectional landscape of the material) and at the nano-scale (FCC austenite lamellae 5 nm wide alternate with 6 nm HCP martensite ones). It should be emphasized once again that the formation of such small microstructural elements does not suppress development of dislocation strain mechanisms even under conditions when grain sizes are about 20 nm. Aspects related to the possibility of strength increasing with toughness retention (strength vs. toughness or hardness vs. ductility) for materials of various classes are discussed in a large number of papers (for example [21,22]) using the concepts of fracture mechanics. In [23], transcrystalline fracture patterns are shown from the perspective of various scale levels as applied to structural steels. At the micro-level, it is proposed to use local fracture stress values to assess the local “work-of-fracture” parameter. Furthermore, at the nano-scale level, ahead of a sharp crack is proposed “to increase the compliance by a cooperative movement of atoms (involving extra work) to allow the crack-tip bond to displace sufficiently”. These representations are consistent with the plastic distortion mechanism considered in V.E. Panin’s paper [24]. The formation of such a hierarchical microstructure, including the lamellar one, can be achieved using the radial–shear rolling procedure. The results of studies of the strain mechanisms for various materials subjected to this processing have been widely published recently [25–29]. In all these papers, the importance of the lamellar microstructure for- mation in the conditions of multilevel structuring is noted for increasing strength and maintaining toughness of rolled bars. A fairly complete review devoted to the formation of multiscale hierarchical mi- crostructures, including hardening and providing ductility and plasticity by forming a laminate structure, is presented in [30]. According to the authors, elastic modulus and yield strength obey the “rule of mean”. However, elongation and fracture toughness cannot be estimated by this principle. It is known that (for laminate composites) elongation is determined by a number of parameters: layer thickness, interface structure, gradient struc- ture, strain hardening exponent, strain rate, etc. They, in turn, depend on crack deflection, blunting, and bridging, as well as stress redistribution, local plane stress deformation, and interfacial delamination cracking. Thus, the ideology of the formation of the laminate composites can be efficiently used to develop the heterogeneous materials with increased strength and toughness. In this case, the laminate microstructure can be formed at various structural levels. In general, the implementation of the laminate microstructure as an effective way to improve fracture toughness can be traced to the example of the Egyptian pyramids [31]. The laminate microstructure has shown particular efficiency in terms of increasing strength and toughness in the manufacture of bladed weapons such as Damascus steel, Japanese combat swords, etc. [32]. In this case, the composites are formed from single base metals, while the layer-by-layer joining of stronger (for example, ultra-high carbon (UHC)) and more ductile (for example, Fe-3Si) steels is considered to be more efficient [33]. A fairly detailed review of research results related to laminate metal composites (LMC) is presented in [34]. It has been shown that one of the effective methods for Metals 2021, 11, 352 4 of 28 the formation of such composites is accumulated roll bonding (ARB), and the highest both strength and toughness values are achieved by forming packets from materials with high and low strength. Numerous examples show that damage critical properties (fracture toughness, fatigue, and impact behavior) are largely depend on local delamination at the layer interfaces. In addition, the key factors that determine the mechanisms of distortion and failure are the parameters of the laminate formation, its ‘architecture’ and the microstructure of the layers and their interfaces. It has been shown in [35] that post-processing heat treatment of a twelve-layer “UHC steel/mild steel” composite significantly improves its impact toughness by enhancing strength at the layer interfaces during tests at a temperature of 79 ◦C. One of the most effective methods of increasing the strength properties of conventional steels (for example, the X52 grade HSLA ones) is to form certain microstructures, phase compositions, and internal stresses in the bulk metal [36]. Varying these parameters at the processing stage enables to obtain the required hardness, ductility, and crack resistance. Most often, processing includes heat treatment to form the optimal phase composition and plastic strain for additional hardening. Low-alloyed low-carbon steels of the pearlite class are best suited to such procedures [37,38]. In this regard, development of complex thermo-mechanical processing technologies is relevant to obtain materials that combine both high strength and sufficient ductility. To control and determine the optimal processing modes, the high-temperature helical rolling (HR) is widely applied, which can strain and heat treat the studied materials to a predetermined degree [39]. The aim of this research has been to investigate the HR parameters (a number of passes) for the microstructural modification and the low-temperature impact toughness improvement of the 09Mn2Si HSLA steel. In order to achieve this purpose, work spent to crack initiation and propagation has been analyzed and compared with patterns of fracture surfaces. The second section describes the details of the experimental studies, including the HR modes. In the third one, the results of studies of the microstructure (optical and SEM micrographs) are shown and impact toughness values are presented in the temperature range from +20 to –70 ◦C. Also, the failure mechanisms in individual regions on the fracture surfaces are discussed. In the fourth section, the authors propose a methodology for computer simulation of the process within the framework of the excitable cellular automata method and its integration with the kinetic theory of fracture (S.N. Zhurkov). A theoretical analysis of the effect of grain shapes and orientations on the strain response patterns of a certain meso-volume simulating the material after the HR processing is carried out in the fifth section. The results presented in Section 5 are not a complete analogy of the full-scale experiment, but are aimed at identifying the role of the processes at the internal interfaces on the propagation of energy flows. Finally, the authors draw conclusions based on the obtained results. It should be noted that the effect of such processing on the resistance to fatigue fracture for this steel has been considered in the previous paper of the authors [40]. 2. Materials and Methods The HR processing of the 09Mn2Si HSLA steel (chemical composition: (C ≤ 0.12; Si = 0.5–0.8; Mn = 1.3–1.7; Ni ≤ 0.3; S ≤ 0.04; p ≤ 0.035; Cr ≤ 0.3; N ≤ 0.008; Cu ≤ 0.3; As ≤ 0.08; Fe = 96–97) was carried out using an RSP “14–40” three-roll mini-mill for screw rolling in the temperature range of 850–500 ◦C with a decrease in the rolling temperature at each subsequent pass (i.e., thermal schedule of the HR were are as follows: first pass: 850 ◦C; second pass: 750 ◦C; third pass: 650 ◦C; fourth pass: 580 ◦C; fifth pass: 500 ◦C.). After each rolling pass, water quenching was performed. Both three- and five–pass HR processing procedures were completed. The obtained results after various numbers of passes were compared to assess the dynamics of changes in the mechanical properties and determine the optimal mode of thermo-mechanical processing. The total degree of true logarithmic strain/accumulated strain of the steel after every pass of helical rolling was Metals 2021, 11, 352 5 of 28 equal to: (first) 0.09/8.8%; (second) 0.23/21.9%; (third) 0.45/41.5%; (fourth) 0.63/57.9%; (fifth) 0.80/73.4%. Due to the metal inhomogeneity along the bar cross-section after the HR processing, the grain microstructure was evaluated within four regions (Figure 1a): the surface layer; area 1 (near the bar surface); area 2 (transition), and area 3 (the bar core). This differentiation is conditional, since the change in the microhardness values was smooth over the cross- section (Figure 1b). A ‘PMT-3′ microhardness tester (LOMO, St. Peterburg, Russia) was used to assess microhardness levels. Metals 2021, 11, 352 5 of 29 Due to the metal inhomogeneity along the bar cross-section after the HR processing, the grain microstructure was evaluated within four regions (Figure 1a): the surface layer; area 1 (near the bar surface); area 2 (transition), and area 3 (the bar core). This differentia- tion is conditional, since the change in the microhardness values was smooth over the cross-section (Figure 1b). A ‘PMT-3′ microhardness tester (LOMO, St. Peterburg, Russia) was used to assess microhardness levels. (a) (b) II I IV III (c) (d) Figure 1. Scheme of the bar volume partitioning into areas (a), microhardness curve over the bar cross-section (b), the scheme of the specimen preparation for the Charpy impact tests (c), and the designation of the fracture areas in them (d): I—crack initiation area; II—crack propagation area; III—shear lips; IV—subcritical crack propagation area (rupture). Specimens for the Charpy impact tests were cut from the bars by electrical discharge machining (Figure 1c). Impact toughness was determined using an ‘Instron 450MPX’ au- tomated pendulum and the Charpy specimens of 55 × 10 × 10 mm3 (actually contained the metal from areas 2 and 3) with a V–notch 2 mm deep (according to ASTM E 23:2007). The tests were carried out at temperatures (T) of +20, 0, –20, –40, and –70 °C. During 10 min immediately before testing, the specimens were cooled in a silicone oil environment using a ‘Lauda rp 870′ installation (chiller). Duration between the specimen removal from the cooling chamber and the test beginning was within 5 s. The lowest test temperature threshold was determined by the environmental operating temperatures in the Arctic re- gion for the tested HSLA steel, as well as the possible temperature range of the ‘Lauda rp 870′ installation. The fracture surfaces of the specimens were observed using a ‘LEO EVO 50′ scanning electron microscope (Zeiss, Oberkochen, Germany) at the ‘NANOTECH’ Center for Collective Use of the ISPMS SB RAS. Figure 1. Scheme of the bar volume partitioning into areas (a), microhardness curve over the bar cross-section (b), the scheme of the specimen preparation for the Charpy impact tests (c), and the designation of the fracture areas in them (d): I—crack initiation area; II—crack propagation area; III—shear lips; IV—subcritical crack propagation area (rupture). Specimens for the Charpy impact tests were cut from the bars by electrical discharge machining (Figure 1c). Impact toughness was determined using an ‘Instron 450MPX’ automated pendulum and the Charpy specimens of 55 × 10 × 10 mm3 (actually contained the metal from areas 2 and 3) with a V–notch 2 mm deep (according to ASTM E 23:2007). The tests were carried out at temperatures (T) of +20, 0, −20, −40, and −70 ◦C. During 10 min immediately before testing, the specimens were cooled in a silicone oil environment using a ‘Lauda rp 870′ installation (chiller). Duration between the specimen removal from the cooling chamber and the test beginning was within 5 s. The lowest test temperature threshold was determined by the environmental operating temperatures in the Arctic region for the tested HSLA steel, as well as the possible temperature range of the ‘Lauda rp 870′ installation. The fracture surfaces of the specimens were observed using a ‘LEO EVO 50′ scanning electron microscope (Zeiss, Oberkochen, Germany) at the ‘NANOTECH’ Center for Collective Use of the ISPMS SB RAS. 3. Results 3.1. Effect of the HR Processing on the Microstructure and Microhardness Figure 1b shows the microhardness (Hµ) values on the bar cross-section after the three-pass HR. The test was performed to a depth of 9 mm, while the bar diameter was 18 mm after the HR processing. Compared to the as-received steel, strengthening had occurred throughout the entire volume of the bar with a gradual Hµ decrease in the core. Metals 2021, 11, 352 6 of 28 The greatest increase in microhardness was observed in the surface layer up to 3 mm thick. Table 1 presents the microhardness values in area 2 according to Figure 1a for the as-received steel specimen and after both three- and five-pass HR processing. The lowest Hµ values were typical for the as-received steel, while the increase in the number of the HR passes resulted in a gradual Hµ enhancing up to 20%. Table 1. Microhardness values in area 2 according to Figure 1a. Steel Conditions Microhardness HV As-received 97 ± 7 After the three-pass HR 116 ± 3 After the five-pass HR 123 ± 7 The as-received steel was characterized by the ferrite-pearlite microstructure with a characteristic equiaxial grain size of 21 ± 2 µm (Figure 2a). Some striping was observed due to the use of rolling procedures for its manufacturing. In the surface layer after the five-pass HR, a fine microstructure was presumably composed of ferrite and fractured pearlite plates (Figure 2b). At a depth of more than 1 mm in the longitudinal cross-section, individual elongated ferrite grains were distinguished (Figure 2c; area 1), preserved after strains and slight deviated from the main axis towards the bar core at an angle of ~20◦. This orientation was due to a significant reduction in the bar diameter as a result of the HR processing, as well as the inclination of the rolls relative to the workpiece axis. An increase in both the size and the number of elongated ferrite grains was observed towards the bar core. At a distance of more than 4 mm from the surface, the microstruc- ture consisted entirely of elongated grains that formed the lamellar steel microstructure (Figure 2d, area 2). It was similar in the bar core (Figure 2e, area 3). Metals 2021, 11, 352 6 of 28 3. Results 3.1. Effect of the HR Processing on the Microstructure and Microhardness Figure 1b shows the microhardness (Hμ) values on the bar cross-section after the three- pass HR. The test was performed to a depth of 9 mm, while the bar diameter was 18 mm after the HR processing. Compared to the as-received steel, strengthening had occurred throughout the entire volume of the bar with a gradual Hμ decrease in the core. The great- est increase in microhardness was observed in the surface layer up to 3 mm thick. Table 1 presents the microhardness values in area 2 according to Figure 1a for the as-received steel specimen and after both three- and five-pass HR processing. The lowest Hμ values were typical for the as-received steel, while the increase in the number of the HR passes resulted in a gradual Hμ enhancing up to 20%. Table 1. Microhardness values in area 2 according to Figure 1a. Steel Conditions Microhardness HV As-received 97 ± 7 After the three-pass HR 116 ± 3 After the five-pass HR 123 ± 7 The as-received steel was characterized by the ferrite-pearlite microstructure with a characteristic equiaxial grain size of 21 ± 2 μm (Figure 2a). Some striping was observed due to the use of rolling procedures for its manufacturing. In the surface layer after the five- pass HR, a fine microstructure was presumably composed of ferrite and fractured pearlite plates (Figure 2b). At a depth of more than 1 mm in the longitudinal cross-section, individ- ual elongated ferrite grains were distinguished (Figure 2c; area 1), preserved after strains and slight deviated from the main axis towards the bar core at an angle of ~20°. This ori- entation was due to a significant reduction in the bar diameter as a result of the HR pro- cessing, as well as the inclination of the rolls relative to the workpiece axis. An increase in both the size and the number of elongated ferrite grains was observed towards the bar core. At a distance of more than 4 mm from the surface, the microstructure consisted entirely of elongated grains that formed the lamellar steel microstructure (Figure 2d, area 2). It was similar in the bar core (Figure 2e, area 3). As-received After the five-pass HR Surface layer Area 1 Area 2 Area 3 (a) (b) (c) (d) (e) Figure 2. Optical images of the 09Mn2Si steel microstructure in the longitudinal cross-section (explanations for a–e are given in the text above). Figure 2. Optical images of the 09Mn2Si steel microstructure in the longitudinal cross-section (explanations for a–e are given in the text above). The most intense elongation of grains had occurred along the bar axis in the longitudi- nal cross-section, while only a decrease in their sizes, comparing with the as-received steel, was observed in the transverse one (Table 2). In areas 2 and 3, the lamellar microstructure consisted of thin elongated ferritic grains oriented predominantly along the bar axis. In the surface layer at a depth of ~1 mm, a finely dispersed microstructure without a pronounced texture had been formed, where the steel had been subjected to the maximum strains during the HR processing. Any clear boundaries between these layers were not observed. Since more pronounced microstructural changes were found after the five-pass HR, quantitative analysis of the grain microstructure was carried out only for these specimens. Metals 2021, 11, 352 7 of 28 Table 2. Grain sizes in the 09Mn2Si steel after the HR processing. Distance from the Bar Surface, mm Grain Sizes, µm As-Received After the Five-Pass HR Transverse Cross-Section Longitudinal Cross-Section Surface layer 21 ± 2 1.8 ± 0.3 2.7 ± 0.3 1 2 ± 0.3 8.5 ± 1.3 2 2.5 ± 0.3 17.3 ± 2.4 3 3.3 ± 0.4 21.6 ± 3.8 4 6.3 ± 0.6 32.7 ± 5.7 5 6.7 ± 0.8 34.5 ± 4.1 3.2. Impact Toughness The specimens for the Charpy impact test were cut from the bar core (actually from areas 2 and 3, according to Figure 1a,c). Due to the smooth nature of the decrease in the microhardness levels over the cross-section, the specimens possessed the hardened surface layers and the ‘softer’ core. The hardest layers (surface one and area 1), the presence of which during the Charpy impact tests could be accompanied by brittle fracture, were removed at the specimen fabrication stage (Figure 1c). Both three– and five–pass HR processing of the 09Mn2Si steel caused its increased impact toughness over the entire test temperature range (Figure 3a, Table 3). Metals 2021, 11, 352 7 of 28 The most intense elongation of grains had occurred along the bar axis in the longitu- dinal cross-section, while only a decrease in their sizes, comparing with the as-received steel, was observed in the transverse one (Table 2). In areas 2 and 3, the lamellar micro- structure consisted of thin elongated ferritic grains oriented predominantly along the bar axis. In the surface layer at a depth of ~1 mm, a finely dispersed microstructure without a pronounced texture had been formed, where the steel had been subjected to the maximum strains during the HR processing. Any clear boundaries between these layers were not observed. Since more pronounced microstructural changes were found after the five-pass HR, quantitative analysis of the grain microstructure was carried out only for these spec- imens. Table 2. Grain sizes in the 09Mn2Si steel after the HR processing. Distance from the Bar Surface, mm Grain Sizes, μm As-Received After the Five-Pass HR Transverse Cross-Section Longitudinal Cross-Section Surface layer 21 ± 2 1.8 ± 0.3 2.7 ± 0.3 1 2 ± 0.3 8.5 ± 1.3 2 2.5 ± 0.3 17.3 ± 2.4 3 3.3 ± 0.4 21.6 ± 3.8 4 6.3 ± 0.6 32.7 ± 5.7 5 6.7 ± 0.8 34.5 ± 4.1 3.2. Impact Toughness The specimens for the Charpy impact test were cut from the bar core (actually from areas 2 and 3, according to Figure 1a,c). Due to the smooth nature of the decrease in the microhardness levels over the cross-section, the specimens possessed the hardened sur- face layers and the ‘softer’ core. The hardest layers (surface one and area 1), the presence of which during the Charpy impact tests could be accompanied by brittle fracture, were removed at the specimen fabrication stage (Figure c). Both three– and five–pass HR processing of the 09Mn2Si steel caused its increased impact toughness over the entire test temperature range (Figure 3a, Table 3). (a) (b) (c) (d) Figure 3. Impact toughness vs. test temperature curves (a); impact deformation curves for the Charpy specimens in the “load-deflection” coordinates (b–d); 1—as-received; 2—after the three-pass HR; 3—after the five-pass HR. Table 3. Impact toughness (KCV) of the 09Mn2Si steel. Test Temperature, °C KCV, J/cm2 As–Received After the Three–Pass HR After the Five–Pass HR −70 196 ± 16 299 ± 17 317 ± 24 −40 226 ± 21 308 ± 19 332 ± 30 −20 285 ± 20 353 ± 25 333 ± 21 0 308 ± 26 355 ± 28 337 ± 25 20 310 ± 25 387 ± 24 336 ± 26 Figure 3. Impact toughness vs. test temperature curves (a); impact deformation curves for the Charpy specimens in the “load-deflection” coordinates (b–d); 1—as-received; 2—after the three-pass HR; 3—after the five-pass HR. Table 3. Impact toughness (KCV) of the 09Mn2Si steel. Test Temperature, ◦C KCV, J/cm2 As–Received After the Three–Pass HR After the Five–Pass HR −70 196 ± 16 299 ± 17 317 ± 24 −40 226 ± 21 308 ± 19 332 ± 30 −20 285 ± 20 353 ± 25 333 ± 21 0 308 ± 26 355 ± 28 337 ± 25 20 310 ± 25 387 ± 24 336 ± 26 Unlike the as-received steel specimens, there was no critical decrease in toughness (KCV) of the specimens after the HR as the test temperature decreased. At the temperature of +20 ◦C, its higher values were observed for the ones after the three-pass HR processing. On the other hand, the difference in the KCV values decreased with reducing the test temperature for the steel after the HR processing using both modes. Thus, the three-pass HR proved to be a transitional one. In this case, both enhancing the impact toughness values (in comparison with the as-received steel specimens) and its reducing with the decrease in the test temperature were observed. The most important obtained result was the stable retention of the impact toughness level throughout the studied temperature range for the steel specimens after the five-pass Metals 2021, 11, 352 8 of 28 HR processing. The absence of this effect for the ones after the three-pass HR could be justified by the following reasons: (i) the homogeneous ‘lamellar’ microstructure had not been formed due to the fewer number of passes; (ii) the softer metal was due to a lower level of accumulated strains; (iii) quenching in water after the third pass had been carried out from the higher temperature, which had caused the formation of the more non-equilibrium microstructure with lower toughness. 3.3. Analysis of the Loading Diagrams for the Charpy Impact Test Specimens According to the data presented in the ‘load–deflection’ graphs (Figure 3b–d), the curves possessed noticeable variations in their shapes even at the temperature of +20 ◦C. In this case, all specimens fractured in the ductile manner, which confirmed fairly smooth changes in the diagrams without abrupt transitions and bifurcation points (Figure 3b; curve 1). This testified to the effective relaxation of stresses at both macro-crack initiation and propagation stages. The shape of the diagram for the steel after the three-pass HR at the pre-fracture stage should be highlighted (Figure 3b, curve 2, constant p value in the deflection range of 30–40 mm). According to the authors, such curve indicated a significant margin of the material ductility even at the final fracture stage. The curvature of the diagram could be caused by a sharp turn (branching) of the fracture front. On the other hand, the steel specimens after the five-pass HR showed diagrams similar to those for the as-received ones, but with a higher energy intensity to fracture, as evidenced by the greater maximum load (Pmax) values (Figure 4a), as well as a larger area under the curves, reflected work spent to the specimen failure. At the temperatures of −40 and −70 ◦C, impact fracture of the as-received steel speci- mens was accompanied by ductile deformation at both the crack initiation and (partially) propagation stages. However, when the critical crack length was reached, the character of its propagation changed from ductile to brittle. This fact was justified by load breakdown in the diagram (Figure 3c,d; curve 1). On the other hand, the diagrams of the specimens after both three- and five-pass HR showed maintaining of the ductile nature of fracture, which indicated the preservation of relaxation processes (at the micro-scale level) and was accompanied by the noticeable increase in fracture toughness (Table 3). One of the parameters reflected the effect of the HR processing on both strength and dynamic crack resistance of the steel was the maximum load (Pmax) that the specimens withstood upon the Charpy impact tests (Figure 4a). This parameter corresponded to the moment of transition from the macro-discontinuity (crack) initiation to its propagation. Its highest values for the specimens after the five-pass HR, primarily at the low test temperatures, showed the important role of the lamellar (fine-grained) microstructure formation in resistance to impact loading. Metals 2021, 11, 352 9 of 28 (a) (b) (c) Figure 4. Maximum load to the striker during the Charpy impact tests (a); work spent to crack initiation (b) and propaga- tion (c); as-received; 2—after the three-pass HR; 3—after the five-pass HR. 3.4. Calculation of Work Spent to Crack Initiation and Propagation under Impact Loading Since impact toughness was determined by the energy levels for both crack initiation and propagation, their values were analyzed separately. It was found that the crack initi- ation energies Ai for the as-received steel specimens were slightly higher at almost all test temperatures than those for the steel after the HR processing. Its values changed in the range of 60–70 J (Figure 4b). The increase in the number of the HR passes was accompa- nied by a slight decrease in the Ai values, which, nevertheless, did not reduce below 55 J. On the other hand, the crack propagation energies Ap were always higher for the steel specimens the HR processing. After the three-pass HR, Ap decreased from ~270 down to ~180 J as the test temperature reduced, while it was about 200 J after the five-pass HR processing (Figure 4c). The greatest drop in the amount of work spent to crack propaga- tion was observed for the as-received steel specimens, which was primarily associated with their brittle fracture at the subzero test temperatures (Figure 3c,d). In order to estimate the ratio of energy spent to crack initiation and propagation to the total amount of energy spent to the specimen fracture, both Ai/A and Ap/A versus tem- perature graphs were plotted. Analysis of the data shown in Figure 5a enabled to conclude that the Ai/A ratio changed from 25.2% at the temperature of +20 °C up to 40.8% at −70 °C. Such a significant specific part of the energy intensity spent to crack initiation for the as- received steel specimen at the temperature of –70 °C indirectly indicated its embrittle- ment. For the ones after the three-pass HR, this was much less pronounced: there was an increase in the Ai/A ratio from 17.0% at +20 °C up to 24.8% at –70 °C. With increasing the number of passes to five, the Ai/A ratio values practically did not change, remaining at a level of ~21.5% both at +20 and –70 °C. Thus, the increase in the number of the HR passes resulted in the exclusion of the sensitivity of the 09Mn2Si steel to the embrittlement effect of temperature (from the point of view of crack initiation in the investigated temperature range). (a) (b) Figure 5. Energy intensity to failure at the main crack initiation (a) and propagation (b) stages: 1—as-received; 2—after the three-pass HR; 3—after the five-pass HR. Figure 4. Maximum load to the striker during the Charpy impact tests (a); work spent to crack initiation (b) and propagation (c); as-received; 2—after the three-pass HR; 3—after the five-pass HR. Metals 2021, 11, 352 9 of 28 3.4. Calculation of Work Spent to Crack Initiation and Propagation under Impact Loading Since impact toughness was determined by the energy levels for both crack initiation and propagation, their values were analyzed separately. It was found that the crack initiation energies Ai for the as-received steel specimens were slightly higher at almost all test temperatures than those for the steel after the HR processing. Its values changed in the range of 60–70 J (Figure 4b). The increase in the number of the HR passes was accompanied by a slight decrease in the Ai values, which, nevertheless, did not reduce below 55 J. On the other hand, the crack propagation energies Ap were always higher for the steel specimens the HR processing. After the three-pass HR, Ap decreased from ~270 down to ~180 J as the test temperature reduced, while it was about 200 J after the five-pass HR processing (Figure 4c). The greatest drop in the amount of work spent to crack propagation was observed for the as-received steel specimens, which was primarily associated with their brittle fracture at the subzero test temperatures (Figure 3c,d). In order to estimate the ratio of energy spent to crack initiation and propagation to the total amount of energy spent to the specimen fracture, both Ai/A and Ap/A versus temperature graphs were plotted. Analysis of the data shown in Figure 5a enabled to conclude that the Ai/A ratio changed from 25.2% at the temperature of +20 ◦C up to 40.8% at −70 ◦C. Such a significant specific part of the energy intensity spent to crack initiation for the as-received steel specimen at the temperature of −70 ◦C indirectly indicated its embrittlement. For the ones after the three-pass HR, this was much less pronounced: there was an increase in the Ai/A ratio from 17.0% at +20 ◦C up to 24.8% at −70 ◦C. With increasing the number of passes to five, the Ai/A ratio values practically did not change, remaining at a level of ~21.5% both at +20 and −70 ◦C. Thus, the increase in the number of the HR passes resulted in the exclusion of the sensitivity of the 09Mn2Si steel to the embrittlement effect of temperature (from the point of view of crack initiation in the investigated temperature range). Metals 2021, 11, 352 9 of 28 (a) (b) (c) Figure 4. Maximum load to the striker during the Charpy impact tests (a); work spent to crack initiation (b) and propaga- tion (c); as-received; 2—after the three-pass HR; 3—after the five-pass HR. 3.4. Calculation of Work Spent to Crack Initiation and Propagation under Impact Loading Since impact toughness was determined by the energy levels for both crack initiation and propagation, their values were analyzed separately. It was found that the crack initi- ation energies Ai for the as-received steel specimens were slightly higher at almost all test temperatures than those for the steel after the HR processing. Its values changed in the range of 60–70 J (Figure 4b). The increase in the number of the HR passes was accompa- nied by a slight decrease in the Ai values, which, nevertheless, did not reduce below 55 J. On the other hand, the crack propagation energies Ap were always higher for the steel specimens the HR processing. After the three-pass HR, Ap decreased from ~270 down to ~180 J as the test temperature reduced, while it was about 200 J after the five-pass HR processing (Figure 4c). The greatest drop in the amount of work spent to crack propaga- tion was observed for the as-received steel specimens, which was primarily associated with their brittle fracture at the subzero test temperatures (Figure 3c,d). In order to estimate the ratio of energy spent to crack initiation and propagation to the total amount of energy spent to the specimen fracture, both Ai/A and Ap/A versus tem- perature graphs were plotted. Analysis of the data shown in Figure 5a enabled to conclude that the Ai/A ratio changed from 25.2% at the temperature of +20 °C up to 40.8% at −70 °C. Such a significant specific part of the energy intensity spent to crack initiation for the as- received steel specimen at the temperature of –70 °C indirectly indicated its embrittle- ment. For the ones after the three-pass HR, this was much less pronounced: there was an increase in the Ai/A ratio from 17.0% at +20 °C up to 24.8% at –70 °C. With increasing the number of passes to five, the Ai/A ratio values practically did not change, remaining at a level of ~21.5% both at +20 and –70 °C. Thus, the increase in the number of the HR passes resulted in the exclusion of the sensitivity of the 09Mn2Si steel to the embrittlement effect of temperature (from the point of view of crack initiation in the investigated temperature range). (a) (b) Figure 5. Energy intensity to failure at the main crack initiation (a) and propagation (b) stages: 1—as-received; 2—after the three-pass HR; 3—after the five-pass HR. Figure 5. Energy intensity to failure at the main crack initiation (a) and propagation (b) stages: 1—as-received; 2—after the three-pass HR; 3—after the five-pass HR. A matching pattern, only in the ‘mirror’ reflection, was characteristic of the Ap/A pa- rameter (Figure 5b). This enabled to assume the similarity of the fracture micro-mechanisms and morphology of the corresponding fractographic surfaces on the steel specimens after the HR processing, but their significant distinctions for the as-received ones. The HR processing had changed the grain microstructure hierarchy, their shapes and sizes, as well as affected the mechanical properties of the steel, including resistance to impact frac- ture. The decrease in the grain sizes caused the significant increase in the microstructure heterogeneity. This resulted in two main effects: • reducing energy consumption to crack initiation, which was expected due to increasing its structural heterogeneity and enhancing the specific length of the grain boundaries (local heterogeneity zones) in micro-volumes of the steel specimens after the HR [41], in comparison with the as-received ones; Metals 2021, 11, 352 10 of 28 • increasing the energy intensity to crack propagation, which was associated with high energy consumption in such a “micro-composite material” with the lamellar mesostructure [42,43]. At the same time, the lamellar microstructure, formed during the HR processing, affected similar both Ai/A and Ap/A ratios. It should be noted the additional ductility depletion for the 09Mn2Si HSLA steel after the five-pass HR, which was especially no- ticeable from the microhardness test results (Table 1). These patterns indicated that this type of steels could possess increased sensitivity to aging, hydrogen absorption, and other degradation factors during long-term operation [44]. Thus, the main reason for the decrease in the steel toughness at the low test tem- peratures was a gradual decrease in resistance to crack propagation. In this regard, the microstructural modification by the HR processing, despite the increase in microhardness, contributed both to the fracture resistance maintenance almost at the initial level during the crack initiation stage and to its significant enhance upon the macro-discontinuity propaga- tion. From the point of view of increasing the specific area of the grain boundaries along the crack propagation path at the second stage, this result was natural for the ‘lamellar’ microstructure of the steel after the HR processing. However, the fact of impact fracture toughness improving, even at the low test temperatures, required a deeper understanding. 3.5. Shapes of Shear Lips The 09Mn2Si HSLA steel was characterized by a sufficient margin of ductility (tough- ness), which was noticeable by the angle of inclination of the shear lips relative to the specimen axis (Figure 1d), which was about 130◦ (Figure 6). It is known from previous papers [45] that shear lips had formed near the macro-band of strain localization, both in the direction of the maximum shear stresses and perpendicular to the tension axis under the greatest normal ones. Shapes and dimensions of the shear lips characterized the effect of lowering the test temperature on the impact resistance of the material. They were bent at the steel specimens after the three-pass HR, which indicated an increase in the energy intensity to failure. At the same time, the rounded shape of the shear lips, even at the test temperatures of −40 and −70 ◦C, corresponded to a significant margin of fracture toughness. After the five-pass HR, discontinuities were at the shear lips due to the formation of strain localization zones (however, they were also observed on the steel specimens after the three-pass HR). This contributed to an additional increase in the energy intensity to failure. Metals 2021, 11, 352 11 of 29 Steel Conditions Test Temperature, °C –70 –40 –20 0 20 As-received After the three-pass HR After the five-pass HR Figure 6. Effect of the Charpy impact test temperature on the 09Mn2Si steel specimens (both as-received and after the HR processing). 3.6. Fractographic Analysis SEM micrographs of the fracture surfaces are shown in Figures 7 and 8. At the tem- perature of +20 °C, all steel specimens after the HR processing fractured in a rather ductile manner. This fact was evidenced by a noticeable change in their shapes over the cross- sections (Figure 7a–c). The fracture surfaces were characterized by developed shapes with strain localization zones, as well as local delamination. Apparently, depth of these layers largely determined the energy intensity to failure, and as mentioned above, was reflected in the shape of the shear lips [46]. • At the micro-level, rounded pits surrounded by shear fracture spots had been formed on the as-received steel specimen in area 1 (crack initiation, Figure 1g)). In general, it had a ductile-brittle pattern (Figure 7a,g). The fracture surface was quite flat, formed according to the ductile-brittle mechanism. This was confirmed by the presence of multiple “chips”, which alternated with areas of ductile fracture of a dimple type. In area 2 (crack propagation, Figure 1m)), the surface morphology was very rough that indicated a clear heterogeneity in the strain process development at the main crack propagation (failure) stage (Figure 7s). It had been formed by large elongated dim- ples oriented across the crack propagation direction, and also included the large zone 4 (Figure 7s), corresponded to the specimen failure according to Figure 1g, which had a smoother surface with a gradual bend towards the point of impact loading. • After the HR, the crack initiation mechanism was predominantly ductile-brittle, with a minor brittle component (Figure 7b,c,h,i). Dimples had deformed profiles oriented along the crack front direction that indicated an increase in the material ductility and its sufficiently high crack resistance. This was confirmed by data on the energy inten- sity to failure at this stage. Apparently, the formation of large dimples could not be realized due to local plastic strains. The proportion of areas with micro-dimples de- creased on the fracture surfaces. At the same time, the fraction of the fracture surface characterized only by flattened cone pits and tunnel ones, as well as their dimensions, increased significantly. Crack propagation was accompanied by significant strains and delamination. At the same time, the fractographic relief of the areas between them looked rather brittle. This indicated the localization of the fracture processes Figure 6. Effect of the Charpy impact test temperature on the 09Mn2Si steel specimens (both as-received and after the HR processing). Metals 2021, 11, 352 11 of 28 3.6. Fractographic Analysis SEM micrographs of the fracture surfaces are shown in Figures 7 and 8. At the temperature of +20 ◦C, all steel specimens after the HR processing fractured in a rather ductile manner. This fact was evidenced by a noticeable change in their shapes over the cross-sections (Figure 7a–c). The fracture surfaces were characterized by developed shapes with strain localization zones, as well as local delamination. Apparently, depth of these layers largely determined the energy intensity to failure, and as mentioned above, was reflected in the shape of the shear lips [46]. • At the micro-level, rounded pits surrounded by shear fracture spots had been formed on the as-received steel specimen in area 1 (crack initiation, Figure 1g)). In general, it had a ductile-brittle pattern (Figure 7a,g). The fracture surface was quite flat, formed according to the ductile-brittle mechanism. This was confirmed by the presence of multiple “chips”, which alternated with areas of ductile fracture of a dimple type. In area 2 (crack propagation, Figure 1m)), the surface morphology was very rough that indicated a clear heterogeneity in the strain process development at the main crack propagation (failure) stage (Figure 7s). It had been formed by large elongated dimples oriented across the crack propagation direction, and also included the large zone 4 (Figure 7s), corresponded to the specimen failure according to Figure 1g, which had a smoother surface with a gradual bend towards the point of impact loading. • After the HR, the crack initiation mechanism was predominantly ductile-brittle, with a minor brittle component (Figure 7b,c,h,i). Dimples had deformed profiles oriented along the crack front direction that indicated an increase in the material ductility and its sufficiently high crack resistance. This was confirmed by data on the energy intensity to failure at this stage. Apparently, the formation of large dimples could not be realized due to local plastic strains. The proportion of areas with micro-dimples decreased on the fracture surfaces. At the same time, the fraction of the fracture surface characterized only by flattened cone pits and tunnel ones, as well as their dimensions, increased significantly. Crack propagation was accompanied by significant strains and delamination. At the same time, the fractographic relief of the areas between them looked rather brittle. This indicated the localization of the fracture processes within individual micro-zones (Figure 7k,l). Area 4 (Figure 7t,u) had a characteristic rough relief, without pronounced distinctions from a similar fracture of the as-received steel specimen. Metals 2021, 11, 352 12 of 28 Metals 2021, 11, 352 12 of 28 rough relief, without pronounced distinctions from a similar fracture of the as-re- ceived steel specimen. As-received After the three-pass HR After the five-pass HR (a) (b) (c) Area 1 (d) (e) (f) (g) (h) (i) Area 2 (j) (k) (l) (m) (n) (o) Figure 7. Cont. Metals 2021, 11, 352 13 of 28 Metals 2021, 11, 352 13 of 28 Area 3 (p) (q) (r) (s) (t) (u) Figure 7. Fracture surfaces of the 09Mn2Si steel after the Charpy impact tests at the temperature of +20 °C. (explanations for a–u are given in the text above). At the temperature of −70 °C, the fracture surface macro-irregularities were observed near the crack initiation area (Figure 8), which indicated the presence of local strains dur- ing specimen failure. The shapes of the specimen transverse cross-sections were deformed at the entire investigated temperature range, which reflected high energy intensity to failure. At the micro-level, crack initiation and propagation areas had been formed due to shear strains for the as-received steel specimens (Figure 8g,m). The fracture surfaces had a smoothed ‘structureless’ appearance with several micro-chips. The fracture mechanism had been mixed with the predominant realization of micro-chips at the boundaries of de- formed micro-zones. The fractographic relief reflected the ductile-brittle pattern of frac- ture. At the failure area (Figure 8s), where the sharp drop in load had been previously found in the load-deflection curves (Figure 3c,d), a brittle fracture surface was observed. The specimens after the HR retained shear fracture pattern. Ductile tear dimples were absent (Figure 8b,c). The surfaces had an even more ‘blurry’ and ‘structureless’ character (Figure 8h,i,n,o). No significant variations were found in the fractographic surface pattern for the specimens after both three- and five-pass HR. This was also consistent with the data presented above on the study of the fracture macro-characteristics. The main differ- ence from the as-received steel specimens was the ductile fracture type at all areas, includ- ing the failure one (Figure 8t,u). As-received After the three-pass HR After the five-pass HR (a) (b) (c) Area 1 Figure 7. Fracture surfaces of the 09Mn2Si steel after the Charpy impact tests at the temperature of +20 ◦C. (explanations for a–u are given in the text above). At the temperature of −70 ◦C, the fracture surface macro-irregularities were observed near the crack initiation area (Figure 8), which indicated the presence of local strains during specimen failure. The shapes of the specimen transverse cross-sections were deformed at the entire investigated temperature range, which reflected high energy intensity to failure. At the micro-level, crack initiation and propagation areas had been formed due to shear strains for the as-received steel specimens (Figure 8g,m). The fracture surfaces had a smoothed ‘structureless’ appearance with several micro-chips. The fracture mechanism had been mixed with the predominant realization of micro-chips at the boundaries of deformed micro-zones. The fractographic relief reflected the ductile-brittle pattern of fracture. At the failure area (Figure 8s), where the sharp drop in load had been previously found in the load-deflection curves (Figure 3c,d), a brittle fracture surface was observed. The specimens after the HR retained shear fracture pattern. Ductile tear dimples were absent (Figure 8b,c). The surfaces had an even more ‘blurry’ and ‘structureless’ character (Figure 8h,i,n,o). No significant variations were found in the fractographic surface pattern for the specimens after both three- and five-pass HR. This was also consistent with the data presented above on the study of the fracture macro-characteristics. The main difference from the as-received steel specimens was the ductile fracture type at all areas, including the failure one (Figure 8t,u). Metals 2021, 11, 352 14 of 29 Area 3 (p) (q) (r) (s) (t) (u) Figure 7. Fracture surfaces of the 09Mn2Si steel after the Charpy impact tests at the temperature of +20 °C. (explanations for a–u are given in the text above). At the temperature of −70 °C, the fracture surface macro-irregularities were ob- served near the crack initiation area (Figure 8), which indicated the presence of local strains during specimen failure. The shapes of the specimen transverse cross-sections were deformed at the entire investigated temperature range, which reflected high en- ergy intensity to failure. At the micro-level, crack initiation and propagation areas had been formed due to shear strains for the as-received steel specimens (Figure 8g,m). The fracture surfaces had a smoothed ‘structureless’ appearance with several micro-chips. The fracture mechanism had been mixed with the predominant realization of micro-chips at the boundaries of de- formed micro-zones. The fractographic relief reflected the ductile-brittle pattern of frac- ture. At the failure area (Figure 8s), where the sharp drop in load had been previously found in the load-deflection curves (Figure 3c,d), a brittle fracture surface was observed. The specimens after the HR retained shear fracture pattern. Ductile tear dimples were absent (Figure 8b,c). The surfaces had an even more ‘blurry’ and ‘structureless’ character (Figure 8h,i,n,o). No significant variations were found in the fractographic surface pattern for the specimens after both three- and five-pass HR. This was also consistent with the data presented above on the study of the fracture macro-characteristics. The main differ- ence from the as-received steel specimens was the ductile fracture type at all areas, includ- ing the failure one (Figure 8t,u). As-received After the three-pass HR After the five-pass HR (a) (b) (c) Figure 8. Cont. Metals 2021, 11, 352 14 of 28 Metals 2021, 11, 352 15 of 29 Area 1 (d) (e) (f) (g) (h) (i) Area 2 (j) (k) (l) (m) (n) (o) Area 3 (p) (q) (r) Figure 8. Cont. Metals 2021, 11, 352 15 of 28 Metals 2021, 11, 352 16 of 29 (s) (t) (u) Figure 8. Fracture surfaces of the 09Mn2Si steel after the Charpy impact tests at the temperature of –70 °C. (explanations for a–u are given in the text above). In this paper, it was shown that the increase in the number of the HR passes from three to five gave a small effect on increasing impact toughness of the 09Mn2Si HSLA steel, which was associated with the exhaustion of the material ductility margin. Analyz- ing the data in Figure 8, it was possible to note the similarity of the fracture mechanisms after the three- and five-pass HR processing. It could be assumed that a further increase in the number of passes would result in embrittlement of the steel and a decrease in its toughness [47]. 4. Justification of the Theoretical Approach to the Analysis of Deformation Processes Under Impact Loading The effects revealed in the experimental section required complex analysis, including the use of computer simulation methods. A multilevel approach for simulation of defor- mation and fracture processes seemed promising [48–50]. However, before solving the issue of combining mechanisms operating at different scale levels (solved, for example, within the framework of molecular dynamics, crystal plasticity, finite element modeling, etc.), it was necessary to analyze the contribution of the processes at each scale. In this regard, it was expedient to apply the excitable cellular automata method at the meso-level, which showed its effectiveness for analyzing the behavior of materials under intense loads, considering rotational strain modes and material fragmentation at the pre-fracture stage. In this paper, the authors presented how this method could be developed to con- sider this period in terms of micro-damage initiation. In the previous paper [51], the authors—within the framework of the excitable cellu- lar automata method [52–54]—implemented a model for the numerical description of the process of material fatigue under cyclic loading. This model used Mott’s theory [55], which postulated the dominance of kinetic (thermal fluctuation) processes when consid- ering the dynamics of sliding of crystallite planes. According to this, the general expres- sion for the shear activation energy E, taking into account the external stress σ, was written as E = nF ± σnωa/2, (1) where a was an interatomic distance; n was the number of atoms in a cluster of the crys- tallographic plane, the area of which was nω; F was the free energy required for shear due to cluster disordering. In this case, the expression for the free energy was written as F = nL (1 − T/TM), (2) where TM was the melting point; L was latent heat of fusion per atom. Mott’s theory described inter-crystalline shear as an oscillatory process directly re- lated to thermal vibrations of the atomic lattice. Thus, for example, the frequencies of clus- ter shifts by distance a could be written as = nU± ⁄ /kT (3) Figure 8. Fracture surfaces of the 09Mn2Si steel after the Charpy impact tests at the temperature of −70 ◦C. (explanations for a–u are given in the text above). In this paper, it was shown that the increase in the number of the HR passes from three to five gave a small effect on increasing impact toughness of the 09Mn2Si HSLA steel, which was associated with the exhaustion of the material ductility margin. Analyzing the data in Figure 8, it was possible to note the similarity of the fracture mechanisms after the three- and five-pass HR processing. It could be assumed that a further increase in the number of passes would result in embrittlement of the steel and a decrease in its toughness [47]. 4. Justification of the Theoretical Approach to the Analysis of Deformation Processes under Impact Loading The effects revealed in the experimental section required complex analysis, including the use of computer simulation methods. A multilevel approach for simulation of defor- mation and fracture processes seemed promising [48–50]. However, before solving the issue of combining mechanisms operating at different scale levels (solved, for example, within the framework of molecular dynamics, crystal plasticity, finite element modeling, etc.), it was necessary to analyze the contribution of the processes at each scale. In this regard, it was expedient to apply the excitable cellular automata method at the meso-level, which showed its effectiveness for analyzing the behavior of materials under intense loads, considering rotational strain modes and material fragmentation at the pre-fracture stage. In this paper, the authors presented how this method could be developed to consider this period in terms of micro-damage initiation. In the previous paper [51], the authors—within the framework of the excitable cellular automata method [52–54]—implemented a model for the numerical description of the process of material fatigue under cyclic loading. This model used Mott’s theory [55], which postulated the dominance of kinetic (thermal fluctuation) processes when considering the dynamics of sliding of crystallite planes. According to this, the general expression for the shear activation energy E, taking into account the external stress σ, was written as E = nF ± σnωa/2, (1) where a was an interatomic distance; n was the number of atoms in a cluster of the crystallographic plane, the area of which was nω; F was the free energy required for shear due to cluster disordering. In this case, the expression for the free energy was written as F = nL (1 − T/TM), (2) where TM was the melting point; L was latent heat of fusion per atom. Mott’s theory described inter-crystalline shear as an oscillatory process directly related to thermal vibrations of the atomic lattice. Thus, for example, the frequencies of cluster shifts by distance a could be written as ς1 = ς0 e−(nU±σ nω a/2)/kT (3) Metals 2021, 11, 352 16 of 28 where ς0 was the frequency of atomic vibrations. Thus, internal stresses relaxed in a crystal at shear areas of the atomic planes. On the other hand, the mechanical energy of an external load enabled to dissipate due to thermal fluctuation processes associated with these shears. In this paper, the authors, adhering to the line of dominance of the kinetic processes at the micro- and mesoscale levels, supplemented the excitable cellular automata model in the part of the thermo-fluctuation theory of strength, first proposed by S.N. Zhurkov in [56] and greatly developed later [57–60]. Considering the structure-dependent kinetic processes was also very important to establish the dependence of strength characteristics from temperature and loading rate upon material fracture. These tasks were relevant when solving, for example, such a challenge as low-temperature brittleness [61]. In the kinetic concept of S.N. Zhurkov, fracture was considered as a process that began immediately after application of an external load to a body and continuously developed in time. According to this theory, duration to failure τ was determined by the following relation: τ = τ0 e(U0−χσ)/kT, (4) where k was Boltzmann′s constant; σ was tensile stress; T was temperature; χ, τ0, and U0 were material constants that had a clearly defined physical meaning. Thus, τ0 was duration of natural vibrations of solid atoms (about 10−13 s); U0 was the level of the energy barrier that determined the probability of rupture of the interatomic bond. It could be concluded that τ0 was the reciprocal of ς0 in expression (3) by S.N. Zhurkov, who described the frequency of an elementary act of shear strain in Mott’s theory. Thus, frequency of natural vibrations of the crystal lattice was a fundamental quantity both in describing plastic strain and fracture of a solid. The parameter χ corresponded to a wide range of defect microstructures nucleated in a solid due to fracture under external loads and strongly depended on the material. It was found for the case of metallic single crystals that this coefficient was strictly inversely proportional to the density of dislocations generated under the external load. For pure polycrystalline metals, this value was determined by the density of dislocations at the grain boundaries. Within the framework of the developed cellular automata approach, the material ‘torque’ (rotation) occurred under a local moment of force in an active element at a certain angle. The crystal lattice was distorted under this rotation, which was accompanied by the formation of single dislocations, and, then, multiple slip lines. Calculation of the material torsion angle enabled to quantify the induced change in density. For this purpose, the authors applied the approximation of the rod torsion model, considered in detail in [53,54]. Thus, χ, being the reciprocal of the density of dislocations, could be calculated proceeding from the local moment of forces by the formula: χ = klat ∆ρd = klatGπ r SCA M , (5) where klat was the proportionality coefficient associated with the type of crystal lattice; r was the radius of an active element; G was shear modulus; SCA was the cross-sectional area of the cellular automaton active element. The action of local moments of forces in the material was the reason for the transition of a cellular automaton active element from the initial ‘elasticity’ state to the ‘pre-failure’ one. The probability of this transition P1 was calculated for the element with index i as P1 = J ∑ j=1 ∣∣∣∣ →Mi − → Mj ∣∣∣∣2(∣∣∣∣ →Mi ∣∣∣∣+ ∣∣∣∣ →Mj ∣∣∣∣)2 , (6) Metals 2021, 11, 352 17 of 28 where j was the local index of a neighboring element on the first coordination sphere of the element with index i; J was the number of such neighboring elements (in the case of FCC packing J = 12). Before the beginning of the rapid fracture stage, stationary accumulation of cracks occurred when the rate N was constant, and the number of cracks tended to a certain critical (threshold) value N* (the number of cracks per unit volume). The critical value was determined from the concentration law [62] N*−1/3/L = K*, (7) where L was a mean crack length; K* was the dimensionless material parameter in the range from 2.1 to 6.0. It was easy to see that K* had the meaning of the cubic root of the quantity, which was inversely proportional to the critical volume fraction of all cracks in the material. Within the framework of the excitable cellular automata method, the last relation could be used to calculate the probability of the transition of an active element from the ‘pre-fracture’ state to the ‘crack’ one (P2). This was reflected in the laws of interaction with neighboring active elements. This probability was defined as the ratio of the density (amount per unit volume) of accumulated microcracks to the critical density: P2 = N N∗ . (8) Furthermore, in terms of a cellular-automaton evolutionary scheme, the probability was expressed as a sum over time steps of the algorithm as P2 = . N0K3∆t t ∑ 0 e−(U0−χσ)/kT, (9) The authors believed that the modification of the steel microstructure by the HR processing contributed to a longer stage of stationary damage accumulation, delaying the loss of the material’s bearing capacity. The degree of structural ‘imperfection’ of the material was effectively considered in the theory of S.N. Zhurkov by means of the χ parameter, which was determined by the density of dislocations at the crystallite level and by the density and energy state of the internal interfaces at the macro-scale one. Thus, within the framework of the described method, it was possible to consider both the temporal features of the damage accumulation process (according to the theory of S.N. Zhurkov) and the large-scale one (analyzed using the physical mesomechanic postulates), including the mesoscopic consideration of the deformation processes in Mott′s theory. The method explicitly took into account the energy dissipation during the nucleation of defect structures in accordance with the relation [54] ∆Ed = kdissG ∣∣∣∆→γ ∣∣∣2π r3 c 2 , (10) where ∆ → γ was an increment of the vector torsion angle of the material in the active element; G was shear modulus; kdiss was the experimentally measured dissipation coefficient; rc was the radius of the active element. The authors suggested the development of the presented approach of cellular au- tomata modeling in determining the critical tensile strength in order to use its values for simulation processes at higher scale levels. 5. Numerical Simulation Results As noted in the experimental section, resistance to impact fracture was significantly affected by the processes at the macro-discontinuity propagation stage. However, the Metals 2021, 11, 352 18 of 28 microstructure formed upon the HR processing was decisive in terms of energy spent to crack initiation. The HR could have had a dual effect. On the one hand, an increase in yield stress restrained the transition of the material to plastic flow and, accordingly, the moment of crack initiation. A competing process was strain hardening (revealed, for example, by the microhardness tests) that reduced toughness. Together with the retention of strain dislocation mechanisms in the steel after the HR processing, it was possible to ensure the KCV values at the as-received steel level, while there was no significant decrease in impact toughness with the decrease in the test temperature compared to the as-received steel. Thus, it was necessary to study how the modification of the microstructure by the HR processing affected propagation of energy fluxes in the steel specimens under impact loading. This section was devoted to solving the issue. The authors noted that they did not try to achieve a complete analogy with the experimental conditions (as was typically done for the finite element method). Nevertheless, a comparison of the obtained results with the above-described experimental data enabled to both develop the simulation method and out- line ways to further improve the properties of structural materials by their microstructure modification via the HR processing. Implementing the modified excitable cellular automata method, two computer simu- lation series were carried out to study in detail impact loading of the steel specimens with different microstructures characterized by grain sizes of 30 × 6 × 6 (grains were elongated along the x-axis), 6 × 6 × 30 (grains were elongated along the z-axis), 9 × 10 × 12 and 3 × 3 × 3 µm (Figure 9). Each specimen had dimensions of 90 × 30 × 60 µm and a symmetrical V-shaped notch 15 µm wide and 13 µm deep in the middle of the lower face. The size of a cellular automaton element was 1.5 µm, the time step was 1 ns. In the first simulation series, the temperature of the specimens was +20 ◦C, but it was −70 ◦ C in the second ones. During each simulation with a duration of 100 µs, an elliptical area with semi-axes of 30 and 10 µm and a center in the middle of the upper face was under impact loading, and the smooth distribution of compression strain was set so that its maximum value (100 s−1) was reached in the ellipse center, and the minimum (zero) was on its boundary. The smooth load distribution on the lower face simulated the support reaction: the value of the compression strain rate was maximal (10 s−1) at the face edges, but it was minimal (2 s−1) at the face center. The rest of the specimen faces were not loaded and did not exchange energy with environment. Metals 2021, 11, 352 18 of 28 decrease in impact toughness with the decrease in the test temperature compared to the as-received steel. Thus, it was necessary to study how the modification of the microstructure by the HR processing affected propagation of energy fluxes in the steel specimens under impact loading. This section was devoted to solving the issue. The authors noted that they did not try to achieve a complete analogy with the experimental conditions (as was typically done for the finite element method). Nevertheless, a comparison of the obtained results with the above-described experimental data enabled to both develop the simulation method and outline ways to further improve the properties of structural materials by their microstructure modification via the HR processing. Implementing the modified excitable cellular automata method, two computer simulation series were carried out to study in detail impact loading of the steel specimens with different microstructures characterized by grain sizes of 30 × 6 × 6 (grains were elongated along the x-axis), 6 × 6 × 30 (grains were elongated along the z-axis), 9 × 10 × 12 and 3 × 3 × 3 μm (Figure 9). Each specimen had dimensions of 90 × 30 × 60 μm and a symmetrical V-shaped notch 15 μm wide and 13 μm deep in the middle of the lower face. The size of a cellular automaton element was 1.5 μm, the time step was 1 ns. In the first simulation series, the temperature of the specimens was +20 °C, but it was –70 ° C in the second ones. During each simulation with a duration of 100 μs, an elliptical area with semi-axes of 30 and 10 μm and a center in the middle of the upper face was under impact loading, and the smooth distribution of compression strain was set so that its maximum value (100 s–1) was reached in the ellipse center, and the minimum (zero) was on its boundary. The smooth load distribution on the lower face simulated the support reaction: the value of the compression strain rate was maximal (10 s–1) at the face edges, but it was minimal (2 s–1) at the face center. The rest of the specimen faces were not loaded and did not exchange energy with environment. (a) (b) (c) (d) Figure 9. Microstructures of the steel specimens with grain sizes: 30 × 6 × 6 μm (a), 6 × 6 × 30 μm (b), 9 × 10 × 12 μm (c) и 3 × 3 × 3 μm (d). Figures 10–12 show plots of time dependences of specific elastic energy, specific dissipated torsion energy and their proportions in the total energy, as well as mean moduli of specific moment of force for each of the four considered specimens. The expressions, according to which the calculations of these parameters were carried out, had been given in [52]. The behavior of the graphs in Figure 10 enabled to conclude that the rate of accumulation of elastic energy gradually decreased in all specimens, and the minimum values of elastic energy were reached in the specimen with fine grains (3 μm in size). Of the other three specimens, the highest value of the accumulated specific elastic energy was shown by the specimen with almost equiaxed grains of ~10 μm, and the smallest one possessed the specimen with horizontally oriented (elongated) grains. The curves shown in Figure 10b, had very significant distinctions in the specimen energy capacity with lowering temperature. For instance, the accumulated energy in the finest-grained specimen at the temperature of −70 °C was 30% lower than that for the same microstructure at +20 °C. It could also be noted that the specimen with the lamellar horizontal microstructure (30 × 6 × 6 μm3) showed a noticeable sensitivity to the decrease Figure 9. Microstructures of the steel specimens with grain sizes: 30 × 6 × 6 µm (a), 6 × 6 × 30 µm (b), 9 × 10 × 12 µm (c), 3 × 3 × 3 µm (d). Figures 10–12 show plots of time dependences of specific elastic energy, specific dissipated torsion energy and their proportions in the total energy, as well as mean moduli of specific moment of force for each of the four considered specimens. The expressions, according to which the calculations of these parameters were carried out, had been given in [52]. The behavior of the graphs in Figure 10 enabled to conclude that the rate of accumula- tion of elastic energy gradually decreased in all specimens, and the minimum values of elastic energy were reached in the specimen with fine grains (3 µm in size). Of the other three specimens, the highest value of the accumulated specific elastic energy was shown Metals 2021, 11, 352 19 of 28 by the specimen with almost equiaxed grains of ~10 µm, and the smallest one possessed the specimen with horizontally oriented (elongated) grains. The curves shown in Figure 10b, had very significant distinctions in the specimen energy capacity with lowering temperature. For instance, the accumulated energy in the finest-grained specimen at the temperature of −70 ◦C was 30% lower than that for the same microstructure at +20 ◦C. It could also be noted that the specimen with the lamellar horizontal microstructure (30 × 6 × 6 µm3) showed a noticeable sensitivity to the decrease in temperature compared with the ‘coarse-grained’ one. In this case, low temperature, which changed the energy state of the grain boundaries, caused a higher barrier to the elastic energy front propagation. Metals 2021, 11, 352 19 of 28 in temperature compared with the ‘coarse-grained’ one. In this case, low temperature, which changed the energy state of the grain boundaries, caused a higher barrier to the elastic energy front propagation. (a) (b) Figure 10. Graphs of time dependences of specific elastic energy for the specimens with grain sizes 30 × 6 × 6 (blue, 4), 6 × 6 × 30 (red, 2), 9 × 10 × 12 (green, 3) and 3 × 3 × 3 μm (black, 1) at the temperatures of +20 °C (a) and −70 °C (b). The value of the specific torsional energy characterized the ability of the material to nucleate new defect structures according to expression (10). Figure 11 shows a gradual increase in the growth rate of the specific torsional energy for all specimens. At the same time, the maximum values were for the fine-grained material during the entire simulation. From three other specimens, they were for the material with grains oriented across the load direction. The minimum values of the torsional energy were accumulated in the specimen with large equiaxed grains (highly likely, due to the small extent of the boundaries themselves). At the end of the loading stage, the values of the torsional energy fractions stabilized in the range of 85–95%. The difference was noticeable only for the fine- grain specimens. (a) (b) Figure 11. Graphs of time dependences of the proportion of the specific dissipated torsional energy in the total energy for the specimens with grain sizes 30 × 6 × 6 (blue, 4), 6 × 6 × 30 (red, 2), 9 × 10 × 12 (green, 3) and 3 × 3 × 3 μm (black, 1) at the temperatures of +20 °C (a) and −70 °C (b). Analysis of the behavior of the plots in Figures 10 and 11 enabled to conclude that the maximum dissipation of elastic energy occurred in the fine-grained microstructure. Analysis of the data on the specific torsional energy indicated that plastic strain developed according to the ‘shear plus rotation’ mechanism. In this case, the relaxation of excess stresses was in the grain bodies due to interplanar slip (described within the framework of Mott’s theory), while the rotational mechanism had to prevail at the grain boundaries. It should be noted that grain boundaries were not only the barrier preventing the elastic energy transfer and absorbing part of it, but they also significantly changed the stress– strain state patterns in the grain bodies. Of the three specimens with the same mean volume, but having different grain orientations, the highest dissipation of elastic energy was shown by the one with the preferred orientation of the boundaries across the load direction. This was due to the fact that the strain front was forced to overcome the area of the grain boundaries larger than Figure 10. Graphs of time dependences of specific elastic energy for the specimens with grain sizes 30 × 6 × 6 (blue, 4), 6 × 6 × 30 (red, 2), 9 × 10 × 12 (green, 3) and 3 × 3 × 3 µm (black, 1) at the temperatures of +20 ◦C (a) and −70 ◦C (b). The value of the specific torsional energy characterized the ability of the material to nucleate new defect structures according to expression (10). Figure 11 shows a gradual increase in the growth rate of the specific torsional energy for all specimens. At the same time, the maximum values were for the fine-grained material during the entire simulation. From three other specimens, they were for the material with grains oriented across the load direction. The minimum values of the torsional energy were accumulated in the specimen with large equiaxed grains (highly likely, due to the small extent of the boundaries themselves). At the end of the loading stage, the values of the torsional energy fractions stabilized in the range of 85–95%. The difference was noticeable only for the fine-grain specimens. Metals 2021, 11, 352 19 of 28 in temperature compared with the ‘coarse-grained’ one. In this case, low temperature, which changed the energy state of the grain boundaries, caused a higher barrier to the elastic energy front propagation. (a) (b) Figure 10. Graphs of time dependences of specific elastic energy for the specimens with grain sizes 30 × 6 × 6 (blue, 4), 6 × 6 × 30 (red, 2), 9 × 10 × 12 (green, 3) and 3 × 3 × 3 μm (black, 1) at the temperatures of +20 °C (a) and −70 °C (b). The value of the specific torsional energy characterized the ability of the material to nucleate new defect structures according to expression (10). Figure 11 shows a gradual increase in the growth rate of the specific torsional energy for all specimens. At the same time, the maximum values were for the fine-grained material during the entire simulation. From three other specimens, they were for the material with grains oriented across the load direction. The minimum values of the torsional energy were accumulated in the specimen with large equiaxed grains (highly likely, due to the small extent of the boundaries themselves). At the end of the loading stage, the values of the torsional energy fractions stabilized in the range of 85–95%. The difference was noticeable only for the fine- grain specimens. (a) (b) Figure 11. Graphs of time dependences of the proportion of the specific dissipated torsional energy in the total energy for the specimens with grain sizes 30 × 6 × 6 (blue, 4), 6 × 6 × 30 (red, 2), 9 × 10 × 12 (green, 3) and 3 × 3 × 3 μm (black, 1) at the temperatures of +20 °C (a) and −70 °C (b). Analysis of the behavior of the plots in Figures 10 and 11 enabled to conclude that the maximum dissipation of elastic energy occurred in the fine-grained microstructure. Analysis of the data on the specific torsional energy indicated that plastic strain developed according to the ‘shear plus rotation’ mechanism. In this case, the relaxation of excess stresses was in the grain bodies due to interplanar slip (described within the framework of Mott’s theory), while the rotational mechanism had to prevail at the grain boundaries. It should be noted that grain boundaries were not only the barrier preventing the elastic energy transfer and absorbing part of it, but they also significantly changed the stress– strain state patterns in the grain bodies. Of the three specimens with the same mean volume, but having different grain orientations, the highest dissipation of elastic energy was shown by the one with the preferred orientation of the boundaries across the load direction. This was due to the fact that the strain front was forced to overcome the area of the grain boundaries larger than Figure 11. Graphs of time dependences of the proportion of the specific dissipated torsional energy in the total energy for the specimens with grain sizes 30 × 6 × 6 (blue, 4), 6 × 6 × 30 (red, 2), 9 × 10 × 12 (green, 3) and 3 × 3 × 3 µm (black, 1) at the temperatures of +20 ◦C (a) and −70 ◦C (b). Analysis of the behavior of the plots in Figures 10 and 11 enabled to conclude that the maximum dissipation of elastic energy occurred in the fine-grained microstructure. Analysis of the data on the specific torsional energy indicated that plastic strain developed according to the ‘shear plus rotation’ mechanism. In this case, the relaxation of excess stresses was in the grain bodies due to interplanar slip (described within the framework Metals 2021, 11, 352 20 of 28 of Mott’s theory), while the rotational mechanism had to prevail at the grain boundaries. It should be noted that grain boundaries were not only the barrier preventing the elastic energy transfer and absorbing part of it, but they also significantly changed the stress–strain state patterns in the grain bodies. Of the three specimens with the same mean volume, but having different grain orientations, the highest dissipation of elastic energy was shown by the one with the preferred orientation of the boundaries across the load direction. This was due to the fact that the strain front was forced to overcome the area of the grain boundaries larger than in other specimens during the propagation process. On the other hand, the ‘vertical’ grain orientation also increased the degree of energy dissipation compared to the equiaxial structure. In this case, the reason was the elastic energy front propagation to the support areas on the lower face, which passed mainly through the grain boundaries oriented parallel to its direction. Thus, there were two competing mechanisms of elastic energy absorption: (i) the microstructure with horizontally oriented grains realized micro-plastic shears in the grain bodies (according to Mott’s theory); (ii) the microstructure with vertically oriented grains was subjected by the intensive development of rotational wave of defect flows along the grain boundaries (according to the theory of V.E. Panin—V.E. Egorushkin). The general tendency towards stabilization of the proportions of elastic and torsional energy for all types of the grain microstructures and temperatures was due to the fact that there was a gradual decrease in the energy gradients and moments of forces both inside the grains and at their boundaries upon loading because of the establishment of an equilibrium regime for the specimen macro-volume. This effect was not only due to the balancing of elastic macro-stresses, but also because of the achievement of a certain threshold value of the torsional energy, when the material was saturated with dislocation microstructures and exhausted the ductility margin [63]. In this sense, the transition point on the curves, shown in Figure 11, to the horizontal section (‘saturation’) indicated the beginning of the pre-fracture stage, the features of which are discussed below. The magnitude of the modulus of the specific moment of force was analogous to the potential energy of torsion, which later (in the process of increasing the density of dislocations and the formation of slip lines) was transformed into the work spent to plastic strain. The time dependences presented in Figure 12, a show a gradually decaying increase in the mean absolute values of the moment of force. For a more detailed study of the processes of the curvature formation at the grain boundaries, time dependences were also built for the three components of the moment of force for all simulation specimens (Figure 12b–d). Since the relationship of the separate components of the moment of force was characterized exclusively by structural features, these graphs were plotted only for the temperature of +20 ◦C. The maximum mean values of the moduli of the moments of forces and its components were shown by the specimen with the fine-grained microstructure. Again, this fact was due to the largest area of the grain boundaries in comparison with other investigated mi- crostructures. The plots of the mean moduli of the moments of forces were almost identical for all other specimens with the same mean grain volume. Nevertheless, the dependences of the mean values of the moduli of various components of the moments were significantly different from each other. Accordingly, the grain shapes affected the spatial distribution pattern of the torsion zones of the material around one or another coordinate axis. For example, torsion of the material under simulation around the (horizontal) x axis was more intense during deformation if the grains were stretched in the horizontal direction. At the same time, ‘swirls’ around the vertical z axis was most pronounced in the microstructure with vertically oriented grains. This confirmed the above thesis about the competing mechanisms of the elastic energy dissipation. In the case of predominantly horizontally oriented grains (across the load direction), the dissipation was due to the formation of defect microstructures in the grain bodies. In the case of the grain microstructure oriented coaxially to impact load, the predominant development of the rotational wave of defect flows was along the grain boundaries. Metals 2021, 11, 352 21 of 28 Metals 2021, 11, 352 21 of 28 (a) (b) (c) (d) Figure 12. Graphs of time dependences of the mean moduli of the specific moment of force (a) and its components X (b), Y (c), Z (d) for the specimens with grain sizes 30 × 6 × 6 (blue, 4), 6 × 6 × 30 (red, 2), 9 × 10 × 12 (green, 3) and 3 × 3 × 3 μm3 (black, 1) at the temperature of +20 °C. Figure 13 shows maps of the grain microstructure on the front faces of the simulation specimens and the distribution patterns of the specific elastic energy, the specific dissipated torsion energy, the mean moduli of the components of the specific moment of force, and the probabilities of micro-crack initiation on these faces after 20 μs. The elastic energy propagation dynamics, presented in Figure 13, demonstrated a gradual inflow of energy from the load application area at the upper face (as well as from the support side located at the lower face) into the specimens. The grain microstructure had a significant effect on the elastic energy front propagation rate (the grain boundaries slowed it down). For this reason, the ‘saturation of the material with elastic energy’ was slowest in the fine-grained microstructure compared to the other ones. The elastic energy dissipation was most intense both at the point of the impact load application and near the supports at the lower face. It should be noted that the greatest increase in the torsional energy was in the regions of the grain boundaries, behind which this energy propagated deep into the specimen, mainly along the grain boundaries. Such a character of the torsional energy front propagation and its dissipation at the nucleation of the defect microstructures was caused by the action of the moments of forces arose as a result of the elastic energy gradients at the grain boundaries. In this case, the action of the local moments of forces also spread into the grain bodies, which subsequently caused their fragmentation and transcrystalline failure. Figure 12. Graphs of time dependences of the mean moduli of the specific moment of force (a) and its components X (b), Y (c), Z (d) for the specimens with grain sizes 30 × 6 × 6 (blue, 4), 6 × 6 × 30 (red, 2), 9 × 10 × 12 (green, 3) and 3 × 3 × 3 µm3 (black, 1) at the temperature of +20 ◦C. Figure 13 shows maps of the grain microstructure on the front faces of the simulation specimens and the distribution patterns of the specific elastic energy, the specific dissipated torsion energy, the mean moduli of the components of the specific moment of force, and the probabilities of micro-crack initiation on these faces after 20 µs. The elastic energy propagation dynamics, presented in Figure 13, demonstrated a gradual inflow of energy from the load application area at the upper face (as well as from the support side located at the lower face) into the specimens. The grain microstructure had a significant effect on the elastic energy front propagation rate (the grain boundaries slowed it down). For this reason, the ‘saturation of the material with elastic energy’ was slowest in the fine-grained microstructure compared to the other ones. The elastic energy dissipation was most intense both at the point of the impact load application and near the supports at the lower face. It should be noted that the greatest increase in the torsional energy was in the regions of the grain boundaries, behind which this energy propagated deep into the specimen, mainly along the grain boundaries. Such a character of the torsional energy front propagation and its dissipation at the nucleation of the defect microstructures was caused by the action of the moments of forces arose as a result of the elastic energy gradients at the grain boundaries. In this case, the action of the local moments of forces also spread into the grain bodies, which subsequently caused their fragmentation and transcrystalline failure. Metals 2021, 11, 352 22 of 28 Metals 2021, 11, 352 22 of 28 S T R U C T Eel Ed Mx My Mz Pd +20 Pd −70 (a) (b) (c) (d) Figure 13. Specimen grain microstructure maps and the patterns of spatial distributions of specific elastic energy (Eel), specific torsional energy (Ed), components of the moment of force and the probability of micro-crack initiation at the temperatures of +20 °С and −70 °С after 20 μs. (explanations for a–d are given in the text above). In Figure 13, the distribution patterns of the X component absolute values of the specific moment of force M were in many respects similar to the same distributions of the specific torsional energy E. A distinctive feature of these patterns was a more significant difference between the fine-grained and other microstructures in the degree of accumulation of the X moment component values near the support areas. Figure 13. Specimen grain microstructure maps and the patterns of spatial distributions of specific elastic energy (Eel), specific torsional energy (Ed), components of the moment of force and the probability of micro-crack initiation at the temperatures of +20 ◦C and −70 ◦C after 20 µs. (explanations for a–d are given in the text above). In Figure 13, the distribution patterns of the X component absolute values of the specific moment of force M were in many respects similar to the same distributions of the specific torsional energy E. A distinctive feature of these patterns was a more significant Metals 2021, 11, 352 23 of 28 difference between the fine-grained and other microstructures in the degree of accumulation of the X moment component values near the support areas. It was noted above that the Y component values of the moment of force reached their maximum values on the lateral specimen faces at the beginning stage of impact loading. Then, they were gradually increased within the specimen. As a result, quasi-symmetric patterns of these component distributions were found at the front faces of the specimens. Their type was determined by the shapes and the sizes of the grains. In the case of the Z component of the moment of force, the distribution patterns were also quasi-symmetric about the central vertical axis (the axis of the applied impact load). In this case, they gradually ‘stabilized’ as they were loaded. The maximum modulus values were achieved in the grain bodies located near the loaded upper and lower faces. This was due to the fact that the grain bodies accumulated the torsional energy as a result of the action of the moments of forces generated along the grain boundaries and the specimen surface. On the other hand, the grain boundaries were channels for the propagation of rotational wave flows deep into the material. The probability distributions of micro-crack initiation, calculated by formula (6), showed that the shapes and the sizes of the grains also had a significant effect on the nature of micro-cracking (Figure 13; Pd). Thus, there was an extremely effective damping effect of the microstructure with the preferred grain orientation perpendicular to the loading axis. Despite the “potentially in