# New Smith Internal Model Control of Two-Motor Drive System Based on Neural Network Generalized Inverse

Authors: Guohai Liu; Jun Yuan; Wenxiang Zhao; Yaojie Mi
Journal: Journal of Control Science and Engineering
 (2016)

---

## Abstract

Multimotor drive system is widely applied in industrial control system. Considering the characteristics of multi-input multioutput, nonlinear, strong-coupling, and time-varying delay in two-motor drive systems, this paper proposes a new Smith internal model (SIM) control method, which is based on neural network generalized inverse (NNGI). This control strategy adopts the NNGI system to settle the decoupling issue and utilizes the SIM control structure to solve the delay problem. The NNGI method can decouple the original system into several composite pseudolinear subsystems and also complete the pole-zero allocation of subsystems. Furthermore, based on the precise model of pseudolinear system, the proposed SIM control structure is used to compensate the network delay and enhance the interference resisting the ability of the whole system. Both simulation and experimental results are given, verifying that the proposed control strategy can effectively solve the decoupling problem and exhibits the strong robustness to load impact disturbance at various operations.


---

## Body

## 1. Introduction

Multimotor drive system is a multi-input multioutput (MIMO), nonlinear, strongly coupled, and time-varying delay system, which has been widely applied in modern industry. High precision coordinated control performance can improve quality and productivity of products, such as electric vehicles and rail transits [1–4]. In order to meet the requirements of industrial applications, the delay problem of the system must be solved, and then the coupling variables of this system can be decoupled.In a multimotor drive system, decoupling speed and tension is a key issue. Traditional decoupling control methods strongly depend on mathematical model of the system, such as PID control, sliding mode control [5], feed-forward control [6], and adaptive fuzzy control [7]. Moreover, the structure characteristics or parameters are easily influenced by load variation and interference factor. So, it is hard to obtain accurate system mathematical model. A new control strategy based on support vector machine (SVM) theory has been proposed to decouple current and rotational speed of a permanent-magnet (PM) motor system [8]. Recently, neural network generalized inverse (NNGI) is used to decouple strong-coupling controlled variables [9–11]. Torque and flux components of five-phase PM motors were decoupled by multiple-reference-frame transformation, in which the artificial neural network (ANN) controller was trained online to adapt system uncertainties [12]. Based on a single artificial neuron requiring no offline training, an intelligent speed controller for the PM motor was proposed to adapt to various drives without extensive knowledge of motor behavior [13]. Since the analysis model of the traveling-wave ultrasonic motor is difficult to obtain, a generalized regression neural-network-based model is developed to solve the problem and the transfer function identification is no longer required [14]. An ANN-based estimator is implemented to eliminate mechanical sensors and then realize the sensorless control of a PM motor [15]. By using a pseudolinear composite system, NNGI can transform a MIMO nonlinear system into several single-input single-output linear systems, in which an accurate mathematic model is not required.In addition, the demand of real-time transmission is high because a large transmission delay will result in uncontrollable system in practice. So, an effective control strategy is necessary to implement the delay compensation. Traditional Smith predictor control is an effective method to conquer the influence of system delay [16, 17]. However, this control method is sensitive to the errors of estimate model and external disturbance. An inaccurate estimate model or a strong interference will affect control quality and even lead to system instability. According to the strong robustness and anti-interference abilities of internal model control [18], the control structure is widely used in many applications. A strategy based on internal model control was proposed for a matrix converter-based PMSM drive system to reduce the adverse impact on drive performance caused by nonlinear output characteristics of matrix converter in the case of input voltage disturbance [19]. The speed regulation problem for a PM motor drive was settled by a fuzzy adaptive law based internal model control scheme, and the effectiveness of the proposed methods has been verified [20]. An internal model control with a conditional integrator was proposed for the robust output regulation of a DC/DC buck converter [21]. In the work, Smith predictor control is codesigned with internal model control structure, and then the Smith internal model (SIM) control structure is formed. This SIM control incorporates the advantages of both Smith estimate control and internal model control, which can not only reduce the effect of network delay, but also improve the robustness of system. However, since it is hard to obtain the estimate model in multimotor system, the application of SIM control is very limited. On the basis of NNGI method, the accurate model of generalized pseudolinear subsystem will be obtained and the estimate model can be acquired without mismatch. Therefore, the delay compensation can be implemented successfully.The purpose of this paper is to propose a new control method, namely, NNGI-based Smith internal model (NNGI-SIM) control, to improve the operating performance of a two-motor drive system. This NNGI-SIM control can avoid the mismatch between prediction model and original model. The theory of NNGI decoupling method is introduced first; then, the two-motor drive system is decoupled by NNGI method; the accurate model of generalized pseudolinear subsystem will be obtained and the mismatch condition will be conquered. After the theory of SIM control is deduced, the NNGI system is structured and codesigned with SIM control strategy to form the new NNGI-SIM control and then applied in two-motor drive system.The structure of this paper is as follows: in Section1, the multimotor drive system will be introduced. In Section 2, the mathematical model of two-motor drive system will be deduced. In Section 3, the new NNGI-SIM control will be formed and applied in two-motor drive system. In Sections 4 and 5, decoupling effect and robustness under load disturbance of the proposed two-motor drive will be verified by both simulations and experiments. Finally, conclusions will be drawn in Section 6.


## 2. Mathematical Model

A two-motor synchronous system is shown in Figure1; it consists of two three-phase asynchronous induction motors, in which motor 1 is the master motor and motor 2 is the slave one. The belt-pulley is installed on the motor shaft, and both motors are combined by transmission belt on the pulley one. Floating rollers strain the belt to increase the friction between the belt and the driving rollers. Two motors operate synchronously by adjusting the action between speed and tension.Figure 1
Physical model diagram of two-motor drive system.When both motors operate in the vector control mode, the rotor flux is considered to be invariant. Then, the mathematic model of the system can be written as(1)
ω
˙
r
1
=
n
p
1
J
1
ω
1
-
ω
r
1
n
p
1
T
r
1
L
r
1
ψ
r
1
2
-
T
L
1
+
r
1
F
,
ω
˙
r
2
=
n
p
2
J
2
ω
2
-
ω
r
2
n
p
2
T
r
2
L
r
2
ψ
r
2
2
-
T
L
2
-
r
2
F
,
F
˙
=
K
T
1
n
p
1
r
1
k
1
ω
r
1
-
1
n
p
2
r
2
k
2
ω
r
2
-
F
T
,where ω
i and ω
r
i are the synchronous angular speed and the electric angular speed, F is the tension of the belt, n
p
i is the pole-pairs number of the number i motor, J
i, ψ
r
i, and L
r
i are rotor inertia, rotor flux, and rotor self-inductance, T
r
i is the electromagnetic time constant, T
L
i is the load torque, K
=
E
/
v is the transfer coefficient, T
=
L
0
/
A
v is the time constant of tension variation, E is Young’s modulus of elasticity, v is the expected line speed, L
0 is the distance between racks, A is the section area, r
i is the radius of belt-pulley, and k
i is the speed ratio (i
=
1
,
2), respectively.


## 3. Controller Design

### 3.1. NNGI System

