Hindawi Publishing Corporation Journal of Applied Mathematics Volume 2014, Article ID 272391, 18 pages http://dx.doi.org/10.1155/2014/272391

Research Article

Nonlinear Control of an Active Magnetic Bearing System Achieved Using a Fuzzy Control with Radial Basis Function Neural Network

Seng-Chi Chen,1 Van-Sum Nguyen,1 Dinh-Kha Le,1 and Nguyen Thi Hoai Nam2

1 Department of Electrical Engineering, Da-Yeh University, Changhua 51591, Taiwan

2 Department of Electrical Engineering, Hue Industrial College, Hue 47000, Vietnam

Correspondence should be addressed to Seng-Chi Chen; amtf.csg@msa.hinet.net

Received 14 May 2014; Revised 27 August 2014; Accepted 10 September 2014; Published 12 November 2014 Academic Editor: Xinkai Chen

Copyright © 2014 Seng-Chi Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Studies on active magnetic bearing (AMB) systems are increasing in popularity and practical applications. Magnetic bearings cause less noise, friction, and vibration than the conventional mechanical bearings; however, the control of AMB systems requires further investigation. The magnetic force has a highly nonlinear relation to the control current and the air gap. This paper proposes an intelligent control method for positioning an AMB system that uses a neural fuzzy controller (NFC). The mathematical model of an AMB system comprises identification followed by collection of information from this system. A fuzzy logic controller (FLC), the parameters of which are adjusted using a radial basis function neural network (RBFNN), is applied to the unbalanced vibration in an AMB system. The AMB system exhibited a satisfactory control performance, with low overshoot, and produced improved transient and steady-state responses under various operating conditions. The NFC has been verified on a prototype AMB system. The proposed controller can be feasibly applied to AMB systems exposed to various external disturbances; demonstrating the effectiveness of the NFC with self-learning and self-improving capacities is proven.

1. Introduction

Magnetic bearings are electromechanical devices that use magnetic forces to completely levitate a rotor or suspend it in an air gap without physical contact. Because the system undergoes no friction or wear, it requires no lubrication. In addition, magnetic bearings do not pollute the environment, have a long working life, and can be used in a wide range of applications in aerospace, energy, transportation, and other high technology fields, as well as in high speed ultra-precision machine tools [1-3]. By contrast, conventional mechanical rotary bearings are in physical contact with the shaft, resulting in friction between the shaft and bearings. The temperature of the bearings increases greatly at high rotating speeds, resulting in substantial energy loss and eventual wearing down of the bearings. However, it is difficult to design active controls for magnetic bearing systems because of their high nonlinearity and unstable open-loop electromagnetic dynamics. The suspending structure of a magnetic

bearing differs from that of a conventional bearing. Because no contact occurs between the shaft and the bearing, nearly no friction occurs in the AMB system. Moreover, its stiffness and damping can be changed through rotor displacement and velocity feedback control.

In recent years, many intelligent control techniques, such as fuzzy control, adaptive PID control, neural networks control, adaptive fuzzy control, and other control methods, have been developed and applied to the position control for AMB system [4-9]. Although fuzzy control has been successfully applied in several industrial automations, it is not an easy task to obtain an optimal set of fuzzy membership functions and rules in FLC. RBFNNs are used in a variety of applications such as pattern recognition, nonlinear identification, and control time series prediction. In this study, an AMB that supports a ventilator rotor was investigated to ensure the absence of friction and wear between the rotor and stator and improve the life of the ventilator and its rate of rotation. This paper proposes a method for controlling the position

(6) Magnetic bearing (7) Ventilator

Figure 1: Architecture of a ventilator magnetic bearing system.

of the rotor by using the neural fuzzy controller (NFC) approach. The method employs a fuzzy controller system with radial basis function neural network (RBFNN) rotation to identify the AMB system by Jacobian transformation [1014]. The parameters of fuzzy logic controller (FLC) can be optimally tuned to solve the problem of unbalanced vibration in the AMB system by applying the gradient descent method [15] and the real time values according to the AMB system information. The experimental results of the AMB system indicated a satisfactory control performance without overshoot and with no steady-state errors. This control system also satisfies the requirements of real time response and stability of the control AMB system under disturbances.

The remainder of this paper is organized as follows. Section 2 introduces the principle and mathematical model of the AMB system. Section 3 discusses the dynamics of the AMB system. A motion of equation is formulated whereby characteristics of the rotor for rolling and pitching motions and interference between the horizontal and vertical directions can be analyzed. Section 4 presents the structure of a fuzzy controller based on RBFNN identification of the position control loop for mechanical model of the AMB system. Section 5 presents the stability analysis for the AMB system. Finally, Section 6 presents the results and discussion, and Section 7 concludes the paper.

2. Structure and Mathematical Model of the Active Magnetic Bearing System

Many recent studies on magnetic bearings have focused on AMBs. The AMB system proposed in this paper is presented in Figure 1. The system included a ventilator, a rotor shaft, a magnetic bearing, a coupling device, and a driving motor. The drive system of the AMB system included differential driving mode power amplifiers and an analog to digital (A/D) converter, as indicated in Figure 2. Two sensors were

positioned to measure the displacement of the rotor from the reference position horizontally and vertically. The A/D converter converted the analog signal received from the position sensors into a digital signal. The NFC used this signal as an input, generated a control effort according to measurements, and conducted it using power amplifiers. Four actuating magnets held the rotor at the reference position and transformed the control signals into magnetic forces. As indicated in Figure 2, two pairs of electromagnetic coils were installed perpendicular to the E-shaped stators. These coils produced attractive electromagnetic forces in the perpendicular direction in response to direct currents. All coils installed in the AMBs had the same turns and were symmetrical and uncoupled. For two degrees of freedom, two opposing electromagnets operated in a differential driving mode [16].

