Next Article in Journal
Relationship of Surface and Bulk Resistivity in the Case of Mechanically Damaged Fibre Reinforced Red Ceramic Waste Aggregate Concrete
Previous Article in Journal
Thermal Properties of Bayfol® HX200 Photopolymer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Steels Susceptibility to Hydrogen Embrittlement: A Thermal Desorption Spectroscopy-Based Approach Coupled with Artificial Neural Network

Department of Mechanical Engineering, School of Engineering, Aalto University, 02150 Espoo, Finland
*
Author to whom correspondence should be addressed.
Materials 2020, 13(23), 5500; https://doi.org/10.3390/ma13235500
Submission received: 6 November 2020 / Revised: 21 November 2020 / Accepted: 26 November 2020 / Published: 2 December 2020
(This article belongs to the Section Advanced Materials Characterization)

Abstract

:
A novel approach has been developed for quantitative evaluation of the susceptibility of steels and alloys to hydrogen embrittlement. The approach uses a combination of hydrogen thermal desorption spectroscopy (TDS) analysis with recent advances in machine learning technology to develop a regression artificial neural network (ANN) model predicting hydrogen-induced degradation of mechanical properties of steels. We describe the thermal desorption data processing, artificial neural network architecture development, and the learning process beneficial for the accuracy of the developed artificial neural network model. A data augmentation procedure was proposed to increase the diversity of the input data and improve the generalization of the model. The study of the relationship between thermal desorption spectroscopy data and the mechanical properties of steel evidences a strong correlation of their corresponding parameters. A prototype software application based on the developed model is introduced and is openly available. The developed prototype based on TDS analysis coupled with ANN is shown to be a valuable engineering tool for steel characterization and quantitative prediction of the degradation of steel properties caused by hydrogen.

1. Introduction

Hydrogen embrittlement (HE) is often designated as the reason for the unexpected failure of engineering components exposed to hydrogen during manufacturing or in service. Failure of a bolted joint led to an offshore spill of approximately 400 barrels of drilling fluids, in late 2012, that was determined to be induced by environmental HE with the contribution of inappropriate use of zinc plating and improper cathodic protection [1]. Corrosion-induced HE in boiler tubes and water wall tubes was found to be a significant problem in the steam-power cogeneration plants initiated by under deposit corrosion attack [2,3]. High-temperature hydrogen attack (HTHA) is an issue causing heat exchanger failures in the gas and oil industry [4,5] by hydrogen dissociating and dissolving in steels reacting with carbon and forming methane. The HTHA mechanism is different from the one driving the low-temperature hydrogen-induced damage, but it is still controlled by hydrogen diffusion. Control of the degradation of steel properties caused by hydrogen is a complex task, since the hydrogen interaction with steels depends on microstructural properties of the materials, stress state, and environmental conditions controlling the hydrogen uptake [6,7]. Studies of HE in carbon steels have often shown the presence of a certain threshold of hydrogen concentration affecting the steel fracture mode [3,8,9]. This effect, however, depends strongly on the microstructure. For example, it does not show the same influence in duplex stainless steels where the rate of HE increases smoothly with the increase of the hydrogen concentration [10]. Djukic et al. proposed a unified model to assess HE by combining phenomenological research of hydrogen–metal interactions with a predictive model based on in-situ hydrogen concentration measurement [3]. Another way to evaluate the macroscopic mechanical response (namely, ultimate tensile strength, yield stress, and elongation to fracture) was proposed by Thankachan et al. for aluminum alloys using artificial neural network (ANN), creating a strong relationship between the input data (chemical composition, strain rate, and hydrogen charging conditions) and the targeted performance parameters of mechanical response. In addition, a strong correlation between steel hardness/strength and HE susceptibility was defined for individual steels after different heat treatment procedures [11,12,13]. The second-order effect on HE susceptibility was attributed to several factors, including chemical composition and heat treatment, that control the microstructure of the materials. These factors ultimately affect hydrogen diffusion and trapping [11]. Nevertheless, systematic research of hydrogen–metal interactions in different alloys is important for the development of a robust modeling approach for in-situ assessment and prediction of HE in actual engineering applications [3,13,14].
With the use of recent advances in machine learning technologies, we have developed and published a first conceptual approach in the prediction of HE susceptibility of steels correlating the macroscopic mechanical property changes, caused by hydrogen, with hydrogen thermal desorption spectroscopy (TDS) of steels [15]. TDS spectroscopy provides indirect assessment of the material microstructural features controlling the hydrogen diffusivity, solubility, and trapping of hydrogen in the studied steel. The correlation between hydrogen thermal desorption characteristics and material microstructural features has been under the active study in recent years [16,17,18,19]. The assessment of creep damage in heat-resistant ferritic stainless steel using thermal desorption analysis was proposed by Yamashita et al., evidencing a clear overall TDS shape variation with creep, depending on the test conditions [20]. The microstructural change associated with formation of non-metallic inclusion caused by the heat treatment procedure was studied in low-carbon ferritic steel showing the effect of vanadium carbide effective surface area on hydrogen trapping [21]. A model of TDA was proposed based on a correlation between the trap site density as a function of the effective surface area of the vanadium carbides [21]. Investigation of hydrogen trapping in lab cast Fe–C–W alloys using TDS evidenced the presence of the peaks correlated to W carbides of different sizes and compositions [22]. Worth noting is that the impact of dislocation density on hydrogen trapping capacity is a very common observation [18,19,20,21]. The parameters of stress state and environmental conditions controlling the hydrogen uptake were defined as constants at this particular stage of the research (see Figure 1). The approach was tested with different ANN model architectures evidencing a strong correlation between TDS of hydrogen and the target value of the hydrogen sensitivity parameter (HSP) standing as a measure of HE susceptibility of steels [15]. The objective of the present paper is to discuss the benefits of the approach and introduce a novel prototype software application for prediction of HSP of steels based on the recently developed ANN model. The results are evaluated considering the TDS data processing, ANN architecture development, and machine learning process towards quality prediction. The approach will promote a safer structural performance, contributing to a more sustainable operation of key engineering systems with high environmental, societal and economic impact, such as power plants.

