Next Article in Journal
Application of Response Surface Methodology for Optimisation of Simultaneous UHPLC-PDA Determination of Oleanolic and Ursolic Acids and Standardisation of Ericaceae Medicinal Plants
Previous Article in Journal
Potential Application of Fluorescence Imaging for Assessing Fecal Contamination of Soil and Compost Maturity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

L1 Adaptive Control for a Vertical Rotor Orientation System

School of Electrical Engineering, Beijing Jiaotong University, Beijing 100044, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2016, 6(9), 242; https://doi.org/10.3390/app6090242
Submission received: 27 May 2016 / Accepted: 19 August 2016 / Published: 29 August 2016

Abstract

:
Bottom-fixed vertical rotating devices are widely used in industrial and civilian fields. The free upside of the rotor will cause vibration and lead to noise and damage during operation. Meanwhile, parameter uncertainties, nonlinearities and external disturbances will further deteriorate the performance of the rotor. Therefore, in this paper, we present a rotor orientation control system based on an active magnetic bearing with L 1 adaptive control to restrain the influence of the nonlinearity and uncertainty and reduce the vibration amplitude of the vertical rotor. The boundedness and stability of the adaptive system are analyzed via a theoretical derivation. The impact of the adaptive gain is discussed through simulation. An experimental rig based on dSPACE is designed to test the validity of the rotor orientation system. The experimental results show that the relative vibration amplitude of the rotor using the L 1 adaptive controller will be reduced to ∼50% of that in the initial state, which is a 10% greater reduction than can be achieved with the nonadaptive controller. The control approach in this paper is of some significance to solve the orientation control problem in a low-speed vertical rotor with uncertainties and nonlinearities.

Graphical Abstract

1. Introduction

Vertical rotation devices are widely applied in industrial, scientific and civilian fields, e.g., flywheels, centrifuges and washing machines. Many of these applications have a fixed bottom bearing with a free upper side. The one-side bearing structure will lead to vibration during operation and, consequently, cause collision, noise and failure. One of the traditional methods to solve this problem in washing machines is through the addition of a passive balancer [1,2,3]. However, passive control approaches, although having the advantages of simplicity and low cost, are unsatisfactory in vibration reduction because of their open-loop property. For the vertical rotation devices, an active non-contacted vibration reduction system, such as an active magnetic bearing (AMB), placed at the upper side to control the attitude of the rotor, offers a better solution to the vibration problem.
The Active Magnetic Bearing (AMB) system is a nonlinear system, with many uncertainties during its operation, such as parameter uncertainty and external random disturbances. Therefore, the control of the AMB system has been an active research topic. For example, [4] presented a feedback linearization approach to control the position of a horizontal rotor; [5], aiming at a 1-kWh flywheel energy storage device, analyzed the advantages and disadvantages of different approaches, such as decentralized control, Linear Quadratic Regulator (LQR) control and cross-feedback control; [6] proposed an on-line parameter identification method to solve the situation in which the rotor mass is unbalanced; [7], based on an AMB test rig, used radial basis function networks to identify the uncertainties and to synthesize an H controller with the rotor setpoint; [8] presented an adaptive back-stepped controller for a flywheel energy storage system; [9] designed a test rig consisting of a flexible rotor supported by AMBs with a μ synthesis controller.
The research above mostly focused on AMBs for horizontal rotors. Research on a vertical rotor has been rather limited. For example, [10] presented an H controller to regulate the rotor position for vertical AMB systems under 10K rpm. This research, for either horizontal or vertical rotors, has focused on high-speed rotors. However, as mentioned above, vertical low-speed rotating devices are widely used, and their dynamic behaviors will be rather complicated [11]. Therefore, it is necessary to propose a control approach for low-speed vertical rotor position control. Meanwhile, under a lower spinning rate, the influences from parameter uncertainties and external disturbances will be more severe.
The L 1 adaptive control theory [12] provides a powerful tool to overcome the problems of uncertainties and nonlinearities. This theory has been utilized in many applications, e.g., [13] presented a one-output drilling direction adaptive control system in the presence of time-delay, unexpected disturbance and other uncertainties; [14] proposed an underactuated robots control system with unmodeled dynamics. These studies have provided significant theoretical and simulation results for L 1 adaptive controllers. Based on these studies, in this paper, the authors propose an L 1 adaptive controller for a two-input-two-output vertical rotor orientation control system and design an experimental rig based on dSPACE to test the performance of the rotor orientation control system. The work in this paper may also be applied to twin-rotor mechanical system control [15,16,17,18].
This paper is organized as follows. Section 2 establishes the mathematical model of the rotor and presents the control aim of the system. Section 3 proposes the architecture of the adaptive control system. Section 4 analyzes the stability and performance of the adaptive control system. Section 5 gives the simulation of the adaptive control system. Section 6 presents an experimental validation of the performance of the rotor orientation control system. Section 7 provides a summary of the work.

