Next Article in Journal
Failure Analysis and Reliability of Ni–Ti-Based Dental Rotary Files Subjected to Cyclic Fatigue
Next Article in Special Issue
On the Mn–C Short-Range Ordering in a High-Strength High-Ductility Steel: Small Angle Neutron Scattering and Ab Initio Investigation
Previous Article in Journal
Low-Cycle Fatigue Behaviour of AISI 18Ni300 Maraging Steel Produced by Selective Laser Melting
Previous Article in Special Issue
Molecular Dynamics Simulation of Crack Propagation in Nanoscale Polycrystal Nickel Based on Different Strain Rates
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ab Initio-Based Modelling of the Yield Strength in High-Manganese Steels

Steel Institute, RWTH Aachen University, Intzestraße 1, 52072 Aachen, Germany
*
Author to whom correspondence should be addressed.
Metals 2018, 8(1), 34; https://doi.org/10.3390/met8010034
Submission received: 15 November 2017 / Revised: 16 December 2017 / Accepted: 3 January 2018 / Published: 5 January 2018
(This article belongs to the Special Issue First-Principles Approaches to Metals, Alloys, and Metallic Compounds)

Abstract

:
An ab initio-based model for the strength increase by short-range ordering of C-Mn-Al clusters has been developed. The model is based on ab initio calculations of ordering energies. The impact of clusters on the yield strength of high-manganese austenitic steels (HMnS) is highly dependent on the configurational structure of the cells that carbon atoms will position themselves as interstitial atoms. The impact of the alloying elements C, Mn, and Al on the potential and actual increase in yield strength is analyzed. A model for the calculation of yield strengths of HMnS is derived that includes the impact of short-range ordering, grain size refinement, and solid solution strengthening. The model is in good agreement with experimental data and performs better than other models that do not include strengthening by short-range ordering.

1. Introduction

High Mn austenitic steels (HMnS) offer great combinations of hardness and ductility with a wide range of possible applications in various industries, such as the automotive industry [1,2,3]. Depending on the stacking fault energy, which is a result of alloying content, temperature, deformation grade, and other factors, HMnS exhibit a variety of different active deformation modes [4,5,6,7]. Those deformation modes include transformation induced plasticity (TRIP), twinning-induced plasticity (TWIP) or pronounced planar glide in combination with the formation of dislocation cells that is sometimes referred to as microband-induced plasticity (MBIP) [1,8,9,10]. The mode of deformation significantly influences the level of work hardening and plays a major role in determining the mechanical properties [7]. One of those material characteristics is the material’s yield strength (YS). In the past, various attempts have been made to model the YS of HMnS based on the chemical composition using linear regression [5,11,12,13]. However, these models often imply negative solid-solution strengthening by some of the alloying elements, like Mn, that are usually attributed to positive strengthening effects [2,14,15].
Interactions of carbon-metal pairs [16,17,18] and clusters are described by the concept of short-range ordering (SRO) phenomena [19] in HMnS and might offer an explanation for an additional increase in strength that might lead to a defective linear regression when using standard methods to calculate the impact of solid solution strengthening [14,20,21]. Short-range ordered clusters are a result of local differences in ordering energies that make particular unit cell configurations preferential to inhibit a carbon atom [22]. Ab initio calculations employing density-functional theory (DFT) have been shown to be suitable for calculations of ordering energies in HMnS for several chemical compositions. For example, thermodynamic predictions of κ-phase stability [23], interface energies [24], stacking fault energies [25], and hydrogen trap sites [26] showed good agreement with the experimental data.
Cell models [27] have been used to determine the location of carbon atoms in relation to the matrix or substitutional alloying atoms. Other studies show that the interaction of carbon atoms with each other is very repulsive and, thus, influences the diffusion behavior [4,25,28,29,30], which also points to a coordination of carbon atoms in octahedral sites dependent of surrounding solute atoms. Additionally, the position of C also influences the stacking fault energy of the material [31].
Short-range ordering has also been used to explain the pronounced planar glide that HMnS exhibit because of glide plane softening that occurs when short-range ordered C-Mn clusters are destroyed by dislocation stresses reorienting the carbon atom from its octahedral position to a tetrahedral position. This phenomenon also explains the serrated flow curves that can be found during tensile testing of HMnS, as the back stresses of dislocation pile-ups lead to an increase in strength. The reorientation of the carbon atom then leads to a temporary destruction of the short-range ordered cluster and as a result, a rapid drop in strength because of subsequent softening of this glide plane [32]. However, short-range ordering can be restored by short-range jumps of C back from the tetrahedral site to its initial octahedral site [1,2]. Those short-range jumps might be activated by temperature, but already occur at low temperatures where diffusion processes are usually negligible [1,32].
In this study, a basic cell model for the nearest-neighboring atoms of a central octahedral carbon atom is expanded to include the corner atoms of an elementary cell of an austenitic face-centered cubic (FCC) unit cell, which can be described as the second-nearest neighbors of the central octahedral site in an austenitic cell [27]. Results of ab initio calculations are employed to show that the position of Mn and Al being coordinated as a nearest neighbor or second-nearest neighbor to a central carbon atom is of major importance to the energy level of the system [33]. It is assumed that this also holds true in the case of a central interstitial nitrogen atom [2,33,34] but, in this work, the authors will focus solely on carbon. Based on the calculated ordering energies, the yield strength of HMnS is a result of a combination of short-range ordered clusters, regular solid-solution hardening, grain size effects, and the base strength of an austenitic lattice [14,35].

2. Materials and Methods

