Abstract
In the exact linearization of involutive nonlinear system models, the issue of singularity needs to be addressed in practical applications. The approximate linearization technique due to Krener, based on Taylor series expansion, apart from being applicable to noninvolutive systems, allows the singularity issue to be circumvented. But approximate linearization, while removing terms up to certain order, also introduces terms of higher order than those removed into the system. To overcome this problem, in the case of quadratic linearization, a new concept called “generalized quadratic linearization” is introduced in this paper, which seeks to remove quadratic terms without introducing third- and higher-order terms into the system. Also, solution of generalized quadratic linearization of a class of control affine systems is derived. Two machine models are shown to belong to this class and are reduced to only linear terms through coordinate and state feedback. The result is applicable to other machine models as well.
1. Introduction
Control of nonlinear systems is gaining increasing attention in recent years due to its technical importance and its impact on various applications as well. When nonlinearities of systems become significant, linear control techniques generally fail to produce the desired results. However, feedback linearization of nonlinear systems allows linear control methods to be applied effectively to the resultant linearized system.
In the application of exact linearization technique [1–3] to involutive systems, an issue of singularity needs to be addressed in practical situations. The singularity is caused when a function of the system state tends towards zero causing the state feedback to fail. Approximate linearization technique due to Krener [4] based on a Taylor series expansion, applicable to non-involutive systems, has the merit of allowing the problem of singularity to be overcome.
Poincare [5] derived what are known as homological equations for approximate linearization of autonomous differential systems. Kang and Krener [6–8] extended the result for the approximate linearization of control affine systems and derived what are termed generalized homological equations. Tall [9, 10] extended Kang and Krener's normal forms to multi-input systems. Devanathan [11, 12] developed necessary and sufficient conditions for quadratic linearization of control affine systems with a single input.
Using the approximate linearization technique, the problem of singularity of state feedback can be circumvented by appropriate definition of state feedback. But the approximate linearization, while removing nonlinearities up to a certain higher order, introduces nonlinearities into the system of higher order than those removed. Introduction of higher-order nonlinearities into a system which originally may not have such nonlinearities is a matter of concern.
In the case of machine models, for example, the second-order dynamic models as given in [13] seem adequate to describe the normal behaviour of the machine. Applying quadratic linearization based on approximate linearization technique to this model introduces third- and higher-order terms into the model even though the machine does not possess nonlinearities of this order during normal operation.
In this paper, a new concept called “generalized quadratic linearization” which seeks to remove the second-order nonlinearity in the model without introducing third-and higher-order nonlinearities in the process is introduced. This is in contrast to the existing quadratic linearization techniques due to Kang and Krener [6–8] and Devanathan [11, 12] which introduce third- and-higher order terms in the process of removing second-order terms. The generalized quadratic linearization, being a stronger condition than quadratic linearization, imposes additional constraints on the quadratic polynomials of the system than those imposed due to quadratic linearization alone. Sufficient conditions on the quadratic polynomials together with coordinate and state feedback transformations are derived for a class of systems for which generalized quadratic linearization is applicable. In particular, the results are shown to apply to the induction motor and permanent magnet synchronous motor models. The results can be extended to other machine models as well. Application of generalized quadratic linearization to machine models also helps to avoid the issue of singularity which is a drawback attributed to the existing exact linearization of machine models [2, 3]. This is the main contribution of the paper.
To summarise the rest of the paper, in Section 2, homological equations for arbitrary order linearization are introduced. Section 3 helps to simplify the homological equations for quadratic linearization, thus paving the way for the central result on generalized quadratic linearization presented in Section 4. In Section 5, generalized quadratic linearization is applied to two machine models including a real-world numerical example. In Section 6, the paper is concluded.
2. Background
Consider a multiple-input control affine system of the form [7, 14, 15] where are in Brunovsky normal form. That is, let be a sequence of integers such that such that where matrix is partitioned into blocks while matrix is partitioned into blocks with . Each block or is of the form is of order , is of order , and is matrix. and are the vector input such that . and are vector polynomials of order and , respectively, .
In order to linearize the system, change of coordinate and state feedback [6] of the following form is considered: where and is the identity matrix of order . Applying the transformations (4), (1) is reduced to provided the following equations, called homological equations [6], are satisfied for : where for and is expressed in terms of , and . ; , and is expressed in terms of , and .
3. Quadratic Linearization
Consider the specialized case of system (1) containing only quadratic terms in as in Quadratic linearization of (9) involves specialization of the result discussed in Section 2 for and simplifying it by applying . That is, applying the following coordinate transformation and feedback,
(9) is reduced to where corresponds to terms of order 3 or higher, provided the following equations are satisfied:
Remark 3.1. Applying quadratic linearization to system (9) introduces third- and higher-order terms into the system (see (11)) even though originally the system does not possess nonlinearities of this order. To avoid this situation, we introduce a new concept of generalized quadratic linearization which seeks to remove the second-order nonlinearity in the system without introducing third- and higher-order nonlinearities in the process. Sufficient conditions for the generalized quadratic linearization for a class of systems are discussed in the next section.
We next simplify the homological equations (12) and (13) by representing them in partitioned form, using the Brunovsky normal form introduced in Section 2.
Let where putting , where where
Using (2), (14)–(17), and (21), the th block row of (12) can be written as The second homological equation, namely, (13), can be simplified as follows. Since the effect of control input in (9) is dependent on the matrix only, without loss of generality, we can then assume that , because the terms in matrix are nonzero only for the first inputs. So, we can write the revised matrix as
We now put (13) in partitioned form. Using (2) and (21)–(23), the th block row of (13) can be written as
Since (26) is to be an identity in , we can drop and rewrite (26) as Since , are arbitrary, it is clear that a necessary condition for (13) to be satisfied is
4. Generalized Quadratic Linearization
Theorem 4.1. Consider the system where and . are in Brunovsky normal form as given in (2), (3), and (25). is given by (16), where is of order and is given by the r-tuple . Then the transformation (4), where , is given by (14), where given by (17), where reduces the system (29) to
Proof. We first show that the quadratic linearization of (29) is achieved using (10), assuming (30)–(33). That is, it needs to be verified that (12) and (13) are satisfied. To verify that (12) is satisfied, it is sufficient to verify that (24) is satisfied. Substituting (30)–(32) into (24), we need to satisfy the following:
Note that is a function of , only. Hence, . Equation (35) then reduces to the following after simplification:
Clearly (36) is an identity. Hence (24) is satisfied for . For the case , it is easy to verify that (24) is satisfied as is a function of , only.
Equation (28) is a necessary condition for (13) to hold, and it follows easily from (31), where . To satisfy (13), note that, from (25), it follows that
Premultiply (13) by and rearrange to yield
which is (33) for . Substitution of (38) satisfies (13). Hence, the result holds for .
As a result of quadratic linearization (i.e., satisfying (12) and (13)), higher-order terms are left in the system as per (11) even though originally (29) did not have third- and higher-order terms. These are now to be removed to effect generalized quadratic linearization. The higher-order terms are functions of , and and are given for the case , through algebraic computations as
For ,
Partitioning similar to in (16) and considering the th block row,
Using (31) and (32) and that being a function of , only as above, (42) can be reduced to
Hence, . Homological equation (7), for , is satisfied if we put and .
For ,
As and , it follows from (44) that . Hence, one can have and to satisfy (7) for . Proceeding this way, it can be proved that .
Noting that from (40), for , and (8) can be written as
from which using the fact that using (38) is derived. Substitution of satisfies (8) for .
Proceeding this way, one can show that , and we can choose , to satisfy (7). as given in (33) is easily derived so as to satisfy (8), where is given in (40), using the same technique as for successively as above, hence the result.
Remark 4.2. The result of Theorem 4.1 provides a solution to the generalized quadratic linearization of a class of systems defined by (29). Equation (29) may appear to be a strong condition. But this is needed because of the requirement that the third- and all higher-order nonlinearities, which are usually introduced into the system during the course of quadratic linearization, identically vanish. As shown in Section 5, system (29) can represent many electric machine models [13] of which induction motor and permanent magnet synchronous motor are examples. Also, temperature and flow variables in heat transfer (composition and flow variables in mass transfer) processes, for example, combine to produce quadratic nonlinearity [16] of the form (9) which is a general form of (29). Further, even if a quadratic model (9) does not correspond exactly to model (29), one can always tune the coordinate and state feedback on the lines described by the authors in [17] to bring the system to a linear normal form. It can also be pointed out that quadratic system of the form can be brought to the form through state feedback as shown in [6]. Hence the forms (9) and (29), which are a special case of (9), are more general than they appear to be.
Remark 4.3. A specialized result of Theorem 4.1 appears in an earlier paper by the authors [18]. In [18], the concept of generalized quadratic linearization is applied to a class of control affine systems with a specific configuration of Brunovsky form for multiple inputs. A PMSM model is shown to belong to this class and the generalized quadratic linearization of PMSM is derived. In Theorem 4.1, however, we have considered a more general class of systems with no assumption on the configuration of the Brunovsky form for multiple inputs. Hence, the theoretical result of [18] is a special case of the result of Theorem 4.1.
5. Generalized Quadratic Linearization of Machine Models
5.1. Squirrel Cage Induction Motor
The squirrel cage induction motor model can be derived as below [13]: where , where , where , and represent the quadrature and direct axis voltages and currents, respectively, and represent direct and quadrature axis fluxes, respectively. represent the estimated direct and quadrature flux, respectively, and represents the angular velocity of the rotor. and represent the reference values of quadrature and direct axis currents, respectively. The superscript refers to variables in the stationary frame: where , and represent the magnetizing, rotor, and stator inductances, respectively. and represent the rotor and stator resistances, respectively. represents the system moment of inertia, and represents the number of poles. and . The model (46) can be reduced to normal form in a standard way using linear transformation of state variables [19] and linear state feedback as follows, where ,, , and are retained for simplicity of notation: where and given by (2), where and , are where , , , , .
Remark 5.1. can be partitioned using (16) as where As is a function of , and only and hence a function of , one can write . The systems given by (48) and (49) are then in line with conditions of Theorem 4.1 as given in (30). Hence, Theorem 4.1 can be applied to system (48) to obtain the solution of generalized quadratic linearization which is given in Corollary 5.2.
Corollary 5.2. The system (48) can be quadratic linearized in the generalized sense using the transformation where where , and are completely arbitrary, and . The system reduces to
Proof. The result follows directly by applying the result of Theorem 4.1 to (48) where , and , and .
5.2. Permanent Magnet Synchronous Motor
The PM machine model [13] can be expressed as and where , and represent the quadrature and direct axes voltages and currents, respectively, and is stator angular velocity: is the field current equivalent to the permanent magnet. , and are the direct, quadrature, and magnetizing inductances, respectively. is the stator resistance. is the number of poles, and is the system moment of inertia. The model (55) can be reduced to the following normal form in a standard way using linear transformation of state variables [19] and linear state feedback, where , , , , and are retained for simplicity of notation: where and given by (2), where and , are with being constants.
Remark 5.3. can be partitioned using (16) as where Note that is a function of and hence . The system given by (57) is in line with conditions of Theorem 4.1 as given in (30). Hence, Theorem 4.1 can be applied to system (57) to obtain the solution of generalized quadratic linearization which is given in Corollary 5.4.
Corollary 5.4. The system (57) can be quadratic linearized in the generalized sense using the transformation where where is completely arbitrary and . The system then reduces to
Proof. The result follows directly by applying the result of Theorem 4.1 to (57) where , and .
5.3. Real-World Example of PMSM
The parameters of an actual machine are obtained from [20] as ohms, mH, 5.7 mH, 3500 rpm, , kgm2, and 123 mWb turns. The mathematical model of the machine in the normal form is given by where and are as given in (57) and where . By applying Corollary 5.4 and choosing
It is seen by verification that the system (64) becomes through transformation (61), given , and as above. The above system is completely linearized with no higher-order terms left behind.
5.4. Simulation Results
Verification of quadratic linearization of PMSM through simulations is given in an earlier paper by the authors [18]. The simulation results show that the closed-loop dynamic response of PMSM, before linearization, exhibits a behaviour that corresponds to varied responses for variations in reference and load conditions. However, it is verified that the quadratic linearized PMSM system gives a uniform closed-loop response for different reference and load conditions, as is expected of a linear system.
6. Conclusion
In this paper, linearization of machine models, which are predominantly quadratic in nature, is considered. Since the machine models only exhibit higher-order nonlinearities such as core loss, stray loss, and saturation, under extreme conditions, quadratic models can be used to represent the machines during normal operation.
The existing exact feedback linearization techniques introduce singularities in the system. This may result in the the faliure of linearization and the corresponding control technique, during the course of operation of the machine. Application of Poincare's approximate linearization technique due to Kang and Krener can remove the quadratic nonlinearity through a homogeneous transformation. This technique does not suffer from singularity issues but, being an approximate method, introduces higher-order terms in the process of quadratic linearization. Hence, in this paper, the concept of generalized quadratic linearization technique is introduced, wherein the higher-order terms introduced during the process of quadratic linearization are removed even as the quadratic term is removed.
A solution to the problem of generalized quadratic linearization is given for a class of control affine systems. Induction motor and PM motor models involving quadratic nonlinearity are considered and are quadratically linearised in the generalized sense. The proposed method can also be extended to wound rotor and synchronous machine models as well.
The generalized quadratic linearization technique proposed in this paper can also be considered as the approximate feedback linearization equivalent of exact linearization of a class of systems of the form (29). Since the solution of the generalized quadratic linearization involves removing third- and higher-order nonlinearities, the result of the paper can also been seen as a contribution towards a special case of the open problem on arbitrary order linearization [21] of control affine systems.
Finally, as per the result of Theorem 4.1, whose convergence needs to be established. Also, Theorem 4.1 of Section 4 deals only with sufficient conditions for generalized quadratic linearization. The question arises as to what are the necessary conditions for (9) to be generalized quadratically linearizable in general: this is our future work.