2. Problem Statement

The vertical rotor is shown in Figure 1. α and β are the attitude angles of the rotor from O X Z and O Y Z , respectively. γ is the included angle between the rotor axis and O Z .
According to the theorem of angular momentum, the motion of the rotor satisfies:
d d t L = T
where L is the angular momentum of the rotor and T is the external torque.
The torque induced by the gravity is given by:
T x g = P c sin α cos β T y g = P c sin β T z g = 0
where P = m g is the gravitational force on the rotor; c = O O c . Then, the state equation of a spinning vertical rotor can be written as:
x ˙ 1 = x 3 cos x 2 x ˙ 2 = x 4 x ˙ 3 = P c J x y sin x 1 cos x 2 L J x y x 4 + x 3 x 4 tan x 2 + T x e J x y x ˙ 4 = P c J x y sin x 2 x 3 2 tan x 2 + L J x y x 3 + T x e J x y
where x 1 = α and x 2 = β ; L = J z ω z is the rotor angular momentum; J x y and J z are the rotational inertia about the O X / O Y axis and the O Z axis, respectively; ω z is the spinning rate of the rotor; T x e and T y e are the external control torques in the O X and O Y directions, respectively.
If we define the following variables:
h 1 ( x ) = x 3 cos x 2 h 2 ( x ) = x 4 h 3 ( x ) = P c J x y sin x 1 cos x 2 L J x y x 4 + x 3 x 4 tan x 2 h 4 ( x ) = P c J x y sin x 2 x 3 2 tan x 2 + L J x y x 3
and h ( x ) = h 1 h 2 h 3 h 4 , then Equation (2) can be rewritten as:
x ˙ = h ( x ) + B 1 u
where x = x 1 x 2 x 3 x 4 is the state vector, u = T x e T y e is the system input and:
B 1 = 0 0 0 0 1 J x y 0 0 1 J x y
Obviously, h ( 0 ) = 0 .
Equation (4) is a nonlinear equation. Meanwhile, parameter uncertainties also exist, e.g., the rotational inertia about the O X / O Y axis, J x y , is difficult to calculate accurately, and J z in L = J z ω z , as well.
Without loss of generality, we suppose that J x y = J x y 0 + θ 1 , J z = J z 0 + θ 2 , where J x y 0 and J z 0 are nominal moments of inertia of the rotor and θ 1 and θ 2 are unknown constants referred to the inertia errors. These errors include the manufacturing error and accessory connections. By considering the unknown parameters, the two elements h ( x ) and B in Equation (4) can be written as:
h ( x ) = ( A + A θ ) x + f 1 ( x ) B 1 = B + B θ
where:
A + A θ = h x x = 0
represents the linear parts of h ( x ) , f 1 ( x ) = h ( x ) ( A + A θ ) x represents the high-order nonlinear parts of h ( x ) , while B and B θ represent the certain and uncertain part of B 1 . The subscript θ in Equation (5) refers to the uncertain parameters. A , A θ , B and B θ are:
A = 0 0 1 0 0 0 0 1 P c J x y 0 0 0 J z 0 ω z J x y 0 0 P c J x y 0 J z 0 ω z J x y 0 0
A θ = 0 0 0 0 0 0 0 0 P c θ 1 J x y 0 ( J x y 0 + θ 1 ) 0 0 ω z ( J x y 0 θ 2 J z 0 θ 1 ) J x y 0 ( J x y 0 + θ 1 ) 0 P c θ 1 J x y 0 ( J x y 0 + θ 1 ) ω z ( J x y 0 θ 2 J z 0 θ 1 ) J x y 0 ( J x y 0 + θ 1 ) 0
B = 0 0 0 0 1 J x y 0 0 0 1 J x y 0 , B θ = 0 0 0 0 θ 1 J x y 0 ( J x y 0 + θ 1 ) 0 0 θ 1 J x y 0 ( J x y 0 + θ 1 )
Next, we define:
f ( t , x ( t ) ) = f 1 ( x ) + A θ x ( t ) + B θ u ( t )
and
u = u m + u ad
where u m = K m x ( t ) is a linear state feedback that renders A m A B K m Hurwitz. u ad is the adaptive input vector used to compensate for the unknown parameter, system nonlinearity and other disturbances. Then, Equation (4) can be rewritten as:
x ˙ ( t ) = A x ( t ) + B ( u m + u ad ) + f ( t , x ( t ) ) x ( 0 ) = x 0
f ( t , x ) should satisfies the following two assumptions:
Assumption 1. 
(Uniform boundedness of f ( t , 0 ) ).
For any t 0 , there exists F > 0 , such that f ( t , 0 ) F .
Assumption 2. 
(Boundedness of partial derivatives).
For arbitrary δ > 0 , there exist d f x ( δ ) and d f t ( δ ) , such that, for arbitrary x , as long as x δ , the partial derivatives of f ( t , x ) are bounded:
f x d f x ( δ ) , f t 1 d f t ( δ )
In this paper, the system equilibrium x = 0 implies that the rotor is along the vertical direction, which is just the control aim of the system. In other words, the control input must ensure that the equilibrium of Equation (8) is x = 0 . While the rotor lies vertically, the external control torques also need to be zero, i.e., ( u m + u ad ) x = 0 = 0 . Therefore, f ( t , 0 ) = 0 for all t 0 , and Assumption 1 is satisfied naturally. Assumption 2 requires that f ( t , x ) is continuous about x and slowly-varying about t. f ( t , x ) consists of three parts: f 1 ( x ) , A θ x and B θ u ( t ) . It is obvious that f 1 and A θ x satisfy the requirement, and the key issue is u ( t ) . Equation (7) indicates that u ( t ) contains u m and u ad , while u m = K m x satisfies the requirement, and u ad is an artificial adaptive control input, which can be designed to be continuous and slowly varying. Therefore, Assumption 2 can be satisfied.
The control aim is to seek a proper adaptive control signal u ad , such that the state in Equation (2) is stable about the origin.