Many models for the calculation of the yield strength of austenitic steels have been published in the past. Those models generally employ two mechanisms to describe the yield strength: a strengthening by grain size reduction, the Hall-Petch effect, and the alloying elements’ contribution to the material’s strength [15,36,37]. Many of those models have used databases with significant amounts of different alloys to employ a linear regression in order to get parameters that determine the effect of alloying elements on the yield strength of the austenitic steel. Assuming that SRO in HMnS influences the yield strength σ Y significantly, a term describing this effect has to be added to the equation [14], resulting in:
σ Y = σ 0 + σ S S + σ G S + σ S R O
where σ 0 is the base strength of a pure austenitic iron lattice and is assumed to be 90   MPa [14]. σ S S is the contribution of alloying elements to the change in yield strength relative to a pure austenitic lattice. Since this paper focuses on HMnS with the main alloying elements C, Mn, and Al, those elements are assumed to be the main contributors to solid solution strengthening and the contributions of other alloying elements are neglected. It is assumed that this does not lead to a large mistake in calculating σ S S because, apart from Mn and Al, no substitutional alloying element exceeds alloying contents of more than 1   wt .   % .
σ S R O is the change in yield strength stemming from the existence of short-range ordered regions in the material. Based on the work of [14], this contribution to yield strength can be calculated as:
σ S R O   =   M τ S R O = M ( E O r d e r b 3 ) = M ( E r a n d o m E S R O b 3 )
where M = 3.06 is the Taylor-factor for fcc with b = α γ 2 being the magnitude of the Burgers vector and the lattice parameter in austenite in Ȧ depending on the molar fractions x i of the alloying elements and temperature in K [38]:
α γ = ( 3.5780 + 0.033 x C   +   0.00095 x M n     0.0002 x N i +   0.0006 x C r   +   0.0056 x A l   +   0.0031 x M o   +   0.0018 x V ) × ( 1   +   2.065   ×   10 9   ×   ( T     300 ) )
Looking at Equation (3), the temperature dependence of the lattice parameter is assumed to be linear.
σ G S describes the strengthening by grain refinement in accordance with the Hall-Petch effect and can be formulated as:
σ G S   =   k y d 1 / 2
Even though k y depends on the cooling procedure and the carbon content [39], the value k y = 11.3   MPa   mm 1 / 2 seems to be a reasonable assumption for HMnS in general, based on the findings in the literature [1,2,14,39,40,41].
As can be seen in Figure 1, one unit cell consists of 14 atoms situated on the lattice positions of the FCC bracket and one central octahedral interstitial site, where a single carbon atom can be situated.
First-nearest neighbors (1NN) are the six direct neighbors of the carbon atom. In Figure 1, they are shown in green. Second-nearest neighbors (2NN) in this model are the grey atoms on the corners of the FCC lattice. Although third-nearest neighbors (3NN) can be included in first-principle calculations, a recent study has shown that neither 3NN Al atoms nor 3NN Mn atoms have a positive impact on the ordering energy [33]. To reduce the complexity of the calculations and limit the amount of possible cell configurations, the presence of 3NN alloying atoms is not accounted for in the analysis carried out in this paper. Instead, it is assumed that the material is the sum of FCC-unit cells that do not impact each other. This is obviously not true in reality, but for the purposes of this study, which is aimed at providing a basic model to connect ab initio calculations with experimental results, it is suitable, since it leads to reduced mathematical complexities. The term “cell configuration” is used to describe the combination of alloying atoms and matrix atoms in the cell. “Coordination sphere” is used throughout this paper to refer to the position of atoms relative to the central carbon atom.
The configurations of the cells based on the alloying contents for this 14-cell model are calculated according to [27]:
n k l = 14 ! ( 14     k     l ) ! k ! l ! θ F e ( 14     k     l ) θ M n k θ A l l
where n k l is the fraction of all cells that have a configuration with k Mn atoms and l Al atoms with θ i being the molar fraction of Fe, Mn, and Al, denying the existence of other substitutional alloying elements and ignoring C, since carbon is an interstitial alloying element. Since diffusion rates of substitutional atoms in austenitic steels are very low, it is assumed that the distribution of substitutional elements according to Equation (5) is fixed and that the carbon atom is distributed in accordance with a minimum ordering energy. For clusters of nitrogen and chromium, this approach has been shown to be a realistic approximation, since the presence of nitrogen does not negatively influence the homogeneous distribution of chromium [34]. Cells of type n k l that contain a carbon atom are classified as n C k l . It should be noted that the assumption that substitutional atoms are distributed according to minimum ordering energies with respect to the presence of a central carbon atom is a very strong assumption, which is unlikely to be true in reality. However, without this assumption, the calculation of ordering energies Δ E C k l is not possible since ab initio results are presently not available for those unfavorable configurations.
The ordering energy Δ E C k l of a cell with k Mn atoms, l Al atoms and a central carbon atom is the combination of two energy contributions in comparison to a FCC-unit cell consisting only of Fe atoms and a central interstitial carbon atom. First, the energy required to bring a carbon atom in solution Δ E S o l in a FCC-unit cell with respect to the presence of k Mn atoms and l Al atoms in the first two coordination spheres. The second factor is the change in energy level Δ E L a t resulting from a lattice distortion due to the presence of Mn and Al:
Δ E C k l = Δ E S o l   +   Δ E L a t
Ab initio calculations show that Al is preferentially situated in the second coordination sphere and Mn has a positive effect on the ordering energy level if it sits on a lattice position as a 1NN [22,25,33]. According to those studies, this means that the difference in ordering energy is so profound that the likelihood of a cell configuration that contains 1NN Al atoms instead of 2NN Al atoms or 2NN Mn atoms instead of 1NN Mn atoms is very low. Accordingly, for further calculations in this work and to reduce the model’s complexity, it is assumed that for any number of Mn atoms up to six, they will always be 1NN and for any number of Al atoms up to 8, they will always be 2NN.
Using the calculated ordering energies Δ E C k l from [33], the six cells with the lowest ordering energies are listed in Table 1 and shown in Figure 2.
In Equation (2), E r a n d o m describes the energy state of an unordered alloying system. In this case, a single carbon atom is situated in the central octahedral sites of a randomly chosen FCC-unit cell. This distribution is based on the assumption that there are no preferential octahedral sites, so neither the presence of Mn nor Al influences the position of carbon atoms. Thus, it can be assumed that, for the calculation of E r a n d o m , the likelihood of the presence of carbon is the same for each type of cell.
Using the random C distribution for E r a n d o m , E r a n d o m is calculated according to:
E r a n d o m = k = 0 14 l = 0 14 n C k l , r a n d o m Δ E C k l ,   for   k   +   l     14   and   k , l     0
n C k l , r a n d o m   =   n k l × θ C
Since n k l = 1 , this method of determining n C k l , r a n d o m ensures that carbon is distributed homogeneously in the material, as well as ensuring that every available carbon atom is situated in an octahedral site of a cell.
Contrary to the homogeneous distribution of carbon for E r a n d o m , the distribution of C for the calculation of E S R O is based on the calculated ordering energies [33]. Since the differences in ordering energy between different cell configurations are significant, as is shown in Table 1, it is assumed that carbon is always situated in the most favorable cell that has not already been assigned a carbon atom. Assuming that at higher temperatures, for example during recrystallization annealing, the diffusion of carbon as an interstitial alloying element is possible, carbon will diffuse to positions with minimal ordering energy in order to minimize the system’s Gibbs energy. This means that the cells selected for saturation with carbon are the cells with the lowest Δ E C k l until n C k l , S R O   =   θ C .
The calculation of ordering energies in the model presented within this work was performed via ab initio simulations based on quantum-mechanical structure optimizations employing density-functional theory (DFT) within the Vienna ab initio Simulation Package (VASP) [33]. The quantum-mechanical structure optimizations are performed based on calculations of lowest-energy states of the electron structure and thus, are conducted at 0 K. Results from ab initio approaches employing DFT have offered good results that are in line with experimental evidence. For example, Mössbauer spectroscopy [17], atom-probe tomography (APT) [42], in situ synchrotron X-ray diffraction measurements [23], and correlative TEM/APT approaches [43] have shown that results from ab initio simulations employing DFT are applicable to systems containing nano-sized κ-carbides that are usually formed during isothermal annealing at temperatures of 600 °C and above. Generally speaking, the effect of temperature changes on the lattice parameter and on intrinsic material properties that derive from electron effects like electrical conductivity is often approximated to be linear [44]. The serrated flow curves of many high-manganese TWIP steels usually occur at temperatures of up to about 200 °C. At higher temperatures, this dynamic strain aging-like effect is either very small or is not detected at all [2,7].
In order to model the temperature dependence of σ S R O , the ab initio results are fitted using a linear approximation of a decrease in σ S R O between 0 K and 500 K. For this purpose, it is assumed that the fact that no serrated flow curve is expected at 500 K can be interpreted as a result of no significant strength increase by short-range ordered clusters.
Based on those assumptions and equations, σ Y , c a l c in MPa was calculated using Equations (2), (3), and (7) for 13 different HMnS with known average grain sizes d m in μm, listed in Table 2. The alloying contents are given in wt .   % .
Experimental yield strength data was taken from the literature [5,14]. Additionally, quasistatic tensile tests have been carried out using a Zwick Z250 tensile testing machine of the AllroundLine by Zwick Roell (Ulm, Germany).
σ S S has been calculated using a linear approximation to explain the differences between experimental and calculated values for the yield strength using the following equation that is a direct result of Equation (1):
σ S S , e x p   =   σ Y , e x p     σ γ - F e     σ G S     σ S R O
This way, σ S S , e x p was determined for all steels in Table 2. A regression analysis was performed with respect to the alloying contents of C, Mn, and Al.
The regression analysis lead to the following equation to calculate σ S S with respect to the alloying contents χ i in weight percent:
σ S S = 84.55 χ C   +   1.35 χ Mn   +   0.44 χ Al
Overall, the lattice distortion effect of carbon is by far the strongest contribution to solid solution hardening. This is in accordance with calculations of the misfit of C in an octahedral site using the crystallographic correlation r OL r Fe = 0.41 . Comparing the size of the octahedral site to the size of a carbon atom shows that a carbon atom is approximately 8% bigger than the octahedral site, leading to a lattice distortion. Using the lattice parameter calculated according to Equation (3), this size difference can be calculated to approximately 0.03 nm. Comparing this value with the lattice distortion effects by Al and Mn in an austenitic steel according to [45], which are 0.000065   nm wt .   %   Mn and 0.00095   nm wt .   %   Al respectively, the difference in contributions between C and Mn and Al in Equation (10) appear to be reasonable.

