# Linearization Method of Nonlinear Magnetic Levitation System

Authors: Dini Wang; Fanwei Meng; Shengya Meng
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

Linearized model of the system is often used in control design. It is generally believed that we can obtain the linearized model as long as the Taylor expansion method is used for the nonlinear model. This paper points out that the Taylor expansion method is only applicable to the linearization of the original nonlinear function. If the Taylor expansion is used for the derived nonlinear equation, wrong results are often obtained. Taking the linearization model of the maglev system as an example, it is shown that the linearization should be carried out with the process of equation derivation. The model is verified by nonlinear system simulation in Simulink. The method in this paper is helpful to write the linearized equation of the control system correctly.


---

## Body

## 1. Introduction

Magnetic levitation generates eddy-current effect by metal objects under the action of high-frequency electromagnetic fields, which enables metal objects to be suspended stably in the air [1]. Magnetic levitation technology is developing rapidly because it essentially eliminates friction and has advantages unmatched by other traditional technologies [2]. The application of the magnetic levitation technology covers numerous fields, such as high-speed magnetic bearing [3, 4], high-speed levitation train [5–7], and wind tunnel magnetic levitation [8]. Therefore, the control design for the maglev system is particularly important.The maglev system is a nonlinear system, which causes great difficulties in its control design [9–11]. Literature [12] proposes a model-independent control system design method based on U-model. In the linear control design method, linearization is often the first step in the control system design and modeling of the magnetic levitation system [13].It is generally considered that the nonlinear equation is linearized as long as the Taylor expansion is used [14]. However, the Taylor expansion method only applies to the original nonlinear function. The equations for general magnetic levitation systems are usually obtained from many derivations. In the process of building the nonlinear model, the calculation of partial derivative has been done. According to the derivation algorithm, if we continue to take the derivative, then the partial derivative should be treated as a constant. In traditional modeling, this constant has been treated as a variable. And the error is introduced as we continue to differentiate. If linearization is simply carried out through the Taylor expansion method according to the derived equation, the result may be wrong.Based on the analysis of the Lagrange equation [15], this paper points out that the partial derivative of Lagrange equation should be taken according to the operating point to get the correct linearized equation. Simulink simulation results show that the method proposed in this paper is very effective for the linearization of the magnetic levitation system and is worth popularizing.


## 2. Maglev System Model and Its Linearization