Generalized inverse system can realize the inverse mapping relationship from the output to the input of original system. By setting the original order, derivatives and their linear combination of expected output of original system can be used as the input of the generalized inverse system. Then, the output can be used as the input of original system to drive the original system, generating the expected output [22].As a MIMO nonlinear system which hasq-dimensional input vector u
=
u
1
,
u
2
,
…
,
u
q
T
∈
R
q and q-dimensional output vector y
=
y
1
,
y
2
,
…
,
y
q
T
∈
R
q, the differential equation can be expressed as(2)
F
y
ε
T
Y
,
u
σ
T
,
U
=
0
,where(3)
ε
=
ε
1
,
…
,
ε
q
T
,
σ
=
σ
1
,
…
,
σ
q
T
,
y
ε
=
y
1
ε
1
,
…
,
y
q
ε
q
,
u
σ
=
u
1
σ
1
,
…
,
u
q
σ
q
,
Y
=
y
1
ε
1
-
1
,
…
,
y
1
,
…
,
y
q
ε
q
-
1
,
…
,
y
q
T
,
U
=
u
1
σ
1
-
1
,
…
,
u
1
,
…
,
u
q
σ
q
-
1
,
…
,
u
q
T
.
ε
j and σ
j are the highest derivative orders of input u
j and output y
j.If the relative order of the systemα
=
α
1
,
α
2
,
…
,
α
q
T is existent and the sum of relative order ∑
j
=
1
q
α
j
=
n (n is the order of the system), the state equation of the system can be transformed into differential equation. Hence, when the relative order is equal to the original order, the generalized inverse of the system can be directly obtained as(4)
u
=
ϕ
¯
Y
¯
,
v
¯
,where(5)
Y
¯
=
y
1
,
y
1
1
,
…
,
y
1
α
1
-
1
,
…
,
y
q
,
y
q
1
,
…
,
y
q
α
q
-
1
T
,
ν
¯
=
v
1
¯
,
v
2
¯
,
…
,
v
q
¯
T
,
v
j
¯
=
a
j
0
y
j
+
a
j
1
y
j
1
+
⋯
+
a
j
α
j
-
1
y
j
α
j
-
1
+
a
j
α
j
y
j
α
j
.
j
=
1
~
q
.By connecting the generalized inverse system in series, the coupled original system can be decoupled into several generalized pseudolinear composite subsystems and the transfer function of the subsystem can be expressed as(6)
G
s
=
d
i
a
g
G
11
s
,
…
,
G
q
q
s
=
d
i
a
g
1
a
1
α
1
s
α
1
+
a
1
α
1
-
1
s
α
1
-
1
+
⋯
+
a
11
s
+
a
10
,
…
,
1
a
q
α
q
s
α
q
+
a
q
α
q
-
1
s
α
q
-
1
+
⋯
+
a
q
1
s
+
a
q
0
.In practice, the structure of MIMO nonlinear system is complicated and it is hard to obtain its accurate mathematical model. Even if the existence of generalized inverse system has been proved, the expression of generalized inverse of the system cannot be derived for various occasions. To solve this problem, an ANN is used to approximate the original system. This method is independent of the accurate mathematical model, thus enhancing the popularization of inverse system algorithm. The NNGI pseudolinear composite subsystem equivalent chart is shown in Figure2.Figure 2
NNGI pseudolinear composite subsystem equivalent chart.


### 3.2. NNGI Construction for Two-Motor Drive System

According to the mathematic model of two-motor drive system in (1), the existence of generalized inverse system will be proved. Variables are chosen as follows.State variable is(7)
x
=
x
1
,
x
2
,
x
3
T
=
ω
r
1
,
ω
r
2
,
F
T
;control variable is(8)
u
=
u
1
,
u
2
T
=
ω
1
,
ω
2
T
;output variable is(9)
y
=
y
1
,
y
2
T
=
ω
r
1
,
F
T
.Then, the mathematic model of two-motor drive system can be expressed as(10)
x
˙
=
f
x
,
u
=
n
p
1
J
1
u
1
-
x
1
n
p
1
T
r
1
L
r
1
ψ
r
1
2
-
T
L
1
+
r
1
x
3
n
p
2
J
2
u
2
-
x
2
n
p
2
T
r
2
L
r
2
ψ
r
2
2
-
T
L
2
-
r
2
x
3
K
T
1
n
p
1
r
1
k
1
x
1
-
1
n
p
2
r
2
k
2
x
2
-
x
3
T
.According to the generalized inverse theory, the corresponding Jacobi matrix can be described as(11)
J
x
,
u
=
∂
y
1
1
∂
u
1
∂
y
1
1
∂
u
2
∂
y
2
2
∂
u
1
∂
y
2
2
∂
u
2
=
n
p
1
2
T
r
1
ψ
r
1
2
J
1
L
r
1
0
K
r
1
k
1
n
p
1
T
r
1
ψ
r
1
2
T
J
1
L
r
1
-
K
r
2
k
2
n
p
2
T
r
2
ψ
r
2
2
T
J
2
L
r
2
,
Det
J
x
,
u
=
K
r
2
k
2
n
p
1
2
n
p
2
T
r
1
T
r
2
ψ
r
1
2
ψ
r
2
2
T
J
1
J
2
L
r
1
L
r
2
.Whenψ
r
1
≠
0 and ψ
r
2
≠
0, the Jacobi matrix is nonsingular. The relative order of the system is α
=
(
α
1
,
α
2
)
=
(
1,2
), α
1
+
α
2
=
3, and it is equal to the order of the system. In addition, the original order of the system is n
e
=
(
n
e
1
,
n
e
2
)
=
(
1,2
), n
e
1
+
n
e
2
=
3, and it is the same as relative order. So, the inverse of the system is existent and can be expressed as(12)
u
=
ϕ
-
y
1
,
y
˙
1
,
y
2
,
y
˙
2
,
y
¨
2
,
v
-
,where(13)
v
¯
=
v
1
¯
,
v
2
¯
T
,
v
1
¯
=
a
10
y
1
+
a
11
y
1
1
,
v
2
¯
=
a
20
y
2
+
a
21
y
2
1
+
a
22
y
2
2
.In (12), to keep the stability of the open-loop system, a
10
=
1, a
11
=
1; a
20
=
1, a
21
=
1.414, and a
22
=
1. Then, the transfer function of the subsystem is(14)
G
s
=
d
i
a
g
G
11
,
G
22
=
d
i
a
g
1
s
+
1
,
1
s
2
+
1.414
s
+
1
.By connecting the generalized inverse system in series, the original system can be decoupled and simplified into a first-order subsystem and a second-order subsystem. As a result, the composite pseudolinear system is derived and shown in Figure3.Figure 3
Pseudolinear composite system.Simultaneously, by adjusting the parametersa
10, a
11, a
20, a
21, and a
22, the zeros and poles of pseudolinear subsystems can be developed exactly. Two generalized pseudolinear subsystems have neither right half-plane zero, and then the original system will obtain open-loop stabilization after it is in series with the NNGI system. The NNGI open-loop control diagram is shown in Figure 4.Figure 4
NNGI open-loop control diagram.


### 3.3. Profibus-DP Network Delay