As indicated in Figure 2, an electromagnetic force was applied to the rotor along the x- or y-axis to keep the rotor in the center of the AMB. The variable ib is the bias current and ix and iy are control currents along the x- and y-axes, respectively; x1 and yl are the rotor displacements at the magnetic bearing. Following Schweitzer [16-18], the total nonlinear attractive electromagnetic forces along the x- and y-axes are given as follows:

fix fix

fly fly -k\

(ib + iX f (ib - ix)'

{xa -xi) {xa +xi) {ib + iyf {ib - iyf

{yg -7i) {yg + 7i)

Figure 3 presents a three-dimensional plot of the magnetic force, which is nonlinearly related to the current and the rotor displacement.

Position sensor

Power amplifier

Figure 2: Drive system of an AMB.

400 . 300 ,

200 x 100 x 0

-400 0.3

% -0-4-2.5 -2 -1-5 .^U

Co»«oK

Figure 3: Characteristics of electromagnetic force.

A Taylor expansion is performed and linearization is conducted using (1), yielding the following equations that pertain to the x-axis:

kds —

dflX iX1,ix)

dflX iX1,ix)

x1=0,ix=0 xg2

= 4k^,

x1=o,ix=o Xg1

where kis and kds are the displacement and current stiffness parameters of the magnetic bearings, respectively. In this AMB system, the coil on the x- and y-axes circulates the same bias current (ih). Because the nominal air gaps along the x-and y-axes are also the same (xg = yg), the position and current stiffness parameters kis and kds that are obtained from the x-axis are the same as those obtained for the y-axis.

3. Dynamic Model of Active Magnetic Bearing System

Rotor dynamics is a crucial aspect of magnetic bearings, particularly at high speeds. It refers to the results of classical vibration theory and gyromechanics; we can know the phenomena of natural vibration, forward and backward whirl, critical speed, precession, and gyroscopic effect, and the dynamic equation for the horizontal shaft magnetic bearing is studied, and the control schemes are also discussed [1921]. In this section, the dynamic equation of the proposed AMB system is derived. An electric AMB, whose rotor and electromagnet are not in contact, was studied. The shaft was suspended horizontally by magnetic forces at one side while being connected to a driving motor with a flexible coupler at the other side. A coordinate of a fixed frame (x^y^z^) was defined as indicated in Figure 4.

The coordinate Xj is in the coupling horizontal direction, yj is perpendicular to the x^-axis, and the z^-axis coincides with the rotor's axis. Four magnetic attractive forces, f1x, f2x, f1y, and f2y, are exerted on the rotor along the Xj- and yj-axes. The total external forces exerted on the rotor are denoted by fex, fey, and fez, and the variables Fx1 and Fx2 are coupling forces.

Figure 4: Coordinate of a fixed frame.

1 /' ' r >''

1 I I I / /

G (x0,y0,z0)

Figure 5: Coordinates of a rotor.

3.2. The Transformation between Rotor and Fixed Coordinate of Rotor Dynamic. The fixed frames are related to the rotor coordinates, as indicated in Figure 5. The zr-axis is obtained by rotating the Zy-axis by f on the horizontal plane and by d on the vertical plane. The rotating frame with regard to the fixed frame of reference, x0, y0, and z0, shows the location of G with respect to the fixed frame. When the rotor rotates by 0 around the zr-axis, the xr- and yr-axes also rotate. Therefore,

xr = x0 cos <p + y0 sin yr = -x0 sin <p + y0 cos (p.

Figure 6: Rotor in moving frame of reference.

To achieve an effective control result, f and d are assumed to be negligible. Hence,

cos 0 sin 0 0

- sin 0 cos 0 0

cos 0 sin 0 0

- sin 0 cos 0 0

where u, v, and w denote the components of the rotor velocity along the Xj-, y^-, and Zy-axes, respectively.

3.2. Equation Motion of Radial Force for AMB System. The rotor was assumed to be a rigid body and (xj, yj, Zj) was assumed to be fixed to the end of the shaft and rotate with it, as indicated in Figure 6. Therefore, the angular velocity of (xr, yr, zr) is equal to that of the rotor, and it follows that [22]

°Xyz = pi +lj +rk'

Ho = Jxpi +Jyqj +Jzrk'

where co xyz are the angular velocities of the coordinate system

(xr, yr, zr), H0 is the angular momentum around the origin O, and Jx, Jy, Jz are the moments of inertia around the xr-, yr-, and zr-axes, respectively. Calculating (8) with respect to time yields

For a convenient expression, (3) can be rearrang matrix form as follows:

'xr cos 0 sin 0 V

yr. - sin 0 cos 0 Jo .

in the

Ho = (H,

l0 ¡ + wxyzHo ' xyz

= [Jxp - (J y - Jz) qr] 1 + [Jxq + (J y - Jz) pr] 1 (9)

+ [lzr+{jy -Jx)]1-

Assuming that the rotor is a rigid body and symmetrical with respect to the zr-axis and that Jx = Jy, the following equations of motion are obtained:

Deriving (6) leads to

Figure 6 presents a coordinate system of moving frame Gxryrzr from fixed frame (O) to the magnetic bearing. The angular velocity components of the rotor along the xr-, yr-, and zr-axes are denoted by p, q, and r, respectively. The moments of external force along xr-, yr-, and zr-axes are denoted by M, N, and L, respectively. In physical terms, p, q, and r represent the angular velocities of pitching, yawing, and rolling, respectively (Figure 6).

The rotor, suspended by four magnetic forces, is attached to a shaft that is connected to a fixed driving motor with a coupler on the other side (Figure 4). Thus the external torques exerted on the rotor can be described as

*r-axis: {fix - fix - ao) l>

^r-axis: (fiy -f2y -Mg-po)l, (13)

zr-axis: Tm - pr - T0,

where ao and po are the damping and stiffness force exerted by the coupler, Tm is the motor torque, p is the damping torque coefficient, and T0 is the coulomb friction torque. Assuming that the rotor has a small displacement xr and yr in the horizontal and vertical planes, respectively, ao and po are given by

ao = i xo + i Xo

kL cl2 . ßo = - yo + - yo,

where c and k are displacement coefficients and l1 and l2 are the distances between the shaft center and the flexible coupling and the magnetic bearing, respectively. By using (4), the forces exerted on the rotor and the moments around the original point are represented by the following matrix forms in (15) and (16), respectively. Consider

cos 0 sin 0 - sin 0 cos 0

cos 0 sin 0 - sin 0 cos 0

fix fix ao fly - fly -Mg-ßo

■(fly -fly -Mg-ßo )l (fix - fix - ao) 1

P = <P

)qr = M, (10)

Equation (17)

pr = N, (11)

Jzr = L. (12) .f.

cos <p sin <p - sin <p cos <p

cos 0 - sin 0 sin <p cos <p

By using (10) and (11), (18) becomes

p - r ' 0 1 d'

A. -1 0 .Y.

. (17)

9' 1 cos 0 - sin 0 M

.f. = Tx sin 0 cos 0 N

0 1 -1 0

cos 0 - sin 0 sin 0 cos 0

Using (14), by letting x1 = If and y1 = -Id and eliminating p and q in (19) and assuming f and d to be small, (19) becomes a linearizing equation of motion as follows:

" -clli -rJzl

d' Jx Jx

.f. rJz -clli

[ Jx Jx ]

d' klli 1 0' '9'

y. ~TX 0 1 .Y.

-fly + fly + M0 fix fix

From (20), the block diagram can be drawn, as indicated in Figure 7. In this figure, the motion along the x^-axis and the y^-axis is independent if r = 0. The horizontal motion is coupled with the vertical motion if r = 0. This is the so-called gyroscopic effect. Furthermore, the coupling effect increases with r. To decrease the interaction between the x- and y-axes, the Jz/Jx ratio shouldbe decreased [23].

4. Neural Fuzzy Controller Design for Position Control

4.1. Fuzzy Logic Controller (FLC). Recently, the FLC has emerged as an effective tool for stabilizing a nonlinear system, such as an AMB system, a magnetic levitation system, or other electronic devices [24-30]. The FLC is a methodical approach for controlling a nonlinear system and is a heuristic technique for enhancing the operation of a closed loop system. Although FLC performance is a function of its ability to simulate many functions simultaneously, its output results are considerably thorough. The FLC presented in Figure 8 is composed of four main parts: one performs fuzzification; one provides the rule base; one is an inference engine; one performs defuzzification.

The structure of a fuzzy controller, based on RBFNN identification of the AMB system, is presented in Figure 9. It consists of a fuzzy controller, a reference model, an RBFNN,

Figure 7: Diagram radial of pitching, yawing, and rolling of rotor.

Crisp input

Knowledge base

Fuzzification Defuzzification

Fuzzy inference engine /

Fuzzy input Fuzzy output

Crisp output

Figure 8: Fuzzy logic controller block.

and an adjusting mechanism. The operating principle of this controller is based on AMB dynamics and control knowledge can be incorporated into an NFC design based on RBFNN identification [14, 31, 32]. A detailed description of these is shown in Figure 9, where x* and xrm are a step reference and a reference model, respectively. xrbj and f are the outputs of neural and fuzzy neural controller, respectively. xp is the response of displacement; Z-1 is the back shift operator. Uj represents the output of the fuzzy controller. FPID is a fuzzy proportional integral derivative control for the current loop [33]. Fd is the external torque from ventilator. FI and DFI are the fuzzification converted controller inputs (e and de), respectively, and the defuzzification based on product inference rule center average defuzzification the conclusions of the inference mechanism into actual inputs.

Based on Figure 9, the tracking error e and the error change de are defined by the following equations:

e (k) = xrm (k) - xp (k) de(k) = e(k)-e(k-l).

The design procedure of the fuzzy controller algorithm is as follows. First, e and de are taken as the input variables

of the fuzzy controller, and their linguistic variables are defined as E and dE. The linguist values of E and dE are {Ao, Ai,A2, A3, A4,A5,A6] and {B0,Bi,B2,B3,B4,B5,B6}, respectively. Each linguistic value of E and dE is based on the symmetrical triangular membership function, which is presented in Figure 10. Second, the membership degrees of e and de are computed. Figure 10 indicates that only two linguistic values are excited in any input value, and it is the membership degree. The variables (e) and (de) are easily derived using the flowchart presented in the figure, and the membership degree is obtained as follows [34]:

^ (e) = ,

' ei+1 - ei (e) = 1- (e).

Similar results are obtained on computing the membership degree (de). Third, the initial fuzzy controller rules are selected by referring to the dynamic response characteristics, such as

IF e is A; and de is Bj THEN uf is cu, (23)

where i and j range from 0 to 6, A; and Bj are fuzzy numbers, and Cjj is the real number. The graph of fuzzification and the fuzzy rule table is presented in Figure 10. Finally, the fuzzy system Uj(e, de) is constructed using the singleton fuzzifier, product inference rule, and the central average defuzzifier method is presented. Although 49 fuzzy rules which are listed in Figure 10 will be inferred, only four fuzzy rules can be effectively excited to generate a nonzero output. Therefore, if an error e is located between et and ei+1 and an error change de is located between dej and de^+1, only four linguistic values AAi+1, Bj, and Bj+1 and the corresponding consequent

NFC for position control loop

Figure 9: Block diagram neural fuzzy control for AMB system.

Input MFs of e (i = 1)

(e) = 1 - Ma1 (e)

MA1 (e) <-

e de \ Ao Ai A2 A3 A4 A5 A5

Bo ^00 % ^02 ^03 ^04 C05 ^06

Bi c10 cii c12 c13 c14 c15 c16

B2 c20 c21 c22 c23 c24 c25 c26

B3 c30 c31 c32 c33 c34 c35 c36

B4 c40 c41 c42 c43 c44 c45 c46

B5 c50 c51 c52 c53 c54 c55 c56

B6 ^50 c61 c62 c63 c64 c65 ^66

Fuzzy rules table of e and de

Figure 10: The membership functions of e, de, and fuzzy rule table.

values Cj,;, Cj+1,;, c^^+1, and cJ+1> i+1 can be excited, and the output of the fuzzy controller can be calculated using the following equation:

. , . T^i cm,n [va„ (e) * vB„ (de)\ Uf (e, de) =-———-

Tn= i !Lj vA„ (e) * ^Bm (de)

i+1 j+1

= Cm,n * dn m,

n=i m=j

where dnm = (e) x (de). And those cmn are adjustable parameters for fuzzy controller. In addition, by using (22), the value of I'm1! dn,m =1 in (24) can be easily derived. The

Ujn is formulated by the output of fuzzy controller (uj) and the output of the integral (I) controller (ujni) as follows:

ufn (k) = ufni (k-l) + (Kpw + K,w) * uf (k)

Aufn (k) = Ufn (k) - Ufn (k-1).

4.2. Radial Basic Function Neural Network

4.2.1. Radial Basis Function Neural Network Architecture. The RBFNN adopted in this part of the study was a three-layer structure, as shown in Figure 11 and comprised an input layer,

a hidden layer, and an output layer. Each layer is explained in detail below.

Input Layer. the RBFNN has three inputs by (k), xp(k), xp(k - 1) and its vector form is represented by

X=[ufn (k),xp (k),xp (k-1)]T. (26)

Hidden Layer. The hidden layer is composed of n RBFNNs that are connected directly to all the elements in the output layer. A node in the hidden layer produces a greater output when the input pattern presented is closer to its center. The multivariate Gaussian function is used as the activated function in the hidden layer of the RBFNN, which is given by the following equations:

hj {X) = e-{imx-^ih>,

j = 1,2,3,

where bj denotes the node center and node variance of the jth neuron, m is the number of neurons in the hidden layer, - Cj|| is the norm value which is measured by the inputs

and the node center at each neuron, and Cj = [cj1 ,Cj2,..., Cjn]T is the center of the jth RBFNN unit.

Output Layer. The network output xrbj is formed by a linearly weighted sum of the number of basic functions in the hidden layer. The network output in Figure 11 is derived using the following equation:

where Wj and hj are the weights from the jth hidden layer neuron to the output layer neuron and the output of the jth node in the hidden layer, respectively.

4.2.2. Training Algorithm in Radial Basis Function Neural Network. To describe the online learning algorithm of the proposed NFC, the energy function E is defined as follows:

v 1 i \2 1 2

Based on the gradient descent method, the learning algorithm of the weights, node center, and variance are adjusted using the following equations:

Wj (k + 1)

= W: (k) + qenn (k) h: (k)

Cji (k+1)

X, (k) - Cj, (k)

= Cji (k) + qenn (k) Wj (k) hj (k)--

bj (k+1)

Ita (k)-cn (Ml

= b: (k) + nenn (k)W: (k)h: (k) 1 ' " "

b3j (k)

where j = 1,2, ...,m, i = 1,2,3, and q is a learning rate. Further, dxp/dujn is Jacobian transformation and can

be derived from Figure 11 and (27) and (28). Consider the following:

dxp dxrbf ™ cn -ufn (k)

4 =^'—¡¡j— (31)

where Kpw and Kiw are the PI controller gains. Therefore, substituting (38) and (39) into (37) and combining with (36), the parameter cm n of the fuzzy controller described in (24) can be adjusted using m = j,j+l, and n = i,i+l,as indicated by the following equation:

4.3. Reference Model. A second order system with natural frequency wn and damping ratio q is used as a reference model in the adaptive control system, as explained in the following equation:

Xrm (s)

x* (s) s2 + 2qun +

By applying the bilinear transformation, the second order system (32) can be transformed into a discrete model using

a^ + aiZ + a2z

w* (z-1) \+b1z-1 +b2z-2' Furthermore, the difference equation is written as follows: wr (k) = -b1wr (k- 1) - b2wr (k-2) + aQw** (k)

+ a1w** (k - 1) + anw** (k -2).

4.4. Fuzzy Logic Control Parameters Adjusting Mechanism. The gradient descent method is used to derive the FLC control law in Figure 9. The adjusting parameters of the fuzzy controller minimize the square error between the rotor displacement and the output of the reference model. The instantaneous cost function is defined as follows:

P _ 1 2_ _ \

= = ~ (Xrm Xp)

e 2' 2

The parameters of cm n are adjusted according to

Acm,n = >

where a represents the adaptive rate of the system. The chain rule is used, and the partial differential equation for Ee in (35) is written as follows:

dEe dxp duf

duf dcm,n

From (24) and using the Jacobian formulation from (31), the following equation is derived:

duf (k)

dcm,n (k)

It is difficult to calculate dxp/duj because of unknown plant dynamics. To overcome this problem and to increase the online learning rate of the connective weights, a delta adaptation law is proposed as follows:

dxp~(r + KT \dXrhf (Kpw + Kiw)

(IT ^ \V » °)1 -Ufn (k) = (Kpw + K,w) ^Wjhj--2-,

kcm,n (k)

™ Cji -ufn (k) (40)

= ae (k) (Kpw + Kiw) dn^mLwjhj-U-•

5. Stability Analysis for Active Magnetic Bearing System

The stability of the controller of an AMB system is critical to its functioning, particularly during start-up. After the operation of the controller starts, the rotor remains on touchdown surfaces. While the rotor stays at the position, the controller can never estimate the motion of the rotor inside the stator because the electromotive force is proportional to the rotational speed. Hence, movement toward its operating position, the center of the two pair of electromagnets, is necessary for the convergence of controller target. To ensure stable functioning of the system, a stability analysis of the NFC was conducted [35-37]. When the parameters of an NFC move toward infinity, the NFC becomes an unstable controller of the AMB system; when the center of one of the membership functions (MFs) approaches infinity, the NFC is stable. Because the firing of the rule which contains these MFs is equal to zero, other rules can identify the system without instability and when the standard deviations of MFs move toward infinity the values of the MFs are equal to one, and the output is finite; therefore, the areas of MFs that are optimal for the NFC to control an AMB system must be used. The RBFNN adopted in this study comprises three inputs, as indicated in (26). The output is represented by (24), and the objective function is defined by (35). The consequent weights are as follows:

W = [^ W2 W3,...,Wn]

mx(n+1) '

Two types of parameter are defined as the antecedent parameters. The first type is the means of the MFs as follows:

f=[f1 f2 Fj,•••,f^].

The second type is the standard deviations of the MFs as follows:

s=№ sj sj,...,sj].

From (29), the discrete Lyapunov function is defined using (44) as follows:

V(k) = E (k) = -e2 (k) = - (xp (k) - xrbf (k)) . (44)

The change of the Lyapunov function at each iteration of (44) is obtained using

AV(k) = V(k+1)-V(k) = - [e2 (k + 1) - e2 (k)]

= 2[e(k+l)-e(k)][e(k+l)+e(k)] (45)

= Ae (k)

Ae (k) + e (k)

Equations (41) to (45) were applied to achieve an approximate error modification for the NFC as shown in the following equation:

Ae (k) = d(k) + 9(k)+p (k)

9(k) =

(de(k) (dS(k)

AS(k),

T ((UrnT A<™)-

where the tr(-) is the trace of matrices. In this controller, only train the consequent parameters so d(k) = 0; 9(k) = 0. Therefore, substituting (46) and (47) into (45) obtains

AV(k)-p(k)

+ e(k)

- tr((W)T a™

1 'r ((£%T

The updating rule for conclusion parameters is

de (k) _ de (k) dxrbf (k) _ dxrbf (k) = X ------- = - dW(k) '

dW(k) dxrhf (k) dW(k) dxrbf (k)

AW (k) = tje (k)

dW(k) '

where q is the learning rate that is used to adjust consequent parameters. From (48) and (49) can be rewriting as follows:

AV(k) = -qe2 (k)

dxrhf (k)

1 2 2 ns + e (k)

dxrhf (k)

where ||-||F is the Frobenius norm, assuming that dxrhf (k)

^w (k) =

dW(k) '

¥wm„ = max ||¥w (k)\\P, k

Yw =\n(pw (k)l\Pf (2-n(pw (k)l\p)2).

According to the standard Lyapunov theory, the change of the Lyapunov function must be less than zero to ensure stability. This can be achieved using (50) and (51) as follows:

AV (k) < 0,

AV(k) = -Ywe2 (k)=^Yw >0,

0 < q <

Equation (52) is an adaptive constraint and the learning rate stability condition changes at each iteration. Therefore, this constraint can be used easily for online training. By calculating (52), the following equation can be derived:

0 < q <

with = max,

* \\Vw(k)\\P.

Equation (53) is a conservative constraint and cannot be calculated through online identification. The equation can also be written using the chain rule as follows:

= 'r ((Ulf AH

= ne(k)(\Vw (k)l

The learning rate of the consequent part is selected to satisfy (54) so that the identifier is stable at each learning cycle. The adaptive controller functions online during AMB system operation. The results of the AMB system are presented in the next section.

6. Results and Discussions

The experimental setup of this study is presented in Figure 12. The laboratory setup included a horizontal shaft magnetic bearing that was symmetrical and was controlled by two axes.

The system was driven by an induction motor with a flexible coupling to isolate the motor vibration. The magnetic bearing included four identical electromagnets that were equally spaced radially around a rotor composed of laminated stainless steel, as indicated in Figure 13. Each electromagnet included a coil and a laminated core composed of silicon steel. The rotor displacement along the vertical y- and horizontal x-axes of the geometric center of the shaft was measured using a pair of eddy current sensors, as indicated in Figure 14. The conversion time of the 16-bit A/D converter was 10 ^s.

Electromagnet

Figure 12: The experimental setup of the AMB system.

Figure 13: Inside view (a) and layout (b) of magnetic bearing.

Figure 14: Eddy sensor and driver power.

Figure 15: Photograph of the experimental setup.

Figure 16: Photograph of current amplifier and current sensing.

[swre]

From 9

FPID controller

M | /W\ -|->< KU

Derivative

FLC bottom

Analog input

Analog input AI0 Advantech PCI-1716 [auto]

T Slider Goto

gain 1.465

Convert PID regulator

Bottom coil

Analog output

Analog output Advantech PCI-1716 [auto]

Figure 17: Model of real time window target for current-control loop.

3 2.5 2 1.5 1

Right electromagnetic coil of magnetic bearing

3 2.5 2 1.5 1

Left electromagnetic coil of magnetic bearing

- u U

Time (s)

46 Time (s)

- Reference

---Current response

- Reference

---Current response

Top electromagnetic coil of magnetic bearing

Bottom electromagnetic coil of magnetic bearing

2.5 2 1.5 1

l^üflfl® r

u u u u

r r

R r

.-L u J Ü

Time (s)

Time (s)

- Reference

---Current response

- Reference

---Current response

(c) (d)

Figure 18: Current loop control of magnetic bearing with and y-axes.

The sampling rate of the 16-bit converter was 100 kHz and the control cycle was approximately 0.1ms. A photograph of the experimental setup is presented in Figure 15. A current amplifier and current sensing circuit design for AMB system are indicated in Figure 16.

The AMB system and a fuzzy basis neural network were implemented using Matlab software; the parameters are listed in Table 1. The requisite interface was a PCI-1716 card comprising an A/D part with 16 channels and a digital input and output part with 16 channels. Matlab software was used to code the proposed controllers.

Table 1: Parameters of an AMB system.

1 Mass of shaft (m)

2 Nominal length of air gap

3 Transverse moment of inertia of rotor (/)

4 Polar moment of inertia of rotor (Jz)

5 Displacement stiffness (Kx)

6 Current stiffness (KJ

7 Bias currents to be used ()

2.72 kg 0.5 mm 0.013 kgm2 0.008 kgm2 342478 N/m 171N/A 1A

Figure 19: The control surface after adjusted NFC.

-1........... .......... .......... ...........

-1.5........... .......... .......... ...........

-2 -1-1-1-1-1-1-1-

0 0.5 1 1.5 2 2.5 3 3.5 4 Time (s)

Figure 20: The output of RBFNN.

An experiment for the AMB system has been verified by the current-control loop, using FPID control and current amplifier (Figure 16). The model established for real time window target (RTWT) test for current-control loop using FPID controller is shown in Figure 17 (single electromagnetic coil), and the sample time is 0.1ms. The current responses of four electromagnetic in current-control loop at the same time are presented in Figure 18. From the results in Figure 18 we can see that the current response (four electromagnetic) is very close to the reference signal with the setting time about 0.01s.

The system and controller parameters are listed in Table

The initial parameters in fy, and Cj were set as follows: Wj = [0.5 0.5 0.5]; b} = [40 0.5 0.5]; and Cj = [40 0.5 0.5; 40 0.5 0.5; 40 0.5 0.5]. Thediscretemodel of the reference model was obtained as follows by using bilinear transformation at the sampling frequency of 1 kHz.

The AMB system included two pairs of electromagnets on the x- and y-axes. The two pairs of electromagnets were simultaneously controlled using two NFCs. The experimental results are presented in Figures 19 and 20 and in Table 3. Figure 20 indicates that the practical output of the NFC converged from -1to1.

Figure 21 indicates the rotor displacement of x- and y-axes in the AMB system. The rotor displacement is small

y-axis

K -0.2

Q -0.4 -0.6

1.5 2 2.5 Time (s)

Q -0.4 -0.6

1.5 2 2.5 Time (s)

Figure 21: Rotor displacement of an AMB system with x- and 7-axes.

S 0.4 ie

1 0.3 me

4 0.2 Q

PlfHWl

0 0.5 1 1.5 2 2.5 3 3.5 4 Time (s)

0 0.5 1 1.5 2 2.5 3 3.5 4 Time (s)

Figure 22: Step response of rotor position of AMB system.

x-axis

x-axis

about 0.12 mm. In general, the rotor displacement in the horizontal direction is smaller than the rotor displacement in the vertical direction because of the effects of gravity on the y-axis. The response of rotor position of x- and y-axes in AMB system covers the entire feasible region (Figure 22). Figures 23(a) and 23(b) show that the orbit of the rotor center is using an NFC at rotating speeds from 10000 rpm to 15000 rpm. From rotating speed of 10000 rpm, the rotor displacement is small about 0.1 to 0.15 mm (Figure 23(a)). When the rotor rotates at a high speed (15000 rpm), the rotor displacement increases about 0.19 to 0.22 mm (Figure 23(b)), but it is still in the permitted limits of nominal length of air gap (xg = 0.5 mm). To evaluate the performance and

characteristics of the RBFNN, we look at an AMB system controlled by an NFC. As we can see, it is unstable system before the first 1.2 s. The parameters of a FLC are adjusted using a RBFNN. The whole controller is applied to the unbalanced vibration in an AMB system into effect after 1.2 s. The results further demonstrate that a short rise time implies a short settling time, low overshoot, and a small steady-state error with external disturbance.

7. Conclusions

In this study, an NFC was developed to levitate rotor displacement in a highly unstable AMB system. The proposed

10000 rpm

15000 rpm

s-l £ -0.1

m -0.2

rn -0.3

TT/r« % r

-0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 Left < horizontal error (mm) > right

-0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 Left < horizontal error (mm) > right

Figure 23: Orbits of rotor center x- and y-axes.

Table 2: The controller parameters.

Table 3: Rule table after adjustment of neuron fuzzy controller.

Parameters Values de e

Reference model A 0 A1 A 2 A 3 A 4 A 5 A 6

a0 98.10-6 Bo -6 -6 -4 -4 -2 -2 0

196.10-6 Bi -6 -4 -4 -2 -2 0 2

a-2 98.10-6 B2 -4 -4 -2.084 -2.109 0 2 2

h -1.96 B 3 -4 -2 -2.012 0 1.964 2 4

h 0.96 B4 -2 -2 0 2.113 2.065 4 4

C 1 Bs -2 0 2.074 2.058 4.022 4 6

20 rad/s B6 0 2 2.055 4.027 4.024 6 6

PI controller

Kiw NFC a

5.6 0.012

0.07 0.5

method can also be used to improve the control performance of other nonlinear systems. The results indicated that the NFC responds well to a reference signal. The ABM system substantially reduced overshoot, shortened the adjustment time, accelerated the response, produced a more robust system, and improved its dynamic and static performance. The proposed controller can be feasibly applied to AMB systems with various external disturbances, and the effectiveness of the NFC with self-learning and self-improving capacities is

proven with RBFNN learning algorithm. This controller has been verified by the position control loop on a prototype AMB system.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors would like to thank the National Science Council and Bureau of Energy, Ministry of Economic Affairs, for financial support under Contract NSC 102-2221-E-212-007 and no. 102-D0624.

References

[1] D. Johnson, G. V. Brown, and D. J. Inman, "Adaptive variable bias magnetic bearing control," in Proceedings of the American Control Conference (ACC '98), vol. 4, pp. 2217-2223, June 1998.

[2] H. W. Cho, C. H. Kim, J. M. Lee, and H. S. Han, "Design and characteristic analysis of small scale magnetic levitation and propulsion system for maglev train application," in Proceedings of the International Conference on Electrical Machines and Systems (ICEMS '11), pp. 1-5, August 2011.

[3] J. Kaloust, C. Ham, J. Siehling, E. Jongekryg, and Q. Han, "Nonlinear robust control design for levitation and propulsion of a Maglev system," IEE Proceedings—Control Theory and Applications, vol. 151, no. 4, pp. 460-464, 2004.

[4] J. Fei and S. Hou, "Robust adaptive fuzzy control for three-phase active power filter," in Proceedings of the IEEE 13th Workshop on Control and Modeling for Power Electronics (COMPEL '12), pp. 1-6, IEEE, Kyoto, Japan, July 2012.

[5] H.-C. Chen, "Optimal fuzzy pid controller design of an active magnetic bearing system based on adaptive genetic algorithms," in Proceedings of the International Conference on Machine Learning and Cybernetics (ICMLC '08), vol. 4, pp. 2054-2060, Kunming, China, July 2008.

[6] F. Matsumura and T. Yoshimoto, "System modeling and control design of ahorizontal-shaft magnetic-bearing system," IEEE Transactions on Magnetics, vol. 22, no. 3, pp. 196-203,1986.

[7] S.-Y. Chen and F.-J. Lin, "Robust nonsingular terminal sliding-mode control for nonlinear magnetic bearing system," IEEE Transactions on Control Systems Technology, vol. 19, no. 3, pp. 636-643, 2011.

[8] T. Dimond, P. Allaire, S. Mushi, Z. Lin, and S. Y. Yoon, "Modal tilt/translate control and stability of a rigid rotor with gyroscopics on active magnetic bearings," International Journal of Rotating Machinery, vol. 2012, Article ID 567670, 10 pages, 2012.

[9] S.-C. Chen and P.-C. Tung, "Application of a rule self-regulating fuzzy controller for robotic deburring on unknown contours," Fuzzy Sets and Systems, vol. 110, no. 3, pp. 341-350, 2000.

[10] B. Shi, L. Y. Wang, L. Peng, M. Xin, and Y. X. Huu, "A modified particle swarm optimization and radial basis function neural network hybrid algorithm model and its application," Intelligent Systems, vol. 1, pp. 134-138, 2009.

[11] N. Al-Holou, T. Lahdhiri, D. S. Joo, J. Weaver, and F. Al-Abbas, "Sliding mode neural network inference fuzzy logic control for active suspension systems," IEEE Transactions on Fuzzy Systems, vol. 10, no. 2, pp. 234-246, 2002.

[12] W. Wu, S. Zhong, and G. Zhou, "A study on PID intelligent optimization based on radial basis function neural networks," in Proceedings of the 3rd International Conference on Consumer Electronics, Communications and Networks (CECNet '13), pp. 57-60, Xianning, China, November 2013.

[13] S. C. Chen, N. V. Sum, and G. Chang, "Application of self-tuning fuzzy PID controller on magnetic levitation system," in Proceedings of the 11th Taiwan Power Electronics Conference & Exhibition, National Tsing Hua University, September 2012.

[14] P. S. Bhowmik, S. Pradhan, M. Prakash, and S. Roy, "Investigation of wavelets and radial basis function neural network for incipient fault diagnosis in induction motors," in Proceedings of the International conference on Circuits, Controls and Communications (CCUBE '13), pp. 1-5, Bengaluru, India, 2013.

[15] F. Guely and P. Siarry, "Gradient descent method for optimizing various fuzzy rule bases," in Proceedings of the 2nd IEEE

International Conference on Fuzzy Systems, vol. 2, pp. 1241-1246, San Francisco, Calif, USA, April 1993.

[16] G. Schweitzer, H. Bleuler, and A. Traxler, Active Magnetic Bearings, vdf Hochschulverlag AG, Zurich, Switzerland, 1994.

[17] K. Zdenko and B. Stiepan, Fuzzy Controller Design: Theory and Applications, CRC Press, 2005.

[18] S.-C. Chen, V-S. Nguyen, D.-K. Le, and M.-M. Hsu, "ANFIS controller for an active magnetic bearing system," in Proceedings of the IEEE International Conference on Fuzzy Systems (FUZZ '13), pp. 1-8, Hyderabad, India, July 2013.

[19] F. Matsumura and T. Yoshimoto, "System modeling and control design of a horizontal shaft magnetic bearing system," IEEE Transactions on Magnetics, vol. 22, no. 3, pp. 196-203,1986.

[20] M. Fumio, K. Hidehiko, and A. Yuji, "Fundamental equation for horizontal shaft magnetic bearing and its control system design," Fuzzy Sets and Systems, vol. 101, pp. 123-130, 2007.

[21] K. Nonami, "Vibration and control of flexible rotor supported by magnetic bearings," Magnetic Bearing, pp. 177-186,1989.

[22] H. Baruh, Analytical Dynamics, McGraw-Hill International, 1999.

[23] S. Y. Yoon, Z. Lin, and P. E. Allaire, Control of Surge in Centrifugal Compressors by Active Magnetic Bearings: Theory and Implementation, Advances in Industrial Control, Springer, 2013.

[24] J. L. Meza, V. Santibanez, R. Soto, and M. A. Llama, "Fuzzy self-tuning PID semiglobal regulator for robot manipulators," IEEE Transactions on Industrial Electronics, vol. 59, no. 6, pp. 27092717, 2012.

[25] M. Maeda and S. Murakami, "A self-tuning fuzzy controller," Fuzzy Sets and Systems, vol. 51, no. 1, pp. 29-40,1992.

[26] S. C. Chen, L. Y. Jyh, N. V. Sum, and S. M. Mao, "A novel fuzzy neural network controller for maglev system with controlled-PM electromagnets," Intelligent Technologies and Engineering Systems, vol. 234, pp. 551-561, 2013.

[27] J. Y. Chen and H. L. Ying, "A self-tuning fuzzy controller design," in Proceedings of the IEEE International Conference on Neural Networks, vol. 3, pp. 1358-1362,1995.

[28] Z. Q. Wu, P. Z. Wang, H. H. Teh, and S. S. Song, "A rule self-regulating fuzzy controller," Fuzzy Sets and Systems, vol. 47, no. 1, pp. 13-21, 1992.

[29] S.-Z. He, S.-H. Tan, F.-L. Xu, and P.-Z. Wang, "Fuzzy self-tuning of PID controllers," Fuzzy Sets and Systems, vol. 56, no. 1, pp. 3746, 1993.

[30] S.-Z. Gao, J. Yang, and J.-S. Wang, "D-FNN based modeling and BP neural network decoupling control of PVC stripping process," Mathematical Problems in Engineering, vol. 2014, Article ID 681259, 13 pages, 2014.

[31] F. A. Ruslan, A. M. Samad, Z. M. Zain, and R. Adnan, "Modelling flood prediction using Radial Basis Function Neural Network (RBFNN) and inverse model: a comparative study," in Proceedings of the IEEE International Conference on Control System, Computing and Engineering (ICCSCE '13), pp. 577-581, Mindeb, Mich, USA, December 2013.

[32] S. S. A. Ali, M. Moinuddin, K. Raza, and S. H. Adil, "An adaptive learning rate for RBFNN using time-domain feedback analysis," The Scientific World Journal, vol. 2014, Article ID 850189, 10 pages, 2014.

[33] M. Ali, S. Khan, M. Waleed, and Islamuddin, "Application of an intelligent self-tuning fuzzy PID controller on DC-DC buck converter," International Journal of Advanced Science and Technology, vol. 48, pp. 139-148, 2012.

[34] M. Negnevitsky, Artificial Intelligence—A Guide to Intelligent Systems, 2nd edition, 2005.

[35] R.-E. Precup, M. L. Tomescu, and §. Preitl, "Fuzzy logic control system stability analysis based on Lyapunov's direct method," International Journal of Computers, Communications and Control, vol. 4, no. 4, pp. 415-426, 2009.

[36] M. A. Shoorehdeli, M. Teshnehlab, A. K. Sedigh, and M. A. Khanesar, "Identification using ANFIS with intelligent hybrid stable learning algorithm approaches and stability analysis of training methods," Applied Soft Computing,vol. 9, no. 2, pp. 833850, 2009.

[37] M. L. Tomescu and G. Petrov, "A stability analysis method for nonlinear systems with fuzzy logic controller," in Proceedings of the 8th International Symposium on Symbolic and Numeric Algorithms for Scientific Computing (SYNASC '06), pp. 141-150, IEEE, Timisoara, Romania, September 2006.

Copyright of Journal of Applied Mathematics is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.