In general, the linearization of nonlinear systems is through Taylor series expansion, and a linear relation can be obtained after omitting the higher order terms [16]. Taylor series refers to the expansion of a nonlinear function fx,y. For example, formula (1) is the equation of a prime mover [17].(1)Jdωdt=TMω,u−TDω,where TM, the driving torque, is a nonlinear function of the rotation rate ω and the control input u and TD, the load torque, is a function of the rotation rate ω. If TM and TD are expanded as Taylor series and take the first term, the linearized equation of formula (1) can be obtained. This is the general notion of linearization. But for a general nonlinear system, it takes a lot of arithmetical operation to get from the original equation to the nonlinear equation. Using Taylor series to linearize the nonlinear terms in the final result, wrong result is often obtained.The model of the typical electromagnetic suspension (EMS) train control system is shown in Figure1.Figure 1
Model of an electromagnetic suspension system.The motion equation of the maglev system model can be written as [11, 15](2)md2ztdt2=−μ0N2am4itzt2+mg,(3)ditdt=itztdztdt−2μ0N2amztRmit−ut,where μ0 represents the permeability in a vacuum; m represents the mass of the suspension magnet, m=15kg; am represents the effective area of the magnetic pole, am=1.024×10−2m2; N is the number of turns of the coil on the electromagnet, N=280; Rm represents the resistance of the coil, Rm=1.1Ω; ut represents the control voltage at both ends of the coil; it represents the current in the coil; zt is the gap between the electromagnet and the guideway; and g represents the gravitational acceleration. Suppose that the state vector is(4)xt=ztz˙titT=x1x2x3T∈R3.Consider the nonlinear termsit/zt2,it/zt,ztit in (2) and (3) as nonlinear functions. The linearized equation of the system can be obtained by Taylor series expansion [16].(5)x˙=Ax+Bu.The corresponding matrices can be written as(6)A=010μ0amN22mi02z030−μ0amN22mi0z02−2Rmμ0N2am−z˙0z02i0i0z0−2Rmμ0N2amz0+z˙0z0,B=002μ0N2amz0.Set the nominal working point asz0, z0=4.0×10−3m, and the corresponding operating current i0=3.0538A. Then, the linearized equation can be written as(7)x˙=0104900.10−6.4184−6659.4763.45−8.7228x+007.9298u.Formulas (5) and (7) are directly obtained by the Taylor expansion method according to formulas (2) and (3). This is a common linearization method, but it proves (see below) that the resulting linearization equation is incorrect.Next, we start with the establishment of the equation. For the dynamic system withk generalized coordinates q1,…,qk, the Lagrange equation is [2, 4](8)ddt∂T∂q˙i−∂T∂qi=−∂V∂qi−∂R∂q˙i+fit,i=1,2…,k,where T represents the kinetic energy stored by the dynamic system, V represents the potential energy of the system, fit represents the external force varying with time, and R represents Rayleigh’s dissipation function, and its general form can be written as(9)R=12Rmq˙2.Suppose the generalized coordinates of the maglev control system are air gapz and charge q (note q˙=i), then the kinetic energy and potential energy of the system are, respectively,(10)T=12mz˙2+12Lq˙2,(11)V=−mgz,where(12)L=μ0N2am2z.By substituting formulas (10)∼(12) into Lagrange equation (8), the equation of generalized coordinates z and q can be written as(13)ddtmz˙−12q˙2−μ0N2am2z2=mg,(14)ddtμ0N2am2zq˙=−Rmq˙+u.The derivative of (14) on the left is(15)ddtμ0N2am2zq˙=μ0N2am2zditdt−μ0N2ami2z2dztdt.Thus, according to (13) and (14), the equation of motion of the system can be finally obtained as(16)d2ztdt2=−μ0N2am4mitzt2+g,(17)ditdt=itztdztdt−2μ0N2amztRmit−ut.Equations (16) and (17) are equivalent to (2) and (3) when analyzing the maglev system. Since we have taken the partial derivative once when we derived (17), the partial derivative should be a constant to the operating point in linearization, that is, it corresponds to a constant coefficient instead of a variable. Therefore, if linearization is carried out, the corresponding constants of all partial derivatives in (15) should be taken, and (15) can be rewritten as(18)ddtμ0N2am2zq˙=μ0N2am2z0ditdt−μ0N2ami02z02dztdt.The corresponding linearized equation can be obtained by substituting (18) into (14). it/zt2 in (16) is indeed a nonlinear function, and its linearization relation can be obtained by Taylor expansion.(19)μ0N2am4mitzt2=μ0N2am4mi0z02+μ0N2am2mi0z02Δi−μ0N2am2mi02z03Δz.In this way, by substituting (18) and (19) into (14) and (16), the final linearized equation can be obtained as(20)x˙=010μ0amN22mi02z030−μ0amN22mi0z020i0z0−2Rmμ0N2amz0x+002μ0N2amz0u,i.e.,(21)x˙=0104900.10−6.41840763.45−8.7228x+007.9298u.Obviously, (20) and (21) are obviously inconsistent with (5) and (7). The state matrix of (7) has the element −6659.4 in the third row and first column, while the element at this position in (21) is 0. In order to verify this linearization method, we use the nonlinear system simulation in this paper because when the deviation is very small, the nonlinear system results should be consistent with the linearization result.


## 3. Simulation of Nonlinear System