3. Results and Discussion

Figure 3 shows that the calculations for the yield strength σ Y , c a l c using the proposed model show good agreement with the experimental data.
The calculations used in this work lead to differences between σ Y , c a l c and σ Y , e x p that are smaller than 10%. Those small differences between σ Y , c a l c and σ Y , e x p might be a result of a negligence of other alloying elements that participate in solid solution hardening, minor inaccuracies of the measurement of the alloying contents, or d m .
Since σ S R O can also be understood as a deviation from linear solid solution hardening behavior, a close relationship between σ S S and σ S R O seems logical. Higher alloying contents of C, Mn, and Al lead to a higher solid solution strengthening, but also to more carbon atoms integrated in Mn-C, Al-C, and Mn-Al-C clusters, as well as a higher amount of cells with more Al atoms in the second coordination sphere or Mn atoms in the first coordination sphere according to Equation (5).
The individual contributions of solid solution hardening σ S S and σ S R O are illustrated in Figure 4.
The strengthening effect of short-range ordered clusters σ S R O is dependent on the chemical composition. For the steel compositions analyzed within the scope of this work, σ S R O at room temperature is approximately in the range of 20 MPa for X30Mn28, containing 0.28 wt. % C and 28.00 wt. % Mn, and 65 MPa for X60MnAl17-1, containing 0.56 wt. % C, 17.00 wt. % Mn, and 1.35 wt. % Al.
Based on Figure 4, one could think, that the connection between σ S S and σ S R O is approximately linear. However, it needs to be kept in mind that the calculation of σ S S is based on linear approximations as well. Additionally, it should be pointed out that although the contribution of σ S S to σ Y is always greater than that of σ S R O , the proportion between σ S S and σ S R O does not seem to follow a linear trend based on chemical composition. This is illustrated by the line in Figure 4 that represents a relation of σ S S σ S R O = 1 . To further point out the differences in the contributions to σ Y , the individual contributions of σ GS , σ S S , and σ S R O are portrayed in Figure 5.
According to Figure 5, the contribution of σ S R O to σ Y is significant. The proportional contribution of σ S R O to σ Y is in the range of 7.1% for X30Mn28 up to 19.1% for X70Mn24. Since a grain size effect on σ S R O is not assumed in this work, the impact of σ S R O for a given steel composition with larger grain size is higher than for the same steel composition with finer grains. Correspondingly, the proportional contribution of σ S R O that is portrayed in Figure 5 is higher for compositions with larger grains and thus, with smaller grain size strengthening effects σ G S .
The comparison of σ S S and σ S R O for the steels X30Mn13 and X30Mn22 could lead to the conclusion that an increase in Mn content does not significantly influence σ S R O . On the other hand, by comparing X30Mn22 with X30Mn28, it does appear that Mn actually reduces σ S R O . Accordingly, based on these calculations, no clear tendency of the impact of a single alloying element on σ S R O can be found.
For a more detailed analysis, the impact of C and Mn on σ S R O is portrayed in Figure 6 and Figure 7 with regards to the development of σ S R O , σ S S , E r a n d o m , and E S R O for steels with and without Al.
Figure 6 analyzes a steel with a constant Mn content of 22 wt. % with increasing C content from 0.0 wt. % up to 0.8 wt. %. Additionally, seven steels that were studied in this work were selected to illustrate the impact of changes in Mn content and Al content.
According to Figure 6a, σ S R O and σ S S increase approximately linearly with C content. This can be explained by an increasing number of cells with carbon atoms. In theory, this increase in σ S R O will only occur as long as there are cells with a preferential structure available, which leads to a decrease in ordering energy E S R O , as can be seen in Figure 6b. If no preferential cells were available, C atoms would have to be situated in non-preferential cells with Δ E C k l > 0   eV . To a small degree, this effect is included in the calculation of E r a n d o m because C is assumed to be homogeneously distributed among all cells, including those with positive and, thus, unfavorable values of Δ E C k l . This explains the lower rate of decrease of E r a n d o m in Figure 6b compared to E S R O . However, σ S R O will always remain positive because, according to Equation (2), as long as E r a n d o m is bigger than E S R O , σ S R O will be positive. This premise is always fulfilled by using the assumption in this work, that carbon will be situated in the cells of lowest possible ordering energy first and that the distribution of alloying elements for the calculation of E r a n d o m and E S R O is the same.
Looking at Figure 6a, the steels with Al, X30MnAl22-1 and X30MnAl17-1, show positive deviations for the calculated values of σ S R O and X30MnAl17-1 shows a negative deviation for σ S S . In the case of solid solution strengthening for X30MnAl17-1, this can be explained by the low Mn content of 16.80 wt. % compared to the 22.00 wt. % used as a basis of calculation of σ S S in Figure 6. Accordingly, X70Mn24 and X60Mn23 show slightly higher values for than σ S S it would be the case for 22.00 wt. % Mn.
The increase of σ S R O can be explained by the impact of the introduction of new cell configurations containing Al atoms. According to Table 1, n C 08 is the lowest possible energy configuration. However, since the Al content of the steels in this study is very low, only very few cells with the configuration n C 08 are available. The majority of Al atoms is in cells also containing Mn atoms. According to [33], the introduction of Al atoms as a 2NN into cells with Mn atoms as 1NN does not necessarily lead to a lower ordering energy, but is rather dependent of the number of 1NN Mn atoms. The molar fraction of Mn in a steel containing 22 wt. % Mn obviously depends on the wt. % of the other alloying elements. For the X30MnAl22-1, the molar fraction of Mn is approximately 0.23 . According to Equation (5), the most common cell configuration is expected to be n 30 . The ordering energy of the cell configuration n C 30 is Δ E C 30 = 0.0501   eV . This means that this cell configuration is preferential to a pure iron FCC-unit cell. Increasing the number of Al atoms in this cell configuration to 1 and 2 respectively results in Δ E C 31 = 0.0082   eV and Δ E C 32 = 0.0338   eV . Both cell configurations are not preferential with respect to n C 30 which according to previous assumptions means that the cells n C 31 and n C 32 will not contain a carbon atom as long as there are other more preferential cell configurations like n 30 that do not already contain a carbon atom. In turn, E S R O will not be negatively influenced by cell configurations like n C 31 and n C 32 . However, E r a n d o m will include the negative impact of unfavorable cell configurations which will lead to an increase of E r a n d o m and a total increase of σ S R O .
Overall, it can be said that the positive impact of Al on σ S R O is a result of an introduction of favorable cell configurations containing Al atoms that inhibit a carbon atom and the avoidance of cell configurations of unfavorable cell configurations.
Figure 7 analyzes a steel with a constant C content of 0.3 wt. % with increasing Mn content from 10 wt. % up to 30 wt. %. Seven steels with varying Al and C content were added to the diagram.
Contrary to the impact of C that is documented in Figure 6, a local maximum of σ S R O at approximately 16 wt. % Mn is observed. Looking at Figure 7b, this can be explained by the stagnation of E S R O for Mn contents of 16 wt. % and higher. At the same time, E r a n d o m is decreasing, leading to a smaller difference between E r a n d o m and E S R O and accordingly to a lower σ S R O .
The stagnation of E S R O can be explained by considering the method of calculation for σ S R O . For Al-free steels, the energetically optimal cell is n 60 . As long as θ C   < n C 60 , an increase in the amount of C atoms will decrease E S R O . This is due to the underlying assumption that C atoms will be situated in preferential cells of lowest possible ordering energy. At θ C   =   n C 60 , all C atoms are situated in cells of optimal configuration and thus, E S R O is minimized. Due to the way E r a n d o m is calculated in Equation (7), E r a n d o m continues to decrease at θ C   > n C 60 . As a result, σ S R O decreases.
In accordance with Figure 6a, Figure 7a shows that additions of Al and C increase σ S R O significantly.
The contribution of solid solution hardening σ S S increases linearly with the Mn and C content, as it is to be expected considering Equation (10).
To conclude, unlike for the addition of C, Mn additions do not always increase σ S R O . This is a result of a stagnation of E S R O because of a saturation of preferential cells with carbon. In theory, the same effect should be noticeable for C and Al, as well. However, because of the comparably low alloying contents of those elements, it is very unlikely that a point of saturation will be reached in common alloys.
This decrease in σ S R O can be understood as the impact of the deviation from linear solid solution hardening. If the ratio of Mn atoms to Fe atoms becomes increasingly larger and the amount of cells with randomly-distributed Mn increases, the impact of an actual ordering of Mn decreases since fewer cells are brought into short-range ordering to minimize the ordering energy.
Overall, the change of E O r d e r with varying contents of alloying elements can be understood as the effectiveness of an alloy addition of those elements with respect to the increase in strength by short-range ordering. Since it is assumed that σ SS follows a linear trend and σ GS is fixed, | δ E O r d e r δ χ i | is an indicator for the effectiveness of an addition of an alloying element. A higher | δ E O r d e r δ χ i | means a greater efficiency in employing an addition of the alloying element i. Looking at Figure 6 and Figure 7, C and Al both appear to be more efficient than Mn in increasing the yield strength by short-range ordering.
In order to compare the accuracy of the model used in this work to calculate σ Y with other descriptive models, σ Y , m o d e l according to models by Pickering [15], Lorenz [15], Bouaziz [11], De Cooman [2], and Choi [46] has been calculated. The equations used for σ Y , m o d e l are listed in Table 3. The coefficients used for the alloying contents are in MPa wt .   % , coefficients for grain size influence are in MPam0.5 unless stated otherwise. It needs to be pointed out that although all of the models are for austenitic steels, only the models of Bouaziz, De Cooman, and Choi are specifically for HMnS. Additionally, some of the models include coefficients for the Hall-Petch effect while others neglect them. In order to be able to compare the models, the grain size effect has been calculated according to Equation (4) and added to the model calculations of Choi and Bouaziz.
Figure 8 compares σ Y , m o d e l with σ Y , c a l c and σ Y , e x p .
Due to the difference in the steel concepts those models were developed for, significant deviations from σ Y , e x p can be seen. Those deviations are a result of a focus on other alloying elements. For example, Pickering does not consider the impact of Mn on yield strength. Nevertheless, the models of Choi [46] and De Cooman [2] offer reasonable results, but are not as close to the experimental data as the model developed in this work. Additionally, the model by De Cooman does not take compositional variances into account, while the model of Choi considers the impact of Mn on yield strength to be negative.
The models of Pickering and Lorenz are not applicable to HMnS based on the results of the calculations. The model by Bouaziz tends to overestimate the yield strength.

