Next Article in Journal
Numerical Investigation on the Effect of Residual Stresses on the Effective Mechanical Properties of 3D-Printed TPMS Lattices
Previous Article in Journal
Recent Advances and Applications of Machine Learning in Metal Forming Processes
Previous Article in Special Issue
In Silico Contact Pressure of Metal-on-Metal Total Hip Implant with Different Materials Subjected to Gait Loading
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling of Tensile Test Results for Low Alloy Steels by Linear Regression and Genetic Programming Taking into Account the Non-Metallic Inclusions

1
ŠTORE STEEL, d.o.o., Research and Development, 3220 Štore, Slovenia
2
Laboratory for Fluid Dynamics and Thermodynamics, Faculty of Mechanical Engineering, University of Ljubljana, 1000 Ljubljana, Slovenia
3
College of Industrial Engineering, 3000 Celje, Slovenia
4
Faculty of Mechanical Engineering, University of Maribor, 2000 Maribor, Slovenia
*
Author to whom correspondence should be addressed.
Metals 2022, 12(8), 1343; https://doi.org/10.3390/met12081343
Submission received: 30 June 2022 / Revised: 6 August 2022 / Accepted: 9 August 2022 / Published: 12 August 2022

Abstract

:
Štore Steel Ltd. is one of the biggest flat spring steel producers in Europe. The main motive for this study was to study the influences of non-metallic inclusions on mechanical properties obtained by tensile testing. From January 2016 to December 2021, all available tensile strength data (472 cases–472 test pieces) of 17 low alloy steel grades, which were ordered and used by the final user in rolled condition, were gathered. Based on the geometry of rolled bars, selected chemical composition, and average size of worst fields non-metallic inclusions (sulfur, silicate, aluminium and globular oxides), determined based on ASTM E45, several models for tensile strength, yield strength, percentage elongation, and percentage reduction area were obtained using linear regression and genetic programming. Based on modeling results in the period from January 2022 to April 2022, five successively cast batches of 30MnVS6 were produced with a statistically significant reduction of content of silicon (t-test, p < 0.05). The content of silicate type of inclusions, yield, and tensile strength also changed statistically significantly (t-test, p < 0.05). The average yield and tensile strength increased from 458.5 MPa to 525.4 MPa and from 672.7 MPa to 754.0 MPa, respectively. It is necessary to emphasize that there were no statistically significant changes in other monitored parameters.

1. Introduction

