Open access peer-reviewed chapter

Sensorless Control of Induction Motor Supplied by Current Source Inverter

By Marcin Morawiec

Submitted: July 4th 2011Reviewed: June 28th 2012Published: November 14th 2012

DOI: 10.5772/51133

Downloaded: 2148

1. Introduction

The circuits used for power conversion applied in drives with induction motor (IM) are classified into two groups: voltage source inverters (VSI) and current source inverters (CSI). The VSI were used more often than the CSI because of their better properties. Nowadays, the development of power electronics devices has enormous influence on applications of systems based on the CSI and creates new possibilities.

In the 1980s the current source inverters were the main commonly used electric machine feeding devices. Characteristic features of those drives were the motor electromagnetic torque pulsations, the voltage and current with large content of higher harmonics. The current source inverter was constructed of a thyristor bridge and large inductance and large commutation capacitors. Serious problems in such drive systems were unavoidable overvoltage cases during the thyristor commutation, as the current source inverter current is supplied in a cycle from a dc-link circuit to the machine phase winding. The thyristor CSI has been replaced by the transistor reverse blocking IGBT devices (RBIGBT), where the diode is series-connected and placed in one casing with transistor. The power transistors like RBIGBT or Silicon Carbide (SiC) used in the modern CSIs guarantee superior static and dynamic drive characteristics.

The electric drive development trends are focused on the high quality system. The use of current sources for the electric machine control ensures better drive properties than in case of voltage sources, where it may be necessary to use an additional passive filter at the inverter output. The Pulse width modulation (PWM) with properly chosen dc-link inductor and input-output capacitors result in sinusoidal inverter output currents and voltages. Methods of calculating proper inductance in dc-link were proposed in [Glab (Morawiec) M. et. al., 2005, Klonne A. & Fuchs W.F., 2003, 2004]. Properties of dc-link circuit of the Current Source Converter (CSC) force the utilization of two fully-controlled inverters to supply system with electric motor. The first of them – CSI - generates the current output vector to supply the induction motor. The second one – Current Source Rectifier (CSR) - generates a DC voltage to supply a dc-link circuit. The strategy for controlling the output current vector of CSI can be realized in two ways [Glab (Morawiec) M. et. al., 2005, Klonne A. & Fuchs W.F., 2003, Kwak S. & Toliyat H.A., 2006]. First of them is based on changes of modulation index while the value of current in dc-link circuit remains constant [Klonne A. & Fuchs W.F., 2003]. The second method is based on changes of dc-link current. In this case the CSI is working with constant, maximum value of PWM modulation index. Control of modulation index in CSI is used in drive systems, where high dynamic of electromagnetic torque should be maintained [Klonne A. & Fuchs W.F., 2003]. High current in dc-link circuit is a reason for high power losses in CSI. The simplified control method is the scalar control: current to slip (I/s). This method is very simple to implement, but the drive system has average performance (only one controller is necessary, the current in dc-link is kept at constant value by PI controller).

The drive system quality is closely connected with the machine control algorithm. The space vector concept, introduced in 1959 by Kovacs and Racz, opened a new path in the electric machine mathematical modelling field. The international literature on the subject presents drive systems with the CSI feeding an induction motor with the control system based on the coordinate system orientation in relation to the rotor flux vector (FOC – Field Oriented Control). Such control consisted of the dc-link circuit current stabilization [Klonne A. & Fuchs W.F., 2003]. In such control systems the control variables are the inverter output current components. This control method is presented in [Nikolic Aleksandar B. & Jeftenic Borislav I.], where the authors analyze control system based on direct torque control. The control process where the control variable is the inverter output current may be called current control of an induction motor supplied by the CSI.

Another control method of a current source inverter fed induction motor is using the link circuit voltage and the motor slip as control variables. That type of control may be called voltage control of a CSI fed induction motor, as the dc-link circuit voltage and angular frequency of current vector are the control variables. Proposed control strategy bases on nonlinear multi-scalar control [Glab (Morawiec) M. et. al., 2005, ]. The nonlinear control may result in better properties in case if the IM is fed by CSI. To achieve independent control of flux and rotor speed, new nonlinear control scheme is proposed. In this control method the inverter output currents are not controlled variables. The voltage in dc-link and pulsation of output current vector are the controlled variables which can be obtained by nonlinear transformations and are proposed by authors in [Krzeminski Z., 1987, Glab (Morawiec) M. et. al., 2005, 2007]. The multi-scalar model is named the extended one because the mathematical model contained dc-link current and output capacitors equations. This full mathematical model of induction machine with the CSI is used to derive new multi-scalar model. In proposed method the output current vector coefficients are not controlled variables. The output current vector and the flux vector are used to achieve new multi-scalar variables and new multi-scalar model. The control system structure may be supported on PI controllers and nonlinear decouplings or different controllers e.g. sliding mode controllers, the backstepping control method or fuzzy neural controller.

2. Structure of the drive system supplied by CSI or CSC

The simplified configuration of the drive system with the CSI is presented in the Fig. 1. The integral parts of the system are the inductor in dc-link and the output capacitors. In the Fig. 1 the structure with the chopper as an adjustable voltage source is presented.

Figure 1.

The CSI with the chopper

Figure 2.

The Current source converter

The chopper with the small inductor Ld (a few mH) forms the large dynamic impedance of the current source. In the proposed system the transistors forms commutator which transforms DC current into AC current with constant modulation index. The current is controlled by voltage source ed in dc-link. In this way the system with CSI remains voltage controlled and the differential equation for dc-link may be integrated with differential equation for the stator. The inductor limits current ripples during commutations of transistors. The transistors used in this structure are named the reverse blocking IGBT transistors (RBIGBT).

In order to avoid resonance problem the CSI or CSC structure parameters (input-output capacitors and inductor) ought to be properly chosen. The transistor CSI or CSC structures should guarantee sinusoidal stator current and voltage of IM if the parameters are selected by iteration algorithm.

2.1. The iteration algorithm selection of the inductor and input-output capacitors

The inductance in dc-link Ld may be calculated as a function of the integral versus time of the difference of input voltage ed and output voltage ud in dc-link [Glab (Morawiec) M. et. al., 2005, Klonne A. & Fuchs W.F., 2003, 2004]. Calculating an inductance from [Glab (Morawiec) M. et. al., 2005, Klonne A. & Fuchs W.F., 2003] may be not enough because of the resonance problem. The parameters could be determined by simple algorithm.

Two criteria are taken into account:

  • Minimization of currents ripples in the system

  • Minimization of size and weight.

The first criteria can be defined as:

wi=Δidmaxid,E1

where

Δidmax is max(idmax(t1) – idmin(t2)),

id is average value of dc-link current in one period.

The current ripples in dc-link has influence on output currents and commutation process. According to this, the wi factor, THDi stator and THDu stator must be taken into account. Optimal value of inductance Ld and output CM ensure performance of the drive system with sinusoidal output current and small THD. In Fig. 3 the iteration algorithm for choosing the inductor and capacitor is shown. In every step of iteration new values wi, THDi, THDu are received. In every of these steps new values are compared with predetermined value wip, THDip, THDup and:

N is number of iteration,

THDi – stator current total harmonic distortion,

THDu – stator voltage total harmonic distortion. Number of iteration is set for user.

In START the initial parameters are loaded. In block Set Ld inductance of the inductor is set. In Numerical process block the simulation is started. In next steps THDi, THDu and wi coefficient are calculated. THDi, THDu and wi coefficient are compared with predetermined value. If YES then CM is setting, if NO the new value of Ld must be set. Comparison with predetermined value is specified as below:

{THDipmaxTHDi(i)THDipminTHDupmaxTHDu(i)THDupminΔidpmaxΔid(i)Δidpmin}E2

where

Δidpmax – maximum value for wi coefficient,

Δidpmin – minimum value for wi coefficient,

THDipmax, THDupmax – maximum predetermined value of THD for range (ΔLd, ΔCM),

THDipmin, THDupmin – minimum predetermined value of THD (ΔLd, ΔCM),

ΔLd – interval of optimal value Ld,

ΔCM – interval of optimal value CM.

For optimal quality of stator current and voltage in a drive system THDi ought to be about 1%, THDu<2% and wi<15% in numerical process. Estimated CSC parameters by the iteration algorithm are shown in Fig. 4 and 5 or Table 1.

Figure 3.

The iteration algorithm for selection of the inductor Ld an capacitors CM

Figure 4.

Inductor inductance from iteration algorithm, where Pn [kW] is nominal machine power for different transistors switching frequency [kHz]

Figure 5.

Capacitor capacitance from iteration algorithm where Pn [kW] is nominal machine power for different transistors switching frequency [kHz]

The value of AC side capacitors CL ought to be about 25% higher than for CM because of higher harmonics in supply network voltage:

CL1,25CME3
.
Pn [kW]Ld[mH]CM[μF]CL[μF]Pn[kW]Ld[mH]CM[μF]CL[μF]
1,513,21010157,63035
2,212,51212226,25060
411,62020305,56070
5,510,52020454,58090
7,59,42222553120150
118,32225752150200

Table 1.

Estimated a CSC parameters

3. The mathematical model of IM supplied by CSC

3.1. Introduction to mathematical model

Differential equation for the dc-link is as follows

ed=idRd+Lddiddτ+udE4
,

where: ud is the inverter input voltage, Rd is the inductor resistance, Ld is the inductance, ed is the control voltage in dc-link, id is the current in dc-link.

Equation (4) is used together with differential equation for the induction motor to derive the models of induction motor fed by the CSI.

The model of a squirrel-cage induction motor expressed as a set of differential equations for the stator-current and rotor-flux vector components presented in αβ stationary coordinate system is as follows [Krzeminski Z., 1987]:

disαdτ=RsLr2+RrLm2Lrwσisα+RrLmLrwσψrα+ωrLmwσψrβ+Lrwσusα,E5
disβdτ=RsLr2+RrLm2Lrwσisβ+RrLmLrwσψrβωrLmwσψrα+Lrwσusβ,E6
dψrαdτ=RrLrψrαωrψrβ+RrLmLrisα,E7
dψrβdτ=RrLrψrβ+ωrψrα+RrLmLrisβ,E8
dωrdτ=LmJLr(ψrαisβψrβisα)1Jm0,E9

where

Rr, Rs are the motor windings resistance, Ls, Lr, Lm are stator, rotor and mutual inductance, u, u, i, i, ψ, ψ are components of stator voltage, currents and rotor flux vectors, ωr is the angular rotor velocity, J is the torque of inertia, m0 is the load torque. All variables and parameters are in p. u.

3.2. The mathematical model of IM contains full drive system equations

The vector components of the rotor flux together with inverter output current are used to derive model of IM fed by the CSI. The model is developed using rotating reference frame xy with x axis orientated with output current vector. The y component of the output current vector is equal to zero.

The variables in the rotating reference frame are presented in Fig. 6.

The output current under assumption an ideal commutator can be expressed

if=Kid,E10

where

K is the unitary commutation function (K=1).

Figure 6.

Variables in the rotating frame of references

If the commutation function is K=1 than

|if|id.E11

The equation (12) results from (11) taking into account ideal commutator of the CSI, according to equation

pDClink=pAC motor sideudid=usxifxE12
,

where: ud is the input six transistors bridge voltage, usx is the stator voltage component.

The full model of the drive system in rotating reference frame xy with x axis oriented with inverter output current vector is as follows

disxdτ=RsLr2+RrLm2Lrwσisx+RrLmLrwσψrx+ωiisy+ωrLmwσψry+Lrwσusx,E13
disydτ=RsLr2+RrLm2Lrwσisy+RrLmLrwσψryωiisxωrLmwσψrx+Lrwσusy,E14
dψrxdτ=RrLrψrx+(ωiωr)ψry+RrLmLrisx,E15
dψrydτ=RrLrψry(ωiωr)ψrx+RrLmLrisy,E16
dωrdτ=LmJLr(ψrxisyψryisx)1Jm0,E17
diddτ=edLdRdLdidusxLd,E18
dusxdτ=1CM(ifxisx)+ωiusyE19
,
dusydτ=1CMisyωiusxE20
.

where: ωi is angular frequency of vectorif, isx, isy are the capacitors currents.

4. The nonlinear multi-scalar voltage control of IM with PI controllers

4.1. The simplified Multi-scalar control

The Nonlinear multi-scalar control was presented by authors [Krzeminski Z., 1987, Glab (Morawiec) M. et. al., 2005, 2007]. This control in classical form based on PI controllers. The simplify multi-scalar control of IM supplied by CSC for different vector components (ψr,is), (ψs,is), (ψm,is) was presented in [Glab (Morawiec) M. et. al., 2005, 2007]. These multi-scalar control structures give different dynamical and statical properties of IM supplied by CSI. In this chapter the simplified control is presented. The simplification is based on (11) and (12) equations. If the capacity CM has small values (a few μF) the mathematical equations (19) -(20) can be ommitted and the output current vector in stationary state is|if||is|. Under this simplification, to achieve the decoupling between two control paths the multi-scalar model based control system was proposed [Krzeminski Z., 1987, Glab (Morawiec) M. et. al., 2005, 2007]. The variables for the multi-scalar model of IM are defined

x11=ωr,E21
x12=idψry,E22
x21=ψrx2+ψry2,E23
x22=idψrx,E24

where

x11 is the rotor speed, x12 is the variable proportional to electromagnetic torque, x21 is the square of rotor flux and x22 is the variable named magnetized variable [Krzeminski, 1987].

Assumption of such machine state variables may lead to improvement of the control system quality due to the fact that e.g. the x12 variable is directly the electromagnetic torque of the machine. In FOC control methods [Klonne A. & Fuchs W.F., 2003, 2004, Salo M. & Tuusa H. 2004] the electromagnetic torque is not directly but indirectly controlled (the isq stator current component). With the assumption of a constant rotor flux modulus, such a control conception is correct. The inaccuracy of the machine parameters, asymmetry or inadequately aligned control system may lead to couplings between control circuits.

The mathematical model for new state of variables (21) - (24) used (15) - (18) is expressed by differential equations:

dx11dτ=LmJLr(x12)1Jm0,E25
dx12dτ=1Tix12+1LdusxψryRrLmLrisyid+v1,E26
dx21dτ=2RrLrx21+2RrLmLrx22,E27
dx22dτ=1Tix22+RrLmLrisxid+v2E28
.

The compensation of nonlinearities in differential equation leads to the following expressions for control variables v1 and v2 appearing in differential equations (27) - (28):

v1=RrLmLrisyidusxLdψry+1Tim1,E29
v2=RrLmLrisxid+usxLdψrx+1Tim2E30
,

where m1, 2 are the PI controllers output and

1Ti=RsLr+RdLdE31
.

The control variables are specified

ed=Ldψryv1ψrxv2x21E32
,
ωi=ψrxv1+ψryv2idx21+x11E33
,when
x21,id0E34
.

The inverter control variables are: voltage ed and the output current vector pulsation. The multi-scalar control of IM supplied by CSI was named voltage control because the control variable is voltage ed in dc-link.

The decoupled two subsystems are obtained:

  • electromagnetic subsystem

dx21dτ=2RrLrx21+2RrLmLrx22,E35
dx22dτ=1Ti(x22+m2)E36
,
  • electromechanical subsystem

dx11dτ=LmJLrx121Jm0,E37
dx12dτ=1Ti(x12+m1).E38
(37)

4.2. The multi-scalar control with inverter mathematical model

The author in [Morawiec M., 2007] revealed stability proof of simplified multi-scalar control while the parameters of the CSI are optimal selected.

When the capacitance CM is neglected the stator current vector isis about ~5% out of phase to ifwhile nominal torque is set. Then the control variables and decoupling are not obtained precisely. The error is small than 2% because PI controllers improved it.

In order to compensate these errors the capacity CM to mathematical model is applied.

From (19) - (20) in stationary state lead to dependences:

isx=ifx+ωifCMusyE39
,
isy=ωifCMusxE40
.

The new mathematical model of the drive system is obtained from (38) - (39) through differentiation it and used (15) - (16) in xy coordinate system:

disxdτ=RdLdid+1Lded1LduxωifCMisyωif2CM2usxE41
,
disydτ=ωifCMid+ωifCMisxωif2CM2usyE42
,
dψrxdτ=RrLrψrx+(ωifωr)ψry+RrLmLrisxE43
,
dψrydτ=RrLrψry(ωifωr)ψrx+RrLmLrisyE44
,
diddτ=RdLdid+1Lded1LdusxE45
,
dusxdτ=1CM(ifxisx)+ωifusyE46
,
dusydτ=1CMisyωifusxE47
.

Substituting (38) - (39) to multi-scalar variables [Krzeminski Z., 1987] one obtains:

x11=ωrE48
,
x12=ifxψryωifCMx32E49
,
x21=ψrx2+ψry2E50
,
x22=ifxψrx+ωifCMx31E51
,and
x31=ψrxusyψryusxE52
,
x32=ψrxusx+ψryusyE53
.

The multi-scalar model for new multi-scalar variables has the form:

dx11dτ=LmJLrx121Jm0,E54
dx12dτ=1Tix12+1Ldusxψryx11x22+v1E55
,
dx21dτ=2RrLrx21+2RrLmLrx22,E56
dx22dτ=1Tix221Ldusxψrx+RrLmLrid2+x11x12+v2E57
,

where

v1=1Ldedψry+ωif(x22RdLCMx32CMRrLmLrps)E58
,
v2=1Ldedψrx+ωif(x12+RdLCMx31+CMRrLmLrqs)E59
,
qs=isxusyisyusxE60
,
ps=usxisx+usyisyE61
.

The compensation of nonlinearities in differentials equation leads to the following expressions for control variables v1 and v2 appearing in differential equations (56), (58):

v1=1Tim11Ldusxψry+x11x22E62
,
v2=1Tim2+1Ldusxψrxx11x12RrLmLrid2E63
,

and the control variables

ed=LdV2x41V1x42ψrxx41+ψryx42E64
,
ωi=V1ψrx+V2ψryψrxx41+ψryx42E65
,

where

x41=x22RdLCMx32CMRrLmLrpsE66
,
x42=x12+RdLCMx31+CMRrLmLrqsE67
,1Tiis determined in (31).

The decoupled two subsystems are obtained as in (34) - (37).

4.3. The multi-scalar adaptive-backstepping control of an IM supplied by the CSI

The backstepping control can be appropriately written for an induction squirrel-cage machine supplied from a VSI. In literature the backstepping control is known for adaptation of selected machine parameters, written for an induction motor [Tan H. & Chang J., 1999, Young Ho Hwang, 2008]. In [Tan H. & Chang J., 1999, Young Ho Hwang, 2008] the authors defined the machine state variables in the dq coordinate system, oriented in accordance with the rotor flux vector (FOC). The control method presented in [Tan H. & Chang J., 1999, Young Ho Hwang, 2008] is based on control of the motor state variables: ωr – rotor angular speed, rotor flux modulus and the stator current vector components: isd and isq. Selection of the new motor state variables, as in the case of multi-scalar control with linear PI regulators, leads to a different form of expressions describing the machine control and decoupling. The following state variables have been selected for the multi-scalar backstepping control

e1=x11*x11E68
,
e2=x12*x12E69
,
e3=x21*x21E70
,
e4=(2RrLmLrx22)*2RrLmLrx22E71
,

where: x11, x12, x21 and x22 are defined in (47) - (50).

The e4 tracking error is defined in (72), it does not influence on the control system properties and is only an accepted simplification in the format of decoupling variables.

Derivatives of the (69) - (72) errors take the form

e˙1=LmJLre2k1e1m˜0JE72
,
e˙2=k1e2k12JLrLme1+LrLmm^˙0+1Tix121Ldusxψry+x11x22+RrLmLrisyidv1E73
,
e˙3=k3e3+e4E74
,
e˙4=k32e3+k3e44(RrLr)2x21+4(RrLr)2Lmx22+2RrLmLrTix22+2RrLmLrLdusxψrx2(RrLmLr)2isxid+2RrLmLrx11x122RrLmLrv2E75
.

The Lyapunov function derivative, with (73) – (76) taken into account, may be expressed as:

V˙=k12e1k22e2k32e3k42e4+e2(f1v1)+e4(f2a3v2)++m˜0(e1Jk1LrLme2+m^˙0γ)E76
,

where

v1=1Ldedψry+ωifx41E77
,
v2=1Ldedψrx+ωifx42E78
,
f1=limit12e1(LmJLrk12JLrLm)+k2e2+k1e2+LrLmm^˙0+1Tix121Ldusxψry++x11x22E79
,
f2=limit12(e3k32e3)+k4e4+k3e44(RrLr)2x21+4(RrLr)2Lmx22+2RrLmLrTix22++2RrLmLrLdusxψrx2(RrLmLr)2id22RrLmLrx11x12E80
,
a3=2RrLmLrE81
.

limit12 – is a dynamic limitation in the motor speed control subsystem,

limit22 – is a dynamic limitation in the rotor flux control subsystem,

k1…k4 and γ are the constant gains.

The control variables take the form:

ed=Ldx41f2a3x42f1a3(ψrxx41+ψryx42)E82
,
ωi=a3ψrxf1+ψryf2a3(ψrxx41+ψryx42)E83
.

The inverter control variables are: voltage ed and the output current vector pulsation. The two decoupled subsystems are obtained as in (34) - (37).

The load torque m0 can be estimated from the formula:

m^˙0=γ(e1J+k1LrLme2)E84
.(84)

4.4. Dynamic limitations of the reference variables

In control systems with the conventional linear controllers of the PI or PID type, the reference (or controller output) variable dynamics are limited to a constant value or dynamically changed by (85) - (86), depending on the drive working point.

Control systems where the control variables are determined from the Lyapunov function (like in backstepping control) have no limitations in the set variable control circuits. The reference variable dynamics may be limited by means of additional first order inertia elements (e.g. on the set speed signal).

The author of this paper has not come across a solution of the problem in the most significant backstepping control literature references, e.g. [Tan H. & Chang J., 1999, Young Ho Hwang, 2008]. In the quoted reference positions, the authors propose the use of an inertia elements on the set variable signals. Such approach is an intermediate method, not giving any rational control effects. The use of an inertia element on the reference signal, e.g. of the rotor angular speed, will slow down the reference electromagnetic torque reaction in proportion to the inertia element time-constant. In effect a "slow" build-up of the motor electromagnetic torque is obtained, which may be acceptable in some applications. In practice the aim is to limit the electromagnetic torque value without an impact on the build-up dynamics. Control systems with the Lyapunov function-based control without limitation of the set variables are not suitable for direct adaptation in the drive systems. Therefore, a solution often quoted in literature is the use of a PI or PID speed controller at the torque control circuit input.

The set values of thex12*, x22*variables appearing in the e2 and e4 deviations can be dynamically limited and the dynamic limitations are defined by the expressions [Adamowicz M.; Guzinski J., 2005]:

x12lim=Ismax2x21x222E85
,
x22lim=f(Usmax2,Ismax2,x11)E86
,

where

x12lim – the set torque limitation,

x22lim – the x22 variable limitation,

Ismax – maximum value of the stator current modulus,

Usmax – maximum value of the stator voltage modulus.

The above given expressions may be modified to:

x12lim=Ismax2x21x212Lm2E87
,

giving the relationship between the x21 variable, the stator current modulus Ismax, and the motor set torque limitation.

For the multi-scalar backstepping control, to the f1 and f2 variables the limit12 and limit22 variables were introduced; they assume the 0 or 1 value depending on the need of limiting the set variable.

Limitation of variables in the Lyapunov function-based control systems may be performed in the following way:

if (x12*>x12lim) then  {limit12=0,e2=x12limx12}E88
,
if (x12*<x12lim) then {limit12=0, e2=x12limx12}E89
,else
limit12=1E90
,
if (x22*>x22lim) then  {limit22=0,e4=x22limx22}E91
,
if (x22*<x22lim) then {limit22=0, e4=x22limx22}E92
,else
limit22=1E93
.

The dynamic limitations effected in accordance with expressions (85) – (86) limit properly the value of x12*and x22*variables without any interference in the reference signal build-up dynamics.

Fig. 7 presents the variable simulation diagrams. The backstepping control dynamic limitations were used.

Figure 7.

Diagrams of multi-scalar variables in the machine dynamic states, the x12ogr = 1.0 and x22ogr = 0.74 limitations were set for a drive system with an induction squirrel-cage machine supplied from a CSC-simulation diagrams, x*12 – diagram of the machine set electromagnetic torque (without signal limitation), x*22 – diagram of the x*22 set signal (without limitation).

4.5. Impact of the dynamic limitation on the estimation of parameters

The use of a variable limitation algorithm may have a negative impact on the control system estimated parameters. This has a direct connected with the limited deviation values, which are then used in an adaptive parameter estimation. Such phenomenon is presented in Fig. 7. The estimated parameter in the control system is the motor load torquem^0. The set electromagnetic torque is limited to the x12lim = 1.0 value. Fig. 7 shows that the estimated load torque increases slowly in the intermediate states. Limitation of the set electromagnetic torque causes the limitation of deviation e2, which in turn causes limited increase dynamics of the estimated load torque. The m^0value for limit12 = 0 in the dynamic states does not reach the real value of the load torque, which should bem^0x12. A large m˜0estimation error occurs in the intermediate states, which can be seen in Fig. 8. The estimation error in the intermediate states is m˜00because the torque limitation, introduced to the control system, is not compensated. The simulation and experimental tests have shown that the load torque estimation error in the intermediate state has an insignificant impact on the speed control. Omitting m^0in the set torque x*12 expression eliminates the intermediate state speed over-regulation. But absence of m^0in x*12 for a steady state gives the deviation value e10and lack of full control over maintaining the rotor set angular speed. Compensation of the limit12 limitation introduced to the control system is possible by installing a corrector in the rotor angular speed control circuit.

A corrector in the form of an e1 signal integrating element was added to the set electromagnetic torque x*12 signal. In this way a system reacting to the change of machine real load torque was obtained. The introduced correction minimizes the rotor angular speed deviation and the corrector signal may be treated as the estimated load torque value.

The correction element is determined by the expression:

KTL=ke1tk1tke1dτE94
,

where

tk-1…tk is the e1 signal integration range,

KTL – correction element,

ke1– is the correction element amplification.

The gain ke1 should be adjusted that the speed overregulation in the intermediate state does not exceed 5%:

0<ke10,1k1E95
,

The correction element amplification must not be greater than k1, or:

ke1k1E96
.For ke1>k1the KTL signal will become an oscillation element and may lead to the control system loss of stability.

The KTL signal must be limited to the x12lim value.

The x*12 set value expression must be modified:

x12*=JLrLmk1e1+KTLE97
,

where

m^0KTLE98
.

The use of (95) in the angular speed control circuit improves the load torque estimation and eliminates the steady state speed error.

Fig. 9 presents the load torque (determined in (96)) estimation as well as x12 and the limit12 limitations.

Figure 8.

Impact of the electromagnetic torque limitation x12lim on the estimated load torque m^0(84).

Figure 9.

Diagrams of the limit12 variable, KTLload torque and electromagnetic torque x12..

5. The nonlinear multi-scalar current control of induction machine

Conception of the CSI current control is based on forced components of the CSI output current. The dc-link circuit inductor could be modeled as the first order inertia element with the time constant T of a value equal to the dc-link circuit time constant value. The dc-link equation may be introduced to the induction machine mathematical model to obtain the set CSI output current component. The time constant T is equal the inductance Ld, it can be written:

disdt=1T(ifis)E99
,
T=RdLdE100
.

5.1. The simplified multi-scalar current control of induction machine

The simplified version of the CSI output current control it is assumed that the output capacitors have negligibly small capacitance, so their impact on the drive system dynamics is small. Assuming that the cartesian coordinate system, where the mathematical model variables are defined, is associated with the CSI output current vector (which with this simplification is the machine stator current) the mathematical model can be obtained (97) and (42) - (43) equations).

The multi-scalar variables have the form

x11=ωrE101
,
x12=isxψrxE102
,
x21=ψrx2+ψry2E103
,
x22=isxψrxE104
,

where

isxis treated as the output current vector component and|if||is|.

For those variables, the multi-scalar model has the form:

dx11dτ=LmJLrx121Jm0E105
,
dx12dτ=1Tix12+v1E106
,
dx21dτ=2RrLrx212RrLmLrx22E107
,
dx22dτ=1Tix22+RrLmLrid+v2E108
.

Applying the linearization method, the following relations are obtained, where m1 is the subordinated regulator output in the speed control line and m2 is the subordinated regulator output in the flux control line

v1=1Tim1E109
,
v2=1Tim2RrLmLridE110
.

The control variables are modulus of the CSI output current and the output current vector pulsation, given by the following relations:

|if|=Tv2ψrxv1ψryx21ψryE111
,
ωi=v1+v2x21isx+x11E112
.

where: Ld inductance, Ti– the system time constant.

5.2. The multi-scalar current control of induction machine

The current control analysis presented in the preceding sections does not take the CSI output capacitors into account. Such simplification may be applied because of the small impact of the capacitors upon the control variables (the machine stator current and voltage are measured). The capacitor model will have a positive impact on the control system dynamics.

The output capacitor relations have the form:

dusdt=1CM(ifis)E113
,

where: usis the capacitor voltage vector, ifis the current source inverter output current vector, isis the stator current vector.

Using the approximation method, relation (38) may be written as follows:

us(k)us(k1)Timp=1CM[if(k)is(k)]E114
.

Deriving us(k) from (38), the motor stator voltage is obtained as a function of the output current, stator current and stator voltage, in the form:

us(k)=TimpCM[if(k)is(k)]+us(k1)E115
,

where

us(k)is the stator voltage vector at the k-th moment, Timp - sampling period.

In the equations (5) - (9) representing the cage induction motor mathematical model the stator current vector components are appeared but the direct control variables do not. The motor stator current vector components cannot be the control variables because the multi-scalar model relations are derived from them. This is a different situation than with the FOC control. The FOC control is based on the machine stator current components described in a coordinate system associated with the rotor flux and the stator current components are the control variables. Therefore, control variables must be introduced into the mathematical model (5) - (9). The control may be introduced considering the machine currents (5) - (6) and equation (97) written for the αβ components and describing the dc-link circuit dynamics. Adding the respective sides of equations (5) and (97) and equations (6) and (97), where equation (97) must be written with the (αβ) components – the mathematical model of the drive system fed by the CSI is obtained:

disαdτ=RsLr2T+RrLm2T+Lrwσ2LrwσTisα+RrLm2Lrwσψrα+ωrLm2wσψrα+Lr2wσusα+12TifαE116
,
disβdτ=RsLr2T+RrLm2T+Lrwσ2LrwσTisβ+RrLm2LrwσψrβωrLm2wσψrα+Lr2wσusβ+12TifβE117
,

and equations (7) - (8).

The multi-scalar variables are assumed like in [Krzeminski Z., 1987]:

x11=ωrE118
,
x12=ψrαisβψrβisαE119
,
x21=ψrα2+ψrβ2E120
,
x22=ψrαisα+ψrβisβE121
.

Introducing the multi-scalar variables (116) - (119), the multi-scalar model of an IM fed by the CSI is obtained:

dx12dτ=12Tix12x11x22x11x22Lm2wσ+Lr2wσ(usβψrαusαψrβ)+v1E122
,
dx22dτ=12Tix22+x11x12+RrLm2Lrwσx21+RrLm2Lrisα2+Lr2wσ(ψrαusα+ψrβusβ)+v2E123
,

where

v1=12Tψrαifβ12TψrβifαE124
,
v2=12Tψrαifα+12TψrβifβE125
.

Applying the linearization method to (120) - (121), the following expressions are obtained:

v1=12Tim1+x11x22+Lm2wσx11x21+Lr2wσ[usβ(k1)ψrαusα(k1)ψrβ]a1x12E126
,
v2=12Tim2x11x21RrLm2wσLrx21RrLm2Lrisα2-Lr2wσ[usα(k1)ψrα+usβ(k1)ψrβ]+a1x22E127
,

where

v1=a1[ifβψrαifαψrβ]E128
,
v2=a1[ifαψrα+ifβψrβ]E129
.

The control variables take the form

ifα=a2ψrαv2ψrβv1x21E130
,
ifβ=a2ψrαv1+ψrβv2x21E131
,

where

a1=LrTimp2wσCME132
,
a2=a1+12TE133
.

The time constant Ti for simplified for both control method is given

Ti=wσLrTRsLr2T+RrLm2T+LrwσE134
,

5.3. Generalized multi-scalar control of induction machine supplied by CSI or VSI

A cage induction machine fed by the CSI may be controlled in the same way as with the voltage source inverter (VSI). The generalized control is provided by an IM multi-scalar model formulated for the VSI machine control [Krzeminski Z., 1987]. The (116) - (119) multi-scalar variables and additional u1 and u2variablesare used

u1=ψrαusβψrβusαE135
,
u2=ψrαusα+ψrβusβE136
,

which are a scalar and vector product of the stator voltage and rotor flux vectors.

The multi-scalar model feedback linearization leads to defining the nonlinear decouplings [Krzeminski Z., 1987]:

U1*=wσLr[x11(x22+Lmwσx21)+1Tvm1]E137
,
U2*=wσLr[x11x12RrLmLris2RrLmLrwσx21+1Tvm2]E138
.

The control variables for an IM supplied by the VSI have the form [Krzeminski Z., 1987]:

usα*=ψrαU2*ψrβU1*x21E139
,
usβ*=ψrαU1*+ψrβU2*x21E140
.

The controls (137) - (138) are reference variables treated as input to space vector modulator when the IM is supplied by the VSI.

On the other side, when the IM is fed by the CSI, calculation of the derivatives of (133) - (134) multi-scalar variables yields the following relations:

du1dτ=RrLru1x11u2+RrLmLrqs1CMx12+v11E141
,
du2dτ=RrLru2+x11u1+RrLmLrps+1CMx22+v22E142
,

where ps and qs are defined in (61) - (62).

By feedback linearization of the system of equations, one obtains

v11=RrLrvp1RrLmLrqs+1CMx12+x11u2E143
,
v22=RrLrvp2RrLmLrps+1CMx22x11u1E144
,

where

vp1 and vp2 are the output of subordinated PI controllers.

The control variables of the IM fed by the CSI have the form:

ifα=CMv11ψrβv22ψrαx21E145
,
ifβ=CMv11ψrα+v22ψrβx21E146
.

As a result, two feedback loops and linear subsystems are obtained (Fig. 15).

6. The speed observer backstepping

General conception of the adaptive control with backstepping is presented in references [Payam A. F. & Dehkordi B. M. 2006, Krstic M.; Kanellakopoulos I.; & Kokotovic P. 1995]. In [Krstic M.; Kanellakopoulos I.; & Kokotovic P. 1995] the adaptive back integration observer stability is proved and the stability range is given.

Proceeding in accordance with the adaptive estimator with backstepping conception, one can derive formulae for the observer, where only the state variables will be estimated as well as the rotor angular speed as an additional estimation parameter.

Treating the stator current vector components i^sα,βas the observer output variables (as in [Payam A. F. & Dehkordi B. M. 2006, Krstic M.; Kanellakopoulos I.; & Kokotovic P. 1995]) and vα,βas the new input variables, which will be determined by the backstepping method, one obtains:

di^sαdτ=a1i^sα+a5ψ^rα+ω^ra4ψ^rβ+a6usα+vαE147
,
di^sβdτ=a1i^sβ+a5ψ^rβω^ra4ψ^rα+a6usβ+vβE148
,
dψ^rαdτ=RrLrψ^rαω^rψ^rβ+RrLmLri^sαE149
,
dψ^rβdτ=RrLrψ^rβ+ω^rψ^rα+RrLmLri^sβE150
.

In accordance to the backstepping method, the virtual control must be determined together with the observer stabilizing variables. In that purpose, the new ςαand ςβvariables have been introduced and linked with the stator current estimation deviations (the integral backstepping structure [Krstic M.; Kanellakopoulos I.; & Kokotovic P. 1995]):

dς˜adτ=i˜sαE151
,
dς˜bdτ=i˜sβE152
.

The stator current vector component deviations are treated as the subsystem control variables [Payam A. F. & Dehkordi B. M. 2006, Krstic M.; Kanellakopoulos I.; & Kokotovic P. 1995]. Adding and deducting the stabilizing functions, one obtains:

dς˜αdτ=i˜sασα+σβE153
,
dς˜βdτ=i˜sβσα+σβE154
,

where

σα=c1ς˜α,
σβ=c1ς˜βE155
,

by introducing the deviation defining variable, one obtains:

zα=i˜sα+c1ς˜αE156
,
zβ=i˜sβ+c1ς˜βE157
.

Transformation of (154) - (155) leads to:

dς˜αdτ=zαc1ς˜αE158
dς˜βdτ=zβc1ς˜βE159
.

Calculation of the (154) - (155) deviation derivatives gives:

z˙α=a5ψ˜rα+a4[ω^rψ˜rβ+ω˜r(ψ^rβψ˜rβ)]+vα+c1i˜sαE160
,
z˙β=a5ψ˜rβa4[ω^rψ˜rα+ω˜r(ψ^rαψ˜rα)]+vβ+c1i˜sβE161
.

By selecting the following Lyapunov function

V=ς˜α2+ς˜β2+zα2+zβ2+ψ˜rα2+ψ˜rβ2+1γω˜r2E162
,

calculating the derivative and substituting the respective expressions, new correction elements can be determined, treated in the speed observer backstepping as the input variables. The Lyapunov function is determined for the dynamics of theςα,β, zα,βvariables and for the rotor flux components. Calculating the(160) derivative, one obtains:

V˙=c1ς˜α2c1ς˜β2c2zα2c2zβ2RrLrψ˜rα2RrLrψ˜rβ2+zα(a5ψ˜rα+ω^ra4ψ˜rβ+ω˜ra4(ψ^rβψ˜rβ)++vα+c1i˜sα+c2zα+ς˜α)+zβ(a5ψ˜rβω^ra4ψ˜rαω˜ra4(ψ^rαψ˜rα)+vβ+c1i˜sβ+c2zβ+ς˜β)+ψ˜rα(RrLrψ˜rαω^rψ˜rβω˜r(ψ^rβψ˜rβ))+ψ˜rβ(RrLrψ˜rβ+ω^rψ˜rα+ω˜r(ψ^rαψ˜rα)).E163

The input variables vα,β, resulting directly from (161), should include the estimated variables and the estimation deviations:

vα=a5ψ˜rαω^ra4ψ˜rβc1i˜sαc2zας˜αE164
,
vβ=a5ψ˜rβ+ω^ra4ψ˜rαc1i˜sβc2zβς˜βE165
.

Taking (162) - (163) into account, the deviation derivatives may be written in the form:

z˙α=ω˜ra4(ψ^rβψ˜rβ)c2zας˜αE166
,
z˙β=ω˜ra4(ψ^rαψ˜rα)c2zβς˜βE167
.

Using (164) - (165), the Lyapunov function may be written as follows:

V˙=c1ς˜α2c1ς˜β2c2zα2c2zβ2+ω˜ra4[zα(ψ^rβψ˜rβ)zβ(ψ^rαψ˜rα)+1γω^˙r]E168
.

The observer, defined by the (145) - (148) and (156) - (157) equations, is a backstepping type estimator.

In the (162) - (165) expressions the rotor flux deviations appear, which may be neglected without any change to the observer properties (145) - (148). Besides, the ς˜α,βdeviations in (162) - (163) may be zero, thus lowering the observer order. Assuming the simplifications, one obtains

vα=c1i˜sαc2zαE169
,
vβ=c1i˜sβc2zβE170
,and
ω^˙r=γa4(zβψ^rαzαψ^rβ)E171
.

where

c1, c2,γ are constant gains,

a4=Lmwσ, a5=RrLmLrwσ,
a6=LrwσE172
.

In Fig. 10, 11 the backstepping speed observer test is shown. When the load torque is set to ~-0.1 p.u. the rotor speed in backstepping observer is more precisely estimated than e.g. Krzeminski’s speed observer.

Figure 10.

The Speed observer test: the estimated rotor speed x11 is changed from 0.1 to -0.1 p.u., the rotor flux and stator current coefficients are shown

Figure 11.

The Speed observer test: the estimated rotor speed x11 in backstepping observer is changed from 0.1 to -0.1 p.u., the estimated rotor speed ω^r_K by Krzeminski’s speed observer [Krzeminski Z., 1999] and the multi-scalar variable: x12, x21, x22 are shown. The load torque m0 is set to -0.1 p.u.

7. The control system structures

In Fig. 12 and Fig. 14 the voltage and current multi-scalar control system structure is shown. These structures are based on four PI controllers and contain: the modulator, the speed observer and decouplings blocks.

Figure 12.

The voltage multi-scalar control system structure

Figure 13.

The voltage multi-scalar adaptive backstepping control system structure

Figure 14.

The current multi-scalar control system structure.

Figure 15.

Generalized Multi-scalar Control System of Induction Machine supplied by CSI or VSI.

In Fig. 14 the ed_ref value is determined in ed reference block. The ed reference block can be PI current controller or other controller.

In Fig. 13 the voltage multi-scalar adaptive backstepping control system structure is shown.

In Fig. 15 generalized multi-scalar control system structure is presented. This control structure is divided into two parts: the control system of IM fed by VSI and the control system of IM fed by CSI.

8. Experimental verification of the control systems

The tests were carried out in a 5.5 kW drive system. The motor parameters are given in Table 2 and the main per unit values in Table 3. In Fig. 16, 17 motor start-up and reverse for control system presented in chapter 4.1-4.2 are shown. In Fig. 18, 19 motor start-up and reverse for control system presented in chapter 5 are shown. Fig. 20,21 presents the same steady state like previous but for adaptive backstepping control system (chapter 4.3). Fig. 22 presents diagram of stator currents and voltages when motor is starting up for voltage control system (chapter 4.3). In Fig. 23 load torque setting to 0.7 p.u. for current control is presented. In Fig. 24, 25 the id current and the sinusoidal stator voltage and stator current are presented.

Figure 16.

Motor start-up (chapter 4.1- 4.2)

Figure 17.

Motor reverse (chapter 4.1- 4.2)

Figure 18.

Motor start-up (chapter 5)

Figure 19.

Motor reverse (chapter 5)

Figure 20.

Motor start-up (chapter 4.3)

Figure 21.

Motor reverse (chapter 4.3)

Figure 22.

The currents and voltages

Figure 23.

Load torque is set to 0.7 p.u. (chapter 5)

Figure 24.

id current and the stator voltage

Figure 25.

The stator current in stationary state.

where:

x11 is the rotor speed, x12 is the variable proportional to electromagnetic torque, x21 is the square of rotor flux and x22 is the additional variables, id is the dc-link current, usα,β are the capacitor voltage components, KTL is correction element (load torque), isα,β are the stator current components.

9. Conclusion

In this chapter two approaches to control of induction machine supplied by current source converter are presented. The first of them is voltage multi-scalar control based on PI controllers or backstepping controller. The voltage approach seems to be a better solution than the second one: current control, because the control system structure is more simple than the current control structure. The voltage in dc-link is the control variables obtained directly from decouplings. The current in dc-link is not kept at constant value but its value depend on induction machine working point. The current control gives higher losses in dc-link and higher transistor power losses than the voltage control. The power losses can be minimized by modulation index control method but the control system is more complicated. Both control systems lead to decoupling control path and sinusoidal stator current and voltage when space vector modulation of transistors is applied.

PARAMETERVALUE
Pn (motor power)5.5 kW
Un (phase to phase voltage)400 V
In (current)10.9 A
J (interia)0.0045 kgm2
nn (rotor speed)1500 rpm
PARAMETERPER UNIT VALUES
Rs (stator resist.)
Rr (rotor resist.)
Lm (mutual-flux induct.)
Ls (stator induct.)
Lr (rotor induct.)
0.045
0.055
1.95
2.05
2.05
Current Source Converter
C (capacitor in dc-link)
Rd (inductor resist.)
CM,L (input-output caps)
0.1
0.002
0.2

Table 2.

The motor drive system parameters

DEFINITIONDESCRIPTION
Ub=3Unbase voltage
Ib=Inbase current
zb=Ub/Ibbase impedance

Table 3.

Definition of per unit values

© 2012 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

Marcin Morawiec (November 14th 2012). Sensorless Control of Induction Motor Supplied by Current Source Inverter, Induction Motors - Modelling and Control, Prof. Rui Esteves Araújo, IntechOpen, DOI: 10.5772/51133. Available from:

chapter statistics

2148total chapter downloads

3Crossref citations

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

Minimizing Torque-Ripple in Inverter-Fed Induction Motor Using Harmonic Elimination PWM Technique

By Ouahid Bouchhida, Mohamed Seghir Boucherit and Abederrezzek Cherifi

Related Book

Frontiers in Guided Wave Optics and Optoelectronics

Edited by Bishnu Pal

First chapter

Frontiers in Guided Wave Optics and Optoelectronics

By Bishnu Pal

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