Adaptive vibration control of the electromagnettrack coupled high frequency resonance for an urban maglev system
 Authors: Zhou D.^{1}^{,2}, Yu P.^{1}, Li J.^{1}, Cui P.^{1}, Song M.^{1}

Affiliations:
 National University of Defense Technology
 State Key Laboratory of Functional Materials for Informatics Shanghai Institute of Microsystem and Information Technology Chinese Academy of Sciences
 Issue: Vol 4, No 2 (2018)
 Pages: 92106
 Section: Original papers
 URL: https://journals.ecovector.com/transsyst/article/view/10183
 DOI: https://doi.org/10.17816/transsyst20184292106
Cite item
Full Text
Abstract
This paper is a study of the electromagnettrack coupled high frequency resonance that frequently occurs in the urban maglev systems, it includes the following points:
Aim: The purpose of this study is to investigate the principle underlying the high frequency resonance occurs between the maglev train and the track, and to develop an appropriate vibration control algorithm which can be applied in the levitation controller, such that the resonance can be eliminated when the maglev train travels along the track.
Materials and methods of the studies: In this paper, the model of the electromagnettrack coupled system is firstly established, in which some special cases, which correspond to the situations when the screws that fasten the Frail to the sleepers are fatigue, or the stiffness of the rubber plates beneath the sleepers weaken for temperature reasons, are studied; and the reason that leads to the coupled resonance are explained as well. Secondly, an adaptive vibration control algorithm, which consists of a vibration observer and a tunable adaptive filter, is designed to suppress the high frequency electromagnettrack coupled resonance.
Results: Using this algorithm, when the train arrives at the spots where the coupled resonance may occur, the vibration observer will detect the occurring of the vibration and estimates its frequency, and then activate the adaptive filter and tune it to absorb the vibration.
Conclusion: The test indicates that this algorithm is capable of tuning itself to handle the unpredictable coupled resonance that occurs along the track, and it is simple and can be easily integrated into the levitation control code in a digital levitation control system.
Full Text
INTRODUCTION
On December 30th, 2017, the Beijing S1 urban maglev line began its commercial operation after years of argument and construction ever since 2011. Earlier in May, 2016, the first low speed urban maglev line in China – the Changsha Maglev Express, began to operate. The successful operation of these two maglev lines in China, together with the commercial maglev lines in Japan and Korea, indicate that the low speed urban maglev technology is tending to be mature, and its spread around the world is predictable.
However, as the development of test lines and the construction of commercial lines around China, more and more problems that related to the levitation control system of the train have been discovered. For example, the stability of the levitation control system under special track conditions, the reliability of the maglev controller, and the sensors under severe circumstances during the operation, etc. Another common problem for the urban maglev system may have been encountered – the vehicleguideway coupled resonance (or coupled selfexcited vibration in some literature). The coupled resonance occurs between the levitation control system and the track (or girder) when the track conditions are poor or the levitation control system is not well tuned. It is reported that the coupled resonance problem might be a fatal problem if it occurs on some lightweighted girders and cannot be well solved [1].
The maglev vehicleguideway coupled resonance problem can be further divided into two major categories: the vehiclegirder coupled selfexcited vibration, and the electromagnettrack coupled resonance. The former generally occurs on an elevated girder when the train is suspending in a standstill. It appears that the girder begins to vibrate and it can be clearly felt inside the train. The frequency of the girder vibration is generally below 30 Hz, which is determined by the natural frequency of the girder, but the energy of the vibration is large enough to cause levitation failure or even damages to the system structures. These problem has been noticed and discussed by Zhou et al.[2] The electromagnettrack coupled resonance appears as the vibration of the steel rails when the electromagnet is suspending at some specified locations. It may appear not as severe as the vehiclegirder coupled selfexcited vibration, but the frequency is generally above 50 Hz, thus it can be clearly heard by passengers onboard the train or people inside the rooms beneath the tracks. To meet the requirements of the commercial line, such as the Beijing S1 maglev line, the electromagnettrack coupled high frequency resonance must be solved, and it is also the major problem to be discussed in this paper.
The electromagnettrack coupled resonance problem has also been noticed by Zhou et al., as discussed in [3–5], in which several methods to eliminate the coupled resonance has been tried. For example, in [3], an LMS cancellation algorithm has been proposed to generate a signal which is negative to the vibratory component mixed in the closed loop of the levitation control system, and different method for producing reference signals for the LMS algorithm have been proposed in these literature. In [5], the amplitude and frequency of the high frequency resonance can even be precisely computed, which agrees quite well to experiments. More detailed and comprehensive discusses of the coupled vibration problems can be found in the dissertation of Zhou [6]. In this paper, considering that the track conditions may differ as the train travels at different segments of the commercial line, and the amplitude, as well as the frequency, may vary as the train travels; therefore, a more reliable and robust method to eliminate the trackinduced high frequency resonance needs to be developed. In this paper, the principle of the high frequency resonance is firstly introduced; and based on this theory, a new vibration control structure, which consists of a vibration observer and a tunable adaptive FIR (finite impulse response) filter, is designed to suppress the electromagnettrack coupled resonance. This adaptive structure can tune itself to cancel the vibration when the high frequency resonance tends to appear, and it is suitable to work while the train is traveling along the maglev lines.
MODELING OF THE ELECTROMAGNETTRACK COUPLED SYSTEM
During the test runs of the maglev trains, it is discovered that the electromagnettrack coupled high frequency resonance always occurs at the locations where the screws which fasten the Frail to the sleepers are fatigue or slack; another case occurs when the stiffness of the rubber plates beneath the sleepers are weaken for some reasons like temperature. Here, a model with a not well supported track is discussed, as shown in Fig. 1. The latter case mentioned above may also be characterized as a weakly supported track.
Fig. 1. A track model with slack screws connecting a sleeper and the rail
For the low speed maglev tracks, the rail is long and slim, thus it can be deemed as an EulerBernoulli beam with multiple supporters. A simplified model of the rail in each side of a track is shown in Fig. 2, in which each supporter is modeled as a combination of a linear spring and a torsal spring. The linear spring corresponds to the flexibility brought by the sleeper and the rubber plate beneath the sleeper; while the torsal spring characterizes the torsal rigidness of the sleeper since the rail are mounted to each sleeper by four hexagon socket screws. However, when the four screws are loose or fatigue, the stiffness of these springs will decrease or even disappear, leading to a weakly supported beam (say, the kth support shown in Fig. 2). Fig. 2 also demonstrates the relative position of an electromagnet beneath the rail.
Fig. 2. Simplified model of the rail with multiple supporters
In Fig. 2, suppose that E, I, and ρ are the Young’s modulus of the track material, the moment of inertia of the track cross section, and the mass per unit length of the track, respectively; ks_{n} and Ts_{n} represent the equivalent stiffness of these two springs at the nth sleeper, respectively; and L_{n} is the span length of the nth span. Then, for the EulerBernoulli beam, the vertical displacement of the rail, y(x,t), satisfies:
$EI\frac{{\partial}^{4}y\left(x,t\right)}{\partial {x}^{4}}+\rho \frac{{\partial}^{2}y\left(x,t\right)}{\partial {t}^{2}}=f\left(x,t\right)$ (1)
where f(x,t) is the electromagnetic force distributed on the rail. Note that for the urban maglev system, f(x,t) can be regarded as a evenly distributed force. Firstly consider the free vibration of the beam, and temporarily set f(x,t) = 0. Suppose that the free vibration of the beam can be expressed as y(x,t) = f(x)sinωt, then according to Eq.(1), we have
$\frac{EI}{\rho}\frac{{\partial}^{4}\varphi \left(x\right)}{\partial {x}^{4}}{\omega}^{2}\varphi \left(x\right)=0$ (2)
Eq. (2) indicates that the form of the solution for the mode shape,f(x), can be described as a combination of trigonometric and exponential functions, thus for the nth span of the beam, it can be written as:
${\varphi}_{\left(n\right)}\left(x\right)={C}_{n1}\mathrm{sin}\left(\lambda x\right)+{C}_{n2}\mathrm{cos}\left(\lambda x\right)+{\u0421}_{n3}sh\left(\lambda x\right)+{C}_{n4}ch\left(\lambda x\right)$ (3)
where λ is the spatial frequency. Substituting Eq. (3) into Eq. (2), yields:
${\omega}^{2}=EI{\lambda}^{4}/\rho $ (4)
Referring to Fig. 2, it can be seen that for the Nth span beam, due to the effects of the springs, the bending moment and the shear force of the beam at the left terminal of the rail satisfy the following two equations, respectively:
$EI\frac{{\partial}^{2}{\varphi}_{\left(1\right)}\left(0\right)}{\partial {x}^{2}}=T{s}_{1}\frac{\partial {\varphi}_{\left(1\right)}\left(0\right)}{\partial x}$ (5)
$EI\frac{{\partial}^{3}{\varphi}_{\left(1\right)}\left(0\right)}{\partial {x}^{3}}=k{s}_{1}{\varphi}_{\left(1\right)}\left(0\right)$ (6)
For the intermediate supports, the balance equations for the bending moment and the shear force of the beam can be expressed as follows:
$EI\frac{{\partial}^{3}{\varphi}_{\left(n\right)}\left(x\right)}{\partial {x}^{3}}{}_{x={L}_{n}}=k{s}_{n+1}{\Phi}_{\left(n\right)}\left(\u0445\right){}_{x=}{L}_{n}+EI\frac{{\partial}^{3}{\varphi}_{(n+1)}\left(x\right)}{\partial {x}^{3}}{}_{x=0}$ (7)
$EI\frac{{\partial}^{2}{\varphi}_{\left(n\right)}\left(x\right)}{\partial {x}^{2}}{}_{x={L}_{n}}=T{s}_{n+1}\frac{\partial {\varphi}_{(n+1)}\left(x\right)}{\partial x}{}_{x=0}+EI\frac{{\partial}^{2}{\varphi}_{(n+1)}\left(x\right)}{\partial {x}^{2}}{}_{x=0}$ (8)
At the right terminal, the bending moment and the shear force of the beam satisfy:
$EI\frac{{\partial}^{2}{\varphi}_{\left(N\right)}\left(x\right)}{\partial {x}^{2}}{}_{x={L}_{N}}=T{s}_{N+1}\frac{\partial {\varphi}_{\left(N\right)}\left(x\right)}{\partial x}{}_{x={L}_{N}}$ (9)
$EI\frac{{\partial}^{3}{\varphi}_{\left(N\right)}\left(x\right)}{\partial {x}^{3}}{}_{x={L}_{N}}=k{s}_{N+1}{\Phi}_{\left(N\right)}\left(x\right){}_{x={L}_{N}}$ (10)
Note that the continuities of the displacement and slope of the beam at the intermediate supports require
f_{(n)}(L_{n}) = f_{(n}_{+1)}(0) (11)
f'_{(n)}(L_{n}) = f'_{(n}_{+1)}(0). (12)
Combing Eqs. (5)–(12) together, there will be totally 4×N equations, and they can be written in a matrix form:
${\mathbf{A}}_{\mathbf{m}}\mathbf{C}=0$ (13)
where A_{m} is a 4N×4N matrix, and C is the coefficients of the mode shapes described by Eq. (3), namely, . By solving Eq. (13), the unknown coefficients, as well as the spatial frequency of the beam, λ, can be obtained.
The solution of Eq. (1) can be described as a combination of all the modal shapes obtained above, namely:
$y\left(x,t\right)=\sum _{k=1}^{\infty}{q}_{k}\left(t\right){\varphi}_{k}\left(x\right)$ (14)
Here, y_{k}(t) is the generalized displacement of the kth order mode. Substituting Eq. (14) into Eq. (1), multiplying both sides of the resultant equation by f_{k}(x), and integrating both sides along the beam, noticing the mode orthogonality condition for linear beams [7], one obtains
${\ddot{q}}_{k}\left(t\right)+{\omega}_{k}^{2}{q}_{k}\left(t\right)=\frac{{F}_{k}\left(t\right)}{{M}_{k}}$ (15)
where ${M}_{k}=\rho {\int}_{0}^{L}{\varphi}_{k}^{2}\left(x\right)\mathrm{d}x$, ${F}_{k}\left(t\right)={\int}_{0}^{L}f\left(x,t\right){\varphi}_{k}\left(x\right)\mathrm{d}x$, and $L=\sum _{k=1}^{N}{L}_{k}$.
Using the Laplace Transform, Eq. (15) can be rewritten as
${q}_{k}\left(s\right)=\frac{{F}_{k}\left(s\right)}{\left({s}^{2}+{\omega}_{k}^{2}\right){M}_{k}}$ (16)
Therefore, the vertical displacement of the beam at a specified location, x_{0}, can be given by
${y}_{k}\left({x}_{0},s\right)=\frac{{F}_{k}\left(s\right){\varphi}_{k}\left({x}_{0}\right)}{\left({s}^{2}+{\omega}_{k}^{2}\right){M}_{k}}$ (17)
For the electromagnetic force, f(x,t) = F_{m}(t)/l_{0}, where F_{m}(t) is the total electromagnetic force produced by the coil pair, and l_{0} is the length that the electromagnetic force covers. According to Fig. 2, it can be seen that
${F}_{k}\left(s\right)=\frac{{F}_{m}\left(s\right)}{{l}_{0}}\left\{\begin{array}{l}{\int}_{{x}_{0}}^{{x}_{0}+{l}_{0}}{\varphi}_{k}\left(x\right)\mathrm{d}x,\mathrm{if}{x}_{0}+{l}_{0}\le L\\ {\int}_{{x}_{0}}^{L}{\varphi}_{k}\left(x\right)\mathrm{d}x,\mathrm{if}{x}_{0}+{l}_{0}>L\end{array}\right.$ (18)
Therefore, Eq.(17) can be rewritten as another form:
${y}_{k}\left({x}_{0},s\right)=\frac{{F}_{m}\left(s\right){g}_{k}\left({x}_{0}\right)}{\left({s}^{2}+{\omega}_{k}^{2}\right)}$ (19)
Here, g_{k}(x_{0}) can be deemed as the kth order mode gain of the beam, which is a coefficient that reflects the deformation response of the beam at x_{0} when the beam is subjected to a unit force, F_{m}. The larger g_{k}(x_{0}) is, the larger will be the displacement of the rail at x_{0}. As a demonstration, for the Beijing S1 low speed urban maglev system, l_{0} = 1.32 m, L_{n} = 1.2 m, ρ = 130 kg/m, EI = 1.93×10^{6} N m^{2} (the last two parameters of are obtained by calculation); the stiffness of the linear and torsal springs are supposed to be: ks_{n} = 1×10^{8} N/m, and Ts_{n} = 7×10^{10} Nm/rad. For the strongly supported rail, all the sleepers are well fastened to the rail, thus the parameters for all the springs are identical with each other. Using the procedure presented above, the first three order mode shapes can be computed, as shown in Fig. 3, and their corresponding mode frequencies are 122.2 Hz, 129.3 Hz, and 141.5 Hz, respectively. This result is rather close to finite element analysis. Meanwhile, the mode gains for these three vibration modes can be obtained as well, which are shown in Fig. 4.
As a comparison, consider the case when the screws of one sleeper slack. Specifically, suppose that the weak support appears at the middle sleeper (see Fig. 1). In this case, the stiffness of the linear spring and the torsal spring will be prominently weakened, which are supposed to be 1/10 of the stiffness of other well fastened supports. The resultant vibration mode shapes as well as mode gains of the weakly supported rail can be computed, as shown in Fig. 5 and Fig. 6, respectively. Comparing Fig. 5 with Fig. 3, it can be seen that, due to the decrease of support stiffness in the middle of the rail, the shape of the first order mode shape changes dramatically – the displacement of the rail in the middle is more obvious than at other supports, which indicates that when the train moves here, the rail displacement will significantly increase. Also, Fig. 6 shows a prominent increase of mode gain at the middle for the first order vibration mode. According to Eq. (19), this suggests that the vertical response of the rail at the middle will be larger than at other locations when the train passes the rail. This is a fundamental reason that leads to electromagnettrack coupled resonance when the track support weakens.
Fig. 3. The first three order mode shapes of the wellsupported rail
Fig. 4. Mode gains for the first three order vibration modes of the wellsupported rail
Fig. 5. The first three order mode shapes of the wellsupported rail
Fig. 6. The first three order mode shapes of the weaklysupported rail
Another significant change as a result of different rail support conditions is the mode frequency. For comparison, the mode frequencies, as well as mode gains when x_{0} = 3.3 m, for these two cases are listed in Table 1. It can be seen that the fundamental mode frequency drops significantly as the middle support becomes weak, while the other two higher mode frequencies only decrease a bit. This is another reason why the coupled resonance tends to occur when the rail is not well supported, since the lower the mode frequency is, the higher is the open loop gain, as is analyzed by Zhou et al. [2–6].
Table 1. Mode frequencies and mode gains for the normal and the weakly supported rail
Mode order  Mode frequency (Hz)  Mode gain (´10^{–3})  
Normal rail  Weak rail  Normal rail  Weak rail  
1  122.2  76.4  1.82  4.53 
2  129.3  126.1  –0.10  –0.31 
3  141.5  133.2  1.51  0.10 
For the levitation control system shown in Fig. 2, the levitation gap δ, which is the vertical distance between the electromagnet and the track measured by the gap sensor, can be described as
$\delta \left(t\right)={y}_{m}\left(t\right){y}_{\mathrm{T}}\left(t\right)$ (20)
where y_{m} is the vertical displacement of the electromagnet, y_{T} is the vertical deformation of the track at x_{0}, which can be obtained by the superposition of all vibration modes using Eq. (14). To keep the levitation gap at a desired value, say, z_{0}, a controller is needed to eliminate the levitation gap error. The PID cascade controller, which has been employed in [8], is also adopted here, which gives
${i}_{e}\left(t\right)={k}_{p}\left[\delta \left(t\right){z}_{0}\right]+{k}_{d}{\dot{y}}_{m}\left(t\right)+{k}_{i}{\int}_{0}^{t}\left[\delta \left(\tau \right){z}_{0}\right]d\tau $ (21)
where i_{e} is the desired current through the electromagnet coil. To achieve rapid response of the current, a current loop is utilized, and the control voltage is given by
$u\left(t\right)={k}_{c}\left[{i}_{e}\left(t\right)i\left(t\right)\right]$ (22)
Here, k_{p}, k_{d}, k_{i}, and k_{c} are constants which can be chosen to achieve desired dynamic performance of the control system. The current through the electromagnet coil satisfies [3]:
$\frac{u\left(t\right)}{2}=i\left(t\right)R+\frac{{\mu}_{0}A{N}_{c}^{2}}{2\delta \left(t\right)}\dot{i}\left(t\right)\frac{{\mu}_{0}A{N}_{c}^{2}i\left(t\right)}{2{\delta}^{2}\left(t\right)}\dot{\delta}\left(t\right)$ (23)
Here, A is the area of the electromagnetic pole, μ_{0} is the space permeability, R is the resistance of the electromagnet coil, and N_{c} is the number of turns of the coil. The electromagnetic force of a levitation unit is:
${F}_{m}\left(t\right)=\frac{{\mu}_{0}A{N}_{c}^{2}}{2}{\left(\frac{i\left(t\right)}{\delta \left(t\right)}\right)}^{2}$ (24)
The vertical movement of the electromagnet can be described by:
${m}_{1}\frac{{\mathrm{d}}^{2}{y}_{m}\left(t\right)}{\mathrm{d}{t}^{2}}={F}_{m}\left(t\right)+{m}_{1}g$ (25)
in which m_{1} is the equivalent mass of the levitation payload. Detailed discussion of the stability of the electromagnettrack coupled system has been presented in previous published literature [3–5], which will not be repeated in this paper. Instead, as a demonstration, the stability of the coupled system discussed here will be investigated through numerical simulations. For example, if the parameters of the control system are chosen according to Table 2, and suppose that the track is well supported, the responses of the levitation gap and the displacement of the track can be obtained, which are shown in Fig. 7. It can be seen from Fig. 7(a) that, under the chosen parameters, the overdamped transient response of the levitation gap reflects a good performance of the levitation controller; more importantly, in this case, the vibratory response of the track tends to disappear, as shown in Fig. 7(b). This is a desired performance for a commercial maglev system.
Table 2. Parameters of the levitation control system
Parameter  Value  Parameter  Value 
k_{p}  4000  m_{1}  820 kg 
k_{d}  100  μ_{0}  4π×10^{7} H/m 
k_{i}  0.01  A  0.01848 m^{2} 
k_{c}  100  N_{c}  360 
l_{0}  1.32 m  R  0.5 Ohm 
Fig. 7. Responses of the levitation gap and the displacement of the track under the well supported track
However, for the weakly supported track shown in Fig. 1, the fundamental frequency of the track drops dramatically, which greatly changes the stability of the electromagnettrack coupled system. Under the same parameters, the responses of the levitation gap and the track are shown in Fig. 8. Compared with Fig. 7, it is obvious that the electromagnettrack coupled system becomes unstable, and high frequency coupled resonance occurs. As the amplitude of the vibration grows, the noise arise which can be clearly heard by passengers aboard the train. This case has been encountered in the Tangshan low speed maglev test line, and it should be eliminated to prevent potential damage to track infrastructures.
Fig. 8. Responses of the levitation gap and the displacement of the track under the weakly supported track
ADAPTIVE VIBRATION CONTROL SCHEME FOR THE ELECTROMAGNETTRACK COUPLED SYSTEM
The selfexcited vibration caused by weakly supported tracks seems inevitable for a practical maglev track since the fatigue or slack of screws may be unavoidable as a result of seasonal thermal expansions and contractions, as well as ground settlement. To eliminate the encountered selfexcited vibration while the train is traveling, some special aspects should be taken into account: first, the frequency of the selfexcited vibration is not a constant, but is variable due to the support conditions; second, the response of the vibration control scheme should be as rapid as possible; third, the vibration control scheme should not disturb the stability of the levitation controller, or at least it should not distinctly decrease the performance of the levitation control system.
Based on these considerations, an adaptive vibration control scheme, which is capable of tuning its parameters to cope with different track situations, is essential. Here, an adaptive vibration controller is proposed, as shown in Fig. 9. The block "CTL" represents the levitation controller which was described by Eqs. (20)–(22), "MAG" is the electromagnet whose output includes the current and the electromagnetic force; "H(s)" is the transfer function of the track. The mechanism inside the dashed box is the designed adaptive vibration controller which includes two major parts: the frequency estimator which was designated by a block "FE", and the adaptive FIR filter with an adaptive tuning mechanism, "ADJ". The basic idea of this scheme can be generalized as follows: as the selfexcited vibration occurs, the vibratory component mixed in the expected current, i_{e}_{0}, will exceed a certain value, and its frequency can then be estimated by the observer "FE"; a reference signal x, which is a sinusoid whose frequency equals the estimated frequency, ω_{r} , is generated for the FIR queue, X. The output of the FIR filter, y_{f}, is then subtracted from i_{e}_{0}, aiming to cancel the vibratory components in i_{e}_{0}. To achieve precise cancellation, the weights of the FIR filter, W, are adjusted by the adaptive mechanism ADJ, using the Least Mean Square (LMS) algorithm. If well tuned, the output of the cancellation, i_{e}, will not contain vibratory component, which indicates that the electromagnetic force, F_{m}, will not include vibratory excitation, either. Then, it can be inferred that the vibratory energy of the track will be dissipated as time elapses, and the selfexcited vibration can thus be eliminated.
Fig. 9. Adaptive vibration control scheme combined with the maglev control system
From Fig. 9, it can be seen that
${y}_{f}={\mathbf{W}}^{T}\mathbf{X}$ (26)
${i}_{e}={i}_{e0}{y}_{f}$ (27)
The purpose of the FIR filter is to produce a cancellation signal which is capable of eliminating the vibratory signal included in i_{e}_{0}, thus its weights needs to be tuned timely. The LMS algorithm is adopted here, and the weights of the FIR filter are updated according to the following rules:
${w}_{n}\left(k+1\right)=\left(1\alpha {\mu}_{1}\right){w}_{n}\left(k\right)+2{\mu}_{1}e{x}_{n}\left(k\right)$ (28)
where μ_{1} is the convergence coefficient of the LMS algorithm, e is the cancellation error of the vibratory component in i_{e}, and α is a leakage coefficient that ameliorates the potential overflows of the weights caused by the accumulation of quantization errors.
To provide the FIR filter with a correct reference signal, a frequency estimator, FE, is utilized. Here, the frequency estimator developed by Regalia [9] is adopted here. This algorithm has been adopted by some other researchers, as in [10], which can be written as:
$\left\{\begin{array}{l}{\dot{x}}_{1}={x}_{2}\\ {\dot{x}}_{2}=2\zeta {\omega}_{r}{x}_{2}{\omega}_{r}^{2}{x}_{1}+k{u}_{hf}\\ {\dot{\omega}}_{r}=g\left(k{u}_{hf}2\zeta {\omega}_{r}{x}_{2}\right){x}_{1}\end{array}\right.$ (29)
where u_{hf} is the input of the estimator. However, Fig. 9 is only a demonstration of the basic principle of the adaptive vibration controller. In practice, more considerations are required to keep the adaptive FIR filter work properly. It can be noted that i_{e}_{0} contains not only the vibratory component which is associated to the vibration of the track, but also the DC (direct current) component which produces a biased current to levitate the payload; the latter, however, would cause the error signal for the LMS algorithm, e, superimposed with a DC component, which would prevent the LMS algorithm from converging. To solve this problem, high pass filters, as well as a compensation method introduced in [4], may also be applied here. For simplicity, it is not repeated in this paper.
Fig. 10 shows the test result of this adaptive vibration control method. The simulation is undertaken when all the parameters of the levitation system are chosen following the example in the previous section. For comparison, the adaptive vibration controller (AVC) was disabled for the first two seconds. As discussed earlier, due to the weak support of the track, selfexcited vibration occurs, and the amplitude of the vibration grows with time. Fig. 10(b) shows the result of the frequency estimation, from which it can be seen that as the amplitude of the vibration exceeds a small value (corresponding to 0.9s), the frequency estimator begin to work and the result converges to the actual vibration frequency, which is 77 Hz. Also, it can be seen from Fig. 10(c) that, the weights of the FIR filter are updated timely, and they are tuned in different manners as the reference frequency changes.
Fig. 10. Test result of the adaptive vibration control method. (a) the levitation gap; (b) estimated frequency; (c) selected weights of the adaptive FIR filter
Then, when time equals 2 s, the AVC is activated and the cancellation begins to take effect. It can be seen that the vibration begin to attenuate quickly afterwards. The weights of the adaptive filter also attenuate as the amplitude of the vibration decreases. The test result is satisfactory since the vibration of the result can be totally eliminated, and the electromagnettrack coupled system becomes stable.
CONCLUSION
In this paper, the principle of the trackinduced coupled resonance problem in the urban maglev system is discussed. The major reason that lead to the electromagnettrack coupled vibration is analyzed, which reveals that the weak support of tracks may dramatically decrease the fundamental frequency of the track, whilst at the same time increase the mode gain of the fundamental vibration mode. These two factors destroy the conditions for the stability of the coupled system, leading to a selfexcited vibration. Aiming to suppress the selfexcited vibration when a maglev train travels, an adaptive vibration controller, which includes a frequency estimator used to observe the frequency, and an adaptive FIR filter which produces a cancellation signal to neutralize the vibratory component mixed in the control path, is presented. The LMS algorithm is employed to tune the weights of the FIR filter, which achieves a satisfactory performance of convergence. The test result indicates that the proposed adaptive vibration controller is capable of suppressing the electromagnettrack coupled resonance.
ACKNOWLEDGEMENT
This work was supported by the National Key Development Program of the 13th Fiveyear plan [grant number 2016YFB1200601B12]; the Research Plan Program of National University of Defense Technology [grant number 2K160202]; the Opening Foundation of the State Key Laboratory of Functional Materials for Informatics [grant number SKL201707]; and National Nature and Science Foundation of China [grant number 11302252].
About the authors
Danfeng Zhou
National University of Defense Technology; State Key Laboratory of Functional Materials for Informatics Shanghai Institute of Microsystem and Information Technology Chinese Academy of Sciences
Author for correspondence.
Email: zhoudanfeng@maglev.cn
ORCID iD: 0000000157076530
PhD, Lecturer
China, Changsha; ShanghaiPeichang Yu
National University of Defense Technology
Email: lofter@163.com
PhD, Lecturer
China, ChangshaJie Li
National University of Defense Technology
Email: lj@maglev.cn
PhD, Professor
China, ChangshaPeng Cui
National University of Defense Technology
Email: cuipeng@maglev.cn
PhD, Lecturer
China, ChangshaMengxiao Song
National University of Defense Technology
Email: songmengxiao14@nudt.edu.cn
PhD, Student
China, ChangshaReferences
 Zhou DF, Hansen CH, Li J, Chang WS. Review of coupled vibration problems in EMS maglev vehicles. International Journal of Acoustics and Vibration. 2010;15(1):1023. doi: 10.20855/ijav.2010.15.1255.
 Zhou DF, Hansen CH, Li J. Suppression of maglev vehiclegirder selfexcited vibration using a virtual tuned mass damper. Journal of Sound and Vibration. 2011;330:883901. doi: 10.1016/j.jsv.2010.09.018.
 Zhou DF, Li J, Hansen CH. Application of least mean square algorithm to suppression of maglev trackinduced selfexcited vibration. Journal of Sound and Vibration. 2011;330:57915811. doi: 10.1016/j.jsv.2011.07.021.
 Zhou DF, Li J, Zhang K. An Adaptive Control Method to Suppress the Maglev TrackInduced SelfExcited Vibration. Proceedings of the 2011 International Conference on Consumer Electronics, Communications and Networks; 2011; Xian Ning, China. p. 4723–4727. doi: 10.1109/CECNET.2011.5768485.
 Zhou DF, Li J, Zhang K. Amplitude control of the trackinduced selfexcited vibration for a maglev system. ISA Transactions. 2014; 53(5):1463–1469. doi: 10.1016/j.isatra.2013.12.016.
 Zhou DF. Research on SelfExcited Vibration Problems in EMS Maglev VehicleGuideway Coupled System [dissertation]. Graduate School of National University of Defense Technology; 2011 (in Chinese).
 Svetlitsky VA. Engineering Vibration Analysis (Worked Problems 2). Berlin: SpringerVerlag; 2004. doi: 10.1007/9783540409700.
 Zhou DF, Yu PC, Wang LC, Li J. An adaptive vibration control method to suppress the vibration of the maglev train caused by track irregularities. Journal of Sound and Vibration. 2017;408:331350. doi: 10.1016/j.jsv.2017.07.037.
 Regalia PA. An improved latticebased adaptive IIR notch filter. IEEE Trans. Signal Process. 1991;39(9):21242128. doi: 10.1109/78.134453.
 Bodson M, Douglas SC. Adaptive algorithms for the rejection of sinusoidal disturbances with unknown frequency. Automatica. 1997;33(12):22132221. doi: 10.1016/s00051098(97)001490.