4. Conclusions

A basic ab initio-based model for the calculation of σ Y has been developed. The model shows good agreement with the experimental values. Discrepancies between experimental σ Y , e x p and calculation results σ Y , c a l c are smaller than 10%. While many complexity-reducing assumptions have been made that lead to inaccuracies when calculating σ S R O , this model portrays that the general concept of combining ab initio simulations with experimental data can lead to feasible results on a macroscopic scale. This is also shown by the calculated positive impact of Mn and Al on yield strength by solid solution strengthening in a reasonable range.
The impact of the alloying content of C, Mn, and Al on the strengthening by short-range ordering was analyzed. For the given steel compositions with fully austenitic microstructures in this work, the values for σ S R O are as high as approximately 65 MPa. This significant strengthening leads to a proportional contribution to σ Y of nearly 20%.
While an increase in C and Al increases σ S R O as long as a sufficient number of atoms of substitutional alloying elements are available to form preferential cell structures, a limiting factor for the effectiveness of strength increases by Mn was found.
The impact of short-range ordered clusters during deformation needs to be evaluated. A higher amount of more stable, energetically favorable clusters might not only lead to an increase in σ S R O and σ Y , it might also impact the beginning of serrations in the flow curve, as well as their amplitude if the serrations are resulting from short-distance carbon atom jumps associated with the cutting of clusters.

