Scholarly article on topic 'Adaptive complementary fuzzy self-recurrent wavelet neural network controller for the electric load simulator system'

Adaptive complementary fuzzy self-recurrent wavelet neural network controller for the electric load simulator system Academic research paper on "Mechanical engineering"

Share paper
Academic journal
Advances in Mechanical Engineering
OECD Field of science

Academic research paper on topic "Adaptive complementary fuzzy self-recurrent wavelet neural network controller for the electric load simulator system"

Advances in


Research Article Engineering

Adaptive complementary fuzzy self-recurrent wavelet neural network controller for the electric load simulator system

Advances in Mechanical Engineering 2016, Vol. 8(3) 1-12 © The Author(s) 2016 DOI: 10.1177/1687814016639250


Wang Chao, Gao Qiang, Hou Yuanlong, Hou Runmin and Min Hao


Due to the complexities existing in the electric load simulator, this article develops a high-performance nonlinear adaptive controller to improve the torque tracking performance of the electric load simulator, which mainly consists of an adaptive fuzzy self-recurrent wavelet neural network controller with variable structure (VSFSWC) and a complementary controller. The VSFSWC is clearly and easily used for real-time systems and greatly improves the convergence rate and control precision. The complementary controller is designed to eliminate the effect of the approximation error between the proposed neural network controller and the ideal feedback controller without chattering phenomena. Moreover, adaptive learning laws are derived to guarantee the system stability in the sense of the Lyapunov theory. Finally, the hard-ware-in-the-loop simulations are carried out to verify the feasibility and effectiveness of the proposed algorithms in different working styles.


Complementary controller, self-recurrent wavelet neural network, Lyapunov theory, electric load simulator

Date received: 30 July 2015; accepted: 22 January 2016 Academic Editor: Yaguo Lei


The electric load simulator (ELS) is a crucial equipment in the hardware-in-the-loop (HIL) experiments, widely used in the gun, aerospace, and marine engineering,1 3 which is designed to generate the torque to simulate the dynamic loads acting on the actuator system and greatly accelerates the development of the research and product of the actuator system. The ELS of steering engines is commonly studied, which is different from the ELS for the gun control. The former is designed to simulate resisting torque, such as the air and water resistance, the latter is used to simulate the resisting and inertia torque. Apart from the same resisting torque, the inertia torque is also time variable, which can be effectively simulated comparatively, because the variations of the barrel angle and the number of the shells

can be computed. In a word, the nonlinearities and uncertainties are playing great negative effects on the performance of the ELS, like the coupling of the torque motor and the load motor, the friction and backlash.4,5 Therefore, how to improve the ELS torque tracking performance has become a research hotspot recently. The ordinary methods adopt the feed-forward compensation using the actuator's position and velocity signals,6,7 or the dual-loop scheme,8 but the constructed

School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing, China

Corresponding author:

Wang Chao, School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210014, China. Email:

|(ccj \V Creative Commons CC-BY: This article is distributed under the terms of the Creative Commons Attribution 3.0 License

( which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages ( open-access-at-sage).

model cannot wholly describe the friction and backlash nonlinearities, especially the approximation precision shows worse in the practical working conditions, so intelligent controllers with nonlinear characteristics and high precisions have been proposed to substitute the mechanism model.

The feed-forward neural networks (NN) with amazing learning ability, parallel computation, and remarkable generalization ability are widely applied to approximate the nonlinear and uncertain system, but they also have the drawbacks of the local minimum and low convergence rate.9,10 In order to solve these problems, the activation functions and the learning structures of the NN must be selected appropriately, which keep the model with the high computation efficiency and short training time. The wavelet functions with the advantages of time-frequency localization are incorporated into the NN, and the constructed wavelet neural networks (WNN) are widely applied to the engineering research area as the function approximation and signal processing.11,12 Considering the sudden changes and the external disturbance of the ELS, the self-recurrent WNN (SW) with storing the previous data is adopted to maintain the system stability.13,14 Meantime, considering the Takagi-Sugeno-Kang (TSK) fuzzy model with SW has generated a new network named fuzzy SW (FSW), which does well in enhancing the function approximation accuracy and generalization ability.15 However, it is hard to determine the learning structure of the FSW because of the trade-off between the learning performance and the number of neurons. If the number of neurons chosen is too large, the computation loading may be too heavy for the practical applications although the learning performance is very good. To avoid this problem, many structure learning algorithms are proposed for the adaption of the NN,16,17 and the number of neurons can be adjusted dynamically by the significant index or rules. In practical control applications, it is impatient duty to have a systematic method of ensuring the stability of the overall system. Recently, many NN methods are introduced based on the Lyapunov stability the-ory,16 19 but what we need to resolve and focus on is the determination of the assumptions and adaptive learning laws. With the above-mentioned motivations, the radial basis function (RBF) NN with a self-recurrent and variable structure (VS) not only guarantees the stability of the control system,19 but also does not need the constrained conditions and knowledge of the controlled system.20,21 In order to further decrease the approximation error, many compensation methods are proposed aiming at eliminating the external disturbances, which reduce the approximation error and improve the control precision.22,23

In this article, the proposed VSFSWC and complementary controller (CC) is a novel control structure,

which comprises the merits of the fuzzy, self-recurrent, VS, WNN, and compensation method for the ELS system. The VSFSWC used as a main controller is online estimating the dynamic system, and the CC dispels the effect of the approximation error existing in the ELS system. Furthermore, the VS learning algorithm effectively decreases the computation loading, and the adaptive learning laws based on the Lyapunov theory are developed to tune the parameters of the WNN. Not only the fast convergence of the torque tracking error can be achieved but also the system stability can be guaranteed. Finally, the proposed control algorithms are comparatively applied to the HIL simulation system for the ELS, and the experimental results show that the proposed VSFSWC and CC can achieve a favorable control performance.

Problem statement and ideal backstepping controller design

Problem statement

Consider the characteristics of the ELS system which are described as the nth order nonlinear equation in the


yn(t)=f( y, t) + g(y, t) u(t) + d(t) (1)

where y(t) and u(t) denote the control input and output of the system, respectively; yn(t) is the nth order derivative of y(t); f (y, t) and g(y, t) are the system dynamic variables; d(t) is the unknown but bounded external disturbances, which is introduced by the self-coupling, friction, and backlash.

Ideal backstepping controller design

The control objective is to find a control law so that the state y(t) can track the reference command yd (t) closely. The ideal backstepping controller can be constructed step-by-step as follows:

Step 1. The tracking error is defined as follows

ei(t)=yi (t) - yd(t) (2)

Then, a stabilization function is as follows

si(t) = yd(t) - liei(t) (3)

where 1i is a positive constant.

Step 2. The second control input e2 (t) is expressed as follows

e2(t)=y2(t) - si(t) (4)

S2(t) = si(t) - ei(t) - l2e2(t) (5)

where A2 is a positive constant; thus, the derivative of the error e1 (t) is obtained as follows

ei(0 = yi(t) - yd(t) = yi(t) - si(i) - Aiei(i) - e2(t) - liei(t)

x(y, t)- - o-n-iO)

i -■

g(y, t)

-<t)+/gy! (12)

- - x(y, t) - e„-i(t) - A„e„(t)

en(t) - - e„-i(t) - A„e„(t)

Vi (t) -

_ e2(t) + e2(t)


Then the derivative of V (t) will be

Vi(t) - ei(t)ei (t) + e2(t)e2(t)+ ••• + e„(t)e„(t)

- ei(t)(e2(t) - Aiei(t)) + e2(t)(e3(t) - A2ß2(t) - ei(t)) + ••• + e„(t)( - e„-i(t) - A„e„(t))

- - Aiei(t) - A2e2(t)-----A^t) < 0

Step i (3 < i < n — 1). An error e,(t) and a stabilization function si(t) are defined

et (t)^y(t) — si—i(t) (7)

si (t) = s i— i (t) — ei—i (t) — Aiei(t) (8)

where A,- is a positive constant; thus, the derivative of the error ei—1(t) is obtained as follows

ei—i(t) = yi— i(t) — si—2(t) = yi(t) — Si—i(t) — ei—2(t) — Ai—iei—i(t) (9) = ei(t) — a,-— iei—i(t) — ei—2(t)

Step n. An error en(t) and a stabilization function sn(t) are defined

en(t)=yn (t) — Sn—i(t) (10)

e n (t) = y n (t) — s n— 1 (t) = u(t) + x(y, t) (11)

where yn(t) is the nth order derivative of y(t), and the nonlinear function x(y, t) is in the following

Since )V1(t) is negative semi-definite, which means that the equation V1 (t) < V (0) can be achieved, it implies that e1(t), e2(t), ...,en(t) are bounded. Defining the equation S(r)[A1ej (t) + A2e2(t) + ••• + Anen(t) < — ¡V1(t) and integrating S(t) with respect to time, the following equation is obtained

S(r)dr < Vi(t) - Vi(0)

Because V1 (0) is bounded and V (t) is non-increasing and bounded, the result can be obtained as follows

S(r)dr < t

Step n + 1. The ideal backstepping controller can be defined

where An is a positive constant; substituting the ideal backstepping control law u = WjBC into (11), we can obtain the following equation

Step n + 2. Consider the candidate Lyapunov function in the following form

In addition, since E(t) is bounded, by Barbalat's Lemma, lim E(t) = 0.26 As a result, the asymptotically stability can be guaranteed in the sense of the Lyapunov theory. However, the dynamic variables f (y, t) and g(y, t) cannot be entirely known in practical applications. Moreover, the external disturbance d(t) causes an outcome of large chattering phenomena, which deteriorates the control precision and the robustness of the ELS system. Therefore, the ideal backstep-ping controller cannot be achieved precisely.

Adaptive complementary fuzzy self-recurrent WNN controller design

Self-recurrent WNN

The structure of the self-recurrent WNN controller (SWC) is shown in Figure 1, which is comprised of an input layer, wavelet layer, and output layer. The input of the input layer is e(t)= [e1 (t), e2(t), ...,en(t)]T, where T is the transpose and n is the number of dimensions.

Layer 1—Input layer: No function is performed in this layer; the node only transmits input values to layer 2.

Layer 2—Wavelet layer: In this layer, each activation function can be represented27

<x) - 2a=2f (2ax - b)

Figure 1. Structure of the fuzzy self-recurrent WNN.

where a and b are the dilation and translation parameters, respectively, which are self-learning based on the following adaptive laws. The Mexican hat mother wavelet function is defined as follows

U(x) = (1 - ||x||2)exp

where ||x|| = xTx. Thus, the activation function of the jth neuron connected with the input layer is described

f,.;.(e,.) = 2flj/2(l - ||2-

-jv-i / - y- ir -, b</|| jexp| i = 1,2, ..., n; j = 1, 2, ..., m where m is the number of the neurons.

Layer 3—Output layer: With the self-recurrent rules, the output of the SWC is as follows14

Jj = II f j + hoJj-i

where ho 2 [0,1] is the self-recurrent feedback gain. The self-recurrent WNN differs from the typical WNN by the self-feedback link with fixed coefficient h0 in the context node.

Fuzzy self-recurrent WNN

In comparison with the Mamdani, the fuzzy logic systems (FLS), and the TSK, the TSK model is more useful for the application of real-time operations.

However, the TSK cannot wholly reflect the characteristics of the ELS system with limited rules; the SWC is incorporated into the TSK model to improve the control rapidity and accuracy. A typical fuzzy WNN for approximating the nonlinear function can be depicted by the linguistic rule.11

If xi isAj, x2 isA*, ...,x, isAk, then xk = MkJk, k = i, 2, ..., m, where x, is the ith input variable of the system for i = i; Ak is the linguistic term characterized by a Gaussian-type membership function mk for i = i, which is calculated as follows

Ak = fi exP I

-(x, - c,k)2 d2

where cik and dik denote the center and scaling parameters of the membership function associated with rule k, then15

In summary, the output of the fuzzy self-recurrent WNN controller (FSWC) can be obtained

Ak gk = mT gfo b)

where m = ^ ..., Mm^ J = [Jl, J2, ..., Jm^ a = [an, ..., aim, ..., ani, ..., anm], and b =

[¿i^ . . . , bim, . . . , ¿nU . . . , bnm].

Considering the universal function approximation property of the FSWC, there exists an expansion of (25) so that it can uniformly approximate the ideal back-stepping controller28,29

x* = m g*(a*, b*) + D1

where Di is the approximation error; m*, 7*, a*, and b* are the ideal parameter vectors of m, g, a, and b, respectively. Since the ideal parameter vectors cannot be obtained, an estimation FSWC is defined as follows

x = m T g (â, b)

where m, 7, and b are the estimation parameter vectors of m, 7, a, and b, respectively. Subtracting equation (27) from equation (26), an estimation error is defined

X = X* — X = m*T g*(a*, b*) + D1

= m Tg + m T~ + m T~ + a1

m gfo b)

where m = m* — m and g = g* — -y. The linearization technique is employed to transform the membership functions into a partially linear form so that the expansion of ~ in a Taylor series can be obtained30

|fl — a(a*- â) +

9~j 9b

lb = b

(b*- b) +j = gfla + gb~ + j

«M e.W


A/ ' ¿



'WNlt J_l_

Variable Structure

Ti ^ b - -

Learning Laws

Complementary Controller

Torque sensor

where a = a* a and b = b* b; j is a vector of ... _ D1 . ,. ,/t-ct-,A/,-

' Figure 2. Block diagram of the VSFSWC and CC.

higher-order terms. ( 9bj/9a) and ( 9bj/9b) are defined

9"' 9a

(j-1) X n 1

■ ■, da1, fflffl

(m—1) X n

0 0 Q Q

(1—1) X n 1 1 (m—1) X n

gmin = min,. |g|

1 <j < m(k)

where n X m(q) is the number of the existing neurons at the qth control interval. If gmin is even bigger than a pre-specified and positive threshold gth, a new neuron should be added. The parameters of the new neuron layer are selected as follows

am(q + 1) — ao

bm(q + 1) — b0

mm(q + 1)

m(q + 1) — m(q) + 1

where a0 and b0 are pre-specified constants. In order to avoid the infinite growth of the number of neurons, the other significant index I(q) is adopted to evaluate the importance of the jth neuron at the qth control interval and defined as follows24

I(q + 1) —

/j(q)exp( - r)

if r7(q)>ro if r,(q) < go

Substituting equation (28) into equation (29) yields

x = mT (g + ga a + gbb + j) + mT (rfla + gbb + j) + Di = mT g + mT (gfla + gb~) + D2 (32)

where D2 = mTg + mTj + D1. D2 denotes the lump of the approximation error and is limited by |D2| < E0, where E0 is a positive constant.

VS learning algorithm

Considering the trade-off between the computation loading and the learning performance, an online VS learning algorithm is proposed to determine whether or not to add and prune the number of the neurons. The output of the SWC is represented as the index to add a neuron

where the initial value setting of j(q) is 1; g0 is the reduction threshold; t is the reduction speed constant. If I < is satisfied, where Ith is a pre-specified threshold, the jth neuron will be pruned.

Stability analysis

The developed VSFSWC and CC is shown in Figure 2, which is composed of a VSFSWC and a CC. The controller output is given as follows31

uvc = UVSFSWC + UCC = — x — en—1 (t) — Anen (t) + UCC

where the output of VSFSWC xa is designed to approximate the controlled system dynamics x(y, t); the output of CC uCC is designed to eliminate the effect of the approximation error. Substituting equation (36) into equation (11) yields

en(t) — x(y, t) - X - en—1(t) - Anßn(t) + wcc

— X - en-1(t) - Anen(t) + wcc

Then substituting equation (32) into equation (37), the equation (38) can be obtained as follows

en(t) = m T g + m T (7fla + gbb)

+ D2 — en—1(t) — Anen (t) + UCC

In this article, the CC is designed as follows

mcc = - £oSgn(e„)

where EE0 is the estimation bound of D2 and sgn() is a sign function.

In order to guarantee the system stability, the Lyapunov function candidate is given as follows

eft) ^ ~T ~

+ + - + f + 2h + aT ä + bTb + -

2h2 2h3 2^4

2 2 2 T i ä'

= - Aißj - A2ß2-----An^2 + ~ ( e„g + —

+ ä e«gflg +

f* = - ~ = hi^ng ä = - ä = h2engflg

1} = - ~ = h3engbg

V2(t) = - Aie2 - A2e2-----A,,^ + A2e„ - -oN +

< - Aiei - A2e2-----A„e;i + |A2||e„| - Eo|e„| +

< - Aie2 - A2e2-----A„e2 + (-o - -o) |e,| +

1 (40)

- Aie2 - A2e2-----A„e, + -

<|"| + t)

- o-E o

where EE0 = E0 - EE0 is the estimation error and bounded; hi, h2, h3, and h4 are positive learning rates. Differentiating (40) with respect to time and using (39), we can obtain

V2 (t) = eiei + e2e2 + ••• + enen

aT a aTa bT b EE 0EE 0 + i-i- + — + — + -0-0

hi h2 h3 h4

= ei(e2 - liei) + e2(e3 - l2e2 - ei)

+ ••• + en(MT 7 + MT (7a~ + 7bb) + A2 - en-i(t) - Anen + ucc)

aT a aTa bT b EE 0EE 0 + i—i- + - + - +

hi h2 h3 h4

For achieving V2(t) < 0, the error estimation law is designed as follows

EE0 = - EE0 = h41 en | (46)

Then equation (44) can be rewritten as

V2(t)= - Aiei - A2e2-----Anen < 0 (47)

As a result, similar to the discussion in the above section, the asymptotical stability of the proposed control system can be guaranteed in the sense of the Lyapunov theory.


In this section, the proposed control algorithms are applied to validate the feasibility and effectiveness, and the reference commands with different algorithms are applied to the HIL simulation platform for the ELS as shown in Figure 3, which are evaluated by the root mean square error of the output torque. The host control computer with the proposed algorithm to achieve torque tracking control of the ELS mainly analyzes and

~t ( „ b \ _ EE0EE0

+ b en gb 7 + — + en(D2 + ucc) + -

b h3 h4

The adaptive learning laws of the VSFSWC and CC are given in the following so that the convergence of the proposed algorithm is guaranteed

Then taking advantage of equations (39) and (41), we can obtain the following

Figure 3. Components of the ELS.

calculates the feedback torque of the torque sensor, the compensation torque of inertia disks, and the actual torque required. The signal conditioning center is mainly used for processing signals from the sensors and the configuration settings needed for the hardware switches. According to the output of the host control computer, the output of drive provides the appropriate current to drive the torque motor. The torque sensor is used to collect the output of the torque motor, which is fed back to the host control computer to form the closed-loop control for the loading torque. The function of the resolver is to gathering the angle position of the gear reducer, which can be used to calculate the real-time command torque for the torque motor. The function of the angular velocity sensor is to gathering the angular velocity of the gear reducer to provide more actual values for the identification methods and the control strategies. The angular acceleration sensor collects the angular acceleration at the end of the gear reducer to compute the inertia moment produced by the inertia disks. The number of inertia disks depends on the actual demand of the gun control system, which provides various rotational inertia values to simulate the inertia torque. The following simulations are carried out in an Inter core i5 CPU with 3.2 GHz rate, 4 GB RAM, and 64-bit operating system. The parameters of the actual gun system are as follows: the rotary inertia of the turret is 7000 kg m2, the total friction torque Mr is less than 1200 Nm, and the center position of gravity is located in 0.35 m away from the rotation center of the turret, the reduction ratio is 360. Therefore, the relative parameters are selected as gth = 0.65, Ith = 0.12, g0 = 0.1, t = 0.01, ho = 0.22 5, h1 = 8, and h2 = h3 = h4 = 0.5, the center value c of mt is c = [-1.080 -0.225 0 0.225 1.080], the scaling value d of mt is d = [0.50 0.35 0.30 0.35 0.50], the input vector is chosen as [ NB NS Z PS PB ], the rotary inertia of every inertia disk is 2.6 X 10"3kg m2, the control period is 2 ms, and the sample frequency is 10 kHz. Considering the experimental results, the parameters of a, b, and |i are changed with the h1, h2, h3, en(t), and so on, EE0 is used in the CC, the function of the CC is too little to reduce the approximation error when the value of h4 is set more than 1.5. h1 plays a great effect on the coefficients of the xt, the value of h1 is set too small to enhance the influence of xt, which leads to the low convergence rate of the control system. h2 and h3 are related to the dilation and translation parameters, which are set too big or small to achieve the low convergence rate of the control system. The important torque motor parameters are illustrated in Table 1. In order to compare the superiority of all proposed algorithms, the step response reference commands with external disturbances are applied to analyze the rapidity and stability. Next, the static and dynamic loading performance is evaluated by changing the

Table 1. Parameters of the torque motor.

Parameters Values

Rated power pr 2.8 kw

Biggest locked-rotor torque Tmax 19.8 N

Continuous locked-rotor torque Tcont 12 N

Continuous current Icont 80 A

Rated speed vr 1200 r/min

Coefficient of electromagnetic torque Kt 2.6 N m/A

Coefficient of counter electromotive force Ke 9.6 mv/rpm

Equivalent inductance of armature circuit Lm 0.0036 H

Equivalent resistance of armature circuit Rm 1.2 O

Rotational inertiaJm 0.038 kg m2

Viscous friction coefficient Bm 0.22 N m s/rad

fO.i &

10.6 s

10.4 ê 0.2

' , ' . . . .......



if —FSWC and CC


- 1 —VSFSWC and CC -

—Ideal Output

120 150 180 210 240 270 300 Times (ms)

Figure 4. Step response of the proposed algorithms.

motion styles of the load motor and torque motor, respectively. Then, the different connection stiffnesses of the ELS are carried out by the proposed algorithms. In the above experiments, the number of the inertia disks is chosen as 3. Finally, the robustness performance is evaluated by changing the number of inertia disks and the frequencies of the load motor.

Step response with external disturbance

In order to evaluate the proposed algorithms, such as SWC, FSWC, FSWC and CC, VSFSWC, and VSFSWC and CC, the step response with external disturbances is applied to the ELS when the load motor is blocked as shown in Figure 4. When the amplitude of the reference command is selected as 10 N, the VSFSWC and CC shows the best performance among all the proposed algorithms, which reaches the steady state in only 96 ms. Although there always exists torque tracking errors using the SWC, FSWC, and VSFSWC, the VSFSWC with smallest torque errors can meet the control requirements. When the CC is introduced into the FSWC and VSFSWC, the torque tracking errors are decreased effectively. An external disturbance is exerted to the ELS at the moment of 180 ms, the similar conclusion above is obtained, which further verifies the rapidity and stability to some extent.

Figure 5. 1 Hz extraneous torques.

Figure 6. 5 Hz extraneous torques.

Figure 7. 10 Hz extraneous torques.

Static loading

In order to analyze the influences of the load motor, the torque motor is working in the stationary state when the load motor is working in the amplitude of 1 rad and different frequencies, such as 1, 5, and 10 Hz. The control requirement of the load motor is 0.8 mil (1 rad = 954.93 mil). The VSFSWC, FSWC and CC, and VSFSWC and CC algorithms are carried out to show the existing extraneous torques in Figures 5-7. The results indicate that the output of extraneous torque is increasing when the frequencies become bigger, especially the bigger output are happening at the moment of changing the rotary direction of the load motor. The maximum and mean absolute values of the output torque are described in

Figure 8. Sinusoidal tracking with ydl(t) = 5 sin 20pt: (a) VSFSWC and CC means the actual output with the VSFSWC and CC algorithm and (b) the number of neuron variation with the VSFSWC and CC algorithm.

detail in Table 2, which further show that the VSFSWC and CC algorithm has great superiority than the VSFSWC and FSWC and CC. The mean tracking errors of the load motor meet the gun control system requirements, which have nothing to do with the torque motor and the control algorithms when the torque motor works in the stationary state, although the load motor works in different frequencies.

Dynamic loading

To validate the torque tracking performance with the load motor's position disturbance, the torque motor is working in different amplitudes and frequencies when the load motor is working in the motion Adit = sin(0.5pt) rad. Considering the practical working conditions, the reference commands of the torque motor are selected as yd1(t) = 5 sin 20pt, yd2 (t) = 5 sin 10pt, yd3(t) = sin 2pt, yd4(t) = 5 sin 2pt, and yd5(t) = 10 sin 2pt, the experimental results including the torque tracking and variations of the number of neurons are presented in Figures 8-12 and Table 3, respectively.

In Table 3, the mean tracking errors of the load motor meet system control requirements. However, the smaller the reference command amplitudes of the torque motor and the higher the reference command

Table 2. Values of the output extraneous torque and track error.

Frequencies of the load motor (Hz) VSFSWC (N m) Max Mean FSWC and CC (N m) Max Mean VSFSWC and CC (N m) Max Mean Mean tracking error of the load motor (mil)

1 0.5237 0.2315 0.3108 0.2518 0.0827 0.0371 0.2518

5 0.8247 0.3615 0.4982 0.2385 0.2327 0.0954 0.2385

10 1.3830 0.9672 0.5723 0.2741 0.2581 0.1126 0.2741

FSWC: fuzzy self-recurrent WNN controller, CC: complementary controller.

Times /s

Figure 9. Sinusoidal tracking with yd2(t) = 5 sin 10pt: (a) VSFSWC and CC means the actual output with the VSFSWC and CC algorithm and (b) the number of neuron variation with the VSFSWC and CC algorithm.

frequencies of the torque motor (within certain limits), the worse the control performance of the torque motor and the bigger the mean tracking errors of the load motor. In addition, the figures about the number of neurons further show that the better the control performance, the less the neuron number. However, what we need to remind is that the control performance is relatively worse when the direction of the torque motor is changed, which plays a negative effect on the overall performance and need to be further studied in the future.

Different connection stiffnesses

In order to further evaluate the proposed algorithms, we change the connection stiffnesses of the ELS system by different couplings and torque sensor. The different

Figure 10. Sinusoidal tracking with yd3(t) = sin 2pt: (a) VSFSWC and CC means the actual output with the VSFSWC and CC algorithm and (b) the number of neuron variation with the VSFSWC and CC algorithm.

connection stiffnesses play effects on the ELS in Figure 13; Gf is the connection stiffness coefficient.

From Figure 13, the frequency characteristic curve is shown when the connection stiffness is increased by five multiples. The system shear frequency is increasing and the phase characteristic tends to be steep when the connection stiffness is improved, which also broadens the system bandwidth. However, the bigger resonance peak happens when the system bandwidth is overextended, which leads to the system instability. Therefore, aiming at different connection stiffness, the proposed VSFSWC and CC is applied to the ELS when the torque motor is working in the motion 5 sin 2pt. The comparison of different connection stiffnesses is given in Table 4, which further verifies the effectiveness of the proposed controller.

Table 3. Comparison of different reference commands and mean tracking errors.

Reference command (N m)

Mean tracking error of the load -motor (mil)

Torque motor


error (%)

Phase error (0)

Number of neurons


Stable status

5 sin 20pt 5 sin I0pt sin 2pt 5 sin 2pt 10 sin 2pt

0.5205 0.4592 0.3382 0.2762 0.2149

6.43 5.81 4.12 3.38 2.26

-8.9329 -5.6582 -3.8128 -2.0792 -1.0381

24 24 2I I8 I5

I2 I2 I2 I2 9

Table 4. Comparison of different connection stiffnesses.

Connection stiffness coefficient Amplitude error (%) Phase error (0) Number of neurons

Maximum Stable status

10 N m/rad 9.43 -9.1538 45 24

50 N m/rad 6.25 -5.8925 30 18

250 N m/rad 3.38 -2.0792 18 12

1250 N m/rad 3.89 -1.6267 24 15

Times /s

Figure 11. Sinusoidal tracking with yd4(t) = 5 sin 2pt: (a) VSFSWC and CC means the actual output with the VSFSWC and CC algorithm and (b) the number of neuron variation with the VSFSWC and CC algorithm.

Figure 12. Sinusoidal tracking with yd5(t) = 10 sin 2pt: (a) VSFSWC and CC means the actual output with the VSFSWC and CC algorithm and (b) the number of neuron variation with the VSFSWC and CC algorithm.

Robustness performance

The robustness performance is evaluated by changing the number of inertia disks and the frequencies of the

load motor with the proposed VSFSWC and CC algorithms when the torque motor is kept working in the stationary state. From Table 5, it is clear that the

Table 5. Comparison of robustness performance.

Number of inertia disks

Frequencies (Hz)

Maximum extraneous torque (N m)

Mean extraneous torque (N m)

0.0912 0.0827 0.0932 0.0964

0.2617 0.2327 0.2831 0.2781

0.3040 0.2581 0.3106 0.3085

0.0409 0.0371 0.0391 0.0412

0.1176 0.0954 0.1072 0.1204

0.1338 0.1126 0.1286 0.1385

number of inertia disks has no specific relationship to do with the control performance, but with the same parameters of the above controller, there inertia disks show better control performance. When the frequencies of the load motor are becoming bigger, the control performance is relatively worse. However, all the experimental results meet the gun control system requirements, which further show that the proposed controller takes advantages of the robustness and application performance.


In this article, the VSFSWC and CC is proposed to address the complexities in the ELS, which is composed of two terms, one is used to construct equivalent nonli-nearities and uncertainties, the other is designed to eliminate the effect of the approximation error. Moreover, in combination with ideal backstepping controller, the adaptive learning laws are introduced to guarantee the control system stability in the sense of the Lyapunov theory. Finally, extensive comparative simulations and experiments prove that the ability of suppressing the negative influence of the extraneous torque and improving the dynamic torque tracking performance is significantly promoted. In future works, the authors intend to further study the effect of changing the direction of the torque motor.

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.


The author(s) received no financial support for the research, authorship, and/or publication of this article.


1. Wang L, Qian LF, Gao Q, et al. Torque load simulator control research in servo system based on fuzzy PID with gray prediction. Acta Armament 2011; 33: 101-104.

2. Wang CW, Jiao ZX, Wu S, et al. Nonlinear adaptive torque control of electro-hydraulic load system with external active motion disturbance. Mechatronics 2014; 24: 32-40.

3. Dao TL, Dinh QT and Kyoung KA. A torque estimator using online tuning grey fuzzy PID for applications to torque-sensorless control of DC motors. Mechatronics 2015; 26:45-63.

4. Yao JY, Jiao ZX and Han SS. Friction compensation for low velocity control of hydraulic flight motion simulator: a simple adaptive robust approach. Chinese J Aeronaut 2013; 26: 814-822.

5. Yao JY, Jiao ZX and Yao B. Robust control for static loading of electro-hydraulic load simulator with friction compensation. Chinese J Aeronaut 2012; 25: 954-962.

6. Yao JY, Jiao ZX, Shang YX, et al. Adaptive nonlinear optimal compensation control for electro-hydraulic load simulator. Chinese J Aeronaut 2010; 23: 720-733.

7. Wang CW, Jiao ZX and Long Q. Adaptive velocity synchronization compound control of electro-hydraulic load simulator. Aerosp Sci Technol 2015; 42: 309-321.

8. Wang CW, Jiao ZX, Wu S, et al. An experimental study of the dual-loop control of electro-hydraulic load simulator (EHLS). Chinese J Aeronaut 2013; 26: 1586-1595.

9. Man Z, Wu HR and Palaniswami M. An adaptive tracking controller using neural networks for a class of nonlinear systems. IEEE T Neural Networ 1998; 9: 947-955.

10. Calise AJ, Hovakimyan N and Idan M. Adaptive output feedback control of nonlinear systems using neural networks. Automatica 2001; 37: 1201-1211.

11. Mojtaba A and Morteza T. Full-adaptive THEN-part equipped fuzzy wavelet neural controller design of FACTS devices to suppress inter-area oscillations. Neurocomputing 2013; 118: 157-170.

12. Soheil G and Morteza T. Single-hidden-layer fuzzy recurrent wavelet neural network: applications to function approximation and system identification. Inform Sciences 2015; 294: 269-285.

13. Yoo SJ, Park JB and Choi YH. Stable predictive control of chaotic systems using self-recurrent wavelet neural network. Int J Control Autom 2005; 3: 43-55.

14. Lin FJ, Chen SY and Huang MS. Adaptive complementary sliding-mode control for thrust active magnetic bearing system. Control Eng Pract 2011; 19: 711-722.

15. Zekri M, Sadri S and Sheikholeslam F. Adaptive fuzzy wavelet network control design for nonlinear systems. Fuzzy Set Syst 2008; 159: 2668-2695.

16. Hsu CF. Intelligent control of chaotic systems via self-organizing Hermite-polynomial-based neural network. Neurocomputing 2014; 123: 197-206.

17. Hsu CF. Self-organizing adaptive fuzzy neural control for a class of nonlinear systems. IEEE T Neural Networ 2007; 18: 1232-1241.

18. Lan YH, Gu HB, Chen CX, et al. An indirect Lyapunov approach to the observer-based robust control for fractional-order complex dynamic networks. Neurocomputing 2014; 136: 235-242.

19. El-Sousy FFM. Adaptive hybrid control system using a recurrent RBFN-based self-evolving fuzzy-neural-network for PMSM servo drives. Appl Soft Comput 2014; 21: 509-532.

20. Wang YN, Mai TL and Mao JX. Adaptive motion/force control strategy for non-holonomic mobile manipulator

robot using recurrent fuzzy wavelet neural networks. Eng Appl Artif Intel 2014; 34: 137-153.

21. Ali MS. Stability of Markovian jumping recurrent neural networks with discrete and distributed time-varying delays. Neurocomputing 2015; 149: 1280-1285.

22. Lee W, Lee CY, Jeong YB, et al. Friction compensation controller for load varying machine tool feed drive. Int J Mach Tool Manu 2015; 96: 47-54.

23. Gao BW, Shao JP and Yang XD. A compound control strategy combining velocity compensation with ADRC of electro-hydraulic position servo control system. ISA T 2014; 53: 1910-1918.

24. Chen CH. Intelligent transportation control system design using wavelet neural network and PID-type learning algorithms. Expert Syst Appl 2011; 38: 6926-6939.

25. Maryam SK, Farid S and Maryam Z. Design of adaptive fuzzy wavelet neural sliding mode controller for uncertain nonlinear systems. ISA T2013; 52: 342-350.

26. Slotine JJE and Li WP. Applied nonlinear control. Engle-wood Cliffs, NJ: Prentice Hall, 1991.

27. Zhang Q. Using wavelet network in nonparametric estimation. IEEE T Neural Networ 1998; 8: 227-236.

28. Lin CM and Li HY. Self-organizing adaptive wavelet CMAC backstepping control system design for nonlinear chaotic systems. Nonlinear Anal: Real 2013; 14: 206-223.

29. Beliczynski B. A method of multivariable Hermite basis function approximation. Neurocomputing 2012; 96: 12-18.

30. Cheng KH. Self-structuring fuzzy-neural backstepping control with a B-spline-based compensator. Neurocomputing 2013; 117: 138-149.

31. Hsu CF, Lee TT and Tanaka K. Intelligent nonsingular terminal sliding-mode control via perturbed fuzzy neural network. Eng Appl Artif Intel 2015; 45: 339-349.