3. L 1 Adaptive Control Design of the Rotor Orientation System

3.1. Semi-Global Linearization of the System

If the nonlinearity f is subject to the two assumptions in Section 2, then, for all x τ L δ , f ( t , x ( t ) ) can be linearized to a time-varying equation about [12], the norm of x [12], i.e.,
f ( t , x ( t ) ) = K θ ( t ) x ( t ) L + σ ( t )
where K θ ( t ) Θ R 2 × 1 , Θ is a convex set and K θ ( t ) d f x ( δ ) , σ F . d f x ( δ ) and F are defined in Assumptions 1 and 2, respectively. In addition, there exists a constant d θ > 0 , such that K ˙ θ d θ .
In this paper, as f ( t , 0 ) = 0 , σ ( t ) = 0 , and Equation (8) can be written as a linear time varying system:
x ˙ = A m x + B u ad + K θ ( t ) x ( t ) L x ( 0 ) = x 0

3.2. State Predictor

For the linearly-parameterized system Equation (11), consider the following state predictor:
x ^ ˙ ( t ) = A m x ^ ( t ) + B K ^ θ ( t ) x ( t ) + u ad x ^ ( 0 ) = x 0
where x ^ ( t ) R 4 × 1 is the state of the predictor, K ^ θ = K ^ θ 1 ( t ) K ^ θ 2 ( t ) is the estimate of K θ ( t ) and K ^ θ ( 0 ) Θ .

3.3. L 1 Adaptive Control Algorithm

To streamline the following derivation, we need to use the projection operator Proj ( · , · ) , which is defined in [12]. The projection operator has the following property: For any y R n , θ * Ω 0 Ω c and θ Ω c , we have:
θ θ * Proj ( θ , y ) y 0
With the projection operator, the L 1 adaptive control algorithm can be expressed as:
K ^ ˙ θ = Γ Proj K ^ θ , B P e x ( t ) L , K ^ θ ( 0 ) = K ^ θ 0 Θ
where Γ > 0 is the adaptive gain and P = P > 0 is a positive matrix that solves the Lyapunov equation A m P + PA m = Q for any Q = Q > 0 . e ( t ) = x ^ ( t ) x ( t ) is the prediction error. The projection operator ensures that K ^ θ Θ for all t 0 .
The adaptive control signal is:
u ad ( s ) = C ( s ) ( η ^ ( s ) + r g )
where C ( s ) is a low-pass filter, η ^ ( s ) is the Laplace transform of η ^ ( t ) = K ^ θ x ( t ) and r g is the reference input signal. In this paper, the control aim is x = 0 at u = 0 , so r g = 0 0 .
The block diagram of the adaptive control is illustrated in Figure 2.