The tensile testing of metallic materials, as one of the most important methods of testing, which is usually performed at room temperature (i.e., between 10 °C and 35 °C), includes the following steps:
-
taking of samples,
-
preparation of test pieces (usually machining),
-
pulling the test piece by tensile force until the breakage,
-
obtaining mechanical properties based on measurements and calculations (e.g., tensile strength, yield strength, percentage elongation, percentage reduction area).
-
Mechanical properties obtained by tensile testing mostly depend on:
-
material chemical composition,
-
material macrostructure (e.g., segregations, sheets and strips) and microstructure (e.g., grain size and morphology, content of microstructures),
-
sampling (location and orientation of test pieces),
-
preparation of test pieces (machining, heat treatment),
-
tensile test performance (e.g., equipment, testing speed).
On the other hand, non-metallic inclusions should also not be neglected while their presence as foreign substances can disrupt the existing structure as well as cause cracks and fatigue failure.
Below, a few studies are presented where, beside the chemical composition, influences of non-metallic inclusions on mechanical properties of steel were also elaborated. They are generally limited to individual steel grade [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17] and individual type of inclusions [2,5,6,7,8,9,10,11,12,13,14,15,16,17] especially based on rare earth elements [5,11,12,15,16,17]. At the same time, attention is paid to grain morphology, size [12,17,18], microstructure [7,11,12], and heat treatment [1,12,13,18]. Only a few industrial studies are available [9,15].
Qiu and coauthors [1] used two ferritic/martensitic steel grades where Y, Ti, and Zr were transformed into steels through vacuum induction melting. Steel was additionally tempered at 650 °C. Y, Ti, and especially Zr based particles were formed. It was shown that the presence Zr-rich inclusions caused a reduction in mechanical properties. A balance between the formation of inclusions and the tempering temperature could support the improvement of strength and impact toughness.
In the paper [18], the influence of boron and niobium on the microstructure and consequent changes in the mechanical properties of hot-rolled low-carbon bainite steel was investigated. Additionally, hardenability and phase transformation were analyzed. For microstructure analysis and non-metallic inclusion contents, laser scanning confocal microscopy and field emission scanning electron microscopy were used. It was found that the boron does not influence hardness, yield, or tensile strength. The addition of boron generally reduced the impact toughness. The authors explain the outcomes in terms of micro-cleanliness, slight microstructural changes, grain size, and precipitation of coarse (Fe,Cr)23(B,C)6.
Gong and coauthors [2] studied the addition of Mg to Cr12Mo1V1 steel to observe changes in the composition, morphology, size, and distribution of non-metallic inclusions. Added Mg formed MgO·Al2O3 and MgO inclusions as well as MnS–MgS precipitates. It was also found that due to the presence of Mg, the Mg-containing inclusions increase, consequently the density of inclusions decreases, which result in efficient inclusion removal. In addition, strength and plasticity are improved.
The effect of the content of non-metallic inclusions on the impact toughness of 32 CDV 13 steel is analyzed in the article [3]. Based on the fractography of the Charpy impact test samples, it was found that the impact toughness increases with decreasing content of non-metallic inclusions. The significant influence was attributed to micro-cracking at the grow of the inclusion-nucleated voids. Moreover, numerical simulations using COMSOL Multiphysics software were used to calculate the concentrated stresses of the steel matrix around the non-metallic inclusions. Calculated stress concentrations are higher at sulfide inclusions than oxide inclusions.
Similarly, Liu et al. [4] investigated the effect of the size and distribution of non-metal inclusions on mechanical properties of low activation martensitic steel. By optical and scanning electron microscope, the aluminium based non-metallic inclusions content was determined. The aluminium inclusions content was reduced using electroslag remelting. Accordingly, the mechanical properties were improved.
In the paper [5], cerium addition to interstitial free steel was studied in order to analyze the influences on mechanical properties. Cerium based inclusions forms Al2O3 and TiN-Al2TiO5-Al2O3 to Ce2O3, Ce2O2S, CeAlO3, TiN-Al2TiO5-Ce2O3, and TiN-Al2TiO5-Ce2O2S composite inclusions with various thermal coefficients and hardness (i.e., brittleness). It can be concluded that due to small differences between thermal coefficients and the hardness of cerium-based inclusions and iron matrix, higher mechanical properties were achieved.
Wu and coauthors [6] studied the influence of MnS inclusions morphology and spatial distribution on elongation, reduction area, yield and tensile strength. In the experiment, hot rolled bars made from steel containing 0.03% of sulphur were used. The variations of elongation and reduction area increased in the transverse direction while the tensile strength and yield strength remained practically unchanged. The MnS inclusions were classified into three types: parallel elongated MnS inclusions, densely distributed tiny MnS inclusions, and spindle-shaped MnS inclusions. Based on their parameters and probe orientation reduction area linear regression was used for reduction area prediction. Based on linear regression results densely distributed tiny MnS inclusions should be avoided.
Similarly, Lu et al. [9] investigated transverse tensile properties of hot rolled bars from two commercial steel grades while varying Ca, Al, Mg, and Mn contents in steel. Consequently, inclusions with different morphology and distribution formed. It was found that due to agmination of MnS inclusions, the transverse ductility worsened due to higher stress concentrations.
Qiu et al. [7] investigated the effect of Zr addition on the inclusion formation, microstructure, tensile properties, and impact toughness of the low activation martensitic steel. The addition of Zr varied from 0.01% to 0.11%. Scanning electron microscope was used for examination of Zr-based inclusions and Zr-based precipitates It was found out that Zr content deteriorates the mechanical properties. Accordingly, the highest yield and tensile strength values were obtained at 0.01% of Zr where Zr3V3C carbides were precipitated.
Zirconium addition was also investigated using cast heavy section steel [8]. Optical and scanning electron microscopes have been used for examination of inclusions. However, compared with the research mentioned above, the mechanical properties (Tensile and Charpy V-notch impact test) improved due to the addition of Zr. The examinations also show that the Zr addition modifies MnS inclusions.
Similarly, the effect of yttrium low activation martensitic steel was investigated [12]. By adding yttrium, AC3 temperature was increased. Afterward, the tempering process was optimized in order to control carbide precipitations. It was found that based on precipitate dispersion, carbides were smaller. On the other hand, the coarsening of martensite laths and coarsening of grains occurred after tempering. This effect decreased after quenching the low activation martensitic steel twice.
Gupta et al. [19] studied the effect of inclusions on the forming operation of steel sheets using two-dimensional finite element analysis. The focus was to understand the void formation and matrix inclusion interactions. Based on the finite element method and properties of the inclusions, the properties of steel matrix were predicted. Accordingly, the failures during hot rolling could also be determined.
The authors in [11] analysed the influences on cerium addition on carbides formation and mechanical properties of cast high grade knives steel. The size and morphology of carbides and mechanical properties were studied using optical microscopy, scanning electron microscopy, and X-ray diffraction. It was found that cerium refines cast microstructure, as well as changes the type of inclusions and carbide morphology, where rod-shaped carbides turn into lamellae. It can be concluded that cerium acts as nucleator of carbides, reducing the size and distribution of carbides in cast steel and consequently improving mechanical properties.
In the research [14], the influence of TiN inclusions on delayed cracking after flame cutting of high strength, wear resistant steel was studied. Optical microscopy, scanning electron microscopy, X-ray diffraction, transmission electron microscopy, and electron backscattered diffractometry were used. It was found out that the delayed cracking occurred from segregations during solidification and also can be attributed to micro TiN inclusions.
In their industrial study, Jonšta and coauthors [15] analyzed cerium addition to 42CrMo4 steel, which was used for forging. Loss of cerium during melting and due to reoxidation and also mechanical properties were monitored. The expected microstructure refinement and improved mechanical properties were not achieved. Moreover, the loss of cerium during steelmaking reached 50%. The locally segregated cerium inclusions caused initiation of cracks. Furthermore, several types of cerium based inclusions were detected in forged bars.
Similarly, the effect of cerium addition on the grain size, tensile, and impact properties of LDX2101 duplex stainless steel was investigated [16]. For inclusion analyses, optical and scanning electron microscopes were used. The results show the formation of different types, sizes, and distributions of inclusions affecting the grain size and consequently mechanical properties. With the addition of cerium, the energy absorbed during the impact test, tensile strength, yield strength, and elongation improved.
In this industrial study, we covered all available tensile strength data (472 cases–472 test pieces) of 17 steel grades which were ordered and used by the final user in rolled condition. The rolled bars were produced from January 2016 to December 2021 in Štore Steel Ltd., which is one of the biggest flat spring steel producers in Europe. The samples were taken from serial production, with test pieces machined using the same equipment by the same group of machine operators and the tensile testing performed by the same operator and using the same tensile test machine. The main motive was to see the influence of non-metallic inclusions on mechanical properties obtained by tensile testing. Based on the geometry of rolled bars, selected chemical composition, and average size of worst field non-metallic inclusions (sulfur, silicate, aluminium and globular oxides), several models for tensile strength, yield strength, percentage elongation, and percentage reduction area were obtained using linear regression and genetic programming. Based on developed models, the influences of individual parameters were also calculated. The modeling results were validated with improved yield and tensile strength of 30MnVS6 produced in 2022.

2. Materials and Methods