Profibus communication uses the main token polling mechanism; each site on the bus connects with another in the bus way. Master station is composed of control center such as computer and PLC, while slave station consists of actuators such as inverter and sensors. Each site on the Profibus has its own logical address; a loop called logical loop is formed between the master stations. A token generated during the operation of the frame will poll along the logical link between the master stations. One master station is allowed to poll along other slave stations only when a token reaches this master station, and then the data transmission will be finished. While the process from sending task information of site to receiving the task information by target is too time intensive, the Profibus-DP network delay is generated.Profibus-DP network delay is indeterminate in practice. However the maximum delay can be calculated and it can be a reference for the delay of ideal internal modele
-
m
s [23, 24]. The approximate expression can be written as(15)
R
h
=
B
+
i
n
t
n
h
n
+
1
·
T
T
R
+
T
c
h
m
+
τ
+
ξ
h
,where B stands for the longest initial blocking, n
h is the number of high priority request messages added to the outgoing queue at critical instant, T
T
R is the target token rotation time, T
c
h
m is the maximum time of executing a high priority message transmission, τ is the token latency, n
+
1 is the number of high priority messages which are transmitted during the maximum token holding time, and ξ
h is the time between the finish of the last maximum token holding time and the end of target message transmission. There are three possible cases and the expression can be written as(16)
ξ
h
=
-
τ
,
n
h
=
i
n
t
n
h
n
+
1
·
n
+
1
,
T
c
h
m
,
n
h
=
i
n
t
n
h
n
+
1
·
n
+
1
+
1
,
n
h
-
i
n
t
n
h
n
+
1
·
n
+
1
·
T
c
h
m
+
τ
,
n
h
>
i
n
t
n
h
n
+
1
·
n
+
1
+
1
.


### 3.4. SIM Control

According to the previous discussions, two-motor drive system is a time-varying delay system. In order to settle the time delay issue, a SIM control is introduced into the controller design. The SIM control incorporates the advantages of Smith estimate control and internal model control. This control strategy is a delay compensation method, which is to separate the delay part of controlled object and connect in series with the system transfer function. It is an effective way to inhibit the effect of time delay. Moreover, SIM control can also restrain the influence of unmeasured disturbances and improve system robustness.The structure of the SIM control is shown in Figure5, in which G
P
(
s
) and G
m
(
s
) are the actual object model and the prediction internal model, respectively. e
-
P
s and e
-
m
s are the related delay of controlled object and prediction internal model, respectively. G
c
(
s
) is the controller. D
(
s
) is the external disturbance. R
(
s
) and Y
(
s
) are the input and the output of the system, respectively.Figure 5
SIM control diagram.


#### 3.4.1. Delay Compensation Analysis

According to Figure5, the output can be written as(17)
Y
s
=
G
c
s
G
P
s
e
-
P
s
1
+
G
c
s
G
m
s
+
G
c
s
G
P
s
e
-
P
s
-
G
m
s
e
-
m
s
R
s
+
1
+
G
c
s
G
m
s
-
G
c
s
G
m
s
e
-
m
s
1
+
G
c
s
G
m
s
+
G
c
s
G
P
s
e
-
P
s
-
G
m
s
e
-
m
s
D
s
.The closed-loop transfer function of control system is(18)
G
s
=
Y
s
R
s
=
G
c
s
G
P
s
e
-
P
s
1
+
G
c
s
G
m
s
+
G
c
s
G
P
s
e
-
P
s
-
G
m
s
e
-
m
s
.When the prediction internal model is accurate,G
m
s
e
-
m
s is equal to G
P
s
e
-
P
s, and the transfer function can be written as(19)
G
s
=
G
c
s
G
P
s
1
+
G
c
s
G
P
s
e
-
P
s
.Hence, there is no delay part in the denominator of closed-loop transfer function and the delay is totally compensated.


#### 3.4.2. Fixed Value Tracking Ability

According to Figure5, when the external disturbance D
(
s
) is zero, then the error of closed-loop system can be written as(20)
E
s
=
R
s
-
Y
s
=
1
-
G
n
c
s
G
m
s
e
-
m
s
1
+
G
n
c
s
G
P
s
e
-
P
s
-
G
m
s
e
-
m
s
R
s
.In consideration of modeling error of actual system, the controlled object can be described with multiplicative uncertainty:(21)
G
P
s
e
-
P
s
=
G
m
s
e
-
m
s
1
+
l
m
s
,where l
m
(
s
) represents the uncertainty and it is bounded.The structure in dashed box is represented byG
n
c
(
s
). A low pass filter is added to improve the capacity of resisting disturbance. The filter is represented as F
(
s
), and then G
n
c
(
s
) can be written as(22)
G
n
c
s
=
G
m
-
1
s
F
s
.Therefore, the transfer function of low pass filter is(23)
F
s
=
1
λ
s
+
1
n
,where λ describes filter parameter and n stands for the order of filter.Then, the error of closed-loop system can be simplified as(24)
E
s
=
λ
s
+
1
n
-
e
-
m
s
λ
s
+
1
n
+
l
m
s
e
-
m
s
R
s
.The limit of l
m
(
s
) can be expressed as(25)
lim
s
→
0
⁡
l
m
s
=
C
,where C is a constant.When the input signal is step signal,(26)
R
s
=
R
0
s
,then, the limit of E
(
s
) is(27)
e
∞
=
lim
s
→
0
⁡
s
E
s
=
lim
s
→
0
⁡
s
λ
s
+
1
n
-
e
-
m
s
λ
s
+
1
n
+
l
m
s
e
-
m
s
·
R
0
s
=
0
.


#### 3.4.3. Anti-Interference Ability

According to Figure5, when the input R
(
s
) is zero, then the error of closed-loop system can be written as(28)
E
s
=
R
s
-
Y
s
=
1
-
G
n
c
s
G
m
s
e
-
m
s
1
+
G
n
c
s
G
P
s
e
-
P
s
-
G
m
s
e
-
m
s
-
D
s
.Similar to the analysis of fixed value tracking ability, the error of closed-loop system can be simplified as(29)
E
s
=
λ
s
+
1
n
-
e
-
m
s
λ
s
+
1
n
+
l
m
s
e
-
m
s
-
D
s
.When the external disturbance is step signal,(30)
D
s
=
D
0
s
,then, the limit of E
(
s
) is(31)
e
∞
=
lim
s
→
0
⁡
s
E
s
=
lim
s
→
0
⁡
s
-
λ
s
+
1
n
-
e
-
m
s
λ
s
+
1
n
+
l
m
s
e
-
m
s
·
D
0
s
=
0
.In conclusion, the SIM control structure can satisfactorily solve delay problem. What is more, this control structure also has good input signal tracking ability and anti-interference ability.


### 3.5. NNGI-SIM Control

As mentioned above, SIM control can remove the delay of forward network from the closed-loop and set it as a gain block before the output. As a result, the delay on the return network path can be totally eliminated. Meanwhile, SIM control can also enhance robustness and anti-interference ability of the system with the internal model control structure. However, the accuracy of prediction modelG
m
(
s
) is crucial. When there is a mismatch between prediction model and actual model, the control quality will be significantly deteriorated. This problem can be solved by NNGI system. After connecting with NNGI system, the original system is decoupled into two subsystems, namely, a first-order speed one and a second-order tension one. The accurate mathematical model of two subsystems can be obtained after pole-zero allocation of the pseudolinear subsystems, and then the internal model can be predicted. Then, the mismatch condition will be avoided. The diagram of NNGI-SIM control strategy of the two-motor drive system is shown in Figure 6.Figure 6
NNGI-SIM control diagram.In Figure5, the function of the construction in dashed box can be expressed as(32)
G
n
c
s
=
G
c
s
1
+
G
c
s
G
m
s
.According to (22) and (32), G
c
(
s
) can be calculated as(33)
G
c
s
=
G
m
-
1
s
F
s
1
-
F
s
.In a two-motor NNGI system, the original system is separated into two subsystems. For the speed subsystem, the transfer function is G
11
=
1
/
(
s
+
1
) and the order of filter n is one. The expression of G
c
1
(
s
) can be transformed into the structure of classical PI controller, in which K
P
1 is 1
/
λ
1 and T
I
1 is 1. According to (23) and (33), G
c
1
(
s
) can be written as(34)
G
c
1
s
=
G
11
-
1
s
F
1
s
1
-
F
1
s
=
1
λ
1
1
+
1
s
.For the tension subsystem, according to (14), the transfer function is G
22
=
1
/
(
s
2
+
1.414
s
+
1
) and the order of filter n is one; then G
c
2
(
s
) can be written as(35)
G
c
2
s
=
G
22
-
1
s
F
2
s
1
-
F
2
s
=
s
2
+
1.414
s
+
1
λ
2
s
.The expression of G
c
2
(
s
) can be transformed into the structure of classical PID controller:(36)
G
c
2
s
=
1.414
λ
2
1
+
1
1.414
s
+
1
1.414
s
,where K
P
2 is 1.414
/
λ
2, T
I
2 is 1.414, and T
D
2 is 1
/
1.414.


