Mathematical Problems in Engineering
Volume 2020 (2020), Article ID 1849816, 9 pages
https://doi.org/10.1155/2020/1849816

Research Article

A Discrete Nonlinear Tracking-Differentiator and Its Application in Vibration Suppression of Maglev System

Zhiqiang Wang,1,2 Zhiqiang Long,¹ Yunde Xie,³ Jingfang Ding,4,5 Jie Luo,¹ and Xiaolong Li¹

¹Engineering Research Center of Maglev, National University of Defense Technology, Changsha, Hunan 410073, China
²Hunan Provincial Key Laboratory of Electromagnetic Levitation and Propulsion Technology, Changsha, Hunan 410073, China
³Beijing Enterprises Holding Maglev Technology Development Company Limited, Beijing 100000, China
⁴Laboratory of Networked Control Systems, Shenyang Institute of Automation, Chinese Academy of Sciences, Shenyang 110016, China
⁵University of Chinese Academy of Sciences, Beijing 100049, China

Correspondence should be addressed to Yunde Xie; xieyd9999@163.com

Received 8 December 2019; Revised 1 March 2020; Accepted 5 March 2020; Published 19 June 2020

Academic Editor: Filippo Cacace

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

Abstract.

Vibration of the maglev train levitation system is harmful to riding comfort and safety. The signal processing method is effective in vibration control. In this paper, a novel kind of second-order nonlinear tracking differentiator is proposed and applied to suppress the vibration phenomenon. The switching curves of the second-order discrete time optimal control system are presented by the isochronous region method. A synthetic function is acquired depending on the cases whether a point in the phase plane can reach the switching curves within one sample step. The discrete form of tracking differentiator is constructed based on the position relationship between the state point and the characteristic curves. Numerical simulation shows that this discrete tracking differentiator can quickly track an input signal without overshoot and chattering and can produce a good differential signal. A test in a maglev test bench also demonstrates the effectiveness of the tracking differentiator in the suppression of the track-train vibration.

1. Introduction

Maglev train is a new type of modern railway transportation tool. Compared with the conventional wheel-rail train, it possesses the advantages of steady and comfortable riding, low noise, small turning radius, and strong climbing ability. The maglev train has a wide application prospect; the Changsha Maglev Express and the Beijing Maglev S1 are well welcomed by passengers. Making the train levitate above the track steadily is the foundation for the operation of the maglev train. To accomplish this purpose, it is necessary to gain the levitation gap of the train, vertical acceleration of the levitation electromagnet, the current of the electromagnet, and other information through the sensors as the feedback and then calculate the control output based on all the feedback. The quality of the sensor signals and corresponding signal processing technique play an important role in the performance of the levitation system.

To simplify the design of control system and to improve system performance, an effective way to acquire differential signal is of great importance. Tracking differentiator (TD) is firstly proposed for this purpose [1]. TD can also be used to design nonlinear PID controller, estimate parameters, and acquire differential information. Also, TD can be employed as a basic component of ADRC (auto disturbance rejection control) in charge of arranging transition process. Han proposed a discrete tracking differentiator based on second-order time optimal system. The differentiator can acquire the differential signal effectively, with the property of fast tracking input signal and no overshoot and chattering [2]. Lots of research studies are carried out about tracking differentiator [3–11]. Xie and Long [12] proposed a high accuracy fast TD, and this TD is used in signal detection in the maglev control system [13, 14]. Other applications of TD can be seen in [15–18].

Digital control systems are common in modern control engineering, which makes the study of discrete form controllers and signal processors valuable. For discrete form TD, the region of the linear zone needed in switching function has a nonignorable influence on TD’s tracking and differential performance. Literature [2] uses the isochronous region method to specify the linear zone when time optimal control adopts nonextreme value; it presents the normal form of control synthetic function; however, the form is too complicated to be implemented in engineering practice. Another strategy, searching the boundary characteristic points which makes the time optimal control adopt nonextreme value and construct sectionwise linear function to avoid complicated computation, is used in this paper. Numerical simulation shows that this new form TD has the same advantage of the form in [2], and it can fast track input signal without overshoot and chattering and can acquire effective differential signal. However, the new form, which is a piecewise linear function, is much more simplified. The characteristic points are obtained in this strategy when the fastest control function adopts nonoptimal value based on the isochronous region method, and then a piecewise fastest control function is constructed based on the characteristics of the fastest control system.

