Next Article in Journal
Design and Testing of a Fully Automatic Aquatic Plant Combing Machine for Crab Farming
Previous Article in Journal
Design and Experimentation of the Millet Combine Harvester Header
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Fault Diagnosis Method with Adaptive Artificial Gorilla Troops Optimization Optimized Variational Mode Decomposition and Support Vector Machine Parameters

School of Electrical and Information Engineering, Anhui University of Technology, Maanshan 243032, China
*
Author to whom correspondence should be addressed.
Machines 2024, 12(9), 637; https://doi.org/10.3390/machines12090637
Submission received: 8 August 2024 / Revised: 7 September 2024 / Accepted: 9 September 2024 / Published: 12 September 2024
(This article belongs to the Section Machines Testing and Maintenance)

Abstract

:
To address the issue of intelligent optimization algorithms being prone to local optima, resulting in insufficient feature extraction and low fault-type recognition rates when optimizing Variational Mode Decomposition and Support Vector Machine parameters, this paper proposes a fault diagnosis method based on an improved Artificial Gorilla Troops Optimization algorithm. The Artificial Gorilla Troops Optimization algorithm was enhanced using Logistic chaotic mapping, a linear decreasing weight factor, the global exploration strategy of the Osprey Optimization Algorithm, and the Levy flight strategy, improving its ability to escape local optima, adaptability, and convergence accuracy. This algorithm was used to optimize the parameters of Variational Mode Decomposition and Support Vector Machine for fault diagnosis. Experiments on fault diagnosis with two datasets of different sample sizes showed that the proposed method achieved a diagnostic accuracy of no less than 98% for samples of varying sizes, with stable and reliable results.

1. Introduction

Motors and bearings, as important components of modern industrial devices and machinery, are prone to various faults under the combined influence of long-term continuous operation and internal and external environmental factors. These faults can affect normal production, cause economic losses, and pose safety hazards. Therefore, diagnosing faults in motors and bearings is of great significance [1]. Fault diagnosis in industrial equipment such as motors and bearings is primarily achieved through signal analysis and the development of diagnostic models [2]. This process is typically divided into three key stages: vibration signal processing, fault feature extraction, and fault identification and classification. Among these, the fault feature extraction stage is critical, as the adequacy of the extracted features directly determines the accuracy of subsequent diagnostics. Therefore, feature extraction is considered the core of the entire fault diagnosis process.
For vibration signal feature extraction, methods such as Empirical Mode Decomposition (EMD) [3], Local Mean Decomposition (LMD) [4], Empirical Wavelet Decomposition (EWT) [5,6], and Variational Mode Decomposition (VMD) [7,8] are widely used. However, EMD suffers from mode mixing and endpoint effects, and it is sensitive to noise and subtle signal changes, resulting in poor decomposition stability [9,10]. LMD is prone to pseudo-fluctuations when processing signal endpoints, and its decomposition process is complex with high computational demands, limiting its applicability [11,12]. EWT’s decomposition performance depends on the selection and design of wavelet filters and has limited frequency resolution, making it ineffective for handling non-stationary signals [13,14]. By contrast, VMD addresses the issues of mode mixing, endpoint effects, decomposition stability, and applicability found in EMD, LMD, and EWT by using variational optimization and mode bandwidth constraints. VMD demonstrates more stable performance when decomposing nonlinear and non-stationary signals [15].
When using VMD to extract vibration signal features, it is necessary to manually select the total number of mode components, K , and the quadratic penalty factor, α . These two parameters determine the performance of VMD [16]. Since α and K typically take values within a certain range, selecting values based solely on experience or comparing a few values and then drawing conclusions lacks systematic analysis and rationality. In recent years, swarm intelligence optimization algorithms have gained widespread attention from researchers for their application in parameter optimization for Variational Mode Decomposition (VMD). These algorithms are favored for their strong global search capabilities, adaptability, and ability to overcome the limitations of experience-based parameter selection, where the chosen parameter values may not be optimal. Reference [17] proposed using the kernel estimate for mutual information (KEMI) as a fitness function and employing a genetic algorithm (GA) to optimize VMD parameters, achieving the identification of single and multiple defects in constant and variable speed bearings. However, a GA relies on the evolution of the population, requiring numerous iterations to find the optimal solution, significantly increasing computational cost in parameter optimization problems. Reference [18] proposed a hybrid diagnostic method based on vibration signals, using the Gray Wolf Optimization algorithm (GWO) to optimize VMD parameter combinations. However, GWO lacks an adaptive mechanism and cannot dynamically adjust algorithm parameters based on feedback during the search process, leading to insufficient search accuracy. Reference [19] employed the Sparrow Search Algorithm (SSA) to optimize the parameters to be determined in VMD and Maximum Correlated Kurtosis Deconvolution (MCKD), but the SSA also suffers from issues such as initial population random creation, leading to uneven population distribution and a lack of diversity. Additionally, when the warning value is less than the set safety threshold, the search ability of the discoverer gradually decreases as the algorithm iterations increase, resulting in insufficient spatial search area and low convergence efficiency. Therefore, this paper selects the Artificial Gorilla Troops Optimizer (GTO), which has diversified search strategies, and addresses its shortcomings by making improvements. The improved GTO is used to optimize the parameters α and K of VMD, and to extract vibration signal features.
Using vibration signal features for fault diagnosis requires classifying and identifying the processed signals. Choosing an appropriate diagnostic method is crucial. In the field of fault diagnosis using vibration signals, machine learning methods such as k-Nearest Neighbors (k-NN), Random Forest (RF), and Support Vector Machines (SVMs) are widely employed for diagnosing faults in motors and other equipment [20,21,22,23]. Moreover, deep learning approaches, including Convolutional Neural Networks (CNN) and Long Short-Term Memory (LSTM) networks, are also commonly used to achieve fault diagnosis [24]. Reference [25] utilized genetic algorithms to select important features from statistical features extracted from motor signals, employing three machine learning methods—K-Nearest Neighbors, decision trees, and Random Forest—for training and testing to achieve precise fault identification. Reference [26] proposed a multi-signal fault diagnosis method based on a Deep Convolutional Neural Network (DCNN), capable of learning from various types of sensor signals. By extracting features from vibration and current signals, the model achieved accurate induction motor fault diagnosis. Reference [27] applied the MPE threshold method to select appropriate VMD parameters for decomposing motor vibration signals and employed the cuckoo search algorithm to optimize the global optimal solution of the Support Vector Machine, resulting in high fault identification rates. These methods each have their respective advantages and disadvantages. In this paper, considering the background of motor and motor bearing vibration signals, the availability of fault samples, and the nonlinear and non-stationary characteristics of vibration signals, we choose the SVM classifier, which supports nonlinear classification and has good generalization ability, for fault diagnosis.
In summary, to address the issue of insufficient fault vibration signal feature extraction leading to low subsequent diagnostic accuracy, this paper presents an improvement of the Artificial Gorilla Troops Optimization algorithm using four strategies to enhance its search performance and adaptability, resulting in a robust intelligent optimization algorithm. This improved algorithm is employed to optimize the parameters of Variational Mode Decomposition and Support Vector Machine, thereby enhancing the generalization capability of the fault diagnosis model. As a result, the model achieves high diagnostic accuracy across different diagnostic targets, including bearings and motors, and different sample sizes of faults. This work provides a reliable diagnostic method for the field of fault diagnosis.

