Robust nonsingular sliding mode control of the maglev train system: case study

This paper deals with a new approach to explore the precise dynamic response of the maglev system train and its control. Magnetic-suspension systems are characterized by high nonlinearity and open-loop instability which are the core components of maglev vehicles. Firstly, we use the electromagnetics and mechanics laws to derive the mathematical expressions of the proposed maglev system. Analytical investigation and theoretical calculation show that for the specific values of the control system parameters, the maglev system train can be significantly improved. It points out that the inherent nonlinearity, the inner coupling, misalignments between the sensors and actuators, and external disturbances are the main issues that should be considered for maglev engineering. Secondly, a control strategy based on the precise model of a nonsing ular robust sliding mode control is designed to reduce the upper bound of both the uncertainty and interference of the sliding mode controller. This approach presents an added value compared to the new sliding control methods in terms of overshoot and speed of convergence which is designed to control the vertical position of the proposed system. By using rigorous mathematical transformation associated with the adaptation laws in the frequency domain, a sufficient condition is drawn for the stability of the dynamical error based on the Lyapunov theory. This allows us a great possibility for interpreting the operation of the maglev train system. Numerical results are presented to show the effectiveness of our proposed control scheme.


Introduction
Maglev systems attract great interest from engineers over the world due to their basic characteristics (contact-less property between the rail and the train, minimal maintenance cost, faster speed and so on) and their different applications such as high-speed, frictionless bearing and levitation of metal slabs during manufacturing [1]. It has a broad spectrum of applications in different domains in day-to-day life for example in transportation [2,3], bearing-less motors [4,5], bio-medical [6], industries [7,8], launching of rockets [9,10], in levitation of metals etc. A step-by-step design of an efficient suspension system is globally determined by an in-depth understanding of the train dynamics. It is also essential for the development of the control laws for the effective operation of the Maglev system. In order to be commercially viable, the Maglev train must have certain technical requirements to fully satisfy passengers. Therefore, it is essential to consider the dynamic as well as static stability as an integral part of Maglev design even for the control development.
From a viewpoint of the mechanics, the dynamics of the mechanical parts are inescapably affected by the electromagnetic properties. A ferromagnetic material suspended to an air spring, which is subject to the electromagnetic force created by the field of an electromagnet. The electrical current flowing along the magnetic circuit is a nonlinear function of the electromagnetic force. consequently, all motions induced in the ferromagnetic material are nonlinear. As a result, extensive dynamical researches have been carried out to achieve satisfactory performance specifications [11,12]. Most of these works essentially provide the theoretical basis for the development of the electromechanical systems and then give us great inspiration to further improve the Maglev system. Some approaches for the analysis of suspension systems were strongly associated with the vibration phenomena such as hopf bifurcation, limit cycle and so on. Recently, Yougang et al., investigated the modeling and the control of the low-speed maglev using fuzzy sliding-mode technique [13]. In their work, they used the Euler-Lagrange strategy to construct the nonlinear dynamic model. In the qualitative analysis viewpoint, Lee investigated the occurrence of the Hopf bifurcation and the impact of guide-way parameters on the Maglev system behavior [14]. Following the same idea, Haiyan and collaborators presented an appropriate levitation chassis analysis with a relevant robust position control design of the suspension system under nonlinear periodic disturbance [15]. They demonstrated that the stable suspension of the maglev vehicle is achieved by controlling individual EMS modules. More recently, Zhai and coworkers proposed a calculation and evaluation of the load performance of Maglev in medium-low speed [16]. They provided a load performance of the magnetic levitation system under different load using theoretical calculations associated with the finite element. We note that the typical modeling of the maglev system met the unceasing improvement and is more confronted to the basic requirements of engineering applications, but now-a-days, there are still many technical difficulties, especially for the analysis and design of a levitation system. Despite the numerous works on this topic, there is no relevant study that provides a systematic design procedure in the simpler manner of the Maglev system to the best of the author's knowledge. We notice that the improvement of the Maglev train dynamic is of paramount importance and can be achieved through the perfect knowledge of the behavior of the parameters.
It is worthwhile to recall that, the maglev system suffers from various control complexities such as inherent instability, one direction force input, nonlinearity modeling and exogenous disturbances [17,18]. Hence, an external controller is necessary to stabilize the overall system. Abundant control approaches of the maglev systems are based on the linear model which use the linear controller [19][20][21] and the references therein. However, this design model tracking performance deteriorates rapidly as deviation with respect to nominal operating point increases. Some modeling technics were proposed by researchers in order to linearize existing models or to develop a theory around the suppression of vibrations in the system. Most researches in literature deal with appropriate nonlinear controllers that control globally instead of locally the dynamic of the system in a real-time. The principal criterion to select the nonlinear controller is strongly based on its ability to give fast and efficient dynamic response i.e. lesser rise time, fast settling time, less peak value, minimum overshoot/undershoot and negligible steady-state error [22]. In addition, the nonlinear controller is designed to deal with the inherent perturbations on the system dynamic. The majority of the nonlinear control presents poor performance and is less efficient against the disturbances. Therefore, the major motivation behind this work is to design an appropriate nonlinear controller for stabilizing the maglev system with best dynamic and also to allow the guide magnets to keep the rail on track by overcoming the effect of disturbances in an efficient way.
Many contributions have been devoted to the development of control strategies of the maglev train. Yang et al., demonstrated that the magnetic levitation subjected to the uncertainty can be efficiently controlled via a particular dynamic surface proposed a robust position control of a magnetic levitation system via dynamic surface control technique of magnetic levitation system containing parameter uncertainty [23]. The position control strategy is constructed to maintain stable perturbed Maglev system. On the other hand, a robust estimator is designed to estimate the system uncertainties using a DSP-based sliding-mode control for an electromagneticlevitation system by Lee et al., [24]. These techniques present the advantages to be simple and consequently easy to be implemented, but they consider only a particular scenario, which may be conservative. Wiboonjaroen and Sujitjorn [25] exploited the linearized version of the Maglev model to derive a novel state-PID feedback controller. We notice that the results obtained penalize the performance of the controller due to the fact that the linear systems do not actually exist. The linearization of such systems inevitably leads to the loss of some important properties of the original systems. Recently, Bobtsov et al., proposed a state observer for sensorless control of maglev systems [26]. More recently, Sun and coworkers developed a minimum levitation unit via the RBF Neural Network for the Maglev systems [27]. The neural network enabled us to improve significantly the proposed strategy as regards time-delay. Although, the promising results obtained by these authors, it appears boring and bulking in the practical case. In [28], a nonlinear control method for the electromagnetic suspension system of the maglev train was designed by Junqi and Yuan. They prove that the adjustment of time delay parameters plays a guiding role in avoiding vibration. The non-adaptive time-delayed feedbacks controllers are designed to control the overall system which stands as a drawback since the feedback gains cannot be infinite. On the same way, You-gang Sun et al., presented a robust levitation control of maglev vehicles subject to time delay and disturbances: design and hardware experimentation [29]. Their works show an adequate design of controller based on the Riccati method and sliding-mode technology that can handle efficiency the time delays and disturbances in the control scheme. The basic requirement is the enhancement in the dynamic response of the system so that the proposed controller should have characteristics like: minimal overshoot/undershoot, faster convergence rate, minimal or zero steady state error and reduced chattering and robustness against external disturbance. To satisfy these requirements we choose the nonsingular sliding mode for designing our proposed controller.
Amongst huge control approaches, the application of sliding mode technic to control Maglev system are numerous [30,31] and give many advantages such as fast response, insensitivity to the matching external disturbances, robustness against uncertainties and good transient performance [32,33]. It's important to recall that sliding mode structure control has achieved success in various domains. However, from the literature viewpoint, there are only a few works that take into account these procedures magnetic suspension systems area. It is worth to recall that in physical Maglev systems, some significant constraints exist, such as unavailability of all the state variables, matched and/or unmatched uncertainties due to parameter variation, model simplification and/or external disturbances [33][34][35]. Due to these considerations, it is become a titanic task to design an appropriate controller with the specifics performance. An attempt to overcome perturbations in the system is to use the nonsingular sliding mode control which uses a single control inputs and allows the finite-time stabilization using without requiring the information of the uncertainties [36,37]. In addition, the singularity and the chattering are eliminated. Nevertheless, this topology of the controller has some weaknesses such as asymptotic stability with a slow response.
The scope of this paper falls into the electro-magnetic suspension technology. We consider a novel dynamical approach associated with a control problem of the nonlinear Maglev system. Facing the challenges previously established, the main contributions of this paper are three-fold: 1. Design and analyze the Maglev system in terms of laws and including stability theories, 2. Propose a robust control strategy for the Maglev train which can ensure the system to be uniformly bounded. 3. Propose an adaptive law to regulate the total system uncertainty, which was not available earlier.
The designing of the novel adaptive robust nonsingular sliding mode control reveals an important contribution of this paper. To the best of the author's knowledge, not of all these aspects of the Maglev have been simultaneously considered in the literature and then become challenging. In this study, we emphasize the need for improvement of Maglev system in their technical performance and the designing of the efficient control strategy which deserves recognition as a profound contribution.
The rest of the paper is organized as follows. Section 2 gives the problem description and dynamical analysis of the maglev system with some associated simulations. Section 3 deals with the designing of the nonsingular sliding mode control. The simulation investigation is presented in Sect. 4. Section 5 concludes the paper.