The remainder of this paper is organized as follows. First, a description about maglev train levitation system is given in Section 2, and then in Section 3, a detailed procedure of how this tracking differentiator is derived is given. In Section 4, simulation results are given to verify the effectiveness of this tracking differentiator in suppression of maglev train vibration, and Section 5 concludes this paper.

2. Maglev Train Levitation System

Figure 1 is a photo of Changsha Maglev Express, which connects the airport and the railway station. The length of Changsha Maglev Express operation line is 18 km, which is the longest ever for mid-low speed maglev train.

Figure 1: Photo of Changsha Maglev Express.

The interactions between the electromagnets and the tracks play a key role in maglev train. The structure and schematic diagram of maglev train levitation system can be seen in Figures 2(a) and 2(b), respectively. A fundamental maglev system consists of electromagnet, sensor, and controller. In the research of maglev trains, since the mechanical decoupling of the suspension modules on both sides has been achieved by the bogie, a single point is usually used as the research object of the suspension control system. We simplify the single-point suspension system into the model shown in Figure 2(b) and make the following assumptions:

(i)Magnetic leakage and edge effects are ignored, and the magnetic flux is considered to be evenly distributed on the air gap.(ii)The magnetic resistances of guide rails, iron cores, and plates are ignored.(iii)It is considered that the levitation force provided by the electromagnet is concentrated on the geometric center and the geometric center of the electromagnet coincides with the center of mass.(iv)It is believed that there is no dislocation between the magnetic pole surface and the guide rail, that is, the electromagnet has no rolling movement relative to the track, but only vertical movement. The force of the air spring on the electromagnet is also vertical.

(a)

(b)

Figure 2: Maglev train levitation system. (a) Structure of maglev train levitation system. (b) Schematic diagram of maglev train levitation system.

2.1. Reference Coordinates

When the interference of track irregularity is considered, a new absolute reference surface is needed. In Figure 2(b), represents the displacement of the track relative to the reference plane, is the absolute displacement of the electromagnet relative to the absolute reference plane, and is the gap between the track and the electromagnet.

2.2. Kinetic Equation

From Newton’s second law, the kinetic equation can be obtained as:

where is the equivalent mass of electromagnet, is the equivalent mass of carriage, is the electromagnetic force, is the current in the electromagnet coil, is the increment of the force exerted on the bogie by the air spring relative to the static state, and its acting point is the installation position of the air spring. Because the frequency of the air spring itself is far lower than the natural frequency of the suspension system, in normal circumstances, the displacement of the suspension module relative to its balance position is also very small due to the control function. Therefore, compared with the electromagnetic force provided by the electromagnet, can be ignored.

2.3. Electrical Equation

where is the control voltage applied at both ends of the coil windings, is the coil impedance, is the number of turns of the coil windings, is the pole area of the electromagnet, and is the vacuum permeability.

2.4. Electromagnetic Force Equation

2.5. Boundary Condition

where and are the values of the coil current and gap in steady state, respectively, and and are also the initial values of the track and electromagnet relative to the reference plane in steady state.

Making a Taylor expansion of the above model at the boundary equilibrium point and denoting , a linearized model can be acquired.

Since electromagnets fluctuate in a small range near the equilibrium point, the term in the electrical equation can be ignored, and the inductance can be approximated as . Taking the state variable as , the state space representation of the maglev system can be obtained as follows:

3. Second-Order Discrete System Control Synthesis Function

3.1. Preliminaries

For a discrete form double integral time optimal control (TOC) system:

where is the state variable, is the control variable with constrained amplitude , and is the sample step; the corresponding state space form expression can be obtained as:

in which . A universal solution of equation (9) with initial value of is:

If after steps, the state point finally arrives at the origin in the phase plane, that is, , then from (11), the mathematical condition the initial point follows is obtained as:

Substituting and into (11) yields:

From optimal control theory for continuous system, any point in the phase plane needs at most one switch before reaching at the origin along the optimal trajectory. For continuous system, the switching is finished at instance with no delay. However, for discrete control system, the switching process takes place within a sample step after the state point passes the switching line with the continuous control law, which is also called a “bang-bang” control strategy. Define those points which can converge to the origin when or always as or . Define those points which can converge to the origin when first then always or first then always as or , . Define those points which can converge to the origin when first then or always as or , . Define the curve constituted by as , by as , by as , by as , by as , and by as . The next task is to determine the mathematical descriptions of these boundary lines.

3.2. Determination of the Boundaries

In this part, the mathematical descriptions for the boundaries mentioned above are determined in detail. Firstly, the boundary is determined. Suppose that a state point reaches the origin in steps, that is, , and . From (12),

Then,

For simplicity, denote as and as ; then, the following equation can be derived:

Formula (15) is the mathematical description of . In the same way, the mathematical description of can be derived as follows:

Formulas (15) and (16) can be combined into one unity formula as follows:

Then, the boundary line is to be determined. For state points , there are , and , , . Then, the following relationships can be obtained:

For simplicity, denote as and as . From formulas (18) and (19), relation (20) can be derived.

Formula (20) is the mathematical description of ; in the same way, the mathematical description of is as follows:

Combining (20) and (21), the description of is as follows:

Furthermore, the boundary is to be determined. For state points , there are , , , . For state points , there are , , , . For simplicity, denote as and as , and the mathematical description of can be derived in the same way:

Then, all boundaries of region and control characteristic lines are determined by (17), (22), and (23). The boundary curves can now been obtained as shown in Figure 3.

Figure 3: Boundaries of linear regions.

3.3. Determination of Control Synthesis Function

Linear zone (the area between and ) and its boundary lines have already been determined in the above part. For any state point in the phase plane, if , then control value adopts an extremum; if , changes according to some linear law in the range of .

As shown in Figure 4, any point in the phase place converges to the origin going through the switching line at most once, and that sign-changing process takes place in linear zone at the adjacent of switching line . One idea is to confirm the steps it needs to reach the switching line . If is no less than 1, then it means . Under this condition, adopts an extremum . Otherwise, adopts some particular value and . Suppose that is above the switching line and after steps it reaches to the point . From (10), , . That is,

Figure 4: Boundaries and for linear region .

The value of can be obtained by solving these equations as . Since in (24), and , . Since and , can be further rewritten as follows:

In the same way, the steps it takes for point below switching line to reach the switching line can be derived as follows:

These two conditions can be combined into a unified one:

If , it is obvious that . If , can reach within one step. Under this condition, a proper control value is needed to make arrive at exactly. Suppose at control value , will arrive at , then it will reach the origin along under or along under . From (12),

Then, the following relationship can be derived:

The solution of in terms of and can be also derived as follows:

Since , has the same sign with . Then, the following solution is also true: . This form can further be rewritten as:

in which:

One root of (32) is to be excluded. From (28) and (29), , and the sign of the square root is related with the sign of , so (32) can be written as:

If is in the first quadrant or the third quadrant, since it is not within , then:

If can arrive at the origin in two steps, every state and control value must satisfy (8), which are the following relationships:

Setting and , (35) can be written in matrix form as:

Then,

There is no harm in defining the control synthesis function as , and based on all these analyses, a time optimal synthesis function can be obtained by the following theorem.

Theorem 1. For an arbitrary in the phase plane, it can reach the origin quick under the fastest synthesis function obtained in the following procedure:

(1)If , cannot arrive at the origin within two steps, else jump to (5).(2)If , , .(3)Calculate the number of steps it needs for to reach the line , . If , , .(4)If can reach in one step, , .(5) can reach the origin in two steps, .(6)The end.

A nonlinear discrete tracking differentiator based on the previous mentioned synthesis function is described as follows. For a given signal sequence , apply

Then, approaches and approaches the differential of . In (38), is named as filtering coefficient and is named as sampling steps. The tracking differentiator above is a piecewise linear function, its structure is clearer, computation is simpler, and the fastest trajectories are shown in Figure 5.

From Figure 5, when the trajectory enters , it may fall in the area bounded by or ; then, under proper control value according to algorithm, the trajectory will move to the origin along .

Figure 5: Switching curves, boundary lines, and fastest trajectories.

4. Numerical Simulations of Discrete Tracking Differentiator

In this part, numerical simulations are given to verify the effectiveness of the proposed tracking differentiator. Firstly, a simulation test for given signals are performed, then a test in maglev test bench is done.

4.1. Given Signal Tests

For a given signal sequence , a tracking differentiator in the form of (38) is constructed and then the tracking numerical simulation results and differential numerical simulation results of a sinuous signal and a square wave signal are given.

4.1.1. Sinuous Test Signal

4.1.2. Square Wave Signal

The simulation result for a sinuous test signal is shown in Figure 6 while the simulation result for a square wave signal is shown in Figure 7. Both simulation results show that the above tracking differentiator can track the sinuous signal and square wave signal without overshoot and the differential of the given signal can be acquired properly.

Figure 6: Tracking and differential of sinuous signal.

Figure 7: Tracking and differential of square wave signal.

4.2. Application in Vibration Suppression of Maglev System

Vibration is a vital problem for maglev train which originates from the elastic property of the track and the coupling of maglev train and track [19, 20]. Many research studies have been carried out about the suppression methods of train-track coupling vibration. But most of them are too complicated and rely heavily on the precise system model. It is verified in practice that a signal processing method is sometimes simple and effective. The common feedback control method for levitation control in maglev system is as follows:

where is the levitation gap, is the current through the electromagnet windings, is the acceleration of the electromagnet, and is the reference gap value. , , , and are corresponding feedback gain and their values can be obtained via pole placement or the LQR method. This control law works well under the assumption that the track is rigid. However, the track itself vibrates under the effect of electromagnet force, and levitation gap is the difference between the displacement of track and the displacement of electromagnet. To get rid of this coupling, the differential of levitation gap signal can be used to restrain the coupling vibration. Considering this, a tracking differentiator based vibration suppression control scheme is designed in Figure 8. In this scheme, the feedback of gap differential is added to the control output, which results in the following:

where signal is obtained from the tracking differentiator and is the corresponding differential coefficient.

Figure 8: Levitation control diagram with tracking differentiator for vibration suppression.

To verify this, a simulation on a maglev test bench described by the model of (6) is made with the results shown in Figure 9. At first, the maglev system works under control law (41), then from t = 4 s to t = 6 s, control law (42) is applied with the differential signal obtained from an ordinary differentiator, and after t = 6 s, tracking differentiator is applied for extracting differential from noised signal. Their performances are labelled separately in different colors in Figure 9, with the red curve representing the control performance by control law (41), the blue curve representing the control performance by control law (42) without tracking differentiator, and the green curve representing the control performance by control law (42) with the proposed tracking differentiator. Even though control scheme (42) has a better vibration suppression ability, during t = 4 s to t = 6 s, the ordinary differentiator amplified the noise in the signal which makes the vibration still strong, whereas the tracking differentiator can suppress the noise more effectively. From this comparative simulation result, it can be seen that the proposed tracking differentiator is effective and superior in suppressing vibration.

Figure 9: Vibration suppression performance using tracking differentiator.

5. Conclusion

A discrete tracking differentiator is constructed by discrete system fastest control synthesis function, the boundary that control value adopts non-extreme value and the areas in which control value changes linearly are obtained. The fastest system control synthesis function is obtained by regionwise linearization. Compared with other types of tracking differentiator, the computation is simple and its performance is as good as nonlinear tracking differentiator. Numerical simulation and experiment results show that this new type of tracking differentiator has good performance in tracking, differentiating, and eliminating fluttering.

Data Availability

Data sharing is not applicable to this article as all datasets are hypothetical during the current study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the National Key Technology R&D Program of China (nos 2016YFB1200601 and 2016YFB1200602).

References

H. J.-Q. W. Wei, “Nonlinear tracking-differentiator,” Journal of Systems Science and Mathematical Sciences, vol. 2, p. 11, 1994.

H. Jinqing and Y. Lulin, “The discrete form of tracking-differentiator,” Aournal of Systems IciTnce and Mathematical Sciences, vol. 3, p. 2, 1999.

H. Jingqing, “The improvement of PID control law by using nonlinearity,” Information and Control, vol. 6, 1995.

H. Jingqing, “Nonlinear PID controller,” Acta Automatica Sinica, vol. 20, no. 4, pp. 487–490, 1994.

C. Chang, W. Chaozhu, and H. Jingqpng, “A new approach to estimate the moving parameters for maneuvering objects,” Journal of Astronautics, vol. 16, no. 1, pp. 30–34, 1995.

W.-G. Zhang and J.-Q. Han, “The application of tracking differentiator in allocation of zero,” Actc Automatici Sinica, vol. 27, no. 5, pp. 724–727, 2001.

J. ean, “From PID to active disturbance rejection control,” IEEE Transactions on Industrial Electronics, vol. 56, no. 3, pp. 900–906, 2009.

Y. Hvang and W. Zhang, “Development of active disturbance rejection controller,” Control Theory & Applications, vol. 19, no. 4, pp. 485–492, 2002.

X.-H. Wang, N.-Q. Chei, and Z.-Z. Yuav, “Nonlinear tracking-differentiator with high speed in whole course,” Control Theory & Applications, vol. 20, no. 6, pp. 875–878, 2003.

L.-Q. Wu and J.-Q. Han, “Active disturbance rejection controller scheme for the linear inverted pendulum,” Control Theory & Applications, vol. 21, no. 5, pp. 665–669, 2004.

M. Zhang and J. Wu, “Control system of renewable energy connected grid based on the auto-disturbances rejection control technology,” Control Theory & Applications, vol. 22, no. 4, pp. 583–587, 2005.

Y.-d. Xie and Z.-q. Long, “A high-speed nonlinear discrete tracking-differentiator with high precision,” Control Theory & Applications, vol. 26, no. 2, pp. 127–132, 2009.

Z.-z. Zhang, L.-h. She, L.-l. Zwang, and X.-l. Li, “Damp signal acquiring and fusing technology for middle-low speed permanent-electro maglev trains,” Journal of the China Railway Society, vol. 4, p. 8, 2011.

X. Song, L. Zhiqiang, H. Ning, and C. Wensen, “A high precision position sensor design and its signal processing algorithm for Maglev train,” Sensors, vol. 12, no. 5, p. 5, 2012.

J. S.-t. H. Zhong-sheng, “An improved model-free adaptive control for a class of nonlinear large-lag system,” Control Theory & Applications, vol. 4, p. 8, 2008.

Y. Guo, W. Wu, and K. Tang, “A new inertial aid method for high dynamic compass signal tracking based on a nonlinear tracking differentiator,” Sensors, vol. 12, no. 6, pp. 7634–7647, 2012.

Y. Zhou, Q.-L. Wang, and D.-H. Qiu, “Cascade form of tracking differentiator with high-order indiscrimination degree,” Transactions of Beijing Institute of Technology, vol. 31, no. 9, pp. 1058–1061, 2011.

A. Xiong and Y. Fan, “Application if a nonlinear velocity estimator in photoelectronic tracking system,” Electronics Optics & Control, vol. 11, 015 pages, 2010.

J. Li, J. Li, D. Zhou, P. Cui, L. Wang, and P. Yu, “The active control of maglev stationary self-excited vibration with a virtual energy harvester,” IEEE Transactions on Industrial Electronics, vol. 62, no. 5, pp. 2942–2951, 2015.

D. Zhou, P. Yu, L. Wang, and J. Li, “An adaptive vibration control method to suppress the vibration of the maglev train caused by track irregularities,” Journal of Sound and Vibration, vol. 408, pp. 331–350, 2017.