2. Improved Artificial Gorilla Force Algorithm Fusing Osprey and Levy Flight (OLGTO)

2.1. The Original Artificial Gorilla Force Algorithm

Artificial gorilla troop optimization algorithm [28] is a new type of intelligent optimization algorithm proposed by mathematizing the collective social life habits of gorillas. The algorithm mainly simulates the population migration and courtship behavior of gorillas and other life behaviors to find the optimal behavior, and it has the advantages of strong ability to find the optimal and fast convergence speed, etc. The GTO algorithm mainly consists of an exploration phase and a development phase, and the development phase contains two kinds of behaviors, namely, following the silverback gorillas and competing for the adult females.
First, the GTO is initialized by setting parameters such as N , M a x I t , β , p , and W . Here, N represents the number of gorillas in the population; M a x I t is the maximum number of iterations; p controls the migration strategy of the gorillas to unknown locations, with values ranging between 0 and 1; β is set to 0.3 to calculate the intensity of the gorillas’ violent behavior; W determines the two mechanisms in the exploitation phase, with a value of 0.5. Subsequently, the gorilla population is randomly initialized in the search space.
Second, during the exploration phase, the silverback gorilla leads the other gorillas to live in groups in the natural environment. This phase primarily involves three mechanisms for global search: migrating to unknown locations, migrating to known locations, and migrating to other gorillas’ locations. Mechanism 1 allows gorillas to randomly explore the space, with the execution condition r a n d 1 < p ; Mechanism 2 enhances the algorithm’s exploration of the space, with the execution condition r a n d 2 0.5 ; Mechanism 3 strengthens the algorithm’s ability to escape local optima, with the execution condition r a n d 2 < 0.5 . The specific formulas for the exploration phase are as follows in Equations (1)–(6):
G X ( t + 1 ) = ( U B L B ) × r 1 + L B , r a n d 1 < p ( r 2 C ) × X r ( t ) + L × H , r a n d 2 0.5 a n d r a n d 1 p X ( t ) L × ( L × ( X ( t ) G X ( t ) ) ) + r 3 × ( X ( t ) G X r ( t ) ) , r a n d 2 < 0.5 a n d r a n d 1 p
C = F × 1 t M a x I t
F = cos ( 2 × r 4 ) + 1
L = C × l
H = Z × X ( t )
Z = [ C , C ]
In the formulas, X ( t ) and G X ( t + 1 ) represent the current position and the position in the next iteration of a gorilla, respectively, while X r ( t ) and G X r ( t ) are positions of randomly selected gorillas. t denotes the current iteration number. r 1 , r 2 , r 3 , r 4 , l and r a n d are all random numbers; r 1 , r 2 , r 3 , l , r a n d [ 0 , 1 ] , while r 4 [ 1 , 1 ] . U B and L B represent the upper and lower bounds of the variables, respectively. Z is a random number within the range [ C , C ] , where the parameter C has significant variation in the initial stages and gradually decreases later on. L indicates the leadership ability of the silverback gorilla, which might make incorrect decisions due to lack of experience in finding food or managing the group. At the end of the exploration phase, the fitness values of G X ( t ) and X ( t ) are calculated and compared. If the fitness value of G X ( t ) is smaller, the position X ( t ) is replaced by G X ( t ) .
Finally, during the exploitation phase, the algorithm employs two behavioral mechanisms: following the silverback gorilla and competing for mature females. The parameter C is used to switch between these mechanisms. If C W , the gorillas follow the silverback gorilla. The silverback, as the leader of the group, guides the gorillas to food sources and ensures the safety of the group. All gorillas in the group adhere to the decisions made by the silverback. If C < W , competition occurs. The silverback gorilla may age and die, allowing a blackback gorilla to potentially become the leader, or other male gorillas may challenge the silverback in combat to dominate the group.
When C W , the mechanism for following the silverback gorilla is described by the following equations:
G X ( t + 1 ) = L × M × ( X ( t ) X s i l v e r b a c k ) + X ( t )
M = 1 N i = 1 N G X i ( t ) g 1 g
g = 2 L
In the equations, X s i l v e r b a c k represents the position of the silverback gorilla (the optimal position). G X i ( t ) denotes the position of each candidate gorilla during the iteration. N indicates the total number of gorillas.
When C < W , the calculation formula for the competitive adult female mechanism is shown in Equation (10).
G X ( t ) = X s i l v e r b a c k ( X s i l v e r b a c k × Q X ( t ) × Q ) × A Q = 2 × r 5 1 A = β × E E = N 1 , r a n d 0.5 N 2 , r a n d < 0.5
where Q represents the intensity of gorilla competition; A denotes the coefficient of the competition degree; and E represents the impact of violence on the dimensions of the solution. r 5 and r a n d are random numbers within the interval [ 0 , 1 ] . When r a n d 0.5 , E is a random value from a normal distribution and within the problem dimensions; otherwise, E is a random value from a normal distribution.
At the end of the development phase, a population operation is conducted, which involves estimating the fitness values of all G X individuals. If G X ( t ) < X ( t ) , the G X ( t ) individuals replace the X ( t ) individuals. The best solution (the minimum fitness value) obtained from the entire population is considered the silverback gorilla.

2.2. Improved Artificial Gorilla Force Algorithm (OLGTO)

Based on the iterative optimization process of the GTO algorithm, it is known that the GTO algorithm has a good capability to avoid local optima, achieving high convergence accuracy and speed. However, in the early stages of iteration, the large number of formula parameters prevents timely information exchange among gorillas, adversely affecting the algorithm’s ability to escape local optima and impacting convergence speed and accuracy in later stages. To address these issues, the algorithm is improved in four aspects.
(1) We have adopted the Logistic chaotic map to initialize the gorilla population, enhancing the diversity of the population. This modification allows the algorithm to more uniformly cover the search space during the initial phase, resulting in higher convergence precision and reducing the risk of the population falling into local optima.
(2) An improvement in the weight factor W has been made to balance the algorithm’s global search and the ability to escape local optima. This adjustment allows the algorithm to more effectively avoid local optima during the convergence process, thereby enhancing the overall optimization performance, with improvements in both convergence speed and solution accuracy.
(3) The Fish Hawk Optimization Algorithm has been used to replace the second formula in the original GTO algorithm’s global exploration strategy in the first phase. This change aims to improve the global search capability of the algorithm. The modified algorithm demonstrates stronger global search capabilities in complex search spaces, enabling faster convergence to the optimal solution and avoiding the reduction in global search performance due to delayed information exchange among populations caused by excessive parameters.
(4) After the development phase, Levy flight has been applied to update the positions of the gorillas, which enhances the robustness of the algorithm. This application makes the algorithm more effective at escaping local optima in the later stages of iteration.

2.2.1. Logistic Chaotic Mapping Initializes Gorilla Population

