Marine diesel engines are widely used in the domain of ship propulsion [1, 2]. The most crucial task is to have a reliable control system to regulate its speed for safe and efficient operation under the inherent instabilities and disturbances, coupled with the unpredictable external environment [3]. With the development of control technology and the requirement for dynamic positioning systems, several control methods have been already applied [4]. These methods can be mainly divided into two categories, namely, the model-based control methods and the optimisation-based black box methods [5]. The model-based control methods usually include optimal control [6], adaptive control [7], sliding mode control [8] and so on. These methods need to obtain an accurate mathematical model. The performance of model-based control methods greatly depends on the accuracy of the model, but it is hard to achieve in practice. By contrast, the optimisation-based black box methods usually include the classic proportional–integral–derivative (PID) control [9], fuzzy control [10], active disturbance rejection control (ADRC) [11], modeless control and so on. These methods can achieve a good control effect without an accurate mathematical system model.
For the speed controller of a marine engine, the main task is to deal with the uncertain disturbances and the adjustment of system parameters. To handle uncertain disturbances, the ADRC method has introduced an extended state observer (ESO) in the feedback loop to compensate for the unmeasured states and the disturbances at the same time [12,13,14]. In addition to ADRC, the ESO can also be combined with other control methods to deal with uncertain disturbances and be simplified into a linear ESO (LESO) [15,16] for easier parameter adjustment [17]. Although simplification of both the control structure and the bandwidth tuning method facilitates the application of the LESO, non-linear ESO has been rarely studied [18, 19]. Therefore, it is essential to further study the non-linear ESO.
On the other hand, due to its advantages in terms of faster convergence rates, higher accuracies and better disturbance rejection properties, finite-time stability is also an important demand for control systems, especially in practice [20,21]. However, as far as our knowledge goes, the finite-time stabilisation of ADRC has never been attempted in the literature before.
In this paper, a novel ADRC law is proposed for a class of uncertain non-linear systems. The considered system uncertainties arise from the possibly unknown system dynamics, external disturbance and the parameter mismatch of control. To stay close to the basic concept of ADRC, the total uncertainties are treated as an extended state of the plant and are then estimated via non-linear ESO; finally, we compensate for these in the control action, in real time. Based on the output of the ESO, a finite-time sliding mode controller is designed [22]. A switching function is also defined to guarantee that the proposed control law is continuously differentiable. Theoretical analysis based on the Lyapunov stability theory show that the origin of the closed loop is semi-global finite-time stable.
The remainder of this paper is organised as follows. The mathematic model for the speed of a marine diesel engine is introduced in the section ‘Marine diesel engine speed model’. In the section titled ‘ADRC design’, we use an inverse hyperbolic sine function to design a non-linear ESO for the diesel engine speed model. The new non-linear ESO method, which has few adjustment parameters, is proposed to achieve precise control of marine diesel engine speed under uncertain disturbances. We design a finite-time sliding mode controller for ADRC (FT-SADRC), and it is proved that the controller can converge in finite time. In the section ‘Simulation experiments’, several experiments are carried out under the conditions of random load disturbances, manoeuvring the operations respectively, and the performance of the new ADRC method is verified in comparison to the linear ADRC (LADRC) method. In the ‘Conclusion’ section, conclusions and future research directions are presented.
In this paper, the electronically controlled diesel engine 7RT-Flex60C is studied. The maximum sustained power is 16520 kW, and the rated speed is 114 r/min. The relationship between the input (fuel injection) and the output (speed) of the diesel engine is obtained from the mean value engine model [23], which is adopted in this article.
We define the following functions:
The load perturbation of the diesel engine is the disturbance caused by complex sea conditions. The ADRC system with load disturbance characteristics is shown in Fig. 1.
According to the working principle of Fig. 1, the marine diesel engine propulsion control system with load disturbance can be expressed as follows:
The ADRC approach consists of three modules: tracking differentiator (TD), ESO and the non-linear state error feedback (NLSEF) control [24]. The TD generates a transient trajectory to avoid the set point jump. The system’s total disturbance is estimated by the ESO. The NLSEF is a non-linear PID controller.
The modules’ design is independent. Each module can be designed separately according to the controlled object and combined together to acquire the ADRC controller. For the system represented by Eq. (4), the TD in discrete form is as follows:
In order to improve the control performance, we apply the sliding mode error feedback control law to the ADRC. We define the following:
Due to the expansion of the observer convergence in Eq. (7), i.e.
And then, we introduce the following lemma, which is useful in the paper.
Consider the following non-linear systems [28]:
Consider the non-linear system, Eq. (12); suppose that there exist any real numbers
In the following section, the finite-time stability theory is used to analyse the convergence of the sliding mode error feedback control.
Based on the finite-time convergence stability theory and the sliding-mode control theory, we can design a controller for the finite-time convergence of the guidance system represented by Eq. (11) and give the following result.
Consider the system represented by Eq. (4); we design the system control quantity in the form of Eq. (15).
The sliding surface of the system can be chosen as follows [29]:
According to the reaching condition of sliding-mode control,
From Eqs (15) and (16), we have
Equation (17) shows that the system meets the reaching condition, i.e. the system’s state can reach the sliding mode surface.
We take the derivative of Eq. (15) and then connect Eq. (11) and Eq. (16):
We selected the non-linear control law as follows:
Therefore, the control input can be obtained as follows:
Next, we analyse the convergent property of the system based on the finite-time convergence stability theory. The motion of the state can be divided into two stages: (i) the approaching stage; and (ii) the sliding state.
For the approaching stage, the Lyapunov function can be chosen as follows:
The derivative of
From Eqs (11) and (22), we have
From Eqs (21) and (23), we have
Consider the operation of the control system,
Therefore,
According to Lemma 1, the state will converge to the sliding surface in a finite time. The convergent time satisfies the following:
After arriving at the sliding mode surface, the system will move along the surface until it converges to zero. So, the state meets the following requirements:
Consider that
And hence, from Eqs (27) and (28), we have
We define that
Equation (29) can be written as
Choose the Lyapunov function for Eq. (31) as follows:
And then, we have the following time derivative of
From Eqs (30) and (31), we have
According to Lemma 1, the convergent time of the system represented by Eq. (31) satisfies the following:
Therefore, the whole convergent time of the control system can be given as follows:
So, the control system is finite-time convergent.
From Eq. (36), we can know that the convergent time relates to the parameters
In this paper, we use the sine function
In order to verify the effectiveness of the proposed FT-SADRC, the RT-Flex60C low-speed two-stroke marine diesel engine, which is installed on the container ship of Orient Overseas International Limited (OOIL), is chosen for the simulation experiments. The main parameters of the RT-Flex60C engine are shown in Table 1.
Parameters of the RT-Flex60C engine
Number of cylinders | 7 |
Stroke | 2250 mm |
Crank link ratio | 0.489 |
Compression ratio | 18.4 |
Stroke number | 2 |
Cylinder diameter | 600 mm |
Rated speed | 114 r/min |
Rated power | 16520 kW |
Rated torque | 1390 kN·m |
Specific fuel consumption | 177 g/kWh |
To make a comprehensive analysis of the control performance of the proposed FT-SADRC method, a standard LADRC method is used to conduct the comparison. The robustness of the control algorithm under the conditions of a random disturbance load and sudden load dumping needs to be verified.
During the actual navigation, the wind and the waves have an effect on the propeller load, which are represented as random disturbances. In this test, random disturbances with an amplitude of 25 × 103 sin(
Gao [17] proposed a tuning method based on the observer bandwidth for LADRC. In this test, the bandwidth of the observer
The FT-SADRC parameter is selected as follows:
The control results of the proposed FT-SADRC and the LADRC with random disturbances are shown in Figs 2 and 3. The following points can be noted:
Although the control performances of both the methods are good, the tracking error of the FT-SADRC is much less than that of the LADRC. Furthermore, the tracking error of the FT-SADRC and the LADRC are 0.559 and 1.785r/min, respectively. In Fig. 3, we find that control input fluctuation is smaller for the FT-SADRC, which can save energy and reduce actuator wear.
This section describes the simulation analysis of the ship manoeuver navigation when entering and leaving a port. In this test, random disturbances with an amplitude of 25 × 103 sin(
For the LADRC [27], the bandwidth of the observer
It can be observed from the results in Fig. 4 that the maximum speed of the proposed FT-SADRC method is less than that of the LADRC method, and the time required to restore normal speed is also less than that of the LADRC method.
Moreover, it can be seen from Fig. 5 that the control input of FT-SADRC requires only about 0.5 s from the sudden load dumping to restore stability, while the LADRC requires >1 s. Therefore, the adjustment speed of the control input of the proposed method is faster than that of the LADRC.
In the context of the marine diesel engine speed control problem under the random load disturbances and unknown uncertainties, the sliding-mode control method is put forward by combining a new ESO (based on the inverse hyperbolic sine function). Although the new observer is non-linear, the adjustment parameters are the same as for linear ESO, which reduces the difficulty in the application of the non-linear observer. The control law is constructed by sliding-mode control for ADRC, and the method can achieve finite-time convergence and improve the convergence speed and control precision. From the simulation experiment using the model of a large, low-speed two-stroke marine diesel engine RT-Flex60c, we can observe that the proposed method has more advantages than the LADRC method, such as control accuracy and disturbance rejection (both random disturbance and sudden load dumping disturbance).
Thus, the main contributions of this article are summarised as follows. (1) In theory, a new ESO and finite-time convergence control law is proposed, which makes full use of the known non-linear term in the diesel engine model and further improves the adaptability of the controller. The proposed method can be applied to control issues of other systems that have known non-linear terms. (2) In practice, the proposed method can maintain smoother, quicker regulation, with fewer tuning parameters of the diesel engine speed under different external and system disturbances and has more practical significance compared with LADRC. Actually, the fuel injection signal of a marine diesel engine is usually maintained at several fixed positions for safety and economy. Therefore, the proposed method under a certain output constraint could be considered further.