Multi-Parameter Fuzzy-Based Neural Network Sensorless PMSM Iterative Learning Control Algorithm for Vibration Suppression of Ship Rim-Driven Thruster

: Aiming to reduce motor speed estimation and torque vibration present in the permanent magnet synchronous motors (PMSMs) of rim-driven thrusters (RDTs), a position-sensorless control algorithm using an adaptive second-order sliding mode observer (SMO) based on the super-twisting algorithm (STA) is proposed. In which the sliding mode coefficients can be adaptively tuned. Similarly, an iterative learning control (ILC) algorithm is presented to enhance the robustness of the velocity adjustment loop. By continuously learning and adjusting the difference between the actual speed and given speed of RDT motor through ILC algorithm, online compensation for the q-axis given current of RDT motor is achieved, thereby suppressing periodic speed fluctuations during motor running. Fuzzy neural network (FNN) training can be used to optimize the STA-SMO and ILC parameters of RDT control system, while improving speed tracking accuracy. Finally, simulation and experimental verifications have been conducted on the vector control system based on the conventional PI-STA and modified ILC-STA. The results show that the modified algorithm can effectively suppress the estimated speed and torque ripple of RDT motor, which greatly improves the speed tracking accuracy.


Introduction
Electric propulsion systems have undergone significant development in recent years due to their advanced control efficiency and performance.The rim-driven thruster (RDT) is a new type of low-carbon emission electric thruster, which integrates the propeller blade and the motor together by embedding the motor into the duct.This new type of electric thruster can eliminate the need for shaft transmission, mechanical seals, and oil circuit systems in traditional ship propulsion systems.It can save space for cabins on ships, improve motor efficiency and hydrodynamic efficiency, and reduce system installation complexity and manufacturing costs [1].For the development of this new type of electric thruster, the selection of the propulsion motor RDT and its control algorithm are the key technologies.
At present, alternating current (AC) induction motors, brushless direct current (BLDC) motors, permanent magnet synchronous motors (PMSM), and high-temperature superconducting motors are included for the selection of RDT motors.Considering the work At present, alternating current (AC) induction motors, brushless direct current (BLDC) motors, permanent magnet synchronous motors (PMSM), and high-temperature superconducting motors are included for the selection of RDT motors.Considering the work efficiency, manufacturing costs and processes performance of RDT, PMSM is currently considered the best choice for RDT propulsion motor.In order to accomplish the vector control of PMSM, it is necessary to obtain accurate rotor position information using position sensors.However, traditional mechanical position sensors (i.e., photoelectric encoders, resolvers, etc.) need to be mounted on the transmission shaft, which does not exist in RDT.They also have problems regarding large size, high costs, and low reliability, especially in specific engineering applications.RDT has a compact structure without the transmission shaft that can work in a harsh underwater environment, see Figure 1.Therefore, it is of great significance to develop anti-chattering position-sensorless control algorithms for RDT drive motors to realize high-performance control [2,3].The PMSM position-sensorless control algorithms mainly include direct calculation methods, back-EMF integration methods, model reference adaptive methods, and various observer methods [4][5][6][7][8].Of these, the sliding mode observer (SMO) can reduce the dependence on motor parameters with advantages of having a fast response speed, strong robustness, and simple structure.Therefore, it is widely used in the position-sensorless control of PMSM [9].The chattering will occur in the control system when using the SMO algorithm due to the switching of the sliding mode surface [10].Especially in the vector control of PMSM at low speeds, it causes relatively large torque ripples.Therefore, it is necessary to eliminate chattering while maintaining control accuracy when using the SMO algorithm.The current research mainly focuses on the comprehensive improvement of switching function, sliding mode gain, estimation accuracy, etc. [11][12][13][14].The first-order SMO always has chattering problems, which have a negative impact on the stability of the system.The secondorder SMO can effectively eliminate chattering while maintaining the system robustness [15,16] and is well received by many scholars.
Due to the influences of motor design, inverter nonlinearities, and current detection errors, the PMSM control system suffers from torque vibration with current harmonics, which also reduces speed control accuracy.Usually, the vibration of motor torque is suppressed by optimizing motor design or improving the motor control strategy.However, it is difficult to modify the design after completing the original motor design [17].With the development of control theory, many improved control strategies have been proposed to reduce torque vibration suppression, such as harmonic voltage injection, proportionalintegral quasi-resonant compensator, and dead-zone feedforward compensation [18,19].Iterative learning control (ILC) is an approach to improve the performance of a system that executes repetitively and periodically over a fixed time interval by learning from previous executions (trials, iterations, and passes) [20].It does not rely on the mathematical model of systems and can be used in PMSM control to effectively improve tracking performance and compensate the control signal through continuous learning [21].Of these, the sliding mode observer (SMO) can reduce the dependence on motor parameters with advantages of having a fast response speed, strong robustness, and simple structure.Therefore, it is widely used in the position-sensorless control of PMSM [9].The chattering will occur in the control system when using the SMO algorithm due to the switching of the sliding mode surface [10].Especially in the vector control of PMSM at low speeds, it causes relatively large torque ripples.Therefore, it is necessary to eliminate chattering while maintaining control accuracy when using the SMO algorithm.The current research mainly focuses on the comprehensive improvement of switching function, sliding mode gain, estimation accuracy, etc. [11][12][13][14].The first-order SMO always has chattering problems, which have a negative impact on the stability of the system.The second-order SMO can effectively eliminate chattering while maintaining the system robustness [15,16] and is well received by many scholars.
Due to the influences of motor design, inverter nonlinearities, and current detection errors, the PMSM control system suffers from torque vibration with current harmonics, which also reduces speed control accuracy.Usually, the vibration of motor torque is suppressed by optimizing motor design or improving the motor control strategy.However, it is difficult to modify the design after completing the original motor design [17].With the development of control theory, many improved control strategies have been proposed to reduce torque vibration suppression, such as harmonic voltage injection, proportionalintegral quasi-resonant compensator, and dead-zone feedforward compensation [18,19].Iterative learning control (ILC) is an approach to improve the performance of a system that executes repetitively and periodically over a fixed time interval by learning from previous executions (trials, iterations, and passes) [20].It does not rely on the mathematical model of systems and can be used in PMSM control to effectively improve tracking performance and compensate the control signal through continuous learning [21].
This paper presents a RDT position-sensorless control algorithm using an adaptive second-order SMO based on the super-twisting algorithm (STA).An adaptive algorithm is used to tune the sliding mode coefficients.Furthermore, a torque ripple suppression method based on ILC algorithm is proposed.A closed-loop PI-ILC algorithm with a forgetting factor is introduced to compensate for the q-axis given current through iterative learning, thereby reducing the real-time speed error.Considering the error and error rate of motor speed as input, and fuzzy neural network (FNN) training is used to optimize the ILC and STA-SMO parameters to improve the accuracy and stability of the control system.The proposed control algorithm can suppress the motor speed torque vibration and improve the speed tracking accuracy.The main innovations of this paper are as follows: (1) A RDT control strategy using an adaptive second-order SMO based on the STA is designed.
(2) To reduce the torque ripple suppression of RDT motor, a closed-loop PI-ILC algorithm is introduced to compensate the speed real-time error to the q-axis given current.(3) FNN training is used to optimize the ILC and STA-SMO parameters to improve the accuracy and stability of the current loop.
The rest of this paper is organized as follows: Section 2 introduces the materials and methods of RDT motor model and proposed control algorithm.Section 3 introduces the application of multi-parameter FNN for parameter tuning.Section 4 introduces the simulation results and discussions of a RDT motor sensorless control model on MATLAB R2017b/Simulink environment.Section 5 introduces the experimental results and the discussions of a RDT motor sensorless control system on experimental setup.The final section presents the conclusions and subsequent ideas of the topic covered in this paper.

Materials and Methods
The mathematic current model of the RDT motor in the αβ-axis reference system can be expressed as: where i α and i β and u α and u β represent the current and voltage in αβ-axis reference system, respectively.L s and R s represent the inductance and resistance of stator winding.e α and e β represent the back EMFs in αβ-axis reference system and can be expressed as: where ψ f is the flux linkages and ω and θ represent the rotor speed and position of the RDT motor.

Adaptive STA-SMO
The super-twisting algorithm (STA) was first proposed by A. Levant to effectively eliminate the chattering problem of first-order SMO.The basic form of STA with perturbation was designed as [22]: where x 1 , x 2 , and x 1 are the state variables, error between estimated values, and actual values of the state variables.k 1 and k 2 are designed to be the sliding-mode coefficients of STA.z 1 and z 2 are perturbation terms of system.
The stable conditions of STA have been proved in [23,24], demonstrating that the perturbation terms should be globally bounded by: And the sliding mode coefficients k 1 and k 2 should satisfy (5) to stabilize the system: where δ 1 > 0 and is a normal number.
In order to obtain the estimated value of the back EMF, the stator currents of RDT motor are estimated using STA based second-order SMO [15,16]: in which îα and îβ denote the estimated values of i α and i β .By comparing Equation (3) with Equation ( 6), it can be concluded that perturbation terms z 1 and z 2 are designed through: By substituting Equation (7) into Equation (4), we obtain: For a sufficiently large δ 1 , the validity of the above inequality can be met.The state equations of currents errors in αβ-axis reference system can be derived by subtracting Equation (1) from Equation (6): When the system reaches the sliding surface and the estimated error is close to 0 ( i α = 0, i β = 0).The equivalent back-EMFs e α and e β can be obtained by the equivalent control method, it can be concluded that: (10) where k 1 and k 2 are the sliding mode coefficients.It is difficult to determine the actual values of k 1 and k 2 in practice.And different sets of k 1 and k 2 will influence the stability conditions of STA-SMO.Therefore, it is crucial to adaptively adjust the sliding mode coefficients.The adaptive algorithm is given as: (11) in which σ 1 and σ 2 are adaptive coefficients.The stability of the adaptive algorithm has been proved in [16].The rotor position estimation method based on the arctangent function will directly introduce the sliding mode chattering into the division operation, resulting in significant estimation errors.Therefore, the PLL structure is used to extract the rotor position information of the RDT motor.It can suppress chattering and maintain good control performance.The structure block diagram of the adaptive STA-SMO is shown in Figure 2.
The rotor position estimation method based on the arctangent function will directly introduce the sliding mode chattering into the division operation, resulting in significant estimation errors.Therefore, the PLL structure is used to extract the rotor position information of the RDT motor.It can suppress chattering and maintain good control performance.The structure block diagram of the adaptive STA-SMO is shown in Figure 2.

Design of ILC Algorithm
By analyzing the previous control data, ILC can find appropriate control input through real-time online iteration, theoretically obtaining accurate tracking trajectories.
Traditional PI control achieves satisfactory results by selecting appropriate proportional and integral parameters but struggles to meet high-precision control requirements.By introducing the ILC algorithm to form an ILC controller, the q-axis given current of motor can be obtained with ILC controller and traditional PI controller.
The type of PI-ILC algorithm is adopted to acquire the learning law [19]: where , 1 ( ) is the reference signal of the current cycle's q-axis current; , ( ) is the reference signal of the q-axis current in the previous cycle called the "learning experience"; is the current cycle speed deviation signal, which is mainly used to compensate for

Design of ILC Algorithm
By analyzing the previous control data, ILC can find appropriate control input through real-time online iteration, theoretically obtaining accurate tracking trajectories.
Traditional PI control achieves satisfactory results by selecting appropriate proportional and integral parameters but struggles to meet high-precision control requirements.By introducing the ILC algorithm to form an ILC controller, the q-axis given current of motor can be obtained with ILC controller and traditional PI controller.
The type of PI-ILC algorithm is adopted to acquire the learning law [19]: where i q,k+1 * (t) is the reference signal of the current cycle's q-axis current; i q,k * (t) is the reference signal of the q-axis current in the previous cycle called the "learning experience"; e k+1 (t) is the current cycle speed deviation signal, which is mainly used to compensate for the "learning experience" loss caused by the forgetting factor; and η and ξ are the closed-loop learning coefficients.Formula (12) in the z-domain is shown as: The iterative process can be regarded as a delay link, which is shown as: in the above equation, N = f sample / f signal , f sample represents the sampling frequency of the digital control system, while f signal is the system signal frequency.
When N takes a value of 1, the z-domain transfer function of the equation is shown as: The structure block diagram of ILC process and closed-loop PI-ILC are shown in Figures 3 and 4  the "learning experience" loss caused by the forgetting factor; and η and ξ are the closed-loop learning coefficients.Formula (12) in the z-domain is shown as: , The iterative process can be regarded as a delay link, which is shown as: , 1 , ( ) ( ) in the above equation, , sample f represents the sampling frequency of the digital control system, while signal f is the system signal frequency.
When N takes a value of 1, the z-domain transfer function of the equation is shown as: The RDT STA-SMO (1 ) ( )   (1 ) ( ) The design of the ILC controller should not only ensure its stability but also ensure that the iterative algorithm converges to its true value.The convergence conditions of the PI type iterative algorithm are determined as follows.
Combining Formulas ( 12) and ( 15) through , it can be seen that: The design of the ILC controller should not only ensure its stability but also ensure that the iterative algorithm converges to its true value.The convergence conditions of the PI type iterative algorithm are determined as follows.
Combining Formulas ( 12) and ( 15) through e k+1 = ω re f − ω r,k+1 , it can be seen that: where J is the inertia moment, n p is the pole pair of the motor, and k t is the torque coefficient.Therefore, the following formulas are obtained: The definition of the infinite parametric number and the absolute value inequality leads to: When lim k→∞ e k (t) = 0, the algorithm converges by considering the convergence condition of the PI-ILC algorithm: where t ∈ [0, T], T is the system tracking time.

Application of Multi-Parameter Fuzzy-Based Neural Network for Parameter Auto-Tuning
The principle of neural network is to analysis errors based on training and online results, modify weights and thresholds, and iteratively obtain models with the online results as outputs.It consists of input and output layer nodes, as well as one or more hidden layer nodes.The input information must be propagated forward to the hidden layer, and the activated information is passed to each output node after the activation function operation of the unit node to obtain the output results.
Fuzzy neural networks (FNNs) have good self-learning and adaptive abilities with strong robustness that only relies on fuzzy rules [25][26][27][28][29].They can perform logical reasoning through established fuzzy rules which have good approximation abilities for certain complex nonlinear systems.
The proposed algorithm consists of an ILC algorithm, adaptive STA-SMO algorithm, and FNN algorithm.FNNs can generate fuzzy rules and adjust affiliation functions through strong self-learning and adaptive training.The control system transfers the speed deviation e(t) and the deviation rate ec(t) to the FNN, which can adjust the ILC parameters α, η, and ξ and the STA-SMO parameters σ 1 and σ 2 to achieve better speed control performance after the fuzzification, fuzzy inference, and anti-fuzzification of the output.
The structure block diagram of FNN is shown as Figure 5.The four layers are input, fuzzification, fuzzy inference, and output layers in order.The calculation formula for input layer is shown as follows: The calculation formula for the fuzzy inference layer is shown as: The calculation formula for the output layer is shown as follows: in which ij ω is the weight coefficient.
The ILC and STA-SMO control parameters output by the output layer are shown as: (1) , ( 2) , ( 3) , ( 4) , ( 5) The parameters to be learned by the FNN are the affiliation function center value ij c , the width value ij b , and the weight value ij ω , and the objective function is shown as follows: in the formula, ( ) r k , ( ) y k , and ( ) ( ) r k y k − are the expected output, actual output of each network learning, and the control error.
The learning algorithm for the centroid, width, and weight of the affiliation function is shown as: The fuzzification layer considers a Gaussian membership function as the analysis method, which is shown as: in which c ij is the central value of membership function; b ij represents the width of function curve; and n is the number of fuzzy subsets.The calculation formula for the fuzzy inference layer is shown as: The calculation formula for the output layer is shown as follows: in which ω ij is the weight coefficient.
The ILC and STA-SMO control parameters output by the output layer are shown as: The parameters to be learned by the FNN are the affiliation function center value c ij , the width value b ij , and the weight value ω ij , and the objective function is shown as follows: in the formula, r(k), y(k), and r(k) − y(k) are the expected output, actual output of each network learning, and the control error.
The learning algorithm for the centroid, width, and weight of the affiliation function is shown as: where k is the iterative steps of network, λ represents the learning rate, and µ represents the inertia coefficient.
The block diagram of FNN and control system is presented in Figure 6: ( where k is the iterative steps of network, λ represents the learning rate, and μ repre- sents the inertia coefficient.
The block diagram of FNN and control system is presented in Figure 6: (1) Determine neural network framework through number of layers and nodes in each network layer and select initial weights, learning rate, and momentum factor.

Simulation Results and Discussions
The simulation model in this paper is based on MATLAB R2017b/Simulink environment.The structure block diagram of sensorless control model for RDT motor is shown in Figure 7.The main parameters of the RDT motor in simulation are summarized in Table 1.

Simulation Results and Discussions
The simulation model in this paper is based on MATLAB R2017b/Simulink environment.The structure block diagram of sensorless control model for RDT motor is shown in Figure 7.The main parameters of the RDT motor in simulation are summarized in Table 1.0.12 mH To verify the system control performance based on the proposed algorithm, comparisons between conventional PI-STA and modified ILC-STA have been carried out through simulations.The motor starts with the initial given speed of 1200 r/min, the load torque changes to 1.5 N•m at 0.4 s, the given speed changes to 800 r/min at 0.6 s, and the load torque changes back to 0 at 0.9 s.
The simulation results of the RDT motor using conventional PI-STA algorithm are shown in Figure 8. From top to bottom, the estimated speed and actual speed, the speed error, the given torque and electromagnetic torque, the estimated and actual position, and the position error and the three-phase stator currents are presented.It can be seen that the speed error is about −3 and 22 r/min under the 1200 r/min given speed and about −1 and 12 r/min under the 800 r/min given speed.The torque errors are about −1.2 and 1.1 Nm, −1.1 and 1.2 Nm, −1.3 and 1.2 Nm, and −0.9 and 0.8 Nm, respectively.The position error of the RDT motor is about 0.7 rad.To verify the system control performance based on the proposed algorithm, comparisons between conventional PI-STA and modified ILC-STA have been carried out through simulations.The motor starts with the initial given speed of 1200 r/min, the load torque changes to 1.5 N•m at 0.4 s, the given speed changes to 800 r/min at 0.6 s, and the load torque changes back to 0 at 0.9 s.
The simulation results of the RDT motor using conventional PI-STA algorithm are shown in Figure 8. From top to bottom, the estimated speed and actual speed, the speed error, the given torque and electromagnetic torque, the estimated and actual position, and the position error and the three-phase stator currents are presented.It can be seen that the speed error is about −3 and 22 r/min under the 1200 r/min given speed and about −1 and 12 r/min under the 800 r/min given speed.The torque errors are about −1.2 and 1.1 Nm, −1.1 and 1.2 Nm, −1.3 and 1.2 Nm, and −0.9 and 0.8 Nm, respectively.The position error of the RDT motor is about 0.7 rad.
The simulation results of the RDT motor using modified ILC-STA are shown in Figure 9.It can be seen that the speed error is about −2 and 11 r/min under the 1200 r/min given speed and about −1 and 7 r/min under the 800 r/min given speed.The torque errors are about −0.4 and 0.4 Nm, −0.6 and 0.6 Nm, −0.5 and 0.5 Nm, and −0.5 and 0.5 Nm, respectively.The estimated position error of the RDT motor is about 0.7 rad.
From the simulation results, it can be seen that compared to the conventional PI-STA, the modified ILC-STA produces smaller speed and torque errors.The estimated position error is almost the same.The RDT motor speed, torque, and stator currents results display less chattering.This proves that the modified algorithm obviously maintains better control performance.Table 2 shows the performance comparison of different control algorithms.From the simulation results, it can be seen that compared to the conventional P the modified ILC-STA produces smaller speed and torque errors.The estimated po error is almost the same.The RDT motor speed, torque, and stator currents results d less chattering.This proves that the modified algorithm obviously maintains better c performance.Table 2 shows the performance comparison of different control algori The simulation results of the RDT motor using modified ILC-STA are shown i ure 9.It can be seen that the speed error is about −2 and 11 r/min under the 1200 given speed and about −1 and 7 r/min under the 800 r/min given speed.The torque are about −0.4 and 0.4 Nm, −0.6 and 0.6 Nm, −0.5 and 0.5 Nm, and −0.5 and 0.5 N spectively.The estimated position error of the RDT motor is about 0.7 rad.From the simulation results, it can be seen that compared to the conventional P the modified ILC-STA produces smaller speed and torque errors.The estimated po error is almost the same.The RDT motor speed, torque, and stator currents results d less chattering.This proves that the modified algorithm obviously maintains better c performance.Table 2 shows the performance comparison of different control algori

Experimental Results and Discussions
In order to verify the performance of the RDT motor sensorless control algorithm, a TMS320F28335 DSP-based experimental setup was built, as shown in Figure 10.The experimental setup consists of a magnetic powder brake, a torque sensor, and a test RDT motor, which are fixed together on a T-slotted platform and connected by elastic couplings.
The main parameters of the test RDT motor are the same as that in the simulation.And the experimental results under different conditions are shown in Figures 11-13.

Experimental Results and Discussions
In order to verify the performance of the RDT motor sensorless control algorithm, a TMS320F28335 DSP-based experimental setup was built, as shown in Figure 10.The experimental setup consists of a magnetic powder brake, a torque sensor, and a test RDT motor, which are fixed together on a T-slotted platform and connected by elastic couplings.The main parameters of the test RDT motor are the same as that in the simulation.And the experimental results under different conditions are shown in Figures 11-13.Condition 1: The tested given speed is set to 1200 r/min with sudden load torque changes from 0 to 1.5 Nm. Figure 11 shows the experimental results comparison between the conventional PI-STA and the modified ILC-STA in each case of step change in load torque.From top to bottom, the estimated speed and actual speed, the speed error, the given torque and electromagnetic torque, the estimated and actual position, the position error, and the three-phase stator currents are plotted.It can be observed that the stable error of the modified algorithm is smaller than that of the conventional algorithm.The motor speed estimation error ripple is about −2 and 22 r/min, respectively.The motor speed error ripple is about 24 r/min, and the electromagnetic torque ripples are about 2.2 Nm and 3 Nm, respectively.The rotor position error between the estimate and actual value of RDT motor is about 0.7 rad.When the modified ILC-STA is used, the motor speed estimation errors are about −2 and 11 r/min, respectively.The motor speed error ripple is about 13 r/min, and the electromagnetic torque ripples are about 0.9 Nm and 1 Nm, respectively.The rotor position error between the estimate and actual value of RDT motor is also about 0.7 rad.From Figure 11, it can be seen that the motor speed and electromagnetic torque waveforms are more stable, and the ripple is smaller than when using the conventional algorithm.The RDT motor speed error ripple deceases by 46% and the torque ripple deceases by 59% and 67%, respectively.The maximum speed error deceases by 11 r/min.Condition 1: The tested given speed is set to 1200 r/min with sudden load torque changes from 0 to 1.5 Nm. Figure 11 shows the experimental results comparison between the conventional PI-STA and the modified ILC-STA in each case of step change in load torque.From top to bottom, the estimated speed and actual speed, the speed error, the given torque and electromagnetic torque, the estimated and actual position, the position error, and the three-phase stator currents are plotted.It can be observed that the stable error of the modified algorithm is smaller than that of the conventional algorithm.The motor speed estimation error ripple is about −2 and 22 r/min, respectively.The motor Condition 2: The tested given speed decreases from 1200 r/min to 800 r/min, with a continuous load torque of 1.5 Nm. Figure 12 shows the experimental results comparison between conventional PI-STA and modified ILC-STA in case of a step change in given speed.It is obvious that the stable error of the modified algorithm is smaller than that of the conventional algorithm.The motor speed estimation error ripples are about −1 and 12 r/min, respectively.The motor speed error ripple is about 13 r/min, and the electromagnetic torque ripple is about 2.5 Nm.When the modified ILC-STA algorithm is used, the motor speed estimation errors are about −1 and 7 r/min, respectively.The motor speed error ripple is about 8 r/min, and the electromagnetic torque ripple is about 1 Nm.From Figure 12, it can be seen that the motor actual speed and electromagnetic torque waveforms are more stable and the ripple is smaller than when using conventional algorithm.The RDT motor speed error ripple and torque ripple decease by 38% and 60%, respectively.The maximum speed error deceases by 5 r/min.is also about 0.7 rad.From Figure 11, it can be seen that the motor speed and electromagnetic torque waveforms are more stable, and the ripple is smaller than when using the conventional algorithm.The RDT motor speed error ripple deceases by 46% and the torque ripple deceases by 59% and 67%, respectively.The maximum speed error deceases by 11 r/min.Condition 2: The tested given speed decreases from 1200 r/min to 800 r/min, with a continuous load torque of 1.5 Nm. Figure 12 shows the experimental results comparison between conventional PI-STA and modified ILC-STA in case of a step change in given speed.It is obvious that the stable error of the modified algorithm is smaller than that of the conventional algorithm.The motor speed estimation error ripples are about −1 and 12 r/min, respectively.The motor speed error ripple is about 13 r/min, and the electromagnetic torque ripple is about 2.5 Nm.When the modified ILC-STA algorithm is used, the error ripple is about 8 r/min, and the electromagnetic torque ripple is about 1 Nm.From Figure 12, it can be seen that the motor actual speed and electromagnetic torque waveforms are more stable and the ripple is smaller than when using conventional algorithm.The RDT motor speed error ripple and torque ripple decease by 38% and 60%, respectively.The maximum speed error deceases by 5 r/min.Condition 3: The tested given speed is set as the speed of 800 r/min with sudden load torque changes from 1.5 Nm to 0. Figure 13 shows the experimental results comparison between conventional PI-STA and modified ILC-STA in the case of step change in load torque.Similarly to the previous conditions, it is obvious that the stable error of the modified algorithm is smaller than that of the conventional algorithm.The motor speed estimation error ripples are about −1 and 12 r/min, respectively.The motor speed error ripple is about 13 r/min, and the electromagnetic torque ripple is about 1.7 Nm.When the modified ILC-STA algorithm is used, the motor speed estimation errors are about −1 and 7 Condition 3: The tested given speed is set as the speed of 800 r/min with sudden load torque changes from 1.5 Nm to 0. Figure 13 shows the experimental results comparison between conventional PI-STA and modified ILC-STA in the case of step change in load torque.Similarly to the previous conditions, it is obvious that the stable error of the modified algorithm is smaller than that of the conventional algorithm.The motor speed estimation error ripples are about −1 and 12 r/min, respectively.The motor speed error ripple is about 13 r/min, and the electromagnetic torque ripple is about 1.7 Nm.When the modified ILC-STA algorithm is used, the motor speed estimation errors are about −1 and 7 r/min, respectively.The motor speed error ripple is about 8 r/min, and the electromagnetic torque ripple is about 1 Nm.Similarly to the previous conditions, it is obvious that the motor speed and electromagnetic torque waveform are more stable and the ripple is smaller than when using conventional method.The speed error ripple and torque ripple decease by 38% and 41%, respectively.The maximum speed error deceases by 5 r/min.

Host computer
Similarly to the simulation results, it can be seen from the experimental results that compared to the conventional PI-STA, the modified ILC-STA produces lower speed errors and smaller electromagnetic torque errors.The motor speed, electromagnetic torque, and stator currents results also have less chattering.Table 3 summarizes the performance comparison of the different control algorithms.

Conclusions
In this paper, a RDT motor position-sensorless control algorithm that uses an adaptive second-order SMO based on STA is presented to reduce system chattering, and a closedloop PI-ILC algorithm with the forgetting factor is used to suppress torque vibration.In addition, FNN training is used to optimize the ILC and STA-SMO parameters.The simulation and experimental analysis have been used to verify the proposed algorithm.The presented RDT sensorless control algorithm can be used to reduce the chattering generated in the RDT motor speed and torque ripple with improved speed tracking accuracy.The accuracy and complexity have been balanced and can be used to explore future research.
Our future work will apply the presented RDT sensorless control algorithm to a real ship, which will greatly contribute to the realization of smart and low carbon emissionoriented maritime traffic control technology.

Figure 1 .
Figure 1.Physical diagram and geometric structure of an RDT.

Figure 1 .
Figure 1.Physical diagram and geometric structure of an RDT.

Figure 2 .
Figure 2. Structure block diagram of the adaptive STA-SMO.
J. Mar.Sci.Eng.2024, 12, x FOR PEER REVIEW 6 of 18 structure block diagram of ILC process and closed-loop PI-ILC are shown in Figures 3 and 4.

Figure 3 .
Figure 3. Structure block diagram of the ILC process.

Figure 3 .
Figure 3. Structure block diagram of the ILC process.

Figure 4 .
Figure 4. Structure block diagram of the closed-loop PI-ILC.

Figure 4 .
Figure 4. Structure block diagram of the closed-loop PI-ILC.

Figure 6 .
Figure 6.Block diagram of fuzzy network ILC control system.

Figure 7 .
Figure 7. Structure block diagram of sensorless control for RDT motor.

Figure 7 .
Figure 7. Structure block diagram of sensorless control for RDT motor.

Figure 8 .Figure 9 .
Figure 8. Simulation results using conventional PI-STA.The simulation results of the RDT motor using modified ILC-STA are shown ure 9.It can be seen that the speed error is about −2 and 11 r/min under the 1200 given speed and about −1 and 7 r/min under the 800 r/min given speed.The torque are about −0.4 and 0.4 Nm, −0.6 and 0.6 Nm, −0.5 and 0.5 Nm, and −0.5 and 0.5 N spectively.The estimated position error of the RDT motor is about 0.7 rad.

Figure 10 .
Figure 10.Experimental setup for the RDT motor.(a) Block diagram of experimental setup.(b) Diagram of physical experimental setup.

Figure 10 .
Figure 10.Experimental setup for the RDT motor.(a) Block diagram of experimental setup.(b) Diagram of physical experimental setup.

JFigure 11 .
Figure 11.Experimental results under a step change of the load torque condition.(a) Conventional PI-STA.(b) Modified ILC-STA.

Figure 11 .
Figure 11.Experimental results under a step change of the load torque condition.(a) Conventional PI-STA.(b) Modified ILC-STA.

Figure 12 .
Figure 12.Experimental results under a step change of a given speed condition.(a) Conventional PI-STA.(b) Modified ILC-STA.

Figure 12 .
Figure 12.Experimental results under a step change of a given speed condition.(a) Conventional PI-STA.(b) Modified ILC-STA.

Figure 13 .
Figure 13.Experimental results under a step change of the load torque condition.(a) Conventional PI-STA.(b) Modified ILC-STA.

Figure 13 .
Figure 13.Experimental results under a step change of the load torque condition.(a) Conventional PI-STA.(b) Modified ILC-STA.
1) Determine neural network framework through number of layers and nodes in each network layer and select initial weights, learning rate, and momentum factor.(2) Calculate the objective function obtained by running k times.(3) Calculate the input and output values of each network layer.(4) Calculate the output of the self-disturbance rejection controller.(5) Perform backpropagation calculations to correct the weights of output and hidden layers.(6) If successful, end iteration; if not, return to step (2).

Table 1 .
Main parameters of RDT motor.

Table 1 .
Main parameters of RDT motor.

Table 2 .
Simulation performance comparison of different control algorithms.

Table 3 .
Experimental performance comparison of different control algorithms.