4. Adaptive Control System Analysis

4.1. Boundedness of State Error and Parameter Error

To evaluate the performance of the adaptive control algorithm, first, we need to analyze the error boundedness. Consider the following prediction error equation derived from Equations (11) and (12):
e ˙ ( t ) = A m e ( t ) + B θ e ( t ) x ( t ) L e ( 0 ) = 0
where θ e ( t ) = K ^ θ ( t ) K θ ( t ) . For the stability of Equation (16), we have the following proposition:
Proposition 1. 
If the adaptive input u ad is defined as Equation (15), then we have the following result:
e L θ m λ min ( P ) Γ
where:
θ m 4 max ϑ Θ ϑ 2 + 4 λ max ( P ) λ min ( Q ) d θ max ϑ Θ ϑ
Proof. 
First, we consider the following positive definite function:
V e ( t ) , θ e ( t ) = e ( t ) P e ( t ) + 1 Γ θ e ( t ) θ e ( t )
Obviously,
V ( e ( 0 ) , θ e ( 0 ) ) = 1 Γ θ e ( 0 ) θ e ( 0 ) 1 Γ max ϑ Θ ϑ 2 θ m Γ
The derivative of V about t is:
V ˙ ( e , θ e ) = e ˙ P e + e P e ˙ + 2 Γ θ e K ^ ˙ θ 2 Γ θ e K ˙ θ = e Q e + 2 θ e B P e x L + Proj ( K ^ θ , B P e x L ) 2 Γ θ e K ˙ θ
According to the property of the projection operator, it follows that:
2 θ e B P e x L + Proj ( K ^ θ , B P e x L ) 0
Then, the derivative yields:
V ˙ ( e , θ e ) e Q e 2 Γ θ e K ˙ θ e Q e + 2 Γ θ e K ˙ θ
Given that K ^ θ Θ , K θ Θ and K ˙ θ d θ , for any τ 0 , the parameter error yields:
θ e K ˙ θ θ e · K ˙ θ = K ^ θ K θ · K ˙ θ 2 max ϑ Θ ϑ · d θ 1 Γ θ e ( τ ) θ e ( τ ) 4 Γ max ϑ Θ ϑ 2
Therefore, if, for any τ, there exists V ( τ ) > θ m / Γ , the quadratic form of the prediction error has to yield:
e ( τ ) P e ( τ ) > 4 λ max ( P ) Γ λ min ( Q ) max ϑ Θ ϑ d θ
Hence:
e ( τ ) Q e ( τ ) λ min ( Q ) e ( τ ) Q e ( τ ) λ max ( P ) > 4 Γ max ϑ Θ ϑ
which means that, for any τ > 0 , as long as V ( τ ) > θ m / Γ , its derivative yields.
V ˙ ( t , x ( t ) ) 0
Given that V ( t , 0 ) θ m / Γ , therefore, for all t 0 , it follows that:
V ( t , x ( t ) ) θ m Γ
Meanwhile, as λ min ( P ) e ( t ) L 2 e ( t ) P e ( t ) V ( t ) , it follows that:
e ( t ) L 2 θ m λ min ( P ) Γ , t 0
Then Equation (17) is proved. ☐
According to Equation (18), the parameter error yields:
1 Γ θ e θ e = V e P e v θ m Γ
which leads to:
θ e L θ m
Therefore, the parameter error is also bounded.

4.2. Boundedness of State Variables