## 3.1. NNGI System

Generalized inverse system can realize the inverse mapping relationship from the output to the input of original system. By setting the original order, derivatives and their linear combination of expected output of original system can be used as the input of the generalized inverse system. Then, the output can be used as the input of original system to drive the original system, generating the expected output [22].As a MIMO nonlinear system which hasq-dimensional input vector u
=
u
1
,
u
2
,
…
,
u
q
T
∈
R
q and q-dimensional output vector y
=
y
1
,
y
2
,
…
,
y
q
T
∈
R
q, the differential equation can be expressed as(2)
F
y
ε
T
Y
,
u
σ
T
,
U
=
0
,where(3)
ε
=
ε
1
,
…
,
ε
q
T
,
σ
=
σ
1
,
…
,
σ
q
T
,
y
ε
=
y
1
ε
1
,
…
,
y
q
ε
q
,
u
σ
=
u
1
σ
1
,
…
,
u
q
σ
q
,
Y
=
y
1
ε
1
-
1
,
…
,
y
1
,
…
,
y
q
ε
q
-
1
,
…
,
y
q
T
,
U
=
u
1
σ
1
-
1
,
…
,
u
1
,
…
,
u
q
σ
q
-
1
,
…
,
u
q
T
.
ε
j and σ
j are the highest derivative orders of input u
j and output y
j.If the relative order of the systemα
=
α
1
,
α
2
,
…
,
α
q
T is existent and the sum of relative order ∑
j
=
1
q
α
j
=
n (n is the order of the system), the state equation of the system can be transformed into differential equation. Hence, when the relative order is equal to the original order, the generalized inverse of the system can be directly obtained as(4)
u
=
ϕ
¯
Y
¯
,
v
¯
,where(5)
Y
¯
=
y
1
,
y
1
1
,
…
,
y
1
α
1
-
1
,
…
,
y
q
,
y
q
1
,
…
,
y
q
α
q
-
1
T
,
ν
¯
=
v
1
¯
,
v
2
¯
,
…
,
v
q
¯
T
,
v
j
¯
=
a
j
0
y
j
+
a
j
1
y
j
1
+
⋯
+
a
j
α
j
-
1
y
j
α
j
-
1
+
a
j
α
j
y
j
α
j
.
j
=
1
~
q
.By connecting the generalized inverse system in series, the coupled original system can be decoupled into several generalized pseudolinear composite subsystems and the transfer function of the subsystem can be expressed as(6)
G
s
=
d
i
a
g
G
11
s
,
…
,
G
q
q
s
=
d
i
a
g
1
a
1
α
1
s
α
1
+
a
1
α
1
-
1
s
α
1
-
1
+
⋯
+
a
11
s
+
a
10
,
…
,
1
a
q
α
q
s
α
q
+
a
q
α
q
-
1
s
α
q
-
1
+
⋯
+
a
q
1
s
+
a
q
0
.In practice, the structure of MIMO nonlinear system is complicated and it is hard to obtain its accurate mathematical model. Even if the existence of generalized inverse system has been proved, the expression of generalized inverse of the system cannot be derived for various occasions. To solve this problem, an ANN is used to approximate the original system. This method is independent of the accurate mathematical model, thus enhancing the popularization of inverse system algorithm. The NNGI pseudolinear composite subsystem equivalent chart is shown in Figure2.Figure 2
NNGI pseudolinear composite subsystem equivalent chart.


## 3.2. NNGI Construction for Two-Motor Drive System

According to the mathematic model of two-motor drive system in (1), the existence of generalized inverse system will be proved. Variables are chosen as follows.State variable is(7)
x
=
x
1
,
x
2
,
x
3
T
=
ω
r
1
,
ω
r
2
,
F
T
;control variable is(8)
u
=
u
1
,
u
2
T
=
ω
1
,
ω
2
T
;output variable is(9)
y
=
y
1
,
y
2
T
=
ω
r
1
,
F
T
.Then, the mathematic model of two-motor drive system can be expressed as(10)
x
˙
=
f
x
,
u
=
n
p
1
J
1
u
1
-
x
1
n
p
1
T
r
1
L
r
1
ψ
r
1
2
-
T
L
1
+
r
1
x
3
n
p
2
J
2
u
2
-
x
2
n
p
2
T
r
2
L
r
2
ψ
r
2
2
-
T
L
2
-
r
2
x
3
K
T
1
n
p
1
r
1
k
1
x
1
-
1
n
p
2
r
2
k
2
x
2
-
x
3
T
.According to the generalized inverse theory, the corresponding Jacobi matrix can be described as(11)
J
x
,
u
=
∂
y
1
1
∂
u
1
∂
y
1
1
∂
u
2
∂
y
2
2
∂
u
1
∂
y
2
2
∂
u
2
=
n
p
1
2
T
r
1
ψ
r
1
2
J
1
L
r
1
0
K
r
1
k
1
n
p
1
T
r
1
ψ
r
1
2
T
J
1
L
r
1
-
K
r
2
k
2
n
p
2
T
r
2
ψ
r
2
2
T
J
2
L
r
2
,
Det
J
x
,
u
=
K
r
2
k
2
n
p
1
2
n
p
2
T
r
1
T
r
2
ψ
r
1
2
ψ
r
2
2
T
J
1
J
2
L
r
1
L
r
2
.Whenψ
r
1
≠
0 and ψ
r
2
≠
0, the Jacobi matrix is nonsingular. The relative order of the system is α
=
(
α
1
,
α
2
)
=
(
1,2
), α
1
+
α
2
=
3, and it is equal to the order of the system. In addition, the original order of the system is n
e
=
(
n
e
1
,
n
e
2
)
=
(
1,2
), n
e
1
+
n
e
2
=
3, and it is the same as relative order. So, the inverse of the system is existent and can be expressed as(12)
u
=
ϕ
-
y
1
,
y
˙
1
,
y
2
,
y
˙
2
,
y
¨
2
,
v
-
,where(13)
v
¯
=
v
1
¯
,
v
2
¯
T
,
v
1
¯
=
a
10
y
1
+
a
11
y
1
1
,
v
2
¯
=
a
20
y
2
+
a
21
y
2
1
+
a
22
y
2
2
.In (12), to keep the stability of the open-loop system, a
10
=
1, a
11
=
1; a
20
=
1, a
21
=
1.414, and a
22
=
1. Then, the transfer function of the subsystem is(14)
G
s
=
d
i
a
g
G
11
,
G
22
=
d
i
a
g
1
s
+
1
,
1
s
2
+
1.414
s
+
1
.By connecting the generalized inverse system in series, the original system can be decoupled and simplified into a first-order subsystem and a second-order subsystem. As a result, the composite pseudolinear system is derived and shown in Figure3.Figure 3
Pseudolinear composite system.Simultaneously, by adjusting the parametersa
10, a
11, a
20, a
21, and a
22, the zeros and poles of pseudolinear subsystems can be developed exactly. Two generalized pseudolinear subsystems have neither right half-plane zero, and then the original system will obtain open-loop stabilization after it is in series with the NNGI system. The NNGI open-loop control diagram is shown in Figure 4.Figure 4
NNGI open-loop control diagram.


## 3.3. Profibus-DP Network Delay

Profibus communication uses the main token polling mechanism; each site on the bus connects with another in the bus way. Master station is composed of control center such as computer and PLC, while slave station consists of actuators such as inverter and sensors. Each site on the Profibus has its own logical address; a loop called logical loop is formed between the master stations. A token generated during the operation of the frame will poll along the logical link between the master stations. One master station is allowed to poll along other slave stations only when a token reaches this master station, and then the data transmission will be finished. While the process from sending task information of site to receiving the task information by target is too time intensive, the Profibus-DP network delay is generated.Profibus-DP network delay is indeterminate in practice. However the maximum delay can be calculated and it can be a reference for the delay of ideal internal modele
-
m
s [23, 24]. The approximate expression can be written as(15)
R
h
=
B
+
i
n
t
n
h
n
+
1
·
T
T
R
+
T
c
h
m
+
τ
+
ξ
h
,where B stands for the longest initial blocking, n
h is the number of high priority request messages added to the outgoing queue at critical instant, T
T
R is the target token rotation time, T
c
h
m is the maximum time of executing a high priority message transmission, τ is the token latency, n
+
1 is the number of high priority messages which are transmitted during the maximum token holding time, and ξ
h is the time between the finish of the last maximum token holding time and the end of target message transmission. There are three possible cases and the expression can be written as(16)
ξ
h
=
-
τ
,
n
h
=
i
n
t
n
h
n
+
1
·
n
+
1
,
T
c
h
m
,
n
h
=
i
n
t
n
h
n
+
1
·
n
+
1
+
1
,
n
h
-
i
n
t
n
h
n
+
1
·
n
+
1
·
T
c
h
m
+
τ
,
n
h
>
i
n
t
n
h
n
+
1
·
n
+
1
+
1
.


## 3.4. SIM Control

According to the previous discussions, two-motor drive system is a time-varying delay system. In order to settle the time delay issue, a SIM control is introduced into the controller design. The SIM control incorporates the advantages of Smith estimate control and internal model control. This control strategy is a delay compensation method, which is to separate the delay part of controlled object and connect in series with the system transfer function. It is an effective way to inhibit the effect of time delay. Moreover, SIM control can also restrain the influence of unmeasured disturbances and improve system robustness.The structure of the SIM control is shown in Figure5, in which G
P
(
s
) and G
m
(
s
) are the actual object model and the prediction internal model, respectively. e
-
P
s and e
-
m
s are the related delay of controlled object and prediction internal model, respectively. G
c
(
s
) is the controller. D
(
s
) is the external disturbance. R
(
s
) and Y
(
s
) are the input and the output of the system, respectively.Figure 5
SIM control diagram.


### 3.4.1. Delay Compensation Analysis

According to Figure5, the output can be written as(17)
Y
s
=
G
c
s
G
P
s
e
-
P
s
1
+
G
c
s
G
m
s
+
G
c
s
G
P
s
e
-
P
s
-
G
m
s
e
-
m
s
R
s
+
1
+
G
c
s
G
m
s
-
G
c
s
G
m
s
e
-
m
s
1
+
G
c
s
G
m
s
+
G
c
s
G
P
s
e
-
P
s
-
G
m
s
e
-
m
s
D
s
.The closed-loop transfer function of control system is(18)
G
s
=
Y
s
R
s
=
G
c
s
G
P
s
e
-
P
s
1
+
G
c
s
G
m
s
+
G
c
s
G
P
s
e
-
P
s
-
G
m
s
e
-
m
s
.When the prediction internal model is accurate,G
m
s
e
-
m
s is equal to G
P
s
e
-
P
s, and the transfer function can be written as(19)
G
s
=
G
c
s
G
P
s
1
+
G
c
s
G
P
s
e
-
P
s
.Hence, there is no delay part in the denominator of closed-loop transfer function and the delay is totally compensated.


### 3.4.2. Fixed Value Tracking Ability

According to Figure5, when the external disturbance D
(
s
) is zero, then the error of closed-loop system can be written as(20)
E
s
=
R
s
-
Y
s
=
1
-
G
n
c
s
G
m
s
e
-
m
s
1
+
G
n
c
s
G
P
s
e
-
P
s
-
G
m
s
e
-
m
s
R
s
.In consideration of modeling error of actual system, the controlled object can be described with multiplicative uncertainty:(21)
G
P
s
e
-
P
s
=
G
m
s
e
-
m
s
1
+
l
m
s
,where l
m
(
s
) represents the uncertainty and it is bounded.The structure in dashed box is represented byG
n
c
(
s
). A low pass filter is added to improve the capacity of resisting disturbance. The filter is represented as F
(
s
), and then G
n
c
(
s
) can be written as(22)
G
n
c
s
=
G
m
-
1
s
F
s
.Therefore, the transfer function of low pass filter is(23)
F
s
=
1
λ
s
+
1
n
,where λ describes filter parameter and n stands for the order of filter.Then, the error of closed-loop system can be simplified as(24)
E
s
=
λ
s
+
1
n
-
e
-
m
s
λ
s
+
1
n
+
l
m
s
e
-
m
s
R
s
.The limit of l
m
(
s
) can be expressed as(25)
lim
s
→
0
⁡
l
m
s
=
C
,where C is a constant.When the input signal is step signal,(26)
R
s
=
R
0
s
,then, the limit of E
(
s
) is(27)
e
∞
=
lim
s
→
0
⁡
s
E
s
=
lim
s
→
0
⁡
s
λ
s
+
1
n
-
e
-
m
s
λ
s
+
1
n
+
l
m
s
e
-
m
s
·
R
0
s
=
0
.


### 3.4.3. Anti-Interference Ability

According to Figure5, when the input R
(
s
) is zero, then the error of closed-loop system can be written as(28)
E
s
=
R
s
-
Y
s
=
1
-
G
n
c
s
G
m
s
e
-
m
s
1
+
G
n
c
s
G
P
s
e
-
P
s
-
G
m
s
e
-
m
s
-
D
s
.Similar to the analysis of fixed value tracking ability, the error of closed-loop system can be simplified as(29)
E
s
=
λ
s
+
1
n
-
e
-
m
s
λ
s
+
1
n
+
l
m
s
e
-
m
s
-
D
s
.When the external disturbance is step signal,(30)
D
s
=
D
0
s
,then, the limit of E
(
s
) is(31)
e
∞
=
lim
s
→
0
⁡
s
E
s
=
lim
s
→
0
⁡
s
-
λ
s
+
1
n
-
e
-
m
s
λ
s
+
1
n
+
l
m
s
e
-
m
s
·
D
0
s
=
0
.In conclusion, the SIM control structure can satisfactorily solve delay problem. What is more, this control structure also has good input signal tracking ability and anti-interference ability.


## 3.4.1. Delay Compensation Analysis