Acknowledgments

The authors are thankful for the financial support from Deutsche Forschungsgemeinschaft (DFG) within the Project C2: “Material Properties of High Mn-Steels” of the Collaborative Research Centre (SFB) 761 “Steel-ab initio”. The provision of experimental data by Sebastian Wesselmecking (Steel Institute, RWTH Aachen University) is very much appreciated as well as the provision of ab initio calculation results by Dronskowski, Tobias Timmerscheidt, and Dmitri Bogdanovski (Department of Inorganic Chemistry, RWTH Aachen University) of Project A1: “ab initio Quantum chemistry of the Fe-Mn-C system” of SFB 761 “Steel-ab initio”.

Author Contributions

Simon Sevsek analyzed the data. Simon Sevsek and Wolfgang Bleck wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. De Cooman, B.C.; Kwon, O.; Chin, K.-G. State-of-the-knowledge on TWIP steel. Mater. Sci. Technol. 2012, 28, 513–527. [Google Scholar] [CrossRef]
  2. De Cooman, B.C. High Mn TWIP steel and medium Mn steel. In Automotive Steels: Design, Metallurgy, Processing and Applications, 1st ed.; Rana, R., Singh, S.B., Eds.; Woodhead Publishing: Cambridge, UK, 2017; pp. 317–385. ISBN 9780081006382. [Google Scholar]
  3. Haase, C.; Ingendahl, T.; Güvenc, O.; Bambach, M.; Bleck, W.; Molodov, D.A.; Barrales-Mora, L.A. On the applicability of recovery-annealed Twinning-Induced Plasticity steels: Potential and limitations. Mater. Sci. Eng. A 2016, 649, 74–84. [Google Scholar] [CrossRef]
  4. Abbasi, A.; Dick, A.; Hickel, T.; Neugebauer, J. First-principles investigation of the effect of carbon on the stacking fault energy of Fe–C alloys. Acta Mater. 2011, 59, 3041–3048. [Google Scholar] [CrossRef]
  5. Bouaziz, O.; Zurob, H.; Chehab, B.; Embury, J.D.; Allain, S.; Huang, M. Effect of chemical composition on work hardening of Fe–Mn–C TWIP steels. Mater. Sci. Technol. 2011, 27, 707–709. [Google Scholar] [CrossRef]
  6. Tewary, N.K.; Ghosh, S.K.; Chatterjee, S. Effect of Al Content in Low Carbon High Manganese TWIP Steel. Key Eng. Mater. 2016, 706, 16–22. [Google Scholar] [CrossRef]
  7. Bäumer, A.; Jimenez, J.A.; Bleck, W. Effect of temperature and strain rate on strain hardening and deformation mechanisms of high manganese austenitic steels. Int. J. Mater. Res. 2010, 101, 705–714. [Google Scholar] [CrossRef]
  8. Choi, K.; Seo, C.-H.; Lee, H.; Kim, S.K.; Kwak, J.H.; Chin, K.G.; Park, K.-T.; Kim, N.J. Effect of aging on the microstructure and deformation behavior of austenite base lightweight Fe–28Mn–9Al–0.8C steel. Scr. Mater. 2010, 63, 1028–1031. [Google Scholar] [CrossRef]
  9. Gutierrez-Urrutia, I.; Raabe, D. Microbanding mechanism in an Fe–Mn–C high-Mn twinning-induced plasticity steel. Scr. Mater. 2013, 69, 53–56. [Google Scholar] [CrossRef]
  10. Liu, F.; Dan, W.J.; Zhang, W.G. Strain hardening model of TWIP steels with manganese content. Mater. Sci. Eng. A 2016, 674, 178–185. [Google Scholar] [CrossRef]
  11. Bouaziz, O.; Allain, S.; Scott, C.P.; Cugy, P.; Barbier, D. High manganese austenitic twinning induced plasticity steels: A review of the microstructure properties relationships. Curr. Opin. Solid State Mater. Sci. 2011, 15, 141–168. [Google Scholar] [CrossRef]
  12. Bouaziz, O. Strain-hardening of twinning-induced plasticity steels. Scr. Mater. 2012, 66, 982–985. [Google Scholar] [CrossRef]
  13. Ohkubo, N.; Miyakusu, K.; Uematsu, Y.; Kimura, H. Effect of Alloying Elements on the Mechanical Properties of the Stable Austenitic Stainless Steel. ISIJ Int. 1994, 34, 764–772. [Google Scholar] [CrossRef]
  14. Kang, J.-H.; Ingendahl, T.; von Appen, J.; Dronskowski, R.; Bleck, W. Impact of short-range ordering on yield strength of high manganese austenitic steels. Mater. Sci. Eng. A 2014, 614, 122–128. [Google Scholar] [CrossRef]
  15. Sieurin, H.; Zander, J.; Sandström, R. Modelling solid solution hardening in stainless steels. Mater. Sci. Eng. A 2006, 415, 66–71. [Google Scholar] [CrossRef]
  16. Oda, K.; Fujimura, H.; Ino, H. Local interactions in carbon-carbon and carbon-M (M: Al, Mn, Ni) atomic paris in FCC gamma-iron. J. Phys. Condens. Matter 1994, 6, 679–692. [Google Scholar] [CrossRef]
  17. Bentley, A.P. Ordering in Fe–Mn–Al–C austenite. J. Mater. Sci. Lett. 1986, 5, 907–908. [Google Scholar] [CrossRef]
  18. McLellan, R.B.; Dunn, W.W. A quasi-chemical treatment of interstitial solid solutions: It application to carbon austenite. J. Phys. Chem. Solids 1969, 30, 2631–2637. [Google Scholar] [CrossRef]
  19. Pelton, A.D.; Kang, Y.-B. Modeling short-range ordering in solutions. Int. J. Mater. Res. 2007, 98, 907–917. [Google Scholar] [CrossRef]
  20. Honeycombe, R.W.K. The effect of temperature and alloying additions on the deformation of metal crystals. Progress Mater. Sci. 1961, 9, 95–120. [Google Scholar] [CrossRef]
  21. Oda, K.; Kondo, N.; Shibata, K. X-ray Absorption Fine Structure Analysis of Interstitial (C, N)-Substitutional (Cr) Complexes in Austenitic Stainless Steels. ISIJ Int. 1990, 30, 625–631. [Google Scholar] [CrossRef]
  22. Von Appen, J.; Dronskowski, R. Carbon-Induced Ordering in Manganese-Rich Austenite—A Density-Functional Total-Energy and Chemical-Bonding Study. Steel Res. Int. 2011, 82, 101–107. [Google Scholar] [CrossRef]
  23. Song, W.; Zhang, W.; von Appen, J.; Dronskowski, R.; Bleck, W. k-Phase Formation in Fe–Mn–Al–C Austenitic Steels. Steel Res. Int. 2015, 86, 1161–1169. [Google Scholar] [CrossRef]
  24. Raabe, D.; Roters, F.; Neugebauer, J.; Gutierrez-Urrutia, I.; Hickel, T.; Bleck, W.; Schneider, J.M.; Wittig, J.E.; Mayer, J. Ab initio-guided design of twinning-induced plasticity steels. MRS Bull. 2016, 41, 320–325. [Google Scholar] [CrossRef]
  25. Medvedeva, N.I.; Park, M.S.; van Aken, D.C.; Medvedeva, J.E. First-principles study of Mn, Al and C distribution and their effect on stacking fault energies in fcc Fe. J. Alloys Compd. 2014, 582, 475–482. [Google Scholar] [CrossRef]
  26. Timmerscheidt, T.; Dey, P.; Bogdanovski, D.; von Appen, J.; Hickel, T.; Neugebauer, J.; Dronskowski, R. The Role of κ-Carbides as Hydrogen Traps in High-Mn Steels. Metals 2017, 7, 264. [Google Scholar] [CrossRef]
  27. McLellan, R.B. Cell models for interstitial solid solutions. Acta Metall. 1982, 30, 317–322. [Google Scholar] [CrossRef]
  28. Bhadeshia, H.K.D.H. Diffusion of carbon in austenite. Met. Sci. 1981, 15, 477–479. [Google Scholar] [CrossRef]
  29. McLellan, R.B.; Ko, C. The C–C interaction energy in iron-carbon solid solutions. Acta Metall. 1987, 35, 2151–2156. [Google Scholar] [CrossRef]
  30. McLellan, R.B.; Ko, C. The diffusion of carbon in austenite. Acta Metall. 1988, 36, 531–537. [Google Scholar] [CrossRef]
  31. Gholizadeh, H.; Draxl, C.; Puschnig, P. The influence of interstitial carbon on the g-surface in austenite. Acta Metall. 2013, 61, 341–349. [Google Scholar] [CrossRef]
  32. Lee, S.-J.; Kim, J.; Kane, S.N.; DeCooman, B.C. On the origin of dynamic strain aging in twinning-induced plasticity steels. Acta Metall. 2011, 59, 6809–6819. [Google Scholar] [CrossRef]
  33. Timmerscheidt, T.; Dronskowski, R. An Ab Initio Study of Carbon-Induced Ordering in Austenitic Fe–Mn–Al–C Alloys. Steel Res. Int. 2017, 88, 1–10. [Google Scholar] [CrossRef]
  34. Sumin, V.V.; Chimid, G.; Rashev, T.; Saryivanov, L. The Neutron-Spectroscopy Proof of the Strong Cr–N Interactions in Nitrogen Stainless Steels. Mater. Sci. Forum 1999, 318–320, 31–40. [Google Scholar] [CrossRef]
  35. Bouaziz, O.; Barbier, D. Benefits of Recovery and partial Recrystallization of Nano-Twinned Austenitic Steels. Adv. Eng. Mater. 2013, 15, 976–979. [Google Scholar] [CrossRef]
  36. Fiore, N.F.; Bauer, C.L. Binding of solute atoms to dislocations. Progress Mater. Sci. 1968, 13, 85–134. [Google Scholar] [CrossRef]
  37. Butt, M.Z.; Feltham, P. Solid-solution hardening. J. Mater. Sci. 1993, 28, 2557–2576. [Google Scholar] [CrossRef]
  38. Saeed-Akbari, A.; Imlau, J.; Prahl, U.; Bleck, W. Derivation and Variation in Composition-Dependent Stacking Fault Energy Maps Based on Subregular Solution Model in High-Manganese Steels. Metall. Mater. Trans. A 2009, 40, 3076–3090. [Google Scholar] [CrossRef]
  39. Kang, J.-H.; Duan, S.; Kim, S.-J.; Bleck, W. Grain Boundary Strengthening in High Mn Austenitic Steels. Metall. Mater. Trans. A 2016, 47, 1918–1921. [Google Scholar] [CrossRef]
  40. Kang, S.; Jung, Y.-S.; Jun, J.-H.; Lee, Y.-K. Effects of recrystallization annealing temperature on carbide precipitation, microstructure and mechanical properties in Fe–18Mn–0.6C–1.5Al TWIP steel. Mater. Sci. Eng. A 2010, 527, 745–751. [Google Scholar] [CrossRef]
  41. De las Cuevas, F.; Reis, M.; Ferraiuolo, A.; Pratolongo, G.; Karjalainen, L.P.; Alkorta, J.; Gil Sevillano, J. Hall-Petch relationship of a TWIP steel. Key Eng. Mater. 2010, 423, 147–152. [Google Scholar] [CrossRef]
  42. Bartlett, L.N.; van Aken, D.C.; Medvedeva, J.; Isheim, D.; Medvedeva, N.I.; Song, K. An Atom Probe Study of Kappa Carbide Precipitation and the Effect of Silicon Addition. Metall. Mater. Trans. A 2014, 45, 2421–2435. [Google Scholar] [CrossRef]
  43. Yao, M.J.; Welsch, E.; Ponge, D.; Haghighat, S.; Sandlöbes, S.; Choi, P.; Herbig, M.; Bleskov, I.; Hickel, T.; Lipinska-Chwalek, M.; et al. Strengthening and strain hardening mechanisms in a precipitation-hardened high-Mn lightweight steel. Acta Mater. 2017, 140, 258–273. [Google Scholar] [CrossRef]
  44. Gottstein, G. Physikalische Grundlagen der Materialkunde, 3rd ed.; Springer: Berlin/Heidelberg, Germany, 2007; ISBN 978-3-540-71104-9. [Google Scholar]
  45. Li, M.; Chang, H.; Kao, P.; Gan, D. The effect of Mn and Al contents on the solvus of κ phase in austenitic Fe–Mn–Al–C alloys. Mater. Chem. Phys. 1999, 59, 96–99. [Google Scholar] [CrossRef]
  46. Choi, W.S.; de Cooman, B.C.; Sandlöbes, S.; Raabe, D. Size and orientation effects in partial dislocation-mediated deformation of twinning-induced plasticity steel micro-pillars. Acta Mater. 2015, 98, 391–404. [Google Scholar] [CrossRef]