Production in Štore Steel Ltd. (Štore, Slovenia) steel plant consists of melting of the scrap using an electric arc furnace, tapping, ladle treatment (i.e., secondary metallurgy), and continuous casting of the 180 mm × 180 mm billets. The cast billets are reheated and rolled in the rolling plant using three rolling stands. First, two rolling stands are duo reversible stands (800 mm and 650 mm diameter rolls) and the last continuous rolling line (460 mm diameter rolls) has 6 horizontal and 4 vertical stands.
In this industrial study, we covered all available tensile strength data (472 cases–472 test pieces) of 17 steel grades which were ordered and used by the final user in rolled condition. The rolled bars were produced from January 2016 to December 2021 in Štore Steel Ltd. which is one of the biggest flat spring steel producers in Europe. The samples were taken from serial production.
According to ISO 14284 (Steel and iron-sampling and preparation of samples for the determination of chemical composition) immersion probes are inserted manually into the melt in tundish using a lance attached directly to the probe assembly. The chemical composition of obtained disc samples (using immersion probes) was determined according to ASTM E1019 (carbon and sulfur) ASTM E415 (for all other elements). For chemical composition determination ELTRA ONHp analyzer (Eltra GmbH, Haan, Germany), ELTRA CSi analyzer (Eltra GmbH, Haan, Germany) (according to ASTM E1019), and SPECTROLAB LAVMC12A (SPECTRO Analytical Instruments GmbH, Kleve, Germany) (according to ASTM E415) were used.
Non-metallic inclusion content was determined according to ASTM E45. According to standard at least six representative samples per lot (i.e., quantity of steel from the same batch and the same geometry) should be examined using a microscope. The minimum polished surface area of a specimen for the microscopic determination of inclusion content is 160 mm2. Based on standard, the inclusions are classified into four categories (types) based on their morphology and also on two subcategories based on their width or diameter. Categories are:
-
A: sulfide type,
-
B: alumina type,
-
C: silicate type and
-
D: globular oxide type.
All can be subcategorized into heavy and thin describing their thickness. In this study, we considered only the size of the inclusions regardless of thickness. According to standard severity numbers are used representing the size of inclusions. Moreover, method A–worst fields was used. Using method A, all the severity numbers, related with size of inclusions, at the worst field of the samples (at least six representative samples) are averaged.
Samples were obtained from rolled bars with diameter from 18 mm to 110 mm and flat bars with thickness from 8 mm to 73 mm and width from 45 mm to 200 mm from the middle billet of the rolling lot. Preparation of samples and test pieces was conducted according to ISO 377 (Steel and steel products—Location and preparation of samples and test pieces for mechanical testing).
Test probes were made according to ISO 527-2 type 1B.
The tensile testing was performed using the Messphysik Beta 300 tensile test machine (MESSPHYSIK MATERIALS TESTING GmbH, Fürstenfeld, Austria) (Figure 1) according to ISO 6892-1.
Following parameters were collected within the period from January 2016 to December 2021 in Štore Steel Ltd. for 472 cases–472 test pieces of 17 steel grades (Table 1):
-
Carbon (Carbon), silicon (Si), manganese (Mn), silicon (Si), phosphorus (P), sulfur (S), chromium (Cr), molybdenum (Mo), nickel (Ni) and aluminium (Al) content in weight%.
-
Severity number for sulfide type of inclusions (A) determined according to standard ASTM E45, method A. Based on mentioned standard the inclusions can subcategorized into heavy and thin describing their thickness. In the research only maximal severity number was used for categorization of inclusions size.
-
Severity number for alumina type of inclusions (B) determined according to standard ASTM E45, method A.
-
Severity number for silicate type of inclusions (C) determined according to standard ASTM E45, method A.
-
Severity number for globular type of inclusions (D) determined according to standard ASTM E45, method A.
-
Tensile (Rm) and yield (Rp0.2) strength, percentage elongation (Elongation) and reduction area (Reduction).
The chemical composition (average and standard deviation) of steel grades used in the research are gathered in the Table 1.
Table 1. The chemical composition (average and standard deviation) of steel grades used in the research.
Table 1. The chemical composition (average and standard deviation) of steel grades used in the research.
Steel GradeNumber
of Cases
wt % Cwt % Siwt % Mnwt % Pwt % Swt % Crwt % Mowt % Niwt % Al
AvStdAvStdAvStdAvStdAvStdAvStdAvStdAvStdAvStd
15CrNi620.1700.0000.3100.0000.5500.0000.0090.0000.0080.0001.5700.0000.0800.0001.5200.0000.0070.000
16MnCr510.160/0.280/1.100/0.012/0.021/0.890/0.060/0.150/0.018/
16MnCrS5130.1710.0140.2680.0061.1150.0520.0110.0010.0250.0040.9750.0560.0330.0100.1070.0260.0180.002
20CrMo5150.1990.0030.2370.0091.0300.0040.0140.0010.0210.0011.1800.0080.2100.0040.2010.0310.0160.000
20MnCr5310.1760.0080.2880.0120.9030.0140.0130.0010.0070.0011.0650.0080.0210.0030.0750.0080.0170.002
20MnCrS570.1910.0040.2400.0001.2800.0260.0140.0010.0230.0001.1770.0080.0400.0000.1230.0080.0190.000
20MnV6310.1700.0000.1570.0711.5460.0150.0120.0020.0230.0010.1940.0150.0340.0150.1700.0000.0220.003
28MnCrB7100.2880.0080.1580.0061.2600.0280.0130.0040.0090.0030.3970.0150.0470.0130.1330.0220.0180.001
30MnVS6600.2870.0090.6060.0281.4290.0190.0120.0020.0290.0040.1660.0280.0290.0060.0820.0210.0140.002
38MnVS650.3720.0040.5300.0001.2980.0040.0120.0020.0470.0040.1320.0180.0280.0040.0820.0180.0100.001
C22140.1670.0090.0960.0080.3180.0060.0080.0020.0060.0020.1210.0260.0210.0040.0740.0250.0230.004
C35S30.3900.0000.2700.0000.6600.0000.0130.0000.0240.0000.2000.0000.0300.0000.0800.0000.0210.000
C45140.4600.0040.2550.0140.6310.0030.0110.0010.0260.0040.2390.0290.0550.0130.1440.0150.0190.001
C45S330.4520.0040.2910.0200.6280.0040.0100.0010.0260.0020.1820.0040.0260.0120.0880.0170.0210.004
C60S1100.5830.0050.2410.0240.7120.0140.0100.0020.0270.0020.1260.0210.0220.0040.0810.0190.0190.004
P460NH310.1850.0060.0940.0281.5540.0340.0110.0020.0130.0010.1260.0080.0270.0100.1780.0080.0210.001
S355J2920.1780.0160.3170.1061.3160.0800.0100.0020.0070.0020.1370.0410.0280.0050.0780.0150.0270.002
Gathered tensile, yield strength, percentage elongation, and reduction area data are presented in Table 2.

3. Modeling of Tensile, Yield Strength, Percentage Elongation and Reduction Area

On the basis of the collected data, the prediction of mechanical properties was conducted using linear regression and genetic programming. For the fitness function, the average relative deviation between predicted and experimental data was selected. It is defined as:
Δ = i = 1 n | Q i Q i | Q i n
where n is the size of the monitored data and Qi and Q’i are the actual and the predicted mechanical property, respectively. For prediction of mechanical properties linear regression and genetic programming were used.

3.1. Modeling of Tensile, Yield Strength, Percentage Elongation and Reduction Area Using Linear Regression

The next model significantly predicts the tensile strength (p < 0.05, ANOVA) and that only 9.13% of total variances can be explained by independent variables variances (R-square). Significantly influential parameters (p > 0.05) are only chromium (Cr) and aluminium (Al) content.
Rm = 0.36 · WIDTH 1.19 · THICKNESS 61.55 · Carbon 17.14 · Si 20.57 · Mn 3964.00 · P 1366.57 · S 63.01 · Cr + 162.26 · Mo 72.24 · Ni 4137.51 · Al + 30.79 · A 13.31 · B 262.22 · C + 4.79 · D + 782.93
The average relative deviation from experimental data is 14.78%.
In addition, the influences of individual parameters on the tensile strength while separately changing individual parameter within the individual parameter range were calculated (Figure 2). It is possible to conclude that according to linear regression results, the most influential are silicate type inclusions (C).
The next model significantly predicts the yield strength (p < 0.05, ANOVA) but only 7.15% of total variances can be explained by independent variables variances (R-square). Significantly influential parameters (p > 0.05) are hot rolled bar width (WIDTH) and thickness (THICKNESS) and chromium content (Cr).
Rp 0.2 = 1.90 · WIDTH 2.37 · THICKNESS 79.04 · Carbon 11.58 · Si 6.31 · Mn + 2463.53 · P 552.76 · S + 55.44 · Cr + 119.89 · Mo + 47.64 · Ni 2211.32 · Al + 16.89 · A 5.43 · B 156.01 · C + 1.05 · D + 509.09
The average relative deviation from experimental data is 16.02%.
Moreover, the influences of individual parameters on the yield strength while separately changing individual parameter within the individual parameter range were calculated (Figure 3). It is possible to conclude that according to linear regression results the most influential are silicate type of inclusions (C).
The next model significantly predicts the percentage elongation (p < 0.05, ANOVA) but only 7.55% of total variances can be explained by independent variables variances (R-square). Significantly influential parameters (p > 0.05) are hot rolled bar width (WIDTH) and thickness (THICKNESS) and chromium content (Cr).
Elongation = 0.039 · WIDTH 0.024 · THICKNESS + 1.13 · Carbon + 2.05 · Si + 0.32 · Mn + 138.87 · P 9.87 · S + 1.22 · Cr 6.62 · Mo 2.21 · Ni + 112.86 · Al 0.43 · A + 0.89 · B 5.50 · C + 0.33 · D + 14.12
The average relative deviation from experimental data is 25.47%.
Moreover, the influences of individual parameters on the percentage elongation while separately changing individual parameter within the individual parameter range were calculated (Figure 4). It is possible to conclude that according to linear regression results the most influential are alumina and silicate type of inclusions (C).
The next model significantly predicts the percentage reduction area (p < 0.05, ANOVA) but only 14.40% of total variances can be explained by independent variables variances (R-square). Significantly influential parameters (p > 0.05) are hot rolled bar width (WIDTH) and thickness (THICKNESS), aluminium content (Al) and aluminium type of inclusions (B).
Reduction = 0.49 · WIDTH 0.39 · THICKNESS + 8.89 · Carbon + 4.93 · Si 2.02 · Mn + 387.25 · P + 144.12 · S + 4.67 · Cr 53.03 · Mo + 20.40 · Ni + 603.89 · Al + 3.61 · A + 2.25 · B + 41.93 · C + 3.51 · D + 13.47
The average relative deviation from experimental data is 31.84%.
Moreover, the influences of individual parameters on the percentage reduction area while separately changing individual parameter within the individual parameter range were calculated (Figure 5). It is possible to conclude that according to linear regression results the most influential are hot rolled bar width (WIDTH) and thickness (THICKNESS), nickel content, and silicate type of inclusions (C).

3.2. Modeling of Tensile, Yield Strength, Percentage Elongation and Reduction Area Using Genetic Programming

