Abstract

This paper presents a new direct adaptive fuzzy controller and its effectiveness is verified by investigating the damping force tracking control of magnetorheological (MR) fluid based damper (MR damper in short) system. In the formulation of the proposed controller, a model of interval type 2 fuzzy controller is combined with the direct adaptive control to achieve high performance in vibration control. In addition, (H infinity) tracking technique is used in building a model of the direct adaptive fuzzy controller in which an enhanced iterative algorithm is combined with the fuzzy model. After establishing a closed-loop control structure to achieve high control performance, a cylindrical MR damper is adopted and damping force tracking results are obtained and discussed. In addition, in order to demonstrate the effectiveness of the proposed control strategy, two existing controllers are modified and tested for comparative work. It has been demonstrated from simulation and experiment that the proposed control scheme provides much better control performance in terms of damping force tracking error. This leads to excellent vibration control performance of the semiactive MR damper system associated with the proposed controller.

1. Introduction

Recently, magnetorheological (MR) fluids based devices are actively studied and some of devices such as shock absorber for passenger vehicles are commercialized. Many MR devices or MR systems especially are popularized in the field of aerospace, ship, automobile, and civil engineering for effective control of unwanted vibrations. It is well known that MR fluid is considered as one of smart fluids in which the rheological properties can be controlled by external stimuli such as magnetic field. On the other hand, the development of new control system for MR devices or systems is also continuously undertaking in which the objective goal is to guarantee both stability and robustness against disturbances or/and uncertainties. Based on this objective, several types of robust controllers have been developed in both indirect adaptive control and direct adaptive control manners. A model of direct adaptive control for output tracking was studied in [1] where infinity tracking technique was used for adaptive model. An adaptive fuzzy control for MIMO system was presented in [2] in which the Lyapunov function and mean-value theorem were combined with a fuzzy model to build adaptive model. An observer of adaptive fuzzy controller was presented in [3] in which the back-stepping design and the supply changing function technique were applied. An input saturation problem was treated with adaptive feedback fuzzy control in [4] where the back-stepping technique and Lyapunov function were used in building the adaptive model. A direct adaptive control using type 2 fuzzy model was presented in [5] where infinity tracking technique was applied to attenuate the fuzzy approximation error and the system uncertainty. The linear matrix inequality fuzzy Lyapunov method was also studied in [6] which concentrated on modeling errors between real system and fuzzy model. Recently, a novel of adaptive sliding mode control was developed by authors in [7] in which a sliding mode controller and Lyapunov method were used for high control performance. Vibration control of MR damper using an adaptive neurofuzzy inference system was also conducted by authors in [8].

In application of feedback/feed-forward controls for tracking force of MR devices, there are several researches. An idealized hysteretic model for MR was studied in [9] using a nonlinear hysteretic bi-viscous model. For seismic response reduction, a clipped-optimal control using acceleration feedback was studied in [10] in which the properties of MR damper were tested to provide good vibration control performance. A model of MR damper for bridge application was derived in [11] using a semiactive tuned mass damper. A combination of the feed-forward and feedback control was also presented in [12] for torque control of MR brake. On the other hand, in research of interval type 2 fuzzy model (IT2FM), many studies on the development of controller type have been undertaken. A comprehensive research about type-reduction of three groups, Karnik-Mendel (KM), closed-form representation, and combined KM-closed form representation, were presented in [13]. In this study, EIASC (Enhanced Iterative Algorithm with Stopping Condition), Wu-Tan (WT) algorithm, and Nie-Tan (NT) algorithm were used to reduce calculation cost. Another research concentrating on KM algorithm in IT2FM was studied in [14] in which all methods were evaluated based on KM method. In order to solve the bottle-neck problem in calculating of type reduction of IT2FM, a new method of combination of KM and NT method was presented in [15]. It is recognized from the above survey that the development of a new adaptive fuzzy controller which is robust to disturbances or/and uncertainties is a hot issue in many control systems including MR device or MR systems. The development of new IT2FM especially is a new horizon in design of new controllers to guarantee both stability and robustness against disturbances and uncertainties [1315].

