Next Article in Journal
User-Configurable Timing and Navigation for UAVs
Previous Article in Journal
On the Potential of the RST-FLARE Algorithm for Gas Flaring Characterization from Space
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Position Estimation and Compensation Based on a Two-Step Extended Sliding-Mode Observer for a MSFESS

1
Beijing University of Aeronautics and Astronautics, Beijing 100191, China
2
PowerChina Road Bridge Group Company Limited, Beijing 100048, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(8), 2467; https://doi.org/10.3390/s18082467
Submission received: 24 June 2018 / Revised: 12 July 2018 / Accepted: 25 July 2018 / Published: 30 July 2018
(This article belongs to the Section Intelligent Sensors)

Abstract

:
This paper aims to deal with the problem of rotor position estimation and compensation for a magnetically suspended flywheel energy storage system under the consideration of measurement noise and unknown disturbances. First, the flywheel system working principle and description are analyzed and, based on this, the mathematical model as well as the coordinates transformation are introduced. For the purpose of the state estimation, a two-step extended sliding-mode observer is considered to obtain the estimates of the rotor angular position. In this control strategy, a traditional sliding-mode observer is adopted as a first-step original state estimator. After that, the relationship between the angular position and the estimation error is established and a second-step observer is designed to obtain the estimation of the error. The estimated error is then used to compensate the real values of the rotor angular position generated by the first-step observer. To reject the influences of the measurement noise and unknown disturbances, the H optimization strategy is considered to determine the second-step observer structure. Finally, experimental results are presented to demonstrate the effectiveness of the proposed method. It is demonstrated that the proposed two-step observer method has a better estimation accuracy and control performance.

1. Introduction

As a novel physical energy storage technology, the magnetically suspended flywheel energy storage system (MSFESS) has the advantages of rapid charging-discharging process, high density and efficiency of energy and power, long life-period with small environment pollution and etc. [1,2,3]. Due to these priorities, the MSFESS has been widely used in the periods of rail traffic energy recovery, electric power magnitude regulation, uninterruptible power supply (UPS), electromagnetic launch and so on. Probably, the MSFESS contains two working modes: charging mode and discharging mode, which can be guaranteed with the active magnetic bearing (AMB), vacuum equipment, flywheel rotor and the motor/generator stator. In the charging mode, the flywheel rotor is driven by the motor stator for a high-speed rotation so that the electric energy can be stored in the form of kinetic energy. In the discharging mode, the motor stator is actuated by the flywheel rotor so that the kinetic energy can be transferred to the electric energy used by the load, in which the flywheel rotor has to decrease the speed gradually. Owing to the advantages of simple structure, small mass and volume, high efficiency and power factor and etc., the permanent magnet synchronous motor (PMSM) is the ideal selection as the driven device for the high-speed flywheel rotor to accomplish the charging-discharging control process [4,5].
In real practice, the vector control strategy is commonly used in PMSM control to improve the speed and current performance, in which the accurate rotor angular position is needed for coordinate transformation. The traditional method to obtain the rotor position is the sensor detection strategy such as resolver, encoder and hall device, which can be the simplest efficient measure by considering the real operation circumstance. But for magnetically suspended devices like MSFESS, twyer and compressor, the rotor is suspended in the seal house and there is a certain gap between the rotor and the body. In this case, a large displacement variation occurs in both axial and radial directions when the rotor begins to suspend or stays static. In addition, it is a common phenomenon that the rotor displacement can vary in real time, which requires a more difficult installation accuracy for the position sensor and a higher cost and reducer reliability. Therefore, the sensorless control strategy has become a leading issue in the PMSM control field [6,7,8,9].
In general, several studies pay attention to the sensorless control strategy of the PMSM: the first method relies on the magnetic path asymmetry of the permanent magnet motor. In this method, additional voltage and current signals are needed to produce a high-frequency excitation, in which the rotor position can be obtained by measuring the system response. This kind of method includes high-frequency signal injection, motor inductance detection and etc. [10,11,12,13]. However, the high-frequency signal excitation method is appropriate for the salient pole PMSM, which has a high-precision requirement for the physical and magnetic structure. The second one is the model estimation method such as the back electromotive force integration method, model reference self-adaption, sliding mode observer, extended Kalman filter and etc. [14,15,16]. The back electromotive force integration method has a wide application due to its simplicity and reliability, but this method relies on the motor stator voltage and current information to calculate the vector phase. When the current is small or zero, the rotor position cannot be obtained, resulting in a step failing out. The model reference self-adaption method contains an adjustment or reference model and a self-adaption scheme, which has a better capability of disturbance rejection, but this method depends on an accurate mathematical model and easily results in estimation error. The extended Kalman filter is an optimal full-order estimator, which can simultaneously reject the measurement noise and unknown disturbance, but this filter requires a complicated algorithm calculation and needs to be done with a powerful processor.
Among the mode-based estimation methods, the sliding-mode observer method has received more and more attention [17,18,19]. The observer designs the current error as the sliding-mode surface and based on this, the states can be converted to the real values by establishing an appropriate sliding-mode function. Due to the robustness regarding model parameters, the sliding-mode observer has been applied in many real experimental systems. However, a low-pass filter has to be used in this observer since there is a certain noise in the current information, resulting in an estimation phase delay and accuracy decline that need to be compensated for real time. To solve this problem, several modified observer strategies have been proposed to improve the estimation accuracy, such as extended sliding-mode observer and robust non-linear observer and etc., but these observers have not disposed the phase delay radically so that the estimation performance can be affected. In a previous study [20], a kind of phase-locked loop-based sliding-mode observer was used, which can improve the state estimation accuracy in place, but the control scheme is difficult to realize in practice. Above all, for the sensorless control strategy of the PMSM, the first physical estimation strategy can be fixable for the high-speed condition, and the estimation performance cannot be guaranteed, while the observer strategy has a certain phase delay and estimation error. In the whole speed areas, the sensorless control strategy is still a research emphasis, and there is difficulty for PMSM high-precision control.
Motivated by this, this paper proposes a two-step extended sliding-mode observer-based control strategy. The first-step observer adapts a traditional extended sliding-mode observer with a phase compensation to obtain the original estimated states, and the second-step establishes the error model and uses the H optimization strategy to reject the influences of unknown disturbances. To the authors’ knowledge, the H optimization control method has a priority and robustness against unknown disturbances [21,22,23]. Therefore, this paper aims to apply this method to the rotor position estimation and compensation for the MSFESS. This paper is organized as follows. Firstly, in Section 2, the system model and description are introduced. Section 3 gives the basic idea of this work, with the design of the sliding-mode observer proposed in Section 4. Experimental setup and results are presented in Section 5. Finally, Section 6 concludes this paper.