Figure 1. FCC-unit cell with 14 lattice positions for matrix atoms and one central interstitial atom.
Figure 1. FCC-unit cell with 14 lattice positions for matrix atoms and one central interstitial atom.
Metals 08 00034 g001
Figure 2. Energetically most favorable cell configurations.
Figure 2. Energetically most favorable cell configurations.
Metals 08 00034 g002aMetals 08 00034 g002b
Figure 3. Comparison of calculated and experimental values for σ Y .
Figure 3. Comparison of calculated and experimental values for σ Y .
Metals 08 00034 g003
Figure 4. Contributions of σ S S and σ S R O to yield strength.
Figure 4. Contributions of σ S S and σ S R O to yield strength.
Metals 08 00034 g004
Figure 5. Proportional contribution of σ S S , σ S R O , and σ G S to σ Y .
Figure 5. Proportional contribution of σ S S , σ S R O , and σ G S to σ Y .
Metals 08 00034 g005
Figure 6. Impact of C on (a) σ S R O , σ S S ; and (b) E r a n d o m and E S R O .
Figure 6. Impact of C on (a) σ S R O , σ S S ; and (b) E r a n d o m and E S R O .
Metals 08 00034 g006
Figure 7. Impact of Mn on (a) σ S R O , σ S S ; and (b) E r a n d o m and E S R O .
Figure 7. Impact of Mn on (a) σ S R O , σ S S ; and (b) E r a n d o m and E S R O .
Metals 08 00034 g007
Figure 8. Comparison of models of σ Y for the steels used in this paper.
Figure 8. Comparison of models of σ Y for the steels used in this paper.
Metals 08 00034 g008
Table 1. Cell configurations of lowest ordering energy Δ E C k l .
Table 1. Cell configurations of lowest ordering energy Δ E C k l .
Cell Name n C k l k Mn Atomsl Al Atoms Δ E C k l   [ eV ]
n C 08 08−1.2867
n C 18 18−1.1524
n C 28 28−1.0559
n C 07 07−0.9304
n C 68 68−0.9287
n C 38 38−0.9233
Table 2. Chemical composition of HMnS used in this study in wt. %; grain size dm in μm; experimentally determined and calculated yield strength σY,exp and σY,calc in MPa, experimental data of alloys marked with * are taken from [14].
Table 2. Chemical composition of HMnS used in this study in wt. %; grain size dm in μm; experimentally determined and calculated yield strength σY,exp and σY,calc in MPa, experimental data of alloys marked with * are taken from [14].
AlloyCMnAl d m σ Y , e x p σ Y , calc
X60MnAl17-10.5617.001.358328356
X50MnAl15-1 *0.4915.401.3027294281
X30MnAl22-10.3322.461.2118277275
X30MnAl17-10.3216.801.4717255276
X70Mn240.7123.500.0116334334
X60Mn230.5723.210.0110328334
X50Mn30 *0.5030.000.0030276273
X50Mn220.5421.950.0014328311
X50Mn18 *0.5317.900.018356341
X30Mn280.2828.000.009291294
X30Mn230.3222.790.0110284291
X30Mn22 *0.3122.280.003381384
X30Mn13 *0.3012.740.0110266274
Table 3. Solid solution hardening models in the literature.
Table 3. Solid solution hardening models in the literature.
Model σ Y , model
Pickering68 + 493N + 354C + 3.7Cr + 14Mo + 20Si + 0.22d−1/2
Lorenz40 + 450N + 525C + 2Mn + 8.4Cr + 22Mo + 5Cu + 0.77d−1/2
Bouaziz228 + 137C − 2Mn + σGS
De Cooman189 + 413/√(d[μm])
Choi97 + 279C − 1.5Mn + 49.6Si + 20.5Al + σGS

Share and Cite

MDPI and ACS Style

Sevsek, S.; Bleck, W. Ab Initio-Based Modelling of the Yield Strength in High-Manganese Steels. Metals 2018, 8, 34. https://doi.org/10.3390/met8010034

AMA Style

Sevsek S, Bleck W. Ab Initio-Based Modelling of the Yield Strength in High-Manganese Steels. Metals. 2018; 8(1):34. https://doi.org/10.3390/met8010034

Chicago/Turabian Style

Sevsek, Simon, and Wolfgang Bleck. 2018. "Ab Initio-Based Modelling of the Yield Strength in High-Manganese Steels" Metals 8, no. 1: 34. https://doi.org/10.3390/met8010034

APA Style

Sevsek, S., & Bleck, W. (2018). Ab Initio-Based Modelling of the Yield Strength in High-Manganese Steels. Metals, 8(1), 34. https://doi.org/10.3390/met8010034

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop