The invention relates generally to rotating electrical machines. More particularly, the invention relates to method for controlling an electrical machine comprising two or more multiphase stator windings. Furthermore, the invention relates to an apparatus and a computer program for the purpose of controlling an electrical machine comprising two or more multiphase stator windings.
In many applications, high level of availability of an electrical drive system is required. For example, in electrical ship propulsion the availability is critical from the viewpoint of safety. The availability can be increased by adding redundancy to the electrical drive system. Due to economic reasons, redundancy by multiplying the number of complete electrical drive systems is only rarely possible. However, the redundancy can be achieved by providing an electrical machine with two or more galvanically isolated multiphase stator windings each of which being supplied with its own multiphase power stage. The electrical machine may comprise, for example, two star-connected three-phase stator windings which are shifted relative to each other by e.g. 30 electrical degrees.
Accurate control of an electrical machine is typically based on a model which models the behavior of currents and voltages of the electrical machine and sometimes also the produced torque. In conjunction with a synchronous electrical machine which may have a salient pole rotor, the currents, voltages, and flux linkages are preferably expressed in suitable rotation-converted forms in a coordinate system bound to the rotor in order to avoid position dependency of inductance parameters of the model. The rotation-converted stator currents can be controlled on the basis of differences between the rotation-converted stator currents and their target values. The target values of the rotation-converted stator currents can be formed on the basis of e.g. the desired torque. The rotation-converted stator currents and voltages are typically expressed in the d-q coordinate system whose coordinate axes are along the direct and quadrature axes of the rotor. The inherent advantage of the d-q coordinate system is that the d-component of the stator currents does not generate flux-linkage on the q-direction, and correspondingly the q-component of the stator currents does not generate flux-linkage on the d-direction. This de-coupling between the d- and q-directions significantly facilitates the control of the rotation-converted stator currents because the d- and q-components of the stator currents can be regulated with e.g. separate regulators that can be, for example, proportional-integrating “PI” regulators.
In a case of an electrical machine which comprises two or more multiphase stator windings, the situation is more complicated. The two or more multiphase stator windings have mutual magnetic couplings. Hence, for example the d-directional stator flux-linkage of one of the multiphase stator windings is dependent not only on the d-component of the stator currents of this multiphase stator winding but also on the d-components of the stator currents of the other multiphase stator windings. In a simple control principle, the above-mentioned mutual magnetic couplings are neglected and the two or more multiphase stator windings are controlled separately from each other. However, the neglecting of the above-mentioned mutual magnetic couplings weakens the accuracy of the control. On the other hand, the control gets significantly more complicated if the mutual magnetic couplings are taken into account, because this requires cooperation between regulators related to different multiphase stator windings.
The following presents a simplified summary in order to provide a basic understanding of some aspects of various invention embodiments. The summary is not an extensive overview of the invention. It is neither intended to identify key or critical elements of the invention nor to delineate the scope of the invention. The following summary merely presents some concepts of the invention in a simplified form as a prelude to a more detailed description of exemplifying embodiments of the invention.
In accordance with the first aspect of the invention, there is provided a new method for controlling an electrical machine comprising two or more multiphase stator windings. The method according to the invention comprises:
The above-mentioned form of the model in which each of the stator flux-linkages of the model is dependent on only one of the rotation-converted stator currents can be obtained, for example, by diagonalizing the inductance matrix of the traditional d-q model of the electrical machine.
The use of the form where each stator flux-linkage is dependent on only one rotation-converted stator current significantly facilitates the control of the rotation-converted stator currents because the rotation-converted stator currents can be regulated with e.g. separate proportional-integrating “PI” regulators.
In accordance with the second aspect of the invention, there is provided a new apparatus for the purpose of controlling an electrical machine comprising two or more multiphase stator windings. The apparatus according to the invention comprises a processing system configured to:
The apparatus may further comprise, for example but not necessarily, power stages for providing supply voltages to the two or more multiphase stator windings on the basis of the voltage control signals.
In accordance with the third aspect of the invention, there is provided a new computer program for the purpose of controlling an electrical machine comprising two or more multiphase stator windings. The computer program according to the invention comprises computer executable instructions for controlling a programmable processor to:
In accordance with the fourth aspect of the invention, there is provided a new computer program product. The computer program product comprises a nonvolatile computer readable medium, e.g. a compact disc “CD”, encoded with a computer program according to the invention.
A number of exemplifying embodiments of the invention are described in accompanied dependent claims.
Various exemplifying embodiments of the invention both as to constructions and to methods of operation, together with additional objects and advantages thereof, will be best understood from the following description of specific exemplifying embodiments when read in connection with the accompanying drawings.
The verbs “to comprise” and “to include” are used in this document as open limitations that neither exclude nor require the existence of unrecited features. The features recited in depending claims are mutually freely combinable unless otherwise explicitly stated.
The exemplifying embodiments of the invention and their advantages are explained in greater detail below in the sense of examples and with target to the accompanying drawings, in which:
Without limiting generality, we assume for the sake of simplicity that the rotor of the electrical machine is a permanent magnet rotor which does not comprise any excitation and damper windings. In this case, the electrical machine can be modeled with the following voltage equations expressed in the d-q coordinate system:
where ud1 and uq1 are rotation-converted stator voltages of the first three-phase stator windings expressed in the d-q coordinate system, ud2 and uq2 are the rotation-converted stator voltages of the second three-phase stator windings expressed in the d-q coordinate system, id1 and iq1 are rotation-converted stator currents of the first three-phase stator windings expressed in the d-q coordinate system, id2 and iq2 are the rotation-converted stator currents of the second three-phase stator windings expressed in the d-q coordinate system, and Rs is the stator resistance.
Ψd1 and Ψq1 are the d- and q-directional stator flux linkages of the first three-phase stator windings, respectively, and Ψd2 and Ψq2 are the d- and q-directional stator flux linkages of the second three-phase stator windings, respectively. These flux linkages can be modeled with the following equations:
Ψd1=Ldid1+Mdid2+ΨPMd,
Ψq1=Lqiq1+Mqiq2,
Ψd2=Ldid2+Mdid1+ΨPMd, and
Ψq2=Lqiq2+Mqiq1, (2)
where Ld is the direct-axis inductance, Md is the direct-axis mutual inductance between the first and second three-phase stator windings, Lq is the quadrature-axis inductance, Mq is the quadrature-axis mutual inductance between the first and second three-phase stator windings, and ΨPMd is the d-directional stator flux linkage created by the permanent magnet rotor.
The model of the electrical machine defined by equations (1) and (2) can be presented in the following matrix form:
where:
In the above-presented equation (3), vectors are denoted with lower case bold letter and matrices with upper case bold letters. This notation will be used also in other parts of this document.
The vector of the rotation-converted stator voltages udq can be obtained from the vector of the stator phase voltages uabc=[ua1, ub1, uc1, ua2, ub2, uc2]T with the following conversion rule:
udq=T1(θ)uabc, (4)
where θ is the angle expressing the rotational position of the rotor with respect to the stator as shown in
where the sub-matrices Tp are given by:
where C1 is a constant scaling coefficient.
The vector of the rotation-converted stator currents idq can be obtained from the vector of the stator phase currents iabc=[ia1, ib1, ic1, ia2, ib2, ic2]T with the following conversion rule:
idq=T1(θ)iabc.
The above-presented description is related to an electrical machine comprising two three-phase stator windings. It is, however, straightforward to generalize the above-presented equations (1)-(6) to a general case where there are N M-phase stator windings. Also in the general case, the voltage equations in the matrix form would be similar to the equation (3) shown above. The inductance matrix would Ldq be a 2N×2N matrix, the vector of the rotation-converted stator currents idq would have 2N elements, and the vector of the rotation-converted stator voltages udq would have 2N elements.
As can be seen from equations (2), the d- and q-directional stator flux linkages of the first three-phase stator winding Ψd1 and Ψq1 are dependent on also the currents id2 and iq2 of the second three-phase stator winding. Correspondingly, the Band q-directional stator flux linkages of the second three-phase stator winding Ψd2 and Ψq2 are dependent on also the currents id1 and iq1 of the first three-phase stator winding. This cross-dependency complicates the control of the currents of the electrical machine. In the matrix form equation (3), the cross-dependency is manifested by the non-zero off-diagonal elements of the inductance matrix Ldq.
In conjunction with embodiments of the present invention, the problem associated with the above-mentioned cross-dependency is avoided by using such form of the model of the electrical machine where each of the stator flux-linkages of the model is dependent on only one of the rotation-converted stator currents in spite of the mutual inductances between the two or more multiphase stator windings when the inductance parameters of the model are not altered e.g. due to magnetic saturation.
The above-mentioned form of the model of the electrical machine can be derived, for example, from the matrix formulation presented by the equation (3).
The inductance matrix Ldq is symmetric and thus it can be converted into a diagonal form with the aid of a second conversion matrix T2:
LDQ=T2LdqT2T) (7)
where LDQ is the diagonalized inductance matrix, and the second conversion matrix T2 is [g1 g2 . . . gn] where g1 g2 . . . gn are the eigenvectors of the inductance matrix Ldq. For example, in the case where there are two multiphase stator windings, the second conversion matrix T2 can be presented in the form:
where C2 is a constant scaling coefficient. For example, C2=1/√2 makes the second conversion matrix orthonormal so that T2T2T=T2TT2 is a unit matrix diag(1, 1, 1, 1).
The matrix formulation represented by equation (3) can be multiplied from the left side with the second conversion matrix T2:
where the rule that T2TT2 is a unit matrix is utilized too.
Equation (9) can be written as:
where:
As the inductance matrix LDQ is a diagonal matrix, each of the stator flux-linkages of this form of the model is dependent on only one of the rotation-converted stator currents in spite of the mutual inductances between the two or more multiphase stator windings. In conjunction with a current control, the voltage vector caused by the rotation of the rotor:
ωJmLDQIDQ+ωJmΨPMDQ
can be taken into account as a feed-forward term uff that is superposed to an output of a current regulating device. Removing the above-mentioned feed-forward term from the equation (10) and then solving, in the Laplace domain, for the vector iDQ of the rotation-converted stator currents yields the following system to be controlled with the current regulating device:
iDQ(s)=(RsI+sLDQ)−1uDQ(s), (11)
where s is the Laplace variable and I is a unit matrix having the same dimensions as the inductance matrix LDQ. The vector iDQ(S) of the rotation-converted stator currents represents the quantity to be controlled, and the vector uDQ(s) of the rotation-converted stator voltages represents the quantity with the aid of which the control is carried out. The rows of the equation (11) represent mutually independent first order systems because the inductance matrix LDQ is a diagonal matrix. For example, the nth row of the equation (11) is:
where iDQ,n(s) is the nth rotation-converted stator current, i.e. the nth element of the vector iDQ(s), uDQ,n(s) is the nth rotation-converted stator voltage, i.e. the nth element of the vector uDQ(s), and LDQ,n is the diagonal element of the nth row, or column, of the inductance matrix LDQ.
In an exemplifying case where there are two multiphase stator windings, the inductance matrix LDQ is
where C2 is the constant scaling coefficient presented in equation (8), and the uDQ, iDQ, ΨPMDQ and Jm are:
The equation (10), when written in the component form, yields:
where the flux linkages of the model are:
As shown by the equations (14), each of the stator flux-linkages ΨD1, ΨQ1, ΨD2 and ΨQ2 of this form of the model is dependent on only one of the rotation-converted stator currents iD1, iQ1, iD2, and iQ2 when the inductance parameters LD1, LQ1, LD2, and LQ2 of the model are constant, i.e. unaltered. Magnetic saturation may cause that e.g. the rotation-converted stator current iD2 affects the value of e.g. LD1 and thus iD2 may affect indirectly to ΨD1. However, when the inductance parameters are unaltered, each stator flux linkage is dependent on only one rotation-converted stator current.
The voltages caused by the rotation of the rotor are:
Therefore, the mutually independent first order systems to be regulated with a current regulating device are:
The rotation-converted stator currents iD1, iQ1, iD2, and iQ2 are bound to the stator phase currents ia1, ib1, ic1, ia2, ib2, and ic2 via a third conversion matrix T3(θ):
iDQ=T3(θ)iabc=T2T1(θ)iabc, (17)
iabc=T3T(θ)iDQ, (18)
where T1(θ) is the first conversion matrix defined in equations (5) and (6). Correspondingly, the rotation-converted stator voltages uD1, uQ1, uD2, and uQ2 are bound to the stator phase voltages ua1, ub1, uc1, ua2, ub2, and uc2 via the third conversion matrix T3(θ):
uDQ=T3(θ)uabc, (19)
An apparatus according to an exemplifying embodiment of the invention comprises:
The processing system 201 is configured to produce the voltage control signals u1 and u2 on the basis of the rotation-converted stator currents iD1, iQ1, iD2, and iQ2 which are expressed in a coordinate system bound to a rotor of the electrical machine and on the basis of a model of the electrical machine modeling at least inductances of the two multiphase stator windings and mutual inductances between the two multiphase stator windings. The stator phase currents ia1, ib1, ic1, ia2, ib2, ic2 are converted into the rotation-converted stator currents iD1, iQ1, iD2, and iQ2 with a conversion rule corresponding to a form of the model of the electrical machine where each of the stator flux-linkages ΨD1, ΨQ1, ΨD2 and ΨQ2 of the model is dependent on only one of the rotation-converted stator currents iD1, iQ1, iD2, and iQ2 when the inductance parameters of the model are unaltered. A conversion rule of the kind mentioned above is presented with equation (17). A functional block 205 shown in
In the exemplifying embodiment illustrated in
[u1u2]T=C3T3T(θ)[uD1,ref,uQ1,ref,uD2,ref,uQ2,ref]T, (21)
where C3 is a constant scaling coefficient, and the voltage control signals u1 and u2 are assumed to be row-vectors of values proportional to the target phase voltages of the two multiphase stator windings. It is also possible to convert the voltage target values uD1,ref, uQ1,ref, uD2,ref, and uQ2,ref into voltage values expressed in the α-β coordinate system illustrated in
In an apparatus according to an exemplifying embodiment of the invention, the processing system 201 is configured to determine the target values of the rotation-converted stator currents iD1,ref, iQ1,ref, iD2,ref, and iQ2,ref at least partly on the basis of a desired torque Tref. By a straightforward analysis it can be shown that the torque Te produced by the electrical machine can be estimated with the following equation:
Te=CTp(ΨPMD1iQ1+(LD1−LQ1)iD1iQ1+(LD2−LQ2)iD2iQ2), (22)
where p is the number of pole pairs of the electrical machine and CT is a constant scaling coefficient.
In the exemplifying case shown in
where C2 is the constant scaling coefficient defined in the equation (8) and id1 is the direct-axis component of a first current space-vector created by the stator phase currents ia1, ib1, and ic1 of the first one of the two multiphase stator windings, iq1 is the quadrature-axis component of the first current space-vector, id2 is the direct-axis component of a second current space-vector created by the stator phase currents ia2, ib2, and ic2 of the second one of the two multiphase stator windings, and iq2 is a quadrature-axis component of the second current space-vector.
The target values of the rotation-converted stator currents iD2 and iQ2 are advantageously, but not necessarily, selected to be zeros, i.e. iD2,ref=iQ2,ref=0. This selection is actually an attempt to achieve symmetrical loading between the two multiphase stator windings because id1=id2 and iq1=iq2 when iD2,ref=iQ2,ref=0. In this case, the reference value iQ1,ref can be given as:
The reference value iD1,ref can be selected, for example, so that the overall stator current is minimized and the stator voltage has a desired value. A functional block 214 shown in
In an apparatus according to an exemplifying embodiment of the invention, the processing system 201 is configured to use a measured position angle θ of the rotor in the conversion of the stator phase currents of the multiphase stator windings into the rotation-converted stator currents.
In an apparatus according to an exemplifying embodiment of the invention, the processing system 201 is configured to estimate the position angle θ of the rotor on the basis of one or more electrical quantities related to the electrical machine and to use the estimated position angle in the conversion of the stator phase currents of the multiphase stator windings into the rotation-converted stator currents. Methods for estimating the position angle θ of the rotor can be found, for example, from the following publications: T. Halkosaari: “Speed Sensorless Vector Control of Permanent Magnet Wind Power Generator—The Redundant Drive Concept,” in Wind Power. InTech, 2010, 558 p., and A. Piippo, M. Hinkkanen, and J. Luomi: “Sensorless Control of PMSM Drives Using a Combination of Voltage Model and HF Signal Injection” In Conference Record of the 39th IEEE-Industry Applications Society (IAS) Annual Meeting, vol 2, pp 964-970, Seattle, Wash., USA.
The processing system 201 can be implemented with one or more processor circuits, each of which can be a programmable processor circuit provided with appropriate software, a dedicated hardware circuit such as, for example, an application specific integrated circuit “ASIC”, or a configurable hardware circuit such as, for example, a field programmable gate array “FPGA”. For example, some of the functional blocks 205-211, and 214 can be implemented with a dedicated or configurable hardware circuit and some of these functional blocks can be implemented with one or more programmable processor circuits, or all of these functional blocks can be implemented with one or more programmable processor circuits or with dedicated or configurable hardware circuits. The present invention is not limited to any methods of implementation.
In a method according to an exemplifying embodiment of the invention, the stator phase currents of two multiphase stator windings are converted into four rotation-converted stator currents, and inductance coefficients between the stator flux-linkages and the rotation-converted stator currents are proportional to: Ld+Md, Lq+Mq, Ld−Md, and Lq−Mq, where Ld is a direct-axis inductance, Md is a direct-axis mutual inductance between the two multiphase stator windings, Lq is a quadrature-axis inductance, and Mq is a quadrature-axis mutual inductance between the two multiphase stator windings.
In a method according to an exemplifying embodiment of the invention, the four rotation-converted stator currents are proportional to id1+id2, iq1+iq2, iq1−iq2, and id2−id1, respectively, where id1 is the direct-axis component of a first current space-vector created by the stator phase currents of the first one of the two multiphase stator windings, iq1 is the quadrature-axis component of the first current space-vector, id2 is the direct-axis component of a second current space-vector created by the stator phase currents of the second one of the two multiphase stator windings, and iq2 is the quadrature-axis component of the second current space-vector.
In a method according to an exemplifying embodiment of the invention, the target values of the rotation-converted stator currents which are proportional to iq1−iq2 and to id2−id1 are zeroes so as to achieve symmetrical loading between the two multiphase stator windings.
A method according to an exemplifying embodiment of the invention comprises measuring a position angle of the rotor and using the measured position angle in the conversion of the stator phase currents of the two or more multiphase stator windings into the rotation-converted stator currents.
A method according to an exemplifying embodiment of the invention comprises estimating a position angle of the rotor on the basis of one or more electrical quantities related to the electrical machine and using the estimated position angle in the conversion of the stator phase currents of the two or more multiphase stator windings into the rotation-converted stator currents.
In a method according to an exemplifying embodiment of the invention, the controlling of the power stages comprises:
In a method according to an exemplifying embodiment of the invention, the regulating of the rotation-converted stator voltages comprises:
The regulator may comprise, for example, a separate proportional-integrating “PI” or proportional-integrating-derivative “P I D” regulator for each of the differences between the rotation-converted stator currents and their target values.
A method according to an exemplifying embodiment of the invention comprises determining the target values of the rotation-converted stator currents at least partly on the basis of a desired torque.
A computer program according to an exemplifying embodiment of the invention comprises software modules for the purpose of controlling an electrical machine comprising two or more multiphase stator windings. The software modules comprise computer executable instructions for controlling a programmable processor to:
The software modules can be e.g. subroutines or functions implemented with a suitable programming language and with a compiler suitable for the programming language and the programmable processor.
A computer program product according to an exemplifying embodiment of the invention comprises a computer readable medium, e.g. a compact disc “CD”, encoded with a computer program according to an embodiment of invention.
A signal according to an exemplifying embodiment of the invention is encoded to carry information defining a computer program according to an embodiment of invention.
The specific examples provided in the description given above should not be construed as limiting. Therefore, the invention is not limited merely to the embodiments described above.
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/FI2012/050504 | 5/25/2012 | WO | 00 | 1/29/2015 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2013/175050 | 11/28/2013 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
5763976 | Huard | Jun 1998 | A |
20100123426 | Nashiki | May 2010 | A1 |
20120146626 | Bieler | Jun 2012 | A1 |
Number | Date | Country |
---|---|---|
101359888 | Feb 2009 | CN |
2 020 334 | Feb 2009 | EP |
2 412 091 | Jun 2014 | EP |
11-262293 | Sep 1999 | JP |
11262293 | Sep 1999 | JP |
2000-004600 | Jan 2000 | JP |
2000004600 | Jan 2000 | JP |
WO 2012046115 | Apr 2012 | WO |
Entry |
---|
Decoupling algorithm in view of multivariable electrical drives control GE44-LARGE (Laboratoire Atlantique de Recherche en GCnie Electrique) CRTT, BP 406, Bd de l'universitt! 44602 Saint-Nazaire Cedex, France, DOC 2002. |
English language translation of of the Notice of Reasons for Rejection issued in corresponding Japanese Patent Application 2015-513233 on Mar. 22, 2016. |
Benkhoris et al., “Decoupling algorithm in view of multivariable electrical drives control”, Proceedings of the American Control Conference, XP010597696, May 8, 2002, vol. 6, pp. 4861-4866. |
Search report issued in corresponding Chinese Application No. 201280073467.8 on May 30, 2016 (1 page). |
Number | Date | Country | |
---|---|---|---|
20150229261 A1 | Aug 2015 | US |