To analyze the boundedness of the system state x ( t ) , we first consider the following ideal system:
x ˙ id ( t ) = A m x id ( t ) + B K θ x id ( t ) + u id ( t ) x id ( 0 ) = x 0
The ideal control signal is:
u id ( s ) = C ( s ) η ( s )
where η ( s ) is the Laplace transform of η ( t ) = K θ ( t ) x id ( t ) .
We define G ( s ) ( s I 4 A m ) 1 B ( s I 2 C ( s ) ) and L max ϑ Θ ϑ 1 , where I 4 and I 2 are, respectively, fourth-order and second-order identity matrices. According to the definition of K θ in Equation (10), K θ L L . It can be proven that x id in Equation (31) is banded input banded state (BIBS) stable when G ( s ) L 1 L < 1 .
Proof. 
According to the definition of G ( s ) , Equation (31) can be transformed to:
s x id ( s ) x ( 0 ) = A m x id ( s ) + B ( s I 2 C ( s ) ) η ( s )
After performing some algebraic operations, we get:
x id ( s ) = G ( s ) η ( s ) + x in ( s )
where x in ( s ) = ( s I 4 A m ) 1 x 0 . For τ ( 0 , ) and t [ 0 , τ ] , the following bound exists [12]:
x id τ L G ( s ) K θ L 1 · x id ( t ) τ L + x in τ L
According to the definition of the ∞-norm and the L 1 -norm, the following bound exists:
x id τ L x id τ ( t ) L
Meanwhile, since G ( s ) K θ L 1 G ( s ) L 1 K θ G ( s ) L 1 L < 1 , therefore,
x id τ L x in L 1 G ( s ) K θ L 1
According to the definition of x in , if it is bounded, x id will be bounded. Therefore, System (31) is BIBS stable.
The system in Equation (31) is equivalent to the system predictor in Equation (12) while K ^ θ = K θ . Because θ e and x in are both bounded, x ^ is bounded. Meanwhile, since e = x ^ x and e is also bounded, therefore, x is bounded, i.e., System (2) is BIBS stable under the L 1 -controller. ☐

4.3. Performance of the Adaptive System

According to Equation (17), the prediction error e can be reduced by enhancing λ min ( P ) or the adaptive gain Γ. λ min ( P ) depends on the performance of the linear part A m , i.e., the wider the stability margin of the linear part of the system is, the smaller will be the prediction error. Similarly, the increase of Γ will reduce the error, as well.

5. Simulation Results for the Adaptive Control System

5.1. Simulation of the Adaptive Controller

Consider the nominal mechanical parameters in Table 1.
According to Table 1, the nominal rotational inertias of the rotor about the two axes are: J x y 0 = 0 . 01864 kg · m 2 and J z 0 = 0 . 008439 kg · m 2 . Then, the system parameters in Equation (8) are:
A = 0 0 1 0 0 0 0 1 161 . 4 0 0 185 . 1 0 161 . 4 185 . 1 0 , B = 0 0 0 0 56 . 65 0 0 56 . 65
The control input is:
u m = K m x = 0 100 7 0 100 0 0 7 x
Therefore,
A m = A B K m = 0 0 1 0 0 0 0 1 161 . 4 5665 396 . 55 185 . 1 5665 161 . 4 185 . 1 396 . 55
Let the parameter uncertainties be θ 1 = 0 . 2 J x y 0 and θ 1 = 0 . 2 J z 0 , the adaptive gain be Γ = 20 , 000 and the low-pass filter be:
C ( s ) = 3 ω c 2 s + ω c 3 ( s + ω c ) 3 0 0 3 ω c 2 s + ω c 3 ( s + ω c ) 3
The bandwidth of C ( s ) will influence the stability and performance of the system. A detailed discussion of the low-pass filter C ( s ) is in [12].
The simulation model using MATLAB/Simulink is illustrated in Figure 3. (a) is the parameter adaptive law in Equation (14); (b) is the nonadaptive control input u m in Equation (39); (c) is the low-pass filter C ( s ) in Equation (41); (d) is the reference input, which is zero in this paper’s situation; (e) is the state predictor in Equation (12); (f) is the plant system in Equation (2); (g) gives the parameter uncertainties; (h) is the external disturbance.
The simulation results for the prediction error in every component are shown in Figure 4. All of the components of the prediction error get reduced as time progresses.
The variation trends of system parameter prediction, K ^ θ , and one of the prediction error components, e 1 , under different adaptive gains are shown in Figure 5. It can be seen that, when Γ is larger, K ^ θ and e 1 will converge more quickly. Therefore, Γ should be as large as possible, as long as the actuator of the control system can be achieved.

5.2. Simulation of the Rotor Orientation System