Logistic chaotic mapping is a simple yet nonlinear mapping with complex dynamic behavior used to generate chaotic sequences. Compared to other chaotic mappings, it has the advantages of a relatively simple expression, strong adaptability, and the generated sequences tend to stabilize without infinitely increasing or decreasing [29].
Using Logistic chaotic mapping to initialize the population can enhance population diversity, thereby improving the optimization performance and global search capability of the algorithm. The Logistic chaotic mapping is mapped to the search space through Equations (11) and (12).
Y n + 1 = α Y n ( 1 Y n )
Z n = a + ( b a ) Y n
where α [ 0 , 4 ] is the control parameter of the system. When α is within the range [ 3.57 , 4 ] , the system enters a chaotic state [30], and in this study, α is set to 4. Y n is a random number within the interval [ 0 , 1 ] ; Y n represents the chaotic sequence generated by Equation (11); Z n denotes the chaotic sequence mapped to the search space; a and b are the lower and upper bounds of the search space, respectively.
Firstly, a chaotic sequence of length 30 (where the population size is 15 and the search space dimension is 2) is generated using Equation (11) and iterated accordingly. Next, the generated chaotic sequence is mapped to the search space with two-dimensional ranges of [ 100 , 2500 ] and [ 3 , 10 ] using Equation (12). Finally, the initial population is generated as described by Equation (13).
Y = y 1 , 1 y 1 , 2 y 2 , 1 y 2 , 2 y 15 , 1 y 15 , 2

2.2.2. Linearly Decreasing Weight Factors

Due to the GTO algorithm’s inability to effectively control the iteration step size during the iterative process, once an optimal solution is identified, other individuals quickly converge towards this optimal solution, leading to premature convergence to a local optimum and losing the opportunity for further exploration of the global optimum. However, by employing the linear decreasing strategy for the convergence factor as shown in Equation (14), the balance between the exploration and exploitation behaviors of the GTO algorithm can be improved [31]. This approach enables the algorithm to comprehensively explore the search space while also performing fine-tuned exploitation near high-quality solutions, thereby enhancing the algorithm’s adaptability and robustness across different problems.
W ( t ) = W max W max W min T max × t
where W max and W min represent the initial and final weight factors, with values of 1.5 and 0.4, respectively. t and T max denote the current iteration number and the maximum number of iterations. The convergence factor starts with a higher weight in the early stages to maintain strong global search capabilities and linearly decreases as the iteration progresses, thereby enhancing the algorithm’s local exploitation ability.

2.2.3. Integrated Osprey Algorithm Global Exploration Strategy

In the original GTO algorithm, the second formula in Equation (1) of the first-stage exploration strategy contains numerous parameters, which leads to a slower information exchange among gorillas and subsequently affects the global search speed and convergence performance of the algorithm. To address this issue, the more efficient Osprey Optimization Algorithm is used to replace the second formula in Equation (1) during the first-stage global exploration strategy [32]. This introduces randomness into the global search process, thereby enhancing the convergence performance of the GTO algorithm. The global exploration strategy of the Osprey Optimization Algorithm is described by Equations (15)–(17).
x i , j P 1 = x i , j + r i , j ( S F i , j I i , j x i , j )
x i , j P 1 = x i , j P 1 , l b j x i , j P 1 u b j l b j , x i , j P 1 < l b j u b j , x i , j P 1 > u b j
X i = X i P 1 , F i P 1 < F i X i , e l s e
where x i , j P 1 represents the new position information of the i-th osprey in the j-th dimension during the first stage; S F i , j denotes the selected fish for the i-th osprey in the j-th dimension; r i , j is a random number within the interval 0 , 1 ; I i , j is a random number from the set 1 , 2 ; and u b j and l b j denote the upper and lower bounds of the optimization search space, respectively. The updated formula for the gorilla position is as follows:
G X ( t + 1 ) = ( U B L B ) × r 1 + L B , r a n d 1 < p X r ( t ) + r 2 × ( X ( t ) I × X r ( t ) ) , r a n d 2 0.5 a n d r a n d 1 p X ( t ) L × ( L × ( X ( t ) G X ( t ) ) ) + r 3 × ( X ( t ) G X r ( t ) ) , r a n d 2 < 0.5 a n d r a n d 1 p

2.2.4. Levy Flight Strategy

In the exploitation phase of the GTO algorithm, the population leader (the silverback gorilla) is responsible for guiding the population towards food sources. At this stage, the accuracy of the silverback gorilla’s exploration is particularly crucial. Therefore, incorporating the Levy flight strategy [33] into the gorilla position update formula (Equation (7)) during the GTO algorithm’s exploitation phase helps the algorithm utilize known optimal solutions while exploring new potential solutions. The Levy flight strategy employs step lengths with heavy-tailed distributions, as described by Equation (19). This strategy allows for both small and large jumps, where large jumps enable the algorithm to explore new areas of the search space and escape local optima, while small jumps facilitate detailed search in known high-quality regions. This multi-scale search capability enhances the algorithm’s flexibility in handling complex problems.
R L = u v 1 λ
where λ is the power-law exponent of the Levy distribution, typically in a range of [ 1 , 2 ] and is set to 1.5 in this study. v ~ N ( 0 , 1 ) , u ~ N ( 0 , σ u 2 ) , and the standard deviation of u is calculated using Equation (20).
σ u = ( Γ ( 1 + β ) sin ( π β 2 ) Γ ( 1 + β 2 ) β 2 β 1 2 ) 1 β
where Γ is the gamma function.
The improved gorilla position update formula is given by Equation (21), as follows:
G X ( t + 1 ) = R L × ( X s i l v e r b a c k X ( t ) ) + X s i l v e r b a c k
This paper improves the GTO algorithm based on the above strategies, balancing the algorithm’s exploration and exploitation capabilities and enhancing the convergence speed. The improved algorithm is used to optimize the parameters of VMD and SVM. The specific process of the improved GTO algorithm is shown in Figure 1.

2.3. Improved Algorithm Testing

To validate the effectiveness and superiority of the improved GTO algorithm, it is tested on the single-peaked function F6 and the multi-peaked function F12 from the CEC2005 benchmark function set. The results are compared with the Whale Optimization Algorithm (WOA), Dung beetle optimization algorithm (DBO), Grey Wolf Optimizer (GWO), Northern Goshawk Optimization (NGO), Harris Hawks Optimization (HHO), the original GTO algorithm, and the improved GTO algorithm (OLGTO). Functions F6 and F12 are defined by Equations (22) and (23), respectively.
F 6 = Σ i = 1 n ( | x i + 0.5 | ) 2
F 12 ( x ) = π n { 10 sin 2 ( π y 1 ) + i = 1 n 1 ( y i 1 ) 2 [ 1 + 10 sin 2 ( π y i + 1 ) ] + ( y n 1 ) 2 } + i = 1 n u ( x i , 10 , 100 , 4 ) y i = 1 + x i + 1 4 u ( x i , a , k , m ) = k ( x i a ) , x i > a 0 , a < x i < a k ( x i a ) m , x i < a
where the dimensions and optimal minimum values for the F6 and F12 functions are 30 and 0, respectively. The search ranges for these functions are [−100, 100] and [−50, 50], respectively.
The population size and maximum number of iterations for the seven algorithms are uniformly set to 100 and 1000, respectively. In the OLGTO and GTO algorithms, the exploration probability p and control parameter Beta are set to 0.03 and 3, respectively. The weight factor W for GTO is set to 0.5. OLGTO employs a linear decreasing strategy and Levy flight strategy, with parameters W max = 1.5 , W min = 0.4 , and λ = 1.5 . The convergence results of the tested algorithms are shown in Figure 2.
In the entire search space, unimodal functions have only one global optimum and are primarily used to evaluate the algorithm’s exploitation capability, i.e., its ability to converge to the optimal solution. Unimodal functions are relatively simple, allowing for clear measurement of convergence speed and accuracy. On the other hand, multimodal functions contain multiple local optima and are highly suitable for assessing the algorithm’s exploration capability. Complex multimodal functions can test the algorithm’s ability to avoid local optima and find the global optimum while evaluating the balance between exploration and exploitation. From the comparative results in Figure 2, it can be observed that for unimodal functions, the OLGTO algorithm consistently converges at a faster rate and ultimately achieves higher convergence accuracy than the other six algorithms. For multimodal functions, except for DBO and OLGTO, the WOA, GTO, GWO, NGO, and HHO algorithms all experience a significant decline in convergence speed after 200 iterations, falling into local optima and resulting in lower final convergence accuracy. Although the OLGTO algorithm shows a noticeable decrease in convergence speed around the 400th iteration, indicating a tendency to fall into local optima, it quickly resumes its convergence and reaches the optimal solution by the 579th iteration. This demonstrates the algorithm’s strong balance between exploration and exploitation and its ability to escape local optima. In summary, the OLGTO algorithm exhibits strong adaptability compared to other algorithms, effectively avoiding local optima. This proves the effectiveness and superiority of the improved algorithm.

3. Fault Diagnosis Model Construction

3.1. VMD Fault Feature Extraction

Variational Mode Decomposition (VMD) is an improved method of the traditional Empirical Mode Decomposition (EMD), suitable for processing nonlinear and non-stationary signals. The principle of VMD is to decompose the original signal into a set of K intrinsic mode functions (IMFs), each with a central frequency. This decomposition is achieved through variational optimization. Additionally, VMD introduces Lagrange multipliers and a quadratic penalty factor α , and uses the Alternating Direction Method of Multipliers (ADMM) to solve the variational optimization problem to obtain the various modes and their corresponding central frequencies. By appropriately setting the total number of modes K and the value of the quadratic penalty factor α , VMD can overcome issues encountered by EMD when processing complex signals, such as mode mixing and endpoint effects, thereby enhancing VMD’s adaptability to different signal decompositions. VMD is used for fault feature extraction in vibration signals, employing a composite index: permutation entropy/mutual information entropy as the fitness function. Permutation entropy reflects the complexity and irregularity of time series, with lower permutation entropy indicating more regular and periodic decomposed sequences [34]. Mutual information entropy measures the correlation between different information contained in different modes of VMD, with higher mutual information entropy indicating stronger correlation [35]. Therefore, the composite index of permutation entropy/mutual information entropy can balance the regularity and correlation of the decomposition modes. The calculation formulas for permutation entropy, mutual information entropy, permutation entropy/mutual information entropy, and the fitness value are shown in Equations (24)–(27), respectively.
P E = p ( i ) log p i
M I = P X , Y log P X , Y / P X P Y
C I = P E M I
f i t n e s s = min 1 ~ K { C I }
where p ( i ) denotes the probability of the i-th permutation; X and Y are two random variables; P ( X , Y ) represents the joint probability distribution of X and Y ; and P ( X ) and P ( Y ) denote the marginal probability distributions of X and Y , respectively. When decomposing with VMD, the IMF component with the minimum fitness value is considered the current optimal IMF component. The mean, variance, peak value, kurtosis, root mean square, peak factor, impulse factor, waveform factor, and margin factor of the optimal IMF are then calculated, providing a 9-dimensional time-domain feature vector as the characteristic vector of the optimal IMF.

3.2. SVM Classification Model

Support Vector Machine is a flexible supervised learning model widely applied in fault diagnosis due to its high classification performance and ability to handle high-dimensional data. The core concept of SVM is to find an optimal separating hyperplane that divides samples of different classes. To achieve this, SVM solves the optimization problem as shown in Equation (28).
min ω , b , ξ 1 2 ω 2 + C i = 1 N ξ i , s . t . y i ( ω x i + b ) 1 ξ i ,   i = 1 , 2 , , N , ξ i 0 ,   i = 1 , 2 , , N
where ω is the normal vector of the hyperplane, b is the bias term, and ξ i are slack variables used to handle non-separable data. The regularization parameter C controls the trade-off between training error and margin. A larger C imposes a higher penalty on misclassification, potentially leading to overfitting, while a smaller C allows for more misclassification, which may improve generalization ability [36].
Furthermore, SVM utilizes a kernel function to map data from a low-dimensional space to a high-dimensional space, making the sample data linearly separable in the higher-dimensional space. A commonly used kernel function is the Radial Basis Function (RBF) kernel, as shown in Equation (29).
K ( x i , y i ) = exp σ x i y j 2
where σ is the kernel parameter that determines the influence range of a single sample. Setting appropriate values for C and σ can enhance the model’s generalization ability, reduce the complexity of high-dimensional calculations, and thus improve the accuracy of fault diagnosis [37].

3.3. OLGTO-VMD-OLGTO-SVM Fault Diagnosis Model

To obtain suitable values for the number of modes K in the Variational Mode Decomposition (VMD), the quadratic penalty factor α , the regularization parameter C , and the kernel parameter σ , and to effectively extract vibration signal features to improve the accuracy of fault diagnosis models, this paper proposes the use of the OLGTO algorithm. The OLGTO algorithm iteratively selects the optimal values for the critical parameters of the VMD and SVM models, establishing the OLGTO-VMD-OLGTO-SVM fault diagnosis model. The overall design process of the model is shown in Figure 3, with the specific steps as follows:
Step 1: Initialize the settings for the gorilla population size, the maximum number of iterations, and other relevant parameters.
Step 2: Distribute the gorilla population randomly and uniformly in the search space using Logistic chaotic mapping. At this stage, the position vectors of the gorilla individuals represent the initial values of α , K , C , and σ .
Step 3: Use the composite metric of permutation entropy/mutual information entropy as the fitness function to calculate the initial fitness values of the gorilla individuals and sort them in ascending order to determine the optimal fitness value and position.
Step 4: Update the optimal fitness value and position according to Formulas (18), (21), and (10) during the exploration and exploitation phases.
Step 5: Check if the current number of iterations exceeds the maximum number of iterations. If it does, the optimization process ends, and the optimal fitness value and position of the gorilla individuals, representing α , K , C , and σ , are outputted. Otherwise, return to Step 2 and continue the optimization process.
Step 6: Substitute the optimal parameters α and K into the VMD to extract features from the vibration signals. Divide the dataset into training and testing sets.
Step 7: Substitute the optimal parameters C and σ into the SVM to construct the fault diagnosis model and input the dataset to achieve fault diagnosis.

4. Fault Diagnosis Experiment Verification and Analysis

4.1. Data Preparation

In practical applications, it is necessary to consider the potential insufficiency of fault sample data. Therefore, fault diagnosis experiments utilize both the small-sample vibration data from the drive end of motor bearings from Case Western Reserve University and the large-sample vibration data from three-phase asynchronous motors to, respectively, verify the feasibility and superiority of the OLGTO-VMD-OLGTO-SVM fault diagnosis model.
The motor bearing fault diagnosis experiment at Case Western Reserve University (CWRU) involves a 2-horsepower (1.5 kW) motor, a torque sensor, a power meter, and an electronic controller. The test bearings, SKF6205, were subjected to single-point faults using electrical discharge machining, creating defect diameters of 0.1778 mm, 0.3556 mm, and 0.5334 mm. An accelerometer placed above the bearing housing on the drive end collected vibration acceleration signals at a sampling rate of 12 kHz. The experiment analyzed ten states: one normal operating condition (labeled as 1) at 1797 r/min and nine fault conditions (three each for inner race defects, rolling element defects, and outer race defects) corresponding to the different defect diameters (labeled as 2, 5, and 8 for inner race; 3, 6, and 9 for rolling element; and 4, 7, and 10 for outer race). Time series vibration data were segmented using a sliding window of 1000, with each sample containing 2048 data points. Each state type included 120 samples, resulting in a dataset of 1200 samples. Due to space limitations, only the first 1 s samples for the normal state and each defect type with a 0.1778 mm defect diameter are shown in Figure 4.
The data collection for the three-phase induction motor is illustrated in Figure 5. The experiment involves a 1-horsepower (750 W) three-phase induction motor, a piezoelectric accelerometer, and a single-stage gearbox. The motor model is YXVF0M-2, with an adjustable speed range of 0–5000 r/min. The accelerometer model is CK-EU-50B, with a measurement range of 50 g and a sampling frequency of 5 kHz. The gearbox is a single-stage gearbox with a module of 1.5 and a gear ratio of 75:25. This study examines five fault conditions: normal operating state (labeled as 1) at a speed of 1500 r/min, motor rotor imbalance fault (labeled as 2), motor rotor misalignment fault (labeled as 3), motor bearing fault (labeled as 4), motor warped rotor fault (labeled as 5), and motor rotor broken bar fault (labeled as 6). Data for these six states were processed using the same sliding window method, with each sample containing 2048 data points. A total of 2000 sample sets were collected for each state, resulting in a dataset of 12,000 samples. The sampled signals for each state during the first second are shown in Figure 6.

4.2. Motor Bearing Fault Diagnosis Experiment

To validate the feasibility of the OLGTO-VMD-OLGTO-SVM method, four sets of experiments were conducted. (1) The motor bearing vibration signals were directly input into the SVM model for diagnosis. (2) The motor bearing vibration signals were directly input into the OLGTO-SVM model for diagnosis. (3) The motor bearing vibration signals were input into the SVM model for diagnosis after feature extraction using OLGTO-VMD. (4) The motor bearing vibration signals were input into the OLGTO-SVM model for diagnosis after feature extraction using OLGTO-VMD.
For experiments (1) and (2), the raw data were processed into a multi-dimensional sequence matrix of 1200 × 2048 and input as 1 × 2048 one-dimensional sequences into the SVM and OLGTO-SVM models, respectively. In experiment (1), the parameters C and σ were set to 3 and 1.5, respectively.
For experiments (3) and (4), the OLGTO algorithm was first used to optimize the parameters α and K in the VMD for feature extraction from the bearing vibration signals. The initial parameter settings are shown in Table 1. For each bearing state, 120 feature vectors were extracted (10 states in total, resulting in 1200 sets). The 1200 × 2048 data samples were reduced to 1200 × 9 feature vector samples after feature extraction. The OLGTO-VMD results are shown in Table 2.
The iterative convergence curves are shown in Figure 7. Due to space limitations, only the convergence curves for the first four bearing states are displayed.
As shown in Table 2 and Figure 7, the total number of mode decompositions for each vibration signal is not less than 5, but the optimal IMF components are consistently IMF1 and IMF2. This indicates that the main characteristic components are fully extracted in the first few mode components after OLGTO-VMD. Additionally, for all 10 states, the optimal fitness values are achieved before the end of the iterations, demonstrating the algorithm’s effective optimization performance.
Secondly, for each bearing state, the first 80 sets were used as the training set, and the last 40 sets were used as the testing set. The parameters C , σ , and data input method for experiment (3) were the same as those in experiment (1). For experiment (4), the training set samples were first input into the OLGTO-SVM model for training. The initial parameter settings for the OLGTO-SVM model are shown in Table 3. Through iteration, the optimal positions of the gorillas, representing the optimal parameters C and σ of the SVM, were obtained. The ultimate goal of fault diagnosis is to achieve accurate classification of different fault types. Therefore, the diagnostic accuracy as shown in Equation (30) was used as the fitness function. The results of OLGTO optimizing the SVM parameters are shown in Table 4, and the fitness curves are shown in Figure 8.
A c c u r a c y = N c o r r e c t N t o t a l × 100 %
As shown in Figure 8, the fitness values for experiments (2) and (4) reach their optimum at the 16th and 3rd iterations, respectively, and remain stable thereafter. This indicates that the algorithm performs well in optimization. The optimal parameters corresponding to the best fitness obtained through the iterative optimization of the OLGTO algorithm were used as the optimal parameters for the SVM model to test the sample data. The results of optimizing the SVM parameters are shown in Table 4.
In summary, the fault diagnosis results for experiments (1) through (4) are shown in Table 5.
The confusion matrices for the fault diagnosis results are shown in Figure 9.
From the experimental results above, it can be seen that the SVM model performs poorly in identifying vibration signals without feature extraction. Even when optimized using OLGTO, there is no significant improvement in the final results. However, after feature extraction using the OLGTO-VMD method, the SVM achieves high diagnostic accuracy. Specifically, the OLGTO-VMD-OLGTO-SVM model achieves a diagnostic accuracy of 99.5%, which is a five-percentage-point improvement over the 94.5% accuracy of the OLGTO-VMD-SVM model. This validates the feasibility of applying the OLGTO-VMD-OLGTO-SVM method in fault diagnosis.

4.3. Fault Diagnosis Experiment of Three-Phase Asynchronous Motor

To validate the superiority of the OLGTO-VMD-OLGTO-SVM method, this method was compared with five other methods: OLGTO-VMD-CNN-BiLSTM (Convolutional Neural Network-Bidirectional Long Short-Term Memory), OLGTO-VMD-KNN (K-Nearest Neighbors), OLGTO-VMD-RF (Random Forest), OLGTO-VMD-BP (Back Propagation Neural Network), and OLGTO-VMD-ELM (Extreme Learning Machine). These six methods were applied to the fault diagnosis of an actual three-phase asynchronous motor, and their performance was compared. The experimental process was as follows: First, the vibration signals were processed using OLGTO-VMD to extract signal features. Then, the first 1200 sets of each motor state were used as the training set, and the last 800 sets were used as the testing set. These sets were input into OLGTO-SVM, CNN-BiLSTM, KNN, RF, BP, and ELM for training and diagnosis, respectively. The parameter settings for the OLGTO-VMD feature extraction were the same as those in Table 1. For each state of the three-phase asynchronous motor, 2000 feature vectors were extracted (six states in total, resulting in 12,000 sets). The 12,000 × 2048 sample data were reduced to 12,000 × 9 feature vector samples after feature extraction. The convergence curves of the parameter optimization iterations are shown in Figure 10, and the decomposition results of the motor vibration signals are shown in Table 6.
As shown by the decomposition results in Table 6, the total number of mode decompositions K for the three-phase asynchronous motor vibration signals is not less than seven, but the optimal IMF components are consistently IMF1 and IMF2. This indicates that the primary characteristic information of the vibration signals is concentrated in the first few IMF components after OLGTO-VMD. Furthermore, as observed from the convergence curves in Figure 10, the optimal fitness value for each state is achieved before the iterations are completed, demonstrating the algorithm’s effective optimization performance.
After feature extraction, multi-dimensional feature vector samples were obtained. These samples were input as one-dimensional sequences into OLGTO-SVM, CNN-BiLSTM, KNN, RF, BP, and ELM for training and diagnosis. The initial parameter settings for the OLGTO-SVM model were the same as those in Table 3. Using diagnostic accuracy as the fitness function, the parameter optimization results are shown in Figure 10.
As shown in Figure 11, it can be observed that the fitness value reaches its optimum for the first time during the fifth iteration and remains stable thereafter with no fluctuations, indicating an overall stable trend. The optimal parameters C and σ obtained from the iterative optimization process of the OLGTO algorithm are then input into the SVM model to test the samples. The optimization results of the SVM model with the optimal parameters are shown in Table 7.
Additionally, the CNN-BiLSTM model was trained using the Adam optimizer with a maximum of 100 training epochs, an initial learning rate of 0.01, and a learning rate adjustment factor of 0.01 starting after 70 epochs. The regularization parameter was set to 0.001. The KNN model was constructed using the fitcknn function, with the number of nearest neighbors set to 1 and the distance metric type set to Euclidean. The RF model was constructed using the TreeBagger function with a minimum leaf size of 1, no out-of-bag predictions, no surrogate splits, and the number of decision trees set to 10. For the BP model, the optimal number of hidden layer nodes was set to 15. The hidden layer used the tansig activation function, the output layer used the purelin activation function, and the Levenberg–Marquardt algorithm was used as the training function. The maximum number of training iterations was set to 1000, with a learning rate of 0.1, a training goal error of 0.00001, a momentum factor of 0.01, a minimum gradient of 1 × 10−6, and a maximum validation failure count of 6. The ELM model was used in regression mode with a sine activation function and 100 hidden layer nodes. The diagnostic results of the six methods are shown in Table 8.
The confusion matrices for the diagnostic results of the six methods are shown in Figure 12.
From the comparison results above, it can be seen that despite the large data sample size, the OLGTO-VMD-OLGTO-SVM method still achieves a diagnostic accuracy of 98.6458%. This represents a significant improvement compared to the accuracies of 92.6667%, 92.8958%, and 95.0833% achieved by the OLGTO-VMD-KNN, OLGTO-VMD-CNN-BiLSTM, and OLGTO-VMD-ELM methods, respectively. While the improvements over OLGTO-VMD-BP and OLGTO-VMD-RF, which achieved 98.125% and 98.4375% accuracy, respectively, are smaller, it is important to note that a 0.3~0.5% increase in accuracy still represents a considerable number of data samples given the large sample size. Therefore, the method proposed in this paper is more reliable for practical applications, demonstrating the superiority of the OLGTO-VMD-OLGTO-SVM method.

5. Conclusions

(1) To address the deficiency of the Artificial Gorilla Troops Optimization algorithm in the early iteration stages, where the numerous formula parameters hinder timely information exchange between gorillas, reducing the algorithm’s global search capability and ability to escape local optima, and thus degrading convergence accuracy, improvements were made using Logistic chaotic mapping, a linear decreasing weight factor strategy, the global exploration strategy of the Osprey Optimization Algorithm, and the Levy flight strategy. Tests on unimodal and multimodal functions showed that the improved GTO algorithm possesses excellent global search capabilities and effectively avoids local optima.
(2) The selection of parameters α and K in VMD determines the effectiveness of vibration signal decomposition and feature extraction, while the parameters C and σ in SVM determine the final diagnostic accuracy. The improved GTO algorithm was used to optimize the parameters of both VMD and SVM, overcoming the difficulty of manually selecting optimal parameters and enhancing the generalization ability and adaptability of the fault diagnosis model.
(3) Fault diagnosis experiments on motor bearings and three-phase asynchronous motors validated the feasibility and superiority of the OLGTO-VMD-OLGTO-SVM method. This method can accurately identify fault types for both large and small sample data.

Author Contributions

H.Z.: Methodology, Validation, Formal analysis, Writing—original draft. T.F.: Writing—review and editing. L.M.: Investigation, Conceptualization, Project administration. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The National Natural Science Foundation Project (51007002); Anhui Industrial Internet Intelligent application and security engineering laboratory open fund (IASII21-05); Research on data synthesis and image detection methods for appearance defects of power equipment (QZ202108); Anhui University of Technology youth fund (QZ202109); Research on ship target recognition technology based on artificial intelligence (QZ202109).

Data Availability Statement

Datasets are available on request from the authors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Liu, C.Y.; Zou, W.; Hu, Z.; Li, H.; Sui, X.; Ma, X.; Yang, F.; Guo, N. Bearing Health State Detection Based on Informer and CNN + Swin Transformer. Machines 2024, 12, 456. [Google Scholar] [CrossRef]
  2. Gao, Z.W.; Cecati, C.; Ding, S.X. A survey of fault diagnosis and fault-tolerant techniques—Part I: Fault diagnosis with model-based and signal-based approaches. IEEE Trans. Ind. Electron. 2015, 62, 3757–3767. [Google Scholar] [CrossRef]
  3. Zheng, J.D.; Su, M.; Ying, W.; Tong, J.; Pan, Z. Improved uniform phase empirical mode decomposition and its application in machinery fault diagnosis. Measurement 2021, 179, 109425. [Google Scholar] [CrossRef]
  4. Yue, S.Y.; Wang, Y.; Wei, L.; Zhang, Z.; Wang, H. The joint empirical mode decomposition-local mean decomposition method and its application to time series of compressor stall process. Aerosp. Sci. Technol. 2020, 105, 105969. [Google Scholar] [CrossRef]
  5. Li, Y.-X.; Jiao, S.-B.; Gao, X. A novel signal feature extraction technology based on empirical wavelet transform and reverse dispersion entropy. Def. Technol. 2021, 17, 1625–1635. [Google Scholar] [CrossRef]
  6. Liu, Q.Q.; Yang, J.T.; Zhang, K. An improved empirical wavelet transform and sensitive components selecting method for bearing fault. Measurement 2022, 187, 110348. [Google Scholar] [CrossRef]
  7. He, X.Z.; Zhou, X.; Yu, W.; Hou, Y.; Mechefske, C.K. Adaptive variational mode decomposition and its application to multi-fault detection using mechanical vibration signals. ISA Trans. 2021, 111, 360–375. [Google Scholar] [CrossRef]
  8. Gu, J.; Peng, Y.; Lu, H.; Chang, X.; Cao, S.; Chen, G.; Cao, B. An optimized variational mode decomposition method and its application in vibration signal analysis of bearings. Struct. Health Monit. 2022, 21, 2386–2407. [Google Scholar] [CrossRef]
  9. Wang, J.Y.; Li, J.; Wang, H.; Guo, L. Composite fault diagnosis of gearbox based on empirical mode decomposition and improved variational mode decomposition. J. Low Freq. Noise Vib. Act. Control 2021, 40, 332–346. [Google Scholar] [CrossRef]
  10. Gu, J.; Peng, Y.X. An improved complementary ensemble empirical mode decomposition method and its application in rolling bearing fault diagnosis. Digit. Signal Process. 2021, 113, 103050. [Google Scholar] [CrossRef]
  11. Goyal, D.; Choudhary, A.; Sandhu, J.K.; Srivastava, P.; Saxena, K.K. An intelligent self-adaptive bearing fault diagnosis approach based on improved local mean decomposition. Int. J. Interact. Des. Manuf. 2022, 1–11. [Google Scholar] [CrossRef]
  12. Wang, J.Y.; Li, J.; Wang, H.E.J. Research on gearbox composite fault diagnosis based on improved local mean decomposition. Int. J. Dyn. Control 2021, 9, 1411–1422. [Google Scholar] [CrossRef]
  13. Qin, C.R.; Wang, D.; Xu, Z.; Tang, G. Improved empirical wavelet transform for compound weak bearing fault diagnosis with acoustic signals. Appl. Sci. 2020, 10, 682. [Google Scholar] [CrossRef]
  14. Kedadouche, M.; Thomas, M.; Tahan, A.J.M.S. A comparative study between Empirical Wavelet Transforms and Empirical Mode Decomposition Methods: Application to bearing defect diagnosis. Mech. Syst. Signal Process. 2016, 81, 88–107. [Google Scholar] [CrossRef]
  15. Chen, X.J.; Yang, Y.; Cui, Z.; Shen, J. Vibration fault diagnosis of wind turbines based on variational mode decomposition and energy entropy. Energy 2019, 174, 1100–1109. [Google Scholar] [CrossRef]
  16. Wang, Z.J.; He, G.; Du, W.; Zhou, J.; Han, X.; Wang, J.; He, H.; Guo, X.; Wang, J.; Kou, Y. Application of parameter optimized variational mode decomposition method in fault diagnosis of gearbox. IEEE Access 2019, 7, 44871–44882. [Google Scholar] [CrossRef]
  17. Kumar, A.; Zhou, Y.Q.; Xiang, J.W. Optimization of VMD using kernel-based mutual information for the extraction of weak features to detect bearing defects. Measurement 2021, 168, 108402. [Google Scholar] [CrossRef]
  18. Yao, G.; Wang, Y.; Benbouzid, M.; Ait-Ahmed, M.; Yao, G.; Wang, Y.; Benbouzid, M.; Ait-Ahmed, M. A hybrid gearbox fault diagnosis method based on GWO-VMD and DE-KELM. Appl. Sci. 2021, 11, 4996. [Google Scholar] [CrossRef]
  19. Liu, Z.C.; Li, S.; Wang, R.; Jia, X. Research on Fault Feature Extraction Method of Rolling Bearing Based on SSA–VMD–MCKD. Electronics 2022, 11, 3404. [Google Scholar] [CrossRef]
  20. Lei, Y.G.; Yang, B.; Jiang, X.; Jia, F.; Li, N.; Nandi, A.K. Applications of machine learning to machine fault diagnosis: A review and roadmap. Mech. Syst. Signal Process. 2020, 138, 106587. [Google Scholar] [CrossRef]
  21. Lu, J.T.; Qian, W.; Li, S.; Cui, R. Enhanced K-nearest neighbor for intelligent fault diagnosis of rotating machinery. Appl. Sci. 2021, 11, 919. [Google Scholar] [CrossRef]
  22. Hosseinpour-Zarnaq, M.; Omid, M.; Biabani-Aghdam, E. Fault diagnosis of tractor auxiliary gearbox using vibration analysis and random forest classifier. Inf. Process. Agric. 2022, 9, 60–67. [Google Scholar] [CrossRef]
  23. Wang, Z.Y.; Yao, L.G.; Cai, Y.W. Rolling bearing fault diagnosis using generalized refined composite multiscale sample entropy and optimized support vector machine. Measurement 2020, 156, 107574. [Google Scholar] [CrossRef]
  24. Xu, L.Y.; Zhao, G.; Zhao, S.; Wu, Y.; Chen, X. Fault Diagnosis Method for Tractor Transmission System Based on Improved Convolutional Neural Network–Bidirectional Long Short-Term Memory. Machines 2024, 12, 492. [Google Scholar] [CrossRef]
  25. Toma, R.N.; Prosvirin, A.E.; Kim, J.-M. Bearing fault diagnosis of induction motors using a genetic algorithm and machine learning classifiers. Sensors 2020, 20, 1884. [Google Scholar] [CrossRef]
  26. Shao, S.Y.; Yan, R.; Lu, Y.; Wang, P.; Gao, R.X. DCNN-based multi-signal induction motor fault diagnosis. IEEE Trans. Instrum. Meas. 2019, 69, 2658–2669. [Google Scholar] [CrossRef]
  27. Guo, Z.J.; Liu, M.; Wang, Y.; Qin, H. A new fault diagnosis classifier for rolling bearing united multi-scale permutation entropy optimize VMD and cuckoo search SVM. IEEE Access 2020, 8, 153610–153629. [Google Scholar] [CrossRef]
  28. Abdollahzadeh, B.; Gharehchopogh, F.S.; Mirjalili, S. Artificial gorilla troops optimizer: A new nature-inspired metaheuristic algorithm for global optimization problems. Int. J. Intell. Syst. 2021, 36, 5887–5958. [Google Scholar] [CrossRef]
  29. Fadhil, M.S.; Farhan, A.K.; Fadhil, M.N. Designing substitution box based on the 1D logistic map chaotic system. In IOP Conference Series: Materials Science and Engineering; IOP Publishing: Bristol, UK, 2021; Volume 1076. [Google Scholar]
  30. Alawida, M. Enhancing logistic chaotic map for improved cryptographic security in random number generation. J. Inf. Secur. Appl. 2024, 80, 103685. [Google Scholar] [CrossRef]
  31. Choudhary, S.; Sugumaran, S.; Belazi, A.; El-Latif, A.A.A. Linearly decreasing inertia weight PSO and improved weight factor-based clustering algorithm for wireless sensor networks. J. Ambient. Intell. Humaniz. Comput. 2023, 14, 6661–6679. [Google Scholar] [CrossRef]
  32. Dehghani, M.; Trojovský, P. Osprey optimization algorithm: A new bio-inspired metaheuristic algorithm for solving engineering optimization problems. Front. Mech. Eng. 2023, 8, 1126450. [Google Scholar] [CrossRef]
  33. Kaidi, W.; Khishe, M.; Mohammadi, M. Dynamic levy flight chimp optimization. Knowl. Based Syst. 2022, 235, 107625. [Google Scholar] [CrossRef]
  34. Ribeiro, M.; Henriques, T.; Castro, L.; Souto, A.; Antunes, L.; Costa-Santos, C.; Teixeira, A. The entropy universe. Entropy 2021, 23, 222. [Google Scholar] [CrossRef] [PubMed]
  35. He, J.; Qu, L.; Wang, P.; Li, Z. An oscillatory particle swarm optimization feature selection algorithm for hybrid data based on mutual information entropy. Appl. Soft Comput. 2024, 152, 111261. [Google Scholar] [CrossRef]
  36. Cortes, C. Support-Vector Networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar] [CrossRef]
  37. Scholkopf, B.; Alexander, J.S. Learning with Kernels: Support Vector Machines, Regularization, Optimization, and Beyond; MIT Press: Cambridge, MA, USA, 2018. [Google Scholar]
Figure 1. Flowchart of improved Artificial Gorilla Troop Optimization algorithm.
Figure 1. Flowchart of improved Artificial Gorilla Troop Optimization algorithm.
Machines 12 00637 g001
Figure 2. Test algorithm convergence results.
Figure 2. Test algorithm convergence results.
Machines 12 00637 g002
Figure 3. Overall fault diagnosis process.
Figure 3. Overall fault diagnosis process.
Machines 12 00637 g003
Figure 4. Bearing vibration signal waveform.
Figure 4. Bearing vibration signal waveform.
Machines 12 00637 g004
Figure 5. Three-phase asynchronous motor data acquisition.
Figure 5. Three-phase asynchronous motor data acquisition.
Machines 12 00637 g005
Figure 6. Motor vibration signal waveform.
Figure 6. Motor vibration signal waveform.
Machines 12 00637 g006aMachines 12 00637 g006b
Figure 7. Convergence curve of OLGTO-VMD (bearing).
Figure 7. Convergence curve of OLGTO-VMD (bearing).
Machines 12 00637 g007aMachines 12 00637 g007b
Figure 8. Fitness curve of OLGTO-SVM (bearing).
Figure 8. Fitness curve of OLGTO-SVM (bearing).
Machines 12 00637 g008
Figure 9. Troubleshooting results of experiments (1)–(4).
Figure 9. Troubleshooting results of experiments (1)–(4).
Machines 12 00637 g009aMachines 12 00637 g009b
Figure 10. Convergence curve of OLGTO-VMD (Electric motors).
Figure 10. Convergence curve of OLGTO-VMD (Electric motors).
Machines 12 00637 g010
Figure 11. Fitness curve of OLGTO-SVM (electric motors).
Figure 11. Fitness curve of OLGTO-SVM (electric motors).
Machines 12 00637 g011
Figure 12. Results of 6 methods of diagnosis (Electric motors).
Figure 12. Results of 6 methods of diagnosis (Electric motors).
Machines 12 00637 g012aMachines 12 00637 g012b
Table 1. OLGTO-VMD parameter settings.
Table 1. OLGTO-VMD parameter settings.
Gorilla
Population Size
Maximum IterationsOptimization
Parameters
Search Range
1520 α [ 100 , 2500 ]
K [ 3 , 10 ]
Table 2. OLGTO-VMD result (bearing).
Table 2. OLGTO-VMD result (bearing).
Bearing Status α K Fitness ValueOptimal IMF Component
1110100.53519IMF1
22104100.72607IMF2
311090.70307IMF1
4249150.72564IMF2
5198970.75217IMF1
61980100.65067IMF1
72271100.66239IMF1
8208290.68290IMF1
92372100.68985IMF1
102491100.79858IMF1
Table 3. OLGTO-SVM parameter settings.
Table 3. OLGTO-SVM parameter settings.
Gorilla
Population Size
Maximum
Iterations
Optimization
Parameters
Search Range
1030 C [ 10 3 , 10 3 ]
σ [ 10 8 , 10 3 ]
Table 4. OLGTO-SVM parameter optimization results (bearing).
Table 4. OLGTO-SVM parameter optimization results (bearing).
Algorithm ModelRegularization
Parameter C
Kernel
Parameter σ
Experiment
OLGTO-SVM471.6035747.5808(2)
20.262742.4093(4)
Table 5. Comparison of experimental results of motor bearing fault diagnosis.
Table 5. Comparison of experimental results of motor bearing fault diagnosis.
MethodsAccuracy
SVM29.5%
OLGTO-SVM29.75%
OLGTO-VMD-SVM94.5%
OLGTO-VMD-OLGTO-SVM99.5%
Table 6. OLGTO-VMD result (Electric motors).
Table 6. OLGTO-VMD result (Electric motors).
Motor States α K Fitness ValueOptimal IMF Component
12054100.65074IMF1
2117670.67622IMF1
32319100.65835IMF1
414480.93037IMF1
5955100.61136IMF1
6154570.661504IMF1
Table 7. OLGTO-SVM parameter optimization results (electric motors).
Table 7. OLGTO-SVM parameter optimization results (electric motors).
Algorithm ModelAlgorithm Model C Algorithm Model σ
OLGTO-SVM81.4547588.5758
Table 8. Comparison of experimental results of motor fault diagnosis.
Table 8. Comparison of experimental results of motor fault diagnosis.
MethodsAccuracy
OLGTO-VMD-KNN92.6667%
OLGTO-VMD-CNN-BiLSTM92.8958%
OLGTO-VMD-ELM95.0833%
OLGTO-VMD-BP98.125%
OLGTO-VMD-RF98.4375%
OLGTO-VMD-OLGTO-SVM98.6458%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fang, T.; Ma, L.; Zhang, H. Research on Fault Diagnosis Method with Adaptive Artificial Gorilla Troops Optimization Optimized Variational Mode Decomposition and Support Vector Machine Parameters. Machines 2024, 12, 637. https://doi.org/10.3390/machines12090637

AMA Style

Fang T, Ma L, Zhang H. Research on Fault Diagnosis Method with Adaptive Artificial Gorilla Troops Optimization Optimized Variational Mode Decomposition and Support Vector Machine Parameters. Machines. 2024; 12(9):637. https://doi.org/10.3390/machines12090637

Chicago/Turabian Style

Fang, Ting, Long Ma, and Hongkai Zhang. 2024. "Research on Fault Diagnosis Method with Adaptive Artificial Gorilla Troops Optimization Optimized Variational Mode Decomposition and Support Vector Machine Parameters" Machines 12, no. 9: 637. https://doi.org/10.3390/machines12090637

APA Style

Fang, T., Ma, L., & Zhang, H. (2024). Research on Fault Diagnosis Method with Adaptive Artificial Gorilla Troops Optimization Optimized Variational Mode Decomposition and Support Vector Machine Parameters. Machines, 12(9), 637. https://doi.org/10.3390/machines12090637

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