In genetic programming, which is one of the most general evolutionary optimization methods, the organisms that undergo adaptation are in fact mathematical expressions (models) [20]. These models consist of the selected function (e.g., basic arithmetical functions) and terminal genes (e.g., independent input parameters, and random floating-point constants). Typical function genes are: addition (+), subtraction (−), multiplication (*) and division (/), and terminal genes (e.g., x, y, z).
Random computer programs for calculating various forms and lengths are generated by means of the selected genes at the beginning of the simulated evolution. The varying of the computer programs is carried out by means of genetic operations (e.g., crossover, mutation) during several iterations, called generations. After the completion of the variation of the computer programs, a new generation is obtained. Each result obtained from an individual program from a generation is evaluated with a comparison to the experimental data. The process of changing and evaluating organisms is repeated until the termination criterion of the process is fulfilled. Figure 6 shows the pseudo code for genetic programming system.
In-house genetic programming system [21,22,23,24], programmed using AutoLISP, which is integrated into AutoCAD (AutoCAD Release 14, Autodesk, San Rafael, CA, USA) (i.e., commercial computer-aided design software), was used. Its settings were:
-
size of the population of organisms: 1000,
-
maximum number of generations: 100,
-
reproduction probability: 0.4,
-
crossover probability: 0.6,
-
maximum permissible depth in the creation of the population: 5,
-
maximum permissible depth after the operation of crossover of two organisms: 30 and
-
smallest permissible depth of organisms in generating new organisms: 2.
Genetic operations of reproduction and crossover were used. For selection of organisms the tournament method with tournament size 7 was used.
The AutoLISP based in-house genetic programming system was run 100 times in order to develop 100 independent civilizations for each individual mechanical property. Each run for tensile, yield strength, percentage elongation, and reduction area lasted approximately 16 min and 57 s, 14 min and 23 s, 17 min and 27 s, 9 min, and 8 s on an I7 Intel processor with 8 GB of RAM, respectively.
The best mathematical model for prediction of tensile strength obtained from 100 runs of genetic programming system is:
R m = ( A + D ) ( 69.4449 D + WIDTH ) D ( ( 8.33336 ( 8.33336 + A ) + ( C a r b o n + D ) ( 8.33336 ( A + D ) THICKNESS + WIDTH ) ) ) ( A + Mn ) ( 69.4449 + WIDTH A + C + 10.3334 D + THICKNESS ( 8.33336 C + D ) ( C + C a r b o n + D + THICKNESS ) + ( A + D ) ( 69.4449 + 10.3334 D + THICKNESS ) D + WIDTH ) Mn ( 69.4449 D + WIDTH THICKNESS + ( A + C ) 2 ( C 8.33336 D + THICKNESS ) WIDTH + ( D + 8.33336 D ( C + D ) ) ( 69.4449 C D + THICKNESS + WIDTH ) WIDTH ) + D ( D C + THICKNESS + ( A + D ) ( ( ) ) C )
Its average relative deviation from experimental data is 15.73%.
Similarly, as with linear regression, individual parameters influence tensile strength while separately changing individual parameters within the individual parameter range were calculated also (Figure 7). It is possible to conclude that according to genetic programming results, the most influential are silicate type inclusions (C).
The best mathematical model for prediction of yield strength obtained from 100 runs of genetic programming system is:
R p 0.2 = ( 101.14 + 20 A + 3 Mo + 3 . P Mo + ( 67.4269 + 5 A + Mo + P WIDTH ) A + 2 Mo 6 WIDTH + 6.74269 ( 33.7134 A Mo + ( 14.4854 + A + C a r b o n ) A WIDTH + WIDTH ) + C ( 20.2281 2 Mo 1 . P Mo + 1 . ( 67.4269 + 5 A + Mo + P WIDTH ) A + 3 WIDTH ) ) 2 THICKNESS ) )
Its average relative deviation from experimental data is 15.24%.
Moreover, the influences of individual parameters on the yield strength while separately changing individual parameter within the individual parameter range were calculated (Figure 8). It is possible to conclude that according to genetic programming results, the most influential are silicate type of inclusions (C).
The best mathematical model for prediction of percentage elongation obtained from 100 runs of genetic programming system is: Metals 12 01343 i001
Its average relative deviation from experimental data is 22.72%.
Moreover, the influences of individual parameters on the percentage elongation while separately changing individual parameter within the individual parameter range were calculated (Figure 9). It is possible to conclude that according to genetic programming results the most influential are silicon content, sulfide (A) and alumina (B) type inclusions.
The best mathematical model for prediction of percentage reduction area obtained from 100 runs of genetic programming system is:
Reduction = 15.6976 A + B + C A + Mn + A Si + Si WIDTH ( 10.1939 B A Si Si 2 ( A Si + Si 2 ) A SiWIDTH ) ( 5.09697 + B C + Mn MnSi A 5.09697 2 A + Mn MnSi Si ( 5.09697 2 A Si ) ) A B C a r b o n + Mn + A Si + Si WIDTH + WIDTH A 5.09697 A + B + Mn Si ( 5.09697 A + B + Mn MnSi A 5.09697 2 A + Mn MnSi Si ( 5.09697 2 A Si ) ) C A + Mn + C A + Mn + A Si + Si WIDTH + Si WIDTH Si ( A + Mn + Si ( A Mn + MnSi + Si ( 5.09697 2 A Si ) ) + WIDTH ) Si ( 5.09697 A + B C a r b o n + 2 Mn Si A A Si MnSi SiWIDTH ) THICKNESS ) )
Its average relative deviation from experimental data is 28.16%.
Moreover, the influences of individual parameters on the percentage reduction area while separately changing individual parameter within the individual parameter range were calculated (Figure 10). It is possible to conclude that according to genetic programming results, the most influential are hot rolled bar width (WIDTH) and thickness (THICKNESS) and aluminum content (Al).

4. Modeling Results and Validation

