Open access peer-reviewed chapter

Modeling, System Identification, and Control of Electromagnetic Actuators

By Alexandru Forrai

Submitted: October 30th 2017Reviewed: February 8th 2018Published: July 11th 2018

DOI: 10.5772/intechopen.75088

Downloaded: 308

Abstract

This chapter is dedicated to modeling, system identification, and control of electromagnetic actuators with the main focus on the actuators used in magnetic levitation, in fuel injection systems, and in variable valve timing (VVT). These actuators have a simple structure, good reliability, and low manufacturing costs. However, from control viewpoint, they are nonlinear systems and are open-loop unstable. Therefore, mathematical modeling, system identification-based parameter estimation, and control strategies are presented, when the moving armature is controlled around an equilibrium position or is controlled between the two extreme positions of the armature.

Keywords

  • electromagnetic actuator
  • modeling
  • identification
  • gain scheduled control

1. Introduction

Electromagnetic actuators are widely used in the industry, and they transform the electric energy into linear motion. From the large variety of applications, in this chapter we are going to focus on:

  • Magnetic levitation

  • Fuel injection systems and variable valve timing actuators used in internal combustion engines

These applications are relevant from control point of view: in the first case, the moving armature is controlled around an equilibrium position; in the second case, the armature might go under control between the two extreme positions—armature open and armature close.

Nevertheless, magnetic levitation—in particular a magnetically levitated train—is a good example, where closed-loop control plays a key role, since the open-loop system is unstable [1]. The system can be linearized around an operating point, and a linear controller can be designed.

Furthermore, magnetic bearings and their control are from a long time the focus of control system design community. Feedback linearization and asymptotically exact linearization of an active magnetic bearing are presented in [2, 3]. Advanced control strategies are discussed in detail in [4, 5].

Therefore, it makes sense to develop high-accuracy mathematical models, to investigate methods for parameter identification, and finally to apply control strategies to improve performance and reliability of the system.

We will discuss these topics in the next sections, but before that let us focus on applications, where electromagnetic actuators are widely used.

1.1. Magnetic levitation

Magnetic bearings in combination with high-speed electric motors are used across many industries, from oil and gas industry to electric power generation industry (i.e., high-speed electric generators) and from the semiconductor industry to nuclear industry, etc.

The main structure of the magnetic bearing is shown in Figure 1 (reproduced from [6]).

Figure 1.

Applications of magnetic levitation.

Another well-known application of the magnetic levitation is the magnetically levitated high-speed train (Maglev; see Figure 1), having speeds over 500 [km/h] [7, 8].

1.2. Fuel injection and variable valve timing (VVT)

The main purpose of the fuel injection system is to deliver fuel to the cylinders. However, how that fuel is delivered is that it makes the difference in engine performance, emissions, and noise characteristics.

Most notable advances achieved in diesel engines resulted directly from superior fuel injection system designs [9].

Unlike its spark-ignited engine counterpart, the diesel fuel injection system delivers fuel under extremely high injection pressures (e.g., around 2000 [bar]). This means that the system component designs and materials should be selected to withstand higher stresses [10].

The actuators used in diesel fuel injection systems can be either electromagnetic (our focus) or piezoelectric [11]. A diesel fuel injection system using an electromagnetic actuator, from Bosch [12], is shown in Figure 2.

Figure 2.

Diesel fuel injection system from Bosch.

Nowadays, most of the fuel injection systems are electronically controlled. However, it is still not enough to deliver an accurate amount of fuel at the proper time to achieve good combustion. Additional aspects are critical to ensure proper fuel injection system performance, such as [9]:

  • Fuel atomization—ensuring that fuel atomizes into very small fuel particles is a primary design objective for diesel fuel injection systems.

  • Bulk mixing—while fuel atomization and complete evaporation of fuel are critical, ensuring that the evaporated fuel has sufficient oxygen during combustion is equally important to ensure optimum engine performance.

  • Air utilization—effective utilization of the air in the combustion chamber is closely tied to bulk mixing and can be accomplished by dividing the total injected fuel into a number of jets.

While conventional fuel injection systems employ a single injection event for every engine cycle, newer systems can use multiple injection events [12].

Using multiple injections—during every engine cycle—higher engine performance and lower engine noise can be achieved. However, the injector lifetime might be reduced, and therefore advanced control algorithms as well as malfunction detection and fault isolation algorithms can be applied (see next sections).

Another relevant application is the electromechanical valve actuators used in automotive engines, to achieve variable valve timing (VVT). With VVT, larger valve overlap, valve lift, duration, and timing adjustments can be achieved depending on engine speed, load, and temperature.

Variable valve timing leads to improved fuel economy and lower emissions by decoupling the valve timing from the piston motion [13]. This is especially valid in case of advanced combustion technologies, as described in [14, 15].

However, the moving components of the valve actuators create unnecessary wear and excessive noise. The armature landing speed shall be kept, e.g., under 0.1 [m/s]; otherwise, they are excessively loud and are damaging to the actuator and engine valve.

Whenever high-performance and high-accuracy control is required, the electromagnetic actuator is driven by a half H-bridge (see Figure 3), which might be equipped optionally with a current sensing resistor RSENSE.

Figure 3.

Electromagnetic actuator driven by a half H-bridge.

Typical voltage and current waveforms as well as the switching order of the commutation elements T1and T2are shown in Figure 3.

After the electromagnetic armature is pulled up, the actuator current is reduced, and the commutation elements are controlled via pulse-width modulation (PWM).

During armature movement, due to the induced electromotive force (e.m.f.), a small current dip as well as a small current peak might be observed (see Figure 3).

The duty factor of the actuator—specified on the data sheet—is defined as

Duty%=TONTON+TOFF100%E1

In practice, exceeding this value might shorten significantly the lifetime of the actuator.

2. Mathematical modeling

The mathematical model of the electromagnetic actuator is described by the voltage equation and by the motion equation.

The voltage equation is

vin=Ri+∂Ψididt+∂ΨzdzdtE2

where vinis the applied voltage, iis the armature current, Ψis the armature flux, zis the armature position, and Ris the electrical resistance of the coil.

If we note with v=ż, the armature speed and then the equation of the motion can be written as

mz¨=FSFmE3

where mis the moving mass, Fmis the electromagnetic force, and FSis the spring force. Since the armature displacement often is very short, the spring force can be considered constant. In case of magnetic levitation, the spring force is replaced by the weight of the moving mass.

The electromagnetic force can be expressed based on the electromagnetic co-energy Wco:

Wco=0iΨdiE4
Fm=Wcozi=ctE5

Then, we can also write

Fm=0i∂ΨzdiE6

The Ψ=Ψizand Fm=Fmizstatic characteristics can be measured. The flux-linkage characteristic is derived by integration (very often the current decay test is used). Thus, the flux linkage for one fixed position zis calculated by

Ψiz=0vintRitdtE7

where at t=0and the following conditions hold: vin=0, i0, and di/dt=0.

Although the model does not take into account the effect of eddy currents, the numerical model can be very accurate and might be written formally into a nonlinear form:

ẋ=fx+i=1mgxuE8

where x=ivzTrepresents the state of the nonlinear system, u=vinFSTis the input vector, and fxand gxare nonlinear functions of the state x. The output vector yis

y=hxE9

where hxin the most general case is a nonlinear function.

Finally, in the aim to illustrate our investigations, let us consider an electromagnetic actuator with parameters (catalog data) mentioned in Table 1 [16].

TypeSolenoid valve
Stroke length10 [mm]
Operating voltage24 [V] d.c.
Maximum current0.6 [A]
Resistance40 [Ω]
Inductance0.35–1.1 [H]
Number of turns2240

Table 1.

The solenoid parameters.

2.1. Nonlinear model and piecewise linearization

The mathematical model described above is too general and is difficult to handle in analytical form. Therefore, we define an analytical model set, which describes the flux-linkage characteristic as

Ψiz=Ψmax1expic1+c2zE10

where the parameters of the model set are Ψmax, c1, and c2.

Furthermore, the partial derivatives of the flux-linkage are

∂Ψi=Ψmaxc1+c2zexpic1+c2zE11
∂Ψz=Ψmaxc2ic1+c2z2expic1+c2zE12

The approach presented in this section is reproduced from [17].

If we approximate the exponential term by Taylor series, we have

expic1+c2z1ic1+c2z+i22c1+c2z2E13

Thus, the magnetic force—based on the analytical model—can be expressed as

FmΨmaxc2i2c1+c2z212i3c1+c2z+i28c1+c2z2E14

The model set above is validated against the measured static (flux and force) characteristics. The “dots” in Figure 4 represent the measured data, and the solid lines represent the calculated model using the above model set, with Ψmax=0.45Wb, c1=0.4A, and c2=0.375103A/m. The above parameters are derived using nonlinear least squares, fitting the measured data (obtained using the current decay test and force measurements) with the analytical model.

Figure 4.

Flux and force characteristics.

Next, let us introduce the following notations, which help us to rewrite the model in a convenient form: χi=∂Ψ/i, χz=∂Ψ/z, and χf=Fm/i.

Since the magnetic force Fmdepends on the square of the current i2, when the current is zero i=0and then χf=0, there is no division by zero in the model.

Thus, the voltage and motion equations are written as

didt=Rχiiχzχiv+1χivinE15
dzdt=vE16
dvdt=χfmikmz+1mFS0E17

Finally, using a piecewise approximation, the system can be written in state-space form as

ẋ=Aizx+Bizuy=CxE18

where x=izvTis the state-space vector, u=vinFS0Tis the input vector, yis the output vector, the Aizand Bizare current and position dependent matrices and C=110if the armature current and position are sensed. We remark, that in practice sensing the armature speed and/or positions with sensor(s) might be expensive solution. Therefore, often only the armature current can be sensed in a cost-effective manner.

The terms Aizcan be written as

Aiz=R/χi0χz/χi001χf/mk/m0E19

where

χzχi=c2ic1+c2zE20

The term Bizcan be written as

Biz=1/χi00001/mE21

where

1χi2c1+c2z2+2ic1+c2z+i22Ψmaxc1+c2zE22

or a coarser approximation will be

1χic1+c2z+iΨmaxE23

Last but not least, the armature movement is subject to the following constraints:

vt=0ifzzmaxandFsFm00ifz0andFsFm0E24

as well as zminztzmax, where zminand zmaxare the minimum and maximum displacements of the armature.

2.2. Linearized mathematical model

From control engineering viewpoint—in case of some applications (e.g., magnetic levitation)—the piecewise linearized model might be too sophisticated. Therefore, in this section a linearized mathematical model around an operating point is derived [1, 18].

Let us approximate the magnetic force as

Fmγi2c1+c2z2E25

where γis a constant.

The equation of motion can be written as

Mzz¨i=mz¨Fs+γi2c1+c2z2=0E26

The equation above can be linearized around an operating point p0=z0z¨0i0as follows:

Mzz¨i=Mz0z¨0i0+Mzp0zz0+Mz¨p0z¨z¨0+Mip0ii0E27

which can be further written as

Mzz¨i=2γi02c1+c2z03zz0+mz¨z¨0+2γi0c1+c2z02ii0=0E28

If we denote with Δz=zz0and Δi=ii0, we obtain

2γi02c1+c2z03Δz+mΔz¨+2γi0c1+c2z02Δi=0E29

If we divide the equation with the moving mass mand apply the Laplace transform, we obtain

s2a2Δzs+kΔis=0E30
ΔzsΔis=ks2a2=ksas+aE31

where kand aare varying with the equilibrium point i0z0:

k=2γi0mc1+c2z02E32

a2=2γi02mc1+c2z03E33

It means that a family of transfer functions are obtained and the system can be viewed as a linear parameter-varying (LPV) system.

The variation of kand avalues with the equilibrium position z0is shown in Figure 5 and can be well approximated by quadratic functions.

Figure 5.

Plant gain and pole variation with the equilibrium position.

3. System identification

3.1. Clustering-based system identification

In the previous section, we have seen that using the current decay test and the nonlinear least squares method, the parameters of the mathematical model can be identified.

However, the current decay test is time-consuming, since measurements shall be performed for each grid point defined by armature current and position iz. Thus, the obvious question might arise: is there a faster solution to identify the parameters?

During the system identification process, we will note the system’s input and output at time tby utand yt, respectively [19].

For single-input single-output linear systems, we can write

yt=φTtθE34

where θis the parameter vector (unknown) and the φis the recorded (known) input-output data vector:

θ=[a1anb0bm]Tφt=[yt1ytnututm]TE35

To emphasize that the calculation of the ytis from the past data, we will write

ŷt=φTtθE36

Now, suppose for a given system that we do not know the values of the parameters in θ, but we have recorded inputs and outputs over the time interval. If the input signal is persistently exciting—condition described in details in [19, 20]—then the solution can easily be computed by modern software tools.

In this section a clustering-based identification method, proposed by Ferrari-Trecate et al. (2003) is used (see [21, 22]), where the plant is assumed to be described by piecewise linear models having ssub-models, such as

yt=φTtθ1+wt,ifφTtC1φTtθs+wt,ifφTtC1E37

where wtis white noise, θii=1,,sare the parameter vectors, φtis a regression vector, and nis the order of the piecewise ARX (PWARX) model.

It is assumed that the order of each sub-model is the same, and u(t) and y(t) are the input and output, respectively.

Furthermore, it is assumed that Ci=1sare polytopic and they satisfy the well-posed condition: i=1sCi=C, CiCj=, and ij.

An important phase of the system identification experiment is input signal design. In case of nonlinear systems, a multilevel random signal is often used [23, 24], and a bi-level pseudorandom binary signal (PRBS) is not suitable for nonlinear systems.

The generation of the multilevel random signal—using shift registers—is done according to [25]. Figure 6 shows a five-level random signal with maximal length, using four-shift registers with coefficients a1=1,a2=1,a3=1,a4=2[25].

Figure 6.

Multilevel random input and corresponding output signal.

This input signal is applied—when the armature is fixed—in order to identify the dynamic inductance denoted by χiand the electrical resistance R. The output signal—armature current—when the multilevel random signal is applied is shown in Figure 6.

Next, the voltage equation is written in discrete form at the time moment t=tk:

ik=TSχi+RTSvink+χiχi+RTSik1E38

where TSis the sampling time.

The equation above defines the regression space (see also references [21, 22]) having in this case two axis, defined by ik1and vink. The data collected during the system identification experiment is shown in the regression space in Figure 7.

Figure 7.

Regression space—Clustering-based identification.

The regression space is clustered in five different regions for i=0.1,…,0.5A, and the parameters are identified for each case using the least squares method. In each defined cluster, we assume that χiis constant, and basically we use a piecewise linear approximation of χi.

The system identification experiments are repeated around different positions, when the armature is fixed; thus, the function χî=χîizcan be estimated.

Now, using the nonlinear least squares, we can minimize the Jobjective function:

J=minχiizχî(iz)2E39

and we can find out the estimated parameters of the model Ψ̂max=0.437Wb, c1̂=0.37A, and c2̂=0.36103A/m. Values, which are in good accordance with the values, are found via the current decay test.

Having the χi=χiizfunction identified, the parameters of the model set, namely, Ψmax, c1, and c2, are found.

This identification is repeated only around different positions z, when the armature is fixed and thus is much faster than identifying Ψ=Ψizaround different current and position values using the current decay test.

3.2. Identification under closed-loop

In practice, it might be the case that the system is open-loop unstable; thus, system identification experiments have to be performed under closed-loop (for more details see [26, 27]).

Closed-loop identification is a very challenging task. Due to the presence of feedback loop, the input signal might not be persistently exciting. In the aim to achieve a persistent excitation of the system, it is recommended in [19] to switch between different simple controller structures.

First, the system shall be stabilized under feedback around an equilibrium position, as shown in Figure 8; details about the controller design Ksare described in the next section.

Figure 8.

Closed-loop system identification.

Under closed loop, the reference input (armature position) is disturbed by a persistently exciting input signal (i.e., pseudorandom binary signal)—as shown in Figure 9, and three linear transfer functions are identified, which are defined as

Ts=YsRs=PsKs1+PsKsE40
Hs=UsEs=Ks1+PsKsE41
Ps=YsUs=TsHs1TsE42

where we used the well-known identity Ss+Ts=1:

Ss=EsRs=1TsE43

Figure 9.

Input and output signal used to identify T s .

The procedure can be repeated around different equilibrium positions; thus, a family of transfer functions can be obtained.

4. Control of electromagnetic actuators

Let us start with the easier case: the moving armature is controlled around an equilibrium position—magnetic bearings and magnetically levitated high-speed trains are typical applications.

The linearized mathematical model, around an equilibrium position, can be written as

Ps=ΔzsΔis=ksas+aE44

where kand aare strictly positive values, varying with the equilibrium point i0z0.

In this section, we are looking for a linear controller, which can stabilize the plant and can fulfill performance and robustness requirements [18].

4.1. PD controller

A very simple PD controller, which can stabilize the plant, is

Ks=kDs+aE45

The closed-loop transfer function shows that

Ts=kkDs+kkDaE46

The system is stable if kkD>a; however, the steady-state error might be significant, since the controller gain kDcannot be made arbitrarily large.

4.2. PI controller

The next option is to consider a PI controller such as

Ks=kPIs+asE47

In this case, the closed-loop transfer function becomes

Ts=kkPIs2as+kkPIE48

Since ais a positive value, we observe that the PI controller cannot stabilize the plant Ps.

4.3. PID controller

Let us consider a PID controller—having a single tuning parameter KPID—in the form

Ks=kPIDs+a2sE49

The block diagram of the control system is shown in Figure 10.

Figure 10.

The block diagram of the control system.

The closed-loop transfer function is

Ts=kkPIDs+as2+kkPIDas+kkPIDaE50

The closed-loop transfer function has a zero at s=a, which might affect the system response (large overshoot), which can be canceled with a prefilter KPREs=a/s+a.

Then, the closed-loop transfer function becomes

Ts=kkPIDas2+kkPIDas+kkPIDaE51

Next, based on performance and robustness specifications, we would like to find a suitable value for the controller gain KPID. Usually, performance specifications are given in terms of settling time Tsetand percent of overshoot P.O.

For a second-order system

Ts=ωn2s2+2τωns+ωn2E52

where ωnis the natural frequency and the τis the damping factor; we have Tset4/τωnand P.O.=100eτπ/1τ.

Since we have only one tuning parameter KPID, the performance specifications are given only in terms of settling time Tset=1.25s. Therefore, the PID controller gain can be calculated as

kkPIDa=2τωn8TsetE53

Thus, we obtain

kPID=8+aTsetkTsetE54

Next, the stability and robustness in a classical framework can be assessed. We calculate the gain and phase margins, obtaining Gm=0.47m/Aand Pm=37degfor the equilibrium position z0=3103m.

4.4. Gain-scheduled controller

We have seen that the actuator can be stabilized around, and equilibrium point and performance and robustness can be guaranteed.

However, we observed that the plant parameters kand aare varying with the operating point i0z0; thus, for good performance and robustness, the controller should take into account that the plant parameters are varying.

A survey of linear parameter-varying control applications can be found in [28], and control applications validated by experiments are presented in [27, 29] for actuators and for medical X-ray systems in [30]. High-accuracy mathematical modeling and a linear parameter-varying observer for fault detection and fault isolation are presented in [17].

We can design a linear parameter-varying controller having the form

Ksz0=kPIDz0s+az02sE55

where both the controller gain and controller zero are dependent on the equilibrium position.

Finally, stability and robustness (quadratic stability) of such a control system can be analyzed using modern software tools, and details are described in [27].

If a simpler approach is preferred, a gain-scheduled controller might be a good choice, which is easier to implement in real time, and its stability and robustness are easier to analyze.

Since perfect cancelation of the varying plant pole at s=awith a fixed controller zero is not possible, we choose to place the controller fixed zero left to the varying poles, such as z=amax. Then, the gain-scheduled controller can be written as

Ksz0=kPIDz0s+amax2sE56

where the only one tuning parameter is controller gain KPIDz0, defined as

kPIDz0=8+az0Tsetkz0TsetE57

where the variation of the values k=kz0and a=az0were shown already in Figure 5.

Next, the gain and phase margins are calculated for different equilibrium positions z0and shown in Figure 11. We can observe that—due to the gain-scheduled controller—the gain and phase margins do not change significantly with the equilibrium position, and such a robustness is difficult to achieve with a single controller, having fixed parameters.

Figure 11.

Gain and phase margin variation with z 0 .

The system response using the gain-scheduled controller is investigated, considering the following two cases:

  • The moving armature is controlled around an equilibrium position, and the set point (reference position) is changed Δz0=1mm(see Figure 12, left plot). The control system—including the prefilter—exhibits approximately P.O.25%overshoot and settling time Tset=1.5s.

  • The moving armature is controlled between the two extreme positions, armature open and armature close, Δz0=10mm(see Figure 12, right plot). In this case the main goal is to achieve so-called soft landing of the moving armature to reduce wear and noise.

Figure 12.

System response with the gain-scheduled controller.

It is important to highlight that controller design is made based on the linearized plant, but validation of the controller in simulations or during hardware-in-the-loop (HIL) experiments shall be done using the nonlinear plant model.

During our control design investigations, we considered that the armature position can be measured. In practice, there are applications, where the armature position cannot be measured in a cost-effective way.

Therefore, we remark that controlling the moving armature without measuring the armature position (e.g., measuring only the current) remains a challenging research topic, which exceeds the goals and the limits of this chapter.

5. Conclusions

This chapter dealt with mathematical modeling, system identification, and control of electromagnetic actuators. Actuators are often used in industrial applications such as magnetic levitation, electromagnetic bearings, as well as in fuel injectors in the automotive industry.

After a detailed mathematical model was presented, two different parameter identification techniques were described. The first one is based on the classical current decay test, and the second one is a clustering-based system identification approach. Since the actuator is open-loop unstable, the main steps of system identification of the actuators under closed-loop control were presented.

Finally, very simple and easy-to-apply control strategies were discussed, when the armature is controlled around a fixed equilibrium position (PID controller) as well as when the armature is controlled between two extreme positions, armature open and armature closed (gain-scheduled PID controller).

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Alexandru Forrai (July 11th 2018). Modeling, System Identification, and Control of Electromagnetic Actuators, Actuators, Constantin Volosencu, IntechOpen, DOI: 10.5772/intechopen.75088. Available from:

chapter statistics

308total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Quantitative Tactile Examination Using Shape Memory Alloy Actuators for the Early Detection of Diabetic Neuropathy

By Junichi Danjo, Sonoko Danjo, Hideyuki Sawada, Keiji Uchida and Yu Nakamura

Related Book

First chapter

Microassembly Using Water Drop

By Taksehi Mizuno

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us