2. Materials and Methods

2.1. Sample Preparation and Mechanical Testing

Steels with different microstructures were chosen for this study to ensure variation in hydrogen susceptibility. The chemical compositions of the studied steels are shown in Table 1. The chemical compositions of the studied high-strength steels are the same, while the microstructures and mechanical properties are different due to different heat treatment procedures applied during the manufacturing process. Microstructures of the studied ferritic–martensitic high-strength steel (FMHSS) grades and corresponding heat treatment procedures were described in detail by Hickel et al. [13]. Tensile specimens with a gauge length of 32 mm and a width of 5 mm were cut by electrical discharge machining from the 1 mm thickness steel sheets parallel to its rolling direction. The gauge length part of the specimens was mechanically polished with No.1200 emery paper.
Constant extension rate tensile (CERT) tests were performed with a 35 kN MTS desktop machine (MTS Systems Corporation, Eden Prairie, MN, USA) at a strain rate of 10−4 s−1. Tensile specimens were tested in as-supplied condition and during continuous hydrogen charging. The hydrogen charging was performed electrochemically using the 1N H2SO4 solution with 20 mg∙L−1 of CH4N2S at the potential of −1.08 V, −1.2 V, and −1.225 V for austenitic stainless steel (ASS), ferritic stainless steel (FSS), and FMHSS grades, respectively. Specimens tested during continuous hydrogen charging were pre-charged to approach the homogeneous hydrogen distribution through the specimen thickness for 72 h, 2 h, and 1 h for ASS, FSS, and FMHSS grades, respectively.
The hydrogen sensitivity parameter (HSP) was calculated from the reduction of elongation to fracture as follows:
HSP = ( ε ε H ) ε × 100 %
where ε is the elongation to fracture of the specimens tested in the as-supplied condition, and ε H is the elongation to fracture of the H-charged specimen.

2.2. Hydrogen Thermal Desorption Spectroscopy

Hydrogen was measured using a thermal desorption spectroscopy (TDS) (Aalto University, Espoo, Finland) apparatus developed at Aalto University Department of Mechanical Engineering, Finland. The TDS apparatus consists of an ultra-high vacuum (UHV) measurement chamber equipped with a vacuum furnace and mass spectrometer. The UHV chamber was coupled with the air-lock chamber, where the specimen was introduced first, and pre-pumped to an intermediate pressure of 2 × 10−5 mbar. After the intermediate pressure was achieved, the specimen was automatically transported to the UHV chamber using the specimen transport system. Measurements were performed under ultra-high vacuum conditions starting from 8 × 10−8 mbar. The hydrogen desorption rate was measured in the temperature range from room temperature (RT) to 1070 K with a linear heating rate of 10 × K∙min−1. Before measurement, the specimens were cleaned by acetone in an ultrasonic bath for 1 min and dried under helium gas flow to remove the water residuals from the specimen surface. The characteristic size of a TDS specimen is 1 × 5 × 15 mm3. Steel samples were measured in the as-supplied condition. Thus, only the metallurgical hydrogen accumulated into the steels during the manufacturing and storage was considered to address the steel microstructural response using hydrogen as a probe.

2.3. Measurement Data Processing

Fitting of the spectroscopy data was performed using five Gaussian peaks and exponential functions (see Figure 2). The equation of the Gaussian peak was introduced to the fitting algorithm as follows:
f ( x ) = a × exp 1 2 × ( x b ) 2 c 2
where a , b , and c are fitting parameters defining the peak amplitude, peak temperature position, and width, respectively.
Fitting parameters of the Gaussian peaks were collected, successively forming the TDS spectra descriptor. The experimental dataset of the TDS spectra descriptors can be generalized as follows:
x 11 x 21 x 31 x i 1 x 12 x 22 x 32 x i 2 x 13 x 23 x 33 x i 3 x 1 j x 2 j x 3 j x ij
where x ij is the descriptor parameter, i is the descriptor number from 1 to 15, and j defines the number of spectra. Descriptor numbers from 1 to 3 corresponded to the first Gaussian peak parameters defining the peak amplitude, peak temperature position, and width, respectively. The next group of three parameters corresponded to the next Gaussian peak and so on. The experimental dataset comprised 50 descriptor lines ( j = [1, 50]) of 15 parameters each. The diversity of data was increased by the generation of an artificial dataset complementing the experimental dataset. The parameters of the first Gaussian peak were chosen randomly from the uniform distribution between the minimum and maximum parameter values of the experimental dataset. The parameters of the other four Gaussian peaks were generated according to their correlation with the parameters of the first Gaussian peak for each steel grade separately. The algorithm of artificial data generation can be summarized as follows:
i = [ 1 ,   3 ] x ¯ i   ~   ( min ( x i ) , max ( x i ) ) i = [ 4 ,   15 ]   { x ¯ i   ~   ( min ( x i ) , max ( x i ) ) if ,   Corr ( x i   =   [ 1 , 3 ] , x i   =   [ 4 , 15 ] ) 0.6 x ¯ i = F ( x i   =   [ 1 , 3 ] ) if ,   Corr ( x i   =   [ 1 , 3 ] , x i   =   [ 4 , 15 ] ) > 0.6
where x ¯ i is the descriptor parameter of the artificial spectroscopy data. F ( x i ) is a linear function defined by the fitting of the relationship between x i   =   [ 1 , 3 ] and x i   =   [ 4 , 15 ] . It is worth noting that the Pearson correlation coefficients between only the same type of descriptor parameters defining the peak height, temperature position, and width were considered.
The described data augmentation procedure allowed the experimental dataset to be complemented with artificial TDS fitting parameters. A dataset with 20,000 TDS samples was generated comprising 2000 TDS samples for each steel. A normalization procedure is mandatory for the input data of the artificial neural network model [14,23]. A normalization procedure was applied individually for each fitting parameter, dividing each element in the column of the dataset matrix (see Equation (3)) by the maximum value from the column. Then, the full dataset was split into training, validation, and test datasets. The training dataset comprised 6000 TDS samples of six steels (AISI 304, AISI 441, VA1000_TM05, VA1400_TM, and two batches of VA1400_MTM steel) from the list of the studied materials with the corresponding targets of the HSP-value measured experimentally. The validation dataset comprised 4000 data pairs of the remaining four steels (AISI 409, VA1000_M05, VA1200_MTM, and one batch of VA1400_MTM steel). The test dataset comprised 10,000 data pairs of all the studied steels, where each studied steel contained 1000 data pairs.

2.4. Regression Artificial Neural Network Modeling

Deep feedforward architecture was considered for the development of the regression ANN model [15]. The architecture of the developed model is schematically shown in Figure 3. The input of the ANN model was the row vector of the normalized TDS fitting parameters, as described in the “measurement data processing” section. The size of the input vector was 15, which corresponded to the number of fitting parameters. The hidden layer comprised three densely connected layers of 668, 380, and 924 neurons (see Figure 3). A single neuron is defined by the input data multiplied by the weights ( W ), and their summation obtains an addition of bias ( b ). The weights and bias transform the input data linearly. The non-linear transformation is performed by the rectified linear unit (ReLU) activation function ( φ ), defined as f ( x ) = max ( 0 , x ) [24]. The ANN ends with the output layer comprising a single neuron without an activation function [23].
The ANN learning process was performed by the RMSprop adaptive learning rate method, a form of stochastic gradient descent proposed by Geoff Hinton [25]. The mean squared error (MSE) between the experimental and predicted HSP-values was considered as a loss score for application in the learning process of the developed ANN model (see Equation (5)), where n is a number of predictions, Y i is the experimental target data, and Y i is the predicted output by the ANN model.
Mean absolute error (MAE) is a common regression metric that allows the performance of the regression ANN model to be evaluated (see Equation (6)) [14,23]. Hyper-parameters of the ANN model such as the number of neurons in the hidden layer and number of layers were optimized first using MAE-based validation by the trial-and-error method, since the general-purpose method for definition of ANN topology did still not exist. Final tuning of the hyper-parameters was performed using the random search algorithm [26]. The minimum MAE was obtained by training the ANN model with about 150 training iterations (epochs). Best model parameters (weights and bias values) corresponding to the minimum of MAE at the validation dataset were carried out using the early stopping of the learning process.
MSE = 1 n i = 1 n ( Y i Y i ) 2
MAE = 1 n i = 1 n | Y i Y i |
Development, training, validation, and testing of the ANN model were performed in Python programming language using Keras open-source neural network library running on top of TensorFlow software for machine learning applications.

3. Results

The effect of electrochemical hydrogen charging on the elongation to fracture of the studied steels and value of the hydrogen sensitivity parameter (HSP) are summarized in Table 2. The relationship between yield stress (YS), ultimate tensile strength (UTS), and HSP were studied, evidencing an increase of the HSP with increases of the YS and UTS of the studied steels (see Figure 4). Linear regression (LR) was applied to model the relationship between the YS, UTS, and HSP. Pearson coefficient of correlation (R) was calculated between the HSP measured experimentally and those calculated using the linear regression model for YS and UTS data to be about 0.94 and 0.89, respectively. The MAE of the linear regression model was calculated for the YS and UTS data to be about 3.9% and 5.5%, respectively. Standard deviation (STD) was found to be about 2.2% and 3.5% for YS and UTS data, respectively.
Tensile toughness was calculated for the specimens tested in as-supplied ( K Ic ) and H-charged ( K IH ) conditions. The results are summarized in Table 2. In the presence of hydrogen charging, hydrogen-induced cracks developed and propagated well below the K Ic . The subcritical threshold value of tensile toughness in the presence of hydrogen could be calculated from about 5% to 38% of K Ic .
Thermal desorption spectroscopy (TDS) curve-fitting parameters were defined according to the measurement data processing approach, as described in the Methods section. R-value was calculated between the TDS curve-fitting parameters of the studied steels measured experimentally and summarized within the matrix plot, as shown in Figure 5. In addition, the matrix plot was complemented with R-values calculated between YS, UTS, and TDS curve-fitting parameters. The fitting parameters ( x 1 ,   x 2 ,   x 3 ,   ,   x 15 ) corresponded to the amplitude, peak temperature position, and width of the fitted Gaussian curves, successively, as described in the Methods section. According to the rule of thumb, the relationship exists if | R | 2 / n , where n is the number of samples [27]. This threshold value was calculated to be about | R | =   0.63. Figure 5 depicts a good correlation between YS, UTS and x 2 , x 5 , x 14 parameters corresponding to the temperature position of the first, second, and fifth Gaussian peak calculated to be | R | >   0.8. A similar rate of correlation could be observed between some of the fitting parameters such as x 2 and x 5 , x 2 and x 14 , x 5 and x 12 , etc. There were, however, parameters with relatively low strength of the cross-correlation such as x 1 , x 4 , and x 6 . It is worth mentioning that the YS and UTS had the best correlation with parameter x 12 , corresponding to the width of the fourth Gaussian peak. Correlation of the YS and UTS with Gaussian peak width was not as systematic as that observed for temperature peak positions. The correlation analysis of the fitting parameters with the HSP showed similar results, evidencing the best rate of correlation between HSP and x 2 , x 5 , x 14 , and x 12 .
The regression ANN model (depicted in Figure 3) was developed and trained according to the methodology described in the Methods section. The ANN model was trained at the training dataset comprising six steel grades and 6000 data pairs. This provided training data covering a large domain of the targeted structural steels. During training, the developed model was validated at the validation dataset comprising four additional steel grades and 4000 data pairs to find the optimal ANN model parameters corresponding to the minimum of MAE calculated at the validation dataset. After the training was complete, the HSP targets were predicted at the test dataset comprising 1000 data pairs for each studied steel to evaluate the accuracy of the developed ANN model. The test dataset included TDS fitting parameters that were not exposed during ANN model training and validation, that is essential for ANN model accuracy assessment. Worthy of note is that the test dataset contained the artificial TDS fitting parameters produced by the data augmentation procedure using the same algorithm as for training and validation datasets, which may have resulted in some data leakage affecting the accuracy assessment. R-value calculated between the experimental and predicted HSP-values showed a linear correlation of 0.99 at the test dataset (see Figure 6). MAE of the HSP prediction using the developed ANN model at the test dataset was calculated to be about 1.4%. The advantage of the proposed measurement data processing is obvious considering our previous study, where MAE was calculated to be about 2.8% and 4.5% for feed-forward and convolutional neural network architectures, respectively [15]. The detailed study of the prediction error calculated at the test dataset for individual steels showed some increase of MAE calculated for steels included in the validation compared to those used in the model learning process (see Figure 7). The standard deviation of the predicted HSP did not vary significantly between the different steels, evidencing, however, some reduction in deviation of the HSP prediction at steel grades included in the validation dataset.

4. Discussion

Susceptibility of steels to HE increases with an increase in material strength. This relationship was extensively studied evidencing that material conditions such as strength and/or hardness are primary conditions defining the susceptibility of steels to HE [11,12,13]. The results, as presented in Figure 4, confirm the statement showing a certain tendency between UTS, YS, and HSP of the studied steel grades. The linear regression modeling reveals a comparable HSP prediction quality at YS and UTS data for the studied steels, caused apparently by its high cross correlation rate, as shown in Figure 5. Brahimi et al. [11] showed the embrittlement ratio decreasing smoothly with an increase in the hardness of individual high-strength steel grades. However, different alloys can have different hydrogen sensitivities despite the similar properties of hardness/strength that complicate the accurate prediction of changes in the steel properties in the presence of hydrogen [11,13]. Such a phenomenon was attributed to a combination of factors such as chemical composition, heat treatment, aging conditions affecting the microstructure of steels, and ultimately, hydrogen transport and trapping. ANN model input data comprising chemical composition together with hydrogen charging parameters were found to be sufficient for effective prediction of YS, UTS, and ε degradation of aluminum alloys, as shown by Thankachan et al. [14]. In order to evaluate the cumulative effect of hydrogen on both stress and strain of the studied steels, the tensile toughness was calculated for steels tested in as-supplied and hydrogen charged conditions (see Table 2). Threshold tensile toughness was designated through K IH as a percentage of K Ic . Figure 8 evidences the threshold tensile toughness decreasing with the increase of the UTS of the steels. However, the deviation of the results was found to be significant for the studied batch of the specimens. One can conclude that considering the cumulative effect of hydrogen on both plasticity and strength of the studied steels the UTS is insufficient for evaluation of susceptibility of steels to hydrogen embrittlement. Quantitative assessment of the effect of hydrogen on steel toughness is of the utmost importance for the characterization of the material’s resistance to fracture and effective modeling of crack propagation in the presence of hydrogen. The characterization and modeling of a joint effect of microstructure and hydrogen on change of the steel toughness using machine learning technology is the objective of future research.
Ebihara et al. proposed a fundamental model for the reproduction of the total thermal desorption spectra of hydrogen based on the mass conservation of hydrogen in the specimen with the activation energy law for the diffusion, trapping, and detrapping processes [27,28,29,30,31]. This model is applied to simulate complex thermal desorption spectra of hydrogen with more than one peak involved [27,29]. The experimental study of correlation between the steels’ microstructure and thermal desorption spectroscopy of hydrogen reveals a certain hydrogen trapping change as a function of microstructural defects and its density [16,17,18,19,20,21]. From above, one can assume that the hydrogen thermal desorption spectra from steels and alloys incorporate the information on hydrogen diffusion and trapping controlled by the material microstructure. Figure 5 shows a strict correlation between the first, second, and fifth peak temperature positions ( x 2 , x 5 , x 14 ) and fourth peak width ( x 12 ) with YS and UTS properties of the studied steels. From above, one can assume the TDS of hydrogen incorporates the primary parameters defining the HE susceptibility of the studied steels. These fitting parameters shows also a high rate of correlation with HSP and play apparently the most important role in ANN model training. At the same time, some TDS curve-fitting parameters, such as x 1 , x 4 , x 6 , x 7 , x 8 , and x 13 , have a weak relationship with the YS, UTS, HSP, and x 2 , x 5 accordingly. However, these parameters contribute apparently to the HSP prediction quality of the studied steels. It is worth noting that the studied linear correlation between the TDS curve-fitting parameters does not describe the causation completely, and the unrelated parameters may have a complex mathematical relationship.
The regression feedforward ANN model was designed to take a complete dataset of the TDS curve fitting parameters as the input data to predict the value of HSP of the studied steels. The average MAE of the ANN model surpassed markedly that which was calculated for the LR model, despite the average MAE of LR model showing a relatively small error at the experimental dataset. However, the prediction quality of the ANN and LR models can change significantly with an increase of diversity of the experimental data. The detailed analysis of MAE calculated for individual steels (see Figure 7) from the test dataset reveals a perceptible increase of MAE at steel grades chosen for validation of the ANN model during the learning process. The highest absolute error of HSP prediction was calculated to be about 5.9% for AISI 304 and VA1400MTM steel grades representing the steel grades with the lower and higher hydrogen sensitivity parameters, respectively. The observed variation can be decreased apparently by the increase of diversity of steels chosen for the ANN model learning process. In addition, one can observe that the standard deviation increased as the HSP prediction came closer to its target value (see Figure 7). The observed effect relates apparently to the bias-variance problem that describes the challenge to minimize simultaneously both sources of error in a supervised learning algorithm. Thus, minimizing bias is accompanied usually with an increase of variance and vice versa. Proper regularization techniques may considerably improve the machine learning regression algorithm performance. For an ANN, the increase of the hidden layer width may reduce the total error caused by bias and variance [32].
The developed ANN model input data includes all the TDS curve-fitting parameters that are probably excessive due to the relatively high cross correlation of some of them (see Figure 5). Feature selection and optimization of the input data is the aim of future research. Nevertheless, the developed ANN model provides the proper prediction of the HSP of steels of different microstructures, the TDS data of which lies under the probabilistic distribution of the training data. TDS of hydrogen incorporates the information on microstructure, hydrogen transport, and trapping of the studied steels that contribute markedly to the ability of the ANN model to predict the degradation of the mechanical properties caused by hydrogen. Presently human supervision of the prediction is needed to prevent the misuse of the developed approach caused by steel conditions, hydrogen distribution, post-processing of the TDS curves, etc. We believe that with further development of the experimental database and improvement of the prediction quality, the approach will contribute to the fracture mechanics analysis of steel engineering components in the presence of sources of hydrogen. Such an engineering tool can markedly improve safety and reduce the time and labor cost of the steel examination. The developed model was deployed, and it is openly available for testing as a prototype web application with interactive control of the input hydrogen TDS fitting parameters and automatic calculation of the HSP (see Figure 9) [33]. New releases of the web application and ANN model are coming in the future with the development of the research.

5. Conclusions

The artificial neural network model predicting the susceptibility to hydrogen parameter was developed. The model was trained using a dataset of TDS curve fitting parameters measured from different steel grades. TDS was treated as a parameter of the microstructure of the studied steel, while the stress state and environmental conditions were assumed to be constant due to the same approach in tensile testing and hydrogen charging procedures applied for austenitic, ferritic, and ferritic–martensitic steel grades, respectively. The HSP was calculated from the reduction of the plasticity of the studied steels caused by hydrogen charging. The ANN model provides a reliable prediction of the HSP calculated from the reduction of elongation to fracture for steels of different microstructures, the TDS data of which lies under the probabilistic distribution of the trained data. MAE of the HSP prediction using the developed ANN model at the test dataset was calculated to be about 1.4%. Worth noting is that the measurement of plasticity change caused by hydrogen is a widely used approach for HSP calculation; however, considering the effect of hydrogen on both strength and plasticity by calculation of the threshold toughness, a cumulative analysis of the hydrogen embrittlement effect is required for reliable fracture mechanics computational modeling. The approach shows promising results in the prediction of the hydrogen embrittlement effect on elongation at fracture, promoting safer structural performance and contributing to a more sustainable operation of engineering systems.

6. Patents

A patent application for the described approach was filed to the Finnish Patent and Registration Office (FI 20197123).

Author Contributions

Conceptualization, E.M.; methodology, E.M.; software, E.M.; validation, E.M., E.F.; formal analysis, E.M., E.F.; investigation, E.M., E.F.; resources, P.V.; data curation, E.M.; writing—original draft preparation, E.M.; writing—review and editing, E.M., E.F., P.V.; visualization, E.M.; supervision, P.V.; project administration, P.V.; funding acquisition, P.V. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by the School of Engineering of Aalto University, via post-doctoral scholarship no. 9155273; Business Finland, via ISA Aalto HydroSafeSteels project no. 7743/31/2018; and the Academy of Finland, via the EARLY project no. 325108.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Greenslande, J. Available online: https://www.bsee.gov/sites/bsee.gov/files/technical-presentations/standards/greenslade-bsee-astm-f16-may-8-presentation-1.pdf (accessed on 17 November 2020).
  2. Meroufel, A.; Al-Sahari, A.; Al-Hajri, M. Severe under-deposit corrosion including hydrogen embrittlement in water wall tubes. In Proceedings of the NACE Corrosion Conference, New Orleans, LA, USA, 26–30 March 2017. [Google Scholar]
  3. Djukic, M.B.; Bakic, G.M.; Zeravcic, V.S.; Sedmak, A.; Rajicic, B. Hydrogen embrittlement of industrial components: Prediction, prevention, and models. Corrosion 2016, 72, 943–961. [Google Scholar] [CrossRef]
  4. Sanders, R.E. Chemical Process Safety, 4th ed.; Elsevier Inc.: Oxford, UK, 2015; pp. 87–123. [Google Scholar]
  5. US Chemical Safety and Hazard Investigation Board Report. 2010-01-I-WA. Available online: https://www.csb.gov/tesoro-refinery-fatal-explosion-and-fire/ (accessed on 17 November 2020).
  6. Popov, B.N.; Lee, J.-W.; Djukic, M.B. Handbook of Environmental Degradation of Materials, 3rd ed.; Elsevier Inc.: Oxford, UK, 2018; pp. 133–162. [Google Scholar]
  7. Venezuela, J.; Zhou, Q.; Liu, Q.; Zhang, M.; Atrens, A. Influence of hydrogen on the mechanical and fracture properties of some martensitic advanced high strength steels in simulated service conditions. Corros. Sci. 2016, 111, 602–624. [Google Scholar] [CrossRef] [Green Version]
  8. Revie, R.W.; Sastri, V.S.; Hoey, G.R.; Ramsingh, R.R.; Mak, D.K.; Shehata, M.T. Hydrogen-induced cracking of linepipe steels part 1—Threshold hydrogen concentration and pH. Corrosion 1993, 49, 17–23. [Google Scholar] [CrossRef]
  9. Yagodzinskyy, Y.; Malitckii, E.; Ganchenkova, M.; Binyukova, S.; Emelyanova, O.; Saukkonen, T.; Hänninen, H.; Lindau, R.; Vladimirov, P.; Moeslang, A. Hydrogen effects on tensile properties of EUROFER 97 and ODS-EUROFER steels. J. Nucl. Mater. 2014, 444, 435–440. [Google Scholar] [CrossRef]
  10. Zakroczymski, T.; Glowacka, A.; Swiatnicki, W. Effect of hydrogen concentration on the embrittlement of a duplex stainless steel. Corros. Sci. 2005, 47, 1403–1414. [Google Scholar] [CrossRef]
  11. Brahimi, S.V.; Yue, S.; Sriraman, K.R. Alloy and composition dependence of hydrogen embrittlement susceptibility in high-strength steel fasteners. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2017, 375, 20160407. [Google Scholar] [CrossRef] [PubMed]
  12. Brahimi, S.; Sriraman, K.; Yue, S. Hydrogen embrittlement characteristics of two tempered martensitic steel alloys for high-strength bolting. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2016, 231, 3214–3227. [Google Scholar] [CrossRef]
  13. Hickel, T.; Nazarov, R.; McEniry, E.; Zermout, Z.; Yagodzinsky, Y.; Hänninen, H.; Rott, O.; Thiessen, R.; Mraczek, K. Hydrogen Sensitivity of Different Advanced High Strength Microstructures (HYDRAMICROS); Publication Office of the European Union: Luxembourg, 2015. [Google Scholar]
  14. Thankachan, T.; Prakash, K.S.; Pleass, C.; Jothi, S.; Prabakaran, B.; Jothi, S. Artificial neural network to predict the degraded mechanical properties of metallic materials due to the presence of hydrogen. Int. J. Hydrogen Energy 2017, 42, 28612–28621. [Google Scholar] [CrossRef] [Green Version]
  15. Malitckii, E.; Fangnon, E.; Vilaça, P. Study of correlation between the steels susceptibility to hydrogen embrittlement and hydrogen thermal desorption spectroscopy using artificial neural network. Neural Comput. Appl. 2020, 32, 1–12. [Google Scholar] [CrossRef] [Green Version]
  16. Koyama, M.; Rohwerder, M.; Tasan, C.C.; Bashir, A.; Akiyama, E.; Takai, K.; Raabe, D.; Tsuzaki, K. Recent progress in microstructural hydrogen mapping in steels: Quantification, kinetic analysis, and multi-scale characterisation. Mater. Sci. Technol. 2017, 33, 1481–1496. [Google Scholar] [CrossRef] [Green Version]
  17. Yamasaki, S.; Manabe, T.; Hirakami, D. Analysis of Hydrogen State in Steel and Trapping Using Thermal Desorption Method. Nippon Steel & Sumitomo Metal Technical Report; Special Issue on Bar & Wire Rods 2017; No. 116; pp. 38–43. Available online: https://www.nipponsteel.com/en/tech/report/nssmc/pdf/116-08.pdf (accessed on 1 December 2020).
  18. Hirakami, D.; Yamasaki, S.; Tarui, T.; Ushioda, K. Competitive phenomenon of hydrogen trapping and carbon segregation in dislocations introduced by drawing or martensitic transformation of 0.35 mass% and 0.8 mass% C steels. ISIJ Int. 2016, 56, 359–365. [Google Scholar] [CrossRef] [Green Version]
  19. Malitckii, E.; Yagodzinskyy, Y.; Vilaça, P. Role of retained austenite in hydrogen trapping and hydrogen-assisted fatigue fracture of high-strength steels. Mater. Sci. Eng. A 2019, 760, 68–75. [Google Scholar] [CrossRef]
  20. Yamashita, H.; Komazaki, S.-I.; Sato, K.; Kimura, K. Creep damage evaluation of heat resistant ferritic stainless steel by hydrogen thermal desorption analysis. Tetsu-to-Hagane 2016, 102, 630–637. [Google Scholar] [CrossRef] [Green Version]
  21. Turk, A.; San-Martin, D.; Rivera-Diaz-Del-Castillo, P.E.J.; Galindo-Nava, E. Correlation between vanadium carbide size and hydrogen trapping in ferritic steel. Scr. Mater. 2018, 152, 112–116. [Google Scholar] [CrossRef]
  22. DePover, T.; Verbeken, K. Thermal desorption spectroscopy study of the hydrogen trapping ability of W based precipitates in a Q&T matrix. Int. J. Hydrogen Energy 2018, 43, 5760–5769. [Google Scholar] [CrossRef]
  23. Chollet, F. Deep Learning with Python; Manning Publications Co.: Shelter Island, NY, USA, 2018; pp. 1–117. [Google Scholar]
  24. Hoppensteadt, F.C.; Izhikevich, E.M. Weakly Connected Neural Networks; Springer: New York, NY, USA, 1997; pp. 103–110. [Google Scholar]
  25. Tieleman, T.; Hinton, G. Lecture 6.5-rmsprop: Divide the gradient by a running average of its recent magnitude. Coursera Neural Netw. Mach. Learn. 2012, 4, 26–30. [Google Scholar]
  26. Bergstra, J.; Bengio, Y. Random search for hyper-parameter optimization. J. Mach. Learn. Res. 2012, 13, 281–305. [Google Scholar]
  27. Krehbiel, T.C. Correlation coefficient rule of thumb. Decis. Sci. J. Innov. Educ. 2004, 2, 97–100. [Google Scholar] [CrossRef]
  28. Ebihara, K.; Kaburaki, H. Numerical modeling of thermal desorption spectra of hydrogen: A review of thermal desorption models. ISIJ Int. 2012, 52, 181–186. [Google Scholar] [CrossRef] [Green Version]
  29. Turnbull, A.; Hutchings, R.; Ferriss, D. Modelling of thermal desorption of hydrogen from metals. Mater. Sci. Eng. A 1997, 238, 317–328. [Google Scholar] [CrossRef]
  30. Enomoto, M.; Hirakami, D.; Tarui, T. Modeling thermal desorption analysis of hydrogen in steel. ISIJ Int. 2006, 46, 1381–1387. [Google Scholar] [CrossRef] [Green Version]
  31. Ebihara, K.-I.; Kaburaki, H.; Suzudo, T.; Takai, K. A Numerical study on the validity of the local equilibrium hypothesis in modeling hydrogen thermal desorption spectra. ISIJ Int. 2009, 49, 1907–1913. [Google Scholar] [CrossRef] [Green Version]
  32. Neal, B.; Mittal, S.; Baratin, A.; Tantia, V.; Scicluna, M.; Lacoste-Julien, S.; Mitliagkas, I. A modern take on the bias-variance tradeoff in neural networks. arXiv 2018, arXiv:1810.08591. [Google Scholar]
  33. Evaluation of Susceptibility of Steels to Hydrogen Embrittlement. 2020. Available online: https://wiki.aalto.fi/display/HE/HSP (accessed on 30 November 2020).
Figure 1. The overview of the conditions considered as the parameters controlling the hydrogen embrittlement in steels in the developed approach.
Figure 1. The overview of the conditions considered as the parameters controlling the hydrogen embrittlement in steels in the developed approach.
Materials 13 05500 g001
Figure 2. Examples of the thermal desorption spectroscopy of hydrogen from the as-supplied (a) austenitic stainless steel, (b) ferritic stainless steel, and (c) ferritic–martensitic high-strength steel. Spectroscopy curves are fitted with the sum of five Gaussian peaks and exponential background.
Figure 2. Examples of the thermal desorption spectroscopy of hydrogen from the as-supplied (a) austenitic stainless steel, (b) ferritic stainless steel, and (c) ferritic–martensitic high-strength steel. Spectroscopy curves are fitted with the sum of five Gaussian peaks and exponential background.
Materials 13 05500 g002
Figure 3. Deep feedforward artificial neural network architecture. Input layer, output layer, and hidden layer nodes are marked in yellow, red, and green colors, respectively. The corresponding outputs of the layer nodes are shown at the right side of the image. The input data of the nodes was multiplied with the weights ( W ), and their summation obtained an addition of bias ( b ). The weights and bias transformed the input data linearly. The non-linear transformation was performed by the rectified linear unit (ReLU) activation function ( φ ) defined as f ( x ) = max ( 0 , x ) .
Figure 3. Deep feedforward artificial neural network architecture. Input layer, output layer, and hidden layer nodes are marked in yellow, red, and green colors, respectively. The corresponding outputs of the layer nodes are shown at the right side of the image. The input data of the nodes was multiplied with the weights ( W ), and their summation obtained an addition of bias ( b ). The weights and bias transformed the input data linearly. The non-linear transformation was performed by the rectified linear unit (ReLU) activation function ( φ ) defined as f ( x ) = max ( 0 , x ) .
Materials 13 05500 g003
Figure 4. Relationship between yield stress (YS) and ultimate tensile strength (UTS) vs. hydrogen sensitivity parameter (HSP) measured experimentally.
Figure 4. Relationship between yield stress (YS) and ultimate tensile strength (UTS) vs. hydrogen sensitivity parameter (HSP) measured experimentally.
Materials 13 05500 g004
Figure 5. R-value calculated between the thermal desorption spectroscopy (TDS) curve-fitting parameters of the studied steels measured experimentally including yield stress (YS), ultimate tensile strength (UTS), and hydrogen sensitivity parameter (HSP) data.
Figure 5. R-value calculated between the thermal desorption spectroscopy (TDS) curve-fitting parameters of the studied steels measured experimentally including yield stress (YS), ultimate tensile strength (UTS), and hydrogen sensitivity parameter (HSP) data.
Materials 13 05500 g005
Figure 6. Correlation between the experimental and predicted HSP at test dataset. Box plot numbers are associated with the studied steel grades according to Table 1.
Figure 6. Correlation between the experimental and predicted HSP at test dataset. Box plot numbers are associated with the studied steel grades according to Table 1.
Materials 13 05500 g006
Figure 7. MAE and STD of predicted HSP calculated for individual steel. Labels of the plot denote steel type and ordinal number of the studied steels according to the Table 1. Steel grades presented in the training dataset are marked by the letter T. The dimension is in %.
Figure 7. MAE and STD of predicted HSP calculated for individual steel. Labels of the plot denote steel type and ordinal number of the studied steels according to the Table 1. Steel grades presented in the training dataset are marked by the letter T. The dimension is in %.
Materials 13 05500 g007
Figure 8. Threshold tensile toughness ( ( K IH · 100 ) / K Ic ) as a function of the ultimate tensile strength (UTS) calculated for the studied steels. Translucent bands around the regression line define the confidence interval of the regression.
Figure 8. Threshold tensile toughness ( ( K IH · 100 ) / K Ic ) as a function of the ultimate tensile strength (UTS) calculated for the studied steels. Translucent bands around the regression line define the confidence interval of the regression.
Materials 13 05500 g008
Figure 9. Prototype interactive application predicting the HSP using the developed ANN model. Parameters of the spectra can be adjusted using sliders. The value of the HSP changes on-line with the change of the input parameters of the TDS spectra. The developed model was deployed and it is openly available online [33].
Figure 9. Prototype interactive application predicting the HSP using the developed ANN model. Parameters of the spectra can be adjusted using sliders. The value of the HSP changes on-line with the change of the input parameters of the TDS spectra. The developed model was deployed and it is openly available online [33].
Materials 13 05500 g009
Table 1. Chemical composition of the steels selected for training and validation of the artificial neural network (ANN) models (wt.%).
Table 1. Chemical composition of the steels selected for training and validation of the artificial neural network (ANN) models (wt.%).
Steel\Type\ContentCSiMnTiMoCuCrNiNbNPS
AISI 304ASS0.0490.431.49--0.4318.28.2-0.047--
AISI 409FSS<0.03<1<10.26--11.750.5--<0.0450.025
AISI 441FSS0.0140.610.420.1380.0240.1217.70.20.393-0.03<0.015
VA1000_M05FMHSS0.1570.192.240.0020.004-0.46-0.0220.0060.0110.001
VA1000_TM05FMHSS
VA1200_MTMFMHSS
VA1400_TMFMHSS
VA1400_MTMFMHSS
ASS—austenitic stainless steel; FSS—ferritic stainless steel; FMHSS—ferritic–martensitic high-strength steel.
Table 2. The elongation to fracture, yield stress (YS), ultimate tensile strength (UTS), hydrogen sensitivity parameter (HSP), and tensile toughness obtained from constant extension rate tensile testing (CERT). Note: three separate batches of VA1400_MTM steel grade were studied [13].
Table 2. The elongation to fracture, yield stress (YS), ultimate tensile strength (UTS), hydrogen sensitivity parameter (HSP), and tensile toughness obtained from constant extension rate tensile testing (CERT). Note: three separate batches of VA1400_MTM steel grade were studied [13].
#SteelElongation to Fracture (%)YS (MPa)UTS (MPa)HSP (%) K I c
  ( M J   m 3 )
K I H
( M J   m 3 )
ε ε H
1AISI 30483.732.930974560.7558159
2AISI 40933.511.326639766.211730
3AISI 44129.17.730942473.312447
4VA1000 M056.20.841059123386.411024
5VA1000 TM057.010.77955102488.98625
6VA1200 MTM8.30.781119134790.512546
7VA1400 TM5.40.341322145993.7955
8VA1400 MTM3.40.1414601494966926
94.80.171387162696.410316
106.40.131113139197.912013
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Malitckii, E.; Fangnon, E.; Vilaça, P. Evaluation of Steels Susceptibility to Hydrogen Embrittlement: A Thermal Desorption Spectroscopy-Based Approach Coupled with Artificial Neural Network. Materials 2020, 13, 5500. https://doi.org/10.3390/ma13235500

AMA Style

Malitckii E, Fangnon E, Vilaça P. Evaluation of Steels Susceptibility to Hydrogen Embrittlement: A Thermal Desorption Spectroscopy-Based Approach Coupled with Artificial Neural Network. Materials. 2020; 13(23):5500. https://doi.org/10.3390/ma13235500

Chicago/Turabian Style

Malitckii, Evgenii, Eric Fangnon, and Pedro Vilaça. 2020. "Evaluation of Steels Susceptibility to Hydrogen Embrittlement: A Thermal Desorption Spectroscopy-Based Approach Coupled with Artificial Neural Network" Materials 13, no. 23: 5500. https://doi.org/10.3390/ma13235500

APA Style

Malitckii, E., Fangnon, E., & Vilaça, P. (2020). Evaluation of Steels Susceptibility to Hydrogen Embrittlement: A Thermal Desorption Spectroscopy-Based Approach Coupled with Artificial Neural Network. Materials, 13(23), 5500. https://doi.org/10.3390/ma13235500

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