An Intelligent Bio-Inspired Cooperative Decoupling Control Strategy for the Marine Boiler-Turbine System with a Novel Energy Dynamic Model

Abstract: This paper presents an intelligent bio-inspired cooperative decoupling control strategy (IBICDC) for the problems of modeling difficulties and strong coupling in the marine boiler-turbine system (MBTS). First, the model of the main steam pressure control loop is successfully constructed by introducing the Martin-Hou equation, which solves the modeling difficulty caused by the complexity of structure, operation mechanism, and operation conditions, as well as the characteristics of nonlinearity, parameter time-varying, and time-delay in the marine boiler (MB). According to the mathematic method of homeomorphic mapping relationship between the rotational speed and the kinetic energy in the marine steam turbine with propeller (MSTP) and the feedback linearization method, the nonlinear degree of the MSTP rotational speed control loop model is reduced and the infinite point of discontinuity in the rotational acceleration when the rotational speed close to 0 is eliminated. Then, the IBICDC inspired by the internal environment regulation mechanism of human body is applied to the strong coupling problem between the two control loops, namely, to eliminate the large value sudden change of the main steam pressure caused by the change of operation conditions. The conventional decoupling methods are also presented. Finally, detailed numerical simulations are conducted to validate the effectiveness of the IBICDC strategy.


Introduction
The 21st century is a century for human beings to explore the ocean. With the depletion of land resources, people pay an increasing attention to the ocean, with the ship technology playing an increasingly important role in social development. As the source of ship power, the related technology of marine power plant needs to be developed. As a large proportion of marine power plants, many technical problems of marine boiler-turbine system (MBTS) have not been effectively solved. One of the problems that attracts the attention of the academic and engineering circles is the problem of strong coupling between MB and MSTP. This problem has been around for a long time, the main difficulties are as follows: 1. Difficulty in modeling: The structure and operation mechanisms as well as the strong nonlinear and time-delay characteristics of MBTS are extremely complex. Moreover, there are few relevant authoritative literatures, and most of the authoritative literatures are about the modeling of steam power plant in the thermal power station in land [1][2][3][4][5][6][7][8][9][10][11][12][13].
The rest of the paper is organized as follows: Section 2.1 briefly introduces the MBTS and the modeling methods used in the paper; Section 2.2 presents the mathematical analysis and design of IBICDC. In Section 3, the simulation settings and results are summarized, and the effectiveness of the proposed method is verified. Then, in Section 4, conclusions and future work are summarized.

Modeling of MBTS
MBTS is an uncertain, time-varying, and non-linear multi-input and multi-output (MIMO) system [1]. The whole system is highly complex [2] and it can be divided into MB and MSTP parts. MB is responsible for the production of superheated and saturated steam (also called the main steam), the structure of which is shown in Figure 1. MB has ten parts. Furnace is the combustion unit, which burns fossil fuels (heavy petroleum) and releases a lot of heat [3]. The heat converts water in the drum into steam. The mixture of water and steam circulates between the steam drum and the mud drum through the Raiser and Downcomer. The steam in the superheater becomes superheated and saturated which is then used to drive the MSTP. In the MSTP, the main steam expands to work and drive the MSTP rotating or changing its rotational speed. The two components of MBTS are connected closely, and the devices of the two components exhibit different mechanical and control characteristics. The working conditions of MBTS are of high complexity, which need to be switched from one working condition to another frequently. Hence the dynamic characteristics of two components are not synchronous. So, when MBTS switches from one working condition to another, it can make the main steam pressure unstable, which in turn leads to a fluctuation in the rotational speed of the MSTP. This phenomenon is called coupling problem [1].  [3,4] In the furnace of MB, heat emits from the fuel combustion [10,11]. The heat exchanging process mainly include two ways [3]: a certain amount of heat of combustion exchanges in components such as Downcomer, Mud Drum, and Riser in the radiation manners, the left heat is exchanged through high temperature gas in the Superheater. At last, the main steam is generated in the Superheater. The combustion and heat exchanging processes can be expressed by the Equation (4) (1) where m denotes the main steam mass flow; u1 is the fuel mass flow; K1 is the gain of the processes; τ1 is the time delay of the process; T1 and T2 are the time constant of combustion and heat exchanging, respectively.

Combustion and Heat Exchanging
In the furnace of MB, heat emits from the fuel combustion [10,11]. The heat exchanging process mainly include two ways [3]: a certain amount of heat of combustion exchanges in components such as Downcomer, Mud Drum, and Riser in the radiation manners, the left heat is exchanged through high temperature gas in the Superheater. At last, the main steam is generated in the Superheater. The combustion and heat exchanging processes can be expressed by the Equation (4): m(s) u 1 (s) = K 1 (T 1 s + 1)(T 2 s + 1) e −τ 1 s (1) where m denotes the main steam mass flow; u 1 is the fuel mass flow; K 1 is the gain of the processes; τ 1 is the time delay of the process; T 1 and T 2 are the time constant of combustion and heat exchanging, respectively.

The Main Steam Pressure
The pressure of the main steam is a critical variable, and a main steam pressure with a good performance makes a boiler much more excellent than others [10,15]. Steam pressure is a function of temperature and density and the analytic formula of pressure state is very complex [33]. More than hundreds of gas state equations have been proposed [33], but most of them have no practical value and poor accuracy when applied in the practice.
Reference [33] presents a state equation of steam pressure, temperature, and volume as follows: where, R is the universal gas constant, 8.314472J/(K·mol); T denotes temperature, K; V is the molar volume of gases, m 3 /mol; P(V,T) denotes the steam pressure, Pa; b is the specific constant of steam volume, Pa(m 3 /mol) j ; A j , B j , and C j are the characteristic constants. Equation (2) is of high complexity and needs a large computation when applied in the practice. Hence, in this research, a much further simplified analytic formula is delivered [10].
The new equation is further simplified by the following assumption: the temperature of the main steam maintains stably at 529 • C.
where c is the characteristic constant; f s is the steam intake mass flow of MSTP, which is a function of the temperature, pressure of the main steam and the steam turbine governing valve opening u 2 [10,13]: where k v denotes the steam turbine governing valve gain, T sup denotes the temperature of the main steam-the control signal of MSTP. Note: for the detailed derivation of Equation (3a) see the Appendix A. According to the Sections 2.1.1 and 2.1.2, the following steam pressure model can be obtained:

The Model of MSTP
According to the kinetic energy theorem, rotational speed n(t) and the kinetic energy E k (t) of MSTP has the following relationship: Substituting Equation (4b) into Equation (4a), we have: Then, Equation (4b) can be written in this manner: where, J is the rotational inertia of MSTP; n(t) is rotational speed of MSTP; P S is the main steam work power in MSTP; P P is the propeller load power. When the rotational speed is close to 0, there is an infinite point of discontinuity in the rotational acceleration.
There is a homeomorphic mapping relationship between the kinetic energy and the rotational speed n(t). To simplify the model, take n 2 (t) as the controlled variable and we can obtained: The main steam work power has the relationship with the enthalpy value h 1 and h 2 of the steam in the inlet and outlet of MSTP, respectively: where the enthalpy value h 1 = h 1 (P in ,T in ) and h 2 = h 2 (P out ,T out ) of the steam can be calculated according to the IAPWS-IF97 [34]. While, in this paper, the P in , T in , P out , and T out are associated with MSTP loads, and they are functions of the rotational speed of MSTP, which can be achieved by data fitting, namely, P in = f 1 (n), T in = f 2 (n), P out = f 3 (n), T out = f 4 (n). And f 1 (n), f 2 (n), f 3 (n), and f 4 (n) are the three order polynomial functions. f 1 = α 11 n 3 + α 12 n 2 + α 13 n + α 14 f 2 = α 21 n 3 + α 22 n 2 + α 23 n + α 24 f 3 = α 31 n 3 + α 32 n 2 + α 33 n + α 34 f 4 = α 41 n 3 + α 42 n 2 + α 43 n + α 44 The propeller load power can be calculated by [35]: The load torque T P of propeller can be expressed as [35,36]: where, K Q is the dimensionless torque coefficient, ρ denotes the sea water density, n is the rotational speed of the MSTP, D is the diameter of the propeller. The load power of the propeller can be obtained as follows: Note: K Q is a function of the rotational speed of MSTP, namely, K Q = f 5 (n), and f 5 (n) is a 20 order polynomial function. The factors of f 5 (n) is also obtained based on actual data fitting. Substituting Equations (3b), (5), and (9) into Equation (4d), the following MSTP rotational speed model can be obtained: where the function ∆h can be further simplified:

The MBTS Coupling Model
According to the above equations and Figure 2, the following coupling model of MBTS can be obtained: Energies 2019, 12, x FOR PEER REVIEW 6 of 24 where the function Δh can be further simplified:

The MBTS Coupling Model
According to the above equations and Figure 2, the following coupling model of MBTS can be obtained:   According to Equation (1), the analytic relation between the fuel input u 1 (t) and the main steam flow m(t) is as follows: Substituting Equations (13a) into Equation (12), the coupling model of MBTS can be written in Equation (13c): The literatures on the decoupling control of MBTS are also scarce, which can only be seen in some academic conferences. And only a few literatures on decoupling control of steam power plant in the thermal power station in land can be seen in authoritative journals [1,[4][5][6][7][8][10][11][12][37][38][39][40]. In [8,[37][38][39][40], linear model predictive control (MPC) is used for the coupling problem. Most of the models used in [8,[37][38][39][40] come from [2], or identified by neural network and fuzzy methods. Since MBTS in this paper is a non-linear strong coupling multi-operation system, MPC is no longer applicable, and the use of nonlinear model predictive control (NMPC) will produce new problems: (1) nonlinear models and multiple operating conditions in MBTS will make the optimization problem of NMPC a non-convex optimization problem, which will cause the sequential quadratic programming method no longer be applicable; (2) generally, the stability of NMPC usually achieved by making the terminal state of the system into a specified terminal constraint set (which can be equality constraints or inequality constraints), but the current NMPC terminal constraint set lacks an effective design method. In addition, it will increase the amount of calculation and calculation time for online optimization, which makes the optimal control input designing very difficult; (3) the input, output, and state constraints not only lead to a large amount of calculation and long calculation time for MPC to solve constraint optimization online, but also cannot obtain the analytic expression of the optimal control input, and the optimal control input often exceeds the allowable range of input constraints.
Literature [21] proposed a compensatory decoupling method based on the linear time-invariant system, but the MBTS system is a nonlinear system, and the transfer function matrix required by the decoupling method cannot be obtained. Therefore, this decoupling method is not applicable to MBTS. Literature [22] put forward a hierarchical fuzzy sliding-mode (HFSM) decoupling control method. However, because of its high dependence on the structural characteristics of the model, it cannot be applied in this research. In [23], nonlinear models were transformed into linear models by local linearization method (Taylor expansion), and the optimal decoupling method (minimize squared tracking error) was adopted to solve the coupling problem of the system, but there were multiple operating conditions in the MBTS system, which required local linearization of the model at multiple operating points and a switching mechanism. In literature [24] a robust decoupling precompensator was proposed for the transfer function matrix coupling system. Reference [25] presented a decoupled optimization technique using genetic algorithms (GAs), and a stable self-learning PID control scheme was proposed in [27], and the PID parameters were adapted by the Lyapunov method to minimize the squared tracking error. Similarly, reference [28] achieved decoupling by using a support vector machine.
Reference [10,20] proposed a decoupling control algorithm based on the output feedback linearization, but there is a dynamic process when the control input is applied to the MBTS (for details, please refer to Equations (12) and (13)). As a result, this decoupling algorithm cannot be directly applied in this paper. To solve the coupling problem between torque and flux of induction motor, an adaptive state feedback linearized speed control algorithm is designed in [26]. Similar to the literatures [10,20], it is also difficult to apply this control algorithm directly in MBTS. Literature [29][30][31] presented an intelligent bi-cooperative decoupling controller inspired by the modulation mechanism of internal environment in human body, which is very useful for this study. However, like [10,20], this decoupling algorithm cannot be directly applied in MBTS. In order to solve this problem, the steam generation process is regarded as an inner control loop, and its input is fuel and output is steam flow. As shown in Equation (1), the model adopts a PID controller as follows.
The control system block diagram and simulation result of steam generation process is presented below.
Actually, the oxygen required for fuel combustion is supplied by air; in order to ensure sufficient burning, our group sets a definitive ratio β for oxygen entering the furnace, see Figure 3(a). The control input u 1 of the PID controller is the reference signal of the fuel flow control loop, βu 1 is the reference signal of the air flow control loop, and u f means the output of the fuel flow loop (the issues of the fuel flow loop and the air flow loop control are not discussed in the research, so u f = u 1 ). PID controller parameters: K P = 0.7496; K I = 0.01578; K D = 3.3057. The method was used to determine the PID controller parameters are the Ziegler-Nichols parameters tuning rules of PID method. em(t) is the error value between the main steam demand value and the main steam mass flow of MB. Then the coupling model of MBTS can be further expressed as: Note: with a controller for the steam generation process, the item "m(t)" in Equation (14a) is the decoupling control input of the MBTS, and recorded as u3 as below: 3 2 sup The method was used to determine the PID controller parameters are the Ziegler-Nichols parameters tuning rules of PID method. e m(t) is the error value between the main steam demand value and the main steam mass flow of MB. Then the coupling model of MBTS can be further expressed as: Note: with a controller for the steam generation process, the item "m(t)" in Equation (14a) is the decoupling control input of the MBTS, and recorded as u 3 as below: Using the Euler Discrete method, Equation (14c) can be written in discrete-time realization manner (the sample time is T s ): where u 3 (t) denotes the steam mass flow of MB; u 2 (t) is the opening of the steam turbine governing valve; P(t) is the main steam pressure value of MB; N is n 2 (t).

Neuroendocrine Regulation Principles (NERP) in the Human Body
In the human body, the biochemical systems, organs, and physiological metabolisms are not independent but interactive with each other [41]. Among these subsystems in the human body, the nervous system and endocrine system play much more important roles [30]. The nervous system and the endocrine system work cooperatively to establish a stable internal environment in the human body. Therefore, a stable internal environment is extremely critical to the normal physiological metabolisms [31,41].
The nervous system is the main regulating system in the human body. Under the direct or indirect manipulation of the nervous system, all the organs, systems, and all kinds of physiological metabolisms are not working alone. They are interconnected and act in close coordination, which makes the human body a whole and complete organism, and realizes and maintains the human life and behaviors.
The nervous system is divided into central nervous system and peripheral nervous system. The peripheral nervous system is distributed widely all over the body, and binds the brain and spinal cord with the other organs closely. This structure allows the central nervous system easy to acquire information of the inner and outer environments of the human body (afferent neuron transmit sensory information). The conventional control system contains a controller, actuator, and plant. In the Figure 4, the nervous system consists of many organs, glands, sensors corresponding to components of the conventional control system. Hence, there remain similarities between the nervous system and conventional control system, and the nervous system principles can be applied to MBTS.  According to the neuroendocrine regulation principles in the human body, the IBICDC control system structure is shown in Figure 5, in which the decoupling coordination control algorithm IBICDC is divided into two parts [15].
Part I-basic part: This part contains two PID controllers, PID1 is for the main steam pressure P, According to the neuroendocrine regulation principles in the human body, the IBICDC control system structure is shown in Figure 5, in which the decoupling coordination control algorithm IBICDC is divided into two parts [15].
Part I-basic part: This part contains two PID controllers, PID1 is for the main steam pressure P, and PID2 is for the N of MSTP; Part II-cooperative decoupling controller: The function of this part is mainly to decoupling the coupling influence between the main steam pressure control loop and the rotational speed control loop, eliminating the violent oscillation between the two loops.

The Structure and Scheme of IBICDC
According to the neuroendocrine regulation principles in the human body, the IBICDC control system structure is shown in Figure 5, in which the decoupling coordination control algorithm IBICDC is divided into two parts [15].
Part I-basic part: This part contains two PID controllers, PID1 is for the main steam pressure P, and PID2 is for the N of MSTP; Part II-cooperative decoupling controller: The function of this part is mainly to decoupling the coupling influence between the main steam pressure control loop and the rotational speed control loop, eliminating the violent oscillation between the two loops.  According to the Equation (14), it can be obtained that the main steam pressure consists of two parts: 1. the main steam generated by MB increases the main steam pressure value; 2. the main steam consumed by MSTP reduces the main steam pressure value.
e P,k = P re f er,k − P k e N,k = N re f er,k − N k where: U P,k and U N,k are the control signals of PID1 and PID2, respectively.
PID1(e P,k , k) = K P1 (e P,k − e P,k−1 ) + K I1 e P,k + K D1 (e P,k − 2e P,k−1 + e P,k−2 ) PID2(e N,k , k) = K P2 (e N,k − e N,k−1 ) + K I2 e N,k + K D2 (e N,k − 2e N,k−1 + e N,k−2 ) The design of the IBICDC decoupling coordination control algorithm: According to the Equations (14) and (15): where u 3 denotes the main steam of MB; u 2 is the opening of the steam valve of MSTP; u 11 is U P,k−1 ; u 12 is the decoupling coordination to MSTP; u 21 is the decoupling coordination to MB; u 22 is U N,k−1 .
The following focuses on the calculation of u 3 and u 2 ; the main steam pressure P should remain unchanged, so dP/dt = 0 in Equation (19): . P = f 11 u 3 + f 12 u 2 (21) .
As dP/dt = 0, which means P is a constant value, transforming Equation (22) into a time-discrete form: f 11,k (∆u 11,k + ∆u 21,k ) + f 12,k (∆u 12,k + ∆u 22,k ) = 0 When MBTS is running in the steady state operation, the parameters of the whole system can be considered as constant value, which implies P/dt = 0, dN/dt = 0, ∆u 11 = 0, ∆u 21 = 0, ∆u 12 = 0, and ∆u 22 = 0. When the operating condition of the system is changing, the rotational speed of MSTP needs to be changed. The main steam pressure of MB changes with the rotational speed. When the main steam pressure value changes sharply and exhibits large fluctuations, the decoupling coordination controller begins to work, which means dP/dt = 0, dN/dt 0 (∆N 0), ∆u 11 = 0, ∆u 21 0, ∆u 12 = 0, and ∆u 22 0. According to Equation (23), we can have: In the same way, when the main steam pressure of MB changes: dP/dt 0, dN/dt = 0(∆N = 0), we can calculate ∆u 12 : According to the principles of the human internal environment, information on changes of environmental parameters inside and outside the human body is collected by the nerve network. However, the reference trajectory (equivalent to the external parameters of human body) and the variation information of the rotational speed (equivalent to the internal parameters of human body) in the MBTS are not reflected in Equation (15). Therefore, in order to achieve a better decoupling system, drastic changes in the main steam pressure should be reduced, the change information of reference trajectory and the MSTP rotational speed to the decoupling controller need to be added. The decoupling controller is represented by: In the same way, when the main steam pressure of MB changes: dP/dt 0, dN/dt = 0 (∆N = 0), we can calculate ∆u 12 : Therefore, u 21 can offset the coupling influence of the MSTP rotational speed loop on the main steam pressure loop. Moreover, u 12 can offset the coupling influence of the main steam pressure loop on the rotational speed control loop, in turn.

The Flexible Reference Trajectory (FRT) Scheme
A critical mechanism of the human neuroendocrine regulation system is slow regulation; when the external interference alters the human body internal environment parameters. The buffering effect of the internal environment in the human body prevents parameters from changing dramatically in a moment as drastically as a step or pulse signal, but smoothly and slowly deviates from the normal value. Then the human body regulation mechanism can adjust the changing parameters gradually to the normal value. In this process, parameters slowly and steadily move from the deviation value to the normal value in a smooth manner, instead of in a drastic manner (which can result in overshoot and waste body matter and energy). Hence, in order to embody this mechanism in the IBICDC, a FRT is introduced in this research. FRT replaces the step reference signal in the control input, which ensures the output can slowly and smoothly follow FRT during the state switching process. The details are described as follows: According to Equation (14c), initial conditions of the main steam pressure P and rotational speed N are P(0) and N(0), respectively.
In order to achieve a simpler formulary of Equation (29), Then, Equation (29) can be transformed to: In order to prove the coupling influence of the MSTP rotational speed control loop on the main steam pressure control loop can be reduced significantly under FRT. Taking the 2-norm of both sides of Equation (31), we can get: Let f 13,0 = 0, f 23,0 = 0. Substituting (16)- (18) into (30): where Let U 1 = KPID 1 P re f er,0 − P 0 P re f er,1 − P 1 · · · P re f er,k−1 − P k−1 T Similarly, it can be obtained: where Substituting (35) and (36) into (30): P srefer,k = aP k−1 + bP refer,k N srefer,k = aN k−1 + bN refer,k (39a) e sP,k = P sre f er,k − P k e sN,k = P sre f er,k − P k (39b) where P srefer,k , N srefer,k are the FRT of the main steam pressure P and rotational speed N, respectively; P refer and N refer are the step reference signals of the system. a and b are the flexible factors, 0 < a < 1, 0 < b < 1, and a + b = 1. Take the 2-norm to Equations (37), (38), and (39a) respectively: Substituting (37) and (40) into (32): It can be known that the control law obtained by introducing the FRT has a better control performance.

Simulations Results and Discussion
The load of MBTS can be divided into the following two types: steady state conditions and dynamic switching conditions. Steady state conditions: (1) high navigational speed condition; (2) medium navigational speed condition; (3) low navigational speed condition. Dynamic switching conditions: (1) high navigational speed directly reduces to low navigational speed conditions; (2) low navigational speed directly rises to high navigational speed conditions; (3) high navigational speed step-down condition; (4) low navigational speed step-rise condition.
The simulation in this paper is mainly for the dynamic switching conditions of MBTS during the navigational speed switching. The rotational speed of the MSTP at different working conditions in the simulations is 180 r/min in the high navigational speed condition, 120 r/min at medium navigational speed condition, and 60 r/min in the low navigational speed conditions. Therefore, four simulation conditions are set as follows: low navigational speed step-rising condition with MSTP rotational speed of 60 r/min→120 r/min→180 r/min; low navigational speed directly rising to high speed condition with MSTP rotational speed of 60 r/min→180 r/min; high navigational speed step-down condition with MSTP rotational speed 180 r/min→120 r/min→60 r/min; high navigational speed directly reducing to low speed condition with MSTP rotational speed 180 r/min to 60 r/min. In order to verify IBICDC, other methods such as: PID, PIDDC, and Intelligent PIDDC have also been performed. For more details see the Appendix B. All the parameters listed at Table 1. From Figures 6-9, the following conclusions can be drawn: 1.The introduction of the FRT can greatly reduce the amplitude oscillation of the main steam pressure, which indicates that the FRT can buffer the coupling influence of the MSTP on the main steam pressure when the operating condition changes; 2.The PID control method can make the main steam pressure reach to the steady-state value when the operating condition changes, but the dynamic process control performance is very poor, and the main steam pressure changes sharply; 3.Compared with the PID control, the PIDDC method can eliminate the coupling influence on the main steam pressure to some extent, but the control performance is not ideal; 4.Compared with the PIDDC method, Intelligent PIDDC further eliminates the coupling influence of the main steam pressure, but the amplitude of the main steam pressure fluctuation is still relatively large; 5.The IBICDC introduces more information about the internal and external changes of the system, so that it can apply more effective control input, eliminate the coupling influence of the main steam pressure, and has a better dynamic control performance than the previous three methods.
Moreover, the simulation results show that although the PIDDC and Intelligent PIDDC control algorithms can also achieve steady-state response, the variation range of all the parameters are beyond the normal range obviously. So, the IBICDC coordinated control algorithm has a better control performance. beyond the normal range obviously. So, the IBICDC coordinated control algorithm has a better control performance.

Conclusions
In this paper, an IBICDC control algorithm is presented based on the NERP of the human body, and applied to the coupling control problems of MBTS. The simulation results show that IBICDC can respond to the expected output quickly, reduce the coupling influence on the main steam pressure, and improve the MBTS control performance. IBICDC makes a smooth regulation with a smaller overshoot, which means less fuels consumption. Compared with other decoupling control methods, IBICDC has the advantages of simple structure, less computation, and easy implementation, so it is worth being applied in engineering. IBICDC can be applied in complicated multivariable process control systems, but many problems occur when applied to the motion control systems (such as ship track control, unmanned aerial vehicle, and unmanned submersibles, etc.). In this research, the FRT in IBICDC can be applied in other control systems with step reference as tracking trajectory, which is conducive to effectively reducing the overshoot and improving the control performance. In addition, in many cases, the control input of the system will exert constraint, and this will result in a great trouble to the controller design. The FRT method can be adopted to avoid such problem. The future work includes: (1) there is an astern working condition for the MSTP. However, this study does not verify the switching process of astern and forward working conditions. So, the next step is to solve the problem. (2) other control approaches based on optimization such as nonlinear model predictive control for this problem are with high research value, and our group is working on this approach.

Conclusions
In this paper, an IBICDC control algorithm is presented based on the NERP of the human body, and applied to the coupling control problems of MBTS. The simulation results show that IBICDC can respond to the expected output quickly, reduce the coupling influence on the main steam pressure, and improve the MBTS control performance. IBICDC makes a smooth regulation with a smaller overshoot, which means less fuels consumption. Compared with other decoupling control methods, IBICDC has the advantages of simple structure, less computation, and easy implementation, so it is worth being applied in engineering. IBICDC can be applied in complicated multivariable process control systems, but many problems occur when applied to the motion control systems (such as ship track control, unmanned aerial vehicle, and unmanned submersibles, etc.). In this research, the FRT in IBICDC can be applied in other control systems with step reference as tracking trajectory, which is conducive to effectively reducing the overshoot and improving the control performance. In addition, in many cases, the control input of the system will exert constraint, and this will result in a great trouble to the controller design. The FRT method can be adopted to avoid such problem. The future work includes: (1) there is an astern working condition for the MSTP. However, this study does not verify the switching process of astern and forward working conditions. So, the next step is to solve the problem. (2) other control approaches based on optimization such as nonlinear model predictive control for this problem are with high research value, and our group is working on this approach.

Conclusions
In this paper, an IBICDC control algorithm is presented based on the NERP of the human body, and applied to the coupling control problems of MBTS. The simulation results show that IBICDC can respond to the expected output quickly, reduce the coupling influence on the main steam pressure, and improve the MBTS control performance. IBICDC makes a smooth regulation with a smaller overshoot, which means less fuels consumption. Compared with other decoupling control methods, IBICDC has the advantages of simple structure, less computation, and easy implementation, so it is worth being applied in engineering. IBICDC can be applied in complicated multivariable process control systems, but many problems occur when applied to the motion control systems (such as ship track control, unmanned aerial vehicle, and unmanned submersibles, etc.). In this research, the FRT in IBICDC can be applied in other control systems with step reference as tracking trajectory, which is conducive to effectively reducing the overshoot and improving the control performance. In addition, in many cases, the control input of the system will exert constraint, and this will result in a great trouble to the controller design. The FRT method can be adopted to avoid such problem. The future work includes: (1) there is an astern working condition for the MSTP. However, this study does not verify the switching process of astern and forward working conditions. So, the next step is to solve the problem. (2) other control approaches based on optimization such as nonlinear model predictive control for this problem are with high research value, and our group is working on this approach.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
In accordance with achievements by Martin-Hou: In 2016, REN Tieliang, WAN Xue et al. simplified (A1) into: In our study, the main steam pressure range is: [7MPa, 10MPa], the desire value of the main steam pressure is 8 MPa, the molar specific volume V = 1.0683047401 × 10 −3 , b = 6.6454139232 × 10 −6 . The molar specific volume V of steam is much larger than b value, so b in item (V-b) of the Equation is negligible and can be directly regarded as V, obtaining the simplified results in this paper as follows: The derivative of the main steam pressure P to time t can be expressed as: The steam temperature T = 529 • C is a fixed value, so there is dT(t)/dt = 0. Equation (A4) can be further simplified as: dP(V(t), T(t)) dt = ∂P(V(t), T(t)) ∂V(t) According to Equation (A2), we have: In the MB, there is V(t) = M e /ρ(t) that can be substituted into Equation (A6) to obtain: where, M(t) is the steam mass in the steam header connecting MB and MSTP, V t is the volume of steam in the main steam header, and M e is the molar mass of water molecules, ρ is the main steam density. The derivative of the main steam mass in the header is: where, m(t) is the amount of steam entering to the header from MB, and f s (t) is the amount of steam consumed by the MSTP. It shows that: Substituted Equations (A7) and (A11) into Equation (A5): Then let: Then Equation (A12) can be expressed as: As for the ship, within the value range of the steam pressure, the value of c changes little, so it is regarded as a constant value, acquiring the equation below: where, u3,k = PID1(k), u2,k = PID2(k).  where u11,k = PID1(k), u22,k = PID2(k), u3(k) = u11,k + u21,k and u2(k) = u12,k + u22,k. 12 22, 21 11 where, u 3,k = PID1(k), u 2,k = PID2(k).
The proposed IBICDC method, with block diagram as follows:  Figure A3. Intelligent PIDDC block diagram of MBTS.