Hindawi Publishing Corporation Discrete Dynamics in Nature and Society Volume 2012, Article ID 279498,21 pages doi:10.1155/2012/279498

Research Article

Implementation of a Fuzzy TSK Controller for a Flexible Joint Robot

M. E. Akbari,1 M. A. Badamchizadeh,2 and M. A. Poor3

1 Department of Electrical Engineering, Ahar Branch, Islamic Azad University, Ahar 5166614776, Iran

2 Faculty of Electrical and Computer Engineering, University of Tabriz, Tabriz, Iran

3 Department of Electrical and Computer Engineering, University of Alberta, Edmonton, AB, Canada T6G 2V4

Correspondence should be addressed to M. E. Akbari, m-akbari@iau-ahar.ac.ir

Received 26 July 2012; Revised 3 September 2012; Accepted 14 September 2012

Academic Editor: Vimal Singh

Copyright © 2012 M. E. Akbari 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.

This paper proposes a fuzzy TSK controller to control a rotary flexible joint manipulator. The flexibility of joint is performed by means of a solenoid nonlinear spring, which is connected between actuator output and joint input in a bilateral connection form to transfer the produced torque; also the smooth model of frictions is used for modeling the dynamics of flexible manipulator. The effect of coulomb friction and also gearbox backlashes is decreased by a pulsation signal as an extra voltage that is added to the control voltage of actuator. Actuator dynamics is modeled by consideration of saturation mode of armature current. Experimental results demonstrate that the proposed controller has an ability to control flexible joint manipulator with a good performance.

1. Introduction

Robot systems with flexible joints are generally used in the aerospace industry, military, chemical, integrated circuits, and nuclear energy systems in order to perform commands with high sensitivity and accuracy [1]. It is demonstrated that the majority of robots have nonlinear dynamic and performance. Hence a high-accurate controlling requires special design and controllers [2] . Linear controllers are still used in robotics though in most cases this type of controllers may result in low performance responses and also encounters stability problems. Industrial environments, where there are strong electrical noises, have the worse conditions because electronic circuits and especially measuring systems are affected by power noises. Similar to other electromechanical systems, robots have uncertainties on their parameters model due to tolerance in physical quantities. In practice, model uncertainties

and disturbance decrease the efficiency of designed controller and also limit the area of convergence. For years, efforts have been made to solve this problem and several methods such as robust, adaptive, and intelligent controllers were proposed [3].

Fuzzy controllers have been widely proposed and analyzed in the literature over the past decade. The most useful property of the fuzzy system is the ability to approximate the mapping from input to output with arbitrary accuracy. Especially, the Takagi-Sugeno-Kang fuzzy model [4] is an adequate fuzzy model for complex systems. This work proposes the design of a TSK-fuzzy controller, which assigns different linear control laws in order to control a flexible joint robot (FJR). These robots have high nonlinearity on dynamics and must be limited in the vicinity of equilibrium point. Such systems may be unstable in out of operating space because of high nonlinearity error. In these cases, even if designing a controller is being realizable, usually closed loop system has no good performance.

The proposed method has some advantages as follows:

(1) improved performance controller for wide changing of payload,

(2) no need for frequent calculations,

(3) decreasing in controller time delay by eliminating some of the low effect sentences,

(4) with adding a vibration signal we could reduce effects of fraction and backlash.

This paper subtends to four sections. In Section 2, we show how to make an acceptable model of Laboratory FJR (L-FJR), which has rigid arm, flexible joint, DC motor, and gearbox as actuator. Also in this section, a new compensator is proposed to decrease friction and backlash effects in FJR modeling. In Section 3, a controller designing based on fuzzy TSK is discussed. In Section 4, a number of methods are proposed for calculating, collecting, and measuring FJR parameters. This model will be completed and modified by applying parameters uncertainties. Also in this section, simulation results and experimental outputs are presented.

2. Modeling

The scheme of the single link FJR is shown Figure 1, which includes rigid arm, flexible joint, and actuator. Usually actuator is a combination of gearbox and electromagnetic motor as a compact unit fixed on the work-table by a rigid base. Several fixed ball-brings on parallel bases are used to fixing shafts and reducing effect of friction. Gearbox generated torque is transferred by fixed shafts from gearbox-out to input of spring and from output of spring to joint of link. Solenoid-type spring is selected for coupling [5]. In Section 2.4, it will be shown this spring is employed in torsion form and also has weak nonlinearity.