The 30MnVS6 steel grade is used in automotive industry for highly stressed parts because the mechanical properties of forged parts using micro-alloyed steels are obtained with the grain refinement using vanadium addition. Consequently, the precipitation hardening enables the elimination of the additional heat treatment processes [25,26]. It is also the third best-selling steel grade at In Štore Steel Ltd., which is used by the final user in rolled condition.
Based on linear regression and genetic programming results, it is possible to conclude that the hot rolled bar width (WIDTH) and thickness (THICKNESS) and silicate type inclusions (C) are most influential. Unfortunately, the geometrical features are defined by the final customer and only the content of silicate type of inclusions could be changed. Based on modeling results (linear regression and genetic programming) the tensile strength decreases with increasing content of silicate type of inclusions. Deductively, we can conclude that by increasing of the silicon content, we increase the content of silicate type inclusions (C).
In the period from January 2022 to April 2022 five successively cast batches of 30MnVS6 were produced with statistically significant reduction of content of silicon (t-test, p < 0.05). The content of silicate type inclusions (Kruskal-Wallis test, p < 0,05), yield, and tensile strength also changed statistically significantly (t-test, p < 0.05). The average yield and tensile strength increased from 458.5 MPa to 525.4 MPa and from 672.7 MPa to 754.0 MPa, respectively. It is necessary to emphasize that there were no statistically significant changes in other monitored parameters.
The chemical composition of and severity numbers for inclusions determined according to standard ASTM E45, method A and geometry and mechanical properties for five successively cast batches of 30MnVS6 with a statistically significant reduction of the content of silicon are presented in Table 3 and Table 4, respectively.
In addition, the developed linear regression and genetic programming model were used for the prediction of mechanical properties of five successively cast batches of 30MnVS6 with statistically significant reduction of content of silicon.
The average relative deviation from experimental data of the linear regression model for tensile strength, yield strength, percentage elongation, and percentage reduction area is 34.60%, 45.13%, 45.16%, and 128.42%, respectively.
The average relative deviation from experimental data of the genetic programming model for tensile strength, yield strength, percentage elongation, and percentage reduction area is 18.67%, 9.47%, 24.21%, and 14.97%, respectively, which is significantly better than when using the model linear regressions.

5. Conclusions

In this paper, the influence of non-metallic inclusions on mechanical properties obtained by tensile testing was studied. From January 2016 to December 2021, all available tensile strength data (472 cases–472 test pieces) of 17 low alloy steel grades, which were ordered and used by the final user in rolled condition, were gathered.
Samples were obtained from rolled bars with diameter from 18 mm to 110 mm and flat bars with thickness from 8 mm to 73 mm and width from 45 mm to 200 mm from the middle billet of the rolling lot.
Based on the geometry of rolled bars, selected chemical composition, and average size of worst fields non-metallic inclusions (sulfur, silicate, aluminium and globular oxides), determined according to ASTM E45, method A, several models for tensile strength, yield strength, percentage elongation, and percentage reduction area were obtained using linear regression and genetic programming.
According to linear regression and genetic programming results, it is possible to conclude that the hot rolled bar width and thickness and silicate type inclusions are most influential. Unfortunately, the geometrical features are defined by the final customer and only the content of silicate type of inclusions could be changed. Based on modeling results (linear regression and genetic programming), the tensile strength decreases with an increasing content of silicate type inclusions. Deductively, we can conclude that by increasing of the silicon content, we increase the content of silicate type inclusions.
Based on modeling results in the period from January 2022 to April 2022, five successively cast batches of 30MnVS6 were produced with a statistically significant reduction of the content of silicon (t-test, p < 0.05). The content of silicate type inclusions (Kruskal–Wallis test, p < 0,05), yield, and tensile strength also changed statistically significantly (t-test, p < 0.05). The average yield and tensile strength increased from 458.5 MPa to 525.4 MPa and from 672.7 MPa to 754.0 MPa, respectively. It is necessary to emphasize that there were no statistically significant changes in other monitored parameters.
The average relative deviation from experimental data of the linear regression model for tensile strength, yield strength, percentage elongation, and percentage reduction area for five successively cast batches of 30MnVS6 are 34.60%, 45.13%, 45.16%, and 128.42%, respectively.
The average relative deviation from experimental data of the genetic programming model for tensile strength, yield strength, percentage elongation, and percentage reduction area for five successively cast batches of 30MnVS6 is 18.67%, 9.47%, 24.21%, and 14.97%, respectively, which is significantly better than when using the model linear regressions.
In the future, the focus will be on variously heat treated samples and more than 50 numbers different steel grades which are produced in the steel plant. Moreover, the results of Charpy impact testing will be taken into account. In particular, a deep investigation of the mechanisms involved should be considered.

Author Contributions

Conceptualization, M.K. and U.Ž.; methodology, M.K. and U.Ž.; software, M.K.; validation, M.K.; formal analysis, M.K. and U.Ž.; investigation, M.K.; resources, M.K. and U.Ž.; data curation, M.K.; writing—original draft preparation, M.K. and U.Ž.; writing—review and editing, M.K. and U.Ž.; visualization, M.K.; supervision, M.K. and U.Ž. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data are not publicly available due to business confidentiality.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Qiu, G.-X.; Wei, X.-L.; Bai, C.; Miao, D.-J.; Cao, L.; Li, X.-M. Inclusion and mechanical properties of ODS-RAFM steels with Y, Ti, and Zr fabricated by melting. Nucl. Eng. Technol. 2022, 54, 2376–2385. [Google Scholar] [CrossRef]
  2. Gong, W.; Jiang, Z.; Zhang, L.; Chen, C.; Dong, Y. Influence of Mg addition on inclusions and mechanical properties of Cr12Mo1V1 steel under high pressure. Mater. Sci. Eng. A 2020, 791, 139410. [Google Scholar] [CrossRef]
  3. Arreola-Herrera, R.; Cruz-Ramírez, A.; Rivera-Salinas, J.E.; Romero-Serrano, A.; Sánchez-Alvarado, R.G. The effect of non-metallic inclusions on the mechanical properties of 32 CDV 13 steel and their mechanical stress analysis by numerical simulation. Theor. Appl. Fract. Mech. 2018, 94, 134–146. [Google Scholar] [CrossRef]
  4. Liu, S.; Huang, Q.; Li, C.; Huang, B. Influence of non-metal inclusions on mechanical properties of CLAM steel. Fusion Eng. Des. 2009, 84, 1214–1218. [Google Scholar] [CrossRef]
  5. Liu, X.; Yang, J.; Zhang, F.; Fu, X.; Li, H.; Yang, C. Experimental and DFT study on cerium inclusions in clean steels. J. Rare Earths 2021, 39, 477–486. [Google Scholar] [CrossRef]
  6. Wu, M.; Fang, W.; Chen, R.-M.; Jiang, B.; Wang, H.-B.; Liu, Y.-Z.; Liang, H.-L. Mechanical anisotropy and local ductility in transverse tensile deformation in hot rolled steels: The role of MnS inclusions. Mater. Sci. Eng. A 2019, 744, 324–334. [Google Scholar] [CrossRef]
  7. Qiu, G.-X.; Zhan, D.-P.; Cao, L.; Zhang, H.-S. Effect of zirconium on inclusions and mechanical properties of China low activation martensitic steel. J. Iron Steel Res. Int. 2021, 28, 1168–1179. [Google Scholar] [CrossRef]
  8. Aydın, E.; Arısoy, C.F.; Sesen, M.K. The Effect of Manganese Sulfide Inclusions and Zirconium Additions on the Mechanical Properties of Heavy Section Cast Steel. Int. J. Met. 2018, 12, 383–395. [Google Scholar] [CrossRef]
  9. Lu, J.; Cheng, G.; Che, J.; Wang, L.; Xiong, G. Effect of Oxides on Characteristics of MnS and Transverse Mechanical Properties in Commercial Al-Killed Non-Quenched and Tempered Steel. Met. Mater. Int. 2019, 25, 473–486. [Google Scholar] [CrossRef]
  10. Maisuradze, M.V.; Björk, T. Microstructure and Mechanical Properties of Calcium Treated 42CRMO4 Steel with Improved Machinability. Trans. Indian Inst. Met. 2022, 75, 681–690. [Google Scholar] [CrossRef]
  11. Qiao, X.; Han, X.; He, Z.; Zhuang, Z.; Yang, X.; Mao, F. Effect of cerium addition on microstructure and mechanical properties of as-cast high grade knives steel. J. Iron Steel Res. Int. 2022, preprint. [Google Scholar] [CrossRef]
  12. Qiu, G.; Zhan, D.; Li, C.; Yang, Y.; Qi, M.; Jiang, Z.; Zhang, H. Effects of Yttrium and Heat Treatment on the Microstructure and Mechanical Properties of CLAM Steel. J. Mater. Eng. Perform. 2020, 29, 42–52. [Google Scholar] [CrossRef]
  13. Brytan, Z.; Król, M.; Benedyk, M.; Pakieła, W.; Tański, T.; Dagnaw, M.J.; Snopiński, P.; Pagáč, M.; Czech, A. Microstructural and Mechanical Properties of Novel Co-Free Maraging Steel M789 Prepared by Additive Manufacturing. Materials 2022, 15, 1734. [Google Scholar] [CrossRef] [PubMed]
  14. Liu, J.; Liu, D.; Zuo, X.; Liu, L.; Yan, Q. Influence of TiN Inclusions and Segregation on the Delayed Cracking in NM450 Wear-Resistant Steel. Metals 2021, 12, 21. [Google Scholar] [CrossRef]
  15. Jonšta, P.; Jonšta, Z.; Brožová, S.; Ingaldi, M.; Pietraszek, J.; Klimecka-Tatar, D. The Effect of Rare Earth Metals Alloying on the Internal Quality of Industrially Produced Heavy Steel Forgings. Materials 2021, 14, 5160. [Google Scholar] [CrossRef] [PubMed]
  16. Gong, W.; Wang, P.; Zhang, L.; Jiang, Z. Effects of Ce on Microstructure and Mechanical Properties of LDX2101 Duplex Stainless Steel. Metals 2020, 10, 1233. [Google Scholar] [CrossRef]
  17. Chen, R.; Wang, Z.; He, J.; Zhu, F.; Li, C. Effects of Rare Earth Elements on Microstructure and Mechanical Properties of H13 Die Steel. Metals 2020, 10, 918. [Google Scholar] [CrossRef]
  18. Ali, M.; Nyo, T.; Kaijalainen, A.; Javaheri, V.; Tervo, H.; Hannula, J.; Somani, M.; Kömi, J. Incompatible effects of B and B + Nb additions and inclusions’ characteristics on the microstructures and mechanical properties of low-carbon steels. Mater. Sci. Eng. A 2021, 819, 141453. [Google Scholar] [CrossRef]
  19. Gupta, A.; Goyal, S.; Padmanabhan, K.A.; Singh, A.K. Inclusions in steel: Micro–macro modelling approach to analyse the effects of inclusions on the properties of steel. Int. J. Adv. Manuf. Technol. 2015, 77, 565–572. [Google Scholar] [CrossRef]
  20. Koza, J.R. Genetic Programming: On the Programming of Computers by Means of Natural Selection; MIT Press: Cambridge, MA, USA, 1992. [Google Scholar]
  21. Kovačič, M.; Župerl, U. Genetic programming in the steelmaking industry. Genet. Program. Evolvable Mach. 2020, 21, 99–128. [Google Scholar] [CrossRef]
  22. Babič, M.; Marinković, D.; Kovačič, M.; Šter, B.; Calì, M. A New Method of Quantifying the Complexity of Fractal Networks. Fractal Fract. 2022, 6, 282. [Google Scholar] [CrossRef]
  23. Kovačič, M.; Lešer, B.; Brezocnik, M. Modelling and optimization of sulfur addition during 70MnVS4 steelmaking: An industrial case study. Adv. Prod. Eng. Manag. 2021, 16, 253–261. [Google Scholar] [CrossRef]
  24. Kovačič, M.; Salihu, S.; Gantar, G.; Župerl, U. Modeling and Optimization of Steel Machinability with Genetic Programming: Industrial Study. Metals 2021, 11, 426. [Google Scholar] [CrossRef]
  25. Kovačič, M.; Novak, D. Prediction of the chemical non-homogeneity of 30MnVS6 billets with genetic programming. Mater. Teh. 2016, 50, 69–74. [Google Scholar] [CrossRef]
  26. Da Maia, B.I.; Futami, A.H.; De Oliveira, M.A. Vanadium Alloy Steel DIN 30MnVS6 Applied in Cold Forging Process. ISIJ Int. 2018, 58, 2318–2322. [Google Scholar] [CrossRef]
Figure 1. Messphysik Beta 300 tensile test machine.
Figure 1. Messphysik Beta 300 tensile test machine.
Metals 12 01343 g001
Figure 2. The calculated influences of individual parameters on the tensile strength while separately changing individual parameter within the individual parameter range.
Figure 2. The calculated influences of individual parameters on the tensile strength while separately changing individual parameter within the individual parameter range.
Metals 12 01343 g002
Figure 3. The calculated influences of individual parameters on the yield strength while separately changing individual parameter within the individual parameter range.
Figure 3. The calculated influences of individual parameters on the yield strength while separately changing individual parameter within the individual parameter range.
Metals 12 01343 g003
Figure 4. The calculated influences of individual parameters on the percentage elongation while separately changing individual parameter within the individual parameter range.
Figure 4. The calculated influences of individual parameters on the percentage elongation while separately changing individual parameter within the individual parameter range.
Metals 12 01343 g004
Figure 5. The calculated influences of individual parameters on the percentage reduction area while separately changing individual parameter within the individual parameter range.
Figure 5. The calculated influences of individual parameters on the percentage reduction area while separately changing individual parameter within the individual parameter range.
Metals 12 01343 g005
Figure 6. The pseudo code for genetic programming system.
Figure 6. The pseudo code for genetic programming system.
Metals 12 01343 g006
Figure 7. The calculated influences of individual parameters on the tensile strength while separately changing individual parameter within the individual parameter range.
Figure 7. The calculated influences of individual parameters on the tensile strength while separately changing individual parameter within the individual parameter range.
Metals 12 01343 g007
Figure 8. The calculated influences of individual parameters on the yield strength while separately changing individual parameter within the individual parameter range.
Figure 8. The calculated influences of individual parameters on the yield strength while separately changing individual parameter within the individual parameter range.
Metals 12 01343 g008
Figure 9. The calculated influences of individual parameters on the percentage elongation while separately changing individual parameter within the individual parameter range.
Figure 9. The calculated influences of individual parameters on the percentage elongation while separately changing individual parameter within the individual parameter range.
Metals 12 01343 g009
Figure 10. The calculated influences of individual parameters on the percentage reduction area while separately changing individual parameter within the individual parameter range.
Figure 10. The calculated influences of individual parameters on the percentage reduction area while separately changing individual parameter within the individual parameter range.
Metals 12 01343 g010
Table 2. The tensile, yield strength, percentage elongation and reduction area (average and standard deviation).
Table 2. The tensile, yield strength, percentage elongation and reduction area (average and standard deviation).
Steel GradeRm [MPA]Rp0.2 [MPa]A [%]Z [%]
AvStdAvStdAvStdAvStd
15CrNi6998.502.12623.5028.999.050.2126.053.46
16MnCr5605.00/364.00/11.00/25.00/
16MnCrS5559.5469.67370.2360.7319.403.5749.668.47
20CrMo5684.80137.10451.00110.5413.465.1235.4114.49
20MnCr5653.2975.79453.1358.3517.553.5350.8316.55
20MnCrS5685.00119.56435.4368.2815.994.2948.4315.00
20MnV6709.74125.16474.3990.4615.794.8440.3111.56
28MnCrB7661.20103.85416.6083.1017.002.9646.2211.88
30MnVS6672.72112.95457.3287.6017.544.2644.7710.74
38MnVS6684.4059.51465.0059.4713.982.2337.389.81
C22700.86200.85470.07142.4116.387.6439.0914.59
C35S633.0052.14404.6724.5418.403.8645.9316.17
C45651.21118.67400.8675.5017.954.7344.6615.70
C45S644.18162.35437.12126.0916.924.9745.0713.11
C60S665.55134.29432.1386.7316.464.5443.0213.58
P460NH654.0072.68457.6170.2716.363.1947.8210.66
S355J2644.28109.22434.9286.1517.755.0045.1413.52
Table 3. The chemical composition of five successively cast batches of 30MnVS6 and severity numbers for inclusions determined according to standard ASTM E45, method A.
Table 3. The chemical composition of five successively cast batches of 30MnVS6 and severity numbers for inclusions determined according to standard ASTM E45, method A.
Batch Numberwt %
C
wt %
Si
wt %
Mn
wt %
P
wt %
S
wt %
Cr
wt %
Mo
wt %
Ni
wt %
Al
Severity
Number A
Severity
Number B
Severity
Number C
Severity
Number D
10.280.651.440.0130.030.140.030.080.0160.51.51.50
20.270.611.410.010.0220.190.040.110.0150.51.51.50
30.290.651.430.0130.0290.150.030.120.0170.52.51.50
40.30.651.490.0130.0260.20.030.120.0190.521.50
50.270.661.420.0120.0280.120.030.080.0140.5230
Table 4. The geometry and mechanical properties for five successively cast batches of 30MnVS6 with statistically significant reduction of content of silicon.
Table 4. The geometry and mechanical properties for five successively cast batches of 30MnVS6 with statistically significant reduction of content of silicon.
Batch NumberHot Rolled Bar Width [mm]Hot Rolled Bar Thickness [mm]Tensile Strength
[MPa]
Yield Strength
[MPa]
Percentage Elongation [%]Percentage Reduction [%]
1606073551218.950.0
2606072749817.847.7
3303078356011.438.1
460607675341451.4
5323275852314.735.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kovačič, M.; Župerl, U. Modeling of Tensile Test Results for Low Alloy Steels by Linear Regression and Genetic Programming Taking into Account the Non-Metallic Inclusions. Metals 2022, 12, 1343. https://doi.org/10.3390/met12081343

AMA Style

Kovačič M, Župerl U. Modeling of Tensile Test Results for Low Alloy Steels by Linear Regression and Genetic Programming Taking into Account the Non-Metallic Inclusions. Metals. 2022; 12(8):1343. https://doi.org/10.3390/met12081343

Chicago/Turabian Style

Kovačič, Miha, and Uroš Župerl. 2022. "Modeling of Tensile Test Results for Low Alloy Steels by Linear Regression and Genetic Programming Taking into Account the Non-Metallic Inclusions" Metals 12, no. 8: 1343. https://doi.org/10.3390/met12081343

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