According to Figure5, the output can be written as(17)
Y
s
=
G
c
s
G
P
s
e
-
P
s
1
+
G
c
s
G
m
s
+
G
c
s
G
P
s
e
-
P
s
-
G
m
s
e
-
m
s
R
s
+
1
+
G
c
s
G
m
s
-
G
c
s
G
m
s
e
-
m
s
1
+
G
c
s
G
m
s
+
G
c
s
G
P
s
e
-
P
s
-
G
m
s
e
-
m
s
D
s
.The closed-loop transfer function of control system is(18)
G
s
=
Y
s
R
s
=
G
c
s
G
P
s
e
-
P
s
1
+
G
c
s
G
m
s
+
G
c
s
G
P
s
e
-
P
s
-
G
m
s
e
-
m
s
.When the prediction internal model is accurate,G
m
s
e
-
m
s is equal to G
P
s
e
-
P
s, and the transfer function can be written as(19)
G
s
=
G
c
s
G
P
s
1
+
G
c
s
G
P
s
e
-
P
s
.Hence, there is no delay part in the denominator of closed-loop transfer function and the delay is totally compensated.


## 3.4.2. Fixed Value Tracking Ability

According to Figure5, when the external disturbance D
(
s
) is zero, then the error of closed-loop system can be written as(20)
E
s
=
R
s
-
Y
s
=
1
-
G
n
c
s
G
m
s
e
-
m
s
1
+
G
n
c
s
G
P
s
e
-
P
s
-
G
m
s
e
-
m
s
R
s
.In consideration of modeling error of actual system, the controlled object can be described with multiplicative uncertainty:(21)
G
P
s
e
-
P
s
=
G
m
s
e
-
m
s
1
+
l
m
s
,where l
m
(
s
) represents the uncertainty and it is bounded.The structure in dashed box is represented byG
n
c
(
s
). A low pass filter is added to improve the capacity of resisting disturbance. The filter is represented as F
(
s
), and then G
n
c
(
s
) can be written as(22)
G
n
c
s
=
G
m
-
1
s
F
s
.Therefore, the transfer function of low pass filter is(23)
F
s
=
1
λ
s
+
1
n
,where λ describes filter parameter and n stands for the order of filter.Then, the error of closed-loop system can be simplified as(24)
E
s
=
λ
s
+
1
n
-
e
-
m
s
λ
s
+
1
n
+
l
m
s
e
-
m
s
R
s
.The limit of l
m
(
s
) can be expressed as(25)
lim
s
→
0
⁡
l
m
s
=
C
,where C is a constant.When the input signal is step signal,(26)
R
s
=
R
0
s
,then, the limit of E
(
s
) is(27)
e
∞
=
lim
s
→
0
⁡
s
E
s
=
lim
s
→
0
⁡
s
λ
s
+
1
n
-
e
-
m
s
λ
s
+
1
n
+
l
m
s
e
-
m
s
·
R
0
s
=
0
.


## 3.4.3. Anti-Interference Ability

According to Figure5, when the input R
(
s
) is zero, then the error of closed-loop system can be written as(28)
E
s
=
R
s
-
Y
s
=
1
-
G
n
c
s
G
m
s
e
-
m
s
1
+
G
n
c
s
G
P
s
e
-
P
s
-
G
m
s
e
-
m
s
-
D
s
.Similar to the analysis of fixed value tracking ability, the error of closed-loop system can be simplified as(29)
E
s
=
λ
s
+
1
n
-
e
-
m
s
λ
s
+
1
n
+
l
m
s
e
-
m
s
-
D
s
.When the external disturbance is step signal,(30)
D
s
=
D
0
s
,then, the limit of E
(
s
) is(31)
e
∞
=
lim
s
→
0
⁡
s
E
s
=
lim
s
→
0
⁡
s
-
λ
s
+
1
n
-
e
-
m
s
λ
s
+
1
n
+
l
m
s
e
-
m
s
·
D
0
s
=
0
.In conclusion, the SIM control structure can satisfactorily solve delay problem. What is more, this control structure also has good input signal tracking ability and anti-interference ability.


## 3.5. NNGI-SIM Control

As mentioned above, SIM control can remove the delay of forward network from the closed-loop and set it as a gain block before the output. As a result, the delay on the return network path can be totally eliminated. Meanwhile, SIM control can also enhance robustness and anti-interference ability of the system with the internal model control structure. However, the accuracy of prediction modelG
m
(
s
) is crucial. When there is a mismatch between prediction model and actual model, the control quality will be significantly deteriorated. This problem can be solved by NNGI system. After connecting with NNGI system, the original system is decoupled into two subsystems, namely, a first-order speed one and a second-order tension one. The accurate mathematical model of two subsystems can be obtained after pole-zero allocation of the pseudolinear subsystems, and then the internal model can be predicted. Then, the mismatch condition will be avoided. The diagram of NNGI-SIM control strategy of the two-motor drive system is shown in Figure 6.Figure 6
NNGI-SIM control diagram.In Figure5, the function of the construction in dashed box can be expressed as(32)
G
n
c
s
=
G
c
s
1
+
G
c
s
G
m
s
.According to (22) and (32), G
c
(
s
) can be calculated as(33)
G
c
s
=
G
m
-
1
s
F
s
1
-
F
s
.In a two-motor NNGI system, the original system is separated into two subsystems. For the speed subsystem, the transfer function is G
11
=
1
/
(
s
+
1
) and the order of filter n is one. The expression of G
c
1
(
s
) can be transformed into the structure of classical PI controller, in which K
P
1 is 1
/
λ
1 and T
I
1 is 1. According to (23) and (33), G
c
1
(
s
) can be written as(34)
G
c
1
s
=
G
11
-
1
s
F
1
s
1
-
F
1
s
=
1
λ
1
1
+
1
s
.For the tension subsystem, according to (14), the transfer function is G
22
=
1
/
(
s
2
+
1.414
s
+
1
) and the order of filter n is one; then G
c
2
(
s
) can be written as(35)
G
c
2
s
=
G
22
-
1
s
F
2
s
1
-
F
2
s
=
s
2
+
1.414
s
+
1
λ
2
s
.The expression of G
c
2
(
s
) can be transformed into the structure of classical PID controller:(36)
G
c
2
s
=
1.414
λ
2
1
+
1
1.414
s
+
1
1.414
s
,where K
P
2 is 1.414
/
λ
2, T
I
2 is 1.414, and T
D
2 is 1
/
1.414.


## 4. Simulation

### 4.1. Procedure of Simulation

In order to verify the decoupling ability of the proposed NNGI system, feed-forward neural network is adopted to approximate inversion system. The design procedure of NNGI pseudolinear composite system is briefly given as follows.


#### 4.1.1. Sample Collection

Since the parameters of the original system model are unknown or variable, the selected signals should stimulate the dynamic and static characters of original system sufficiently. So, the square wave signal is chosen to verify dynamic performance of the proposed motor drive in each frequency band. Differentiators in the NNGI system can be obtained by using five-point numerical derivative algorithm, in which high computing accuracy is offered. The speed value ranges from 100 r/min to 600 r/min and tension value ranges from 100 N to 600 N.


#### 4.1.2. Neural Network Training

According to (12), y
1
,
y
1
+
y
˙
1
,
y
2
,
y
˙
2
,
y
¨
2
+
1.414
y
˙
2
+
y
2 are chosen as the input and u is chosen as the output to train the NNGI system. Before training neural network, the sampled data should be normalized and break off both ends to improve convergence precision. Furthermore, the neural network adopts 5-15-2 3-layer structure; the activation functions of hidden layer and output layer choose tansig and pureline functions, respectively. The weight and the threshold of the neural network should be adjusted by using the Levenberg-Marquardt algorithm until convergence precision attains the expectant effect. The output of neural network should be antinormalized back to the original range. Then, the training of NNGI system is over.


#### 4.1.3. Integration of Pseudolinear Composite System