2.1. Base Model

Dynamic models for link and actuator sections can be created by two Lagrange differential equations as follows [6]:

d SLi SLi

dt~6&i SWi

i = 1; link of system, i = 2; actuator of system,

Figure 1: Scheme of a single link FJR.

where Li represents difference between the kinetic and potential energy in link and actuator as follows:

Li^ti) = Ki(diA) + prm,

xir pi in link and actuator are expressed by the following equations:

pi(6i) = 2mL • g • i • cos(0i) + mp • g • i • cos(6i),

1 • 2 1 • 2 Ki(Qi,Qi) = 2Il(0i)2 + 2mp(id0 ,

Ui = Uc - Ufl - Uei, ( P2(02) = 0,

actuator

i 2 M 62,62) = 2 J • {62)1,

U2 = Uin - UC - UFa - UE2,

where Uei, and UFa, UFL are axles position shaft-encoders torques and friction torques of link and actuator. Also UC and Uin are the transfer torques from actuator to link joint and gearbox out-torque. In these equations, mL, mP, IL, and J are mass of link, payload, inertia moment of arm, and actuator, respectively. Also i, g, Q\, and d2 are the length of arm, the gravity acceleration constant, the rotation angle of link, and actuator-shaft to system base ({Ao}), respectively. By replacing L\ and L2 in the Lagrange equations in (2.1), FJR dynamics is described by

(Il + mp(£)2)0i + ^^mL + m^jg • i • sin(6i) = Uc - Ufl - Uei, J • 02 + B • ^2 = Uin - Uc - UFa - UE2-

-1 -2 -3 -4 -5

W Fri ction torque

Coulomb

- - Cf Viscous

Speed (rad/s)

Figure 2: Friction characteristic of a compact motor and gearbox.

2.2. Friction Modeling

As it was mentioned UFa, UFL are friction torques in the link and actuator which are generated by viscous and Coulomb frictions. Dynamic of friction torques can be written as [7]

UFL = FVL • 01 + Ucl,

. (2.6) UFa = Fva • 02 + Uca-

Generally Fvl, FVa are constants which are related to type of surfaces and quality of lubrication. Usually, using lubrication and ball-brings can decrease Coulomb friction on the joint and arm sections. But lubrication cannot be used in brushing of DC motor, so Coulomb friction has always high value in gearbox and motor section. In this paper, UCa is selected as follows [8]:

[Cf ■ sgn(02),

Un< W ■ S

¿2 / 0,

& = 0,

Sf ■ sgn(Un), otherwise,

where Cf, Sf, Uin, and W are coefficient of kinetic and static frictions, applied torque and proportional weight respectively. It is also proved that Coulomb friction is generated by kinetic and static frictions and generally Cf is less than Sf (Cf < Sf ). As an example, Figure 2 shows the modified figure of [9] for friction characteristic of a compact motor and gearbox.

Often, designing a nonlinear controller for nonlinear systems will be conducted to use of numerical methods, so dynamic equations should have Taylor expansion. There is no Taylor series for equations similar to (2.7), because it has some nonderivativeable points on

defined domain. This problem can be solved by using the following approximate equation (2.8) instead of (2.7). This equation is derivativeable in all defined domain

Uca = ^tanh(C ■ +

A ■ 02

1 + b ■ é\t

02 e K.

Also Sf and Cf of (2.7), which have positive values, are functions of A, B, ¡5 as follows:

Sf = Max

Cf = P,

tanh(C ■ 02)

A ■ 02 1 + B ■ Q\

where C is a positive number which represents the slope of UCa in low speeds. Lubrication on ball-brings and gearbox is usually used in order to decrease Cf and Sf. But lubrication is not possible between collector and coals in DC and universal motors, because it may cause electrical conductivity problem. Based on (2.7), Coulomb friction torque has high value at low speeds, so if rotor rotates without stopping, then starting torque is eliminated and Coulomb friction is decreased, which is shown experimentally in Section 4.

2.3. Motor Modeling

In most papers, the model of electromagnetic motor, which is used as driver, is described by a linear and low-order differential equation. In spite of simplicity, this model is used in most usual applications yet. But certainly linear model cannot be used in robots because the direction of rotor is always changed between CW and CCW status [10]. So nonlinearity and saturation must be considered in motor model especially in lightweight robots. Also in this paper, robot is actuated by a DC motor. Therefore, robot input torque is a function of armature current of DC motor as follows [10]:

Uout-motor = Km ' ta '

\ ia(Max) /

where ia, Km, and ia(Max) are armature current, positive value as armature current to out-torque conversion constant, and maximum allowable current in armature, respectively. If the armature current is equal or greater than the maximum allowable current, saturation would happen in rotor and stator cores. In this paper, saturation is modeled by smooth function tanh(-) as an approximation of saturation dynamics. Uout-motor is applied to input of coupling spring as Uin. Also, armature current has a relation with the motor's electrical parameters as follows [15]:

ia = 7-Ô (Vm - KbQm), (2.11)

Ls + ^

Fuzzy TSK controller

Pulse from vibration source

Figure 3: Model of DC motor and input mixed voltage.

Figure 4: Tensional and spastic modes.

where R, L, Vm, 6m, and Kb are armature resistance, armature inductance, input voltage, motor speed at gearbox out, and speed to back e.m.f voltage conversion constant, respectively Figure 3 shows a model of DC motor that is used as robot actuator. In this model, motor input voltage is a mixture of nonlinear controller output and symmetric pulse vibration voltages.

In order to actuate link in any position, it is necessary that motor out-torque must be always more than required torque.

2.4. Coupling Model

In FJRs, flexible gearbox out-torques is applied to robot joint as actuator. But in this paper, flexible gearbox was replaced by a set of rigid gearbox and tortuous solenoid spring for having high flexibility. In this working mode, spring coefficient is very low. Other problem, which complicates designing process, is nonlinearity of spring. As it is shown in Figure 4, effective diameter of spring can be changed in Tensional and Spastic modes, because distance of spring bases is always fixed.

Increasing and decreasing in spring diameter on Tensional and Spastic modes can cause higher nonlinearities in FJRs specifications and also model may become asymmetric according to the following equations [11]:

U = j -J

U+ = a+y + a+y3 + a+y5, y> 0,

U- = a-y + a-y3 + a5^5, y < 0, (2.12)

a+ = a- > 0 or a+ = a- < 0 i = 1,3, and 5,

where y is defined as follows:

y = {02 - Q1 | Q1 e R, Q2 e R}. (2.13)

All a+ and a- are Tensional and Spastic modes' polynomial coefficients. In (2.12) it is assumed that coefficients of even powers of U+ and U- polynomials are all zero, because of symmetric structure to origin. Generally (2.12) is not derivable on the origin; therefore, it will not have Taylor expansion. In order to eliminate conflicting, the following equation can be used instead of (2.12) in defined domain as a continues and derivable approximation

Uc = Si^ + a3y3 + a5y5, (2.14)

where Sj are new polynomial coefficients with uncertainty that are bounded between lower and upper size of a*. Therefore, aLi can be defined based on a+ and a- as follows:

Si = ai (1 + Wai),

a- + a

L | (a- & a+) e R+ or (a- & a+) e R-

(2.15)

I a - a

a- + a

where wai represents maximum deviation a+ and a- from aj. Figure 5 shows how to connect spring to bases. In this paper, a new method is proposed for symmetrical transferring of torque from actuator to spring where 50% of coupled torque is applied to the spring at A by OA link and the rest is applied to A' by OA' link and the fixed arc of spring (AA'). Spring properties may change if connection points (A, A') are fixed by high-temperature methods such as soldering. Therefore, to avoid this problem spring and bases were fastened by screws and spanners in laboratory experiments. Also applied force on connection points (A, A') is equal to U-1|OA|, which can be reduced by incrementing the bases length (OA, OA'). This structure has some advantages as follows:

(1) symmetrical transmission for torque by two sides of spring bases,

(2) decrease in applied forces on connection points,

(3) decrease in friction forces of boll brings.

3. Uncertainty and TSK Fuzzy Controller

In this section, the new model of FJR is presented, which has uncertainty as tolerance of parameters and disturbance inputs. It is also assumed that states of system are full available and system equations in state space are described as follows [12]:

x = f (x) + g1(x)w + g2(x)u,

z = h(x) + l(x)u,

where x e Rn is the state vector, w e Rmi is the disturbance input vector, u e Rm2 is the control input vector, and z e Rp is the penalty variable vector. Also, it is assumed that functions f (x), g1(x), g2(x), h(x), and l(x) are defined and are smooth in the neighborhood of X on Rn.

3.1. Uncertainty

By using the discussions of previous section, new state space model of FJR can be written as follows:

x1 = x2, x2 = f2(x1,x2,x3 ,W),

x3 = x4,

x4 = f4(x1,x3,x4, x5, W), x5 = f5(x4,x5,W) + g25(x) ■ Vm, xx*6 — xd x1-

These equations are written by using (2.5), (2.6), (2.10), (2.11), and (2.14), where wd is the input disturbance vector. The state vector on (3.2) is defined as:

X1 = 01, X3 = 02, X5 = ia, x2 = 0i, x4 = 02, x6 = e ■ dt.

Also the models of shaft encoders were eliminated for simplicity, because of their low effects. Equation for X6 is added to FJR describing equations to decrease the tracking error. Where xd and x6 are the input command as desired position and integral of tracking error. In this paper, tolerances of mechanical parameters such as mass and length of link are also considered, because these parameters can be changed easily. It is assumed that measuring of fixed mechanical and electrical parameters in compact motor and gearbox has been done with enough precision. Tolerance of parameters can be defined as follows:

mL = m1(1 + a1w1),

mP = m2(1 + a2w2), (3.4)

e = e1(1 + a3W3),

where, a1, a2, and a3 are positive numbers as maximum deviations of mL, mP, and e respectively. Also w1, w2, and w3 are external input signals, which are limited by |wi| < 1. Therefore, wd is defined as follows:

w = [w1 ,w2,w3,wa1,wa3,wa5,xd ]

3.2. TSK Controller

T-S fuzzy models are universal approximations [13,14]. The form of a nonlinear system is

x = F(x) + b(x)u,

where x = [x1, x2,...,x„]T is the state vector. F(x) and b(x) are the nonlinear functions of the plant and u is the control signal u(t) = [u1(t), u2 rule in the fuzzy rule are in the following form: Ri: If Z1(t) is Mi1 and ... and Zg(t) is Mig

Then <|x(t) = Ax(t) + Bu(t), (3.7)

y(t) = Cix(t).

In the ith implication, r is the number of fuzzy rules, Mj are the fuzzy sets, Ai and Bi are the n x n and n x p constant state matrixes, and x(0) is the initial state vector:

. (t) h Wi[Z(t)][AiX(t) + Bju(t)]

X(t)~ sr_i Wi [Z(t)] '

r (3-8)

^i[Z(t)] = r^^L _* X(t) _ ¿^i[Z(t)][AiX(t) + Biu(t)],

Xi_i Wi[Z(t)]

Wi is its corresponding membership function, wi[Z(t)] = Tlg^Fj[Zj(0], and Fj[Zj(t)] is the grade of membership of Zj(t) in Fj (fuzzy sets).

The final state of fuzzy systems and the output is represented as

x (t) = ^ fiiAix(t) + ^ ¡A.iBiu(t), i=1 i=1

y(t) = Tj ^iCix(t)-i=1

Considering (3.7), design of the controller via parallel distributed pole placement results in another fuzzy system with r rules in the form of the following: if Z1 (t) is Mi1 and Z2(t) is Mi2 and ... and Zg(t) is Mig

Then u(t) = -Kix(t), i = 1,2,...,r. (3.10)

Then the output of the fuzzy state feedback controller is u(t) = - ^ri=1 fii[Z(t)]Kix(t), where the eigenvalues of Ai - BiKj are in the left half-plane. By substation (3.10) in to (3.7),

X(t) _ X - BiKj)x(t). (3.11)

i_1 j_1

The nonlinearity in system dynamic is

z1 = cos(x1), (3.12)

T-S fuzzy modeling of FJR results in two rules as follows

Z1(t) = cos(x1(t)) = ^ g Mi(z1(t))b^jx1(t), (3.13)

-1 0 1 " Figure 6: Membership function M1(z1(i)) and M2(z1(i)).

\ ^ — \

\ Ma( Za(t))= 0N / N/ /\

'"■M2 (za (t))/ 0 \ ' N \

, • \ __ . — ■ \

where the membership functions are

Ml(zl(t)) =

Z1(t) - (2/n)cos-1(z1(t)) (1 - 2/n)cos-1(z1(t)) '

Z1(t)/ 0, otherwise,

(3.14)

M2(Z1(t)) =

cos-1 (Z1(t)) - Z1(t) (1 - 2/n)cos-1(Z1(t))' 0,

Z1(t) = 0, otherwise,

(3.15)

b is slope of sectors and defined as

Z1(t) = bx(t), b = 1 or ^-2),

b1 = 1 (Zero),

= — (Not Zero) n

(3.16)

Membership functions of Z\ are titled as "Zero," "Not Zero." These membership functions are shown in Figure 6. Hence, the proposed two rules of T-S fuzzy modeling can be defined as follows.

Rule 1 (If Z(t) is zero). Then x(t) = Aix(t) + B ■ uc

mglbi + k k

0 0 0 -J

B Kb Km

(3.17)

Rule 2 (If z(i) is not zero). Then x(i) = A2x(i) + B ■ uc

mglb2 + k I

B Kb Km

(3.18)

The state feedback controllers are obtained as follows:

Ki = [-1069.5 - 75.988 574.84 73.981], K2 = [-913.34 - 44.88 618.26 73.981].

The final state feedback gain is K = ¡1 [Z(t)]Ki + ¡2 [Z(t)]K2.

(3.19)

4. Experimental Results

In order to analyze the proposed controller performance, a single flexible joint robot test-bed was prepared in Islamic Azad University of Ahar. Figure 7 shows the photo of this robot. This FJR has two main sections, electromechanic and control sections. Electro-mechanic section is split to compact gear-boxed DC motor as actuator, solenoid spring as torque transformer, and plastic rigid link as arm. Also, in order to avoid undesired movements, main shaft from gearbox output to joint of arm is fastened by fixed parallel bases on main foundation.

Controller section has four high speed microcontrollers. Gearbox-shaft and link positions are measured by high precision shaft encoders as measured output signals. Shaft encoders signal wires are connected to two separable and independent microcontrollers as inputs. These outputs are wired to the interrupt inputs of microcontrollers. Outset, the pulses of shaft-encoders are decoded and then rotational angles are separately calculated by a special algorithm in microcontrollers.

In the next phase, decoded information is sent to main computer by RS232 protocol as serial communication and serial to USB converter. The main controller is a software controller programmed in Visual Basic. The controller on main computer gets gearbox-shaft and link positions by serial line and control law coefficients of a file that was already created by the other program in MATLAB. Special numerical code as DC motor voltage is created by use of control law and positions information in controller.

The controller on main computer gets gearbox-shaft and link positions by serial line and control law coefficients of a file that was already created by the other program in MATLAB. Special numerical code as DC motor voltage is created by use of control law and positions information in controller.

Then this code is sent to the fourth microcontroller by serial line. In order to create apposite commands for generating PWM form DC voltage as motor power supply, DC voltage is generated by four IGBTs as a controllable bridge and PWM method. According to the controller commands, produced voltage can be altered in range from -24 V to +24 V.

Figure 7: Single link laboratory FJR test bed.

MEMS tilt sensors

Figure 8: Scheme of single link FJR and controller test-bed.

4.1. System Mechanical and Electrical Parameters

Table 1 shows values of link and spring mechanical parameters.

Also, Table 2 shows values of compact gear-boxed DC motor set electrical and mechanical parameters. Figure 8 shows scheme of single link FJR and controller test-bed.

Coulomb and Viscose friction coefficients are presented for link and compact gear-boxed DC motor in Table 3. This table also shows coefficients at two states: without vibration and with vibration signal. Of course vibration signal had proper quantities during experimentally tests.

Discrete Dynamics in Nature and Society Table 1: Link and spring properties.

Parameters Nominal values Tolerance

Joint stiffness a1 = 3.35 a3 = -0.16 1.5% 97%

Link Length L = 0.53 6%

Link mass M = 0.12 9%

Gravity coefficient g = 9.8 1%

Inertia I = 0.083 4%

Table 2: Compact DC motor and gearbox set properties.

Parameters Nominal values Tolerance

Coil resistance r = 9 ohm 4%

Coil inductance l = 0.1 H 5%

Torque const. K^K-m — 3.3 3%

Back emf const. Kb = 5.62 3%

Gearbox ratio N = 100 : 1 —

Motor inertia J = 2.15 3%

Operating Max. Voltage V = 36 —

Max. armature current ia (Max) =4A —

Table 3: Values of friction coefficients.

Parameters No vibration With vibration

Link Coulomb const. Fc-link = 0.05 —

Link Viscose const. F^-link = 0.1 —

Motor coulomb const. ß = 0.83, A = 3.6 ß = 9, c = 1 ß= ß 0.3, A = 1.5 = 3.1, c = 1

Motor Viscose const. Fc-motor = 1.43 —

Table 4: Sf for some vibration signals.

Sf (N. Vib) = 0.98 2vp-p 4vp-p 6Vp-p

5 Hz 0.68 0.27 0.37

10 Hz 0.52 0.29 0.4

15 Hz 0.63 0.35 0.39

Table 5: Cf for some vibration signals.

Cf (N. Vib) = 0.82 2vp-p 4vp-p 6Vp-p

5 Hz 0.55 0.2 0.27

10 Hz 0.46 0.23 0.24

15 Hz 0.52 0.3 0.15

4.2. Compensated Model for Frictions

The friction model can be compensated by applying a pulsation form signal (vibration signal) with special frequency and amplitude. Tables 4 and 5 show values of Sf, Cf by some vibration signals. These tables were gathered by various experimentally tests on test-bed FJR.

-11-1-1-1-1-1-1-1-1-1-

0 2 4 6 8 10 12 14 16 18 20

Time (s)

- Motor speed (rad/s)

...... Motor volatge (V)

Figure 9: Sf, Cf measuring test.

Figure 9 shows a test result for measuring of Sf, Cf. This test was done with square wave signal as a vibration signal with f = 15 Hz and vP-P = 4v.

4.3. Effect of Vibration Signal on Friction and Backlash

In order to show advantages of proposed controller with vibration signal, a number of tests were done on FJR for various input commands. These tests are repeated with commands 45', 120', and 150'.

Figure 10 shows that response to command 45' with vibration signal is more accurate than without it and also shows that link angle was sat to input command at t = 8 s by vibration signal, but it has ringing form on response for without vibration signal even after t = 16 s. Figures 11 and 12 show response of FJR to commands 120' and 150' that are in unstable domain. It means that FJR has been stabled even in unstable domain by proposed controller. Certainly it is predictable that link response to these commands would have high settling time overshoot.

4.4. Wide Changing in Payload

In order to show advantages of proposed controller law to traditional. The other tests were done on FJR by use of three different payloads.

Figures 13 and 14 show the test results of controller k(x,mp) and k(x) in the same conditions. Similar to Section 4.3, optimal vibration signal has been applied in these tests.

It can be concluded from Figures 13 and 14 that new controller has more compatibility with mp. Also, in addition to solving backlashes problem and decreasing of friction effect, the proposed controller has more accuracy and higher performance compared to traditional controllers.

0 5 10 15 20

Time (s)

- Vib. (on)

- Vib. (off)

Figure 10: Response to Xd = 45' with vibration signal.

- Vib. (on)

- Vib. (off)

Figure 11: Response to Xd = 120' with vibration signal.

4.5. Uncertainties and Disturbance

Robustness of the proposed controller is verified by several experimentally tests. The tests have been conducted by changing the link length. Figure 15 shows the step response of link position in controlled system by proposed control law for link uncertainties. Also Figure 16 shows the response of link position to step command 60', where external disturbances are applied as impact forces to link.

The results of tests show that the proposed controller can guarantee the system stability against parameter tolerances and external disturbances. The effect of Q6 on new controller is verified by some tests, the results are shown in Figure 17.

(g 150

jn 100

--——^

Time (s)

- Vib. (on)

- Vib. (off)

Figure 12: Response to Xd = 150' with vibration signal.

^ 40 in

10 Time (s)

- mp = 500 gr

mp = 600 gr - mp = 700 gr

Figure 13: Effect of mv changing with k(x, mv).

Results of Figure 17 show operation of controlled system can be improved by fine adjusting of Q6. Inappropriate value for Q6 can drive system toward unstable mode. Motor current curves are shown in Figure 18. They are measured in two modes: with vibration signal and without vibration signal.

Figure 18 shows, if motor control voltage is mixed with vibration signal, armature current will have smooth function with low variations after t = 6 s. It is needed to mention that in all tests, position of link and gearbox shaft have been measured by use of shaft encoder and data acquisition board.

Time (s)

- mp = 500 gr

mp = 600 gr - mp = 700 gr

Figure 14: Effect of changing with k(x).

Vib (on)

Time (s)

— -%20 len. of link +%0 len. of link

— +%20 len. of link

Figure 15: Response to Xd = 60' with link uncertainty.

5. Conclusion

In this paper, Takagi-Sugeno fuzzy approach for modeling nonlinear systems is represented. Using this method, the nonlinear elements of the system can be replaced by linear subsystems and result can be achieved in a simpler model and less complicated equations. Also saturation problem in actuator and Coulomb friction were modeled by derivable and smooth models

-20 -1-1-1-1-1-1-1-

0 5 10 15 20 25 30 35 40

Time (s)

- Vib. (on)

- Vib. (off)

Figure 16: Response to Xd = 60' with external disturbance.

Vib (on)

- Q = 80 - Q = 100

- Q = 200

Figure 17: The effect of Q6.

instead of classical methods in this paper. Moreover, experimental results of proposed controller on FJR showed that TSK fuzzy-based controller has an ability to stabilize and control FJR with acceptable performance. Moreover, it was shown that backlashes and vast starting torque problems in gearbox can be solved by use of vibration signal. Also proposed approach can be applied to any multilink flexible on stable domain with a good performance.

Time (s)

- Vib. (on)

- Vib. (off)

Figure 18: Applied voltage to DC motor for Xd = 120'.

Also, it can be expanded as an extension of the technique to the more challenging cases of systems with unstable dynamics with minimum performance, such as flexible link robots.

References

[1] B. Siciliano, "Control in robotics: open problems and future directions," in Proceedings of the IEEE International Conference on Control Applications, vol. 1, no. 1, pp. 81-85, September 1998.

[2] M. W. Spong, K. Khorasani, and P. V. Kokotovic, "An integral manifold approach to the feedback control of FJRs," IEEE Journal of Robotics and Automation, vol. 3, no. 4, pp. 291-300,1987.

[3] K. Khorasani, "Adaptive control of FJR," IEEE Journal of Robotics and Automation, vol. 8, no. 2, pp. 250-267,1992.

[4] W. H. Ho and J. H. Chou, "Design of optimal controllers for takagi-sugeno fuzzy-model-based systems," IEEE Transaction on Systems, Man, and Cybernetics A, vol. 37, no. 3, pp. 329-339, 2007.

[5] M. W. Spong, "Modeling and control of elastic joint robots," Journal of Dynamic Systems, Measurement, and Control, vol. 109, pp. 310-319,1987.

[6] J. D. Zhang and S. T. Zhang, "Controller design of T-S fuzzy systems with standard fuzzy partition inputs," in Proceedings of the 4th IEEE Conference on Industrial Electronics and Applications, pp. 31013106, May 2009.

[7] M.-Y. Hsiao, C.-H. Liu, S.-H. Tsai, K.-L. Tsai, P.-S. Chen, and T.-T. Chen, "A Takagi-Sugeno fuzzy-model-based modeling method," in Proceedings of the IEEE International Conference on Fuzzy Systems, pp. 1-6, July 2010.

[8] J. Chiasson, Modeling and High-Performance Control of Electric Machines, Wiley Interscience, 2005.

[9] S. Ozgoli and H. D. Taghirad, "A survey on the control of flexible joint robots," Asian Journal of Control, vol. 8, no. 4, pp. 332-344, 2006.

[10] G. Palli, C. Melchiorri, and A. De Luca, "On the feedback linearization of robots with variable joint stiffness," in IEEE International Conference on Robotics and Automation, pp. 1753-1759, 2008.

[11] J. Huang, "An algorithm to solve the discrete HJI equation arising in the L2 gain optimization problem," International Journal of Control, vol. 72, no. 1, pp. 49-57,1999.

[12] W. M. Lu and J. C. Doyle, "Hœ control of nonlinear systems: a convex characterization," IEEE Transactions on Automatic Control, vol. 40, no. 9, pp. 1668-1675,1995.

[13] W. H. Ho and J. H. Chou, "Design of optimal controllers for Takagi-Sugeno fuzzy-model-based systems," IEEE Transactions onSystems, Man and Cybernetics A, vol. 37, pp. 329-339, 2007.

[14] W. H. Ho, J. T. Tsai, and J. H. Chou, "Optimal controller design of time-varying TS-fuzzy-model-based systems via Chebyshev series and genetic algorithm," in IEEE International Conference on Systems, Man and Cybernetics, vol. 3, pp. 2168-2173, 2005.

Copyright of Discrete Dynamics in Nature & Society 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.