Next Article in Journal
Performance Study on a Soft X-ray Betatron Radiation Source Realized in the Self-Injection Regime of Laser-Plasma Wakefield Acceleration
Next Article in Special Issue
Acoustic Properties Comparison of Ti6Al4V Produced by Conventional Method and AM Technology in the Aspect of Ultrasonic Structural Health Monitoring of Adhesive Joints
Previous Article in Journal
Monitoring of Structures and Infrastructures by Low-Cost GNSS Receivers
Previous Article in Special Issue
A Spatial Model for Repairing of the Dam Safety Monitoring Data Combining the Variable Importance for Projection (VIP) and Cokriging Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Classification of Grain-Oriented Electrical Steel Sheets by Magnetic Barkhausen Noise Using Time-Frequency Analysis and Selected Machine Learning Algorithms

Center for Electromagnetic Fields Engineering and High-Frequency Techniques, Faculty of Electrical Engineering, West Pomeranian University of Technology, ul. Sikorskiego 37, 70-313 Szczecin, Poland
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(23), 12469; https://doi.org/10.3390/app122312469
Submission received: 5 October 2022 / Revised: 30 November 2022 / Accepted: 1 December 2022 / Published: 6 December 2022
(This article belongs to the Special Issue Structural Health Monitoring: Latest Applications and Data Analysis)

Abstract

:
In this paper, a combination of Magnetic Barkhausen Noise (MBN) and several classical machine learning (ML) methods were used to evaluate both the grade and the magnetic directions of conventional and high grain oriented electrical sheets subjected to selected surface engineering methods. The presented analysis was conducted to compare the performance of two machine learning approaches, classical ML and deep learning (DL), in reference to the same MBN examination problem and based on the same database. Thus, during the experiment, 26 classical ML algorithms were used including decision trees, discriminant analysis, support vector machines, naïve Bayes, nearest neighbor, artificial neural networks and ensemble classifiers. The experiments were carried out considering a different number of recognized magnetic directions and hence the number of determined classes as well. The results of classification accuracy of the applied ML methods were compared with those obtained for the DL model presented in a previous paper. The highest accuracy was obtained for ML models based on artificial neural networks and ensemble bagged trees. However, the accuracy did not reach 89% in the best case—for the smallest number of determined classes. Nevertheless, the achieved results generally indicated an approx. 10 percent advantage of the deep learning model over the classical ones in terms of accuracy in each of the considered cases.

1. Introduction

Grain-Oriented (GO) electrical steels are frequently used for industrial electrical engineering applications. They are particularly important in shaping electromagnetic field distributions during the production of machines, such as transformers, in order to optimise these devices’ operational properties [1,2,3]. The magnetic parameters of these materials have been significantly improved due to technological progress in their manufacturing process [4]. Various surface engineering methods applied to GO electrical sheets, such as laser or mechanical scribing, have made it possible to reduce losses. Consequently, these techniques have an impact on the increase in efficiency of the electrical machines. GO steels are characterized by preferable magnetic performance depending on the direction of magnetization and this effect is known as magnetic anisotropy. Thus, the material contains axes of easy and hard magnetization. When an external magnetic field is applied to a ferromagnetic material in the direction parallel to the axis of easy magnetization rather than the axis of hard magnetization, the reorganization of the domain structure occurs with greater activity and dynamics. Thus, to obtain the desired properties of the designed electrical machines during production, the ability to carry out quick non-destructive verification tests of the angular magnetic properties of the steel sheet becomes crucial.
To determine the directions of the axes of easy and hard magnetization, various methods [5,6,7,8,9] have been used for example: induced magnetic field [7], torque [5] or magnetization curves [5] measurement methods. A promising method for measuring axes of easy and hard magnetization is the magnetic Barkhausen noise method. It is related to the Barkhausen effect occurring during the process of magnetic domain structure reorganization which takes place when the AC magnetic field is applied to the ferromagnetic material. The Barkhausen effect has a stochastic, non-stationary nature and its course strictly depends on the microstructure and mechanical stress of the tested material [10]. The rearrangement of the domain structure causes local abrupt changes of the magnetization in the sample. When a coil is applied to such a material, a corresponding voltage will be induced in the coil, called the MBN signal [11]. Considering the close relationship between the mechanical and magnetic properties of ferromagnetic materials, it becomes possible both to influence the magnetic properties of steel using engineering methods and to test those steels subsequently as well [10]. Due to these properties, the MBN method is also applicable to study the anisotropy of the magnetic properties of steel [12,13,14,15,16,17,18,19]. The subjects of the research are often steel elements used in pipelines services. Strongly directional magnetic properties in these elements, e.g., can have a large impact on the effectiveness of the non-destructive method used to test these components, such as the magnetic flux leakage [10,12,14,15]. In [12], the authors used the MBN method to investigate the origin of the magnetic easy axis and characterize the anisotropy of pipeline API X70 steel. In paper [12], the authors proposed a method for determining the magnetic easy axis affected by the roll magnetic anisotropy in API-5L steel, while in [14,15], the authors presented the results of a study on the identification of various processes distinguished in magnetization dynamics based on the analysis of an MBN course. In [17], the authors presented the results of research on the utilization of MBN for the monitoring of variation in magnetic anisotropy affected by crystallographic texture in cold rolled low carbon steel sheets. Anisotropy was also tested in grain-oriented steels. Exemplary, in [13] the authors demonstrated the possibility of conducting quick tests and showing the relationship between the initial permeability and the root mean square of the MBN signal.
Nevertheless, given the complexity of Barkhausen noise signals, the application of data processing techniques becomes a critical step in obtaining often hidden and complex components of knowledge about the state of the material under study. The majority of MBN analyses are performed in the time domain T, however the frequency domain F is also considered [20]. However, it should be noted that the MBN signal changes both in time and frequency, and the nature and dynamics of the Barkhausen effect are dependent on the type of material and the homogeneity of its structure. As a result, it becomes crucial to track changes in the dynamics of the MBN course over time. Thus, a variety of analytic approaches based on time-frequency (TF) transformations of the MBN signal have recently been developed [21,22,23,24,25]. The article [21] presents an analysis of the MBN TF representations obtained for the determination of the influence of case hardening on MBN characteristics. Several different transformation algorithms were applied and qualitatively compared. Studies [22] and [23] analyze MBN spectrograms to investigate the processes taking place in electrical steel sheets during magnetization and provide a quantitative description of these phenomena as a function of the measurement angle. In [23], a resultant anisotropy of the steel was analyzed based on the TF spectrograms of the entire MBN burst, while in [22] the possibility of determining the influence of various factors on the magnetization dynamics was considered by extraction of bursts subsequences. In [23] and [24] the authors discussed the effect of signal processing parameters and the influence of measurement conditions. For that reason, various magnetization currents, frequencies and shapes of the excitation signal on the obtained TF spectrograms were considered to achieve the greatest distinguishability between the states of the examined steel. In article [25], the time-frequency analysis along with the J-A model was used to study the effect of stresses. The discussed works showed that the analysis of combined time and frequency (TF) domain representations can provide information that is additional to knowledge gathered through traditional approaches. It should be noted, however, that no single feature of the signal expressed in any domain adequately explains the influence of all physical factors affecting a ferromagnetic material. As a result, typical signal analysis approaches may not be enough. Therefore, there is an increasing necessity to apply techniques of searching for complicated rules and relationships within measurement data to build knowledge about the tested material based on multiple features analyses. Two issues are important in this area. The first concerns the search for human-definable features in the measured signals. The second relates to the use of methods based on automatic search for features that do not have straight definition. Currently, methods based on artificial intelligence are very popular in research on finding patterns from data. Various methods, such as machine learning algorithms, shallow and deep artificial neural networks, are used. In the context of MBN, several works [26,27,28,29] using artificial intelligence methods to extract information from signals can be found in the literature. In [26], a Deep Convolutional Neural Network (DCNN) was used to classify various CGO and HGO electrical sheets subjected to a different method of surface engineering. In addition, the proposed deep network model allowed for the evaluation of the angular deviation of a magnetic easy axis from the rolling direction. In [27], the MBN signals expressed in the time domain and autoencoders were used to estimate the magnetocrystalline energy of the examined steel, whereas in [28] a convolutional neural network for the quantitative description of plastic deformation was proposed. Last but not least, paper [29] presented an exploratory data analysis and the application of machine learning methods to describe the surface quality of a material.
Additionally, considering the degree of the measured signal processing, the algorithms can be fed with data representing different abstraction levels. In case of the first representation type, the original raw signals are used (data of low abstraction level). The second representation refers to the features vector (data of medium abstraction level) extracted from those measured raw signals (i.e., characteristic or statistical value such as median, mean, maximum or minimum). Under the first representation type, the full information gathered within the measured signal is transferred to the processing algorithms. In the case of the second representation type, the computational effort is much lower, however, it carries a possible loss of information. In the previous paper [26], the authors presented the first type of approach, where the DCNN was utilized. Classification was performed on time-frequency spectrograms of measured MBN signals. This allowed the authors to obtain good classification results. Deep learning makes it possible to search for complex relationships between data, and especially in the case of a very large data set, it shows very good performance. In contrast to the direct feature engineering used in classical ML, the DL procedure is far more complicated, and the results are more challenging to interpret. Therefore, in this paper, to verify the performance and assess the possible advantage of the first type approach to MBN analysis, the second type of measured data representation was used. For this need, first, the TF features were extracted from the TF spectrograms obtained in [26]. Then, considering feature vectors for each measured case, a further classification approach was performed using several machine learning methods. For the broad spectrum of analysis, the classification task was realised using several frequently applied standard (classical) machine learning methods, i.e., decision trees, discriminant analysis, support vector machines, naïve Bayes, nearest neighbor, artificial neural networks and ensemble classifiers. The work was carried out on the same TF representation of measured data, thanks to which it is possible to directly compare the performance of both approaches, classical ML and DL. This type of research has not been conducted so far for the purpose of analyzing MBN signals; thus, such a comparison may be beneficial in the future development of expert systems based on this method.

2. Experimental Setup

The measurements were carried out for samples made of several grades of 3% SiFe grain-oriented electrical steel sheets. Figure 1 shows the rotation machine operating in four axes (X, Y, Z and A) with the attached transducer placed over the sample. The transducer combines two sections: an excitation section and a measurement section. The transducer is made up of a C-shaped core with an excitation coil and a coil for measuring the magnetic flux density (B-coil). The MBN signal measurement section consists of two coils wound on a ferromagnetic core (MBN-coil). The details of the transducer and the measuring system were given in [22,23,30].
Conventional (CGO) and high grain-oriented (HGO) steel grades were utilized for sample preparation. The steel sheets had different thicknesses of 0.23 mm, 0.27 mm and 0.30 mm, and were subjected to various surface engineering methods, including coating with and decoating of oxide layer, laser scribing of the surface or nitriding. A detailed description of the tested samples, along with the metallographic assessment of their properties, can be found in the previous paper [26]. The set of the inspected samples with their basic description is given in Figure 2.
For each sample, the measurements were made for different angular positions of the transducer between 0° and 180° using an angular step of 1.8° (for 0° position the main axis of the transducer is along the axis transverse (TD) to the rolling direction of steel; for 90° position the main axis of the transducer is along the rolling direction of steel—RD). The excitation signal was sinusoidal with a frequency of 10 Hz. The field strength obtained during the measurements was about 1.6 kA/m. MBN signals were collected via a data acquisition (DAQ) card (NI-PCI-6251) using a sampling frequency of 250 kS/s. Prior to the acquisition, first, the signals were bandpass filtered in the frequency range of 2 kHz to 100 kHz. Next, the gathered MBN signals were transformed into TF representation using short-time Fourier transform (STFT). It allows for the determination of the content of sinusoidal components of different frequency and phase in local sections of the base signal as it changes over time [22,23,30]. In the transformation process, the Kaiser function with a size of 512 points was used, which resulted into a resolution in time of 512 µs and a frequency of 488 Hz.

3. MBN Signals Time-Frequency Transformation and Feature Extraction

Further, the MBN TF representations SBN(t, f) were expressed in the form of spectrograms |SBN(t, f)|2. Figure 3 shows the exemplary spectrograms, obtained for MBN signals acquired for four different samples and the orientation of the transducer along the RD (90°): the CGO steel sample with thickness of 0.23 mm (C23#1), the CGO steel sample with thickness of 0.27 mm after decoating of the oxide layer (C27P#2), the CGO steel sample with thickness of 0.27 mm after nitriding process (C27N#2) and the HGO steel sample after laser scribing (H23LS#1).
The results were presented in common range and colour map was expressed in linear scale. In the case of the CGO C23#1 sample (Figure 3a), there are two clearly visible areas of increased activity. The first is in the vicinity of 0.02 s, and the next is in the subperiod from about 0.03 to about 0.04 s. Figure 3b shows a spectrogram achieved for the laser scribed HGO H23LS#1. In this case, on the obtained TF characteristics, three clearly distinguished areas of activity can be noted, which also cover a much wider time segment of the magnetic field half-period than in the previous cases. Figure 3c shows the spectrogram obtained for CGO sample C27P#2 subjected to decoating of the oxide layer. As it can be seen, the presented characteristic of activity in the TF space is similar to that presented in Figure 3a in the sense that it also shows two distinguished areas of activity. Nevertheless, in this case, the areas of this activity cover a smaller space of the TF distribution, which is especially noticeable in the case of the second activity area. In addition, the amplitudes over the spectrogram, in general, take lower values (it should be emphasized that both sheets come from different manufacturers thus are characterized by a different value of magnetic permeability and have different thickness as well, which affects the differences in received power density levels). In general, at least the three key phases of the magnetization process can be distinguished for GO steels (visible especially when magnetization direction is close to the easy magnetization axis): the first is related to the reversed domain nucleation process bonded with the magnetocrystalline anisotropy, while the second is related to the movement of the 180° domains occurring close to the main MBN peak, and the third is related to the movement of the 90° domains [22]. While in the case of HGO three areas that may refer to the described magnetization phases are clearly visible, in the other two cases (Figure 3a,c) they may be hidden. In the case of C23#1, the first visible activity subperiod is relatively compact in time while the second one is more stretched. This may indicate the overlapping of the subperiods of the second and third magnetization phases in that situation. Similarly, in the case of C27P#2, around 0.04 s, one more area with a much lower level of value can be noticed. Figure 3d shows the spectrogram of nitrided CGO steel which was first subjected to the oxide layer decoating process. This spectrogram is characterized by only a single, but wide, area of activity, which proves the significant influence of nitriding on the original structure of the material. A broader presentation of the spectrograms obtained during the research and the discussion is presented in the previous paper. Nevertheless, based on the presented exemplary distributions of |SBN(t, f)|2, changes in the course of the instantaneous dynamics and the degree of activity of the MBN can be clearly noticed depending on the steel grade and machining process.
Next, further processing is required. Thus, the feature extraction procedures were performed to quantify information expressed by the TF spectrograms and track the changes in the analysed distributions. These features [30] refer to the description of instantaneous, as well as general, changes in the time-frequency representation. They allow for the definition of 66 features which can be divided into three major groups: describing statistical properties (such as minimum, maximum, mean, median, centroid, variance, standard deviation, skewness or kurtosis, etc.); describing the dynamics of variance, uniformity of distribution or degree of disorder (i.e., the TF distribution’s shape, the energy or entropy); describing the characteristics’ values referring to the symmetry, centre shift in the TF plane, flatness, homogeneity or monotonicity, etc. The definitions and the properties of the used TF features were discussed in previous paper by the authors [30]. Considering only the features calculated from the entire spectrogram, the vector dimension was reduced to 13 selected features.
Feature extraction was performed for all spectrograms. Thus, a feature vector was obtained for each examined steel grade and each measured angular position of the transducer. Pictorially, the features can be presented using a polar chart. Such a presentation makes it possible to observe how the value of a given feature changes with the angular position of the transducer, which allows for the evaluation of the angular distribution of the magnetic properties of a given steel. Exemplary polar charts are presented in Figure 4.
The figure shows the obtained angular characteristics of the four selected features for the single selected sample (the CGO steel after laser scribing C27LS#1). Each chart was normalized to the (0–1) range. As it can be seen, these distributions have a different shape depending on the given TF feature. Figure 4a presents the distribution of the concentration measure feature BNCM. As it can be seen, from 0° to about 75° there is a monotonic increase in the value, and then after reaching maximum value, a slight decrease in the value of the feature is noticeable until an angle of 90° is reached. The BNCM feature assumes higher values for spectrograms with an even energy distribution across the entire TF plane, but at the same time it is not sensitive to small values. Therefore, as the level of MBN activity increases throughout the spectral band while approaching the vicinity of the easy magnetization axis, this feature also takes higher values. Figure 4b presents the distribution of the spectral entropy feature BNSE. The value of this feature in the angular interval between 0° and 60° maintains almost a constant close to 1, while for angle between 60° and 90°, there is a sharp decrease in the value to 0. Because BNSE can assess the degree of disorder in a distribution, when the feature reaches values close to 0 for the direction close to the axis of easy magnetization, it can be interpreted as leading the distribution to attain a higher degree of order and to accommodate the existing energy states. Figure 4c shows the coefficient of variation BNCoV feature. The value increases from 0° to about a 90° angle, where at about 90° it obtains the highest value. The BNCoV expresses the ratio of the mean and the standard deviation values and the observed monotonic increase of its value from TD to RD direction confirms previous observations. The last presented feature distribution (Figure 4d) relates to the statistical quantity of the spectrograms—the median value BNMedian. The BNMedian characteristic is well corelated with the concentration measure BNCM feature. This supports prior findings that at angles congruent with or near the RD direction, spectrogram’s value essentially grows over the whole TF space. The resulting distributions allow for the same conclusions to be made. A general indication of the directions of high and low MBN activity can be seen; these directions are also consistent with the directions of easy and hard magnetization, respectively.
Figure 5 presents the angular characteristics of a single selected feature, BNCM, for three selected samples of different grades: 0.30 mm thick CGO steel C30#1, 0.30 mm thick laser scribed CGO C30LS#1 and 0.27 mm thick nitrided CGO steel C27N#2. All distributions were normalized to the greatest value from all distributions. At first, one can notice the significant difference in the BNCM value levels between samples. Furthermore, for each sample, a different course of characteristic is obtained. In the case of C30#1 (a red-dotted chart), one can observe a monotonic, generally constant, growth of the feature’s value, starting from the TD direction and heading into the RD direction, with the maximum close to the easy magnetization axis. The characteristic of the BNCM for the C30LS#1 laser scribed CGO (yellow-dotted chart) takes much smaller values compared to the results obtained for other presented in the figure; however, as in the first case, the distribution presents a relatively narrow angular range where the greatest feature’s value occurred. In case of the nitrided CGO C27N#2 (purple-dotted line), a great increase in value of the feature in the angular range from 0° to about 60° can be noticed, while after reaching the greatest value, a gentle decrease in the value is further observed while approaching the 90° angle.
Due to the various characteristics of the features, it is necessary to apply an appropriate decision-making method. It is possible to use an expert to make this decision but then there is a human factor that can cause erroneous predictions (emotions, not-well physical condition, random happenstance). As a substitute for this, an algorithmized decision-making mechanism in the form of a machine learning model can be used.
The next part of the paper will present the methodology for conducting research using machine learning algorithms, which will be described, and a summary of the results obtained will be presented.

4. Machine Learning (ML) Algorithms and Results

Artificial intelligence algorithms are becoming increasingly powerful nowadays through the increasing computing capabilities of computers. One of the most accurate machine learning algorithms is artificial neural networks (ANN). Nowadays, considering the last decade of advancement, one can split ANN into shallow and deep networks. The difference between the two groups refers to the number of layers between the input and output, called the hidden layers. The shallow ANN are made up of much fewer layers completing a specific task, while the deep networks can have a much more complex structure, consisting of layers that can perform different functions. In this paper, the shallow ANN were used. The ANN can be categorized considering the information transferring direction, into: feed forward and backpropagation [31]. The first type is characterized by the fact that the input data goes forward after each neuron and returns the result. In the second case, the error (the difference between the expected output and the received one) is calculated and then the weights are updated to minimize the error. Over the years, the ANN models have been refined to produce good classification and regression results. The inspiration for such an algorithm originates from the connections of biological neurons. They are used in various issues, such as in ocean engineering [32], to forecast performance of solar still [33], as well as in medical applications [34]. Artificial neural networks have also been used in assessing the degradation state of a material in contact with hydrogen, which can change the mechanical properties of the material [35]. The ANN network, along with the application of an algorithm to optimize its layers, has been shown to be an effective tool for predicting product decomposition in solid waste gasification [36]. ANNs have also found use in the analysis of MBNs, for example, to assess the stress values [37,38] or hardness and case hardening depth [39]. It should be mentioned that various classical machine learning techniques, such as support vector machines (SVM), decision trees, Bayesian approaches or regression models also have been used recently as a result of improved computational capabilities. Each algorithm is based on different methods for extracting patterns from data. For example, SVM can be used to implement linear or non-linear classification. The main part of which is to set an appropriate margin between the separating hyperspace and learning examples. Meanwhile, in the case of decision trees, the decision-making process is carried out by adding more nodes and branches to the tree. This process proceeds from the perspective of information increment towards the selected feature. It is repeated until each side reaches the leaf nodes. The algorithm evolved into the random forest algorithm which is a blackbox type algorithm [40,41]. One of the simplest classifiers that uses a measure of the distance between objects in multidimensional space is k-nearest neighbor (kNN). Since it is assumed that items near to one another are similar, objects are classified in the database according to the class of their nearest neighbors. The naïve Bayesian classifier assumes that features in a class are not dependent to any other feature. New data are classified based on the highest probability of belonging to a specific class. Discriminant analysis classifies data by finding linear combinations of features. Machine learning is used nowadays in many applications, e.g., rock type recognition using self-organizing maps [42], in logistics [43] or medical applications [44]. These methods have been refined over the years. When the use of a single model becomes insufficient, the use of ensemble learning methods is applied, which refer to combining several weaker models into a stronger ensemble. In the case of this article, e.g., a bagged decision tree algorithm is used. It allows one to select individual trees and then choose the class that received the best result [40], but the algorithm is complex and difficult to interpret. These methods are supported by appropriate modifications, for example, bagging (which is short for Bootstrap AGGregatING), boosting or stacking. In the case of MBN, the machine learning algorithms have been used to assess surface quality [29], predict residual stress [20,45], detect Barkhausen large jumps [46] and predict hardness and case depth [45,47]. All these applications show the great importance of the applied ML approach on the performance of the MBN method, which confirms the need to consider the route in further applications.

4.1. Discription of Considered Set of Classical ML Algorithms

In this paper, a wide variety of classifiers were employed in the task, and several versions were evaluated in each case to provide the possibility of broad comparisons of performance between classical ML and DL approaches. Considering the properties discussed in the main section of the paper, frequently applied types of classifiers have been used to conduct experiments, that is, decision trees, discriminant analysis, SVM, naïve Bayes, kNNs, ensemble classifiers and ANNs. The decision trees and discriminant analysis are relatively fast and are not characterized by high complexity of interpretation. kNN has similar properties. On the other hand, they may not provide satisfactory efficiency of the prediction process. The naïve Bayes classifier, although it is relatively simple, often outperforms more sophisticated classifiers. Models based on SVM and ANN allow for high accuracy, but they are much more complex. Ensemble learning combining multiple base classifiers (e.g., decisions tress or nearest neighbor) enables classification rules to be strengthened to achieve higher predictor accuracy. Thus, the wide experiment was introduced to search for a promising classifier type and its configuration. The following factors were taken into consideration while choosing the settings for each variant: training speed, memory utilization, prediction accuracy and response complexity. Table 1 lists seven classifier types considered during the experiment together with the most crucial variables for each of their particular versions. In connection with the theory about the lack of free lunches in machine learning [48] and the impact of multiple elements related to the dataset, multiple classifiers were analyzed in the context of the set of features created from the spectrograms.
The experiments were realized using an application built into MATLAB software called Classification Learner [49]. The experiment was carried out for the set of configurations listed in the Table 1. The data transformation scheme is presented in Figure 6. As it was described in the previous section, first the spectrograms were parameterized, then, the feature vectors were used as the input to the ML algorithms.

4.2. Presentation of the Results Obtained for Classical ML Models

From all calculated classifiers, only those that delivered the best outcome for the test data were presented. During the learning process, the data were subjected to 10-fold cross-validation and 20% of the feature vector cases (records) in the database were used to obtain the final verification. During the experiments, the performance of the successively analysed ML algorithms was investigated. The effect of varied transducer angular step size on the accuracy obtained for given the ML algorithm was assessed in these studies. Although the basic transducer rotation step used for the measurement was 1.8°, it is possible to reduce the angular measurement resolution and increase the rotation step by only considering a portion of the measurements (portion of basic steps) during the learning process. Thus, for this reason, the analysed rotation steps that the transducer took between successive measurements were equal to 1.8°, 3.6°, 9° and 21.6°. The angular step size influences both the number of unique instances in the training set and the ML algorithm’s resolution of differentiating unique cases. Using this method, the accuracy %Acc of proper prediction for different angular resolutions of measurements was calculated according to:
% A c c = P R N o R × 100
where: PR is the number of properly predicted records and NoR is the number of all records in the database.
Table 2 presents the best result obtained for all trained ML models for four angular resolutions: 1.8° (there were 51 angles, thus 51 classes as well considered between TD and RD), 3.6° (26 classes), 9° (11 classes), 21.6° (5 classes). Finally, to determine the total number of classes for all samples, the above values of the number of angle classes should be multiplied by ten, that is by the total number of the inspected steel grades.
In the Table 2, in the row named “Method” by successive columns referring to the results achieved for different angular resolutions, the algorithm for which the best result was obtained was indicated. In the case of angular resolution 1.8°, the greatest accuracy was reached for the ensemble bagged tree algorithm, while for other resolutions, it was artificial neural network (with different setup) that was the most accurate. Generally, in case of the highest number of classes (referring to the angular resolution of 1.8°—the highest angular resolution case) the Ensembles Bagged Trees models showed an advantage of 5–10% over the naïve Bayes and kNNs methods. The ANN and SVM structures did not reach enough flexibility and failed to obtain the solution in this case. Along with the reduction of the class numbers (by reducing the angular resolution to 3.6), ANNs (especially narrow and medium neural network structures, for which the accuracies differed from each other by approx. 1% maximum) obtained 6–8% advantage over the Ensembled Subspace Discriminant and Bagged Trees, leaving kNN models far behind (which reached close to 20% of difference in accuracy in reference to the ANNs’ result). Further reduction of the class number, referring to the reduction of the angular resolution to 9 degrees, sustained the percentage advantage of the ANN models (narrow, medium and bilayered, reaching similar accuracies within 1% of each other) over the ensembled ones (Subspace Discriminant and Bagged Trees). However, in this case, the SVM models (especially the linear and quadratic ones) were up to just 3 percent behind the ANN models. Finally, in the last case (of the smallest number of classes and the angular resolution of 21.6 degrees), the order of the best models did not change, but the difference of the accuracy between them decreased and did not exceed 3% this time. Moreover, as it can be seen, the best accuracy result was obtained for an angular resolution of 21.6° and the worst for an angular resolution of 1.8°, which is as expected. Utilizing greater angular steps increases the proportion of correctly identified steel grades and magnetic directions (the angles of the transducer alignment with respect to TD and RD). This is because the classification process has been obviously simplified and the number of potential classes has been decreased (there were 510 classes for 1.8° step size, but only 50 for 21.6°). However, it is crucial to achieve both the highest possible classification accuracy and the greatest number of detected directions, which, in this case, define the resolution of the classification operation. An accuracy of 50.9% was obtained for an angular resolution of 3.6°, which was stated in the previous study [26] as the most optimal for identifying grades of GO sheets and their angles (according to the average misalignment angle between the RD and easy magnetic axis of the GO steels).

4.3. Performance Comparison of the Obtained Classical ML Models and the DL Model along with a Discussion

Next, the performance of the obtained ML models was validated by comparison to that of the DCNN previously attained [26]. The detailed configuration of utilized DCNNs have been shown in Table 3, while Table 4 presents the chosen training algorithm along with its parameters, which were selected through optimization methods. The details of the process of adjustment of the DCNN structure and training algorithm parameters were presented in [26].
Table 5 presents the accuracies (%Acc DCNN) of the DCNN model obtained in [26] for the same datasets (of TF representations of MBN measured signals) as those used in this paper. The table also includes the differences in the accuracies obtained for the DCNN and the best ML models (%Acc DCNN%Acc) presented in Table 2. As it can be seen, the greatest difference was obtained for an angle resolution of 9° and the smallest for 21.6°. Nevertheless, in each of the analyzed cases, there is a clear advantage (oscillating within 10%) of using a solution based on deep networks, i.e., based on a full representation of (raw) measurement data.
When analyzing the obtained results of approaches of both types in the context of their efficiency for the analysis of MBN signals, several aspects should be discussed. The deep learning approach requires a much bigger representation of data. Thus, it is time-consuming but allows one to capture all information, to automate the process of search and definition of features and to obtain subtle sensitivity to changes in MBN signals. On the other hand, the feature extraction process leads to a higher abstraction level of the analyzed TF data. Further, a representation with a higher level of abstraction leads to generalization of data but also to the loss of some information, which may be reflected in the loss of accuracy of the applied ML algorithms. In the case of a stochastic-type signal, such as the MBN one, this may be of much importance for the final capability of proper material assessment, especially if observed differences in signals between successive considered material states are not obvious. However, simultaneously, this approach reduces the complexity of the problem (due to generalization of data via representation of features) and shortens the analysis time as well. The performance of the ML models can be improved by combining the ML model with the feature selection procedure allowing for definition of the set of relevant features for the analysis. One should also consider optimizing the hyperparameters of the selected ML algorithms through the possible use of parameter grid search algorithms or with the help of swarm algorithms to obtain better results. Thanks to this, it is possible to reduce the distance between the two discussed approaches. Additionally, it is possible to further increase accuracy by implementing an update algorithm supporting the predictions of models (such as the one presented in the previous publication [26]).
On the other hand, the advantage of the deep machine learning approach in terms of accuracy is clearly visible, especially in the context of the growing number of analyzed base cases and labeled classes. Classic ML models reach a state of accuracy saturation faster (along with the decreasing number of considered classes). It can be seen on the basis of the growing difference between the results achieved in favor of the DL model for the angular resolution of 1.8, 3.6 and 9 degrees, respectively (see Table 5). The ML models slightly reduce this difference only after the DL model approaches the limit of 100 percent accuracy (for the resolution of 21.6° because of the closing to maximum limit).

5. Conclusions

This article presents the results of work on the application of machine learning to the analysis of the time-frequency representation of Barkhausen’s noise signals for the purpose of evaluating the angular properties of selected grain-oriented electrical steel sheets. The paper refers to two possible approaches based on a different form of representation of input signals, relating to the degree of their postprocessing. For the requirements of the first approach, in the previous work [26], the authors presented the results of a procedure for analyzing the basic representation of MBN signals (raw signals transformed into time-frequency domain) based on deep learning. This paper presents a second type: a classic approach based on the extraction of unambiguously definable features of signals, and the search for intricate rules and relationships between them. For this purpose, a group of classic machine learning algorithms with varying degrees of model complexity and computational requirements was considered. Ensembled classifiers (Bagged Trees along with the Subspace Discriminant) showed relatively good behavior in reference to the others in each analyzed case. This showed a potential in this classifier type. However, it must be underlined that the ANN and SVM models performed better in the case of a smaller number of labelled classes, but the difference was not significant. Finally, both approaches, referring to the classical ML models and the DL one, were compared. The obtained results make it possible to compare the accuracy of both approaches to MBN analysis by carrying out research on a common database. Such a comparison was not the subject of any previous work and thus is a value added to future work on this topic.
Presented in this article, the comparison of machine learning methods (based on the parameter vector) with deep learning (based on the full representation of measurement signals) shows the superiority of the latter approach in the context of the achieved accuracy. In the case of machine learning, the results in accuracy were 10.1% to 13.9% worse than those obtained from deep learning methods. This leads to the need to pay special attention to the deep learning methods that were presented in the previous paper [26]. Nevertheless, it seems that considering the high availability of computer units with relatively high computing capabilities, the possibility of building systems based on deep learning seems to have significant potential for use in measurement data analysis systems, especially of a stochastic nature, such as magnetic Barkhausen noise.

Author Contributions

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

Funding

The work was funded by Research Funds of Faculty of Electrical Engineering of WestPomeranian University of Technology, Szczecin.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to express special thanks to Paweł Kochmański from the West Pomeranian University of Technology in Szczecin in particular for valuable discussions and cooperation in the realization of metallographic research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Korzunin, G.S.; Sysolyatina, I.P.; Chstyakov, V.K.; Bulychev, O.A. The Effect of the Anisotropy of Magnetic Properties of Electrical Steel on Magnetic Losses in Power Transformer Cores. Russ. J. Nondestruct. Test. 2010, 46, 632–637. [Google Scholar] [CrossRef]
  2. Elleuch, M.; Poloujadoff, M. Anisotropy in Three-Phase Transformer Circuit Model. IEEE Trans. Magn. 1997, 33, 4319–4326. [Google Scholar] [CrossRef]
  3. Lee, J.-J.; Kwon, S.-O.; Hong, J.-P.; Ha, K.-H. Cogging Torque Analysis of the PMSM for High Performance Electrical Motor Considering Magnetic Anisotropy of Electrical Steel. World Electr. Veh. J. 2009, 3, 365–369. [Google Scholar] [CrossRef] [Green Version]
  4. Wodzyński, A.; Suliga, M.; Chwastek, K. Blachy elektrotechniczne o ziarnach zorientowanych—Wybrane zagadnienia. Pr. Inst. Elektrotechniki 2014, 267, 159–173. [Google Scholar]
  5. Cullity, B.D.; Graham, C.D. Introduction to Magnetic Materials, 2nd ed.; IEEE: Hoboken, NJ, USA; Wiley: Hoboken, NJ, USA, 2009; ISBN 978-0-471-47741-9. [Google Scholar]
  6. Guo, V.W.; Lu, B.; Wu, X.; Ju, G.; Valcu, B.; Weller, D. A Survey of Anisotropy Measurement Techniques and Study of Thickness Effect on Interfacial and Volume Anisotropies in Co/Pt Multilayer Media. J. Appl. Phys. 2006, 99, 08E918. [Google Scholar] [CrossRef]
  7. Silva, E.; Paula, A.; Leite, J.; Leite, J.; Andrade, L.; de Albuquerque, V.; Tavares, J. Detection of the Magnetic Easy Direction in Steels Using Induced Magnetic Fields. Metals 2016, 6, 317. [Google Scholar] [CrossRef] [Green Version]
  8. Sievert, J. The Measurement of Magnetic Properties of Electrical Sheet Steel—Survey on Methods and Situation of Standards. J. Magn. Magn. Mater. 2000, 215–216, 647–651. [Google Scholar] [CrossRef]
  9. Tumański, S. Modern Methods of Electrical Steel Testing—A Review. Electrotech. Rev. 2021, 1, 164–169. [Google Scholar] [CrossRef]
  10. Tumański, S. Handbook of Magnetic Measurements; CRC Press: Boca Raton, FL, USA, 2016; ISBN 978-1-4398-2952-3. [Google Scholar]
  11. Santa-aho, S.; Laitinen, A.; Sorsa, A.; Vippola, M. Barkhausen Noise Probes and Modelling: A Review. J. Nondestruct. Eval. 2019, 38, 94. [Google Scholar] [CrossRef] [Green Version]
  12. Clapham, L.; Heald, C.; Krause, T.; Atherton, D.L.; Clark, P. Origin of a Magnetic Easy Axis in Pipeline Steel. J. Appl. Phys. 1999, 86, 1574. [Google Scholar] [CrossRef]
  13. Caldas-Morgan, M.; Padovese, L.R. Fast Detection of the Magnetic Easy Axis on Steel Sheet Using the Continuous Rotational Barkhausen Method. NDT E Int. 2012, 45, 148–155. [Google Scholar] [CrossRef]
  14. Martínez-Ortiz, P.; Pérez-Benitez, J.A.; Espina-Hernández, J.H.; Caleyo, F.; Hallen, J.M. On the Estimation of the Magnetic Easy Axis in Pipeline Steels Using Magnetic Barkhausen Noise. J. Magn. Magn. Mater. 2015, 374, 67–74. [Google Scholar] [CrossRef]
  15. de Campos, M.F.; Campos, M.A.; Landgraf, F.J.G.; Padovese, L.R. Anisotropy Study of Grain Oriented Steels with Magnetic Barkhausen Noise. J. Phys. Conf. Ser. 2011, 303, 012020. [Google Scholar] [CrossRef]
  16. Pérez-Benítez, J.A.; Espina-Hernández, J.H.; Man, T.L.; Caleyo, F.; Hallen, J.M. Identification of Different Processes in Magnetization Dynamics of API Steels Using Magnetic Barkhausen Noise. J. Phys. D Appl. Phys. 2015, 48, 295002. [Google Scholar] [CrossRef]
  17. Chávez-Gonzalez, A.F.; Martínez-Ortiz, P.; Pérez-Benítez, J.A.; Espina-Hernández, J.H.; Caleyo, F. Comparison of Angular Dependence of Magnetic Barkhausen Noise of Hysteresis and Initial Magnetization Curve in API5L Steel. J. Magn. Magn. Mater. 2018, 446, 18–27. [Google Scholar] [CrossRef]
  18. Kleber, X.; Vincent, A. On the Role of Residual Internal Stresses and Dislocations on Barkhausen Noise in Plastically Deformed Steel. NDT E Int. 2004, 37, 439–445. [Google Scholar] [CrossRef]
  19. Akçao, F.Ü.; Lu, Z.; Gür, C.H. Monitoring Magnetic Anisotropy Variations in Cold-Rolled Steels by Magnetic Barkhausen Noise Method. In Proceedings of the 11th European Conference on Non-Destructive Testing (ECNDT 2014), Prague, Czech Republic, 6–10 October 2014. [Google Scholar]
  20. Sorsa, A.; Santa-aho, S.; Vippola, M.; Lepistö, T.; Leiviskä, K. Utilization of Frequency-Domain Information of Barkhausen Noise Signal in Quantitative Prediction of Material Properties. AIP Conf. Proc. 2014, 1581, 1256–1263. [Google Scholar] [CrossRef]
  21. Padovese, L.; Martin, N.; Millioz, F. Time-Frequency and Time-Scale Analysis of Barkhausen Noise Signals. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2009, 223, 577–588. [Google Scholar] [CrossRef] [Green Version]
  22. Maciusowicz, M.; Psuj, G. Time-Frequency Analysis of Barkhausen Noise for the Needs of Anisotropy Evaluation of Grain-Oriented Steels. Sensors 2020, 20, 768. [Google Scholar] [CrossRef] [Green Version]
  23. Maciusowicz, M.; Psuj, G. Use of Time-Frequency Representation of Magnetic Barkhausen Noise for Evaluation of Easy Magnetization Axis of Grain-Oriented Steel. Materials 2020, 13, 3390. [Google Scholar] [CrossRef] [PubMed]
  24. Maciusowicz, M.; Psuj, G.; Chudzik, P. Time-Frequency Representation of Magnetic Barkhausen Noise under Various Measurement Conditions. In Proceedings of the 2019 International Interdisciplinary PhD Workshop (IIPhDW), Wismar, Germany, 15–17 May 2019; pp. 24–29. [Google Scholar]
  25. Liu, J.; Tian, G.; Gao, B.; Zeng, K.; Liu, Q.; Zheng, Y. Magnetic Barkhausen Noise Transient Analysis for Microstructure Evolution Characterization with Tensile Stress in Elastic and Plastic Status. Sensors 2021, 21, 8310. [Google Scholar] [CrossRef] [PubMed]
  26. Maciusowicz, M.; Psuj, G.; Kochmański, P. Identification of Grain Oriented SiFe Steels Based on Imaging the Instantaneous Dynamics of Magnetic Barkhausen Noise Using Short-Time Fourier Transform and Deep Convolutional Neural Network. Materials 2022, 15, 118. [Google Scholar] [CrossRef] [PubMed]
  27. Perez-Montes, F.; Ortega-Labra, O.; Manh, T.L.; Perez-Benitez, J.A. Enhancing the Precision of Magnetocrystalline Anisotropy Energy Estimation from Barkhausen Noise Using a Deep Neural Network. Mater. Today Commun. 2020, 24, 101145. [Google Scholar] [CrossRef]
  28. Yan, Z.; Sun, G.; Liu, X.; Li, Y.; He, C.; Xing, Z.; Wang, X.; Zhang, Y.; Ning, M. FilterNet: A Deep Convolutional Neural Network for Measuring Plastic Deformation from Raw Barkhausen Noise Waveform. J. Magn. Magn. Mater. 2022, 555, 169330. [Google Scholar] [CrossRef]
  29. Griffin, J.M.; Mathew, J.; Gasparics, A.; Vértesy, G.; Uytdenhouwen, I.; Chaouadi, R.; Fitzpatrick, M.E. Machine-Learning Approach to Determine Surface Quality on a Reactor Pressure Vessel (RPV) Steel. Appl. Sci. 2022, 12, 3721. [Google Scholar] [CrossRef]
  30. Maciusowicz, M.; Psuj, G. Use of Time-Dependent Multispectral Representation of Magnetic Barkhausen Noise Signals for the Needs of Non-Destructive Evaluation of Steel Materials. Sensors 2019, 19, 1443. [Google Scholar] [CrossRef] [Green Version]
  31. Farizawani, A.G.; Puteh, M.; Marina, Y.; Rivaie, A. A Review of Artificial Neural Network Learning Rule Based on Multiple Variant of Conjugate Gradient Approaches. J. Phys. Conf. Ser. 2020, 1529, 022040. [Google Scholar] [CrossRef]
  32. Portillo Juan, N.; Negro Valdecantos, V. Review of the Application of Artificial Neural Networks in Ocean Engineering. Ocean Eng. 2022, 259, 111947. [Google Scholar] [CrossRef]
  33. Immanual, R.; Kannan, K.; Chokkalingam, B.; Priyadharshini, B.; Sathya, J.; Sudharsan, S.; Raghu Nath, E. Performance Prediction of Solar Still Using Artificial Neural Network. Mater. Today Proc. 2022; in press. [Google Scholar] [CrossRef]
  34. Fu, J.; Cai, W.; Zeng, B.; He, L.; Bao, L.; Lin, Z.; Lin, F.; Hu, W.; Lin, L.; Huang, H.; et al. Development and Validation of a Predictive Model for Peripherally Inserted Central Catheter-Related Thrombosis in Breast Cancer Patients Based on Artificial Neural Network: A Prospective Cohort Study. Int. J. Nurs. Stud. 2022, 135, 104341. [Google Scholar] [CrossRef]
  35. Thankachan, T.; Prakash, K.S.; David Pleass, C.; Rammasamy, D.; 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]
  36. Sun, C.; Ai, L.; Liu, T. The PSO-ANN Modeling Study of Highly Valuable Material and Energy Production by Gasification of Solid Waste: An Artificial Intelligence Algorithm Approach. Biomass Conv. Bioref. 2022. [Google Scholar] [CrossRef]
  37. Wang, P.; Ji, X.; Zhu, L.; Tian, G.; Yao, E. Stratified Analysis of the Magnetic Barkhausen Noise Signal Based on Wavelet Decomposition and Back Propagation Neural Network. Sens. Actuators A Phys. 2013, 201, 421–427. [Google Scholar] [CrossRef]
  38. Wang, P.; Zhu, L.; Zhu, Q.; Ji, X.; Wang, H.; Tian, G.; Yao, E. An Application of Back Propagation Neural Network for the Steel Stress Detection Based on Barkhausen Noise Theory. NDT E Int. 2013, 55, 9–14. [Google Scholar] [CrossRef]
  39. Jedamski, R.; Epp, J. Non-Destructive Micromagnetic Determination of Hardness and Case Hardening Depth Using Linear Regression Analysis and Artificial Neural Networks. Metals 2021, 11, 18. [Google Scholar] [CrossRef]
  40. Géron, A. Hands-on Machine Learning with Scikit-Learn, Keras, and TensorFlow: Concepts, Tools, and Techniques to Build Intelligent Systems, 2nd ed.; O’Reilly Media, Inc.: Sebastopol, CA, USA, 2019; ISBN 978-1-4920-3264-9. [Google Scholar]
  41. Raschka, S.; Mirjalili, V. Python Machine Learning: Machine Learning and Deep Learning with Python, Scikit-Learn, and TensorFlow 2, 3rd ed.; Expert Insight; Packt: Birmingham, UK; Mumbai, India, 2019; ISBN 978-1-78995-575-0. [Google Scholar]
  42. Hussain, M.; Liu, S.; Ashraf, U.; Ali, M.; Hussain, W.; Ali, N.; Anees, A. Application of Machine Learning for Lithofacies Prediction and Cluster Analysis Approach to Identify Rock Type. Energies 2022, 15, 4501. [Google Scholar] [CrossRef]
  43. Tsolaki, K.; Vafeiadis, T.; Nizamis, A.; Ioannidis, D.; Tzovaras, D. Utilizing Machine Learning on Freight Transportation and Logistics Applications: A Review. ICT Express, 2022; in press. [Google Scholar] [CrossRef]
  44. Shehab, M.; Abualigah, L.; Shambour, Q.; Abu-Hashem, M.A.; Shambour, M.K.Y.; Alsalibi, A.I.; Gandomi, A.H. Machine Learning in Medical Applications: A Review of State-of-the-Art Methods. Comput. Biol. Med. 2022, 145, 105458. [Google Scholar] [CrossRef] [PubMed]
  45. Tam, P.L.; Hammersberg, P.; Persson, G.; Olavison, J. Case Depth Evaluation of Induction-Hardened Camshaft by Using Magnetic Barkhausen Noise (MBN) Method. Nondestruct. Test. Eval. 2021, 36, 494–514. [Google Scholar] [CrossRef]
  46. Alimi, R.; Ivry, A.; Fisher, E.; Weiss, E. Machine Learning Detection Algorithm for Large Barkhausen Jumps in Cluttered Environment. IEEE Magn. Lett. 2021, 10, 1–5. [Google Scholar] [CrossRef]
  47. Sorsa, A.; Santa-aho, S.; Aylott, C.; Shaw, B.A.; Vippola, M.; Leiviskä, K. Case Depth Prediction of Nitrided Samples with Barkhausen Noise Measurement. Metals 2019, 9, 325. [Google Scholar] [CrossRef]
  48. Wolpert, D.H. The Lack of A Priori Distinctions Between Learning Algorithms. Neural Comput. 1996, 8, 1341–1390. [Google Scholar] [CrossRef]
  49. The MathWorks, Inc. MATLAB and Statistics and Machine Learning ToolboxTM; The MathWorks, Inc.: Natick, MA, USA, 2022. [Google Scholar]
Figure 1. View of the rotation machine with sample and transducer (with dimensions given in mm).
Figure 1. View of the rotation machine with sample and transducer (with dimensions given in mm).
Applsci 12 12469 g001
Figure 2. Samples and measuring procedure: (a) designation of the samples, (b) testing procedure; note: #1 and #2 in indexes refers to the producer no. 1 and no. 2 of the utilized steel sheets.
Figure 2. Samples and measuring procedure: (a) designation of the samples, (b) testing procedure; note: #1 and #2 in indexes refers to the producer no. 1 and no. 2 of the utilized steel sheets.
Applsci 12 12469 g002
Figure 3. Exemplary TF spectrograms obtained for alignment of the transducer along the rolling direction RD for four different test samples: (a) CGO steel sample with thickness of 0.23 mm (C23#1), (b) high GO steel sample after laser scribing (H23LS#1) (c) CGO steel sample with thickness of 0.27 mm after decoating of the oxide layer (C27P#2) and (d) CGO steel sample with thickness of 0.27 mm after nitriding process (C27N#2).
Figure 3. Exemplary TF spectrograms obtained for alignment of the transducer along the rolling direction RD for four different test samples: (a) CGO steel sample with thickness of 0.23 mm (C23#1), (b) high GO steel sample after laser scribing (H23LS#1) (c) CGO steel sample with thickness of 0.27 mm after decoating of the oxide layer (C27P#2) and (d) CGO steel sample with thickness of 0.27 mm after nitriding process (C27N#2).
Applsci 12 12469 g003
Figure 4. Polar charts (normalized to (0–1) range) of four exemplary time-frequency features extracted from MBN TF spectrograms obtained for selected sample C27LS#1: (a) BNCM, (b) BNSE, (c) BNCoV, (d) BNMedian.
Figure 4. Polar charts (normalized to (0–1) range) of four exemplary time-frequency features extracted from MBN TF spectrograms obtained for selected sample C27LS#1: (a) BNCM, (b) BNSE, (c) BNCoV, (d) BNMedian.
Applsci 12 12469 g004
Figure 5. Polar charts of the selected time-frequency feature (the concentration measure BNCM) extracted from MBN spectrograms obtained for various samples; the presented angular characteristics were normalized according to the greatest and smallest value in the set.
Figure 5. Polar charts of the selected time-frequency feature (the concentration measure BNCM) extracted from MBN spectrograms obtained for various samples; the presented angular characteristics were normalized according to the greatest and smallest value in the set.
Applsci 12 12469 g005
Figure 6. Scheme of data workflow.
Figure 6. Scheme of data workflow.
Applsci 12 12469 g006
Table 1. Classifiers with selected hyperparameters.
Table 1. Classifiers with selected hyperparameters.
TreeFineMediumCoarse
Maximum number of splits100204
Split criterionGini’s diversity indexGini’s diversity indexGini’s diversity index
DiscriminantLinearQuadratic
PresetLinearQuadratic
SVMLinearQuadraticCubicGaussian(with scaling)
Kernel functionLinearQuadraticCubic3xGaussian(with scaling) Kernel scale (1, 3.6, 14)
Multi-class methodOne vs. OneOne vs. OneOne vs. OneOne vs. One
KNNFineMediumCoarseCosineCubicWeighted
Number of neighbours110100101010
Distance metricEuclideanEuclideanEuclideanCosineMinkowski (cubic)Euclidean
Distance weightEqualEqualEqualEqualEqualSquared inverse
Naïve BayesGaussian Naïve Bayes
PresetGaussian Naïve Bayes
Neural NetworkNarrowMediumWideBilayeredTrilayered
PresetNarrow neural networkMedium neural networkWide neural networkBilayered Neural NetworkTrilayered Neural Network
Number of fully connected layer11123
First layer size10251001010
Second Layer sizeXXX1010
Third Layer sizeXXXX10
ActivationReLUReLUReLUReLUReLU
Iteration limit10001000100010001000
Regularization strength00000
Standardize dataYesYesYesYesYes
EnsembleBoosted TreesBagged TreesSubspace Discriminant Subspace KNNRUSBoostedTrees
PresetBoosted TreesBagged TreesSubspace Discriminant Subspace KNNRUSBoostedTrees
Ensemble method AdaBoostBagSubspaceSubspaceRUSBoost
Learner TypeDecision TreeDecision TreeDiscriminantNearest neighboursDecision Tree
Maximum number of splits2025,999XX20
Number of learners3030303030
Learning rate0.1XXX0.1
Subspace dimensionXX77X
Note: The red frames depict the algorithms and their base configuration for which the highest accuracy was obtained during the computational experiment.
Table 2. The highest accuracy score obtained for all trained ML models for four angular resolutions.
Table 2. The highest accuracy score obtained for all trained ML models for four angular resolutions.
Angular Resolution1.8°3.6°21.6°
MethodEnsemble
Bagged Tree
Neural Network—Medium Neural Network—NarrowNeural network—Narrow
%Acc [%]26.350.976.388.8
Table 3. Structure of the DCNN [26].
Table 3. Structure of the DCNN [26].
Lp.LayersSizeStrideNumber of Filters
1Image Input145 × 94 × 1N/AN/A
2Convolution60 × 601 × 132
3Batch NormalizationN/AN/AN/A
4ReLUN/AN/AN/A
5AveragePooling8 × 85 × 5N/A
6Convolutuon4 × 43 × 3128
7Batch NormalizationN/AN/AN/A
8ReLUN/AN/AN/A
9AveragePooling1 × 12 × 2N/A
10Fully ConnectedN/AN/AN/A
11SoftmaxN/AN/AN/A
12Classification OutputN/AN/AN/A
Table 4. Training algorithm and its parameters [26].
Table 4. Training algorithm and its parameters [26].
ParameterValue/Name
Training algorithmRMSProp
Initial learning rate0.005
Mini Batch Size32
Number of epoch75
Table 5. Accuracy for DCNN and accuracy difference between two approaches.
Table 5. Accuracy for DCNN and accuracy difference between two approaches.
Angle Res.%Acc DCNN [%] [26](%Acc DCNN%Acc) [%]
1.8°39.510.2
3.6°62.311.4
90.213.9
21.6°98.910.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Maciusowicz, M.; Psuj, G. Classification of Grain-Oriented Electrical Steel Sheets by Magnetic Barkhausen Noise Using Time-Frequency Analysis and Selected Machine Learning Algorithms. Appl. Sci. 2022, 12, 12469. https://doi.org/10.3390/app122312469

AMA Style

Maciusowicz M, Psuj G. Classification of Grain-Oriented Electrical Steel Sheets by Magnetic Barkhausen Noise Using Time-Frequency Analysis and Selected Machine Learning Algorithms. Applied Sciences. 2022; 12(23):12469. https://doi.org/10.3390/app122312469

Chicago/Turabian Style

Maciusowicz, Michal, and Grzegorz Psuj. 2022. "Classification of Grain-Oriented Electrical Steel Sheets by Magnetic Barkhausen Noise Using Time-Frequency Analysis and Selected Machine Learning Algorithms" Applied Sciences 12, no. 23: 12469. https://doi.org/10.3390/app122312469

APA Style

Maciusowicz, M., & Psuj, G. (2022). Classification of Grain-Oriented Electrical Steel Sheets by Magnetic Barkhausen Noise Using Time-Frequency Analysis and Selected Machine Learning Algorithms. Applied Sciences, 12(23), 12469. https://doi.org/10.3390/app122312469

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