After putting the NNGI system in front of the original system, the pseudolinear composite system is constructed. A strong coupling and nonlinear system is separated into two pseudolinear composite subsystems. By adopting the SIM control structure, two independent control loops are formed to achieve the aim of decoupling control.A MIMO nonlinear system is decoupled by NNGI method, and then the accurate model of generalized pseudolinear subsystem can be obtained, in which the controllers’ parameters depend on the accurate mathematical model. The robustness of the controllers is mainly determined by the filters, and the first-order filter is considered as the best controller. After a lot of simulation experiments, the filter parametersλ
1 and λ
2 are chosen to be 0.2 and 1.2, respectively, and then the filters can meet the system requirements. Hence, the controller parameters of two pseudolinear subsystems can be optimized. According to the analysis of previous discussions, the expressions of G
c
1 and G
c
2 can be written as the classical PID controller. According to (34) and (36), K
P
1 is 5 and T
I
1 is 1 in the speed subsystem, K
P
2 is 1.2, T
I
2 is 0.8, and T
D
2 is 0.8 in the tension subsystem. The simulation model is developed according to Figure 6. Based on the NNGI-SIM control strategy, the coupling issue can be settled and the delay of Profibus-DP network can also be compensated.


### 4.2. Comparative Discussions

In order to evaluate the proposed NNGI-SIM control system,ω
r
1 is set to track square wave response and F is set to constant 300 N. ω
r
1 increases from 200 r/min to 250 r/min at 80 s and decreases from 250 r/min to 200 r/min at 140 s. As shown in Figure 7, due to the effect of system delay, the responses of traditional PID control have high frequency shiver during the whole simulation. Large fluctuations of the tension responses occur at the constant of speed sudden changing. Therefore, it cannot effectively decouple speed and tension. Moreover, the response speed is slow and transient time is too long. By contrast, with the NNGI-SIM control structure, the shiver problem has been solved. Furthermore, the tension fluctuation caused by the sudden change of speed is restrained to a great extent. The disturbance caused by speed impact is significantly reduced. Hence, this two-motor system is effectively coupled; the speed of motors and tension will not be influenced by each other.Figure 7
Simulated responses whenω
r
1 suddenly changes.
(a)
PID control strategy
(b)
NNGI-SIM control strategyIn order to evaluate the restraining disturbance ability of the proposed control strategy,ω
r
1 and F are set to be constants of 300 r/min and 300 N, respectively. Figure 8 is the sudden load increase from 2 Nm to 5 Nm at 100 s, while Figure 9 is the sudden load decrease from 5 Nm to 2 Nm at 100 s. At the PID control operation, speed and tension have large fluctuations and both of them are significantly affected by load jump. In contrast, at the proposed control operation, the fluctuation of speed is smoothed and tension is immune to this load change, indicating that the proposed system has a strong antidisturbance ability. Furthermore, the fluctuation issue caused by transmission delay has been settled completely. Since SIM control structure can separate the delay part of controlled object and connect in series with the system transfer function, the operation of system is immune to delay. So, the tension responses by NNGI-SIM control are relatively smooth and the delay of SIM control structure is obviously compensated.Figure 8
Simulated responses for sudden increased load impact.
(a)
Speed responses
(b)
Tension responsesFigure 9
Simulated responses for sudden decreased load impact.
(a)
Speed responses
(b)
Tension responses


## 4.1. Procedure of Simulation

In order to verify the decoupling ability of the proposed NNGI system, feed-forward neural network is adopted to approximate inversion system. The design procedure of NNGI pseudolinear composite system is briefly given as follows.


### 4.1.1. Sample Collection

Since the parameters of the original system model are unknown or variable, the selected signals should stimulate the dynamic and static characters of original system sufficiently. So, the square wave signal is chosen to verify dynamic performance of the proposed motor drive in each frequency band. Differentiators in the NNGI system can be obtained by using five-point numerical derivative algorithm, in which high computing accuracy is offered. The speed value ranges from 100 r/min to 600 r/min and tension value ranges from 100 N to 600 N.


### 4.1.2. Neural Network Training

According to (12), y
1
,
y
1
+
y
˙
1
,
y
2
,
y
˙
2
,
y
¨
2
+
1.414
y
˙
2
+
y
2 are chosen as the input and u is chosen as the output to train the NNGI system. Before training neural network, the sampled data should be normalized and break off both ends to improve convergence precision. Furthermore, the neural network adopts 5-15-2 3-layer structure; the activation functions of hidden layer and output layer choose tansig and pureline functions, respectively. The weight and the threshold of the neural network should be adjusted by using the Levenberg-Marquardt algorithm until convergence precision attains the expectant effect. The output of neural network should be antinormalized back to the original range. Then, the training of NNGI system is over.


### 4.1.3. Integration of Pseudolinear Composite System

After putting the NNGI system in front of the original system, the pseudolinear composite system is constructed. A strong coupling and nonlinear system is separated into two pseudolinear composite subsystems. By adopting the SIM control structure, two independent control loops are formed to achieve the aim of decoupling control.A MIMO nonlinear system is decoupled by NNGI method, and then the accurate model of generalized pseudolinear subsystem can be obtained, in which the controllers’ parameters depend on the accurate mathematical model. The robustness of the controllers is mainly determined by the filters, and the first-order filter is considered as the best controller. After a lot of simulation experiments, the filter parametersλ
1 and λ
2 are chosen to be 0.2 and 1.2, respectively, and then the filters can meet the system requirements. Hence, the controller parameters of two pseudolinear subsystems can be optimized. According to the analysis of previous discussions, the expressions of G
c
1 and G
c
2 can be written as the classical PID controller. According to (34) and (36), K
P
1 is 5 and T
I
1 is 1 in the speed subsystem, K
P
2 is 1.2, T
I
2 is 0.8, and T
D
2 is 0.8 in the tension subsystem. The simulation model is developed according to Figure 6. Based on the NNGI-SIM control strategy, the coupling issue can be settled and the delay of Profibus-DP network can also be compensated.


## 4.1.1. Sample Collection

Since the parameters of the original system model are unknown or variable, the selected signals should stimulate the dynamic and static characters of original system sufficiently. So, the square wave signal is chosen to verify dynamic performance of the proposed motor drive in each frequency band. Differentiators in the NNGI system can be obtained by using five-point numerical derivative algorithm, in which high computing accuracy is offered. The speed value ranges from 100 r/min to 600 r/min and tension value ranges from 100 N to 600 N.


## 4.1.2. Neural Network Training

According to (12), y
1
,
y
1
+
y
˙
1
,
y
2
,
y
˙
2
,
y
¨
2
+
1.414
y
˙
2
+
y
2 are chosen as the input and u is chosen as the output to train the NNGI system. Before training neural network, the sampled data should be normalized and break off both ends to improve convergence precision. Furthermore, the neural network adopts 5-15-2 3-layer structure; the activation functions of hidden layer and output layer choose tansig and pureline functions, respectively. The weight and the threshold of the neural network should be adjusted by using the Levenberg-Marquardt algorithm until convergence precision attains the expectant effect. The output of neural network should be antinormalized back to the original range. Then, the training of NNGI system is over.


## 4.1.3. Integration of Pseudolinear Composite System

