Research Article  Open Access
Xuemin Li, Yufei Liu, Haoyu Shu, Runzhi Wang, Yunlong Yang, Chunyue Feng, "Disturbance ObserverBased Discrete SlidingMode Control for a Marine Diesel Engine with Variable Sampling Control Technique", Journal of Control Science and Engineering, vol. 2020, Article ID 2456850, 17 pages, 2020. https://doi.org/10.1155/2020/2456850
Disturbance ObserverBased Discrete SlidingMode Control for a Marine Diesel Engine with Variable Sampling Control Technique
Abstract
This paper proposes a disturbance observerbased discrete slidingmode control scheme with the variable sampling rate control for the marine diesel engine speed control in the presence of system uncertainties and disturbances. Initially, a slidingmode controller based on the fast power reaching law is employed, which has a good dynamic quality of the arrival stage and can suppress chattering. To satisfy the practical requirements in the digital controller and the crank anglebased fuel injection in engine speed control, the proposed method is discretized under the variable sampling rate condition. A disturbance observer based on the secondorder slidingmode control is designed to compensate the system uncertainties and disturbances, by doing such the requirement of the parameters of the slidingmode controller to be reduced significantly. In addition, a cylinderbycylinder mean value engine model (MVEM) is built by restructuring the combustion torque model, based on which numerical simulations are carried out by comparing the proposed method with PID and the extended state observer (ESO)based sliding mode controllers. The common operation situations of the marine diesel engines are taken into account, including starting process, acceleration and deceleration, load variation, and varied propulsion system parameters. The results demonstrate that the proposed disturbance observerbased slidingmode controller has prominent control performance and strong robustness.
1. Introduction
In the ship domain, diesel engines have been widely used as marine main engines for propulsion [1]. These engines are usually controlled to work at the specified rotating speed to keep the ship sailing at the appropriate speed, although this adds unknown load disturbances to the main engines [2]. The speed fluctuation and overshoot will cause harm to the whole propulsion system and even paralyze the entire system. Hence, for ship propulsion, good speed control performance of the marine main engines must be guaranteed even in the face of the strong nonlinearity in diesel engines and harsh marine environments. In addition, the emergence of new technologies in recent years also presents new challenges for speed control. Some new techniques greatly introduce additional dynamic nonlinearity for diesel engines, and some put forward higher requirements for speed control [3]. Therefore, it is necessary to design an accurate and robust speed controller for marine diesel engines.
Since the implementation of electroniccontrolled fuel control system, proportionintegrationdifferentiation (PID) has always been the most widely used algorithm for the speed control of diesel engines because of its simplicity and reliability [4]. Nevertheless, the PID controller can hardly maintain good control performance in all working conditions for the nonlinear system [5]. To overcome such problem, the PID controller is often studied to combine with other algorithms, such as fuzzy logic control [5], genetic algorithm [6], H∞ algorithm [7], and model predictive control (MPC) [8]. These composite algorithms improve the control performance by optimizing the PID parameters in real time. Meanwhile, some scholars try to use some other control methods to replace PID. The active disturbance rejection control (ADRC) is introduced to improve the performance of marine main engine under complex operating conditions in [9–11]. Ouladsine et al. established a neural controller with pollution constraints which realizes the multivariable control of engine speed and emission [12]. Slidingmode control (SMC), known as a robust control method [13–15], also has been studied in engine speed control. A simple SMC is designed by Zhang for speed control of generators, but the control performance is not satisfactory [16]. Yuan and his collaborators propose discrete slidingmode variable structure control in the literature [17] and verify the control performance for diesel engine through simulation and experiment, and the results show that its control effect is better than the traditional PID. Afterwards, multisliding surface control is proposed in his literature [18], and different sliding surfaces are selected in the acceleration and steadystate conditions to achieve better control performance in the various working conditions of the marine diesel engine. However, both of these methods rely on the accurate model. However, it should be noted that the simplified and accurate diesel engine model is hard to build, which severely limits the use of these controllers. In [19], the diesel engine model is approximated by the RBF neural network, based on which a slidingmode controller is developed. Khan [19] and Li [18] use the supertwisting scheme to control engine speed, where the accurate diesel engine model is no longer required but with the sacrifice of robustness.
It is well known that the traditional SMC requires the switching gain larger than the total disturbance, which sacrifices its control performance. Meanwhile, the chattering problem is hard to avoid [20]. Faced with this situation, one way is used to compensate the system disturbance by means of the disturbance observer. By doing such, the switch gain is only required to be larger than the estimated disturbance error rather than the total disturbance. On the basis of such idea, Wei proposes a controller via combining the disturbance observer and terminal SMC for a class of multipleinputmultipleoutput continuous nonlinear systems subject to disturbances [21]. In [22], a disturbance observerbased SMC is designed for a class of fractionalorder nonlinear systems. By combining the disturbance observer with SMC for the staringmode spacecraft, the stability of the spacecraft platform is improved and the chattering is reduced [23]. Motived by the statements above, in this study, the disturbance observerbased SMC is employed to regulate the marine diesel engine speed.
There is another issue should be considered in engine speed control, although the output speed of the diesel engine is continuous, its working process is not continuous with inherent eventtriggered property [11, 24–26], making the output of the controller can only be executed once per cylinder per cycle. For this reason, the disturbance observerbased SMC employed in this paper should be discrete. To avoid the chatting phenomena caused by highfrequency switching, the fast power reaching law is adopted by referring to [27]. The disturbance observer is structured based on a secondorder slidingmode algorithm which can significantly reduce the chattering in the discrete system. Considering that the controller output can only be executed at a fixed segment on the crankangle domain, variable sampling rate control is adopted instead of traditional timing control, which guarantees that the controller calculates once per working cycle for each cylinder at any engine speed (refer to [11] and the references therein for more details). In addition, a cylinderbycylinder mean value engine model (MVEM) is built for the validation of the proposed controller. Compared with the conventional MVEM, the adopted engine model is closer to the actual diesel engine, making the verification of the control algorithm more practical and reasonable.
The rest of the paper is organized as follows. First, the cylinderbycylinder MVEM is introduced. Then, a simplified discrete mathematical model of the diesel engine is described, based on which a discrete slidingmode controller and a slidingmode disturbance observer are designed. After that, based on the proposed cylinderbycylinder MVEM, simulation verification is carried out with vast comparative results to demonstrate the superior control performance of the proposed method. Finally, the conclusion about the whole work is summarized, and further research of this paper is discussed.
2. CylinderbyCylinder Mean Value Model of Marine Diesel Engine
Diesel engine is a complicated system with thermodynamic system, hydraulic system, electronic system, rotary and reciprocation machinery system, etc. Hence, it is difficult to establish an accurate mathematical model. As a compromise, MVEMs composed with primary empirical formula are widely used in correlation control research because of their easier implementation and higher calculating speed compared with detailed models [28]. Besides, it is very convenient to analyse the whole system and design qualified controllers with its intuitive system equations. However, as mentioned above, the actual working process of the diesel engine is periodic and cylinderbycylinder, while the traditional MVEMs only consider the average effect of the periodic work. For a fourstroke diesel engine, its working process includes the following features:(1)The working process is cylinderbycylinder(2)Each cylinder completes a cycle every two revolutions, including the suction stroke, the compression stroke, the power stroke, and the exhaust stroke(3)The control input, i.e., fuel injection quantity per cylinder per cycle, produces output only after a cylinder has completed its power stroke
It can be seen that no matter for the individual cylinder or the whole diesel engine, the working process is discontinuous. Therefore, although the study of continuous control algorithm based on continuous MVEMs can realise the rigorous theoretical analysis, there are still some differences with the actual control of diesel engine. To improve the situation, a cylinderbycylinder MVEM of a marine diesel engine is presented in this study, which can reflect the practical working characteristics mentioned above. This cylinderbycylinder engine model is based on the MVEM with the remodel of the combustion torque model which is the most critical part of the diesel engine model. Figure 1 shows the main structure of the developed model.
2.1. Mean Value Engine Model
As shown in Figure 1, the engine model is mainly composed of the engine body (cylinders, crankshaft, etc.), turbocharger (turbine and compressor), exhaust manifold, and intake manifold. Except for the indicated torque in the engine body, the rest parts are modelled according to the empirical formula as MVEMs.
The main function of a diesel engine is to output the effective torque through the crankshaft and then drive the load rotation. Here, the loads are the gearbox, the spindle system, and its connected propeller. For the sake of analysis, assume that there is only one propeller and that the gear box has a reduction ratio of 1. Thus, the rotational dynamics of the total propulsion system including the diesel engine and the loads can be given aswhere is the engine speed, i.e., the rotation speed of crankshaft, is the total rotational inertia of the moving parts of the propulsion system, is the total frictional loss torque, is the load torque, and and denote the total indicated torque and the total pump loss torque for all cylinders, respectively.
For MVEMs, and are usually calculated by empirical formulas:where is the fuel injection quantity per cylinder per cycle, is the fuel low calorific value, is the gross indicated efficiency, is the number of cylinders, is the total displacement of the engine, and and , respectively, denote the pressure in intake manifold and exhaust manifold. One can see that only cyclic averages are obtained. In this paper, they are simulated via the cylinderbycylinder combustion torque model which will be described in the next subsection.
When a diesel engine runs, friction loss is caused by all moving parts, especially crankshaft and pistons. It is generally believed that these losses are mainly related to engine speed. Here, the frictional loss torque is assumed to be a quadratic polynomial of engine speed as follows [29]:where , , and are the experimental coefficients.
Ignoring the relatively small friction loss of the gearbox and the spindle system, the load torque is considered to be mainly from the propeller torque which can be approximated bywhere is a coefficient depending on the propeller torque coefficient , the density of water , and the propeller diameter . The relationship between them can be described by as [30]. For a propeller, is invariable, but varies with the advance coefficient and propeller pitch. So, is also varying. Here, two form variations will be taken into account during controller validation as follows.
Remark 1. As equation (2) shows, the gross indicated efficiency is an important index for diesel engines. Meanwhile, it is affected by many factors among which the excess air coefficient and engine speed are the most important. Referring to [31], it can be described aswhere are coefficients. If an accurate indicated torque is desired, the excess air coefficient is also needed to be accurately modelled. For this purpose, the turbocharger model, exhaust manifold model, and intake manifold model should be considered well. The air mass flow into the cylinders is calculated according to the air pressure and temperature in the intake manifold and the volumetric efficiency which is a function of the pressure and the engine speed [29]. The turbocharger model mainly provides boundary conditions and inlet gas flow to the intake manifold. The exhaust manifold model determines the rotation speed and working state of the turbocharger. For simplicity, the detailed modelling approach is not described here and more details are available in [32].
2.2. CylinderbyCylinder Combustion Torque Model
As mentioned above, only averages of and are calculated by empirical formulas for MVEMs. It is essential that the torque applied to the crankshaft by the cylinder pressure should be through the crank link mechanism. The force diagram of a single crank link mechanism is shown in Figure 2.
From the figure, the following relationship formulas can be obtained:where is the connecting rod length, is the crank arm length, and is the piston diameter. It should be noted that denotes the tilt angle of the connecting rod and denotes the angle at which the crankshaft rotates from the topdeadcenter (TDC). Meanwhile, it is supposed that , , and mean that the crankshaft is at the TDC of the compression phase. According to equation (7), the torque produced by a certain cylinder on the crankshaft can be obtained by
The model of the cylinder pressure is the core of the cylinderbycylinder engine model. According to the working process of the individual cylinder, it is divided into three stages to simulate for a working cycle. The first stage is the working stage, where only the intake valve is open. The second stage is the most complex stage, where both the intake and exhaust valves are closed. The third stage refers to the working stage where only the intake valve is open. In the first and third stages, is considered to be approached to the pressure in the intake manifold model and exhaust manifold model, respectively. The second stage includes the compression phase and the power phase, and the cylinder pressure can be treated as the interpolation of that in two ideal processes, namely, the pure compression process and ideal combustion process.
The pure compression process refers to the second stage process with no fuel supply, and the cylinder pressure can be measured by the motoring test method. For simulation, it can be regarded as a polytropic process. The pressure and temperature can be given aswhere , , and , respectively, represent the cylinder pressure, temperature, and volume at the intake valve closing, denotes the instantaneous volume of the cylinder, and is the polytropic exponent.
The ideal combustion process refers to the second stage process with a negligible time from the start of combustion to the end of combustion. Before combustion, the pressure and temperature of this ideal process are the same as the pure compression process, and after combustion, they can be modelled as another polytropic process as follows:where is the polytropic exponent and , , and denote the cylinder pressure, temperature, and volume of the ideal combustion process at the end of combustion, respectively. According to the description, the combustion process can be regarded as a constant volume combustion, and the and are determined bywhere denotes the at the start of combustion, is the specific heat at constant volume of mixture in the cylinder. And it is important to note that keeps the value at the start of combustion during the whole second stage and is updated only once in the next cycle. is the fuel conversion efficiency, and it is considered to be proportional to the gross indicated efficiency of the MVEM. That is,
The pressure ratio proposed by Matekunas is defined as [33]
Eriksson has found that the pressure ratio curve is similar to the proportion of fuel already burned [27]. If represents the pressure ratio when it has completely combusted at , which can be given aswhere is the percentage of fuel that has been burned which can be described by the wellknown Vibe function. According to the ideal combustion process, one can get
Combining equations (16)–(18), it yields to
As known, until the start of combustion and since the end of combustion, which leads to and , respectively [34].
In summary, the cylinder pressure of a working cycle can be described as
For the blowdown phase and the transitions between these three stages, the pressure can be calculated by interpolation using a cosine function.
Above all, for a sixcylinder diesel engine, the total torque generated by the cylinders can be obtained as
Combined with equation (21), one can see the in MVEMs is the cyclic average of . Therefore, the rotational dynamics (1) can be rewritten as
Thus, together with other parts of the original MVEM, the new cylinderbycylinder MVEM is formed. Some simulation results of this model are presented in Figures 3 and 4.
(a)
(b)
The cylinder pressure based on equation (19) is shown in Figure 3, and the three stages described above can be observed. This cylinder pressure produces cyclical fluctuation for the engine speed as shown in Figure 4, which cannot be given by the MVEMs. Moreover, by choosing the appropriate coefficient in equation (15), the speed responses are also similar except for the cyclic fluctuation when the fuel injection of the two models is the same. As can be seen from the acceleration at 3 s, the cylinderbycylinder model does not respond to the change of the fuel injection until the next cycle. This also causes the delay in the acceleration process. All these make the model closer to the actual working process of the diesel engine.
Remark 2. It should be noted that the original MVEM is borrowed from the work [32], which has been verified by extensive experiments. Hence, as shown in Figure 4, the mean effect of engine speed in the proposed engine model is almost the same as that in the original MVEM, which can be regarded as the verification of the proposed engine model.
3. Disturbance ObserverBased Discrete SlidingMode Controller
It can be seen from the model description that the diesel engine is a discontinuous system, and the control input will only be valid after the current working cylinder completes the power stroke. Thus, the general used constant sampling rate control will have problems; hence, a variable sampling rate control is employed in this study. The controller calculates only at a few fixed crank angles, which means that the sampling interval varies with the engine speed. Therefore, in this section, the control problem is simplified and discretized with the variable sampling control theory, based on which the controller and observer are designed.
3.1. Control Problem Simplification
Although the cylinderbycylinder model is closer to the reality engine, it is not convenient for the design of the controller. On the contrary, the MVEM is not realistic but simple in form, which can be regarded as a simplification of the controlled object. According to equations (1), (2), (4), and (5), the control problem equation can be described aswhere
Obviously, is not accurate and difficult to obtain for the actual diesel engine or cylinderbycylinder model. However, its variation range is limited, so it is regarded as a known constant with some unknown deviations. For the actual situation, both and can be measured experimentally, but there is always some offset. Besides, both of them have all kinds of interference in practical work, especially is strongly affected by the complex marine environment. As a result, some parts of are also unknown and unmodelable. Assume that is also unknown and treat all these unknown parts as unknown disturbance , given state variable and control input , and equation (23) can be rewritten aswhere is the known part of and .
Equation (26) can be discretized with a variable sampling interval as follows:where is the variable sampling interval which is related to the engine speed. This function can be regarded as a simplified discrete mathematical model of the diesel engine and the control problem equation when designing the controller in the below section.
3.2. Discrete SlidingMode Controller
Most slidingmode algorithms achieve their control objectives and maintain good robustness by fast switching, resulting in chattering often, especially when the control frequency is not fast enough. The control frequency of a diesel engine is variable and limited by the number of cylinders and the engine speed. Therefore, a suitable slidingmode controller is needed to suppress chattering. In this study, a slidingmode controller based on the reaching law proposed in [35] is used for the speed control of a diesel engine. The fast power reaching law is selected, which is helpful to improve the dynamic quality of the arrival stage and suppress chattering [36]. It can be described aswhere the coefficient and . When the initial value , the solution of this firstorder nonlinear equation can be expressed aswhich shows that equation (28) converges in finite time.
Considering the variable sampling rate control of the diesel engine, equation (28) can be written in discrete form asThat is,
Obviously,
is required to ensure the convergence of equation (30). Therefore, the convergence accuracy of the equation in discrete form is . It can be seen that the smaller coefficients are, the higher the steadystate convergence accuracy. However, it also reduces the transient control performance. This should be taken into account when selecting coefficients.
For system 27, let , and its discrete time derivative can be obtained aswhere denotes the reference speed. Combining equations (27) and (30), it yields to
Hence, the control law can be set to bewhere the disturbance is unknown. If the disturbance is ignored in equation (35), the system equation becomesthe convergence domain of which depends on the bound of the disturbance. This can have serious adverse effects on control performance. Therefore, an SMDO is designed to compensate for it.
3.3. SlidingMode Disturbance Observer
Treating as an additional state variable and assuming , the original firstorder system can be described as a secondorder system:where .
The disturbance observer can be given aswhich is often used by scholars [37–40]. and are the items that need to be designed.
Define observation error . Considering equations (37) and (38), an error system can be obtained as
Therefore, the design of disturbance observer becomes design and to make the error system converge. If , and the observer 29 becomes a firstorder system. Only needs to be designed so that converges to zero and , as in [40]. However, many scholars usually design both and , e.g., the extended state observer (ESO) proposed by Han [37], the nonlinear disturbance observer in [39], and the linear observer with optimization [38]. The purpose of all these observers is to achieve and , which is the same as a secondorder slidingmode controller (SOSMC). There are designs and analysis of secondorder sliding modes for continuous systems, such as twisting algorithm, suboptimal algorithm, and supertwisting algorithm [41, 42]. However, high control frequency is often required in these algorithms to ensure their control performance and robustness. In this study, the SOSMC controller described in Theorem 1 and the similar controller for discrete control are considered as follows.
Theorem 1 (see [43]). Considering the secondorder system,
Suppose there are and such that the disturbance and the variable structure control lawcan suppress the disturbance and make system (40) reach the origin after switching at most once in finite time.
Like many slidingmode controllers, it is difficult to switch accurately on the slidingmode surface when used for discrete control， which easily causes chattering. To overcome the problem, a similar control law for discrete control is proposed based on the principle of control law (41) by Han [29]. It is defined aswhere , , and are coefficients corresponding to the and in control law (41) and the sampling step of the discrete system, respectively. The two control laws converge system (28) in an almostuniform trajectory, if the step size is small enough and , , and .
Figure 5 gives the phase trajectories and the convergence trajectories of these two controllers where the corresponding parameters are the same except for . The phase trajectory A shows the secondorder slidingmode characteristics of controller (41). However, it can be seen from phase trajectory B and convergence trajectory B that the control performance becomes worse and chattering occurs easily when the control law is discontinuous. Meanwhile, controller (42) can still guarantee good control performance except for the smaller fluctuation due to discrete calculation. This fluctuation can also be weakened by appropriately increasing , and the convergence rate will be slower as a cost as convergence trajectory D shows.
According to the above description, controllers (41) and (42) can be used for continuous and discrete observer design, respectively. In this paper, the observer is proposed aswhere is a coefficient so that will change with the calculated step . In order to reduce the overshoot and steadystate fluctuation, a part of the convergence rate needs to be sacrificed when selecting the parameters of . Here, is introduced to improve the convergence rate. The error system can be rewritten as
The can be taken as a new disturbance, and it is bounded in diesel engine systems. According to Theorem 1, there are parameters of to make error system (44) converge in finite time.
In conclusion, based on this slidingmode disturbance observer (SMDO), the designed discrete slidingmode controller (DSMC) can be described as
Combing with system function (27), it yields to
Comparing with equation (36) shows that the convergence accuracy will be improved if which is guaranteed by the SMDO.
(a)
(b)
4. Results and Discussion
To verify the effectiveness of the proposed SMDO and DSMC, simulation on the cylinderbycylinder diesel engine model is carried out and the results are analysed in this section. Moreover, the ESO proposed in Ref. [37] and a PID controller are employed to compare with the proposed observer and controller, respectively.
In Table 1, the primary parameters of the modelled diesel engine are given. Similar to the general operating conditions of a marine diesel engine, the simulations include the starting process, steady process, acceleration process, deceleration process, and loadchanging process. To simulate the starting process in a real engine, the engine is driven by a starting torque until it reaches the set speed 350 rpm, and then it is driven by combustion torque since the controller starts working. After the starting process, the reference speed is increased from 700 rpm to 1900 rpm by 400 rpm each time and then decreased by 600 rpm to 1300 rpm. A slope limit for the reference speed is employed to ensure the safety of the propulsion system as the actual engine runs [30], which is set as 400 rpm/s. When a ship is running in the sea, the propeller torque is affected periodically by waves, which can be simulated by periodic disturbance torque in 13 to 15 seconds. In harsh sea conditions, the propeller inevitably emerges and reenters the water, which would cause dramatic load changing. Given such situation, the load torque is discharged and increased by 50% at 9 s and 10 s to simulate the effect of emerging and reentering. Especially, varied inertia moment is considered to evaluate the robustness of the proposed controller. Detailed simulation results and analysis are presented below.

4.1. SlidingMode Disturbance Observer Performance
As described above, the ESO is employed for comparison, and the simulation results of the disturbance observers with optimized parameters are given in Figure 6 and 7.
(a)
(b)
(a)
(b)
(c)
(d)
Since equation (27) is the average and simplification of the cylinderbycylinder model, the observed disturbance is actually an average. However, the combustion torque of the cylinderbycylinder model varies periodically with the crank angle; hence, the actual disturbance torque of the model is difficult to calculate. Therefore, the comparison between the observed disturbance and the actual disturbance cannot be given. The following results of the two observers for the engine speed are shown in Figure 6. It can be seen that, the speed tracking performance of the ESO is better than that of the SMDO in all the compared situations. However, as the speed estimation error is the input of the observers, too fast speed following is sometimes at the expense of disturbance observation effect. It is obvious that the observation of the SMDO is better than that of the ESO at initial stage (as shown in the downside subplot of Figures 6 and 7(a) due to the faster convergence in the SMDO). More detailed disturbance observation comparison is given in Figure 7.
As can be seen from Figure 7(a), the observer starts operation at nearly 0.8 seconds where the engine speed reaches 350 rpm. Then, both observers converge quickly, but the ESO converges significantly slower and has a brief oscillation. After that, the two estimated results are close, as shown in Figure 7(c). However, the SMDO is significantly faster than the ESO in response to sudden loading and unloading as present in Figures 7(b) and 7(d), and this will have a great effect on the speed control which will be discussed later. Overall, the SMDO performs better than the ESO. A more intuitive comparison can be seen from the speed tracking effect of the DSMC based on the two observers as shown below.
4.2. Discrete SlidingMode Controller Performance
For providing better comparison, three controllers are tested; they are SMDObased DSMC, ESObased DSMC, and PID. All these controllers perform the calculation of control input at a preset crank angle before fuel injection occurs in the individual cylinder. In order to avoid introducing meaningless periodic fluctuations which is shown in Figure 4, the engine speed is averaged between two adjacent calculations. Moreover, the disturbances are added to the engine speed calculation and the load torque separately for simulating the actual situation more closely. Therefore, the average speed is fluctuating even in steadystate condition, which is more likely to cause chattering of SMC.
The overall results of engine speed and control signal response curves for the compared controllers are shown in Figure 8, and some detailed results are shown in Figures 9–13.
(a)
(b)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(e)
(f)
As can be seen from Figure 8, there is no significant chattering under DSMC whether it is based on the ESO or SMDO, and all the compared controllers can achieve satisfactory control performance under the most operation conditions. Detailed presentation and analysis are described below.
Case 1. Starting process
As mentioned above, the diesel engine model is firstly driven by a starting torque until the speed reaches to 350 rpm, and then the controller takes over the control and starts the closedloop control. Furthermore, to avoid overregulation, the desired speed is set to always be 50 rpm higher than the actual speed before the actual speed approaches the idle speed. The simulation results after the controllers start working in the starting process are shown in Figure 9. Visibly, the overshoot and stabilization time of the PID controller reach to 140 rpm and 1.2 seconds, while the overshoots of the DSMC based on the ESO and SMDO are only 31 rpm and 20 rpm, respectively, and the stabilization time is less than 1.0 seconds. Moreover, since the SMDO converges faster and does not oscillate like the ESO, the starting process of the SMDObased DSMC is significantly smoother. After the speed stabilizes at idle condition, the load torque is added at 2 s. During loading, the three controllers have similar control effects, while the SMDObased DSMC is slightly better. Hence, it can be concluded that the SMDObased DSMC has a better starting performance than PID controller and the ESObased DSMC.
Case 2. Acceleration and deceleration process
Figure 10 shows the acceleration and deceleration processes. It can be seen that DSMCs based on the ESO and SMDO have a better tracking effect in all processes, except that there are slight overshoots but almost negligible. Meanwhile, the two DSMC methods maintain similar performance during all three acceleration processes, while the performance of the PID controller becomes worse in the processes in Figures 10(a) and 10(c).
Case 3. Load variation process
Figures 11 and 12, respectively, show the engine speed responses of two kinds of load variation mentioned above. During the sudden loading and unloading, one can see that the overshoot of the ESObased DSMC is higher than that of the other two controllers due to the slow response of the ESO which is shown in Figure 7. The convergence time of the PID controller is almost twice than that of DSMCs. In contrast, the SMDObased DSMC has the best control performance with the minimum overshoot and the shortest stabilization time. Facing the periodic disturbance, the control performance of the SMDObased DSMC is still the best with only 31 rpm speed fluctuation, while the PID controller is the worst with 52 rpm speed fluctuation. The DSMC based on the ESO has slightly worse performance than that based on the SMDO, in which the speed fluctuation is 39 rpm. Therefore, it can be concluded that the SMDObased DSMC is more capable of handling the load variation because of the harsh marine environment.
Case 4. Robustness verification
In practice, a diesel engine may be used to different propulsion systems with different parameters such as inertia moment and torque coefficient. This requires the controller to have good robustness, otherwise the controller parameters need to be recalibrated for different propulsion systems. Here, the simulation is carried out to verify the robustness of the controllers, in which the inertia moment of the propulsion system is increased by 25% and the parameters of the controllers remain unchanged. Some simulation results are given in Figure 13.
As can be seen that, during the starting process, the performance of the PID controller is still not satisfactory and the performance of the ESObased DSMC gets worse with overshoot increasing to 51 rpm. During the first loading process, the control effect of the three controllers becomes better than their original conditions as shown in Figure 9 while DSMCs are obviously better than the PID controller. In Figures 13(c) and 13(d), the DSMCs based on the ESO and SMDO maintain the good control performance, but the PID controller gains the worse control performance. Due to the increase of rotational inertia, the control performance of the three controllers gets better when suddenly loading and unloading. However, the DSMC based on the SMDO still has the best performance, and the PID controller is still significantly worse than the DSMCs. Facing the periodic disturbance, the speed fluctuation of the DSMCs based on the ESO and SMDO is 25 rpm and 20 rpm, respectively, which means that, for the DSMCs, the control performance under such condition is better than the original condition as shown in Figure 12. On the contrary, the control effect of the PID controller becomes worse, and the speed fluctuation increases to 58 rpm.
Summarizing the above analysis and comparison, the DSMCs show superior performance over the PID controller. While ensuring the robustness, the chattering is almost eliminated. As can be seen from the above simulation results, the steadystate fluctuation is basically the same as the PID controller. The robustness and control effect of the SMDObased DSMC are better than that of the ESObased DSMC, which further indicates that the performance of the designed SMDO is better than that of the ESO. Therefore, control performance of the marine diesel engine speed system can be further improved by the proposed SMDObased DSMC.
5. Conclusions
In this article, a disturbance observerbased discrete slidingmode control scheme is proposed for the speed control of the marine diesel engine. The system uncertainties and external disturbances are estimated by the designed disturbance observer based on the secondorder slidingmode control. With such estimate to compensate the uncertainties and disturbances, a discretized slidingmode controller based on the fast power reaching law is employed as the speed control law. In order to verify the controller, a cylinderbycylinder MVEM of a marine diesel engine is established by remodelling the combustion torque model, which can simulate the periodic operation of the diesel engine. The numerical simulations are carried out on this model. To show the superiority of the proposed method, the PID controller and ESObased DSMC are employed for comparison. Moreover, the common operation situations of the marine diesel engines are taken into account, including starting process, acceleration and deceleration, load variation, and varied propulsion system parameters. The results show that the designed SMDO is more stable and has better tracking performance than the ESO, and the proposed SMDObased DSMC has better control performance and stronger robustness under various working conditions than the PID and ESObased DSMC. Therefore, the speed control performance of a marine diesel engine can be significantly improved by using the proposed controller.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
References
 N. Xiros, Robust Control of Diesel Ship Propulsion, Springer, London, UK, 2002.
 R. D. Geertsma, R. R. Negenborn, K. Visser, and J. J. Hopman, “Design and control of hybrid power and propulsion systems for smart ships: a review of developments,” Applied Energy, vol. 194, pp. 30–54, May 2017. View at: Publisher Site  Google Scholar
 S. Roy, O. P. Malik, and G. S. Hope, “Adaptive control of speed and equivalence ratio dynamics of a diesel driven powerplant,” IEEE Transactions on Energy Conversion, vol. 8, no. 1, pp. 13–19, 1993. View at: Publisher Site  Google Scholar
 S.H. Lee, J.S. Yim, J.H. Lee, and S.K. Sul, “Design of speed control loop of a variable speed diesel engine generator by electric governor,” in Proceedings of the 2008 IEEE Industry Applications Society Annual Meeting, Institute of Electrical and Electronics Engineers Inc, Edmonton, Canada, 2008. View at: Google Scholar
 D. J. McGowan, D. J. Morrow, and B. Fox, “Integrated governor control for a dieselgenerating set,” IEEE Transactions on Energy Conversion, vol. 21, no. 2, pp. 476–483, 2006. View at: Publisher Site  Google Scholar
 F. A. Mohamed and H. N. Koivo, “Diesel engine systems with genetic algorithm self tuning PID controller,” in Proceedings of the International Conference on Future Power Systems, Amsterdam, Netherlands, November 2005. View at: Google Scholar
 N. I. Xiros, “PID marine engine speed regulation under full load conditions for sensitivity H∞norm specifications against propeller disturbance,” Journal of Marine Engineering & Technology, vol. 3, no. 2, pp. 3–11, 2004. View at: Publisher Site  Google Scholar
 R. Wang, X. Li, J. Zhang et al., “Speed control for a marine diesel engine based on the combined linearnonlinear active disturbance rejection control,” Mathematical Problems in Engineering, vol. 2018, 18 pages, 2018. View at: Publisher Site  Google Scholar
 C. Wang, Y. Zhou, W. Pan, Y. Han, and F. Zhou, “The application of ADRC in the ship main engine speed controller based on genetic algorithm,” in Proceedings of the 2011 International Conference on Materials Science and Computing Science, vol. 327, pp. 17–22, Trans Tech Publications, Wuhan, China, 2011. View at: Google Scholar
 R. Wang, X. Li, J. Zhang et al., “Speed control for a marine diesel engine based on the combined linearnonlinear active disturbance rejection control,” Mathematical Problems in Engineering, vol. 2018, Article ID 7641862, 18 pages, 2018. View at: Publisher Site  Google Scholar
 R. Wang, X. Li, Y. Liu et al., “Variable sampling rate based active disturbance control for a marine diesel engine,” Electronics, vol. 8, no. 4, p. 370, 2019. View at: Publisher Site  Google Scholar
 M. Ouladsine, G. Bloch, and X. Dovifaaz, “Neural modelling and control of a Diesel engine with pollution constraints,” Journal of Intelligent & Robotic Systems, vol. 41, no. 23, pp. 157–171, OctNov 2004. View at: Google Scholar
 C. Mu and H. He, “Dynamic behavior of terminal sliding mode control,” IEEE Transactions on Industrial Electronics, vol. 65, no. 4, pp. 3480–3490, 2018. View at: Publisher Site  Google Scholar
 C. Mu, W. Xu, and C. Sun, On Switching Manifold Design for Terminal Sliding Mode Control, Elsevier, Amsterdam, Netherlands, 2016.
 C. Mu and C. Sun, A New Finite Time Convergence Condition for Supertwisting Observer Based on Lyapunov Analysis, Wiley, Hoboken, NJ, USA, 2015.
 G.C. Zhang, “Speedfrequency controller design based on sliding mode for marine dieselgenerator,” in Proceedings of the 2010 2nd WRI Global Congress on Intelligent Systems, vol. 2, pp. 31–34, IEEE Computer Society, Wuhan, China, December 2010. View at: Google Scholar
 Y. Yuan, Y. He, L. Cai, and X. Mao, “Discrete sliding mode variable structure control over the rotating speed of marine diesel engines,” Proceedings of the Institution of Mechanical Engineers, Part I: Journal of Systems and Control Engineering, vol. 231, no. 5, pp. 367–379, 2017. View at: Publisher Site  Google Scholar
 Y. Yuan, M. Zhang, Y. Chen, and X. Mao, “Multisliding surface control for the speed regulation system of ship diesel engines,” Transactions of the Institute of Measurement and Control, vol. 40, no. 1, pp. 22–34, 2018. View at: Publisher Site  Google Scholar
 J. Zhang, H. Lan, and Y. Sun, “Study on the speed control of diesel engine based on neural sliding mode control,” in Proceedings of the International Conference on Material and Manufacturing, ICMM 2011, vol. 299300, pp. 1190–1193, Trans Tech Publications, Jinzhou, Liaoning, China, 2011. View at: Google Scholar
 J. Yang, S. Li, and X. Yu, “Slidingmode control for systems with mismatched uncertainties via a disturbance observer,” IEEE Transactions on Industrial Electronics, vol. 60, no. 1, pp. 160–169, 2013. View at: Publisher Site  Google Scholar
 X. Wei and L. Guo, “Composite disturbanceobserverbased control and terminal sliding mode control for nonlinear systems with disturbances,” International Journal of Control, vol. 82, no. 6, pp. 1082–1098, 2009. View at: Publisher Site  Google Scholar
 S. Shao, M. Chen, and Q. Yang, “Sliding mode control for a class of fractionalorder nonlinear systems based on disturbance observer,” in Proceedings of the IEEE International Conference on Industrial Technology, ICIT 2016, pp. 1790–1795, Institute of Electrical and Electronics Engineers Inc, Taipei, Taiwan, May 2016. View at: Google Scholar
 S. N. Wu, X. Y. Sun, Z. W. Sun, and X. D. Wu, “Slidingmode control for staringmode spacecraft using a disturbance observer,” Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, vol. 224, no. 2, pp. 215–224, 2010. View at: Publisher Site  Google Scholar
 E. De Santis, M. D. Di Benedetto, and G. Pola, “Digital idle speed control of automotive engines: a safety problem for hybrid systems,” Nonlinear Analysis: Theory, Methods & Applications, vol. 65, no. 9, pp. 1705–1724, 2006. View at: Publisher Site  Google Scholar
 S. Yurkovich and M. Simpson, “Comparative analysis for idle speed control: a crankangle domain viewpoint,” in Proceedings of the 1997 American Control Conference, vol. 1, pp. 278–283, Albuquerque, Mexico, June 1997. View at: Google Scholar
 M. Mosleh and A. AlAli, “Discrete speed controller design of a marine diesel engine including sampling effects due to fuel injection,” Journal of Vibration and Control, vol. 8, no. 5, pp. 659–671, 2002. View at: Publisher Site  Google Scholar
 L. Eriksson, “Requirements for and a systematic method for identifying heatrelease model parameters,” in Proceedings of the 1998 SAE International Congress and Exposition, SAE International, Detroit, MI, USA, February 1998. View at: Google Scholar
 G. Theotokatos, C. Guan, H. Chen, and I. Lazakis, “Development of an extended mean value engine model for predicting the marine twostroke engine operation at varying settings,” Energy, vol. 143, pp. 533–545, Jan 15 2018. View at: Google Scholar
 J. Han, Active Disturbance Rejection Control Techniquethe Technique for Estimating and Compensating the Uncertainties, National Defense Industry Press, Beijing, China, 2008.
 P. Mizythras, E. Boulougouris, and G. Theotokatos, “Numerical study of propulsion system performance during ship acceleration,” Ocean Engineering, vol. 149, pp. 383–396, 2018. View at: Publisher Site  Google Scholar
 L. Guzzella and A. Amstutz, “Control of diesel engines,” IEEE Control Systems, vol. 18, no. 5, pp. 53–71, 1998. View at: Google Scholar
 J. Wahlstrom and L. Eriksson, “Modelling diesel engines with a variablegeometry turbocharger and exhaust gas recirculation by optimization of model parameters for capturing nonlinear system dynamics,” Proceedings of the Institution of Mechanical Engineers Part DJournal of Automobile Engineering, vol. 225, no. D7, pp. 960–986, 2011. View at: Publisher Site  Google Scholar
 F. A. Matekunas, “Engine combustion control with ignition timing by pressure ratio management,” 1986, US Patent. View at: Google Scholar
 R. Anjum, A. Yar, Q. Ahmed, and A. I. Bhatti, “Cyclic torque imbalance detection in gasoline engines using second order sliding mode,” in Proceedings of the 2nd IEEE Conference on Control Technology and Applications, CCTA 2018, pp. 1046–1051, Institute of Electrical and Electronics Engineers Inc, Copenhagen, Denmark, 2008. View at: Google Scholar
 W. Gao and J. C. Hung, “Variable structure control of nonlinear systems. A new approach,” IEEE Transactions on Industrial Electronics, vol. 40, no. 1, pp. 45–55, 1993. View at: Google Scholar
 S. Yu, X. Yu, B. Shirinzadeh, and Z. Man, “Continuous finitetime control for robotic manipulators with terminal sliding mode,” Automatica, vol. 41, no. 11, pp. 1957–1964, 2005. View at: Publisher Site  Google Scholar
 J. Han, “From PID to active disturbance rejection control,” IEEE Transactions on Industrial Electronics, vol. 56, no. 3, pp. 900–906, 2009. View at: Publisher Site  Google Scholar
 Z. Gao, “Scaling and bandwidthparameterization based controller tuning,” in Proceedings of the2003 American Control Conference, vol. 6, pp. 4989–4996, Institute of Electrical and Electronics Engineers Inc, Denver, CO, USA, 2003. View at: Google Scholar
 Y. B. Shtessel, I. A. Shkolnikov, and A. Levant, “Smooth secondorder sliding modes: missile guidance application,” Automatica, vol. 43, no. 8, pp. 1470–1476, 2007. View at: Publisher Site  Google Scholar
 W. Liu and P. Li, “Disturbance observerbased faulttolerant adaptive control for nonlinearly parameterized systems,” IEEE Transactions on Industrial Electronics, vol. 66, no. 11, pp. 8681–8691, Nov 2019. View at: Publisher Site  Google Scholar
 G. Bartolini, A. Ferrara, and E. Usai, “Chattering avoidance by secondorder sliding mode control,” IEEE Transactions on Automatic Control, vol. 43, no. 2, pp. 241–246, 1998. View at: Publisher Site  Google Scholar
 A. Levant, “Sliding order and sliding accuracy in sliding mode control,” International Journal of Control, vol. 58, no. 6, pp. 1247–1263, 1993. View at: Publisher Site  Google Scholar
 J. Han, “Invariability of time optimal feedback control,” Journal of Systems Science and Mathematical Sciences, vol. 25, no. 4, pp. 498–506, 2005. View at: Google Scholar
Copyright
Copyright © 2020 Xuemin Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.