Advances in Power Electronics

Advances in Power Electronics / 2011 / Article

Research Article | Open Access

Volume 2011 |Article ID 273081 |

Parvathy Ayalur Krishnamoorthy, Kamaraj Vijayarajan, Devanathan Rajagopalan, "Application of Quadratic Linearization for the Control of Permanent Magnet Synchronous Motor", Advances in Power Electronics, vol. 2011, Article ID 273081, 11 pages, 2011.

Application of Quadratic Linearization for the Control of Permanent Magnet Synchronous Motor

Academic Editor: Henry S. H. Chung
Received30 Jun 2011
Accepted08 Aug 2011
Published16 Nov 2011


Many of the existing control methods for the permanent magnet synchronous motor (PMSM) either deal with steady state models or consider dynamic models under particular cases. A dynamic model of the PM machine allows powerful control-theoretic techniques such as linearization to be applied to the system. Existing exact feedback linearization of dynamic model of PMSM suffers from singularity issues. In this paper, we propose a quadratic linearization approach for PMSM based on the approximate linearization technique which does not introduce singularities. A MATLAB simulation is used to verify the effectiveness of the linearization technique proposed. Also, to account for higher-order and unmodelled dynamics of PMSM, tuning of the linearizing transformation is proposed and verified using simulation.

1. Introduction

Permanent magnet (PM) machines, particularly at low power range, are widely used in the industry because of their high efficiency. They have gained popularity in variable frequency drive applications. The merits of the machine are elimination of field copper loss, higher power density, lower rotor inertia, and a robust construction of the rotor [1].

Steady state models, such as in [2], have been used to develop control strategies for the permanent magnet synchronous motor. Vector control [3] enables independent control over the magnitude and angle of the current with respect to the rotor such that instantaneous control over torque is possible. But vector control implemented, in practice [4], assumes a constant flux in the closed loop and uses static models to effect the control scheme. The steady state models used to describe the machine behaviour do not capture explicitly the dynamics of the machine involved. Decoupled control [5], while considering the dynamics of the machine, corresponds to application of particular control strategies for machine control.

A dynamic model of a PM machine using direct and quadrature axis variables [1] allows more powerful general control theories to be brought to bear on the problem of the control of PM machine. One such control-theoretic approach is quadratic linearization of the PM machine model. Linearization of a nonlinear system allows a simple fixed controller to be applied for the control of the linearized system and yet get a uniform closed loop response for different reference and load conditions.

Bodson and Chiasson [6] have designed a controller using differential geometric method based on exact feedback linearization. The main drawback with this method and its variations [7] is that, even if the system is linearizable, the linearizability is subject to certain function of the state being nonsingular in a region of operation of the machine. This puts a constraint on the practical implementation of the control strategy. Zhu et al. [8] provide a static feedback linearization of a PM model (see [8, equation (25)]) which has a simpler quadratic term compared to our model. With the assumption that šæš‘‘=šæš‘ž, the solution of quadratic linearization becomes trivial. Our method is more general and is applicable to PMSM where šæš‘‘ā‰ šæš‘ž.

Krener [9] formulated an approximate feedback linearization technique based on Taylor series expansion; Kang and Krener [10, 11] extended the work of Poincare [12] on the lines of Krener to consider quadratic linearization of control affine systems. Quadratic linearization can be applied to a PMSM model [1] as the model has a predominant quadratic nonlinearity. But being an approximate technique, higher-order terms are introduced. Consideration of core loss involves fourth-order terms involving product of squares of both current and angular velocity. These together with the stray losses and unmodelled dynamics coupled with the third- and higher-order nonlinearities introduced due to quadratic linearization are best accounted for by tuning the transformations against an actual PM machine on the lines similar to those proposed by Levin and Narendra [13].

In this paper, input and state transformations are derived for a 4-dimensional PM machine model in order to linearize PMSM machine model. The PMSM model is quadratic linearized using the approximate technique. Tuning rules are derived for the linearizing transformations to account for higher-order terms, by back propagation of error between the outputs of quadratic linearized system with a normal form output.

The linearization technique is verified using SIMULINK model which is developed for interior permanent magnet (IPM) machine. The core loss which consists of higher-order terms is included in the SIMULINK model of PMSM, and tuning rules are also simulated. The closed loop response and open loop gain for the system before and after tuning are obtained.

The simulation results after linearization indicate a uniform closed loop response for different reference and load conditions, thus verifying the theory. The simulation results after tuning also verify the effectiveness of tuning.

To summarise the rest of the paper, in Section 2, background required for quadratic linearization is given. In Section 3, the model of PM synchronous motor is derived and is reduced to normal form. In Section 4, PMSM model is linearized. In Section 5, a SIMULINK model of PMSM is constructed and simulation is presented to show the effectiveness of the proposed method. In Section 6, tuning of the transformation coefficients is derived by including core loss in the PMSM model and also simulation results are given to verify the effectiveness of tuning. In Section 7, the paper is concluded.

2. Background

Consider a single input control affine system of the form [10]Ģ‡š‘„=š“š‘„+šµš‘¢+š‘“(2)(š‘„)+š‘“(3)(š‘„)+ā‹Æ+š‘“(š‘š)(š‘„)+ā‹Æ+š‘”(1)(š‘„)+š‘”(2)(š‘„)+ā‹Æ+š‘”(š‘šāˆ’1)(š‘„)+ā‹Æ,(1) where š“ and šµ are matrices in the controller normal formāŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£00ā‹…ā‹…1āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦š“=010ā‹Æ0001ā‹…ā‹…0ā‹…ā‹…ā‹…ā‹…ā‹…00ā‹…ā‹…ā‹…ā‹…1000ā‹…ā‹…0;šµ=,(2)š“ is an š‘›Ć—š‘› constant matrix and šµ is an š‘›Ć—1 constant matrix.

š‘„=[š‘„1š‘„2ā‹Æš‘„š‘›]š‘‡ and š‘¢ is a scalar input. š‘“(š‘š)(š‘„), š‘”(š‘šāˆ’1)(š‘„) are vector homogeneous polynomials of order š‘š and (š‘šāˆ’1), respectively, š‘š=2,3,ā€¦.

In order to linearize the system, a change of coordinates and feedback [10, 11] of the following form is considered:š‘¦=š‘„+šœ™(š‘„),(3)š‘¢=(1+š›½(š‘„))š‘£+š›¼(š‘„),(4) where š‘¦=[š‘¦1š‘¦2ā‹Æš‘¦š‘›]š‘‡, š‘£ is a scalar, and š›¼(š‘„), š›½(š‘„), and šœ™(š‘„) are polynomials given by šœ™(š‘„)=šœ™(2)(š‘„)+šœ™(3)(š‘„)+ā‹Æ+šœ™(š‘š)(š‘„)+ā‹Æ,š›¼(š‘„)=š›¼(2)(š‘„)+š›¼(3)(š‘„)+ā‹Æ+š›¼(š‘š)š›½(š‘„)+ā‹Æ,(š‘„)=š›½(1)(š‘„)+š›½(2)(š‘„)+ā‹Æ+š›½(š‘šāˆ’1)(š‘„)+ā‹Æ.(5)

Applying the transformations (3) and (4), (1) is reduced to Ģ‡š‘¦=š“š‘¦+šµš‘£(6) provided that the following equations, called generalized homological equations [10], are satisfied:āˆ’š“šœ™(š‘š)(š‘„)+šµš›¼(š‘š)(š‘„)+š‘“ī…ž(š‘š)(š‘„)+šœ•šœ™(š‘š)(š‘„)šœ•š‘„š“š‘„=0,šµš›½(š‘šāˆ’1)(š‘„)š‘£+šœ•šœ™(š‘š)(š‘„)šœ•š‘„šµš‘£+š‘”ī…ž(š‘šāˆ’1)(š‘„)š‘£=0;āˆ€š‘£,(7) where š‘“ī…ž(š‘š)(š‘„)=š‘“(š‘š)(š‘„) for š‘š=2 and š‘“ī…ž(š‘š)(š‘„) is expressed in terms of š‘“(š‘šāˆ’š‘–)(š‘„), š‘–=0,1,2ā€¦(š‘šāˆ’2) and šœ™(š‘šāˆ’š‘—)(š‘„), š‘—=1,2ā€¦(š‘šāˆ’2), š‘š>2. š‘”ī…ž(š‘šāˆ’1)(š‘„)=š‘”(1)(š‘„), š‘š=2 and š‘”ī…ž(š‘šāˆ’1)(š‘„) is expressed in terms of š‘”(š‘šāˆ’š‘–)(š‘„), š‘–=1,2ā€¦(š‘šāˆ’1) and šœ™(š‘šāˆ’š‘—)(š‘„), š‘—=1,2ā€¦(š‘šāˆ’2),š‘š>2.

Remark 1. Quadratic linearization involves specialization of the above result for š‘š=2. Consider Ģ‡š‘„=š“š‘„+šµš‘¢+š‘“(2)(š‘„)+š‘”(1)(š‘„)+š‘‚(3)(š‘„,š‘¢),(8) where š“ and šµ are in controller normal form [10] and š‘‚(3)(š‘„,š‘¢) corresponds to terms of order 3 or more. Equations (3) and (4) for š‘š=2 can be written as š‘¦=š‘„+šœ™(2)ī€·(š‘„),(9)š‘¢=1+š›½(1)(ī€øš‘„)š‘£+š›¼(2)(š‘„).(10) Applying (9) and (10), (8) can be written as Ģ‡š‘¦=š“š‘¦+šµš‘£+š‘‚ī…ž(3)(š‘¦,š‘£)(11) provided that the homological equations āˆ’š“šœ™(2)(š‘„)+šµš›¼(2)(š‘„)+š‘“(2)(š‘„)+šœ•šœ™(2)(š‘„)šœ•š‘„š“š‘„=0,šµš›½(1)(š‘„)š‘£+šœ•šœ™(2)(š‘„)šœ•š‘„šµš‘£+š‘”(1)(š‘„)š‘£=0;āˆ€š‘£(12) can be solved for šœ™(2)(š‘„),š›¼(2)(š‘„), and š›½(1)(š‘„).

3. Machine Model

The PM machine model can be derived [1, 14, 15] as Ģ‡š‘„=š“š‘„+šµš‘¢+š‘“(2)ī‚ƒš‘„(š‘„),(13)š‘„=1š‘„2š‘„3š‘„4ī‚„š‘‡=ī‚ƒšœƒšœ”š‘’š‘–š‘žš‘–š‘‘ī‚„š‘‡,ī‚ƒš‘¢š‘¢=1š‘¢2ī‚„š‘‡=ī‚ƒš‘£š‘žš‘£š‘‘ī‚„š‘‡,(14) where š‘£š‘ž,š‘£š‘‘,š‘–š‘ž, and š‘–š‘‘ represent the quadrature and direct axis voltages and currents, respectively, and šœƒ and šœ”š‘’ represent rotor position and angular velocity, respectively.āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£š“=0100001.5š‘šœ†š½00āˆ’šœ†š‘šæš‘žāˆ’š‘…šæš‘ž0000āˆ’š‘…šæš‘‘āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦,āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£001šµ=šæš‘ž00001šæš‘‘āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦,š‘“(2)āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£0ī€·šæ(š‘„)=1.5š‘š‘‘āˆ’šæš‘žī€øš‘–š‘‘š‘–š‘žš½āˆ’šæš‘‘š‘šœ”š‘’š‘–š‘‘šæš‘žšæš‘žš‘šœ”š‘’š‘–š‘žšæš‘‘āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦,(15)Ī» is the flux induced by the permanent magnet of the rotor in the stator phases. šæš‘‘ and šæš‘ž are the direct and quadrature inductances, respectively. š‘… is the stator resistance, š‘ is the number of pole pairs, and š½ is the system moment of inertia.

The model (13) can be reduced, in a standard way using linear change of coordinates and feedback [16], to the Brunovsky form [17] for two inputs as in (16) below (where š‘„,š‘¢,š“, and šµ are retained for simplicity of notation):Ģ‡š‘„=š“š‘„+šµš‘¢+š‘“(2)(š‘„),(16) where āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦š“=0100001000000000,šµ=00001001,š‘“(2)āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£0š‘˜(š‘„)=1š‘„3š‘„4š‘˜2š‘„2š‘„4š‘˜3š‘„2š‘„3āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦,š‘˜1=ī€·šæ1.5š‘š‘‘āˆ’šæš‘žī€øš‘Ž4š½š‘Ž1,š‘˜2=āˆ’šæš‘‘š‘š‘Ž1š‘Ž4šæš‘ž,š‘˜3=šæš‘žš¶12š‘Ž1šæš‘‘š‘Ž4,š‘Ž1=1.5š‘šœ†š½,š‘Ž4=āˆ’š‘…šæš‘‘,š¶1=1šæš‘ž.(17)

4. Linearization of Interior Permanent Magnet (IPM) Synchronous Motor

The quadratic linearization given in Section 2 is extended to the case of two inputs in a straightforward way and applied to a PM machine model in this section [14, 15].

Theorem 2. Given the 4-dimensional model of a PM synchronous motor (IPM model) of the form (16), the system can be linearized using the following transformations: š‘¦=š‘„+šœ™(2)ī€·š¼(š‘„),(18)š‘¢=2+š›½(1)(ī€øš‘„)š‘£+š›¼(2)(š‘„),(19) where š‘¢=[š‘¢1š‘¢2]š‘‡, š‘£=[š‘£1š‘£2]š‘‡, šœ™(2)āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£00š‘˜(š‘„)=1š‘„3š‘„40āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦š›¼,(20)(2)āŽ”āŽ¢āŽ¢āŽ£(š‘„)=āˆ’š‘˜2š‘„2š‘„4āˆ’š‘˜3š‘„2š‘„3āŽ¤āŽ„āŽ„āŽ¦š›½,(21)(1)ī‚µšµ(š‘„)=(āˆ’1)š‘‡šœ•šœ™(2)(š‘„)šµī‚¶āŽ”āŽ¢āŽ¢āŽ£š‘˜šœ•š‘„=āˆ’1š‘„4š‘˜1š‘„3āŽ¤āŽ„āŽ„āŽ¦00,(22) where š¼2 is the identity matrix of order 2.
The system then reduces to Ģ‡š‘¦=š“š‘¦+šµš‘£+š‘‚(3)(š‘¦,š‘£),(23) where š‘‚(3)(š‘¦,š‘£) represents third- and higher-order nonlinearities.

Proof. Applying transformations (18) and (19) (which correspond to a natural extension of (9) and (10) to two inputs) to (16), the homological equations to be considered can be written as āˆ’š“šœ™(2)(š‘„)+šµš›¼(2)(š‘„)+š‘“(2)(š‘„)+šœ•šœ™(2)(š‘„)šœ•š‘„š“š‘„=0,(24)šµš›½(1)(š‘„)+šœ•šœ™(2)(š‘„)šœ•š‘„šµ=0.(25) By choosing šœ™(2)(š‘„) and š›¼(2)(š‘„) as in (20) and (21), (24) reduces to āˆ’āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£0š‘˜1š‘„3š‘„400āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦+āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£00āˆ’š‘˜2š‘„2š‘„4āˆ’š‘˜3š‘„2š‘„3āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦+āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£0š‘˜1š‘„3š‘„4š‘˜2š‘„2š‘„4š‘˜3š‘„2š‘„3āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦+āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£0000000000š‘˜1š‘„4š‘˜1š‘„3āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£š‘„00002š‘„300āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦=0.(26)
Premultiplying (25) by šµš‘‡ and noting that šµš‘‡šµ=š¼2, (25) reduces to š›½(1)ī‚µšµ(š‘„)=(āˆ’1)š‘”šœ•šœ™(2)(š‘„)šµī‚¶.šœ•š‘„(27) Also, substitution of (27) satisfies (25). Hence, the homological equations are satisfied and quadratic linearization is achieved, hence the proof.

5. Experimental Simulation Results

Given the parameters š‘…=2.875ā€‰Ī©, šæš‘ž=9ā€‰mH, šæš‘‘=7ā€‰mH, šœ”š‘’=3500ā€‰rpm, š‘=4, š½=0.0008ā€‰kgm2 and šœ†=175ā€‰mWb turns, of an actual PM machine, Figure 1 shows the SIMULINK model of the PMSM which is constructedby using speed and torque blocks and control circuit. š‘£š‘ž and š‘£š‘‘ are taken as inputs to the motor. The model is reduced to normal form in a standard way by using the linear transformations. The PM model in Figure 1 is especially configured for the IPM where šæš‘‘ā‰ šæš‘ž. š‘1 and š‘2 blocks in Figure 2 include the above linear transformation and the nonlinear transformation (18) and (19) as well.

Prior to linearization, the open loop steady state gain of šœ”š‘’ versus š‘£š‘ž of the PMSM model is investigated and the results are given in Table 1. In this table, it is observed that the open loop steady state gain of šœ”š‘’ versus š‘£š‘ž (keeping š‘£š‘‘ constant) is not constant because of the system nonlinearity. To verify the linearity of the system after linearization, we investigated the variation of its gain of š‘¦2 (a scaled version of šœ”š‘’ as can be seen from (18) and (20)) with input š‘£1 (see Figure 2) and the results are given in Table 2. The table reveals that the gain of the system is nearly constant, thus verifying that by applying the homogeneous linearizing transformation, the PMSM model is made nearly linear for the given set of inputs.

š‘£ š‘ž šœ” š‘’ š‘˜ = š‘‘ šœ” š‘’ / š‘‘ š‘£ š‘ž


š‘£ 1 š‘¦ 2
Ɨ 1 0 āˆ’ 6
š‘˜ = š‘‘ š‘¦ 2 š‘‘ š‘£ 1
Ɨ 1 0 āˆ’ 6


We now proceed to show that the nearly constant gain of the linearized model results in a uniform closed loop response on a range of set point and load inputs with a fixed controller. This is in contrast to the case before linearization under the corresponding conditions.

Figures 3 and 4 show the time response of angular velocity šœ”š‘’ by closing the loop around PMSM model before linearization when š‘£š‘ž=5 and 30 units, respectively. It is observed that the dynamic response for š‘£š‘ž=5 is more oscillatory compared to the case of š‘£š‘ž=30 with a fixed controller of proportional gain = 50 and integral constant = 2. This is to be expected since the loop gain is higher in the former case with a higher static gain in the motor as can be seen from Table 1.

Figures 5 and 6 show the time response of š‘¦2 of the transformed PMSM system (Figure 2) in closed loop when š‘£1=5 and 30 units, respectively. It is observed that a uniform output response is obtained in the closed loop after linearization when the reference is varied. Since the static gain in Table 2 is nearly uniform, the loop gain is also nearly constant for the extreme points in the operating range, thus resulting in the uniform dynamic responses as in Figures 5 and 6.

Figures 7 and 8 show the time response of angular velocity šœ”š‘’ of the PMSM model before linearization in closed loop when š‘£š‘ž=15 units with the load torque š‘‡ varied as š‘‡=0 and 1 unit, respectively. Note that the dynamic responses in Figures 7 and 8 are different, with the response of Figure 8 corresponding to an inverse response (i.e., speed reduces first before increasing). The latter corresponds to a response of a nonminimum phase system.

Figures 9 and 10 show the time response of š‘¦2 of the transformed PMSM system in closed loop when š‘£1=15 units and š‘‡=0 and 1 unit, respectively. It is observed that a uniform output response is obtained in the closed loop after linearization when load torque is varied. Hence, it is verified that a linearized system gives a uniform closed loop response for the different reference and load conditions used for testing as above.

6. Tuning

6.1. Core Loss

The core loss or iron loss, caused by the permanent magnet (PM) flux and armature reaction flux, is a significant component in the total loss of a PMSM, and, thus, it can have a considerable effect on the PMSM modeling and performance prediction.

The net core loss š‘ƒš‘™š‘ [2] for the machine is computed as follows: š‘ƒš‘™š‘=īƒ¬1.5šœ”š‘Ÿ2ī€·šæš‘žš¼š‘žī€ø2š‘…š‘+1.5šœ”š‘Ÿ2ī€·šœ†š‘Žš‘“+šæš‘‘š¼š‘‘ī€ø2š‘…š‘īƒ­,(28) where š‘…š‘ represents core loss resistance, šœ†š‘Žš‘“ represents magnet flux linkage and šœ”š‘Ÿ denotes rotor electrical speed.

The mechanical torque equation including core losses is given byš‘‡š‘’=ī‚€2š‘ƒī‚š½š‘‘šœ”š‘’š‘‘š‘”+š‘‡š‘™+š‘‡š‘™š‘,š‘‡(29)š‘™š‘=š‘ƒš‘™š‘šœ”š‘Ÿ.(30)

6.2. Tuning Formula

To account for the core loss (28), unmodelled dynamics, and higher-order nonlinearities, tuning of the quadratic linearization transformations š‘1 and š‘2 can be done [13]. Figure 11 shows the block diagram for tuning.

Error (šø) can be calculated as ī€·šœ€šø=š‘‡šœ€ī€ø=ī€ŗ(š‘¦āˆ’Ģ‚š‘¦)š‘‡ī€»(š‘¦āˆ’Ģ‚š‘¦)1/2,(31) where šœ€=(šœ€1šœ€2šœ€3šœ€4)š‘‡.

The error can be written as ī€·šœ€šø=21+šœ€22+šœ€23+šœ€24ī€ø1/2.(32) Since šœ™(2)(š‘„) and š›½(1)(š‘„) are both functions of š‘˜1, we shall redefine šœ™(2)āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£00š‘˜(š‘„)=1š‘„3š‘„40āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦,š›½(1)āŽ”āŽ¢āŽ¢āŽ£š‘˜(š‘„)=āˆ’ī…ž1š‘„4š‘˜ī…ž1š‘„3āŽ¤āŽ„āŽ„āŽ¦00(33) so that šœ™(2)(š‘„) and š›½(1)(š‘„) can be independently tuned by tuning š‘˜1 and š‘˜ī…ž1, respectively. š›¼(2)(š‘„) is not varied.

6.2.1. Updation of š‘2 Transformation Coefficients

Tuning of š‘2 transformation implies the tuning of šœ™(2)(š‘„). As šœ™(2)(š‘„) is a function of only š‘˜1š‘„3š‘„4, the coefficient š‘˜1 has to be updated based on the error between the outputs of quadratic linearized system and normal form. The updation law is derived as follows: Ī”š‘˜1=šœ•šøšœ•š‘˜1=šœ•šøšœ•š‘¦šœ•š‘¦šœ•š‘˜1.(34) From (31), it is seen that šœ•šø/šœ•š‘¦š‘–=šœ€š‘–/šø, š‘–=1,2,3,4.

Hence, šœ•šø/šœ•š‘¦=[šœ€1/šøšœ€2/šøšœ€3/šøšœ€4/šø],āˆ“Ī”š‘˜1=ī‚øšœ€1šøšœ€2šøšœ€3šøšœ€4šøī‚¹āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£00š‘„3š‘„40āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦=šœ€3šøš‘„3š‘„4.(35) Updation of šœ™(2)(š‘„) is done by using the formula š‘˜1(š‘š)=š‘˜1(š‘šāˆ’1)āˆ’šœŒ1Ī”š‘˜1(š‘š);0<šœŒ1<1,(36) where š‘š corresponds to the updating step and šœŒ1 corresponds to the accelerating factor.

6.2.2. Updation of š‘1 Transformation Coefficients

Tuning of š‘1 transformation is achieved by tuning of š›½(1)(š‘„). As š›½(1)(š‘„) is a function of š‘˜ī…ž1š‘„3 and š‘˜ī…ž1š‘„4, the coefficient š‘˜ī…ž1 has to be updated based on the error between the outputs of quadratic linearized system and normal form. The updation law is derived as follows: Ī”š‘˜ī…ž1=šœ•šøšœ•š‘˜ī…ž1=šœ•šøšœ•š‘¦šœ•š‘¦šœ•š‘„šœ•š‘„šœ•š‘¢1šœ•š‘¢1šœ•š‘˜ī…ž1,(37) where šœ•šø=ī‚øšœ€šœ•š‘¦1šøšœ€2šøšœ€3šøšœ€4šøī‚¹,šœ•š‘¦=āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£ī€·šœ•š‘„10000100001+š‘˜ī…ž1š‘„4ī€øš‘˜ī…ž1š‘„3āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦,0001šœ•š‘„šœ•š‘¢1=āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£0āˆ’1š‘˜2š‘„40āˆ’1š‘˜2š‘„2āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦,(38) assuming that the steady state of the SIMULINK model is reached within the tuning period: šœ•š‘¢1šœ•š‘˜ī…ž1=āˆ’š‘£1š‘„4āˆ’š‘£2š‘„3.(39) Thus, Ī”š‘˜ī…ž1=((š‘£1š‘„4+š‘£2š‘„3)/šø)(šœ€2/š‘˜2š‘„4+(šœ€3š‘˜ī…ž1š‘„3+šœ€4)/š‘˜2š‘„2).

Updation of š›½(1)(š‘„) is done by using the formula š‘˜ī…ž1(š‘š)=š‘˜ī…ž1(š‘šāˆ’1)āˆ’šœŒ2Ī”š‘˜ī…ž1(š‘š),0<šœŒ2<1,(40) where š‘š corresponds to the updating step and šœŒ2 corresponds to the accelerating factor.

6.3. Tuning Simulation Results

PMSM model including core loss is given in Figure 12 where the block due to loss torque (30) is also included. The tuning of the transformations is done using the updations for šœ™(2)(š‘„) and š›½(1)(š‘„) as per formulae (36) and (40) as given in Figure 13. The simulation diagram for tuning is done by using memory blocks to store the updated values of šœ™(2)(š‘„) and š›½(1)(š‘„). It is seen that the error after tuning is reduced to 0.01. Table 3 shows the variation of š‘¦2 versus š‘£1 of the PMSM model after linearization including core loss prior to tuning of the coefficients. It is observed that the open loop steady state gain is not constant due to the effect of core loss. Table 4 shows the variation of š‘¦2 versus š‘£1 for the linearized system after incorporating tuning of the transformation coefficients. The table reveals that the gain of the system is nearly constant thus verifying the effectiveness of tuning. Figures 14 and 15 show the closed loop time response of š‘¦2 of the linearized system including core loss before tuning when š‘£1=5 units and 30 units, respectively; š‘˜š‘=50; š‘˜š‘–=2. Figures 16 and 17 show the closed loop time response of š‘¦2 of the linearized system including core loss after tuning when š‘£1=5 units and 30 units, respectively; š‘˜š‘=50; š‘˜š‘–=2. It is observed that the dynamic response before tuning for š‘£1=5 units is more oscillatory compared to the case of š‘£1=30 units. This is to be expected since the loop gain is higher in the former case with a higher static gain in the plant or motor as can be seen from Table 3. Since the static gain in Table 4 is nearly uniform, the loop gain is also nearly constant for the extreme points in the operating range, thus resulting in the uniform dynamic responses in Figures 16 and 17. Thus, the effectiveness of tuning the linearizing transformations to overcome the unmodelled core loss is verified.

š‘£ 1 š‘¦ 2
Ɨ 1 0 āˆ’ 6
š‘˜ = š‘‘ š‘¦ 2 / š‘‘ š‘£ 1
Ɨ 1 0 āˆ’ 6


š‘£ 1 š‘¦ 2
Ɨ 1 0 āˆ’ 6
š‘˜ = š‘‘ š‘¦ 2 / š‘‘ š‘£ 1
Ɨ 1 0 āˆ’ 6


7. Conclusion

As the PMSM is inherently nonlinear, to design a controller that can provide a predictable uniform performance of the drive under varying operating conditions, it is necessary to linearize the PMSM. In this paper, the dynamic model of a PM synchronous motor involving quadratic nonlinearity is linearized. The technique used is based on the control input extension of Poincareā€™s work due to Kang and Krener which is in line with the approximate linearization technique of Krener [9, 18]. The existing techniques of exact feedback linearization introduce singularities in the system, which may cause difficulties in the implementation of the closed loop control. In the proposed method, the problem of singularities does not arise.

The PMSM machine model, together with the state and input transformations, are simulated using SIMULINK. The simulation results show that the quadratic linearizing transformations effectively linearize the system thus, supporting the theory. The simulation verifies that a uniform response under a fixed controller is obtained for the linearized system for variations of reference speed and load conditions, in contrast to the case before linearization.

Further, to account for the core loss, unmodelled dynamics and third- and higher-order nonlinearities, tuning of the transformation parameters is proposed by comparing the output of the linearized system with a normal form output. After tuning, it is shown that the linearized system shows improved linearity in terms of static gain compared to the condition before tuning. Closed loop system response for the linearized system is also shown to be uniform due to the effect of tuning of linearization transformations, under varying reference inputs.

The proposed linearization method can be extended to induction motor and wound synchronous motor models as well [19].


  1. B. K. Bose, Modern Power Electronics and AC Drives, Pearson Education, 2002.
  2. R. Monajemy, Control strategies and parameter compensation of permanent magnet synchronous motor drives, Ph.D. thesis, Virginia Polytechnic Institute and State University, Blachsberg, Va, USA, 2000.
  3. J. K. Seok, J. K. Lee, and D. C. Lee, ā€œSensorless speed control of nonsalient permanent-magnet synchronous motor using rotor-position-tracking PI controller,ā€ IEEE Transactions on Industrial Electronics, vol. 53, no. 2, pp. 399ā€“405, 2006. View at: Publisher Site | Google Scholar
  4. P. Pillay and R. Krishnan, ā€œModeling, simulation, and analysis of permanent-magnet motor drives. I. The permanent-magnet synchronous motor drive,ā€ IEEE Transactions on Industry Applications, vol. 25, no. 2, pp. 265ā€“273, 1989. View at: Google Scholar
  5. H. Zhu, X. Xiao, and Y. Li, ā€œPI type dynamic decoupling control scheme for PMSM high speed operation,ā€ in Proceedings of the 25th Annual IEEE Applied Power Electronics Conference and Exposition (APEC '10), pp. 1736ā€“1739, Palm Springs, Calif, USA, February 2010. View at: Publisher Site | Google Scholar
  6. M. Bodson and J. Chiasson, ā€œDifferential-geometric methods for control of electric motors,ā€ International Journal of Robust and Nonlinear Control, vol. 8, no. 11, pp. 923ā€“954, 1998. View at: Google Scholar
  7. J. Chiasson, ā€œNonlinear controllers for an induction motor,ā€ Control Engineering Practice, vol. 4, no. 7, pp. 977ā€“990, 1996. View at: Publisher Site | Google Scholar
  8. G. Zhu, A. Kaddouri, L. A. Dessaint, and O. Akhrif, ā€œA nonlinear state observer for the sensorless control of a permanent-magnet AC machine,ā€ IEEE Transactions on Industrial Electronics, vol. 48, no. 6, pp. 1098ā€“1108, 2001. View at: Publisher Site | Google Scholar
  9. A. J. Krener, ā€œApproximate linearization by state feedback and coordinate change,ā€ Systems and Control Letters, vol. 5, no. 3, pp. 181ā€“185, 1984. View at: Google Scholar
  10. W. Kang and A. J. Krener, ā€œExtended quadratic controller normal form and dynamic state feedback linearization of nonlinear systems,ā€ SIAM Journal on Control and Optimization, vol. 30, no. 6, pp. 1319ā€“1337, 1992. View at: Google Scholar
  11. A. J. Krener and W. Kang, ā€œExtended normal forms of quadratic systems,ā€ in Proceedings of the 29th IEEE Conference on Decision and Control, pp. 2091ā€“2096, IEEE, New York, NY, USA, December 1990. View at: Google Scholar
  12. V. I. Arnold, Geometric Methods in the Theory of Ordinary Differential Equations, Springer, New York, NY, USA, 1983.
  13. A. U. Levin and K. S. Narendra, ā€œControl of nonlinear dynamical systems using neural networks. Controllability and stabilization,ā€ IEEE Transactions on Neural Networks, vol. 4, no. 2, pp. 192ā€“206, 1993. View at: Publisher Site | Google Scholar
  14. A. K. Parvathy, V. Kamaraj, and R. Devanathan, ā€œA new linearisation technique for permanent magnet synchronous motor model,ā€ in Proceedings of the Joint International Conference on Power System Technology and IEEE Power India Conference (POWERCON '08), pp. 1ā€“5, New Delhi, India, October 2008. View at: Publisher Site | Google Scholar
  15. A. K. Parvathy, R. Devanathan, and V. Kamaraj, ā€œApplication of quadratic linearization to control of Permanent Magnet synchronous motor,ā€ in Proceedings of the 1st International Conference on Electrical Energy Systems (ICEES '11), pp. 158ā€“163, 2011. View at: Publisher Site | Google Scholar
  16. B. C. Kuo, Automatic Control Systems, Prentice-Hall, New Delhi, India, 2001.
  17. P. Brunovsky, ā€œA classification of linear controllable systems,ā€ Kybernetika, vol. 6, no. 3, pp. 173ā€“188, 1970. View at: Google Scholar
  18. G. S. Cardoso and L. Schnitman, ā€œAnalysis of exact linearization and aproximate feedback linearization techniques,ā€ Mathematical Problems in Engineering, vol. 2011, Article ID 205939, 17 pages, 2011. View at: Publisher Site | Google Scholar
  19. A. K. Parvathy, V. Kamaraj, and R. Devanathan, ā€œComplete quadratic linearisation of machine models,ā€ in Proceedings of the IEEE International Conference on Control Applications, pp. 1130ā€“1133, Singapore, October 2007. View at: Publisher Site | Google Scholar

Copyright © 2011 Parvathy Ayalur Krishnamoorthy et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.