After putting the NNGI system in front of the original system, the pseudolinear composite system is constructed. A strong coupling and nonlinear system is separated into two pseudolinear composite subsystems. By adopting the SIM control structure, two independent control loops are formed to achieve the aim of decoupling control.A MIMO nonlinear system is decoupled by NNGI method, and then the accurate model of generalized pseudolinear subsystem can be obtained, in which the controllers’ parameters depend on the accurate mathematical model. The robustness of the controllers is mainly determined by the filters, and the first-order filter is considered as the best controller. After a lot of simulation experiments, the filter parametersλ
1 and λ
2 are chosen to be 0.2 and 1.2, respectively, and then the filters can meet the system requirements. Hence, the controller parameters of two pseudolinear subsystems can be optimized. According to the analysis of previous discussions, the expressions of G
c
1 and G
c
2 can be written as the classical PID controller. According to (34) and (36), K
P
1 is 5 and T
I
1 is 1 in the speed subsystem, K
P
2 is 1.2, T
I
2 is 0.8, and T
D
2 is 0.8 in the tension subsystem. The simulation model is developed according to Figure 6. Based on the NNGI-SIM control strategy, the coupling issue can be settled and the delay of Profibus-DP network can also be compensated.


## 4.2. Comparative Discussions

In order to evaluate the proposed NNGI-SIM control system,ω
r
1 is set to track square wave response and F is set to constant 300 N. ω
r
1 increases from 200 r/min to 250 r/min at 80 s and decreases from 250 r/min to 200 r/min at 140 s. As shown in Figure 7, due to the effect of system delay, the responses of traditional PID control have high frequency shiver during the whole simulation. Large fluctuations of the tension responses occur at the constant of speed sudden changing. Therefore, it cannot effectively decouple speed and tension. Moreover, the response speed is slow and transient time is too long. By contrast, with the NNGI-SIM control structure, the shiver problem has been solved. Furthermore, the tension fluctuation caused by the sudden change of speed is restrained to a great extent. The disturbance caused by speed impact is significantly reduced. Hence, this two-motor system is effectively coupled; the speed of motors and tension will not be influenced by each other.Figure 7
Simulated responses whenω
r
1 suddenly changes.
(a)
PID control strategy
(b)
NNGI-SIM control strategyIn order to evaluate the restraining disturbance ability of the proposed control strategy,ω
r
1 and F are set to be constants of 300 r/min and 300 N, respectively. Figure 8 is the sudden load increase from 2 Nm to 5 Nm at 100 s, while Figure 9 is the sudden load decrease from 5 Nm to 2 Nm at 100 s. At the PID control operation, speed and tension have large fluctuations and both of them are significantly affected by load jump. In contrast, at the proposed control operation, the fluctuation of speed is smoothed and tension is immune to this load change, indicating that the proposed system has a strong antidisturbance ability. Furthermore, the fluctuation issue caused by transmission delay has been settled completely. Since SIM control structure can separate the delay part of controlled object and connect in series with the system transfer function, the operation of system is immune to delay. So, the tension responses by NNGI-SIM control are relatively smooth and the delay of SIM control structure is obviously compensated.Figure 8
Simulated responses for sudden increased load impact.
(a)
Speed responses
(b)
Tension responsesFigure 9
Simulated responses for sudden decreased load impact.
(a)
Speed responses
(b)
Tension responses


## 5. Experimental Verification

The experimental platform of multimotor drive system is designed as shown in Figure10, by which the number of the test motors can be up to three. The whole system includes industrial computer, Siemens S7-300 PLC, Siemens micromaster vector inverters, 2.2 kW three-phase asynchronous motors, photoelectric encoder, magnetic powder brakes, and tension sensor.Figure 10
Experimental platform.In the experimental platform, S7-300 PLC is the main controller which sends the frequency signal to the inverters. In addition, the inverters drive motors in vector control mode and the motors connect conveyor belt with roller.The PLC controls two inverters through Profibus-DP communication which is the main way of message transmission delay. According to (15) and (16), this study sets the data transfer rate of Profibus-DP network as 1.5 Mbps, and then the data bit period is equal to 0.667 us. Besides, n
h is chosen to be 20 and T
T
R is fixed to 8 ms, and then B is equal to 2.18 ms, T
c
h
m is equal to 1.27 ms, τ is approximately equal to 0.04 ms, nis equal to 6, and ξ
h is equal to 7.66 ms in the third case.Hence, according to (15), the maximum delay of Profibus-DP network is(37)
R
h
=
B
+
i
n
t
n
h
n
+
1
·
T
T
R
+
T
c
h
m
+
τ
+
ξ
h
=
28.46
m
s
.The hardware experiment condition is similar to simulation one, of which the procedure can be divided into three parts, namely, sample collection, training offline, and control implementation. However, differently, the experimental motor will be affected by external disturbances and mechanical vibration, hence resulting in the fluctuations of response speed and tension. To eliminate these harmful effects, antipulse interference average filter is applied to filter the real-time sampled signals. This filter method collectsN sets of date continually, removes the maximum and the minimum, and takes the average of the remaining date as the output of the filter. The feedback signal after filtering has no significant fluctuation and is smooth. It is conducive to improve curve-fitting condition and convergence precision in neural network training.This work implements the feed-forward neural networks with instruction list of Siemens PLC. After transferring the trained weights and thresholds to the corresponding register units, the NNGI system is accomplished in S7-300 PLC platform. Then, by adopting the SIM control structure to deal with the problem of Profibus-DP network delay, the whole closed-loop system is formed.In order to verify the decoupling effectiveness, the tension supply of 300 N is chosen and the speed supply is set to suddenly change from 300 r/min to 400 r/min. The responses of the experiment system with PID control and with NNGI-SIM control are compared in Figure11. As shown in Figure 11(a), the tension response fluctuates dramatically when the speed suddenly changes. Thus, the strong coupling problem cannot be solved by the traditional PID control. Moreover, the transmission delay leads to the fluctuation of responses, which agrees with the simulated one. Comparatively, in Figure 11(b), with the sudden change of speed from 300 r/min to 400 r/min at 80 s, the jitter has been reduced, so the sudden change of speed has less influence on tension response. It is obvious that the tension and the speed are successfully decoupled by using the proposed control strategy. Hence, the NNGI-SIM control is superior to PID control in the two-motor drive system with network delay.Figure 11
Experimental results whenω
r
1 suddenly changes.
(a)
PID control strategy
(b)
NNGI-SIM control strategyThe measured responses of speed and tension under sudden load impact are shown in Figures12 and 13. In Figure 12, the motor starts at no load and then increases to 3 Nm suddenly. By contrast, in Figure 13, the load decreases from 3 Nm to 0 Nm.Figure 12
Experimental results for sudden increased load impact.
(a)
Speed responses
(b)
Tension responsesFigure 13
Experimental results for sudden decreased load impact.
(a)
Speed responses
(b)
Tension responsesAs shown in Figures12 and 13, the proposed control strategy reduces the impact on speed and tension caused by the sudden application of load. The fluctuation of tension is mitigated effectively under NNGI-SIM control. However, the improvement of speed response is relatively less because the magnetic powder brake has direct effects on motor shaft. Overall, the NNGI-SIM control strategy achieves the good dynamic and static performances, and it has strong robustness to the load disturbances.


## 6. Conclusion

In this paper, a new control strategy has been proposed to solve the Profibus-DP network delay problem and decouple the speed and the tension of the two-motor drive system by incorporating NNGI and SIM control. The NNGI strategy has been used to transform the MIMO two-motor drive system to several single-input single-output systems, and then the speed and tension are decoupled validly. By adopting the SIM control structure to compose a closed-loop system, the delay issue has been compensated. The designed system can provide good starting and tracking capabilities. The utilities of the proposed control strategy have been verified by both the simulation and the experiments at various operations. The proposed control strategy has a great application prospect in the field of multimotor control system.



---

*Source: 5794087-2016-02-28.xml*