Figure 4 and Figure 5 imply that the adaptive control algorithm in Equations (14) and (15) is effective. On this basis, to further simulate the performance of the orientation system, we introduce more uncertainties. Besides the parameter uncertainty, consider the following uncertain system:
x ˙ = A m x + B ( u ad + K θ x + σ ) , x ( 0 ) = x 0
where σ refers to an external random disturbance. Unlike Equation (11), here σ 0 .
Let the adaptive controller activate at t = 4 s , then we obtain the result of the rotor angular drift, γ d = α 2 + β 2 , from the vertical direction, as shown in Figure 6. It is obvious that, after the controller actuates, the angular displacement of the rotor reduces from more than > 1 to < 0 . 4 . Therefore, the orientation control is effective.

6. Experimental Results for the Rotor Orientation System

6.1. Experimental Devices

To demonstrate the feasibility of the adaptive control algorithm, we designed an experimental platform as shown in Figure 7. The iron rotor lies vertically in the center. It is connected to the drive motor at the bottom with a universal bearing, which enables the rotor to swing freely while rotating with the motor coaxially. Four electromagnets lie perpendicularly around the rotor to provide the orientation control forces. Two gap sensors are installed at a 90 interval to detect the attitude of the rotor and provide feedback signals to the control system.

6.2. Experimental Process

The block diagram is illustrated in Figure 8. In the experiment, we choose dSPACE as the digital controller.
The experimental program is as follows:
1.
Measurement of the rotor drift angle.
The drift angles of the rotor in the X / Y directions a n g l e x and a n g l e y , respectively, are converted from the two gap values g a p 1 and g a p 2 . The two angles are the state variables x 1 and x 2 in Equation (2).
2.
Calculation of the electromagnetic force.
According to Equations (15) and (39), the control torque is:
u = T x e T y e = u m + u ad
The torque is given by:
T x e = ( F y F y + ) h m T y e = ( F x + F x ) h m
where F x + , F x , F y + and F y are the forces of each electromagnet in the X and Y directions and h m is the installation height of the electromagnets.
Because the electromagnetic force can only be attractive, to make them continuous and smooth, the forces of each direction are realized via a bias control, i.e.,
F x + = F n 0 + M y e 2 h m F x = F n 0 M y e 2 h m F y + = F n 0 M x e 2 h m F x = F n 0 + M x e 2 h m
where F n 0 is the reference value of each electromagnet, about which each force varies. The control force signal for each electromagnet is determined through Equations (43) and (45).
3.
Calculation of the current in every electromagnet.
The electromagnetic force of each electromagnet yields:
F = A I 2 ( B + g ) 2
where A and B are factors corresponding to the electromagnet. Substituting Equations (46) into the control force Equations (45) gives the current signal of each electromagnet.
I x + = ( B + g x + ) F n 0 + M x c 2 h m / A I x = ( B + g x ) F n 0 M x c 2 h m / A I y + = ( B + g y + ) F n 0 + M y c 2 h m / A I y = ( B + g y ) F n 0 M y c 2 h m / A
where g x ± and g y ± are the gaps of each electromagnet. The current signal in Equation (47) can be loaded onto each electromagnet via a power amplifier.

6.3. Experimental Results

The experimental results for the orientation control system are shown in Figure 9. To clearly show the performance of the adaptive orientation control system, a result of the nonadaptive orientation control, that u ad = 0 and only u m activates, is also illustrated, as a contrast. γ d 0 is the average angular displacement without control, which is about 0 . 8 of the adaptive control and 0 . 7 of the nonadaptive control. The control system starts to actuate at t = 4 s .
It can be seen from Figure 5 that, after the controller activates, the average angular displacement of the rotor from the vertical direction reduces to less than ∼50% of the adaptive controller and about ∼60% of the nonadaptive controller. The adaptive controller has a better performance in vibration reduction.
Furthermore, by comparing Figure 9 with Figure 6, it can be seen that the trend of the angular displacement is basically the same.

7. Conclusions

In the paper, an adaptive orientation control system for a vertical low-speed rotor is presented. The nonlinearity and parameter uncertainty are transformed into an equivalent linear time-varying form about the norm of the system state. Then, an adaptive algorithm is established based on the state predictor, and an adaptive position control method is presented to overcome the parameter uncertainty, nonlinearity and random disturbance. The stability and boundedness of the prediction error has been analyzed. The analysis indicates that the prediction error is uniformly bounded and that the bound is inversely proportional to the square root of the adaptive gain. Meanwhile, the state of the adaptive control system is BIBS stable. Finally, the control system is verified via simulation and experiment. Simulation results demonstrate the boundedness of the prediction error and the parameter estimation. The experimental result shows that the position control method reduces the amplitude of rotor vibration effectively during rotation. The work in this paper is of significance for control methods and the design of vertical AMB systems and similar rotor mechanical systems.

