2. Mechanics of Plastic Flow during Deformation
The mechanical properties of dual phase ferritic-martensitic steel depend on the chemical composition, the volume fraction of martensite, the volume fraction of ferrite, the carbon content in martensite, the grain size of martensite, and their strength [
17,
18]. To describe the behavior of dual phase ferritic-martensitic steel under plastic deformation, various constitutive equations were proposed [
9,
15,
19,
20,
21,
22]. Increasing the intercritical temperature increases the amount of austenite generated and this is transformed to martensite during rapid cooling. Thus, the strength and hardness of the material increases as well. The carbon content in martensite is larger for dual phase steels with low volume fractions of martensite. Otherwise, the carbon content in martensite decreases when the volume fraction of martensite increases. The carbon content in martensite controls the phase hardness and influences the final properties of the material. By controlling the metallurgical processes, it is possible to reach ferritic-martensitic structures with volume fractions of martensite from 35 to 50% with a wide combination of strength and plastic properties [
17,
18].
The effective use of dual phase steels in the automotive industry requires a better understanding of how they behave in crashes, as well as how they behave when processed by stamping to the structural parts of the safety zones. Nowadays, numerical simulations of crash tests and metal forming processes, based on the Finite Element Method, are widely used to predict the deformation behavior of materials. Thus, to describe the material behavior under deformation, the following constitutive equations are used [
23]:
where
σs is the true stress, K is the material constant,
n is the strain-hardening exponent that expresses the intensity of the strain-hardening and the ability of the material to deform uniformly,
ϕi is true strain
and
ϕ0 is the pre-strain [
24].
These models can be used to prepare the production of DP steel grades with precisely defined “tailor-made” properties for the components of the vehicle’s deformation zones at the front and the side impact. When selecting material for the car-body safety zones, the main criterion is resistance to deformation (that is, deformation work) that is consumed at the crash. This can be determined by the tensile test record
σs–
ϕ (
Figure 1):
where
V0 is the specimen volume on the initial length
L0,
ϕUE is the uniform true strain (true strain at tensile strength), and
ϕ0.002 = 0.002 is the true strain at yield strength.
The parts of the body deformation zones are elastically and plastically deformed during impact and during their production. However, crash tests are only concerned with plastic deformation. After we insert Equation (1) into Equation (3) and make adjustments, we get the following:
After the integration and adjustment Equation (4) we get the following:
The values of the material constant K and the strain-hardening exponent n can be determined from the tensile test record by regression analysis. However, to gain a better understanding of the mechanics of the deformation process, the strain-hardening exponent
n, and the material constant K can be determined from the mutual bonds of the mechanical properties of metallic materials. If Equation (1) (or, by analogy, Equation (3)) is subjected to a logarithmic operation, we get the following linear dependence:
and we express the contribution to strain-hardening at tensile strength R
m (R
m refers to the ultimate tensile strength) with respect to the yield strength, depending on the uniform deformation, in the interval from
ϕ = 0.002 to
ϕUE [
25] as follows:
which yields the exponent of the strain-hardening
n:
The
n value is not constant throughout the uniform deformation, so it is necessary to expect a certain uncertainty in the calculation of the deformation work and the actual strength, especially in the case of minor strains. The exact determination of the strain-hardening exponent requires the division of an even deformation region into several intervals and the expression of the strain-hardening exponent in terms of deformation:
where
n0 is the strain-hardening exponent found in the first interval (for example,
ϕi is the true strain between 0.002 and 0.02), p is the constant determined by the approximation of the dependence of the strain-hardening exponent on the deformation at individual intervals.
From the Equation (6), it follows that the material constant K will be
Upon adjustment, we get the following:
The above-mentioned mechanical properties of materials (the yield strength Re, the tensile strength Rm, the material constant K, the strain-hardening exponent n, the maximum value of uniform deformation ϕUE, and so forth) are given by their internal structure—the structure of the material, which, in turn, depends on the chemical composition of steel and on its production technology. The production of “tailored” or “customized” steel grades, with exactly defined properties, requires knowledge of not only the above-mentioned relationships that determine the mechanical properties but also knowledge of the relationships between the structural parameters and the mechanical properties of the metallic materials.
Due to the fact that the structural parts of the safety zones are deformed at higher strain rates when a car crashes, the influence of the strain rate needs to be included in constitutive equations. In Reference [
26] the authors included the influence of strain rate and temperature into these equations. Authors from References [
27,
28] included the strain rate influence into the constitutive equations when predicting the deformation work. It has been found that the influence of the strain rate was low at quasistatic strain rates [
28], but a notable effect was found at higher strain rates and that it is connected to the evolution of the dislocation density.
In the literature, the structural nature of the material properties of ferritic-martensitic steel grades is given a great deal of attention [
29,
30,
31]. Based on the dislocation theory, founded on the motion of dislocations and their interaction with various obstacles (grain boundaries, precipitates, interstitial atoms, fractions of different phases, as well as other dislocations), the actual stress necessary for plastic deformation flow can be expressed in terms of the individual contributors to hardening:
where
σ0 is Peierls stress necessary to overcome the lattice friction stress, the resistance of alloying elements dissolved in solid solution, the precipitation matrix resistance, and the lattice defects [
32]; ∆
σg is the hardening effect depending on the size of the ferritic grain; ∆
σS is the effect of substitute hardening; ∆
σIN is the effect of interstitic hardening; ∆
σP is the effect of precipitation hardening; ∆
σPR is the effect of perlite hardening; ∆
σSG is the effect subgrain hardening (also possible to be expressed as ∆
σFMaB—the hardening through bainitic or martensitic fractions or plates); ∆
σD is the dislocation density hardening effect, and so forth [
33,
34].
For dual phase ferritic-martensitic (DP) steel grades, Equation (1) can be adjusted as follows:
The Peierls stress
σ0 [
32,
33]
The hardening effect of ferritic grain size:
The hardening effect of the martensitic or bainitic fractions:
The hardening effect of the dislocation density:
where
dα is the mean grain size of ferrite,
ky is the strengthening coefficient, α is a material constant, G is the shear modulus (80,000 MPa), b is Burger’s vector, and
ρD is the dislocation density.
The aims of the experimental research were to prepare materials with different volume fractions of martensite up to 50% from commercial steels, to describe the relationship between the mechanical properties and the temperature of the intercritical annealing, and to describe the relationship between properties that are sensitive to changes of the sub-structural parameter when cold deformed.
3. Materials and Methods
The deformation behavior of the dual phase steel types (Equations (13)–(17)) depends mainly on the chemical composition, the volume of martensite, the morphology and distribution in the ferrite matrix, as well as the ferrite grain size d
α. The aim of the experimental research was to prepare materials (states) with a martensite volume of up to 50% from commercially produced low carbon steel types of 3–3.3 mm thickness: A (X60), B (S460MC), and C (DP600), whose chemical composition and values of carbon equivalent C
E calculated from Equation (18) [
35] are listed in
Table 1.
The microstructures of the initial materials A, B, and C used are shown in
Figure 2. The low carbon steel microstructure (A) is ferritic-pearlite (
Figure 2a). The low-carbon micro-alloyed steel (B) microstructure is ferritic-pearlite with a low perlite content (
Figure 2b). The microstructure of steel C is a ferritic-martensitic one with a martensitic volume of 24% (
Figure 2c). As can be seen from
Table 1, the carbon content and the average size of the ferritic grain
dα are approximately equal for as-received steels A, B, and C. The dispersion of the mean ferrite grain size under the surface and in the middle of the sheet thickness of the as-received A, B, and C materials was ±10%.
Prior to the heat treatment, the proper starting and final temperatures of the transformation of ferrite to austenite A
C1, A
C3, A
r1, and A
r3 were set according to Andrews [
36] (
Table 2). The non-recrystallization temperature Tnr and the critical cooling time between 800 and 500 °C for the beginning of the perlite precipitation were calculated according to the equations listed in Reference [
37].
The as-received materials samples were prepared by single-step annealing in a flowing cantalum furnace REH-B-10-60 (Linn High Therm GmbH, Bad Frankenhausen, Germany) with a protective argon atmosphere. The samples made out of material A (marked as DPA) were annealed at temperatures 740, 790, and 840 °C. Then, considering the results reached, the samples made out of materials B and C (marked as DPB and DPC) were annealed at temperatures of 750 and 820 °C (which lie between the temperatures A
C1–A
C3) with the same steady-state 10 min for each temperature, followed by cooling in water with a cooling rate of 30 °C/s [
38].
Samples for metallographic analysis were hot mounted in dentacrylate, wet grinded (sandpaper 220–1200), and polished by diamond grit in suspension. Then, the samples were etched in 2% Nital.
The grain size
dα was identified by the linear method according to the Slovak standard STN 42 0462 on the microscope Olympus GX71. The volume fraction of secondary phases (
VFSP) was measured by the grid method (square foil 15 × 15 cm with grid 1 × 1 cm) and by the image analysis method using the image analyzer Image J at a magnification 1000× [
38].
Mechanical properties of as-received materials A, B, and C, and the samples after annealing DPA, DPB, and DPC were measured by static tensile tests according to STN EN ISO 6892-1 at room temperature on a testing machine TIRAtest 2300. These are shown in
Table 3. The transversal feed was 1 mm·min
−1 and the corresponding quasistatic strain rate was 0.003 s
−1. Five specimens for each material and annealing state were tested. The specimen’s shape is shown in
Figure 3 and the dimensions were as follows:
L0 = 20 and 35 mm,
LC = 50 mm,
a = 14 mm,
h = 15 mm,
b = 8 ÷ 10 mm [
38].
4. Results and Discussion
The range of annealing conditions within the temperatures of A
C1–A
C3 applied to commercially available low carbon steel X60 (A), low-alloyed steel S460MC (B), and dual phase steel DP 600 (C) were allowed to reach a wide range of microstructure states, with a martensite volume between 20.4% and 68.2% and a ferritic grain size between 4.7 and 7.7 µm, as seen in
Figure 4,
Figure 5 and
Figure 6, respectively. Hereinafter, these phases are designated as DPA
740, DPA
790, DPA
820, DPB
750, DPB
820, DPC
750, and DPC
820.
In determining the volume fraction of ferrite and martensite, the ferrite fraction was evaluated as the dominant phase, while the sum of all other phases (martensite, residual austenite, and bainite) represents the fraction of the secondary phase particles (FSP). This means that the fraction of the purely martensitic phase is slightly overestimated. However, the fraction of bainite, cementite, and residual austenite in the analyzed states of DPA, DPB, and DPC was ±3% within the distribution of the volume fraction of martensite. In the samples of the DPA
740, DPA
790, and DPA
840 states obtained by the heat treatment from the initial material A, the volume fraction of the martensite ranged between 23.4% and 68.2% and the ferritic grain size ranged from 4.7 to 7.7 µm; in the samples of the DPB
750 and DPB
820 states obtained by the heat treatment from the initial material B, the volume fraction of the martensite ranged between 22.4% and 58.6%, the ferritic grain size ranged from 7.7 to 3.1 µm; and in the samples of the DPC
750, DPC
820 states obtained by the heat treatment from the initial material C, the volume fraction of the martensite ranged between 23.8% and 64.4% and the ferritic grain size ranged from 5.2 to 7.3 µm. Thus, the assumption that the states obtained from the A and C materials with higher values of carbon equivalent would result in greater percentages of the secondary phase fractions—shown in
Table 2 and
Figure 7a—depending on the annealing temperature in the range between 740 and 820 °C, has been confirmed.
The dependence of the volume fractions of the secondary phase (
Figure 7a) can be described by the regression model:
Independent of the observed structural differences, increased annealing temperature T
IA resulted in refinement of the mean size of the ferritic grain in examined states—
Figure 7b. This tendency of grain size refinement in the annealing temperature interval between 740 and 820 °C has been described by the regression model:
The interaction between the fractions of the secondary phase and the size of the ferrite grains as it is shown in
Figure 8 is described by the regression model:
According to References [
38,
39], the increase in the secondary phase fractions of the dual phase ferritic-martensitic steel grades is mostly due to the number of grains of the secondary phase fractions rather than due to their volume in the structure. Regardless of the annealing temperature, it is possible to further increase the martensite fraction in the volume and the refinement of ferrite grain by increasing the rate of cooling [
40].
The decisive criterion for the choice of steel sheets for the structural parts of the deformation zone of the car-body is the deformation work, which expresses the absorption capacity during a crash. The deformation work can be determined with greater uncertainty from the conventional stress-strain diagram as seen in
Figure 1.
or with less uncertainty (more precisely) from the true stress–true strain diagram by application of Equation (3).
For the obtained states, the attention was focused on the analysis of the relationships between the mechanical properties (the yield strength Re, the tensile strength Rm, the uniform elongation
Ag, the total elongation A, the material constant K, and the strain-hardening exponent n) and the parameters of the dual phase ferritic-martensitic steels structure. It follows from the measured results (
Table 2 and
Figure 9) that the increase in the volume fraction of the secondary phase resulted in the increased yield strength values by about 80 MPa and the increased tensile strength by about 78 MPa in the DPA states. For the DPB states, the yield strength increased by 92 MPa and the tensile strength increased by 81 MPa; for the DPC, the yield strength increased by 85 MPa and the tensile strength increased by 96 MPa. Higher strength properties (yield strength, tensile strength) of the obtained states are mainly related to the volume fraction of the secondary phase. The results obtained indicate a linear dependence of the tensile strength on the volume fraction of the secondary phase, described by the following regression equation:
However, the dependency of the tensile strength on the volume fraction of the secondary phase was found to be a linear, ranging between 20% and 50%, described by the regression equation as follows:
No increase in tensile strength has been observed for the volume fraction of the secondary phase greater than 50%. Rather, a decrease in tensile strength values has been noted. We assume that on the one hand, martensite contributes to an increase in tensile strength due to the increased volume of the harder phase (martensite), on the other hand, the carbon content of martensite decreases with the increasing volume of martensite. As is well known, the strength of martensite is mainly determined by its carbon content. Another reason for the reduced tensile strength values may be due to the size of martensite islands and martensite distributions in the ferrite matrix [
11,
32,
39]. At lower annealing temperatures
TIA 740 or 750 °C, martensite was dispersed along the borders of ferritic grains, however, at the annealing temperature of 820 °C, the size of martensitic fractions was greater in comparison to the states obtained under the annealing temperature of 740 or 750 °C (
Figure 4,
Figure 5 and
Figure 6).
It follows, from
Figure 9b, that the lower the size of the ferritic grain d
α−0.5, the larger the yield strength and tensile strength values that have been observed. These curves have been described by the regression models in Equations (27) and (28).
and similarly, in the 20% to 50% interval for the tensile strength
However, it should be noted that the tensile strength dependency on the size of the ferritic grain dα–0.5 is more of a tendency because the residual dispersion value R2 = 0.192 was low.
The influence of strengthening contributors in terms of the individual structural parameters on the strength properties of dual phase ferritic-martensitic materials cannot be assessed separately. For this reason, attention was focused on expressing the summary influence of the structural parameters on the yield strength according to Equation (13) by the unit sum
Sj of the individual parameters of the structure. The unit sum
Sj of the parameters of the structure was determined as the ratio of yield strengths expressed by the Equations (25) and (27) with respect to the yield strength value of the reference material DPA
750 (R
e = 406 MPa) as follows:
where
i is the number of parameters of the structure (
i = 3), R
e,ref is the yield strength value of the reference material, and k
yαi is a constant expressing the influence of the ferritic grain size.
The dependence of the yield strength on the unit sum of the parameters of the structure S
iRe is given in
Figure 10 and expressed through the following regression model:
Then, after inserting Equation (29) into Equation (30), while taking into account the relationships given by Equations (25) and (27) and after the subsequent adjustment, we arrive at the following result:
If we insert Equation (22) into Equation (31) as d
α and we insert Equation (21) as V
FSP, we obtain the relation from which we can predict the yield strength dependent on the annealing temperature:
To determine the deformation work, it is necessary to know the value of the uniform elongation
Ag.
Figure 11 shows that the increasing volume fraction of the secondary phase causes the values of the total elongation and the uniform elongation to drop. A similar trend was found in Reference [
32]. The uniform deformation values
Ag ranged from 8.8 to 16.3 and total deformation A ranged from 17.1 to 25.3. At the annealing temperature of 840 °C, the lower
Ag and A values were recorded compared to the states obtained at 750 °C. We assume that this tendency may be related in particular to the morphology and the distribution of secondary phase fractions.
It should be noted that in most metallic materials, the dependence force on the elongation, or the conventional strength on the deformation, is flat in the area of the maximum uniform elongation
Ag. If we determine the elongation ∆
Li from the conventional diagram (
Figure 1) at the moment when
The
Ag value will be determined by the following relation:
Then, the value of the uniform deformation does not allow for the precise deformation work to be determined by Equation (3). The deformation work in the interval from
A = 0.2% to the maximum uniform deformation
Ag, does not express the overall deformation work of the material, as the material resistance to the deformation increases even with a greater deformation than
Ag (
Figure 1). For this reason, we recommend using a reduced elongation value
Ags, to be determined from the tolerance range as 1/3 of the difference of the total elongation A and the uniform elongation
Ag relative to the standard quadratic deviation of the measured values of
A and
Ag (STDEVA,
Ag), according to the six-sigma method:
Then the reduced value of the uniform elongation will be:
and the true strain (or real deformation) will be:
Trend analyses of the dependence of the immediate stress value on deformation in the interval from
A = 0.002 to the maximum uniform deformation of
Ags (
Figure 1) allow the designers in the automotive industry to understand the differences of mechanical behavior in conventional and advanced high-strength steel grades in crashes and it allows designers to optimize the choice of the materials for individual “tailored” parts of deformation zones. When compared to the initial material A, the (gradient) curve directions (
Figure 12) under the DPA, DPB, and DPC states show that the obtained states of DPA, DPB and DPC exhibit a more favorable course of material resistance to deformation and thus, the course of the deceleration at a crash in comparison to the curve of the steel A.
The direction or the slope of the curves express the degree of steel deformability and also the intensity of the deformation resistance upon deformation.
Figure 12 shows that the states obtained at lower annealing temperatures, whose strain-hardening exponent values are higher than in the states obtained at higher temperatures, exhibit the greatest deformation resistance.
Figure 13 and
Figure 14 show that the material constant K and the strain-hardening exponent n depend on the thermomechanical history of steel, with the strain-hardening exponent n being more sensitive to changes in the parameters of the structure than the material constant K. In terms of physics, the strain-hardening exponent n determines the ability of the steel to distribute the stress along the tensile specimen. For low-carbon steel grades, used in the production of complex car-body shapes, the required value of
n is >0.22. The higher the
n-value, the more uniform the deformation distribution, the greater the steel’s resistance to deformation, and the better its formability [
39,
40,
41].
Equation (13) allows for the estimation of the real strength of the material using empirical models based only on the structural parameters of ferritic-martensitic steel grades. This model describes the deformation behavior dependent on the volume fraction of the secondary phase, the ferrite grain size, and the lattice friction stress required for the dislocation motion. Another important parameter that affects the deformation behavior of metallic materials is the density of dislocations. The density of the dislocations is different for each material. The combined effects of the structural parameters of model (13) differ in the density of dislocations. We can express the density of dislocations from the actual difference in tensile strength:
and the actual strength at the yield strength:
After adjustment, we obtain the contribution coming from the density of dislocations upon deformation at the tensile strength limit:
The high values of the dislocation densities in the dual phase steel grades listed in
Table 4 are not surprising since the transformation of austenite into martensite is the cause of great stress in ferrites. Near the martensite fractions, the dislocation density may be even higher (
Figure 15). For example, Reference [
34] states that, depending on the deformation in DP 500, the dislocation density values can range from 1.5 × 10
14 m
−2 up to 1.7 × 10
15 m
−2.