Consequently, the technical originality of this work is to propose a new direct adaptive fuzzy controller using both IT2FM and adaptation method based on -infinity tracking technique for the robust damping force tracking control of a cylindrical MR damper. Thus, a final goal is to demonstrate superior control performances of the proposed new direct adaptive fuzzy controller compared with two existing adaptive fuzzy controllers [1, 2]. In order to clarify the difference from the previous work done by the authors [7], it is noted here that the previous study deals with the indirect adaptive fuzzy controller, while current study treats the direct adaptive fuzzy controller whose structure is much different from the indirect method. In order to achieve the final goal of this work, as a first step, a new interval type 2 fuzzy control system is formulated. Then, a new type of adaptive controller is designed based on infinity tracking technique and Lyapunov method. Based on the mathematical model, a closed-loop control structure is built and experimentally realized for a small-sized MR damper. MR fluid used in this work is manufactured in the laboratory using plate-like iron particles. Its field-dependent rheological characteristics are different from general MR fluid containing spherical iron particles. In order to demonstrate superior control performance of the proposed controller, both tracking control of desired damping force of MR damper and vibration control performance are evaluated and compared with existing two controllers.

2. Interval Type 2 Fuzzy Logic System and NT Algorithm

The model of interval type 2 fuzzy logic system (IT2FLS) has been developed from the type 1 fuzzy logic system (T1FLS). The structure of IT2FLS includes five components: fuzzifier, rule base, inference engine, type reducer, and defuzzifier. The rule base of IT2FLS can be expressed as follows: where are fuzzy sets and is the centroid of a consequent IT2FLS.

In this study, NT algorithm is utilized for finding output of IT2FM. In general, NT algorithm shows the best algorithm in calculation with minimum computational cost. Hence, the application of NT algorithm will save time of fuzzy calculation. The relationship between NT algorithm and IT2FM is shown in Figure 1. The defuzzification of the IT2FLS-NT is calculated from several steps as follows.

Step 1. Set up input vector . This vector is also vector state of the system whose component is defined as .

Step 2. Find , where is lower value and is upper value of IT2FLS .

Step 3. Calculate lower and upper firing values of IT2FLS where is lower firing value and is upper firing value.

Step 4. Use NT algorithm for calculating of . Here is lower interval type 1 value and is upper interval type 1 value given by

Step 5. Calculate defuzzification of IT2FLS as follows:where ,  , and . It is remarked that   is consequent vector and , are lower and upper consequent membership vectors, respectively.

3. New Direct Adaptive Fuzzy Controller

Consider a nonlinear system is governed by the following equation:where and are two unknown nonlinear function vectors, is control function, is an external disturbance vector, , where is upper bound of , and is the state vector of the system. The function and can be expressed in two parts. One is a nominal function and the other is unknown bounded uncertainty as follows:where , , , and . and   are two positive vectors. Hence, the system (4) can be rewritten aswhere is uncertain part (disturbance) which can be redefined as . Then, the relationship between system (6) and IT2FLS is determined based on (3) as follows:where and , . Note that , are the centroids of consequent vectors and , are consequent membership vectors of , , respectively. The difference between a desired output and the measure output is defined by . Thus, the error vector for the system is defined as , and a vector is the chosen coefficients such that all of the roots of the polynomial are located in the open left-half complex plane [1, 2].

In relation with fuzzy control, by assuming the disturbance of a direct adaptive fuzzy control is expressed as follows:where is a fuzzified value of . Using (6) and (8), the derivative of is expressed as follows:whereNow, define the minimum approximation error as follows:Then, substituting (11) into (9) yields the following equation:Let . From (12), the equivalence control is defined aswhere is the estimates of . It is noteworthy that the equation is established without the minimum approximation error . Due to the presence of fuzzy approximation error and equivalent control term, using only is not sufficient to ensure the stability of the closed-loop system. Therefore, it is necessary to add a robust compensator to deal with it as follows:where is an adaptive parameter, and is a solution of the following Riccati-like equation:where , is a prescribed attenuation level, and , is consequent membership value of the IT2FLS. The resulting fuzzy control is determined as follows:By substituting (16) into (12), (16) can be rewritten aswhere .

