Next Article in Journal
Enhancing Energy Power Quality in Low-Voltage Networks Integrating Renewable Energy Generation: A Case Study in a Microgrid Laboratory
Previous Article in Journal
Adsorption Air-Conditioning for Automotive Applications: A Critical Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Intelligent Backstepping Control of Permanent Magnet-Assisted Synchronous Reluctance Motor Position Servo Drive with Recurrent Wavelet Fuzzy Neural Network

1
Department of Electrical Engineering, National Central University, Taoyuan 32001, Taiwan
2
Department of Electrical Engineering, National Taipei University of Technology, Taipei 10608, Taiwan
*
Author to whom correspondence should be addressed.
Energies 2023, 16(14), 5389; https://doi.org/10.3390/en16145389
Submission received: 14 June 2023 / Revised: 10 July 2023 / Accepted: 13 July 2023 / Published: 14 July 2023
(This article belongs to the Topic Advanced Electrical Machines and Drives Technologies)

Abstract

:
An intelligent servo drive system for a permanent magnet-assisted synchronous reluctance motor (PMASynRM) that can adapt to the control requirements considering the motor’s nonlinear and time-varying natures is developed in this study. A recurrent wavelet fuzzy neural network (RWFNN) with intelligent backstepping control is proposed to achieve this. In this study, first, a maximum torque per ampere (MTPA) controlled PMASynRM servo drive is introduced. A lookup table (LUT) is created, which is based on finite element analysis (FEA) results by using ANSYS Maxwell-2D dynamic model to determine the current angle command of the MTPA. Next, a backstepping control (BSC) system is created to accurately follow the desired position in the PMASynRM servo drive system while maintaining robust control characteristics. However, designing an efficient BSC for practical applications becomes challenging due to the lack of prior uncertainty information. To overcome this challenge, this study introduces an RWFNN as an approximation for the BSC, aiming to alleviate the limitations of the traditional BSC approach. An enhanced adaptive compensator is also incorporated into the RWFNN to handle potential approximation errors effectively. In addition, to ensure the stability of the RWFNN, the Lyapunov stability method is employed to develop online learning algorithms for the RWFNN and to guarantee its asymptotic stability. The proposed intelligent backstepping control with recurrent wavelet fuzzy neural network (IBSCRWFNN) demonstrates remarkable effectiveness and robustness in controlling the PMASynRM servo drive, as evidenced by the experimental results.

1. Introduction

The permanent magnet synchronous motor (PMSM) is an electric motor type that has found extensive applications in diverse fields, including electric vehicles, industrial automation, robotics, and aerospace plane [1,2]. Compared to conventional induction motors, PMSM exhibits high efficiency exceeding 90%, resulting in substantial energy savings. It also possesses a high magnetic field strength and low internal resistance, enabling PMSM to achieve a higher power density and greater output power within a relatively smaller volume. Moreover, PMSM offers excellent torque-speed performance and sensitive current control ability, making it ideal for precision applications such as automatic control systems [3,4,5]. Despite its advantages, rare-earth elements such as NdFeB can be costly and raise concerns regarding the supply chain monopoly and trade wars. Additionally, high-performance applications using PMSM face challenges such as operating at the region of flux-weakening control with large direct-axis current and the potential for uncontrolled generator mode due to the flux linkages generated by permanent magnets (PM). These drawbacks limit the potential benefits of PMSM and should be taken into account when considering their use in specific applications. Hence, there is a need to decrease the utilization of rare-earth PMs.
Due to their rapid dynamic response, wide speed range, affordability, and high efficiency, the synchronous reluctance motor (SynRM) has gained significant popularity and is now utilized in numerous applications. Unlike induction motors and PMSM, SynRM offers inherent advantages by eliminating the need for permanent magnetic materials, cages, and excitation windings. This not only enhances its robustness but also reduces its overall cost, making it an attractive alternative [6,7,8,9]. However, SynRM has limited overload capability compared to other types of motors, meaning that it may not be suitable for applications that require high levels of torque for short periods of time. Based on the aforementioned reasons, a relatively new machine, known as the permanent magnet-assisted synchronous reluctance motor (PMASynRM), has been developed to address the challenges associated with the scarcity of rare-earth PMs. One approach to mitigating these difficulties involves reducing the amount of rare-earth PMs used in the rotor or replacing them with ferrite magnets [10,11,12,13]. Maximum torque per ampere (MTPA) has been widely utilized as a control strategy in PMSMs and SynRMs control to maximize output torque [14,15,16,17,18,19]. In [17,18], the MTPA control strategy has been implemented in the control of IPMSMs and SynRMs to determine the ideal current angle for maximizing the output torque based on a given stator current. Moreover, in [19], an MTPA control with nonlinear simultaneous equations was derived from the Lagrange multiplier method, which could be solved by numerical algorithms. The MTPA control technology enables the production of equivalent torque with minimal current. It achieves this by identifying the ideal current angle that maximizes output torque at a given stator current while also minimizing copper loss during the process. However, PMASynRM has inherent drawbacks, such as nonlinear and time-varying control characteristics, which make achieving high-performance servo applications and the traditional MTPA quite challenging [20].
The backstepping control (BSC), despite its advantages in providing a recursive and systematic design methodology for nonlinear feedback control, may encounter undesired chattering phenomena due to the presence of a sign function. Several approaches, such as adaptive control [21] and intelligent control [22], have been proposed for integration with BSC to overcome this issue and enhance control performance. In [21], a PMSM drive system was targeted, and an adaptive backstepping (ABS) control approach was introduced as a solution. The purpose of this method is to achieve precise tracking responses by utilizing the robustness properties of the ABS control. In [22], to overcome the limitations posed by the nonlinear and time-varying control natures of a SynRM, a robust position controller was devised for a SynRM servo drive system. This was achieved through the introduction of an intelligent BSC approach, employing a recurrent feature selection fuzzy neural network. In addition, to improve control performance and facilitate model-free controller design, intelligent control methods, including fuzzy mechanisms, NNs, and FNNs, have been widely used as universal approximators in various studies. Among these, the FNN has been particularly popular due to its combined advantages of neural networks and fuzzy logic and has been applied in various fields, such as photovoltaic systems, robotics, and motor control [23,24,25]. In fact, the FNN’s robustness and convenience have garnered attention for controlling permanent magnet linear synchronous motors [25]. Moreover, the recurrent neural network (RNN) is capable of mapping and storing temporal information dynamically [26,27], utilizing time delays from earlier states, and approximate information can be obtained from internal feedback states. This makes RFNN a better option for dynamic performance than the pure feedforward FNN. Furthermore, the wavelet transform is an influential and effective technique used for the analysis of intricate time-varying signals. It offers numerous advantages and capabilities, making it a valuable tool in signal processing [28], and has been extensively studied for its applications that combine the learning capabilities of artificial neural network (ANN) and wavelet decomposition. Recently, researchers have proposed integrating wavelet functions into FNNs to create the wavelet fuzzy neural network (WFNN) with the goal of improving adaptive and learning capabilities for complex engineering problems [29]. By analyzing non-stationary signals to identify local details, reducing data complexity and handling uncertainty through fuzzy logic, and leveraging NNs’ self-learning characteristics to improve model accuracy, the WFNN is capable of describing nonlinear systems with uncertainties and possesses a fast learning capability. Additionally, this study proposes an intelligent control system using the capabilities of FNN, RNN, and WNN, where an online trained recurrent wavelet fuzzy neural network (RWFNN) [30] is utilized to enhance control performance.
High-performance applications of PMASynRM are limited due to nonlinear and time-varying control features. Therefore, the purpose of this article is to develop a high-performance PMASynRM servo drive that simultaneously achieves robust position control and high energy efficiency by using intelligent backstepping control recurrent wavelet fuzzy neural network (IBSCRWFNN) control with MTPA. To address the issue of MTPA, a Maxwell 2D simulation tool is employed in the design process of the PMASynRM. Then, the optimal current angle command for MTPA control is subsequently analyzed using finite element analysis (FEA), and apply the result by a lookup table (LUT) method to ensure proper functionality. Moreover, the PMASynRM servo drive is controlled by using the BSC to solve the presence of unavoidable uncertain system dynamics in the PMASynRM servo drive system. However, the bound of lumped uncertainty in the BSC is difficult to determine in real-life situations. To overcome this issue, the suggested method involves approximating the BSC using the RWFNN. Furthermore, this research incorporates an adaptive compensator to account for potential deviations resulting from the approximation of the RWFNN. In addition, the utilization of the Lyapunov stability method to generate online learning algorithms [22] for the IBSCRWFNN is proposed, ensuring robust performance in position control. The rest of this study is organized as follows: in Section 2, the focus will be on describing the modeling of the position servo drive system for PMASynRM with MTPA control based on the results of FEA. In Section 3, the PMASynRM servo drive is controlled by using the BSC to solve the presence of unavoidable uncertainties. To overcome the difficulty of the BSC, the RWFNN is discussed in Section 4. In Section 5, the Lyapunov stability method is proposed to generate online learning algorithms for the IBSCRWFNN. Section 6 will present the experimental results of the PI control, BSC, and the proposed IBSCRWFNN control. Finally, the research findings are thoroughly discussed in Section 7, presenting the conclusive remarks.

2. Modeling of PMASynRM Position Servo Drive System

2.1. Modeling of PMASynRM Servo Drive System

The stator voltage equations of PMASynRM in the d-q reference frame can be formulated as follows:
v d = R s i d + d d t λ d ω e λ q
v q = R s i q + d d t λ q + ω e λ d
where R s is the stator resistance; v d and v q are the d-axis and q-axis voltage; i d and i q are the d-axis and q-axis stator currents; λ d and λ q are the d-axis and q-axis flux linkages; ω e is the rotor electrical angular velocity. Furthermore, the flux linkage equations in the d-q reference frame are represented in the following:
λ q = L q i q
λ d = L d i d + λ m
The equations above denote the d-axis and q-axis inductances as L d and L q , respectively, and λ m as rotor PM flux. Neglecting magnetic saturation, the resulting electromagnetic torque in the d-q reference frame of PMASynRM can be represented by the following equation:
T e = 3 2 P 2 [ λ m i q + ( L d L q ) i q i d ]
where P is the pole number. Both PMASynRM and PMSM have electromagnetic torque consisting of reluctance torque and PM torque. However, in the case of PMASynRM, the reluctance torque plays the main driving role owing to the substantial disparity between the inductance values of the d-axis and q-axis inductances. The expression for the mechanical dynamic equation of the PMASynRM is as follows:
T e = J d ω r d t + B ω r + T L
J represents the inertia coefficient; B represents the damping coefficient; ω r is the speed response; T L represents the external load and friction torque.

2.2. PMASynRM Position Servo Drive System

The PMASynRM utilized in this study is a 4-pole, 36-slot motor with a power rating of 4.5 kW, voltage rating of 214 V, current rating of 9.4 A, a rated speed of 1500 rpm, and a maximum torque of 25 Nm. The mechanical designed parameters for the PMASynRM and the specifications of the servo drive system are outlined in detail in Table 1 and Table 2.
A typical PMASynRM position servo drive with coordinate transformation includes proportional position control, PI speed control, LUT for current angle, and PI current controllers, as illustrated in Figure 1. θ r is the position command; θ r is the position response; θ e is the rotor electrical position; ω r is the speed command. The d-q axis current commands are represented by i d and i q respectively. The three-phase currents are denoted by i a , i b , and i c . Similarly, the d-q axis voltage commands are represented by v d and v q respectively. By utilizing Hall current sensors with a transformation ratio of 1 V/6.67 A with an analog-to-digital converter (ADC), it is very effective to measure the three-phase currents. By implementing the space vector pulse width modulation (SVPWM) technique, control of the voltage source inverter (VSI) can be achieved. The VSI operates at a switching frequency of 20 kHz. Moreover, the VSI employs a silicon carbide (SiC) power MOSFET, which has a voltage rating of 900 V and a current rating of 36 A. Furthermore, using a load driver in torque control mode to control an industrial PMSM as the load for the PMASynRM, the performance of the position servo drive system for the PMASynRM is evaluated. In addition, a torque meter is connected to measure the torque output. The QEP interface is used to connect the encoder to the DSP, as illustrated in Figure 1.

2.3. Speed and Position Controllers Design

The speed and position controllers discussed in this article utilize PI controllers. These controllers are designed based on the controlled plants’ small signal models. The speed controller depicted in Figure 2, where s is the Laplace operator, assumes ideal conditions where T L = 0 and sets i d = 5 A . By referring to Equations (7) and (8), as well as Table 2, the torque constant can be deduced as K t = 1.2267 .
T e = K t i q
K t = 3 2 P 2 [ λ m + ( L d L q ) i d ]
Moreover, the bode diagrams of the speed-controlled plant and the loop gain are shown in Figure 3. The design requirements for the speed controller include a desired bandwidth (BW) f = 20   Hz and phase margin (PM) ϕ m = 70 . The graphical representation in Figure 3 reveals that the controlled plant necessitates compensation of −3.01 dB and −20.1°. The obtained result of the PI speed controller is
K P + K I s = 0.664 + 30.5385 s
where K P is a proportional gain; K I is an integral gain; s is the Laplace operator. Furthermore, the position controller depicted in Figure 4 utilizes a proportional controller K P P . In addition, the bode diagrams of the position-controlled plant and the loop gain are shown in Figure 5. The design requirements for the speed controller include a desired bandwidth (BW) f = 2 Hz and phase margin (PM) ϕ m = 70 . The graphical representation in Figure 5 reveals that the controlled plant necessitates compensation of 21.8 dB and −19.5°. Since the position controller consists of a single constant term, it can only compensate for the gain through a proportional gain K P P as follows:
K P P = 12 . 3

3. BSC System

To rewrite the ideal dynamic equation using Equations (5) and (6), it can be expressed as follows:
ω ˙ r = B ¯ J ¯ ω r + 3 P [ λ ¯ m + L ¯ d L ¯ q i d ] 4 J ¯ i q T L J ¯ = A m ω r + B m i q + C m T L
where A m = B ¯ J ¯ ; B m = 3 P [ λ ¯ m + L ¯ d L ¯ q i d ] 4 J ¯ ; C m = 1 J ¯ ; The symbol “ ¯ ” represents the nominal value. Considering the presence of uncertainties necessitates a rewrite of the dynamic Equation (11) as follows:
ω ˙ r = ( A m + Δ A m ) ω r + ( B m + Δ B m ) U + ( C m + Δ C m ) T L = A m ω r + B m U + F
where U = i q is the torque current command, the time-varying parameter variations are indicated by Δ A m , Δ B m , and Δ C m . Moreover, F represents the lumped uncertainty, which is defined as follows:
F = Δ A m ω r + Δ B m U + ( C m + Δ C m ) T L ,   F F b
and F b is defined as lumped uncertainty bound. Furthermore, the following definitions are used for the error in position tracking and its derivative:
e 1 = θ r ( t ) θ r ( t )
e ˙ 1 = θ ˙ r ( t ) θ ˙ r ( t )
The term θ ˙ r ( t ) = ω r ( t ) can be regarded as a virtual control input, and the following stabilizing function λ 1 is defined:
λ 1 = c 1 e 1 ω r ( t )
The constant c 1 is a positive value, and the first Lyapunov function is selected as
V 1 = 1 2 e 1 2 > 0
The function V 1 is positive definite. In addition, the definition of the virtual control error is as follows:
e 2 = ω r ( t ) + λ 1 = ω r ( t ) c 1 e 1 ω r ( t )
Obtaining the derivative of V 1 can be performed as follows:
V ˙ 1 = e 1 e ˙ 1 = e 1 ( e 2 c 1 e 1 ) = e 1 e 2 c 1 e 1 2
Assuming that e 2 = 0 is satisfied, then the derivative of V 1 will be negative. Additionally, the derivative of e 2 can be obtained as
e ˙ 2 = ω ˙ r ( t ) + λ ˙ 1 = ω ˙ r ( t ) c 1 e ˙ 1 ω ˙ r ( t )
The replacement of Equation (12) with Equation (20) yields the following equation:
e ˙ 2 = A m ω r ( t ) + B m U + F c 1 e ˙ 1 ω ˙ r ( t )
Then, the selection of the second Lyapunov function is performed as follows:
V 2 = 1 2 e 1 2 + 1 2 e 2 2 = V 1 + 1 2 e 2 2 > 0
where V 2 is a positive-definite function. The derivative of V 2 can be obtained by
V ˙ 2 = V ˙ 1 + e 2 e ˙ 2 = e 1 e 2 c 1 e 1 2 + e 2 e ˙ 2 = c 1 e 1 2 + e 2 ( e 1 + e ˙ 2 ) = c 1 e 1 2 + e 2 ( A m θ ˙ r ( t ) + B m U + F c 1 e ˙ 1 ω ˙ r ( t ) e 1 )
To ensure system stability based on Lyapunov’s condition, V ˙ 2 must be negative semidefinite. Consequently, utilizing Equation (23), a BSC control law is proposed as follows [22]:
U B S C = B m 1 [ A m θ ˙ r ( t ) + c 1 e ˙ 1 + ω ˙ r ( t ) + e 1 c 2 e 2 F b sgn ( e 2 ) ]
where c 2 is a positive constant; sgn ( ) is a sign function. The dynamic Equation (12) of the PMASynRM position servo drive system indicates that the implementation of the BSC law, as outlined in Equation (24), ensures system stability. By substituting Equation (24) into Equation (23), we can derive the resulting equation as follows:
V ˙ 2 = c 1 e 1 2 + e 2 [ c 2 e 2 F b sgn ( e 2 ) + F ] = c 1 e 1 2 c 2 e 2 2 e 2 F b + e 2 F c 1 e 1 2 c 2 e 2 2 e 2 F b + e 2 F c 1 e 1 2 c 2 e 2 2 e 2 ( F b F ) c 1 e 1 2 c 2 e 2 2 0 .
Thus, parametric uncertainty and external torque disturbance do not affect the stability of the BSC system. Figure 6 illustrates the control system’s capability to maintain stability even in the presence of disturbances. However, it is worth noting that the use of a sign function can lead to chattering phenomena. A boundary layer approach can be employed to mitigate this issue by substituting the sign function with a saturation function. This substitution helps reduce the occurrence of chattering phenomena.
s a t ( e 2 ) = 1 e 2 > ϕ e 2 ϕ , ϕ e 2 ϕ 1 e 2 < ϕ
The saturation function is denoted as s a t ; the boundary layer is set as ϕ > 0 . Thus, the BSC control law (24) is modified as follows:
U B S C = B m 1 [ A m θ ˙ r ( t ) + c 1 e ˙ 1 + ω ˙ r ( t ) + e 1 c 2 e 2 F b s a t ( e 2 ) ]

4. IBSCRWFNN System

The BSC system can ensure system stability when F F b . However, the lumped uncertainty is unknown in the real world, making it challenging to determine the upper bound F b . Moreover, asymptotic stability is a crucial requirement for position servo drives. In order to overcome the limitations associated with the BSC law described in Equation (27), an RWFNN controller [30] is proposed. The primary objective of designing the RWFNN controller is to achieve improved performance by providing an effective approximation of the BSC law. The control block diagram of the IBSCRWFNN system is shown in Figure 7. The control law for the IBSCRWFNN system is designed as follows to achieve asymptotic stability in position servo drives.
U = U ^ R W F N N + U ^ c
The RWFNN controller, represented by U ^ R W F N N , plays a crucial role in learning the BSC law to handle unknown system dynamics. Simultaneously, the compensator, denoted as U ^ c , is specifically designed to minimize the approximated error introduced by the RWFNN controller. This combination of U ^ R W F N N and U ^ c effectively addresses the unknown system dynamics and improves the overall performance of the control system.
Furthermore, the network structure is represented in Figure 8.
Figure 9 illustrates the flowchart outlining the proposed RWFNN controller. The detailed description of the operational mechanisms in the proposed RWFNN is as follows:
  • Measuring:
Utilizing the eQEP module in the DSP, the position response is measured with the assistance of an incremental encoder that has a resolution of 2500 counts/rev. Then, the RWFNN controller receives and utilizes e 1 = θ r ( t ) θ r ( t ) and e 2 = ω r ( t ) + λ 1 for generating control signals.
2.
RWFNN Input Layer:
Two input signals are fed into this layer of the proposed RWFNN controller: the tracking error of the rotor position x 1 1 = e 1 and the virtual control error x 2 1 = e 2 . To describe the input and output of each node i in this layer, the following expression is used:
n e t i 1 ( N ) = x i 1 ( N )
y i 1 ( N ) = f i 1 ( n e t i 1 ( N ) ) = n e t i 1 ( N ) ,   i = 1 ,   2
The network inputs are represented by n e t i 1 ( N ) , where the superscript and subscript correspond to the layer and node numbers, respectively. N denotes the sampling iteration number, while y i 1 ( N ) is the output of node ith. The unity function is denoted as f i 1 ( ) .
3.
RWFNN Membership Layer:
Layer 2 takes the outputs of layer 1 as its inputs. Additionally, the membership function utilized in this layer is the Gaussian function. The following elucidate the correlation between the input and output of each node in a comprehensive manner:
n e t j 2 ( N ) = ( x i 2 ( N ) m j 2 ) 2 ( σ j 2 ) 2
y j 2 ( N ) = f j 2 ( n e t j 2 ( N ) ) = exp ( n e t j 2 ( N ) ) ,   j = 1 , 2 , ,   6
The input is denoted by x i 2 ( N ) = y i 1 ( N ) ; the mean and standard deviation of the Gaussian function for node jth are represented by m j 2 and σ j 2 respectively; the output of node jth is denoted by y j 2 ( N ) ; f j 2 ( ) is an exponential function.
4.
RWFNN Wavelet Layer:
The propagation of signals in the wavelet layer is illustrated below:
ϕ i k 3 ( x ) = 1 σ i k 3 1 ( x i 1 ( N ) m i k 3 ) 2 ( σ i k 3 ) 2 × exp ( x i 1 ( N ) m i k 3 ) 2 ( σ i k 3 ) 2 ,   k = 1 , 2 , ,   9
ψ k 3 ( N ) = w i k 3 ϕ i k 3 ( x )
The input to node ith from layer 1, directed towards the wavelet function of node kth, is represented as ϕ i k 3 , and the connective weight as w i k 3 . In the wavelet layer, the output of node kth is represented by ψ k 3 . The dilation and translation variables of the wavelet function are expressed as m i k 3 and σ i k 3 , respectively.
5.
RWFNN Rule and Recurrent Layer:
The layer comprises a rule layer and a recurrent layer. Each node l , denoted as ∏, performs a multiplication operation on the input signals and outputs their product. Additionally, the nodes in the rule and recurrent layer employ a multiplication operation on the output signals derived from the membership layer, wavelet layer, and recurrent layer. This dynamic mapping process enhances the overall mapping capability of the system. The nodes are summarized as follows:
y j l 4 ( N ) = w j l 4 y j 2 ,   l = 1 , 2 , ,   9
n e t l 4 ( N ) = y j l 4 ψ k 3 w r 4 y l 4 ( N 1 )
y l 4 ( N ) = f l 4 ( n e t l 4 ( N ) ) = n e t l 4 ( N )
The output of the lth node in this layer is represented by y l 4 ( N ) . The calculation involves the utilization of the connecting weight, denoted as w j l 4 , between layer 2 and layer 4, and the recurrent weight is denoted by w r 4 . y l 4 ( N 1 ) denotes the previous output of node lth in this layer, and the unity function is denoted as f l 4 ( ) . Each node in the network incorporates a feedback loop using the recurrent technique to achieve dynamic mapping and higher sensitivity to previously obtained data.
6.
RWFNN Output Layer:
The inputs to layer 5 are obtained from the outputs of layer 4 and compute the final output by summing them up. In this layer, the output y o 5 ( N ) is mathematically expressed as follows:
n e t o 5 ( N ) = l 9 w l o 5 x l 5 ( N ) ,   o = 1
y o 5 ( N ) = f o 5 ( n e t o 5 ( N ) ) = n e t o 5 ( N )
The output of the rule layer is represented as x l 5 ( N ) = y l 4 ( N ) . The connective weight is represented by w l o 5 . The final output of the RWFNN is depicted as y o 5 ( N ) . The unity function is denoted as f o 5 ( ) .
7.
Online Network Parameters Learning:
All the adaptation laws of online network parameters learning will be given in Theorem 1 in the following section.

5. Stability Analysis of IBSCRWFNN System

The structure of the five-layer RWFNN, as illustrated in Figure 8, can be expressed as
U R W F N N   ( e 1 ,   e 2 ,   W ,   m j ,   σ j ,   m i k ,   σ i k ,   R ) W Γ W = [ w 1 5   w 2 5   w 3 5   w 4 5   w 5 5   w 6 5   w 7 5   w 8 5   w 9 5 ] R 1 × 9                     Γ = [ x 1 5   x 2 5   x 3 5   x 4 5   x 5 5   x 6 5   x 7 5   x 8 5   x 9 5 ] T R 9 × 1 m j = [ m 1 2   m 2 2   m 3 2   m 4 2   m 5 2   m 6 2 ] T R 6 × 1                                             σ j = [ σ 1 2   σ 2 2   σ 3 2   σ 4 2   σ 5 2   σ 6 2 ] T R 6 × 1 m i k = [ m 1 3   m 2 3   m 18 3 ] T R 18 × 1                                                           σ i k = [ σ 1 3   σ 2 3   σ 18 3 ] T R 18 × 1 R = [ w r 1 4   w r 2 4   w r 3 4   w r 4 4   w r 5 4   w r 6 4   w r 7 4   w r 8 4   w r 9 4 ] R 1 × 9
The universal approximation property guarantees the existence of an optimal U R W F N N for any nonlinear function. Consequently, a designed optimal U R W F N N is employed to learn the BSC law U B S C in order to achieve the following:
U B S C = U R W F N N ( e 1 ,   e 2 ,   W ,   m j ,   σ j ,   m i k ,   σ i k ,   R ) + ε = W Γ + ε
The reconstructed error is represented by ε , which is the minimum value; W , m j , σ j , m i k , σ i k , and R are the optimal values of W , m j , σ j , m i k , σ i k , and R respectively. Additionally, the control law illustrated in Equation (28) can be expressed as
U = U ^ R W F N N ( e 1 ,   e 2 ,   W ^ ,   m ^ j ,   σ ^ j ,   m ^ i k ,   σ ^ i k ,   R ^ ) + U ^ c = W ^ Γ ^ + U ^ c
where W ^ , m ^ j , σ ^ j , m ^ i k , σ ^ i k , and R ^ represent the estimated values of W , m j , σ j , m i k , σ i k , and R correspondingly. The equation below is obtained by subtracting (41) from (42):
U ˜ = U B S C U = U B S C U ^ R W F N N ( e 1 ,   e 2 ,   W ^ ,   m ^ j ,   σ ^ j ,   m ^ i k ,   σ ^ i k ,   R ^ ) U ^ c = U R W F N N ( e 1 ,   e 2 ,   W ,   m j ,   σ j ,   m i k ,   σ i k ,   R ) + ε U ^ R W F N N ( e 1 ,   e 2 ,   W ^ ,   m ^ j ,   σ ^ j ,   m ^ i k ,   σ ^ i k ,   R ^ ) U ^ c = W Γ + ε W ^ Γ ^ U ^ c = W ˜ Γ + W ^ Γ ˜ + ε U ^ c
where W ˜ = W W ^ and   Γ ˜ =   Γ   Γ ^ . A linearization technique is employed to convert the RWFNN into a partially linear form. This technique involves obtaining the Taylor series expansion of   Γ ˜ , which can be expressed as
Γ ˜ = Γ m j T m ˜ j + Γ σ j T σ ˜ j + Γ m i k T m ˜ i k + Γ σ i k T σ ˜ i k + Γ R T R ˜ + N h
where m ˜ j = m j m ^ j , σ ˜ j = σ j   σ ^ j , m ˜ i k = m i k m ^ i k , σ ˜ i k = σ i k   σ ^ i k , R ˜ = R R ^ ; the high-order term is represented by N h . In addition,
Γ m j T = x 1 5 m 1 x 1 5 m 2 x 1 5 m 6 x 2 5 m 1 x 2 5 m 2 x 2 5 m 6 x 9 5 m 1 x 9 5 m 2 x 9 5 m 6 R 9 × 6 Γ σ j T = x 1 5 σ 1 x 1 5 σ 2 x 1 5 σ 6 x 2 5 σ 1 x 2 5 σ 2 x 2 5 σ 6 x 9 5 σ 1 x 9 5 σ 2 x 9 5 σ 6 R 9 × 6 Γ m i k T = x 1 5 m 1 x 1 5 m 2 x 1 5 m 18 x 2 5 m 1 x 2 5 m 2 x 2 5 m 18 x 9 5 m 1 x 9 5 m 2 x 9 5 m 18 R 9 × 18 Γ σ i k T = x 1 5 σ 1 x 1 5 σ 2 x 1 5 σ 18 x 2 5 σ 1 x 2 5 σ 2 x 2 5 σ 18 x 9 5 σ 1 x 9 5 σ 2 x 9 5 σ 18 R 9 × 18 Γ R T = x 1 5 w r 1 x 1 5 w r 2 x 1 5 w r 9 x 2 5 w r 1 x 2 5 w r 2 x 2 5 w r 9 x 9 5 w r 1 x 9 5 w r 2 x 9 5 w r 9 R 9 × 9
Rewriting (44), Γ can be calculated as follows:
Γ = Γ ^ + Γ ˜ = Γ ^ + Γ m j T m ˜ j + Γ σ j T σ ˜ j + Γ m i k T m ˜ i k + Γ σ i k T σ ˜ i k + Γ R T R ˜ + N h
Substituting (44) and (45) into (43), expressing the estimated error in Equation (46) can be performed in the following manner:
U ˜ = W ˜ Γ * + W ^   Γ ˜ + ε U ^ c = W ˜ ( Γ ^ + Γ m j T m ˜ j + Γ σ j T σ ˜ j + Γ m i k T m ˜ i k + Γ σ i k T σ ˜ i k + Γ R T R ˜ + N h ) + W ^ ( Γ m j T m ˜ j + Γ σ j T σ ˜ j + Γ m i k T m ˜ i k + Γ σ i k T σ ˜ i k + Γ R T R ˜ + N h ) + ε U ^ c = W ˜ Γ ^ + W ˜ Γ m j T m ˜ j + W ˜ Γ σ j T σ ˜ j + W ˜ Γ m i k T m ˜ i k + W ˜ Γ σ i k T σ ˜ i k + W ˜ Γ R T R ˜ + W ˜ N h + W ^ Γ m j T m ˜ j + W ^ Γ σ j T σ ˜ j + W ^ Γ m i k T m ˜ i k + W ^ Γ σ i k T σ ˜ i k + W ^ Γ R T R ˜ + W ^   N h + ε U ^ c = W ˜ Γ ^ + W ^ Γ m j T m ˜ j + W ^ Γ σ j T σ ˜ j + W ^ Γ m i k T m ˜ i k + W ^ Γ σ i k T σ ˜ i k + W ^ Γ R T R ˜ + W ˜ Γ m j T m ˜ j + W ˜ Γ σ j T σ ˜ j + W ˜ Γ m i k T m ˜ i k + W ˜ Γ σ i k T σ ˜ i k + W ˜ Γ R T R ˜ + W * N h + ε U ^ c = W ˜ Γ ^ + W ^ Γ m j T m ˜ j + W ^ Γ σ j T σ ˜ j + W ^ Γ m i k T m ˜ i k + W ^ Γ σ i k T σ ˜ i k + W ^ Γ R T R ˜ U ^ c + H
The H , which is named as the uncertain term, can be expressed as follows:
H = W ˜ Γ m j T m ˜ j + W ˜ Γ σ j T σ ˜ j + W ˜ Γ m i k T m ˜ i k + W ˜ Γ σ i k T σ ˜ i k + W ˜ Γ R T R ˜ + W N h + ε
Theorem 1.
Given the PMASynRM servo drive system described in (12), the proposed IBSCRWFNN achieves absolute asymptotic stability under the following condition.
1.
Implementation of the IBSCRWFNN control as illustrated in (28);
2.
Adoption of the RWFNN adaptation law as described in (48)–(53);
3.
The compensators, illustrated in Equations (54) and (55), are developed with an adaptive law.
W ^ ˙ T = η w e 2 Γ ^
m ^ ˙ j T = η m j e 2 W ^ Γ m j T
σ ^ ˙ j T = η σ j e 2 W ^ Γ σ j T
m ^ ˙ i k T = η m i k e 2 W ^ Γ m i k T
σ ^ ˙ i k T = η σ i k e 2 W ^ Γ σ i k T
R ^ ˙ T = η r e 2 W ^ Γ R T
U ^ c = ψ ^
ψ ^ ˙ = γ e 2
where η w , η m j , η σ j , η m i k , η σ i k , η r are positive constant learning rate parameters; ψ ^ represents the value of the estimated online approximated error; and γ is a positive constant.
Proof. 
The proposed IBSCRWFNN is designed with a Lyapunov function given by
V 3 ( e 1 ( t ) ,   e 2 ( t ) ,   ψ ˜ ( t ) ,   W ˜ ,   m ˜ j ,   σ ˜ j ,   m ˜ i k ,   σ ˜ i k ,   R ˜ ) = 1 2 e 1 2 + 1 2 e 2 2 + B m 2 η w W ˜ W ˜ T + B m 2 η m j m ˜ j T m ˜ j + B m 2 η σ j σ ˜ j T σ ˜ j + B m 2 η m i k m ˜ i k T m ˜ i k + B m 2 η σ i k σ ˜ i k T σ ˜ i k + B m 2 η r R ˜ T R ˜ + B m 2 γ ψ ˜ 2 > 0
The function V 3 is chosen to be positive-definite, and ψ ˜ = ψ ψ ^ . ψ is the symbol used to denote the approximated error, and it is defined by ψ = H F / B m . Furthermore, the approximated error ψ is assumed to be bounded by ψ F b . Given that the sampling interval in the experiment is considerably shorter than the fluctuations observed in H and F , the approximated error ψ is treated as a constant during the estimation process. However, it is difficult to know the upper bound F b . Hence, a proposed adaptation law is put forth to modify the value of the online estimated approximated error ψ ^ within the compensator. Differentiating ψ ˜ with respect to time yields ψ ˜ ˙ = ψ ^ ˙ . By utilizing Equation (12) and taking the derivative of V 3 with respect to time, the following expression can be derived:
V ˙ 3 = e 1 e ˙ 1 + e 2 e ˙ 2 B m η w W ˜ W ^ ˙ T B m η m j m ^ ˙ j T m ˜ j B m η σ j σ ^ ˙ j T σ ˜ j B m η m i k m ^ ˙ i k T m ˜ i k B m η σ i k σ ^ ˙ i k T σ ˜ i k B m η r R ^ ˙ T R ˜ B m γ ψ ˜ ψ ^ ˙ = e 1 e 2 c 1 e 1 e 2 B m U B S C U e 1 + c 2 e 2 F B m η w W ˜ W ^ ˙ T B m η m j m ^ ˙ j T m ˜ j B m η σ j σ ^ ˙ j T σ ˜ j B m η m i k m ^ ˙ i k T m ˜ i k B m η σ i k σ ^ ˙ i k T σ ˜ i k B m η r R ^ ˙ T R ˜ B m γ ψ ˜ ψ ^ ˙ = c 1 e 1 2 c 2 e 2 2 B m η w W ˜ W ^ ˙ T e 2 B m W ˜ Γ ^ B m η m j m ^ ˙ j T m ˜ j e 2 B m W ^ Γ m j T m ˜ j B m η σ j σ ^ ˙ j T σ ˜ j e 2 B m W ^ Γ σ j T σ ˜ j B m η m i k m ^ ˙ i k T m ˜ i k e 2 B m W ^ Γ m i k T m ˜ i k B m η σ i k σ ^ ˙ i k T σ ˜ i k e 2 B m W ^ Γ σ i k T σ ˜ i k B m η r R ^ ˙ T R ˜ e 2 B m W ^ Γ R T R ˜ + e 2 B m U ^ c e 2 B m ψ B m γ ψ ˜ ψ ^ ˙ = c 1 e 1 2 c 2 e 2 2 B m η w W ˜ W ^ ˙ T + e 2 B m W ˜ Γ ^ B m η m j m ^ ˙ j T m ˜ j + e 2 B m W ^ Γ m j T m ˜ j B m η σ j σ ^ ˙ j T σ ˜ j + e 2 B m W ^ Γ σ j T σ ˜ j B m η m i k m ^ ˙ i k T m ˜ i k + e 2 B m W ^ Γ m i k T m ˜ i k B m η σ i k σ ^ ˙ i k T σ ˜ i k + e 2 B m W ^ Γ σ i k T σ ˜ i k B m η r R ^ ˙ T R ˜ + e 2 B m W ^ Γ R T R ˜ + e 2 B m U ^ c e 2 B m ψ ^ B m γ ψ ˜ [ ψ ^ ˙ + γ e 2 ]
Moreover, by substituting (48)–(55) into (57), it can be concluded that
V ˙ 3 ( e 1 ( t ) ,   e 2 ( t ) ,   ψ ˜ ( t ) ,   W ˜ ,   m ˜ j ,   σ ˜ j ,   m ˜ i k ,   σ ˜ i k ,   R ˜ ) = c 1 e 1 2 c 2 e 2 2 0
Since V ˙ 3 ( e 1 ( t ) ,   e 2 ( t ) ,   ψ ˜ ( t ) ,   W ˜ ,   m ˜ j ,   σ ˜ j ,   m ˜ i k ,   σ ˜ i k ,   R ˜ ) 0 is negative semidefinite,
V ˙ 3 ( e 1 ( t ) ,   e 2 ( t ) ,   ψ ˜ ( t ) ,   W ˜ ,   m ˜ j ,   σ ˜ j ,   m ˜ i k ,   σ ˜ i k ,   R ˜ ) < V ˙ 3 ( e 1 ( 0 ) ,   e 2 ( 0 ) ,   ψ ˜ ( 0 ) ,   W ˜ ,   m ˜ j ,   σ ˜ j ,   m ˜ i k ,   σ ˜ i k ,   R ˜ )
which implies that e 1 ( t ) ,   e 2 ( t ) ,   ψ ˜ ( t ) ,   W ˜ ,   m ˜ j ,   σ ˜ j ,   m ˜ i k ,   σ ˜ i k , and R ˜ are all bounded. By defining Ω ( t ) = c 1 e 1 2 + c 2 e 2 2 = V ˙ 3 ( t ) and integrating with respect to time, one can obtain the following equation:
0 t Ω ( τ ) d τ = V 3 ( 0 ) V 3 ( t )
Since V ˙ 3 ( e 1 ( 0 ) ,   e 2 ( 0 ) ,   ψ ˜ ( 0 ) ,   W ˜ ,   m ˜ j ,   σ ˜ j ,   m ˜ i k ,   σ ˜ i k ,   R ˜ ) is bounded and V ˙ 3 ( e 1 ( t ) ,   e 2 ( t ) ,   ψ ˜ ( t ) ,   W ˜ ,   m ˜ j ,   σ ˜ j ,   m ˜ i k ,   σ ˜ i k ,   R ˜ ) is also bounded and nonincreasing, thus
lim t 0 t Ω ( τ ) d τ <
Furthermore, the boundedness of Ω ˙ ( t ) implies that Ω ( t ) is uniformly continuous. Applying Barbalat’s Lemma, it can be demonstrated that lim t Ω ( t ) 0 . As a result, both e 1 and e 2 will approach zero as t . Consequently, the proposed IBSCRWFNN system exhibits asymptotic stability [22]. □

6. Experimentation

The experimental setup, depicted in Figure 10, comprises various components such as the PMASynRM servo drive, DSP TMS320F28075 board, and a SiC-based VSI with 4.5 kW. An industrial 7.5-kW PMSM drive is operated in torque control mode as the load. Moreover, two load torques 10 Nm (case 1) and 20 Nm (case 2), are set in the experimentation. The control of position and speed of the PMASynRM is determined using an incremental encoder, which interfaces with a quadrature encoder pulse (QEP) interface and has a sampling interval of 1 ms. The control of current operates with a sampling interval of 0.1 ms. The PMASynRM servo drive is then controlled by delivering switching commands for space vector pulse width modulation (SVPWM) to the voltage source inverter (VSI).
The suggested position control system is evaluated based on three performance metrics: T M , T a v e r , and T s d . These metrics represent the maximum tracking error, average tracking error, and standard deviation of the tracking error, respectively. They are utilized to assess and validate the control performance.
T M = max N ( T e r r o r ( N ) )
T a v e r = N = 1 h T e r r o r ( N ) h
T s d = N = 1 h ( T e r r o r ( N ) T a v e r ( N ) ) 2 h
Given that T e r r o r ( N ) = θ r ( N ) θ r ( N ) and h represents the total number of iterations, the control performance of the system is demonstrated by measuring the responses of periodical step and sinusoid commands. To model the periodical step reference input, a second-order transfer function with a rise time of 0.6 s is utilized as the reference model in the following:
ω n 2 s 2 + 2 ξ ω n s + ω n 2 = 30 s 2 + 11 s + 30
In Equation (64), ξ and ω n represent the damping ratio and undamped natural frequency, correspondingly. Moreover, the control performance of the proposed IBSCRWFNN position controller is compared with that of the BSC position controller through experimental results analysis. Furthermore, to compare the control performance, experimental results of the PI control, BSC, and the proposed IBSCRWFNN control are presented and analyzed. The parameters of PI control have been designed in Section 2.3, and the parameters of BSC and the proposed IBSCRWFNN control are provided as follows:
c 1 = 1.2 ,   c 2 = 0 . 45 ,   F b = 10 ,   η w = 0 . 05 ,   η m j = 0 . 15 ,   η σ j = 0 . 3 ,   η m i k = 0 . 01 ,   η σ i k = 0 . 01 ,   η r = 0.02
The parameters are iteratively adjusted to achieve optimal transient control performance while ensuring stability using a trial-and-error process. In addition, in order to strike a balance between computational resources and control performance, the network structure of the RWFNN has been designed with specific numbers of neurons in each layer: 2 in the input layer, 6 in the membership layer, 27 in the wavelet layer, 18 in the rule layer, and 1 in the output layer. Additionally, for the 32-bit floating-point DSP with 120 MHz using the “C” program, the total operation cycles and execution time for the PI controller are 60 and 0.0005 ms; the proposed BSC controller are 393 and 0.003275 ms; the proposed IBSCRWFNN controller are 9437 and 0.0786 ms. Consequently, the total execution time of the proposed IBSCRWFNN controller remains below 1 ms, which aligns with the sampling interval of the speed control loop.
In the experimentation, the objective is the control of the rotor position of PMASynRM to periodically track step and sinusoid position commands with minimum tracking errors. The test scenarios are outlined in Table 3 for the assessment of the robustness of various controllers under different operating conditions.
Figure 11 illustrates the experimental results of command tracking using periodical step commands for both case 1 and case 2 of the PI control system. The position command, response, and error are shown in Figure 11a,d; the current commands are shown in Figure 11b,e; the speed command and response are shown in Figure 11c,f. Figure 12 and Figure 13 illustrate the experimental results of command tracking using periodical step and sinusoid commands for both case 1 and case 2 of the BSC control system. The position command, response, and error are shown in Figure 12a,d and Figure 13a,d; the current commands are shown in Figure 12b,e and Figure 13b,e; the speed command and response are shown in Figure 12c,f and Figure 13c,f. In addition, Figure 14 and Figure 15 illustrate the experimental results of command tracking using periodical step and sinusoid commands for both case 1 and case 2 of the IBSCRWFNN control system. The position command, response, and error are shown in Figure 14a,d and Figure 15a,d; the current commands are shown in Figure 14b,e and Figure 15b,e; the speed command and response are shown in Figure 14c,f and Figure 15c,f.
From the experimental results, it can be observed that the d-axis current command is effectively generated using the FEA-based look-up table (LUT) for MTPA control. Moreover, the BSC position controller performs better than the PI controller, and the proposed IBSCRWFNN controller outperforms the BSC controller. The rotor response of the PMASynRM is significantly enhanced by the proposed IBSCRWFNN position controller, resulting in reduced tracking errors under different reference inputs. This improvement can be attributed to the parallel processing and online learning capabilities of the RWFNN used in the control network. In other words, the robustness of the position control is improved by employing the suggested IBSCRWFNN position controller. Furthermore, the quantified results of maximum tracking error and transient response time of all experiments are also presented in Table 3. In addition, the performance measurements of PI, BSC, and the proposed IBSCRWFNN position controllers are compared in Figure 16, considering two operating cases with periodical steps and sinusoid reference commands. The proposed IBSCRWFNN position controller exhibits lower maximum, average, and standard deviation tracking errors thanks to its faster convergence rate and improved generalization performance.

7. Conclusions

In this study, an IBSCRWFNN control was proposed for a high-performance PMASynRM servo drive system. First, the dynamic model of the PMASynRM servo drive was analyzed using ANSYS Maxwell-2D capabilities. The FEA results were utilized to generate a LUT for the MTPA current angle command. Subsequently, a BSC position tracking system was developed to confront the existing lumped uncertainty of the motor drive. Moreover, the proposed RWFNN was employed as an alternative to the BSC law to address the challenges associated with the dynamic model of the motor required by the BSC in the PMASynRM servo drive. Furthermore, the Lyapunov stability method was employed to derive online learning algorithms for the RWFNN, ensuring asymptotical stability. Finally, the experimental results demonstrate that the proposed IBSCRWFNN exhibits excellent control performance in terms of position tracking control for the PMASynRM servo drive. This study presents several significant contributions, which include: (1) the successful creation of the IBSCRWFNN specifically designed for a high-performance PMASynRM position servo drive system; (2) the successful development of an online learning algorithm that allows for the real-time training of the RWFNN using the Lyapunov stability theorem; (3) the effective implementation of the IBSCRWFNN in a floating point DSP, ensuring robust position control performance for the high-performance PMASynRM.
The future works of this study are as follows: (1) in the experiment, an optical encoder was used to obtain the position of the motor rotor. Subsequently, a sensorless control method, which eliminates the need for sensors, can be further incorporated to reduce system costs. (2) In this study, the online learning rates of the intelligent control algorithm are adjusted through trial and error. In the future, it is possible to explore a self-adjusting network learning rate to optimize the intelligent control algorithm. (3) In addition to reducing motor copper losses by maximizing torque per ampere, the development of optimal efficiency control is also possible.

Author Contributions

F.-J.L. designed and developed the main parts of research work, including theory derivation and analyses of the obtained results. M.-S.H. were mainly responsible for the preparation of the paper. Y.-C.C. and S.-G.C. contributed to the DSP-based control platform and writing parts. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science and Technology Council of Taiwan under Grant Most 110-2221-E-008-054-MY3.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lara, J.; Xu, J.; Chandra, A. Effects of Rotor Position Error in the Performance of Field Oriented Controlled PMSM Drives for Electric Vehicle Traction Applications. IEEE Trans. Ind. Electron. 2016, 63, 4738–4751. [Google Scholar] [CrossRef]
  2. Yang, Y.; Sandra, M.C.; Rong, Y.; Berker, B.; Anand, S.; Hossein, D.; Ali, E. Design and Comparison of Interior Permanent Magnet Motor Topologies for Traction Applications. IEEE Trans. Transp. Electrif. 2017, 3, 86–97. [Google Scholar] [CrossRef]
  3. Fang, S.; Liu, H.; Wang, H.; Yang, H.; Lin, H. High Power Density PMSM with Lightweight Structure and High-Performance Soft Magnetic Alloy Core. IEEE Trans. Appl. Supercond. 2019, 29, 0602805. [Google Scholar] [CrossRef]
  4. Kong, Y.; Lin, M.; Jia, L. A Novel High Power Density Permanent-Magnet Synchronous Machine with Wide Speed Range. IEEE Trans. Magn. 2020, 56, 7505206. [Google Scholar] [CrossRef]
  5. Park, S.-H.; Lee, E.-C.; Park, J.-C.; Hwang, S.-W.; Lim, M.-S. Prediction of Mechanical Loss for High-Power-Density PMSM Considering Eddy Current Loss of PMs and Conductors. IEEE Trans. Magn. 2021, 57, 6300205. [Google Scholar] [CrossRef]
  6. Bedetti, N.; Calligaro, S.; Petrella, R. Stand-Still Self-Identification of Flux Characteristics for Synchronous Reluctance Machines Using Novel Saturation Approximating Function and Multiple Linear Regression. IEEE Trans. Ind. Appl. 2016, 52, 3083–3092. [Google Scholar] [CrossRef]
  7. Nikmaram, B.; Davari, S.A.; Naderi, P.; Garcia, C.; Rodriguez, J. Sensorless Simplified Finite Control Set Model Predictive Control of SynRM Using Finite Position Set Algorithm. IEEE Access 2021, 9, 47184–47193. [Google Scholar] [CrossRef]
  8. Taghavi, S.M.; Pillay, P. A Mechanically Robust Rotor with Transverse Laminations for a Wide-Speed-Range Synchronous Reluctance Traction Motor. IEEE Trans. Ind. Appl. 2015, 51, 4404–4414. [Google Scholar] [CrossRef]
  9. Ferrari, M.; Bianchi, N.; Doria, A.; Fornasiero, E. Design of synchronous reluctance motor for hybrid electric vehicles. IEEE Trans. Ind. Appl. 2015, 51, 3030–3040. [Google Scholar] [CrossRef]
  10. Kerdsup, B.; Takorabet, N.; Nahidmobarakeh, B. Design of Permanent Magnet-Assisted Synchronous Reluctance Motors with Maximum Efficiency-Power Factor and Torque per Cost. In Proceedings of the 2018 XIII International Conference on Electrical Machines (ICEM), Alexandroupoli, Greece, 3–6 September 2018; pp. 2465–2471. [Google Scholar]
  11. Wu, W.; Zhu, X.; Quan, L.; Du, Y.; Xiang, Z.; Zhu, X. Design and Analysis of a Hybrid Permanent Magnet Assisted Synchronous Reluctance Motor Considering Magnetic Saliency and PM Usage. IEEE Trans. Appl. Supercond. 2018, 28, 5200306. [Google Scholar] [CrossRef]
  12. Nasiri-Zarandi, R.; Karami-Shahnani, A.; Toulabi, M.S.; Tessarolo, A. Design and Experimental Performance Assessment of an Outer Rotor PM-Assisted SynRM for the Electric Bike Propulsion. IEEE Trans. Transp. Electrif. 2023, 9, 727–736. [Google Scholar] [CrossRef]
  13. Liu, Z.; Hu, Y.; Wu, J.; Zhang, B.; Feng, G. A Novel Modular Permanent Magnet-Assisted Synchronous Reluctance Motor. IEEE Access 2021, 9, 19947–19959. [Google Scholar] [CrossRef]
  14. Li, K.; Wang, Y. Maximum Torque Per Ampere (MTPA) Control for IPMSM Drives Based on a Variable-Equivalent-Parameter MTPA Control Law. IEEE Trans. Power Electron. 2019, 34, 7092–7102. [Google Scholar] [CrossRef]
  15. Lin, F.J.; Liu, Y.T.; Yu, W.A. Power Perturbation Based MTPA with an Online Tuning Speed Controller for an IPMSM Drive System. IEEE Trans. Ind. Electron. 2018, 65, 3677–3687. [Google Scholar] [CrossRef]
  16. Li, K.; Wang, Y. Maximum Torque per Ampere (MTPA) Control for IPMSM Drives Using Signal Injection and an MTPA Control Law. IEEE Trans. Ind. Inform. 2019, 15, 5588–5598. [Google Scholar] [CrossRef]
  17. Yousefi-Talouki, A.; Pescetto, P.; Pellegrino, G. Sensorless Direct Flux Vector Control of Synchronous Reluctance Motors Including Standstill, MTPA, and Flux Weakening. IEEE Trans. Ind. Appl. 2017, 53, 3598–3608. [Google Scholar] [CrossRef] [Green Version]
  18. Dianov, A.; Tinazzi, F.; Calligaro, S.; Bolognani, S. Review and Classification of MTPA Control Algorithms for Synchronous Motors. IEEE Trans. Power Electron. 2022, 37, 3990–4007. [Google Scholar] [CrossRef]
  19. Kim, H.S.; Lee, Y.; Sul, S.K.; Yu, J.; Oh, J. Online MTPA Control of IPMSM Based on Robust Numerical Optimization Technique. IEEE Trans. Ind. Appl. 2019, 55, 3736–3746. [Google Scholar] [CrossRef]
  20. Lin, F.J.; Huang, M.S.; Hung, C.Y.; Chien, Y.C. Intelligent Computed Torque Control with Recurrent Legendre Fuzzy Neural Network for Permanent-Magnet Assisted Synchronous Reluctance Motor. IEEE Access 2023, 11, 54017–54028. [Google Scholar] [CrossRef]
  21. Liu, T.-H.; Pu, H.-T.; Lin, C.-K. Implementation of an Adaptive Position Control System of a Permanent-Magnet Synchronous Motor and Its Application. IET Electr. Power Appl. 2010, 4, 121–130. [Google Scholar] [CrossRef]
  22. Lin, F.J.; Chen, S.G.; Hsu, C.W. Intelligent backstepping control using recurrent feature selection fuzzy neural network for synchronous reluctance motor position servo drive system. IEEE Trans. Fuzzy Syst. 2019, 27, 413–427. [Google Scholar] [CrossRef]
  23. Mansouri, M.M.; Hadjeri, S.; Brahami, M. New method of detection, identification, and elimination of photovoltaic system faults in real time based on the adaptive Neuro-fuzzy system. IEEE J. Photovolt. 2021, 11, 797–805. [Google Scholar] [CrossRef]
  24. Dong, C.; Yu, Z.; Chen, X.; Chen, H.; Huang, Y.; Huang, Q. Adaptability Control Towards Complex Ground Based on Fuzzy Logic for Humanoid Robots. IEEE Trans. Fuzzy Syst. 2022, 30, 1574–1584. [Google Scholar] [CrossRef]
  25. Chen, S.Y.; Liu, T.S. Intelligent tracking control of a PMLSM using self-evolving probabilistic fuzzy neural network. IET Electr. Power Appl. 2017, 11, 1043–1054. [Google Scholar] [CrossRef]
  26. Lin, F.J.; Sun, I.F.; Yang, K.J.; Chang, J.K. Recurrent Fuzzy Neural Cerebellar Model Articulation Network Fault-Tolerant Control of Six-Phase Permanent Magnet Synchronous Motor Position Servo Drive. IEEE Trans. Fuzzy Syst. 2016, 24, 153–167. [Google Scholar] [CrossRef]
  27. Fei, J.; Chen, Y.; Liu, L.; Fang, Y. Fuzzy Multiple Hidden Layer Recurrent Neural Control of Nonlinear System Using Terminal Sliding-Mode Controller. IEEE Trans. Cybern. 2022, 52, 9519–9534. [Google Scholar] [CrossRef] [PubMed]
  28. Khan, M.A.; Uddin, M.N.; Rahman, M.A. A Novel Wavelet-Neural-Network-Based Robust Controller for IPM Motor Drives. IEEE Trans. Ind. Appl. 2013, 49, 2341–2351. [Google Scholar] [CrossRef]
  29. Huang, W.; Oh, S.; Pedrycz, W. Fuzzy Wavelet Polynomial Neural Networks: Analysis and Design. IEEE Trans. Fuzzy Syst. 2017, 25, 1329–1341. [Google Scholar] [CrossRef]
  30. Lin, F.J.; Tan, K.H.; Luo, W.C.; Xiao, G.D. Improved LVRT performance of PV power plant using recurrent wavelet fuzzy neural network control for weak grid conditions. IEEE Access 2020, 8, 69346–69358. [Google Scholar] [CrossRef]
Figure 1. PI Controlled PMASynRM position servo drive.
Figure 1. PI Controlled PMASynRM position servo drive.
Energies 16 05389 g001
Figure 2. PI speed controller.
Figure 2. PI speed controller.
Energies 16 05389 g002
Figure 3. Bode diagrams of speed controller. (a) Magnitude; (b) Phase.
Figure 3. Bode diagrams of speed controller. (a) Magnitude; (b) Phase.
Energies 16 05389 g003
Figure 4. Proportional position controller.
Figure 4. Proportional position controller.
Energies 16 05389 g004
Figure 5. Bode diagrams of position controller. (a) Magnitude; (b) Phase.
Figure 5. Bode diagrams of position controller. (a) Magnitude; (b) Phase.
Energies 16 05389 g005
Figure 6. Control block diagram of BSC controlled PMASynRM position servo drive.
Figure 6. Control block diagram of BSC controlled PMASynRM position servo drive.
Energies 16 05389 g006
Figure 7. IBSCRWFNN controlled PMASynRM position servo drive.
Figure 7. IBSCRWFNN controlled PMASynRM position servo drive.
Energies 16 05389 g007
Figure 8. Network structure of RWFNN.
Figure 8. Network structure of RWFNN.
Energies 16 05389 g008
Figure 9. RWFNN flowchart for position control.
Figure 9. RWFNN flowchart for position control.
Energies 16 05389 g009
Figure 10. Experimental setup.
Figure 10. Experimental setup.
Energies 16 05389 g010
Figure 11. Experimental results of PI control with periodical step command. (a) Position command, response, and error at case 1; (b) d-q axis current commands at case 1; (c) Speed command and response at case 1; (d) Position command, response, and error at case 2; (e) d-q axis current commands at case 2; (f) Speed command and response at case 2.
Figure 11. Experimental results of PI control with periodical step command. (a) Position command, response, and error at case 1; (b) d-q axis current commands at case 1; (c) Speed command and response at case 1; (d) Position command, response, and error at case 2; (e) d-q axis current commands at case 2; (f) Speed command and response at case 2.
Energies 16 05389 g011
Figure 12. Experimental results of BSC control with periodical step command. (a) Position command, response, and error at case 1; (b) d-q axis current commands at case 1; (c) Speed command and response at case 1; (d) Position command, response, and error at case 2; (e) d-q axis current commands at case 2; (f) Speed command and response at case 2.
Figure 12. Experimental results of BSC control with periodical step command. (a) Position command, response, and error at case 1; (b) d-q axis current commands at case 1; (c) Speed command and response at case 1; (d) Position command, response, and error at case 2; (e) d-q axis current commands at case 2; (f) Speed command and response at case 2.
Energies 16 05389 g012
Figure 13. Experimental results of BSC control with periodical sinusoid command. (a) Position command, response, and error at case 1; (b) d-q axis current commands at case 1; (c) Speed command and response at case 1; (d) Position command, response, and error at case 2; (e) d-q axis current commands at case 2; (f) Speed command and response at case 2.
Figure 13. Experimental results of BSC control with periodical sinusoid command. (a) Position command, response, and error at case 1; (b) d-q axis current commands at case 1; (c) Speed command and response at case 1; (d) Position command, response, and error at case 2; (e) d-q axis current commands at case 2; (f) Speed command and response at case 2.
Energies 16 05389 g013
Figure 14. Experimental results of IBSCRWFNN control with periodical step command. (a) Position command, response, and error at case 1; (b) d-q axis current commands at case 1; (c) Speed command and response at case 1; (d) Position command, response, and error at case 2; (e) d-q axis current commands at case 2; (f) Speed command and response at case 2.
Figure 14. Experimental results of IBSCRWFNN control with periodical step command. (a) Position command, response, and error at case 1; (b) d-q axis current commands at case 1; (c) Speed command and response at case 1; (d) Position command, response, and error at case 2; (e) d-q axis current commands at case 2; (f) Speed command and response at case 2.
Energies 16 05389 g014
Figure 15. Experimental results of IBSCRWFNN control with periodical sinusoid command. (a) Position command, response, and error at case 1; (b) d-q axis current commands at case 1; (c) Speed command and response at case 1; (d) Position command, response, and error at case 2; (e) d-q axis current commands at case 2; (f) Speed command and response at case 2.
Figure 15. Experimental results of IBSCRWFNN control with periodical sinusoid command. (a) Position command, response, and error at case 1; (b) d-q axis current commands at case 1; (c) Speed command and response at case 1; (d) Position command, response, and error at case 2; (e) d-q axis current commands at case 2; (f) Speed command and response at case 2.
Energies 16 05389 g015
Figure 16. Maximum, average, and standard deviation of tracking errors used for the PI, BSC, and IBSCRWFNN. (a) Periodical step command at case 1; (b) Periodical step command at case 2; (c) Periodical sinusoid command at case 1; (d) Periodical sinusoid command at case 2.
Figure 16. Maximum, average, and standard deviation of tracking errors used for the PI, BSC, and IBSCRWFNN. (a) Periodical step command at case 1; (b) Periodical step command at case 2; (c) Periodical sinusoid command at case 1; (d) Periodical sinusoid command at case 2.
Energies 16 05389 g016
Table 1. Mechanical designed parameters of PMASynRM.
Table 1. Mechanical designed parameters of PMASynRM.
ParametersValues
Pole number4 pole
Slot number36 slot
Air gap length0.3 mm
Rotor inner diameter31 mm
Rotor outer diameter94.4 mm
Stator inner diameter95 mm
Stator outer diameter160 mm
Stack length150 mm
Table 2. Specifications of PMASynRM servo drive.
Table 2. Specifications of PMASynRM servo drive.
ParametersValues
Power rating4.5 kW
Phase voltage rating214 V
Phase current rating9.4 Arms (delta)
Speed rating1500 rpm
Torque rating25 Nm
d-axis inductance19.6 mH
q-axis inductance84.3 mH
Stator resistance1.01 Ω
Inertia coefficient0.0069 Nm/(rad/s2)
Damping coefficient0.0013 Nm/(rad/s)
Magnetic flux0.0854 Wb
DC-link540 V
Sampling time of current/speed and position loop0.1 ms/1 ms
Switching frequency20 kHz
Encoder2500 counts/turn
Table 3. Test scenarios and quantified results of experiment.
Table 3. Test scenarios and quantified results of experiment.
CaseControllerPosition Command
(Degree)
Maximum Tracking Error
(Degree)
Transient Response Time
(s)
Load Torque
(Nm)
Figure
Case 1PIperiodical step44.51.2510Figure 11a
Case 2PIperiodical step50.51.3120Figure 11c
Case 1BSCperiodical step28.31.1510Figure 12a
Case 2BSCperiodical step341.1920Figure 12c
Case 1BSCperiodical sinusoid13.20.8210Figure 13a
Case 2BSCperiodical sinusoid150.8520Figure 13c
Case 1IBSCRWFNNperiodical step14.50.6210Figure 14a
Case 2IBSCRWFNNperiodical step15.30.5920Figure 14c
Case 1IBSCRWFNNperiodical sinusoid4.90.3210Figure 15a
Case 2IBSCRWFNNperiodical sinusoid5.50.3720Figure 15c
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

Lin, F.-J.; Huang, M.-S.; Chien, Y.-C.; Chen, S.-G. Intelligent Backstepping Control of Permanent Magnet-Assisted Synchronous Reluctance Motor Position Servo Drive with Recurrent Wavelet Fuzzy Neural Network. Energies 2023, 16, 5389. https://doi.org/10.3390/en16145389

AMA Style

Lin F-J, Huang M-S, Chien Y-C, Chen S-G. Intelligent Backstepping Control of Permanent Magnet-Assisted Synchronous Reluctance Motor Position Servo Drive with Recurrent Wavelet Fuzzy Neural Network. Energies. 2023; 16(14):5389. https://doi.org/10.3390/en16145389

Chicago/Turabian Style

Lin, Faa-Jeng, Ming-Shi Huang, Yu-Chen Chien, and Shih-Gang Chen. 2023. "Intelligent Backstepping Control of Permanent Magnet-Assisted Synchronous Reluctance Motor Position Servo Drive with Recurrent Wavelet Fuzzy Neural Network" Energies 16, no. 14: 5389. https://doi.org/10.3390/en16145389

APA Style

Lin, F. -J., Huang, M. -S., Chien, Y. -C., & Chen, S. -G. (2023). Intelligent Backstepping Control of Permanent Magnet-Assisted Synchronous Reluctance Motor Position Servo Drive with Recurrent Wavelet Fuzzy Neural Network. Energies, 16(14), 5389. https://doi.org/10.3390/en16145389

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