System modeling
Maglev train structure has two suspension systems: the upper car body between wheels and bogie and the lower between the bogie and car body. The suspension is necessary not only for isolating the car from vibrations and bumps but also to reduce the forces between wheels and rail. The bogie modeling is the key to the Maglev vehicle model [38]. In order to gain more insight into the behavior of the system, two majors factors are taken into account, the magnetic damping and stiffness of the air spring in the modeling which are the important parameters that determine the dynamic response and the stability of the Maglev system. These quantities in the equation of motion allow to gain more knowledge about the (nonlinear) dynamical behavior of the system for certain parameters and also enable the evaluation of their effects on the operation of the maglev [39,40]. For the sake of accuracy and for a judicious design of the Maglev train, various forms of active and passive damping mechanisms must be incorporated into the design to make the no steady-state operation stable [41]. The predominant approaches in the literature consider stiffness coefficient similar to the mechanical spring. The stability exists in the direction of positive stiffness if the reaction force acts to oppose perturbative displacements.
In order to analyze the dynamic of the maglev system, we consider the standard model structure presented in Fig. 1. This model is similar to those encountered in literature which consists of three subsystems: the spring system, guideway, and electromagnet system [29,42]. The specificity of our approach compared to other existing works is that the influence of the spring system and the friction are not ignored. For more details, we refer the reader to [43][44][45]. We neglected the magnetic flux leakage and the edge effect which cannot significantly affect the dynamical behavior of the Maglev. We also consider only the vertical motion of the electromagnet.
Remark 1 Keep the levitation system steady is more difficult since the maglev force is inherently unstable, which leads to the dynamic characteristics performing diversity and complexity. Two approaches have been presented in literature such as the dynamics characteristics of the maglev vehicle structure and other way is focusing on controller design [46,47]. The big challenge is to study the impact of each parameter on the dynamic behavior of the device and to detect those whose influence is significant for the Maglev design. Unfortunately, by referring to the theoretical and experimental aspects, most strategies although satisfactory encountered in the literature remain qualitative (Hopf, bifurcation …) which does not allow the rigorous establishment of the performance concepts of the Maglev system. An appropriate mathematical approach appears interesting and significantly helpful for us to understand the dynamical behavior of Maglev trains and to design relevant control laws.