Acknowledgments

The work in this paper was supported by the National Science Fund of China (Grant Number 51077003); and the Fundamental Research Funds for the Central Universities of China (Grant Number 2016YJS143).

Author Contributions

Sijia Liu and Yu Fan conceived of and designed the control algorithm. Sijia Liu and Jun Di designed the experiments. Sijia Liu and Mingming Ji performed the experiments and analyzed the data. Sijia Liu wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bae, S.; Lee, J.; Kang, Y.; Kang, J.; Yun, J. Dynamic analysis of an automatic washing machine with a hydraulic balancer. J. Sound Vib. 2002, 257, 3–18. [Google Scholar] [CrossRef]
  2. Chen, H.W.; Zhang, Q.J. Stability analyses of a vertical axis automatic washing machine with a hydraulic balancer. Mech. Mach. Theory 2011, 46, 910–926. [Google Scholar] [CrossRef]
  3. Chen, H.W.; Zhang, Q.J.; Fan, S.Y. Study on steady-state response of a vertical axis automatic washing machine with a hydraulic balancer using a new approach and a method for getting a smaller deflection angle. J. Sound Vib. 2011, 330, 2017–2030. [Google Scholar] [CrossRef]
  4. Smith, R.D.; Weldon, W.F. Nonlinear control of a rigid rotor magnetic bearing system: Modeling and simulation with full state feedback. IEEE Trans. Magn. 1995, 31, 973–980. [Google Scholar] [CrossRef]
  5. Ahrens, M.; KuCera, L.; Larsonneur, R. Performance of a Magnetically Suspended Flywheel Energy Storage Device. IEEE Trans. Control Syst. Technol. 1996, 4, 494–502. [Google Scholar] [CrossRef]
  6. Lum, K.Y.; Coppola, V.; Bernstein, D. Adaptive autocentering control for an active magnetic bearing supporting a rotor with unknown mass imbalance. IEEE Trans. Control Syst. Technol. 1996, 4, 587–597. [Google Scholar]
  7. Gibson, N.S.; Choi, H.; Buckner, G.D. H control of active magnetic bearings using artificial neural network identification of uncertainty. In Proceedings of the IEEE International Conference on Systems, Man and Cybernetics, Washington, DC, USA, 5–8 October 2003; Volume 2, pp. 1449–1456.
  8. Sivrioglu, S. Adaptive backstepping for switching control active magnetic bearing system with vibrating base. Control Theory Appl. IET 2007, 1, 1054–1059. [Google Scholar] [CrossRef]
  9. Mushi, S.E.; Lin, Z.; Allaire, P.E. Design, construction, and modeling of a flexible rotor active magnetic bearing test rig. IEEE/ASME Trans. Mechatron. 2012, 17, 1170–1182. [Google Scholar] [CrossRef]
  10. Tsai, N.; Kuo, C.; Lee, R. Regulation on radial position deviation for vertical AMB systems. Mech. Syst. Signal Process. 2007, 21, 2777–2793. [Google Scholar] [CrossRef]
  11. Shi, M.; Wang, D.; Zhang, J. Nonlinear dynamic analysis of a vertical rotor-bearing system. J. Mech. Sci. Technol. 2013, 27, 9–19. [Google Scholar] [CrossRef]
  12. Hovakimyan, N.; Cao, C. L1 Adaptive Control Theory: Guaranteed Robustness with Fast Adaptation; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2010; Volume 21. [Google Scholar]
  13. Sun, H.; Li, Z.; Hovakimyan, N.; Başsar, T.; Downton, G. L1 Adaptive Control for Directional Drilling Systems. IFAC Proc. Vol. 2012, 45, 72–77. [Google Scholar] [CrossRef]
  14. Nguyen, K.D.; Dankowicz, H. Adaptive control of underactuated robots with unmodeled dynamics. Robot. Auton. Syst. 2014, 64, 84–99. [Google Scholar] [CrossRef]
  15. Tao, C.W.; Taur, J.S.; Chang, Y.H.; Chang, C.W. A novel fuzzy-sliding and fuzzy-integral-sliding controller for the twin-rotor multi-input–multi-output system. IEEE Trans. Fuzzy Syst. 2010, 18, 893–905. [Google Scholar] [CrossRef]
  16. Allouani, F.; Boukhetala, D.; Boudjema, F. Particle swarm optimization based fuzzy sliding mode controller for the twin rotor MIMO system. In Proceedings of the 2012 16th IEEE Mediterranean Electrotechnical Conference, Hammamet, Tunisia, 25–28 March 2012; pp. 1063–1066.
  17. Nejjari, F.; Rotondo, D.; Puig, V.; Innocenti, M. Quasi-LPV modelling and non-linear identification of a twin rotor system. In Proceedings of the 2012 20th Mediterranean Conference on Control & Automation (MED), Barcelona, Spain, 3–6 July 2012; pp. 229–234.
  18. Rădac, M.B.; Precup, R.E.; Petriu, E.M.; Preitl, S. Iterative data-driven tuning of controllers for nonlinear systems with constraints. IEEE Trans. Ind. Electron. 2014, 61, 6360–6368. [Google Scholar] [CrossRef]