Now, in order to prove the stability of the control system, consider the Lyapunov function candidate as follows:Then, the derivative of (18) is expressed as follows:From (19), adaptation laws are established as follows:Applying (20), (19) is determined as follows:By integrating (21) from to , the following equation is obtained:where . The group is always more than or equal to zero value. Hence, from (22), the remained group is analyzed as follows:Since the value is set by and , the error of matrix and the gain of must be equal to or less than . Hence, the infinity tracking form is satisfied [2, 16], and then (21) is always less than zero value (Lyapunov stability). Now, from the boundedness of parameters of , is guaranteed by closed sets defined as follows:where , , , are the choosing parameters. Hence the adjusted adaptation laws are redefined as

To observe the states of system, an observer for the nonlinear system is needed to converge the estimated state of . In this work, the observer is constructed by the Luenberger observer [17] which is defined as follows:where is a finite gain vector, , is the observability matrix of the system. From the above analysis, the closed-loop of the proposed controller is shown in Figure 2. In Figure 2, adaptation laws use output of IT2FLS-NT model and error of system for calculating adaptive values of and . These values are input of fuzzy controller which is input of observer and plant modules. Output of plant is input of IT2FLS-NT model and base for finding error of system which is also input of IT2FLS-NT model. Value of output observer is also base for calculating error of system.

4. MR Fluid and MR Damper

4.1. New MR Fluid

As mentioned in Introduction, MR fluid is suspension of microsized magnetic particles dispersed in a nonmagnetic carrier liquid [1820]. It is smart fluid whose properties can be controlled by the magnetic field intensity. In the absence of magnetic field, the properties of the fluid are isotropic, while in the presence of magnetic field the magnetized particles form chain aligned in the direction of the field which results in the appearance of the field-dependent yield stress. Reversible transition from solid to liquid is also possible. Most of MR fluids developed so far have spherical iron particles. Recently, it has been reported that MR fluid containing plate-like iron particles provides better performance than MR fluid featuring spherical iron particle in terms of sedimentation characteristic [21, 22]. Therefore, in this work, a new type of bidisperse MR fluid is used for the application to MR damper. In order to undertake this work, two different sizes of plate-like iron particles are prepared. It consists of iron plate-like particles by weight of two different average diameters suspended in a heavy paraffin oil (64cSt) using mechanical stirrer. Iron microsized particles (Industrial Metal Powder, Pune, India) having average particle size 2 μm (small size) and 19 μm (large size) are used to prepare MR fluid. The structural property of micron-sized magnetic particle is characterized by an X-ray diffractometer, D2-phaser (Bruker XRD, Germany). Bidisperse MR fluid sample is prepared with a variable weight fraction (weigh of small size particles ()/weight of large size particles ()) of 4.00. The particle density is about 7.8 g . A small amount of stabilizer is then added, stirred, and homogenized for a long time till good dispersion is achieved. Figure 3 shows the X-ray diffraction pattern of plate-like iron particles recorded at room temperature. The XRD result matches with the reported value of Fe given in the literature [23]. The magnetization measurement is carried out using a home-built magnetometer. The magnetization saturation of the particles used in this work is about 1250 ± 15 kA m−1. In a field gradient, alignment of the magnetic moment of large size particle is faster compared with the small-sized particle. Under this condition, initial susceptibility and saturation magnetization of larger particles is higher compared with the smaller particle. As a result, the fraction magnetization rises as the weight fraction increases.

Rheological properties in static mode for the MR sample are measured using a Physica MCR 301 (Anton Paar, GmbH, Austria) parallel plate rheometer coupled with a commercial magnetorheological device (MRD 180/1T magneto-rheological cell). A diameter of plate is 20 mm and a gap between parallel disks is kept to 1 mm during whole experiment. The magnetic circuit is designed so that the magnetic flux lines are normal to the parallel disk. Magnetic field is perpendicular to the velocity direction. When the magnetic field is applied the shear stress increases quickly with increasing shear rate as shown in Figure 3(a). The shear stress increases with increasing magnetic field strength due to enhancement of magnetic dipole-dipole interactions between the particles. This means that the MR fluid behaves as a Bingham fluid, where a yield stress is needed to initiate flow in the fluid. The Bingham model only applies for the linear portion of the stress-shear rate curve at high shear rate where the suspension flows like a Newtonian liquid [24]. According to this model shear stress is obtained bywhere is the shear stress, is the Bingham yield stress caused by the applied magnetic field strength , and is the field independent plastic viscosity defined as the slope of the shear-stress curve at higher shear rate. The dynamic yield stress value increases with increasing magnetic field strength as shown in Figure 3(b). The data fitted in the second order polynomial equation is expressed as follows:In the case of bidisperse suspension of micron-sized magnetic particles due to low applied magnetic field, a pair of short-chain segments forms from the large size particles first. When the external field increases, the substitution of small particles by large particles in a chain would have the effect of weakening the chain structure formed under shearing. The main reason to use bidisperse MR fluid is because it produces a higher yield stress compared with monodisperse suspension. The maximum yield stress generated is 38 kPa at 228 kA m−1. In application point of view, the yield stress is the most important rheological parameter. The sedimentation rate of this MR fluid is 0.3% per day. The characteristics of bidisperse MR fluids are listed in Table 1.

4.2. Model of MR Damper

In this work, a simple dynamic system model installed with MR damper shown in Figure 4 is considered. From Figure 4, the equation of motion can be derived bywhere is the mass of piston inside damper, is viscous damping, is the accumulator stiffness, is external force, and is displacement. The value of belongs to the damping force of MR damper given bywhere is signum function, . Substituting (31) into (30) and by defining the displacement as a state variable yields the following equation of motion in a state space form:From (32), control of the system is expressed by , , and . To change the system (32) into direct controller, the value is chosen by 1. In this application, the semiactive control is applied in damping force control of MR damper. Hence, the damping force due to the yield stress of MR fluid can be expressed as follows [9]:where is the postyield viscous damping, is the preyield viscous damping, is the yield force, is velocity, is the decelerating velocity, and is the accelerating velocity. Values of and are found by parameters of , , . The can be seen as the desired damping force of the system. From (32), the control is the semiactive controller for damping force control of MR damper. On the other hand, the observer of the system to estimate the states can be designed byThe values of , are properly chosen by considering both stability and robustness of MR damper control system.

Now, the closed-loop control system for damping force control of MR damper is established as shown in Figure 5. Adaptation laws are used as the outputs of IT2FLS-NT model and the error of system to find adaptive values related function and in Riccati-like equation. The values of output of adaptation laws are input of the fuzzy controller module. The result of the fuzzy controller is input of damper controller, plant, and observer. It is remarked that damper controller supports the plant (MR damper) in control action. The measured force from the output of plant is used as the input for damper controller and IT2FLS-NT modules. The dynamical responses such as displacement and velocity of plant are used for IT2FLS-NT modules and for calculating error system. On the other hand, the output of the observer is also used for calculating error system and the value of error system is used for the input of IT2FLS-NT model.

5. Results and Discussions

5.1. Simulation Results

In order to verify high control performance of the proposed direct adaptive fuzzy controller, in this work both computer simulation and experimental realization are undertaken. The fuzzy rules are shown in Tables 2(a) and 2(b). The centroid of four rule bases of IT2FLS is formed by experiment tests based on trial and error method with considering maximal limit of damping force of MR damper. It is here remarked that the relationship between the centroids and rule bases in Table 2(b) has been clearly analyzed in [25, 26] in which the objective of controlled velocities is equal to or less than the exciting velocity, and the error force is also equal to or less than the limited force. The fuzzy models used in this work are shown in Figures 6(a) and 6(b). There are two parameters for these models: velocity of system and error of damping force. The values of fuzzy sets for velocity of , and error force of , are chosen as follows: In simulation of the proposed controller, the value is chosen as 1 for direct controller and vector . The initial value of is chosen as 10, and the matrix    of Riccati-like equation is chosen as . The constants , are chosen as 20, 10, respectively. The excited frequency is set by 3 Hz. The matrix for observer of the system is used by . From the relation of damping force and control in (32), the stability and robustness of the system will be evaluated by observing tracking performance of actual damping force to the desired force.

In this simulation, the desired force function for the system is chosen as , where is desired magnitude of force set by  N, and is angular velocity where the frequency of the system is set by  Hz. Based on the tracking performance of the damping force, the proposed controller is evaluated in terms of stability and robustness to disturbances. In progress of simulation, the output of control is always evaluated with the desired value of the damping force. This guarantees that the output strictly appears according to the desired value. The external force is chosen as where the frequency is also used as 3 Hz. The initial states of the system are used as for both dynamic and observer states. The values of , are chosen as 20, 20, respectively. The values of , are used as 0.05, 0.05, respectively. The simulation results on the damping force tracking of the proposed controller are shown in Figure 7. Actual force of the system tracks well the desired damping force as shown in Figures 7(a) and 7(b) which show in 1 cycle. The stability of the system after nearly 60 s is shown in Figure 7(c). In this figure, the time of stability is long because the function becomes a constant of the system, and the variable control is applied directly into the system. This progress will take a longer time than indirect adaptive control. The stability is shown in Figure 7(d) in which the velocity of the system becomes small and stable. The applied current of the proposed control in 1 cycle and 100 s is shown in Figures 7(e) and 7(f), respectively. It is seen from Figure 7(a) that the input current is applied in the first time of cycle and holds its value throughout the control process.

In order to highlight good control performance of the proposed controller, a comparative work between the proposed and existing controllers is undertaken. In this work, two existing controllers which are similar to the proposed one are adopted and modified: Tong and Li [1] and Liu and Wang [2]. Adaptation laws in [1, 2] are modified by considering the dynamic model of MR damper which is used in the proposed control system. In simulation for Tong et al. controller, the control vector is used as , value of is set by 10, and the initial states of the system are used by for both dynamic and observer states. The values of , are chosen as 20, 0.05, respectively. The control results of Tong et al. controller are shown in Figure 8. It is seen from Figures 8(a) and 8(b) that the actual force is tracking well to the desired force. But the dynamic parameters such as displacement and velocity are unstable as shown in Figures 8(c) and 8(d). This phenomenon is from the performance of adaptation laws of Tong et al. controller in which all parameters of Riccati-like equation are constant and are not changed through the control process. The applied current of Tong et al. controller is shown in Figures 8(e) and 8(f). It is observed that the results are different from the results shown in Figures 7(e) and 7(f) in which the current input is held and changed a little throughout the control process. The second controller for comparison, Chen et al. controller, is adopted and modified in a same way of Tong et al. controller. It is noted that all setup parameters of Chen et al. controller are similar to Tong et al. controller. The control results of Chen et al. controller are shown in Figure 9. It is observed from results that control results of Chen et al. controller are almost the same as the Tong et al. controller. The applied currents of Tong et al. controller and Chen et al. controller are similar and the current magnitudes are lower than the proposed controller. This directly indicates that the energy consumption of the proposed controller is larger than the others. However, the proposed controller can guarantee the robustness of the stability during control action, while Tong et al. controller and Chen et al. controller are not sufficient to retain the stability of the system due to the lack of the control robustness. From the above observations, the proposed controller provides much better tracking control performance than Tong et al. controller and Chen et al. controller. The proposed controller especially strongly and effectively ensures the robustness against the disturbances resulting in high control performance.

5.2. Experiment Results

In order to more clearly investigate superior control performances of the proposed controller, an experimental apparatus is set up as shown in Figure 10. Force signal is collected by a load cell 20 kgf, and signal force is changed by an amplifier built in dSpace box (DS1104). A wired-type LVDT (Linear Variable Differential Transformer) sensor is used for measuring dynamical response such as displacement and velocity. Exciting vibration is generated by DC motor, and control input from computer system through current amplifier box is applied to MR damper. From (4), the values of , , , , are determined as 941.37 Ns/m, 946.15 Ns/m, 0.073 cm/s, 0.101 cm/s, and 7.436 cm/s, respectively. The control objective in experiment is also to control damping force to follow the desired value.

Control results obtained from experimental realization of the controllers are shown in Figures 11 and 12. The desired force is chosen as 70 N. It is noted that the uncontrolled force (no applied current) and maximal force (maximal applied current) of damper are measured as 21 N and 120 N, respectively. In Figure 11, control results of the proposed controller are shown at four different exciting voltages: 6 V (Voltage), 8 V, 10 V, and 12 V. The exciting frequencies correlative with four voltages are 1.6 Hz, 2.1 Hz, 2.6 Hz, and 3.1 Hz, respectively. It is noted that the exciting voltages indicate the input to be directly applied to the motor. It is seen that controlled force takes nearly 0.04 s to reach the desired force at on-state and 0.08 s at off-state. These times are not changed at four different excited voltages. The positive maximal controlled damping force in Figures 11(a1), 11(b1), 11(c1), and 11(d1) are obtained by 70.998 N, 70.997 N, 70.992 N, and 70.987 N, respectively. The positive minimal controlled damping force in Figures 11(a1), 11(b1), 11(c1), and 11(d1) are 70.098 N, 70.088 N, 70.266 N, and 70.577 N, respectively. The enlarged view of damping forces and control input currents at four excited voltages are shown in Figures 11(a2), 11(a3), 11(b2), 11(b3), 11(c2), 11(c3), 11(d2), and 11(d3). In Figures 11(a2), 11(b2), 11(c2), and 11(d2), it is obviously observed that the stability of the system is not good and hence causes poor tracking performance. It is also seen that the variation of damping force is larger than the proposed controller. From control signals shown in Figures 11(a3), 11(b3), 11(c3), and 11(d3), it is observed that the variation of input current is not large and the maximal control current applied to MR damper is identified by 1.86 A.

For comparison, Tong et al. controller is experimentally realized at same conditions as the controller proposed in this work. The control results are shown in Figure 12. The excited voltages are the same as the proposed controller as 6 V, 8 V, 10 V, and 12 V. The positive maximal controlled damping force in Figures 12(a1), 12(b1), 12(c1), and 12(d1) are obtained by 74.941 N, 73.916 N, 73.964 N, and 73.914 N, respectively. The positive minimal controlled damping force in Figures 12(a1), 12(b1), 12(c1), and 12(d1) are identified by 62.051 N, 63.490 N, 61.758 N, and 62.458 N, respectively. The enlarged view of damping forces and control input currents for four excited voltages are shown in Figures 12(a2), 12(a3), 12(b2), 12(b3), 12(c2), 12(c3), 12(d2), and 12(d3). In Figures 12(a2), 12(b2), 12(c2), and 12(d2), it is obviously observed that the stability of the system is not good, and hence this result in poor tracking performance. It is also seen that the variation of damping force is larger than the proposed controller. This indicates that Tong et al. controller cannot sufficiently control the disturbances and hence the robustness of the system is not fully guaranteed. From the input histories shown in Figures 12(a3), 12(b3), 12(c3), and 12(d3), it is seen that the variation of the control current is not large and the maximal control current is 1.83 A. Figure 13 presents tracking performance of damping force and corresponding applied control input. The positive force and negative force controlled with two different controllers are nearly 27 N and −27 N, respectively. It is clearly seen that the proposed controller is much better than Tong et al. controller in term of tracking performance of damping force. This is because the proposed controller has adaptation law to overcome external disturbances, while Tong et al. controller does not have it.

6. Conclusion

In this work, a new direct adaptive fuzzy controller was developed and its superior control performance was verified through both simulation and experiment associated with MR damper system. The proposed controller was developed based on model of interval type 2 fuzzy and tracking technique. For performance verification, MR damper containing plate-like iron particles was adopted and its damping force tracking was evaluated. In addition, two existing controllers of Tong et al. and Chen et al. were adopted and modified for comparative work with the proposed controller. Control results achieved from simulation and experiment demonstrate that the proposed controller is much better than the others in terms of stability and robustness. This directly indicates the enhancement of damping force tracking performance of the filed-dependent damping force of MR damper of the proposed controller and hence excellent vibration control performance has been achieved. It has been also observed that Tong et al. controller and Chen et al. controller are not converged well and hence cannot sufficiently guarantee the robustness against the disturbance. It is finally remarked that the proposed direct adaptive fuzzy controller can be effectively utilized to numerous MR application control systems subjected to uncertainties or/and disturbances without significant modification.

Conflict of Interests

The authors declare that there is no conflict of interests.

Acknowledgment

This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean Government (MEST) (no. 2010-0015090). This financial support is gratefully acknowledged.