2. System Modeling and Description

As shown in Figure 1, the MEFESS is composed of the seal house, PMSM, AMB, electric pole mechanical bearing, etc. The PMSM is the crucial component for energy transmission in the form of electric and kinetic, which is divided into a motor stator and coaxial permanent magnet. The AMB possesses a 5-degree freedom in the axial and radial directions, which is able to maintain the high-speed rotor suspending in the vacuum body. The high-speed rotor is used as the kinetic energy storage component, and the mechanical bearing provides the strength protection once the AMB is running out. Outside of them, the seal house is reinforced to fix all of the components and thus provides a vacuum.
Considering the problem of rotor strength during high-speed rotation, the PMSM adopts a mechanism of a surface-mounted permanent magnet with a high-intensity sheath to fasten it. The minimum power scheme of the PMSM driver is shown in Figure 2. The six insulated gate bipolar transistors (IGBTs) VT1~VT6 are used as the power switches of the rectifier and inverter. The resistor R0 and relay KT compose the power soft-start module, while the resistors R1~R3 and capacitors C1~C3 are able to stabilize the DC-bus voltage Udc. The resistor RL is traded as the simulate load. In the three-phase interface of the motor electric pole, three current sensors U/V/W are connected to measure the stator current.
Based on this, the kinematic model for the surface-mounted PMSM (SPMSM) in the coordinate of three-phase static a/b/c can be described as follows:
{ u s = R s i s + L s d i s d t + d ψ s d t , φ s = φ f [ sin ( ω r t + θ r ) sin ( ω r t + θ r + 2 π / 3 ) sin ( ω r t + θ r + 4 π / 3 ) , ]  
where u s = [ u a u b u c ] T , i s = [ i a i b i c ] T are stator three-phase voltage and current, respectively. L s , R s are armature inductance and resistor. θ r , ω r are angular position and speed.
Generally, the PID-based vector control strategy was adopted for the SPMSM system, including speed loop during charging (voltage loop during discharging), double current loop and position loop, as shown in Figure 3. V r e f , ω r e f , I q r e f and I d r e f are voltage, speed and current reference inputs, respectively. G v ( s ) , G ω ( s ) , G i q ( s ) and G i d ( s ) are PID controllers.
In the vector control loop as shown in Figure 3, the coordinate transformations need to be disposed, including CLARKE transformation (three-phase a/b/c to two-phase static α/β coordinate, 3s/2s), PARK transformation (two-phase static α/β to two-phase rotating d/q coordinate, 2s/2d) with its inverse IPARK transformation. The kinematic model for the SPMSM in the coordinate of two-phase static α/β can be described as follows:
{ u α = R s i α + L s d i α d t + e α , u β = R s i β + L s d i β d t + e β , e α = ω r φ f sin θ r , e β = ω r φ f cos θ r ,  
where u α , u β , i α , i β are motor stator two-phase static voltage and current, respectively. e α , e β are motor back electromotive force in the α/β coordinate. φ f is the excitation flux linkage of the rotor permanent magnet.
According to the IPARK transformation, the kinematic model in the coordinate of two-phase rotating d/q can be formulated as:
{ u d = R s i d + L s d i d d t ω r L s i q , u q = R s i q + L s d i q d t + ω r L s i d ,  
where u d , u q , i d , i q are motor stator two-phase rotating voltage and current, respectively. In this formulation, the variation of the permanent magnet can be ignored.
It is known that the rotor angular position θ r is the crucial parameter in the coordinates transformation. Usually, the position sensor resolver or encoder is applied to obtain the angular position and speed. However, for the MSFESS, the rotor is suspended in the geometric center and hence there is a certain large distance between the motor stator and rotor. Accordingly, there are nearly no position sensors to satisfy this distance circumstance. In addition, the special sensors are usually expensive, which cannot be tolerated due to economic consideration. Therefore, the sensorless vector control becomes more and more important in practical engineering fields. In this paper, we focus on the real-time position estimation and compensation for the SPMSM based on the two-step extended slide-mode observer.

3. Basic Idea of This Work

The basic idea of this work is to estimate and compensate the angular position for the SPMSM by using the sensorless control strategy. Generally, the sensorless estimation approaches can be summarized as back electromotive force method, high-frequency signal injection, observers, etc. Among these methods, it has been shown that the observer-based sensorless estimation (SLE) is proved to be one of the most effective methods. For this purpose, the scheme of the traditional extended sliding-mode observer is first given as follows.
According to the relationships in Equation (2), the angular position and speed can be formulated as:
{ θ r = a r c tan e α e β ω r = d θ r / d t , ,  
where it can be seen that the value of the position signal θ r depends on the back electromotive forces e α and e β in the α/β coordinate.
On the other hand, the PMSM system of the MSFESS is subjected to multi-source disturbances such as mass imbalance torque, magnetic bias force, gyroscope effect, etc. The disturbance torques have been introduced in [24] and are omitted here. The feedback current signals contain noise, which can also be disposed as disturbances. To achieve satisfactory control performance, the SLE system should be sensitive to position signal and simultaneously be robust to the unknown disturbances.
Set the state variable as x = [ i α , i β , e α , e β ] T , control input as u = [ u α , u β ] T , and two-phase static current after CLARKE transformation as measured output y = [ i α , i β ] T . For the PMSM, the mentioned disturbance torques act directly on the control forward loop, which can be equivalent to a disturbance voltage in the control input. Therefore, the unknown disturbance is composed of the disturbance voltage v d and measurement noise n d , i.e., d = [ v d , n d ] T . In the case of unknown disturbance, the following state-space model is obtained based on the kinematic model in the α/β coordinate:
{ x ˙ = A x + B u + B d d , y = C x + D d d ,  
where
A = [ R s L s 0 1 L s 0 0 R s L s 0 1 L s 0 0 0 ω r 0 0 ω r 0 ] , B = [ 1 L s 1 L s 0 0 ] , B d = [ 1 L s 0 1 L s 0 0 0 0 0 ] , C = [ 1 0 0 0 0 1 0 0 ] , D d = [ 0 0.1 0 0.1 ] ,  
with the variance of the noise set to 0.1 according to the circuit condition. It is easy to know the plant is completely observable and hence the observer is excitable.
For the purpose of SLE, the following extended slide-mode observer is considered [25,26,27]:
{ x ^ ˙ = A ^ x ^ + B ( u u s m o ) y ^ = C x ^ u s m o = K s l i d e s a t ( y y ^ )  
where x ^ , y ^ are the estimates of x and y , respectively. The observer matrix is A ^ = [ R s L s 0 1 L s 0 0 R s L s 0 1 L s 0 0 ω ^ r 0 0 ω ^ r 0 ] , with the estimated speed. The goal of the saturation control z is to drive current estimation error to zero. This is achieved by proper selection of the slide-mode gain K s l i d e and correct formation of the estimated back emf.
Provided that the evaluation function is z = [ e e α , e e β ] T with the estimated error e e α = e α e ^ α , e e β = e β e ^ β . Then we have
z ( s ) = G z u ( s ) u ( s ) + G z d ( s ) d ( s ) .  
It is expected that the transfer matrix gain from the disturbance d ( s ) to evaluation function z ( s ) is minimal. It is the so-called H optimization control or sub-optimization control, i.e.,
G z d ( s ) min , or G z d ( s ) < γ .  
Due to this, the impact of the unknown disturbance can be rejected within an adaptable range, even with the worst disturbance circumstance. In this paper, we propose an optimal SLE method based on the extended sliding-mode observer in order to accomplish SLE in the whole speed range.
Remark 1.
It is seen that the α/β current mode has two basic states iα and iβ as described in (2). To improve the state estimation accuracy, the back electromotive forces eα and eβ can be extended new states as shown in (5). Therefore, it is usual to call this sliding mode observer of (6) as an extended observer and some references like [27] have proposed this observer.

4. Design of Sensorless Estimation System

4.1. H Optimization-Based Sensorless Estimation

In this section, we start with the design of SLE system by applying the traditional extended sliding-mode observer method for the SPMSM.
By defining e = x x ^ = [ e x α e x β e e α e e β ] T , it follows from (5)–(7) that:
{ e ˙ = A ^ e B u s m o + B d d , z = C z e , u s m o = K s l i d e s a t ( y y ^ ) ,  
where the evaluation matrix C z = [ 0 0 1 0 0 0 0 1 ] . The estimation error e x α = x α x ^ α , e x β = x β x ^ β , e e α = e α e ^ α , e e β = e β e ^ β .
In Equation (9), the saturation function is denoted as follows:
s a t ( y y ^ ) = { y ˜ t h , ( y y ^ ) < y ˜ t h , ( y y ^ ) , y ˜ t h < ( y y ^ ) < y ˜ t h , y ˜ t h , ( y y ^ ) > y ˜ t h ,  
where y ˜ t h = y y ^ is a preset constant.
According to (9) and (10):
u s m o = { K s l i d e y ˜ t h , ( y y ^ ) < y ˜ t h , K s l i d e ( C e + D d d ) , y ˜ t h < ( y y ^ ) < y ˜ t h , K s l i d e y ˜ t h , ( y y ^ ) > y ˜ t h ,  
Then, the error equation can be rewritten as:
( y y ^ ) < y ˜ t h , { e ˙ = A ^ e + B K s l i d e y ˜ t h + B d d , z = C z e , | y y ^ | < y ˜ t h , { e ˙ = ( A ^ B K s l i d e C ) e + ( B d B K s l i d e D d ) d , z = C z e , ( y y ^ ) > y ˜ t h , { e ˙ = A ^ e B K s l i d e y ˜ t h + B d d , z = C z e ,  
From Equation (12), it can be seen that the observer stability depends on the system matrix A ^ and ( A ^ B K s l i d e C ) . This is the typical field of the modern control theory. In this paper, we focus on the transfer matrix from unknown disturbance d to evaluation function z is minimal, i.e., H optimization control. The case in | y y ^ | < y ˜ t h is considered with another special cases.
According to the traditional extended slide-mode observer in (6), the state estimates x ^ = [ i ^ α i ^ β e ^ α e ^ β ] can be obtained. Owing to the measurement of current i α and i β , the values of the errors e x α and e x β can also be determined with e x α = i α i ^ α , e x β = i β i ^ β .
Then, considering the errors e x α and e x β as the measured output y e = [ e x α e x β ] , the error equation in (12) can be rewritten as
| y y ^ | < y ˜ t h , { e ˙ = ( A ^ B K s l i d e C ) e + ( B d B K s l i d e D d ) d , y e = C e e , z = C z e ,  
where C e = [ 1 0 0 0 0 1 0 0 ] .
Generally, the sliding mode observer gain K s l i d e is designed as follows [17,18,19]:
K s l i d e > max ( | e α | , | e β | ) ,  
where the sliding mode plane exists only if the gain K s l i d e is larger than the back electromotive force in the α/β coordinate.
It is obvious that the formulation in (14) just provides the upper bound of the observer gain. The detailed value of the gain is not determined so that the states’ estimation accuracy cannot be guaranteed. In this way, the value of the gain can be approximately confirmed based on the formulation in (14).
Therefore, the error estimation equation can be further formulated as:
  { e ^ ˙ = A e e ^ + L e ( y e y ^ e ) , y ^ e = C e e ^ , z ^ = C z e ^ ,  
where the matrices A e = A ^ B K s l i d e C , B e = B d B K s l i d e D d are constant matrices.
Combining Equations (13) and (15), it is obtained that:
  { e ˜ ˙ = ( A e L e C e ) e ˜ + B d d , z ˜ = C z e ˜ ,  
where the evaluation output error z ˜ = z z ^ .
Based on this, the evaluation output is further expressed as:
z ˜ ( s ) = G z d ( s ) d ( s ) ,  
where the transfer function matrix is:
G z d ( s ) = C z ( s I A e + L e C e ) 1 B d .  
Lemma 1.
Considering the following linear system [23]
x ˙ ( t ) = A x ( t ) + B w ( t ) , z ( t ) = C x ( t ) + D w ( t ) ,  
Given a positive coefficient γ > 0 , if there exists a positive matrix P > 0 with the following linear matrix inequality satisfied:
[ A T P + P A P B C T B T P γ I D T C D γ I ] < 0  
Then, the original linear system is robustly stable and the following H-norm is met:
G z w ( s ) < γ ,   w .  
By using Lemma1, the error equation in (16) can be transferred as:
γ min ,   s . t . [ ( A e L e C e ) T P + P ( A e L e C e ) P B d C z T B d T P γ I 0 C z 0 γ I ] < 0 , P > 0  
The formulation in (18) can be considered as a traditional H optimization problem which is able to be disposed by the linear matrix inequality and objective function. Here, the detailed process is omitted.
Remark 2.
It is certain that the H-optimization control method has been widely used in reality since the worst disturbance case has been considered. The influence of the disturbance on the estimation performance can be reduced effectively by using the H-optimization technique. But this approach has a little conservativenss, and is especially subject to a large estimation error. Due to this, this paper adapts a traditional observer as the first step to obtain the basic states estimations, and the second one is further designed to compensate the errors based on H optimization control.

4.2. Back Electromotive Force Error Compensation

After designing the observer gain L e by using the linear matrix inequality, the error e = [ e x α e x β e e α e e β ] T can be obtained. Then, putting the estimation error into the observer (13), we get
e ^ α 1 = e ^ α + e e α ,   e ^ β 1 = e ^ β + e e β ,  
where e ^ α 1 , e ^ β 1 are estimates of back electromotive forces after compensation.
According to the relationship in (4), the estimates of the angular position and speed can be configured as:
{ θ ^ r = a r c tan e ^ α 1 e ^ β 1 , ω ^ r = d θ ^ r d t .  
Due to the measurement noise in the current i α and i β , it is inevitable that the estimates of the back electromotive forces e ^ α and e ^ β contain much high-frequency disturbances. For this problem, the traditional method is to design a low-pass filter to reject the noise.
θ ^ r = a r c tan e ^ α 1 e ^ β 1 f f i l t e r ( s ) .  
Usually, the low-pass filter is designed as a one-order filter as follows:
f f i l t e r ( s ) = ω c s + ω c ,  
where ω c is the cut-off frequency.
It is known that there is some phase delay in the low-pass filter, resulting in a postponed angle estimation error. In some references, the phase compensation method has been used to make up the angle error as follows:
Δ θ r = arctan ( ω ^ r ω c ) .  
Based on this, the angular position after compensation is determined as:
θ ^ r = arctan ( e ^ α 1 e ^ β 1 ) + arctan ( ω ^ r ω c ) .  
In this way, the two-step angle error estimation and compensation strategy can be accomplished. The control scheme block is shown in Figure 4. The overall angular position estimation and compensation algorithm can be summarized in Algorithm 1.
Algorithm 1. Given the original system (5) and a sampling period Ts, set computational number k = 0 and initial time t 1 = t 2 = 0 . The observer gain matrix K s l i d e and L have been calculated based on formulation (14) and Lemma 1.
Step 1. Sample the Hall sensor to refresh the measurements i a i b i c V a V b V c .
Step 2. Substitute i a i b i c V a V b V c into (2) to calculate the parameters i α i β V α V β in the α/β coordinate.
Step 3. Put i α i β V α V β into (6) to calculate the estimation parameters i ^ α i ^ β e ^ α e ^ β .
Step 4. Substitute i ^ α i ^ β e ^ α e ^ β into (15) to calculate the error estimation parameters e ^ i α e ^ i β e ^ e α e ^ e β .
Step 5. Combine e ^ e α e ^ e β into (19) to obtain the error compensation parameters e ^ α 1 e ^ β 1 .
Step 6. Design the low-pass filter in (22) to compensate the error compensation parameters e ^ α 1 e ^ β 1 to accomplish the position estimation and compensation according to (24).
Remark 3.
For the traditional application, the first-step observer is suitable for a medium speed or frequency circumstance. In this case, it is not essential for the second-step observer combination. The existence of the proposed observer is always true since the original system is robustly stable. For the MSFESS, the design of the second-step observer can further improve the angle estimation accuracy.

5. Experimental Results

5.1. Experimental Setup

The experimental setup is shown in Figure 5. The MSFESS consists of an aluminium outermost shell, inside PMSM and 5-freedom AMB, charging-discharging power unit and controller, master control screen and the outside black cabinet. The input power is 540 VDC voltage with the bus link, while the assistant 3-phase 380 VAC power is provided for AMB, vacuum and fan systems. The protection earth is configured near to the front door for prevention of the electric leakage. All of the required parameters are shown in Table 1.
Accordingly, the single gimbal system (4) can be governed with the following matrices:
A = [ 333.3 0 3333.3 0 0 333.3 0 3333.3 0 0 0 2198 0 0 2198 0 ] , B = [ 3333.3 3333.3 0 0 ] , B d = [ 3333.3 0 3333.3 0 0 0 0 0 ] , C = [ 1 0 0 0 0 1 0 0 ] , D d = [ 0.1 0.1 ] .  
Based on the rated voltage of the MSFESS in Table 1, the back electromotive forces can be near to 600 / 3 = 346.4 ( V A C ) , i.e., e a ,   e b ,   e c 346.4 ( V A C ) . After CLARKE transformation 3s/2s, the back electromotive forces in coordinate can be estimated as max ( e α ,   e β ) max ( e a ,   e b ,   e c ) . Therefore, the first-step observer gain K s l i d e can be designed as K s l i d e = 346.4 .
After that, by using the LMI box tool in MATLAB software, the solution of the inequality in Lemma 1 can be obtained exactly. Based on that, the positive coefficient and matrix can be calculated as:
γ = 0.01 ,   P = [ 333.3 0 3333.3 0 0 333.3 0 3333.3 0 0 0 2198 0 0 2198 0 ] .

5.2. Angular Position Estimation Results

For the MSFESS, the flywheel can be accelerated to a high speed when the magnetically suspended system can function properly since the flywheel axis will strike the bearing without suspension. However, when the motor is controlled in a low speed, it can also be safe even if the AMB system is powered off and in that case, the flywheel rotor will stay connected with the bearing body. For this working condition, a kind of position sensor can be installed at the terminal of the rotor axis and applied as the angle benchmark to evaluate the observer estimation result. Figure 6 shows the real picture of the resolver, consisting of a stator and rotor, while the rotor is fixed at the terminal of the axis and the stator is connected to the flywheel body. However, if the flywheel axis is suspended by the AMB system, the clearance between the stator and rotor will exceed the maximum tolerance value of the resolver, and in this case, the resolver cannot be used any more.
For the sensorless control, the SPMSM will be started in an open-loop mode since there is no initial current in the motor stator and the observer has no angle estimation. In the open-loop start-up mode, a simulated angle function is generated by the program code and the motor will rotate gradually according to the frequency value of the angle function. The experiment is developed to compare the estimated performance by using one-step and the proposed two-step observer. Figure 7 and Figure 8 show the resolver angle and its estimation by the proposed observer. In this condition, the flywheel speed is set to 1200 rpm and 3000 rpm with the angle frequency 40 Hz and 100 Hz, respectively. Actually, the motor can be started in the position sensor mode but it is not practical for the sensorless control strategy. In these curves, the resolver angle is located in the field of 0~1, i.e., 0~360 deg (red line) and a constant of 1.5 has been added to the estimated angles for simplicity (blue line), while the green line is the estimated angle error.
It is seen that at both 1200 rpm and 3000 rpm, the estimated angular position can track the simulated and resolver angles. At a low speed, the simulated angle can be accurate with respect to the real resolver angle but it is not correct in a high-speed range. Additionally, it is not obvious in the estimation performance by using the two methods since the common sliding-mode observer is accurate enough for the low speed circumstance.
In order to evaluate the estimation performance of the proposed observer, the resolver sensor is removed and the AMB system is powered for flywheel full-suspension. Due to this, the flywheel speed is set to 9000 rpm with the angle frequency 300 Hz. Also, the open-loop start-up method is used to achieve the typical V/F control algorithm. Figure 9 shows the estimation results by using the simulated and estimation strategy. It is demonstrated that the estimated value by using the one-step observer is not perfect with a lagging phase and fluctuant magnitude; there is much noise in the estimation results. By using the proposed two-step observer, the estimated angular position can also be appropriate for the simulated value, and the estimated results are smoother, which are effective for the control application.

5.3. Control Results

In order to test the control performance by using the sensorless control strategy, the estimated angular position is fed back to the control loop in Figure 3 to accomplish the flywheel acceleration process. The flywheel speed is set to the medium-speed 9000 rpm and the two steps observer are experimented. In the closed-loop control structure, the Iq and Id currents are the main control objects. When the currents in the d/q coordinate are controlled perfectly, the control performance of the outer speed loop or voltage loop will be better. Therefore, the control results of the Iq and Id current loops are able to indicate the control performance improvement by using the two steps for the sliding-mode observer. In the charging mode with the speed from 0 to 9000 rpm, Figure 10 shows the current measurements of Iq and Id, while the yellow line is the reference Iqref with the blue one Iq the feedback value, and the green one is the Id feedback value with its reference Idref zero.
In the charging mode of the MSFESS, the constant torque strategy is used to accelerate the rotor speed. In this case, the effective charging current remains constant about Iqref = 100 A (as shown in Figure 10) so that the ideal Iq feedback value could be used to track the reference value Iqref. The Iqref = 0 strategy is adapted to minimize the invalid current. It can be seen from Figure 10 that both the Iq and Id current measurements have larger fluctuations, while the maximum of the Iq feedback value can reach 150 A with 50% overshot and the maximum overshot of the Id feedback value can also be worse than 55%. These control conditions cannot be applied in the flywheel charging-discharging process. However, by using the proposed two-step observer, both the fluctuations of the Iq and Id current measurements can be largely reduced, and the maximum of the Iq feedback value can be less than 110 A with 10% overshot, while maximum overshot of the Id feedback value can also be better than 15%. These control conditions can satisfy the flywheel rapid charging-discharging requirements. Above all, the proposed sensorless control algorithm can improve the control performance of the MSFESS in the high-speed field through the two-step sliding-mode observer and H control method.

6. Conclusions

In this paper, the problem of angular position estimation in the sensorless control strategy has been studied for the MSFESS. Through analysis of the system modeling and principle description, it is obvious that the sensorless control has to be essential for the MSFESS based on the rotor angular position estimation strategy. Compared with other estimation methods, the sliding-mode observer algorithm has the advantages of robustness, reliability and simplicity. Based on this, a kind of two-step sliding-mode observer has been designed, including a traditional first-step sliding-mode and an extended second-step observer. In this scheme, the first step is used to guarantee the stability and convergence of the observer, while the second one can improve the estimation accuracy of the states. For the second-step observer design, the H optimization approach has been used based on the linear matrix inequality solving such that the influence of the unknown disturbances on the state estimation accuracy can be rejected within a low-level range. Finally, real experiments have been developed to demonstrate the effectiveness of the proposed method. It is shown that the angular position estimation accuracy can be further improved by using the two-step sliding-mode observer compared with the one step. Thus, the control results in both charging and discharging processes can also satisfy the desired performance.

Author Contributions

Conceptualization, S.L. and Y.F.; Methodology, S.L.; Software, S.L.; Validation, S.L., Y.F. and P.L.; Writing-Original Draft Preparation, S.L. and P.L.; Writing-Review & Editing, S.L. and P.L.; Visualization, X.X.; Supervision, X.X.; Project Administration, Y.F.; Funding Acquisition, Y.F.

Funding

This work was supported in part by the 61th Postdoctoral Science Foundation of China (No. 2017M610735), the National Nature Science Foundation of China (No. 6151001034) and the Key Project of Ministry of Science and Technology of China (No. 2018YFB0905500).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dai, X.J.; Deng, Z.F.; Liu, G. Review on advanced flywheel energy storage system with large scale. Trans. China Electrotech. Soc. 2011, 26, 133–140. [Google Scholar]
  2. Wang, C.D. Development of flywheel battery. Power Technol. 2014, 37, 732–745. [Google Scholar]
  3. Chen, Y.; Gann, S.; Zhou, J. Flywheel energy storage technology. Power Technol. 2016, 40, 1718–1721. [Google Scholar]
  4. Li, S.S.; Fu, Y.Y.; Liu, P.; Wang, Z.Q. Research on Integrated Application and Charging-discharging Control Method for the Magnetically Suspended Flywheel Storage-based UPS System. Proc. CSEE 2017, 37, 170–176. [Google Scholar]
  5. Li, S.S.; Liu, P.; Fu, Y.Y. Research on problem of the transient charging-discharging switching for the magnetically suspended flywheel-based PMSM system. Small Spec. Electr. Mach. 2018, accepted. [Google Scholar]
  6. Silva-Saravia, H.; Pulgar-Painemal, H.; Mauricio, J.M. Flywheel energy storage model, control and location for improving stability: The Chilean case. IEEE Trans. Power Syst. 2017, 32, 3111–3119. [Google Scholar] [CrossRef]
  7. Zhang, X.; Yang, J.Q. A robust flywheel energy storage system discharge strategy for wide speed range operation. IEEE Trans. Ind. Electron. 2017, 64, 7862–7873. [Google Scholar] [CrossRef]
  8. Müller, T.; See, C.; Ghani, A.; Bati, A.; Thiemann, P. Direct flux control—Sensorless control method of PMSM for all speeds—Basics and constraints. Electron. Lett. 2017, 53, 1110–1111. [Google Scholar] [CrossRef]
  9. Paitandi, S.; Sengupta, M. Analysis, design and implementation of sensorless V/f control in a surface-mounted PMSM without damper winding. Sadhana Acad. Proc. Eng. Sci. 2017, 42, 1317–1333. [Google Scholar] [CrossRef]
  10. Chaoui, H.; Okoye, O.; Khayamy, M. Current sensorless MTPA for IPMSM drives. IEEE-ASME Trans. Mechatron. 2017, 22, 1585–1593. [Google Scholar] [CrossRef]
  11. Almarhoon, A.H.; Zhu, Z.Q.; Xu, P.L. Improved pulsating signal injection using zero-sequence carrier voltage for sensorless control of dual three-phase PMSM. IEEE Trans. Energy Convers. 2017, 32, 436–446. [Google Scholar] [CrossRef]
  12. Ni, R.; Xu, D.; Blaabjerg, F.; Lu, K.; Wang, G.; Zhang, G. Square-wave voltage injection algorithm for PMSM position sensorless control with high robustness to voltage errors. IEEE Trans. Power Electron. 2017, 32, 5425–5437. [Google Scholar] [CrossRef]
  13. Tian, L.S.; Yu, Q. Wide-speed-range sensorless control scheme for IPMSMs based on super-twisting observer and HF signal injection. Int. J. Appl. Electromagn. Mech. 2017, 54, 367–387. [Google Scholar] [CrossRef]
  14. Urbanski, K. A new sensorless speed control structure for PMSM using reference model. Bull. Pol. Acad. Sci. Tech. Sci. 2017, 65, 489–496. [Google Scholar] [CrossRef] [Green Version]
  15. Jung, T.U.; Jang, J.H.; Park, C.S. A back-EMF estimation error compensation method for accurate rotor position estimation of surface mounted permanent magnet synchronous motors. Energies 2017, 10, 1160. [Google Scholar] [CrossRef]
  16. Feng, G.D.; Lai, C.Y.; Kar, N.C. Expectation-maximization particle-filer-and Kalman-filter-based permanent magnet temperature estimation for PMSM condition monitoring using high-frequency signal injection. IEEE Trans. Ind. Inform. 2017, 13, 261–1270. [Google Scholar]
  17. Chang, X.J.; Liu, L.; Ding, W. Novel nonsingular fast terminal sliding mode control for a PMSM chaotic system with extended state observer and tracking differentiator. J. Vib. Control 2017, 23, 2478–2493. [Google Scholar] [CrossRef]
  18. Liang, D.L.; Li, J.; Qu, R.H. Sensorless control of permanent magnet synchronous machine based on second-order sliding-mode observer with online resistance estimation. IEEE Trans. Ind. Appl. 2017, 53, 3672–3682. [Google Scholar] [CrossRef]
  19. Sheng, L.C.; Li, W. Sensorless control of a shearer shor-range cutting interior permanent magnet synchronous motor based on a new sliding mode observer. IEEE Access 2017, 5, 18439–18450. [Google Scholar] [CrossRef]
  20. Bariša, T.; Sumina, D.; Pravica, L.; Čolović, I. Flying start and sensorless control of permanent magnet wind power generator using induced voltage measurement and phase-locked loop. Electr. Power Syst. Res. 2017, 152, 457–465. [Google Scholar] [CrossRef]
  21. Zhang, X.M.; Han, Q.L. Network-based H∞ filtering using a logic jumping-like trigger. Automatica 2013, 49, 1428–1435. [Google Scholar] [CrossRef]
  22. Zhang, X.M.; Han, Q.L.; Yu, X.H. Survey on Recent Advances in Networked Control Systems. IEEE Trans. Ind. Inform. 2016, 12, 1740–1752. [Google Scholar] [CrossRef]
  23. Zhong, M.Y.; Li, S.S. Robust optimal Hi/H∞ fault detection and its applications in a TGSCM system. Trans. Inst. Meas. Control 2018, 40, 719–732. [Google Scholar] [CrossRef]
  24. Li, S.S.; Fu, Y.Y.; Liu, P. Research on disturbance analysis and detection method for the magnetically suspended flywheel-based PMSM system. Energy Storage Sci. Technol. 2018, 7. [Google Scholar] [CrossRef]
  25. Zhang, B.L.; Han, Q.L.; Zhang, X.M. Recent advances in vibration control of offshore platforms. Nonlinear Dyn. 2017, 89, 755–771. [Google Scholar] [CrossRef]
  26. Zhang, B.L.; Han, Q.L.; Zhang, X.M.; Yu, X. Sliding mode control with mixed current and delayed states for offshore steel jacket platforms. IEEE Trans. Control. Syst. Technol. 2014, 22, 1769–1783. [Google Scholar] [CrossRef]
  27. Mohamed, G.; Sofiane, A.A.; Nicolas, L. Adaptive super twisting extended state observer based sliding mode control for diesel engine air path subject to matched and unmatched disturbance. Math. Comput. Simul. 2018, 151, 111–130. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of a magnetically suspended flywheel energy storage system (MSFESS).
Figure 1. Schematic diagram of a magnetically suspended flywheel energy storage system (MSFESS).
Sensors 18 02467 g001
Figure 2. Minimum power scheme of permanent magnet synchronous motor (PMSM) driver. IGBT: insulated gate bipolar transistor.
Figure 2. Minimum power scheme of permanent magnet synchronous motor (PMSM) driver. IGBT: insulated gate bipolar transistor.
Sensors 18 02467 g002
Figure 3. PID-based vector control block diagram. SPMSM: surface-mounted PMSM; SVPWM: space vector pulse width modulation.
Figure 3. PID-based vector control block diagram. SPMSM: surface-mounted PMSM; SVPWM: space vector pulse width modulation.
Sensors 18 02467 g003
Figure 4. Two-step angular position estimation and compensation.
Figure 4. Two-step angular position estimation and compensation.
Sensors 18 02467 g004
Figure 5. Experimental setup. UPS: unterruptible power supply.
Figure 5. Experimental setup. UPS: unterruptible power supply.
Sensors 18 02467 g005
Figure 6. Resolver in the MSFESS.
Figure 6. Resolver in the MSFESS.
Sensors 18 02467 g006
Figure 7. The simulated angle and its estimation at 1200 rpm by using one-step observer (a) and two-step observer (b).
Figure 7. The simulated angle and its estimation at 1200 rpm by using one-step observer (a) and two-step observer (b).
Sensors 18 02467 g007
Figure 8. The simulated angle and its estimation at 3000 rpm by using one-step observer (a) and two-step observer (b).
Figure 8. The simulated angle and its estimation at 3000 rpm by using one-step observer (a) and two-step observer (b).
Sensors 18 02467 g008
Figure 9. The simulated angle and its estimation at 9000 rpm by using one-step observer (a) and two-step observer (b).
Figure 9. The simulated angle and its estimation at 9000 rpm by using one-step observer (a) and two-step observer (b).
Sensors 18 02467 g009aSensors 18 02467 g009b
Figure 10. The close-loop control performance at 9000 rpm by using one-step observer (a) and two-step observer (b).
Figure 10. The close-loop control performance at 9000 rpm by using one-step observer (a) and two-step observer (b).
Sensors 18 02467 g010
Table 1. MSFESS parameters.
Table 1. MSFESS parameters.
ParameterIndexSymbol
Electric capacity3 kWhEN
Rated speed10,500 rpmnN
Rated power250 kVAPN
Rated voltage600 VVN
Rated current255 AIN
Stator resistance0.1 ΩRs
Stator inductance300 μHLs
Pole-pairs2p
Rotor moment of inertia16.2 kg·m2JN
Bearing typeAMB-
Vacuum value<10 Pa-

Share and Cite

MDPI and ACS Style

Li, S.; Fu, Y.; Liu, P. Position Estimation and Compensation Based on a Two-Step Extended Sliding-Mode Observer for a MSFESS. Sensors 2018, 18, 2467. https://doi.org/10.3390/s18082467

AMA Style

Li S, Fu Y, Liu P. Position Estimation and Compensation Based on a Two-Step Extended Sliding-Mode Observer for a MSFESS. Sensors. 2018; 18(8):2467. https://doi.org/10.3390/s18082467

Chicago/Turabian Style

Li, Shusheng, Yongling Fu, and Ping Liu. 2018. "Position Estimation and Compensation Based on a Two-Step Extended Sliding-Mode Observer for a MSFESS" Sensors 18, no. 8: 2467. https://doi.org/10.3390/s18082467

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