Next, we will simulate and analyze the nonlinear equations of the maglev system. In the nonlinear simulation of Simulink, a custom function module Fcn is mainly used. The function of this module is to process the input by the rules conforming to C language. The corresponding formulas may consist of one or more components.u is the input of the module, which supports vectorized input. If u is a vector, u [i] represents the ith element of the vector, and u [1] or a single u represents the first element. The function module can perform both mathematical operations (such as addition, subtraction, multiplication, and division) and relational operations (such as identically equal, not equal, greater than, and equal), as well as three kinds of logical operations. The precedence of these operators conforms to the rules in C language.The simulation model of the magnetic levitation system based on differential (2) and (3) is shown in Figure 2. In Figure 2, the function of Mux module is to organize several signals into one signal. Here, variables are expressed in increments so as to compare with linear results. In the figure, f1u, f2u, and f3u are all custom function modules.Figure 2
Simulation model of the nonlinear system.Adopt the state feedbackk=k1k2k3=24843.50249.7285−20.9242 in the simulation. And the feedback control signal should be elicited from Δz,Δz˙,Δi. In the case of module f3u, it describes the operation on the right side of formula (2).(22)f3u=−μ0N2am4m△i+i0△z+z02+g.The input is two variables:Δz,Δi. The setting of f3u is shown in Figure 3.Figure 3
The dialog box of Fcn3.In Figure3, 1.6814∗10−5 corresponds to the value of μ0N2am/4m, 9.8 is the value of gravitational acceleration g, and u1, u2 are input signals arranged in a certain order by Mux module in front of the module f3u. Here, u1 corresponds to Δi and u2 corresponds to Δz, and i0 and z0 are values of corresponding nominal operating points. The values of i0 and z0 at the equilibrium point are calculated in advance according to the following equation:(23)mz¨t=−μ0N2am4i0z02+mg=0.The functionf1u describes the first term on the right side of (3), and f2u describes the second term on the right side of (3). The problem of the nominal operating point is also set as follows:(24)f1u=△i+i0△z+z0△z˙,f2u=2μ0N2am△z+z0RmΔi−Δu.


## 4. Verification of Linearized Equations

Two linearization methods are proposed above. One is the simple Taylor series method ((5) and (7)), and the other is the method of taking the partial derivative as a constant when taking the column equation ((20) and (21)). In order to further verify the linearization method in this paper, the simulation checking is adopted, and the simulation results of nonlinear systems with small deviations are compared.Figure4 shows the simulation block diagram of the linearized system. The dotted box in the figure shows the difference between the two linearized results (7) and (21). The same state feedback is used in the simulation.(25)k=k1k2k3=24843.50249.7285−20.9242.Figure 4
Simulation model of the linearized system.Figure5 shows the simulation results of the step response. The short line in the figure is the simulation result of system (7), and the dotted line is the simulation result of system (21) of this method.Figure 5
Comparison of the simulation results.The simulation block diagram of the nonlinear system is shown in Figure2. The solid line in Figure 5 is the simulation result of the nonlinear system, and the amplitude of the input step signal is 0.1mm, 0.5mm, and 0.1mm, respectively. When a small signal (0.1mm) is input, the dotted line of the system (21) is basically coincident with the solid line of the nonlinear system. The simulation result of system (7) is obviously inconsistent with the nonlinear result under small signals, especially the static differences are not equal. In this case, the static error of the nonlinear system is relatively easy to calculate. The static error between the system output and the step input is also shown in Figure 5 (the calculation process of the static error is abbreviated). The simulation results show that the static error of the system (21) is consistent with that of the nonlinear system.


## 5. Conclusions

The linearization of nonlinear magnetic levitation cannot be carried out simply by Taylor expansion according to the derived nonlinear equation. In the process of derivation, the partial derivative is included in the equation, and the linearized equation can be obtained by taking the partial derivative at the operating point. The method is simple and the physical concept is clear. On the contrary, if we linearize the equation by the Taylor series based on the derived nonlinear equation, wrong result may be obtained. In this paper, the simulation method of the nonlinear magnetic levitation system is presented, which separates the operating point (setting) from the incremental change, and is worthy of promoting.



---

*Source: 9873651-2020-06-22.xml*