Remark 2
The order of the Maglev will be high if all factors are taken into consideration which will be boring and the handing very difficult. For the sake of simplicity, we mention some important hypothesis considered in the analysis of this work: (1) Only vertical motion of the electromagnet is considered, namely, there is no displacement between the surface and the guideway. (2) the edge effect and the magnetic flux leakage are neglected, consequently, the flux is uniformly distributed in the air gap.
(3) the electromagnet develops an electromagnetic force which the origin is considered to be on the center of mass.
We derive the electromagnetic force as follows: where E m represents the magnetic field energy of volume V; e m is the magnetic field energy density; A is the magnetic pole area of the electromagnet; 0 is air permeability and F denotes the electromagnetic attraction force. It is well known that the magnetic flux density can be expressed as B = 0 Ni(t) 2z(t) and taking into account the Eq. (1), we finally derive the expression of the electromagnetic force generated by the electromagnet as follows: where z(t) is the magnet gap or the magnet displacement, N the number of turns and i(t) the current of the solenoid coil. We suppose that components of the Maglev train are rigid bodies. Using the electrical and mechanical laws, it is found that the mechanical and the electrical dynamics are described by the following equations Δi(t) and Δz(t) are the fluctuations of current and airgap when the Maglev moves on the guideway. At any point of the Maglev operation, the linearized electromagnetic force can be express as follows [48,49] where are the equivalent stiffnesses. Neglecting the fluctuation of the current in the electromagnetic force for the sake of simplicity, Eqs. (3) becomes Because very low damping is considered, the system is simplified to a constant damping coefficient: (̇z(t), z(t)) = . However, the equation is still nonlinear .1 Configuration of the Maglev system because of the nonlinear spring forces. u(t) is the control voltage acting on both ends of the electromagnet. m is the equivalent mass of the train body and the levitation chassis. K is the stiffness of the air spring. z e (t) the designed rating magnet gap z(t) , the mechanical resistance and a m the area of magnetic the pole while g is the gravity constant. Δz(t) = z(t) − z e is the magnet gab deviation. The simplified model of the slow medium maglev train can be obtained by introducing the new variables y and as z = y and i = y which lead to the following coupled systems: where ℏ = 0 N 2 A 2 denotes the force transmission coefficient.

Remark 3
It's easy to obtain the simplified model in [42,50] by setting ≈ 0 . In the following analysis, we recall that the mechanical resistance of the air spring is known and remains a fundamental characteristic of the Maglev described in [51][52][53].

Equilibrium point and stability analysis
The equilibrium points are solutions of the equations
If y 0 = 0 , then, 0 = ± √ −2 mg ℏ + K ℏ y e . We then obtain two equilibrium points Ω 1 0, On the other hand, if 0 = 0 , the equilibrium point is therefore Ω 3 y e + mg K , 0 . The Jacobian matrix of the corresponding system is By considering the linearization around the Ω y 0 , 0 , the eigenvalues are solutions of the following characteristic equation, namely For the equilibrium points Ω 1 and Ω 2 , the following inequalities are verified 1 ˃0, 2 ˃0 and 1 2 ˃ 3 . According to the Routh-Hurwitz criteria, hence, these points are all unstable. And the other hand, considering Ω 3 , we have 1 ˃0, 2 ˃0 and 1 2 ˂ 3 . According to the Routh-Hurwitz criteria, the corresponding equilibrium point is stable.
Note: According to the stiffness criterion, it's easy to verify that the stability of the Maglev system can be expressed as follows: a) Stable: K > 0 b) Unstable: K < 0

Case u(t) = 0
Many authors have proven that numerous dynamical systems encountered in literature can depict a variable fixed point called AC equilibrium points [54] (the reader can also find an excellent text in Ref. [55]). This will allow a supplementary stringent aspect of many mechanicals systems which is submitted to small vibrations at a rest point.
Taking from (6) We have two equilibrium points as follows We obtain three equilibrium points With These results are relaxed and proven at the end (see the appendix).

Amplitude response
In the aims to observe the evolution of the amplitude response of the Maglev system, an analytical expression of displacement z(t) and the current i(t) flowing through the electromagnet are for high interest. In order to solve the set of Eqs. 6 and 7, let us make the change of variables: � y e + mg K Setting (t) weak, then the product (t) (t) is too weak. Letting u(t) = U m sin(Ωt) , it's obvious to obtain from (22) that If we seek for harmonic solution, we set C = 0 and derive the following solutions and substituting (24) Finally, the analytical solution of the electromechanical system (6)-(7) is as follows: where Comparing (20) and (29), we identify the constant parameter defined as Taking into consideration Eq. (29), the solution can be defined as z 0 is the initial position of the maglev which stands as a mean value of z(t) . z m is the amplitude of z(t) . I 1 and I 2 are the amplitude of the current for the fundamental and for the third harmonic. z , 1 and 2 denote the corresponding phase for the position z(t) and the current i(t) respectively. It's straightforward to obtain after some developments the following results: where (a 1 , a 2 ) and (a 01 , a 02 ) denotes the constants coefficients of the trigonometric function (cos(2Ωt), sin(2Ωt)) and (sin(2Ωt), cos(2Ωt)) referring to the solutions y(t) and (t) respectively.
From Eq. (33), we can highlight two essential observations: (a) the presence of harmonic currents is one of factors which have the most influence on the technical design of the electromagnet or the associated generator, more particularly with reference to the waveform. These currents cause distortions that move the wave profile away from the ideal sine reference. Accordingly, in this article, the system's parameters should properly choose to eliminate the third harmonics or to render this higher order component negligible and thus attenuate their effect. (b) The mechanical frequency vibration is probably twice that of the source voltage.

Simulations validation of the theory
Let consider the physical parameters of Maglev system as presented in Refs. [21,36,38] which are listed in Table 1.
For the sake of clarity, we choose the control input voltage as u(t) = 220 √ 2 sin(2 ft) and the radian frequency is = 2 f where f is chosen as 50Hz . Now, let us analyze the dynamical behavior of the displacement and the current flowing in the coil. The numerical method used to find the desired solutions is fixed step size fourthorder Runge-Kutta algorithm with step 67 × 10 -5 . The initial conditions considered are z(0) = z 0 ;̇z(0) = 0 and i(0) = z(0) (0) = 0 . We derive the corresponding  showed on Fig. 2a, b. We clearly observe that the oscillation frequency is twice the one of the current. In order to investigate the effects of each parameter on the behavior of the system, let us depict Fig. 3. The red curve represents the numerical solution obtained from Eq. (19) and the blue one is the analytical solution. We notice that for weak values of the vibration frequency, the two curves are similar. A disparity appears in the gap [1.5; 1.85] and reveals an irregular behavior of the electromagnet in the dynamics of the overall system. The mathematical approximation of the Eq. (19) does not take into consideration this discontinuity. This behavior is observed for m = 520Kg. Figure 3 shows the same curves for m = 620Kg. We obviously point out that the solution obtained remains rigorously true in a sufficient time interval to be able to make a complete analysis of the Maglev. From Fig. 2a, it appears that the great variation of the displacement may destroy the performance of the maglev system. The current that is given to the electric coils of the guideway will be alternating in nature. Thus the polarity of the coils will be changing in period. Figures 3 and 4  respectively, the numerical solution of the displacement and the analytical solution from Eq. (19). We observe a perfect agreement between these curves which confirms the performance of the general proposed solution within a suitable interval. In addition, one of the important aspects of this study is also to know the parameter which would considerably reduce the harmonics in the signals of the Maglev, in fact, without loss of generality, from Eqs. (27) and (31), it clearly appears that the relevant design of the value of air spring (K) of the electromagnet assigns considerably the dynamical properties of the vehicle and then becomes fundamental to eliminate the harmonic in the Maglev system. Hence, the change causes a pull force for the train in the front and when it is added to this force, the magnetic field behind the train adds more forward thrust. Figures 5a, b, present the numerical and analytic solution of the displacement. These results are obtained for the small values of the vibrating frequencies (10-30 Hz), when the frequency increases, the maglev system become boring and the control problem takes place. One important issue to improve technologically the Maglev system is the selection of a device which enables slow variation of the displacement of the air gap, Fig. 6 shows the zone in which the oscillations appears.

Control scheme description and objective
To deal with the solution of the Maglev control, the electromechanical Eqs. (6)-(7) can be put in the canonical form of differential equation. By letting y(t) = x 1 (t),ẋ 1 (t) = x 2 (t), x 3 (t) = , g 0 = K m x 1e + g and x 1e = y e , we obtain Let consider x 1s , x 2s and x 3s be the target values of x 1 , x 2 and x 3 respectively. It obvious to note that from Eqs. (39), the equilibrium point is x 0 = x 10 , 0, x 30 T , where x 30 satis- The control objective is to design a robust control u(t) which will drive the states x 1 , x 2 and x 3 to its desired constant values x 1s , x 2s and x 3s respectively which can be expressed in the mathematical form by

Problem formulation
We first introduce the change coordinators as follows: It's obvious to remark that when 1 , 2 , 3 are driven to zero as t → ∞ , then x 1 will converge to x 10 , x 2 will converge to zero, and x 3 will converge to x 30 = √ −2 K ℏ x 10 + g 0 ℏ . The first derivative of Eq. (41) leads to the following (39)

Robust sliding mode control design
Let consider Eqs. (42) of uncertain nonlinear maglev system as follows where Λ(t) is the bounded disturbance that affects the dynamic of the maglev system. We assume that

Remark 2.
The diagram overview of the control is illustrated in Fig. 7. The full order adaptive nonsingular sliding mode control is considered in the scheme under the exogenous disturbances.
A sliding surface is chosen as 1 , 2 and 3 are constants which are selected such that the polynomial p 3 + 3 p 2 + 2 p + 1 is Hurwitz. The aim of this part is to develop an adaptive nonsingular sliding mode controller such that the closed loop (44) states converge to zero in limited time.
We design the control law as u eq can be designed as u n can be taken as where v is designed as where u n (0) = 0, ′ is a small positive number and the positive constant which is automatically update by the designer according to the following law: where 0 ≻ 0 is adaptation gain Theorem The system (44) reaches and converges to the sliding surface in limited time if the selected sliding surface (45) and the control law as above are considered.

Remark 3
In order to obtain the expression of u n , let consider the Laplace transform of Eqs. (48) We also assume that u n (0) = 0 and for t = 0 , we have

Numerical simulation investigating robust control of Maglev train
The simulations are based on Fig. 7. The values on the Tab.1 are considered in the numerical investigations. Matlab simulations are performed to verify the performance of the proposed control scheme. The external disturbance representing time-varying parameter uncertainties is considered as Λ(t) = 5.10 −3 sin(0.5t) . The initial conditions of the differential maglev Eqs. (44) are 10 −3 , 0.8, 0.4, 0.6, 10 −3 . The results of the proposed scheme are validated using Routh Hurwitz criteria. We first select the parameters of the polynomial as follows = 35, = 10 −3 , 1 = 100, 2 = 121, 3 = 151, 0 = −0.01 . It's obvious to mention for this case that the Hurwitz criteria is not satisfied. With this condition the maglev train remains uncontrollable. Figures 8 a, b, c present the surface and the states 1 (t) and 2 (t) when the control is not effective.
For the sake of simplicity, we choose i = 1, (i = 1, 2, 3) which confers to the stability of the maglev system (Hurwitz criteria is verified). For convenience, from Fig. 6, we select arbitrarily x 10 close to y e in order to remain in the oscillatory region such that Let consider the initial position as y e = 10mm , the target air gap is taken as y = 3.7mm . The simulation results can be observed in Fig. 9 a, b, c. From Fig. 9a, we can observe that the displacement of the rail and the air gap reaches stability. From Fig. 9 b, c, one can observe that, after a transient period, the current and the speed become stable and converge to zero which is in good agreement with our theoretical results. It is important to recall that the convergence to zero of the state variables describing the current through the coil doesn't imply that the quantity vanishes, but it also shows the great reduction of the harmonic in the maglev system. Figures 10 and 11 depict the sliding surface s(t) and the control gain k(t) respectively. It is observed that the control effort which is measured by ∫ ∞ 0 u(t)dt is weak and therefore can be easily implemented in real time. From the above results, the response system of the Maglev system has a good control performance. We can clearly observe that the proposed sliding mode control helps to stabilize the maglev on the rail at around 9.5 s despite perturbations. Let's point out some promising results obtained in the literature using the relevant control (60) | | y e − x 10 | | ≤ 10%mm.
techniques for the Maglev [56,57]. Although, findings are satisfactory, they remain very conservative and present a real problem of energy consumption. Furthermore, the results obtained in this research are in perfect agreement with the theoretical results and provide a significant contribution to the Maglev control domain.

Conclusion
This work explores a new and simple concept of the critical analysis for the motion stability of a general class of Maglev train incorporated into differential equations. This paper discussed the AC equilibrium point and its stability. It is found that the initial value of the displacement varies between unstable and stable intervals with the amplitude of the control input. To investigate numerically the complex Maglev dynamic, the derivation of the equilibrium point and the corresponding stability analysis is performed and the behavior depending on the control input have revealed the behavior of the system's states with respect to the vibrating frequencies. The analytical studies and numerical ones show perfect agreement and can allow predicting efficiently the displacement's amplitude of the electromagnet as well as the current flowing in the coil. These features are of great interest in the design of maglev vehicles. The adaptive robust nonlinear sliding mode control for the maglev system is proposed in the second part of this work. The bound of uncertainty is assumed to be known. To deal with the difficulties face in Refs. [38], the robust nonsingular sliding mode controller is designed in the frequency domain which guarantees the uniform boundedness of the system regardless of the disturbances. From this study, we can point out a major drawback such as an accurate analytical expression of the displacement and the current due to the high nonlinearity found in the Maglev system. Our next investigation will be on the test of the present model when the edge effect and the magnetic flux leakage are not neglected, consequently, the flux is not uniformly distributed in the air gap. We will propose to use the fuzzy logic issue to reduce the complexities in analysis and control design. We think also that the analysis of the proposed system with fractional damping is another good idea.

Compliance with ethical standards
Conflict of interest The authors declare that there are no conflicts of interest, regarding this paper.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.