Figure 1. Rotor scheme.
Figure 1. Rotor scheme.
Applsci 06 00242 g001
Figure 2. L 1 adaptive control system algorithm.
Figure 2. L 1 adaptive control system algorithm.
Applsci 06 00242 g002
Figure 3. Block diagram in Simulink.
Figure 3. Block diagram in Simulink.
Applsci 06 00242 g003
Figure 4. Prediction error simulation results for: (a) the x-direction component e 1 ; (b) the y-direction component e 2 ; (c) the x-direction angular error e 3 ; and (d) y-direction angular error e 4 .
Figure 4. Prediction error simulation results for: (a) the x-direction component e 1 ; (b) the y-direction component e 2 ; (c) the x-direction angular error e 3 ; and (d) y-direction angular error e 4 .
Applsci 06 00242 g004
Figure 5. Simulation comparison for Γ values: (a) comparison of parameter K θ ; (b) comparison of error e 1 .
Figure 5. Simulation comparison for Γ values: (a) comparison of parameter K θ ; (b) comparison of error e 1 .
Applsci 06 00242 g005
Figure 6. Simulation result for the orientation control system.
Figure 6. Simulation result for the orientation control system.
Applsci 06 00242 g006
Figure 7. Experimental platform for the orientation control system: (a) electromagnet; (b) rotor; (c) universal joint; (d) motor.
Figure 7. Experimental platform for the orientation control system: (a) electromagnet; (b) rotor; (c) universal joint; (d) motor.
Applsci 06 00242 g007
Figure 8. Block diagram of the orientation experiment: (a) gap sensor; (b) electromagnet; (c) rotor; (d) gap calculation; (e) angle calculation; (f) controller; (g) control current calculation; (h) power amplifier.
Figure 8. Block diagram of the orientation experiment: (a) gap sensor; (b) electromagnet; (c) rotor; (d) gap calculation; (e) angle calculation; (f) controller; (g) control current calculation; (h) power amplifier.
Applsci 06 00242 g008
Figure 9. Experimental results for the angular drift of the rotor orientation system.
Figure 9. Experimental results for the angular drift of the rotor orientation system.
Applsci 06 00242 g009
Table 1. Rotor parameters.
Table 1. Rotor parameters.
ParameterValueUnit
MaterialIron
Outer diameter150mm
Inner diameter142mm
Height150mm
Mass1.737kg

Share and Cite

MDPI and ACS Style

Liu, S.; Fan, Y.; Di, J.; Ji, M. L1 Adaptive Control for a Vertical Rotor Orientation System. Appl. Sci. 2016, 6, 242. https://doi.org/10.3390/app6090242

AMA Style

Liu S, Fan Y, Di J, Ji M. L1 Adaptive Control for a Vertical Rotor Orientation System. Applied Sciences. 2016; 6(9):242. https://doi.org/10.3390/app6090242

Chicago/Turabian Style

Liu, Sijia, Yu Fan, Jun Di, and Mingming Ji. 2016. "L1 Adaptive Control for a Vertical Rotor Orientation System" Applied Sciences 6, no. 9: 242. https://doi.org/10.3390/app6090242

APA Style

Liu, S., Fan, Y., Di, J., & Ji, M. (2016). L1 Adaptive Control for a Vertical Rotor Orientation System. Applied Sciences, 6(9), 242. https://doi.org/10.3390/app6090242

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop