Next Article in Journal
Propagating Particle Tracking Uncertainty Defined by Fuzzy Numbers in Spatially Variable Velocity Fields
Next Article in Special Issue
Heat Integration of Liquid Hydrogen-Fueled Hybrid Electric Ship Propulsion System
Previous Article in Journal
Insights on the Extreme Storm Surge Event of the 22 November 2022 in the Venice Lagoon
Previous Article in Special Issue
Improvement and Optimization Configuration of Inland Ship Power and Propulsion System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Active Disturbance Rejection Control Method for Marine Permanent-Magnet Propulsion Motor Based on Improved ESO and Nonlinear Switching Function

1
School of Marine Engineering, Dalian Maritime University, Dalian 112026, China
2
School of Navigation and Naval Architecture, Dalian Ocean University, Dalian 112026, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2023, 11(9), 1751; https://doi.org/10.3390/jmse11091751
Submission received: 6 August 2023 / Revised: 26 August 2023 / Accepted: 28 August 2023 / Published: 7 September 2023
(This article belongs to the Special Issue New Advances on Energy and Propulsion Systems for Ship)

Abstract

:
In the control of marine permanent-magnet propulsion motors, active disturbance rejection control has attracted much attention because it can deal with external load disturbances and uncertainties of motor parameters at the same time. However, the conventional second-order ADRC has the problem of slow disturbance observation speed. To this end, this paper proposes an improved third-order extended state observer using the proportional–integral disturbance update law to improve the tracking performance and anti-external disturbance ability of the motor control system. Then, aiming at the problem that the structure does not effectively use the current information, resulting in large speed fluctuations when the load changes, the measured value of the q-axis current is used as the disturbance feedforward compensation item to further improve the load disturbance suppression ability of the motor. Finally, in order to suppress the influence of the current periodic disturbance caused by unmodeled dynamics on the steady-state accuracy of the motor, a nonlinear switching function with bounded gain and an IIR low-pass filter are designed to suppress the periodic disturbance without affecting the dynamic performance of the system. Combined with the established ship propeller load model, the effectiveness of the method is verified on the motor experimental platform: When suddenly changing 75% of the propeller load, the motor speed decreases by about 20%, and the adjustment time is 0.1 s, which improves the performance by more than 70% compared to PI control and conventional ADRC methods.

1. Introduction

Under the background of carbon peaking and carbon neutrality goals, the issue of ship emission reduction and energy efficiency optimization is becoming more and more important. Among the different means of transportation, global maritime transport is responsible for 2–3% of global greenhouse gas (GHG) emissions, and it is predicted to increase to 17% by 2050 if no changes are adapted. Hence, the international maritime organization (IMO) has aimed to reach a 50% reduction in GHG emissions by 2050 compared to 2008. Hence, to meet these strict GHG emissions reduction targets, ships must use alternative power sources [1,2]. These strict regulations have brought new challenges to ship designers and ship owners, and have had a significant impact on the energy efficiency design indicators and energy efficiency operation indicators of ships. It is worth noting that by the end of the 20th century, with the rapid development of fast-switching semiconductor devices and the higher utilization of power converters, ship propulsion systems began to explore electrification [2,3]. In this stage, generators are used together with internal combustion engines, and the electrical energy obtained from the generator is used to run the propulsion motors. Therefore, regardless of the speed of internal combustion engine, the speed of propeller can be freely adjusted, thereby improving the overall energy efficiency [3]. Now, under the strict policy background, to further reduce ship emissions and achieve a more flexible energy distribution, full-electric propulsion systems using lithium batteries and fuel cells have begun to emerge; so, electric propulsion technology will be more advantageous. Therefore, electric propulsion ships are becoming the trend of future ships with their advantages of low emission, high efficiency and low noise. As the power source of an electric propulsion ship, the control performance and efficiency of the motor and its variable frequency drive system determine the maneuverability and energy efficiency of the ship; so, it is necessary to study its control strategy in depth.
Permanent-magnet synchronous motors (PMSMs) are widely used in propulsion devices for large ships and marine equipment due to their advantages of high efficiency and high power density. However, the permanent-magnet synchronous motor is a typical nonlinear and strongly coupled system [4,5], which is difficult to described with an accurate mathematical model. The permanent-magnet synchronous motor is susceptible to various load changes, external sea state changes and parameter changes, and the long-term operation of the motor system in a high-salt, high-temperature, and high-humidity marine environment will lead to problems such as the reduction of sensor accuracy [6,7]. The traditional PI control method cannot meet the requirements of high-performance control of marine propulsion motors; so, it is necessary to use nonlinear controllers to improve their tracking and anti-disturbance capabilities.
In recent years, relevant researchers have proposed a variety of control methods and applied them to the motor control, such as Active Disturbance Rejection Control (ADRC) [8,9,10] Model Predictive Control (MPC) [11,12,13], sliding mode control (SMC) [14], etc. Ref. [8] proposed the nonlinear ADRC controller, but due to the complexity of parameter setting and stability analysis of the nonlinear ADRC controller, this method was not widely used until [9] proposed a linear ADRC parameter-tuning method based on observer bandwidth. Ref. [10] used ADRC in permanent-magnet synchronous motor control and showed that the method had better anti-disturbance capabilities and parameter robustness. Extended State Observer (ESO) is the core link of ADRC, which can observe and estimate the system state and lumped disturbances. Ref. [13] adopted the extended state observer to solve the problem of steady-state error and poor robustness of parameters in the motor model predictive control. Ref. [14] combined the sliding mode control with the extended state observer and updated the permanent-magnet synchronous motor speed sliding mode control law in real time by using the lumped disturbances estimated by the observer to improve the speed tracking performance and reduce the steady-state speed tracking error. These control methods start from different practical problems, and use ESO’s effective state observation methods that do not depend on accurate models to improve the dynamic response performance, harmonic suppression ability, and stability and robustness of the system. These studies also show that ADRC has good disturbance suppression performance and parameter robustness, and is suitable for ship electric propulsion working conditions. For electric propulsion ships, permanent-magnet synchronous motors are usually used as propulsion motors to drive the propellers. During navigation, they often face complex ocean environments and disturbances such as wind, waves, and currents. At this time, the torque load of permanent-magnet synchronous motors is constantly changing, and a set of PI parameters is very difficult to adapt well to all working conditions. Hence, it is necessary to consider a more complete control strategy. If the changing load and friction are regarded as disturbances, the control problem of time-varying systems is equivalent to the suppression of disturbances. The problem to be solved by ADRC is the estimation and compensation of disturbances. In addition, there are unmodeled dynamics (such as Flux Harmonics, Dead-Time Effects, and Measurement Error Effects) in the motor itself [4], which cause torque and speed to fluctuate. Therefore, the control performance of disturbances and unmodeled dynamics determines the propulsion efficiency of the ship, and it is necessary to study the disturbance suppression strategy of the ship propulsion motor in depth.
Although ADRC has good control performance, it can still be improved for different control problems. To further improve the anti-disturbance capabilities of ADRC, researchers have proposed many improvement methods. First of all, although the linear ADRC parameter setting is simple and easy to use, and the suppression effect is better when there is a large disturbance, the nonlinear active disturbance rejection controller can achieve more accurate state estimation, and has a better suppression effect on small disturbances. To solve this problem, relevant scholars have proposed a series of linear/nonlinear active disturbance rejection switching control method to improve the system’s anti-disturbance capabilities. Ref. [15] analyzed the effective operating range and steady-state observation error of the linear/nonlinear ADRC controller, proposed a linear/nonlinear switching control method, and explained its necessity, but the direct switching method may cause current oscillation in the motor control. Therefore, Ref. [16] proposed a linear/nonlinear switching active disturbance rejection control method applied to permanent-magnet synchronous motor speed regulation and designed a hysteretic switching strategy, which achieved better anti-disturbance capabilities. Ref. [17] improved the nonlinear function and proposed a class of linear/nonlinear switching active disturbance rejection controllers. By designing the switching transition strategy, the smooth switching of the control variable under different disturbance amplitudes and the improvement of the anti-disturbance performance were achieved. However, the parameter design and stability analysis of the nonlinear method are complicated. Ref. [18] designed a new nonlinear function to combine the advantages of linear and nonlinear ADRC and apply it to rotor position estimation, which improved the system’s anti-disturbance capabilities and robustness. Ref. [19] adopted the observation disturbance optimization controller switching strategy and used it in the propeller propulsion motor control of the aircraft to enhances its speed response and anti-disturbance performances. Refs. [16,17,18,19] suggestions are all in the motor system, by switching the linear/nonlinear ADRC controller to improve the anti-disturbance capabilities, but this kind of method mainly combines the advantages of the linear/nonlinear ADRC controller; the parameter configuration and the stability analysis are still relatively complicated, and the control performance of the linear/nonlinear active disturbance rejection controller itself has not been improved in essence.
In terms of the structure of ADRC itself, some scholars have proposed some improvement methods for the problem of poor disturbance suppression performance in the middle frequency band of the controller, which suppresses the problems of speed, torque ripples, and motor noise caused by periodic disturbances in the motor. Ref. [20] proposed a control method that combines an improved vector resonant controller and an active disturbance rejection control controller for suppressing the current harmonics of permanent-magnet synchronous linear motors. Ref. [21] combined the quasi-resonant controller with the active disturbance rejection control law to suppress the first and second harmonics of the motor speed caused by current measurement errors, but this method is prone to stability problems under high-speed conditions. Ref. [22] designed an extended harmonic state observer (EHSO) for the estimation and attenuation of the selective periodic disturbance, which addresses the torque-ripple reduction in PMSM drives for smooth speed control. Moreover, a pole-placement strategy is proposed and analyzed through sensitivity function to improve the disturbance-rejection capability and the relative stability. Ref. [23] proposed a two-degree-of-freedom ADRC current-control method based on the improved extended state observer to simultaneously improve the current control’s dynamic response speed and steady-state control accuracy, eliminate the influence of ADRC tracking performance and disturbance suppression coupling on harmonic suppression. And, it is worth noting that the observed disturbance update law of the improved ESO is designed as a proportional–integral–repetitive control structure. Ref. [24] proposed an adaptive ADRC for the current loop of PMSM to suppress the uncertain periodic and aperiodic disturbances simultaneously, and it does not require the disturbance frequency information.
In summary, ADRC is a control method with strong anti-disturbance ability and good parameter robustness, which is suitable for the working conditions of ship electric propulsion [25]. However, for time-varying disturbances, the ADRC controller also has the problem of poor anti-disturbance capabilities, which needs to be further solved.
In the ship propulsion motor speed control system, the ADRC can be used to achieve direct speed control [16], and accurately estimate and compensate the system state variables in the presence of internal and external disturbances. It requires the observer to quickly and accurately deal with the lumped disturbances that deviate from the established standard integral-series controlled object [9], but in the actual system, the bandwidth of the controller is limited by the switching frequency of the inverter. The state tracking and speed estimation of ADRC based on integral observation structure are limited, and it is difficult to take the advantages of active disturbance rejection controllers. As a result, when the ADRC is used to control the speed in the permanent-magnet synchronous motor of the ship, it has a poor suppression effect on the speed fluctuation caused by the sudden change in the load torque, the first and second harmonics caused by the current measurement error and sixth harmonic caused by inverter dead-time effects and non-sinusoidal flux linkage.
In order to solve the above problems, this paper proposes a PMSM active disturbance rejection control method based on improved observer and current disturbance feedforward compensation. First, in order to improve the tracking performance and the disturbance suppression performance of ADRC, the state and disturbance observation laws of the extended state observer are improved to a proportional–integral control structure, which improves the dynamic response speed of the motor control system. Second, in order to further improve the system’s ability to suppress sudden changes in load torque and reduce speed fluctuations in the dynamic process, the measured current information is used to construct a disturbance feedforward compensation item to make the adjustment process faster and more stable. Third, in order to suppress the first, second, and sixth harmonics of the speed caused by the current measurement error and inverter nonlinearity, the continuous nonlinear function with bounded gain and the IIR low-pass filter are used to design dynamic filtering method, which can adjust the filtering weight according to the observation error, quickly track the current in the dynamic process and filter the current harmonics in the steady state, taking into account the dynamic response speed and steady-state control accuracy. Finally, this paper uses MATLAB/Simulink 2021a software to model the propeller load model of the ship’s electric propulsion system. On the permanent-magnet synchronous motor experimental platform, experiments such as motor load mutation, sudden current measurement error, ship staged start, and propeller torque mutation are carried out. Experimental results on motor speed and current verify the effectiveness and superiority of the proposed method. There are two main highlights in this paper:
(1)
In this study, the tracking performance and anti-disturbance performance of the controller are improved by reconstructing the structure of the extended state observer, which significantly enhances the dynamic response speed of the system.
(2)
On the basis of the original controller, a current disturbance compensation item with dynamic filtering is designed to reduce the speed fluctuation when the motor torque changes suddenly. The dynamic filtering link effectively filters the current harmonics, suppresses the influence of periodic disturbances such as current measurement errors of the motor system, and improves the current-control accuracy.
The remainder of the paper is organized as follows: Section 2 presents the mathematical model of the electrical propulsion system. Section 3 presents the Improved Active Disturbance Rejection Controller (IADRC) and the smooth switching strategy. Section 4 presents the stability and the performance analysis of the improved structure. Section 5 discusses the overall experimental results of the system. Finally, the conclusions are given in Section 6.

2. Ship Propulsion Motor and Its Propeller Load Model

2.1. Motor State Equation under the Framework of ADRC

Assuming that in the core saturation, eddy current and hysteresis losses are ignored, using the vector control strategy with id = 0, the surface-mounted permanent-magnet synchronous motor (SPMSM) motion equation, and voltage equation in the synchronous rotating coordinate system are expressed as follows:
T e T L = J ω ˙ m + B ω m T e = K t i q = 1.5 n p λ r i q u d = L s i ˙ d + R s i d L s n p ω m i q u q = L s i ˙ q + R s i q + n p ω m ( L s i d + λ r )
where Te and TL are electromagnetic torque and load torque, respectively; J is the moment of inertia; ωm is mechanical angular speed; B is friction coefficient; Kt is torque coefficient; id and iq are d-axis and q-axis current, respectively; ud and uq are the d-axis and q-axis voltages; Ls and Rs are the stator inductance and resistance, respectively; np is the number of pole pairs of the PMSM; and λr is the magnetic flux of the permanent magnet.
The second-order equation of the mechanical angular speed in PMSM can be obtained from Equation (1):
ω ¨ m = K t J L s u q T ˙ L J B ω ˙ m J n p K t λ r J L s ω m K t R s J L s i q
Let
b = K t J L s
f ω m , i q , T L = b b 0 u q T ˙ L J B ω ˙ m J n p K t λ r J L s ω m K t R s J L s i q
where b0 is the estimated value of b; then, Equation (2) can be rewritten as follows:
y ¨ = b 0 u + f ω m , i q , T L
where y and u represent ωm and uq, respectively; and f(ωm, iq, TL) represents the lumped disturbances of PMSM. By estimating and compensating f(ωm, iq, TL) in the control system, the influence of various disturbances and unmodeled dynamics on the system can be effectively suppressed, and the original complex permanent-magnet synchronous motor system control problem is transformed into the control of a unit gain double integrator. For this problem, the PD control method can achieve a good control performance and has strong robustness [9].

2.2. Mathematical Model of the Propeller Load

In the working condition of ship electric propulsion, the propulsion motor drives the propeller to overcome the seawater resistance moment to rotate, and at the same time, the seawater generates thrust on the propeller and transmits it to the hull to push the ship to sail [25,26,27]. This process requires the motor to provide electromagnetic torque for the propeller, and the motor works under complex conditions where the propeller interacts with the hull and seawater. The thrust P (unit is N) and torque Tp (unit is N·m) generated when the propeller rotates are
P = K p ρ D 2 v p 2 + n 2 D 2
T p = K T ρ D 3 v p 2 + n 2 D 2
where ρ is the density of seawater, generally 1025 kg/m3; D is the diameter of the propeller and the unit is m; and vp is propeller velocity relative to water and the unit is m/s. n is the rotation speed of the propeller and the unit is r/s; Kp is the dimensionless coefficient of propeller thrust; and KT is the dimensionless coefficient of drag torque, which are the functions of the advance ratio L′. According to the derivation process in Appendix A, when n > 0, Kp and KT can be expressed approximately as an 8-order Chebyshev polynomial [28]:
K p = 0.348 0.051 L 1.224 L 2 + 0.007 L 3 + 3.618 L 4 1.273 L 5 4.670 L 6 + 1.155 L 7 + 1.944 L 8
K T = 0.047 0.012 L 0.139 L 2 + 0.068 L 3 + 0.407 L 4 0.306 L 5 0.571 L 6 + 0.228 L 7 + 0.261 L 8
When the propeller rotates at speed n to propel the ship forward, and the advance speed of the propeller is vp, then the distance traveled by the propeller per revolution is
h p = v p n
The ratio of the hp to the propeller diameter D is called the advance ratio L of the propeller, also known as the advance speed coefficient. To avoid the situation that the values of Kp, KT, and L are infinite when the propeller speed is close to 0 [27], a bounded form of advance speed ratio is adopted in this paper. The expression of the advance ratio L′ in bounded form is [28]
L = v p v p 2 + n 2 D 2
When the propeller rotates, it can cause additional resistance to the ship, and the impact of this resistance increase can be expressed by the thrust deduction coefficient t; so, the effective thrust Pe (unit is N) generated by the propeller is
P e = 1 t P = 1 t K p ρ D 2 v p 2 + n 2 D 2
As the ship sails through the water, the water around the ship is affected by the ship, which creates currents that flow with the ship. This phenomenon is called a wake. The wake will cause a difference between propeller forward speed vp and ship speed vs. The relationship between the two is
v p = 1 w v s
where w is the wake coefficient; and the total resistance suffered by the ship is
R s = c v s 2
where c is the resistance coefficient; and the motion equation of the ship–machine–propeller system can be expressed as
m + Δ m d v s d t = P e P s
where m is the mass of the hull and Δm is the mass of the attached water, which is generally taken as 5–15% of the mass of the hull. The experimental part of this paper adopts the ship parameters in Table 1.
According to the data in Table 1 and Equations (6)–(15), the propeller load model is built as shown in Figure 1:
Therefore, according to the motor speed and the propeller load model, the load torque generated by the propeller to the propulsion motor can be obtained.

3. Improved Active Disturbance Rejection Controller (IADRC)

In the ADRC method, the tracking differentiator (TD) generates tracking signals and differential tracking signals based on the input signals. The extended state observer (ESO) is used to observe the system state and compensate for lumped disturbances. The control law combines the deviation between the observation value and the given value and the disturbance compensation value, and the three parts work together to control the plant.

3.1. Tracking Differentiator (TD)

The mathematical model of the tracking differentiator is expressed as
ω ˙ 1 = ω 2 ω ˙ 2 = f h a n ( ω 1 ω ref , ω 2 , r , h )
where ωref is the reference input speed; ω1 is the tracking signal of the reference speed; ω2 is the differential of the tracking signal; r is the speed factor; h is the control step; and fhan(f1,f2,r,h) is the optimal tracking function. They are given by
d = r h 2 y = f 1 + h f 2 a 1 = d ( d + 8 | y | ) a 2 = h f 2 + sign ( y ) ( a 1 d ) / 2 a 3 = ( sign ( y + d ) sign ( y d ) ) / 2 a 4 = ( h f 2 + y a 2 ) a 3 + a 2 a 5 = ( sign ( a 4 + d ) sign ( a 4 d ) ) / 2 f h a n ( f 1 , f 2 , r , h ) = r ( a 4 d sign ( a 4 ) ) a 5 r sign ( a 4 )
For the fhan function, the selection of r and h affects the performance of the tracking differentiator, and the differential signal helps improve the speed of motor speed regulation. Taking r = 50,000 and h = 0.0001, the TD tracking characteristics are shown in Figure 2:
Figure 2 shows that the traditional fhan function has different tracking performance for different rotational speed tracking values, but after a simple transformation,
ω ˙ 1 = ω 2 , ω r 1 = ω 1 / λ , ω r 2 = ω 2 / λ ω ˙ 2 = f h a n ( ω 1 λ ω ref , ω 2 , r , h )
For different given signals, the same tracking performance can be obtained, and the transformed function tracking characteristics are shown in Figure 3:
The tracking speed can be adjusted by adjusting λ to obtain the desired tracking effect.

3.2. Extended State Observer (ESO)

According to the motor motion equation, a third-order extended state observer (as shown in Figure 4) is designed, and then an ADRC controller is designed. According to the controller design method [16], the Linear ADRC controller directly outputs the reference voltage to adjust the motor speed, and the q-axis current feedback link is canceled, reducing the dependence of the motor control system on the current sensor.
According to Equation (5), the state equation of PMSM can be expressed as follows:
y = x 1 = ω m x ˙ 1 = x 2 = ω ˙ m x ˙ 2 = x 3 + b 0 u x ˙ 3 = f ˙ ω m , i q , T L
According to the state Equation (19) of the motor, the third-order extended state observer is designed as follows:
e = ω m z 1 = y y ^ z ˙ 1 = z 2 + β 1 e z ˙ 2 = z 3 + β 2 e + b 0 u z ˙ 3 = β 3 e
where z1, z2 and z3 are the observed values of x1, x2 and x3, respectively. And the control law u can be designed as
u 0 = k 1 ω r 1 z 1 + k 2 ω r 2 z 2
u = u 0 z 3 b 0
When the tracking transition process ends, ωr1 becomes a given signal value ωref, ωr2 is 0, and Equation (22) becomes:
u = k 1 ω r 1 z 1 k 2 z 2 z 3 b 0
Extended State Observer (ESO) is the core of the Active Disturbance Rejection Controller (ADRC), which is used to estimate the state variables and lumped disturbance f(ωm, iq, TL) of the system. Therefore, the operating characteristics of the ESO play a key role in the control performance of the ADRC.

3.3. Improved Extended State Observer (IESO)

Refs. [23,24,25,26,27,28,29] show that the use of the proportional–integral disturbance update law can improve the anti-disturbance performance of the second-order extended state observer, which is a part of the first-order ADRC controller. In this paper, to improve the state tracking and estimation capabilities of the third-order observer, this section reconstructs the structure of the third-order extended state observer. The state equation of the improved extended state observer is as follows:
e = ω m z 1 = y y ^ z ˙ 1 = z 2 z 2 = z 3 + b 0 u s + β 1 e z 3 = β 2 e + β 3 e s
where s is the Laplace operator. Compared with the structure in Figure 4, the z2 and z3 observe update laws of the IESO are designed as a proportional–integral control structure. The structure of the improved ADRC controller is shown in Figure 5:
The control law u can be designed as
u 0 = k 1 ω r 1 z 1 + k 2 ω r 2 z 2
u = u 0 z 3 b 0
When the tracking transition process ends, Equations (25) and (26) become
u 0 = k 1 ω r 1 z 1 k 2 z 2 u = k 1 ω r 1 z 1 k 2 z 2 z 3 b 0

3.4. q-Axis Current Disturbance Compensation Item

It can be seen from Figure 4 and Figure 5 that the original active disturbance rejection controller regards the motor current disturbance item as a part of the lumped disturbance, only estimates compensation through the disturbance observation link, and cancels the q-axis current feedback, which can avoid current measurement errors for the impact. However, when the load torque changes greatly, the proportion of this item in the aggregate disturbance is larger. At this time, it is difficult to obtain excellent anti-disturbance performance by only relying on the observation link to deal with these disturbances. In addition, according to Equations (3)–(5), the coefficient of the current disturbance term and the ADRC gain coefficient b only differ by the resistance value RS. This is a parameter value that is convenient to obtain; so, it is necessary to introduce a current disturbance compensation term in the system.
Rewrite Equation (5) as
ω ¨ m = b 0 u q + f ω m , i q , T L b 0 R s i q
Then, the lumped disturbances become
f ω m , i q , T L = b b 0 u q T ˙ L J B ω ˙ m J n p K t λ r J L s ω m b b 0 R s i q
In Equations (28) and (29), the q-axis current disturbance item in the total disturbance is separated and compensated into the disturbance observation item z3. Through feedforward compensation, this item can quickly compensate the disturbance caused by the change in the motor load torque, and further improve the system tracking and anti-disturbance performance. The designed compensation structure (IADRC + iq, “iq” represents the q-axis current disturbance compensation term, IADRC with the q-axis current disturbance compensation item) is shown in Section 3.5 (except for the smooth switching strategy).

3.5. Design of Smooth Switching Filtering Strategy

After introducing the q-axis current disturbance item, although the system anti-disturbance performance is further improved, the harmonics in the current also directly participate in the system control, resulting in the increase in the first, second, and sixth harmonic content of the motor current. Since the low-pass filtering effect of the system’s rotational inertia, the sixth harmonic has little effect on the motor speed.
At present, the notch filter is mainly used to suppress the current harmonics, but the increased resonance term makes it easy to reduce the system phase margin, resulting in system instability. Since these harmonics are introduced by compensation, the first, second, and sixth harmonics existing in the q-axis current can be filtered out by using an IIR low-pass filter, which fundamentally solves this problem. However, when the filter cutoff frequency is low, it will affect the effect of the current compensation item in the dynamic process (filter delay), and then affect the anti-disturbance performance of the motor. Therefore, it is necessary to adjust the filtering effect according to the working condition of the motor.
To solve the above problems, this paper obtains a nonlinear switching function that is convenient for parameter adjustment by modifying the fal function. Equation (30) of the fal function is as follows:
f a l e , α , δ = e δ 1 α | e | δ | e | α s i g n e | e | > δ
Figure 6 is the fal function characteristic curve when α = 0.01; δ takes different values and δ = 0.01, and α takes different values.
Because the fal function gain item is a continuous nonlinear function with bounded gain, and by adjusting the parameter values α, δ can easily obtain different nonlinear characteristics; so, it is suitable to be transformed into a smooth transition function for the switching process. Equation (31) of the designed transition function f is as follows:
f e , α , δ = 1 | e | δ | e | α s i g n e δ 1 α e | e | > δ
Figure 7 is the f function characteristic curve when α = 0.01, δ takes different values and δ = 0.2, while α takes different values.
The output of the f function is the weight of the motor current item after filtering, and 1 − f is the weight of the original current item. Scaling the observed deviation e in Equation (24) can match the gain range of the switching function. Through experiments, the fluctuation range of the error e during the steady-state operation of the motor can be easily obtained. At this time, the output of the f function should be 1. During the operation of the motor, the weight is dynamically adjusted according to the working conditions. Since the f function is continuous, there will be no sudden change in the current during the adjustment process, and the dynamic and steady-state performance of the system will be taken into account. The effectiveness of the method is verified by experiments.
The improved controller structure (IFADRC + iq, IADRC+iq controller with dynamic filtering method) is shown in Figure 8:
Figure 8 shows that as the observe error e changes, the current filter weights change accordingly. In the dynamic adjustment process, the current in the compensation item is closer to the real current to ensure the response speed; in the steady state process, all the compensation items are the filtered current, and the current harmonic content is very small when the cut-off frequency is very low, which avoids the current and speed ripples in the steady state.

4. IADRC Stability and Frequency Domain Performance Analysis

4.1. IADRC Stability Analysis

The state equation of the motor in Equation (19) can be written as
x ˙ = A x + B u + H y = C x
where
A = 0 1 0 0 β 1 1 0 β 2 0 , B = 0 b 0 0 , C = [ 1 0 0 ] , H = 0 0 1 f ˙ ω m , i q , T L + 0 0 0 0 β 1 0 0 β 2 0 x
The improved observer state equation in Equation (24) can be written as
z ˙ = A z + B u + L y y ^ + W y ^ = C z
where
L = 0 0 β 3 , W = 0 β 1 β 2 ω ˙ m
Let ei = xizi, I = 1, 2, 3. and subtract Equation (33) from Equation (32), the error equation can be written as
e ˙ = A e e + H W
where
A e = A L C = 0 1 0 0 β 1 1 β 3 β 2 0
The characteristic equation of the matrix Ae is
s I A e = s 1 0 0 s + β 1 1 β 3 β 2 s = s 3 + β 1 s 2 + β 2 s + β 3
when H-W is bounded and the root of the characteristic equation is on the left side of the complex plane, the IESO is bounded-input bounded-output (BIBO) stable [9]. Let s3 + β1s2 + β2s + β3 = (s + ωo)3. When β1 = 3ωo, β2 = 3ωo2, β3 = ωo3, and ωo > 0, the observer is stable and has a faster convergence speed.
Theorem 1. 
The IADRC in Figure 5 yields a BIBO stable closed-loop system if the observer in Equation (33) and the control law Equation (27) for the plant are stable.
Proof of Theorem 1. 
Equation (27) can be written in a state feedback form of u = (1/b0)[− k 1  −  k 2 − 1]z = Fz, where F = (1/b0)[− k 1 k 2 − 1]. The closed-loop system is then represented by the state-space equation of
x ˙ z ˙ = A B F L C A L C + B F x z + G G ω r 1 + H W
where G = k 1 B/b0. Equation (36) is BIBO stable if its eigenvalues are in the left half plane. By applying elementary row and column operations, it is obvious that the closed-loop eigenvalues satisfy
e i g A B F L C A L C + B F = e i g A + B F B F 0 A L C = e i g A + B F e i g A L C = { root of s 2 + k 2 + β 1 s + k 2 } { root of s 3 + β 1 s 2 + β 2 s + β 3 } Q . E . D .
Since ωr1, as the reference signal, is always bounded, the only nontrivial condition on the plant is that H and W are bounded. In other words, the system state x1, x2, and disturbance f(ωm,iq,TL) must be differentiable, which is a reasonable assumption. □

4.2. IADRC Frequency Domain Performance Analysis

Figure 9 is the frequency domain model of permanent-magnet synchronous motor using IADRC controller (the IADRC structure is presented in Figure 5 and Figure 8). From Figure 9, it can be seen that there are various disturbances and unmodeled dynamics in the motor-speed control system, these factors will affect the motor control effect. For the controller used in this paper, unmodeled dynamics such as these external loads and current measurement errors, flux harmonics, etc., can be regarded as disturbance items in the feedforward channel of the control loop; so, it is very important to improve the anti-disturbance capabilities of the controller.
Next, the frequency domain characteristics of the control system are analyzed, and Equations (20) and (23) are regarded as a dual-input single-output system with two inputs y, ωr1 and output u, and the transfer function from y and ωr1 to u is derived as follows [30]:
u s ω m s = β 1 k 1 + β 2 k 2 + β 3 s 2 + β 2 k 1 + β 3 k 2 s + β 3 k 1 s 3 + ( β 1 + k 2 ) s 2 + β 1 k 2 + β 2 + k 1 s 1 b 0
  u s ω r 1 s = k 1 s 3 + k 1 β 1 s 2 + k 1 β 2 s + k 1 β 3 s 3 + ( β 1 + k 2 ) s 2 + β 1 k 2 + β 2 + k 1 s 1 b 0
Neglecting the sampling delay, the frequency domain model of the permanent-magnet synchronous motor speed control system under ADRC is shown in Figure 10:
where
  G f s = ω m s u s u s ω r 1 s = k 1 s 3 + k 1 β 1 s 2 + k 1 β 2 s + k 1 β 3 β 1 k 1 + β 2 k 2 + β 3 s 2 + β 2 k 1 + β 3 k 2 s + β 3 k 1   G c s = b 0 u s ω m s = β 1 k 1 + β 2 k 2 + β 3 s 2 + β 2 k 1 + β 3 k 2 s + β 3 k 1 s 3 + ( β 1 + k 2 ) s 2 + β 1 k 2 + β 2 + k 1 s
when using IADRC,
G f I s = k 1 s 3 + k 1 β 1 s 2 + k 1 β 2 s + k 1 β 3 β 1 k 2 + β 2 s 3 + β 1 k 1 + β 2 k 2 + β 3 s 2 + β 2 k 1 + β 3 k 2 s + β 3 k 1   G c I s = β 1 k 2 + β 2 s 3 + β 1 k 1 + β 2 k 2 + β 3 s 2 + β 2 k 1 + β 3 k 2 s + β 3 k 1 s 3 + ( β 1 + k 2 ) s 2 + k 1 s
Figure 10 shows that the disturbances in the system are mainly concentrated in two places: current loop feedforward channel disturbances (Dead-Time Effects, etc.), and speed loop feedforward channel disturbances (load torque, flux harmonics, etc.). These disturbances include both DC disturbances and periodic disturbances; so, it is necessary to analyze the frequency domain anti-disturbance capabilities of the control system for different disturbance locations.
When using ADRC and IADRC, respectively, the anti-disturbance transfer function of the current loop feedforward channel is:
  G d i s = G i s G s s 1 + G c s G i s G s s / b 0 + G i s G s s n p λ r
    G d i I s = G i s G s s 1 + G c I s G i s G s s / b 0 + G i s G s s n p λ r
Therefore, the amplitude–frequency curve can be obtained as in Figure 11.
When using ADRC and IADRC, respectively, the anti-disturbance transfer function of the speed loop feedforward channel is
G d s s = G s s 1 + G c s G i s G s s / b 0 + G i s G s s n p λ r
G d s I s = G s s 1 + G c I s G i s G s s / b 0 + G i s G s s n p λ r
The amplitude–frequency curve can be obtained as in Figure 12.
The amplitude–frequency characteristic curve in Figure 11 and Figure 12 shows that the improved control structure can effectively improve the anti-disturbance capabilities of the motor in the low-frequency range.

5. Experimental and Simulation Analysis

To verify the effectiveness of the proposed method, this paper conducts experiments on the Teknic-2310P motor platform produced by the TEXAS INSTRUMENTS company and uses a TMS320F28379D dual-core microprocessor to execute the researched control algorithm. Another PMSM is mechanically coupled with the tested motor as a load. The experimental platform is shown in Figure 13:
In the experiment, under the FOC framework, the improved active disturbance rejection direct speed control method proposed in this paper is used to drive the power motor, and the torque output of the load motor is controlled according to the load model of the ship electric propulsion system. The overall structure of the ship’s propulsion motor vector control model is shown in Figure 14:
The propeller load parameters, motor parameters, and controller parameters used in the experiment are shown in Table 1, Table 2 and Table 3, respectively. Among them, the PI controller adopts the parameters set by Matlab example for the platform, and the parameters of the extended state observer adopt the bandwidth method (in Section 4) configuration.

5.1. IADRC Disturbance Suppression Experiment

To verify the superiority of the improved method of ADRC proposed, firstly, the three control methods of PI, ADRC, and IADRC are used to test the anti-disturbance capabilities of the speed control system when the external load is suddenly added. The experimental results are shown in Figure 15:
In Figure 15a, when PI is used, the speed drops by 21.2 rpm, and the adjustment time is 0.14 s; when ADRC is used, the speed drops by 20.2 rpm, and the adjustment time is 0.15 s; and when IADRC is used, the speed drops by 10.4 rpm, and the adjustment time is 0.04 s. In Figure 15b, when PI is used, the speed drops by 49.6 rpm, and the adjustment time is 0.2 s; when ADRC is used, the speed drops by 48 rpm, and the adjustment time is 0.22 s; and when IADRC is used, the speed drops by 24.2 rpm, and the adjustment time is 0.05 s. In Figure 15c, when PI is used, the speed drops by 93.5 rpm, and the adjustment time is 0.23 s; when ADRC is used, the speed drops by 92 rpm, and the adjustment time is 0.26 s; and when IADRC is used, the speed drops by 48 rpm, and the adjustment time is 0.07 s. In Figure 15d, under the same load of 1A, IADRC has a faster start-up speed than ADRC. To show the experimental results more clearly, Table 4 summarizes these data.
The results in Table 4 indicate that when IADRC is used, compared with PI and ADRC, the speed fluctuation of the motor is reduced by about 50%, the adjustment time is shortened by more than 60%, and the start-up time is shortened by 64%. These results show that the tracking performance is significantly improved after the observer structure is improved, and the disturbance suppression ability is enhanced.

5.2. Experimental Analysis after Introducing q-Axis Current Disturbance Compensation Item

From the analysis in Section 3.4, it can be known that the current feedback value contains the information of the q-axis current disturbance item, which is missing in the original structure, which greatly limits the anti-disturbance capabilities of the system, and this problem is more serious when the load changes greatly. To solve this problem, this paper introduces the current disturbance term in the observation value z3. To verify the effectiveness of this method, the anti-disturbance experiment of the traditional ADRC and IADRC before and after the introduction of the current disturbance item, when the system suddenly adds an external load; the dynamic response performance of the four methods is shown in Figure 16:
In Figure 16a, when using ADRC, the speed drops by 47 rpm, and the adjustment time is 0.21 s; and when using traditional ADRC + iq, the speed drops by 24 rpm, and the adjustment time is 0.09 s. When using IADRC, the speed drops by 20 rpm, and the adjustment time is 0.04 s; and when using IADRC + iq, the speed drops by 7 rpm, and the adjustment time is 0.02 s. In Figure 16b, the speed drops by 89 rpm when ADRC is used, and the adjustment time is 0.25 s; and when ADRC + iq is used, the speed drops by 43 rpm, and the adjustment time is 0.12 s. When using IADRC, the speed drops by 41 rpm, and the adjustment time is 0.06 s; and when using IADRC + iq, the speed drops by 15 rpm, and the adjustment time is 0.03 s. To show the experimental results more clearly, Table 5 summarizes these data.
The results in Table 5 show that introducing the q-axis current disturbance item in the ADRC and IADRC controller can reduce the speed fluctuation caused by the external load change, which can further improve the load-disturbance-suppression ability of the motor.
However, this is only an improvement on the disturbance suppression performance, and does not indicate the impact on the steady-state performance of the system. To further analyze the influence of the introduction of the q-axis current disturbance compensation item on the system, the following experiment is carried out. After the motor is started, the phase current offset error (iA_offset = 0.1A, iB_offset 0.04A) is added. The steady-state phase current waveform is shown in Figure 17a, and the response waveforms of IADRC and IADRC + iq (speed and phase current) are shown in Figure 17b.
In Figure 17a, when using PI, the phase current measurement errors will cause significant offset in motor phase current. However, phase current measurement error has little effect on ADRC and IADRC, because these two control structures do not introduce q-axis current, thus avoiding this problem.
Then, Figure 17b shows that when IADRC+iq is adopted, the phase currents are distorted and the motor speed ripples increase because the current harmonics will be introduced into the control system along with the current. Next, FFT analysis is carried out on the motor speed and phase current. The results show that the 1st, 2nd, 6th, and 12th harmonics in the speed and the 2nd, 3rd, 5th, and 7th harmonics in the phase current have content increase; the harmonic content values are shown in Table 6.
The FFT (Fast Fourier Transform Algorithm) analysis results in Table 6 show that when IADRC + iq is used (compared with IADRC), the motor speed and current harmonic total distortion increase by 1.03% and 5.33%, respectively. Which indicate that the q-axis current disturbance compensation item will have an adverse effect on the steady-state accuracy of the system. To overcome this problem, this paper adopts the method of dynamic filtering to deal with the q-axis current harmonics, and it will be verified experimentally in the next section.

5.3. Experimental Analysis of Dynamic Filtering Method

From the comprehensive Figure 8, Figure 16 and Figure 17, it can be seen that after the introduction of the q-axis current disturbance term, although the system’s anti-disturbance capabilities are further improved, the harmonics in the current are also feedforward-compensated to the control law, directly participating in the calculation of the voltage vector, and will lead to an increase in the harmonic content in the current. According to the analysis in Section 3.5, using a low-pass filter and designing a smooth switching strategy can solve this problem better; the weights are dynamically adjusted according to the working conditions when the motor is running, taking into account both dynamic and steady-state performance. After the current measurement error is introduced, the disturbance-suppression experiment results and phase current waveforms under the IADRC and IFADRC + iq methods are shown in Figure 18:
Figure 18a shows that after adopting the smooth switching filter strategy, the anti-disturbance capabilities of the motor are still significantly improved compared with IADRC, and the regulation process is smoother. The speed drops of the IADRC and IFADRC (IFADRC + iq) are 41 rpm and 16 rpm, respectively. The FFT analysis results of the phase current (Figure 18b) is shown in Table 7.
The FFT results in Table 7 show that when IFADRC + iq is used, in other words, the dynamic filtering method, the third, fifth, seventh and eleventh harmonics are reduced by 0.16%, 3.62%, 3.48%, and 0.13%, respectively. In addition, the phase current waveforms (Figure 18b) show a significant improvement in current distortion.
This section completes the design of the controller proposed in this paper, which can be divided into three parts:
(1)
IADRC: Active Disturbance Rejection Controller with Improved Observation Structure;
(2)
The q-axis current disturbance compensation item;
(3)
Dynamic filtering strategy.
This method has the advantages of fast response, strong anti-disturbance capabilities, and low steady-state ripples.

5.4. Ship–Machine–Propeller Speed Control System Experiment

To further verify the starting performance of the proposed method under the condition of ship electric propulsion, the motor starting experiment under the propeller load was carried out (multiply the calculated propeller torque by 1/8000 to match the experimental motor load range); when the actual ship sailed, to avoid excessive torque, the method of staged starting is usually adopted. The experimental results of the staged starting of the propulsion motor under the three control methods are shown in Figure 19.
To make the comparison effect more obvious, this paper also conducted a set of direct starting experiments, and the results are shown in Figure 20:
The experimental results in Figure 20 show that the propulsion motor has a faster starting speed when using IFADRC (IFADRC + iq), which shows that the method has a stronger anti-disturbance ability when the load torque is constantly changing, and the motor starting performance and ship maneuverability are better. When PI, ADRC and IFADRC are used, the motor starting times are 0.19 s, 0.39 s and 0.39 s, respectively.
To further verify the control performance of the proposed method for the step load disturbance under the ship’s electric propulsion condition, the motor-speed control experiment under the propeller load was carried out. Ref. [31] indicated that the propeller may be exposed to the water surface in wind and waves, and the propeller torque of the ship will change significantly. Therefore, by changing 75% of the propeller torque to verify the anti-disturbance performance of the control system when the propeller is exposed to the sea due to severe sea conditions, the experimental results are shown in Figure 21:
The experimental results in Table 8 show that when PI control is used, the speed fluctuations are 278 rpm and −186 rpm, respectively; and when ADRC is used, the speed fluctuations are 279 rpm and −158 rpm, respectively. It is difficult to quickly suppress the disturbance effect when the propeller torque changes greatly, and the speed fluctuation is large. It is not conducive to ship maneuvering and stable navigation. When using IFADRC, the speed fluctuations are 64 rpm and −58 rpm, respectively, and the adjustment time is shorter, and the current response is fast without fluctuation. Compared with PI and ADRC, the total speed fluctuation of IFADRC is reduced by 73.7% and 72.1%, respectively. This method can still maintain good anti-disturbance capabilities under large disturbances.
To sum up, when the improved ADRC is used to propel the motor, the speed and torque can be adjusted quickly and smoothly. It can overcome various uncertainties and the influence of internal and external disturbances, and is suitable for ship electric propulsion conditions with high mobility requirements and time-varying disturbances.

5.5. Simulation Analysis of Another Ship Propulsion Motor

In order to better demonstrate the effectiveness of the method proposed in this paper, the motor parameters of another electric propulsion ship are used for simulation. The ship has a total length of 21.1 m, an empty ship weight of 29 t, and is equipped with an AC propulsion motor with a power level of 100 kw. The simulation was built in the MATLAB/simulink environment and carried out according to the parameter design method mentioned in the paper. Table 9 shows the motor parameters used in the simulation.
To show the simulation results more clearly, the results of Figure 22 are shown in Table 10.
The simulation results in Table 10 show that the method proposed in this paper is also suitable for high-power motor control. The control method proposed in this paper can still improve the tracking speed and disturbance suppression ability of the high-power propulsion motor system.

6. Conclusions and Outlook

This article proposes an improved ADRC method for ship-propulsion motors, which effectively improves the tracking performance, disturbance rejection performance, and steady-state control accuracy of the propulsion motor system.
(1) In order to improve the disturbance observation speed of the third-order ESO, the observer structure is reconstructed, and then the IADRC method is proposed, which has the advantages of fast-tracking speed and strong disturbance-suppression ability. In the 200 rpm load start-up and 3A anti-load disturbance experiments, the index of IADRC is about 0.14 s and 44 rpm smaller than ADRC, respectively.
(2) To utilize the measured current information, a current disturbance compensation item is built in the controller, which further improves the anti-disturbance performance of the system. The results of anti-disturbance experiments show that adding current compensation items to ADRC can reduce about 50% of the speed drop and 53% of the adjustment time; adding current compensation items to IADRC can reduce about 64% of the speed drop and 50% of the adjustment time. However, current experiments show that the current compensation term introduces current harmonics, and in IADRC, the fifth and seventh harmonics cause significant current distortion and degrade the current quality.
(3) To solve the above problems, this paper designs a dynamic filtering method, including IIR low-pass filter and nonlinear smooth switching function. While ensuring the dynamic performance, the harmonic content of the steady-state current is reduced. Current FFT analysis shows that the fifth and seventh harmonic content are reduced by about 70% after using the filtering method. The phase current waveforms show that the current distortion is significantly improved.
(4) Combined with the built propeller load, the effectiveness of the proposed method is verified. The results of the start-up experiment with propellers show that the start-up time of IFADRC is 0.2 s shorter than that of ADRC. In the experiment where the propeller torque changes by 75%, the speed fluctuation of IFADRC is more than 70% smaller than that of PI and ADRC, and the adjustment time is shortened by more than 50%.
Our future work will verify the effectiveness of the proposed algorithm in actual ships. This will significantly improve the anti-disturbance performance of propulsion motors, and enhance ship maneuverability and stability.

Author Contributions

Conceptualization, H.G., T.X. and Y.L.; methodology, H.G., T.X. and Y.L.; software, H.G. and T.X.; validation, H.G., T.X. and Q.Z.; formal analysis, H.G. and T.X.; investigation, H.G. and T.X.; resources, H.G.; data curation, H.G. and T.X.; writing—original draft preparation, H.G., T.X. and B.G.; writing—review and editing, H.G., T.X., S.L. and B.G.; visualization, H.G., T.X. and S.L.; supervision, H.G., T.X. and Q.Z.; project administration, H.G., T.X., Y.L. and S.L.; funding acquisition, H.G., Y.L. and Q.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the following: (a) National Natural Science Foundation of China, China (Grant No. 52201354, 51979021); and (b) Fundamental Research Funds for the Central Universities (Grant No. 3132023621).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

The Ship
PThrust of the propeller.
TpTorque of the propeller.
ρDensity of seawater.
DDiameter of the propeller.
vpPropeller velocity (relative to water).
nRotation speed of the propeller.
KpThe dimensionless coefficient of propeller thrust.
KTThe dimensionless coefficient of drag torque.
LAdvance ratio.
LAdvance ratio in bounded form.
hpDistance traveled by the propeller per revolution.
tThrust deduction coefficient.
PeEffective thrust.
vsShip speed.
wWake coefficient.
cResistance coefficient.
mMass of the hull.
ΔmMass of the attached water.
H/DPropeller screw pitch ratio.
A/AdPropeller disc square ratio.
The inverter and PMSM
SPMSMSurface-mounted permanent-magnet synchronous motor.
TeElectromagnetic torque.
TLLoad torque.
JMoment of inertia.
ωmMechanical angular speed.
BFriction coefficient.
KtTorque coefficient.
id, iqd-axis and q-axis current.
ud, uqd-axis and q-axis voltage.
Ls, RsStator inductance and resistance.
npPole pairs of the PMSM.
λrMagnetic flux of the permanent magnet.
bControl gain of the PMSM.
b0The estimated value of b.
f(ωm, iq, TL), f’(ωm, iq, TL)The lumped disturbances of the PMSM.
Proposed control structure
x1,x2,x3System state.
z1,z2,z3Estimation of system state.
ADRCActive disturbance rejection control
ADRC + iqADRC with the q-axis current disturbance compensation item
TDTracking differentiator
fhan(f1,f2,r,h)Optimal tracking function of TD.
rSpeed factor of TD
hControl step of TD
λTracking characteristic adjustment coefficient
ωrefReference speed
ω1, ωr1Tracking value of the reference speed.
ω2, ωr2Differential value of the reference speed.
ESOExtended state observer
eMeasured speed x1 minus observed speed z1.
ωoBandwidth of the observer.
β1,β2,β3Gains of the observer.
k 1 ,   k 2 ;   k 1 ,   k 2 Gains of the control law.
IADRCImproved Active Disturbance Rejection Controller
IESOImproved Extended State Observer
IADRC + iqIADRC with the q-axis current disturbance compensation item
IFADRC/IFADRC + iqIADRC+iq controller with dynamic filtering method.
f functionNonlinear switching function
α, δParameters of the f function

Appendix A

The thrust P (unit is N)and torque Tp1 (unit is N·m) are generated by the propeller, and the thrust coefficient Kp1, torque coefficient KT1 and advance ratio L for n 0 are defined as
K p 1 = P / ρ D 4 n 2 K T 1 = T p / ρ D 5 n 2 L = v p / ( D n )
where ρ is the density of seawater, generally 1025 kg/m3; D is the diameter of the propeller, and the unit is m; vp is propeller velocity relative to water, and the unit is m/s; n is the rotation speed of the propeller, and the unit is r/s; Kp1 is the dimensionless coefficient of propeller thrust; and KT1 is the dimensionless coefficient of drag torque, which is the function of the advance ratio L′. For a given screw pitch ratio H/D, the relation between Kp1, KT1 and L is expressed as
K p 1 = K p 1 ( L ) , K T 1 = K T 1 ( L )
On the Kp1-L plane or KT1-L plane, Kp1 and KT1 may stride across four quadrants generally and, Kp1, KT1, and L may stretch to infinite. Thus, Equation (A2) is named as a boundless form of the propeller properties. But, it is more convenient to use following form in the propeller’s all-around dynamic analysis. When n and vp are not equal to zero at the same time, they are defined as
K P = P / ρ D 2 v p 2 + D 2 n 2 = K p 1 / ( 1 + L 2 ) K T = T P / ρ D 3 v p 2 + D 2 n 2 = K T 1 / ( 1 + L 2 ) L = v p / v p 2 + D 2 n 2
From Equation (A3), it is deduced that
L = v p D N v p / D N 2 + 1
Thus, the following transformation relation is obtained:
L = L / 1 + L 2 n > 0 L = L / 1 + L 2 n < 0
It can be seen that L ( ,   ) is mapped into L ( 1 ,   1 ).
Generally, a function which is continuous for x ( 1 ,   1 ) can be expressed approximately as an n-order Chebyshev polynomial, that is
f x 1 2 a 0 T 0 ( x ) + a 1 T 1 ( x ) + + a n T n ( x )
where T0(x) = 1, T1(x) = x, T2(x) = 2x2 1, T2(x) = 4x3 3x, ….
The general recursive formula for Tk(x) is Tk+l − 2xTk(x) + Tk−1 = 0, (k = 1, 2, …). The Chebyshev polynomial has the following features. First of all, Tj, and Tk (j, k = l, 2, … (j  k)) are orthogonal to each other. The polynomial coefficients do not change due to the number of terms selected. Secondly, the fitting result of its finite items is the best approximation in the sense of minimum mean square error. Finally, the expression from a given Chebyshev polynomial can be easily replaced by an ordinary polynomial. As long as the Chebyshev polynomial coefficients a0~an are substituted into Equation (A6), and T0~Tn are expanded according to the definition, the corresponding ordinary polynomial coefficients b0~bn can be obtained.
0 n b i x i b 0 + b 1 x + + b n x n
The main disadvantage of using ordinary polynomial fitting directly is the non-generalizability of the fitting results. When the selected order n is different, in order to ensure a high fitting accuracy, the entire set of coefficients must be changed, and the fitting must be performed again. If the coefficients remain the same, simply increasing or decreasing the number of terms will reduce the fitting accuracy. When fitting with ordinary polynomials, in order to meet different precision requirements, many sets of coefficients must be given, and there are many tables, which is cumbersome. There is no such problem when using Chebyshev polynomial fitting.
The Chebyshev polynomial fitting is performed on the characteristics of the propeller, and the result is as follows:
K P ( L ) = 1 2 a 0 P T 0 ( L ) + a 1 P T 1 ( L ) + + a n P T n ( L ) K T ( L ) = 1 2 a 0 T T 0 ( L ) + a 1 T T 1 ( L ) + + a n T T n ( L )
The coefficients are calculated according to the following equations:
a i P = S P a i P ( 0.965 ) = 0.88 ( 0.35 a i P ( 0.965 ) + 0.65 a i P ( 1 ) ) i = 0 , 1 , , 8 a i T = S T a i T ( 0.965 ) = 0.918 ( 0.35 a i T ( 0.965 ) + 0.65 a i T ( 1 ) ) i = 0 , 1 , , 8
where aip, aiT is obtained from the attached Table 1 and attached Table 2 in Ref. [28] according to H/D = 0.9 and H/D = 1. The method of combining interpolation and proportional correction is used to determine aip and aiT, where 0.35 and 0.65 are interpolation coefficients, and SP and ST are proportional coefficients. Table A1 shows the final calculation results of the coefficients of Equation (A8).
Table A1. Chebyshev polynomial coefficients of propeller thrust and torque characteristics (n > 0).
Table A1. Chebyshev polynomial coefficients of propeller thrust and torque characteristics (n > 0).
a0a1a2a3a4a5a6a7
Kp0.3753−0.2378−0.1606−0.01930.00210.0532−0.02780.0205
KT0.0494−0.0302−0.0214−0.00420.00100.0064−0.00170.0039
a8a9a10a11a12a13a14a15
Kp0.01730.00190.00000.0141−0.0021−0.00190.0048−0.0036
KT0.00220.00120.00020.0020−0.0003−0.00030.0007−0.0004
Therefore, the 8th-order Chebyshev polynomial fitting polynomial with a better accuracy and time efficiency can be obtained. In addition, in order to facilitate the use in the propeller simulation, Chebyshev polynomials are often converted into ordinary polynomials, the conversion method is mentioned in Equations (A6) and (A7):
K P = 0.348 0.051 L 1.224 L 2 + 0.007 L 3 + 3.618 L 4 1.273 L 5 4.670 L 6 + 1.155 L 7 + 1.944 L 8
K T = 0.047 0.012 L 0.139 L 2 + 0.068 L 3 + 0.407 L 4 0.306 L 5 0.571 L 6 + 0.228 L 7 + 0.261 L 8

References

  1. Bai, H.; Yu, B.; Gu, W. Research on Position Sensorless Control of RDT Motor Based on Improved SMO with Continuous Hyperbolic Tangent Function and Improved Feedforward PLL. J. Mar. Sci. Eng. 2023, 11, 642. [Google Scholar] [CrossRef]
  2. Haxhiu, A.; Abdelhakim, A.; Kanerva, S.; Bogen, J. Electric Power Integration Schemes of the Hybrid Fuel Cells and Batteries-Fed Marine Vessels-An Overview. IEEE Trans. Transp. Electrif. 2022, 8, 1885–1905. [Google Scholar] [CrossRef]
  3. Reusser, C.A.; Young, H.A.; Osses, J.R.P.; Perez, M.A.; Simmonds, O.J. Power Electronics and Drives Applications to Modern Ship Propulsion Systems. IEEE Ind. Electron. Mag. 2020, 14, 106–122. [Google Scholar] [CrossRef]
  4. Yang, J.; Chen, W.H.; Li, S.; Guo, L.; Yan, Y. Disturbance/Uncertainty Estimation and Attenuation Techniques in PMSM Drives—A Survey. IEEE Trans. Ind. Electron. 2017, 64, 3273–3285. [Google Scholar] [CrossRef]
  5. Bu, F.; Yang, Z.; Gao, Y.; Pan, Z.; Pu, T.; Degano, M.; Gerada, C. Speed Ripple Reduction of Direct-Drive PMSM Servo System at Low-Speed Operation Using Virtual Cogging Torque Control Method. IEEE Trans. Ind. Electron. 2021, 68, 160–174. [Google Scholar] [CrossRef]
  6. Chen, H.; Wang, X.; Benbouzid, M.; Charpentier, J.-F.; Aϊt-Ahmed, N.; Han, J. Improved Fractional-Order PID Controller of a PMSM-Based Wave Compensation System for Offshore Ship Cranes. J. Mar. Sci. Eng. 2022, 10, 1238. [Google Scholar] [CrossRef]
  7. Jia, D.; Wang, T.; Amirat, Y.; Tang, Y. An Inequality Indicator for High-Resistance Connection Fault Diagnosis in Marine Current Turbine. J. Mar. Sci. Eng. 2023, 11, 97. [Google Scholar] [CrossRef]
  8. Han, J. Extended state observer for a class of uncertain plants. J. Control. Decis. 1995, 10, 85–88. (In Chinese) [Google Scholar] [CrossRef]
  9. Gao, Z. Scaling and bandwidth parameterization based controller tuning. In Proceedings of the 2003 American Control Conference, Denver, CO, USA, 4–6 June 2003; pp. 4989–4996. [Google Scholar]
  10. Li, S.; Liu, Z. Adaptive Speed Control for Permanent-Magnet Synchronous Motor System With Variations of Load Inertia. IEEE Trans. Ind. Electron. 2009, 56, 3050–3059. [Google Scholar] [CrossRef]
  11. Liu, M.; Chan, K.W.; Hu, J.; Xu, W.; Rodriguez, J. Model Predictive Direct Speed Control With Torque Oscillation Reduction for PMSM Drives. IEEE Trans. Ind. Inform. 2019, 15, 4944–4956. [Google Scholar] [CrossRef]
  12. Liu, Y.; Zhang, Z.; Guo, H.; Zhang, Q.; Yang, P.; Zhu, P.; Zhuang, X. An Active DC-Link Stabilization Strategy for Permanent Magnet Synchronous Motor Drive System Based on Inverter Voltage Compensation. IEEE-ASME Trans. Mechatron. 2022, 27, 2786–2795. [Google Scholar] [CrossRef]
  13. Zhang, Z.; Liu, Y.; Liang, X.; Guo, H.; Zhuang, X. Robust Model Predictive Current Control of PMSM Based on Nonlinear Extended State Observer. IEEE J. Emerg. Sel. Top. Power Electron. 2023, 11, 862–873. [Google Scholar] [CrossRef]
  14. Qu, L.; Qiao, W.; Qu, L. Active-Disturbance-Rejection-Based Sliding-Mode Current Control for Permanent-Magnet Synchronous Motors. IEEE Trans. Power Electron. 2021, 36, 751–760. [Google Scholar] [CrossRef]
  15. Li, J.; Xia, Y.; Qi, X.; Gao, Z. On the Necessity, Scheme, and Basis of the Linear–Nonlinear Switching in Active Disturbance Rejection Control. IEEE Trans. Ind. Electron. 2017, 64, 1425–1435. [Google Scholar] [CrossRef]
  16. Hao, Z.; Yang, Y.; Gong, Y.; Hao, Z.; Zhang, C.; Song, H.; Zhang, J. Linear/Nonlinear Active Disturbance Rejection Switching Control for Permanent Magnet Synchronous Motors. IEEE Trans. Power Electron. 2021, 36, 9334–9347. [Google Scholar] [CrossRef]
  17. Lin, P.; Wu, Z.; Liu, K.Z.; Sun, X.M. A Class of Linear–Nonlinear Switching Active Disturbance Rejection Speed and Current Controllers for PMSM. IEEE Trans. Power Electron. 2021, 36, 14366–14382. [Google Scholar] [CrossRef]
  18. Zhu, L.; Zhang, G.; Jing, R.; Bi, G.; Xiang, R.; Wang, G.; Xu, D. Nonlinear Active Disturbance Rejection Control Strategy for Permanent Magnet Synchronous Motor Drives. IEEE Trans. Energy Convers. 2022, 37, 2119–2129. [Google Scholar] [CrossRef]
  19. Zhang, C.; Zhang, C.; Li, L.; Liu, H. An Enhanced Nonlinear ADRC Speed Control Method for Electric Propulsion System: Modeling, Analysis, and Validation. IEEE Trans. Power Electron. 2023, 38, 4520–4528. [Google Scholar] [CrossRef]
  20. Wang, Z.; Zhao, J.; Wang, L.; Li, M.; Hu, Y. Combined Vector Resonant and Active Disturbance Rejection Control for PMSLM Current Harmonic Suppression. IEEE Trans. Ind. Inform. 2020, 16, 5691–5702. [Google Scholar] [CrossRef]
  21. Wang, B.; Tian, M.; Yu, Y.; Dong, Q.; Xu, D. Enhanced ADRC with Quasi-Resonant Control for PMSM Speed Regulation Considering Aperiodic and Periodic Disturbances. IEEE Trans. Transp. Electrif. 2022, 8, 3568–3577. [Google Scholar] [CrossRef]
  22. Hu, M.; Hua, W.; Wang, Z.; Li, S.; Wang, P.; Wang, Y. Selective Periodic Disturbance Elimination Using Extended Harmonic State Observer for Smooth Speed Control in PMSM Drives. IEEE Trans. Power Electron. 2022, 37, 13288–13298. [Google Scholar] [CrossRef]
  23. Lin, S.; Cao, Y.; Li, C.; Wang, Z.; Shi, T.; Xia, C. Two-Degree-of-Freedom Active Disturbance Rejection Current Control for Permanent Magnet Synchronous Motors. IEEE Trans. Power Electron. 2023, 38, 3640–3652. [Google Scholar] [CrossRef]
  24. Tian, M.; Wang, B.; Yu, Y.; Dong, Q.; Xu, D. Adaptive Active Disturbance Rejection Control for Uncertain Current Ripples Suppression of PMSM Drives. IEEE Trans. Ind. Electron. 2023, 1–11. [Google Scholar] [CrossRef]
  25. Zhao, K.Q.; Yang, L.; Hu, H.X.; Shuang, Z. Ship PMSM Nonlinear ADRC Parameter Self-tuning Based on Neural Network. In Proceedings of the 25th International Conference on Electrical Machines and Systems, Chiang Mai, Thailand, 29 November–2 December 2022; pp. 1–5. [Google Scholar]
  26. Bai, H. Position Estimation of a PMSM in an Electric Propulsion Ship System Based on High-Frequency Injection. Electronics 2020, 9, 276. [Google Scholar] [CrossRef]
  27. Sunarsih; Izzuddin, N.; Priyanto, A. Thrust and Torque Estimation Scheme Based on Chebyshev Polynomial for Ship Manoeuvring Simulation. In Proceedings of the 2nd International Materials, Industrial, and Manufacturing Engineering Conference, Bali, Indonesia, 4–6 February 2015; pp. 359–363. [Google Scholar] [CrossRef]
  28. Li, D.; Wang, Z.; Chi, H. Chebyshev fitting way and error analysis for propeller atlas across four quadrants. J. Mar. Sci. Appl. 2002, 1, 52–59. [Google Scholar] [CrossRef]
  29. Zuo, Y.; Zhu, X.; Quan, L.; Zhan, C.; Du, Y.; Xiang, Z. Active Disturbance Rejection Controller for Speed Control of Electrical Drives Using Phase-Locking Loop Observer. IEEE Trans. Ind. Electron. 2019, 66, 1748–1759. [Google Scholar] [CrossRef]
  30. Jin, H.; Song, J.; Lan, W.; Gao, Z. On the characteristics of ADRC: A PID interpretation. Sci. China Inf. Sci. 2020, 63, 258–260. [Google Scholar] [CrossRef]
  31. Taskar, B.; Yum, K.K.; Steen, S.; Pedersen, E. The effect of waves on engine-propeller dynamics and propulsion performance of ships. Ocean Eng. 2016, 122, 262–277. [Google Scholar] [CrossRef]
Figure 1. Ship propeller load model structure.
Figure 1. Ship propeller load model structure.
Jmse 11 01751 g001
Figure 2. Transition process tracking curves with different tracking values.
Figure 2. Transition process tracking curves with different tracking values.
Jmse 11 01751 g002
Figure 3. Transition process tracking curves with different tracking values (after transformation).
Figure 3. Transition process tracking curves with different tracking values (after transformation).
Jmse 11 01751 g003
Figure 4. Active disturbance rejection controller structure (ADRC).
Figure 4. Active disturbance rejection controller structure (ADRC).
Jmse 11 01751 g004
Figure 5. Improved active disturbance rejection controller (IADRC) structure.
Figure 5. Improved active disturbance rejection controller (IADRC) structure.
Jmse 11 01751 g005
Figure 6. Outputs of the function fal with different α or δ.
Figure 6. Outputs of the function fal with different α or δ.
Jmse 11 01751 g006
Figure 7. Outputs of the function f with different α or δ.
Figure 7. Outputs of the function f with different α or δ.
Jmse 11 01751 g007
Figure 8. Improved active disturbance rejection controller with current disturbance compensation.
Figure 8. Improved active disturbance rejection controller with current disturbance compensation.
Jmse 11 01751 g008
Figure 9. Distribution of disturbance/unmodeled dynamics in PMSM systems.
Figure 9. Distribution of disturbance/unmodeled dynamics in PMSM systems.
Jmse 11 01751 g009
Figure 10. Equivalent feedback system of ADRC.
Figure 10. Equivalent feedback system of ADRC.
Jmse 11 01751 g010
Figure 11. Disturbance transfer function amplitude–frequency curve of q-axis current closed-loop system with ADRC and IADRC, respectively.
Figure 11. Disturbance transfer function amplitude–frequency curve of q-axis current closed-loop system with ADRC and IADRC, respectively.
Jmse 11 01751 g011
Figure 12. Disturbance transfer function amplitude–frequency curve of the speed closed-loop system with ADRC and IADRC, respectively.
Figure 12. Disturbance transfer function amplitude–frequency curve of the speed closed-loop system with ADRC and IADRC, respectively.
Jmse 11 01751 g012
Figure 13. PMSM experimental platform.
Figure 13. PMSM experimental platform.
Jmse 11 01751 g013
Figure 14. PMSM active disturbance rejection control system for ship propulsion.
Figure 14. PMSM active disturbance rejection control system for ship propulsion.
Jmse 11 01751 g014
Figure 15. Dynamic performances at 900 rpm under different load torque. (a) 1A. (b) 2A. (c) 3A. (d) Motor starting performance under 1A load.
Figure 15. Dynamic performances at 900 rpm under different load torque. (a) 1A. (b) 2A. (c) 3A. (d) Motor starting performance under 1A load.
Jmse 11 01751 g015aJmse 11 01751 g015b
Figure 16. Dynamic performances at 300 rpm under different load torque. (a) 2A. (b) 3A.
Figure 16. Dynamic performances at 300 rpm under different load torque. (a) 2A. (b) 3A.
Jmse 11 01751 g016
Figure 17. Performances at 300 rpm with current measurement error. (a) Phase current. (b) Speed and current.
Figure 17. Performances at 300 rpm with current measurement error. (a) Phase current. (b) Speed and current.
Jmse 11 01751 g017
Figure 18. Dynamic and steady-state performances at 300 rpm under different load torque. (a) Anti-disturbance performance under 3A load torque. (b) Phase current under 0.8A load torque.
Figure 18. Dynamic and steady-state performances at 300 rpm under different load torque. (a) Anti-disturbance performance under 3A load torque. (b) Phase current under 0.8A load torque.
Jmse 11 01751 g018
Figure 19. Experimental results of step-by-step starting of propulsion motor under three control methods.
Figure 19. Experimental results of step-by-step starting of propulsion motor under three control methods.
Jmse 11 01751 g019
Figure 20. Experimental results of direct starting of propulsion motor under three control methods.
Figure 20. Experimental results of direct starting of propulsion motor under three control methods.
Jmse 11 01751 g020
Figure 21. Anti-disturbance performance of three methods when propeller load torque changes suddenly.
Figure 21. Anti-disturbance performance of three methods when propeller load torque changes suddenly.
Jmse 11 01751 g021
Figure 22. Ship propulsion motor simulation results. (a) Start-up simulation under 50% load. (b) Disturbance-suppression performance simulation.
Figure 22. Ship propulsion motor simulation results. (a) Start-up simulation under 50% load. (b) Disturbance-suppression performance simulation.
Jmse 11 01751 g022
Table 1. Parameters of a ship.
Table 1. Parameters of a ship.
ParametersValue
Ship tonnage m92 t
Propeller diameter D0.9 m
Number of blades4
Propeller screw pitch ratio H/D0.965
Propeller disc square ratio A/Ad0.45
Thrust deduction coefficient t0.145
Wake coefficient w0.157
Resistance coefficient c694.2
Table 2. Main parameters of PMSM.
Table 2. Main parameters of PMSM.
ParametersValue
motor pole pairs np4
stator resistance Rs0.36 Ω
stator inductor Ls0.2 mH
permanent magnet flux λr0.0064 Wb
inverter input voltage Vdc30 V
motor rated current Ir7 A
Table 3. Parameters of the controller.
Table 3. Parameters of the controller.
ParametersADRC/IADRC
r50,000
h0.0001
λ0.8
k1400,000
k24000
k 2 20,000
k 2 1000
b01.36 × 10−7
ωo4000 rad/s
β112,000 rad/s
β24.8 × 107 rad/s
β36.4 × 1010 rad/s
Table 4. Dynamic performance comparison of the PMSM (Figure 15).
Table 4. Dynamic performance comparison of the PMSM (Figure 15).
Performance MetricsPIADRCIADRC
Δωm (rpm)1A21.220.210.4
2A49.64824.2
3A93.59248
Δt (s)1A0.140.150.04
2A0.20.220.05
3A0.230.260.07
Start-up time t(s)-0.250.09
Table 5. Dynamic performance comparison of the PMSM (Figure 16).
Table 5. Dynamic performance comparison of the PMSM (Figure 16).
Performance MetricsADRCADRC + iqIADRCIADRC + iq
Δωm (rpm)2A4724207
3A89434115
Δt (s)2A0.210.090.040.02
3A0.250.120.060.03
Table 6. Harmonic analysis of Figure 17b.
Table 6. Harmonic analysis of Figure 17b.
Harmonics of Speed1st2nd6th12thTHD
IADRC1.22%0.37%0.03%0.02%1.45%
IADRC + iq2.35%0.43%0.28%0.04%2.48%
Harmonics of Current2nd3rd5th7thTHD
IADRC2.68%0.65%1.44%1.36%3.59%
IADRC + iq6.27%0.95%5.13%4.92%8.92%
Table 7. Harmonic analysis of the phase current in Figure 18b.
Table 7. Harmonic analysis of the phase current in Figure 18b.
FFT of Phase Current3rd5th7th11th
IADRC + iq0.95%5.13%4.92%0.18%
IFADRC + iq0.79%1.51%1.44%0.05%
Table 8. Motor dynamic response performance in Figure 21.
Table 8. Motor dynamic response performance in Figure 21.
Performance MetricsPIADRCIFADRC + iq
Δωm (rpm)up27827964
drop18615858
total464437122
Δt (s)up0.240.260.10
drop0.340.410.15
Table 9. Parameters of the PMSM in simulation.
Table 9. Parameters of the PMSM in simulation.
ParametersValue
motor pole pairs np12
stator resistance Rs0.36 Ω
stator inductor Ls3.074 H
permanent magnet flux λr0.305 Wb
Moment of inertia J0.8051 kg·m2
Table 10. Simulation results.
Table 10. Simulation results.
Performance MetricsADRCADRC + iqIADRCIFADRC + iq
Δωm (rpm)50%-rated load1751367142
100%-rated load33624315096
Start-up time (50% load) t(s)0.12-0.26-
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

Guo, H.; Xiang, T.; Liu, Y.; Zhang, Q.; Liu, S.; Guan, B. Active Disturbance Rejection Control Method for Marine Permanent-Magnet Propulsion Motor Based on Improved ESO and Nonlinear Switching Function. J. Mar. Sci. Eng. 2023, 11, 1751. https://doi.org/10.3390/jmse11091751

AMA Style

Guo H, Xiang T, Liu Y, Zhang Q, Liu S, Guan B. Active Disturbance Rejection Control Method for Marine Permanent-Magnet Propulsion Motor Based on Improved ESO and Nonlinear Switching Function. Journal of Marine Science and Engineering. 2023; 11(9):1751. https://doi.org/10.3390/jmse11091751

Chicago/Turabian Style

Guo, Haohao, Tianxiang Xiang, Yancheng Liu, Qiaofen Zhang, Siyuan Liu, and Boyang Guan. 2023. "Active Disturbance Rejection Control Method for Marine Permanent-Magnet Propulsion Motor Based on Improved ESO and Nonlinear Switching Function" Journal of Marine Science and Engineering 11, no. 9: 1751. https://doi.org/10.3390/jmse11091751

APA Style

Guo, H., Xiang, T., Liu, Y., Zhang, Q., Liu, S., & Guan, B. (2023). Active Disturbance Rejection Control Method for Marine Permanent-Magnet Propulsion Motor Based on Improved ESO and Nonlinear Switching Function. Journal of Marine Science and Engineering, 11(9), 1751. https://doi.org/10.3390/jmse11091751

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