Method and apparatus for operating an electric machine, electric machine
10333449 ยท 2019-06-25
Assignee
Inventors
Cpc classification
H02P23/14
ELECTRICITY
International classification
H02P23/14
ELECTRICITY
H02P21/14
ELECTRICITY
Abstract
A method for operating an electric machine comprises ascertaining an electric stator current of the electric machine, calculating a magnetic flux which is generated based on the ascertained electric stator current, and controlling a torque generated by means of the electric machine, as a function of the calculated magnetic flux.
Claims
1. A method for operating an electric machine, comprising the steps: ascertaining an electric stator current of the electric machine, calculating a magnetic flux which is generated based on the ascertained electric stator current, the calculating of the magnetic flux comprising approximating the magnetic flux by means of a polynomial which is of at least the second degree, and controlling a torque generated by means of the electric machine, as a function of the calculated magnetic flux.
2. The method as claimed in claim 1, wherein the polynomial is of the third order and the absolute term of the polynomial is equal to zero.
3. The method as claimed in claim 1, wherein a differential inductance resulting from the corresponding stator current is measured in each case at different stator currents, wherein the coefficients of the polynomial are ascertained based on the measured differential inductances.
4. The method as claimed in claim 3, wherein for ascertaining the coefficients, it is specified that the linear term of the polynomial is equal to 1 or is equal to the torque constant of the electric machine.
5. A method for operating an electric machine, comprising the steps: ascertaining an electric stator current of the electric machine, calculating a magnetic flux which is generated based on the ascertained electric stator current, the ascertained stator current being divided into a torque-generating current and into a field-generating current by means of the Park transformation, wherein for calculating the magnetic flux, only the torque-generating current is used, and controlling a torque generated by means of the electric machine, as a function of the calculated magnetic flux.
6. A method for operating an electric machine, comprising the steps: ascertaining an electric stator current of the electric machine, calculating a magnetic flux which is generated based on the ascertained electric stator current, and controlling a torque generated by means of the electric machine, as a function of the calculated magnetic flux, wherein an actual torque corresponding to the calculated magnetic flux is calculated, which, is compared to a specified setpoint torque, wherein the control of the torque generated by means of the electric machine is carried out based on the comparison.
7. An apparatus for operating an electric machine, comprising: an ascertainment device for ascertaining an electric stator current of the electric machine, a processor for calculating a magnetic flux which is generated based on the ascertained electric stator current, the processor being configured to calculate the coefficients of a polynomial which approximates the magnetic flux which is generated based on the ascertained electric stator current, based on inductances measured at different stator currents, and a control device for controlling a torque generated by means of electric machine, as a function of the calculated magnetic flux.
8. The apparatus for operating an electric machine as claimed in claim 7, wherein the ascertainment device comprises a current-measuring device for measuring the electric stator current.
9. The apparatus for operating an electric machine as claimed in claim 7, comprising an inductance-measuring device for measuring the inductance, which results from the electric stator current.
10. The apparatus for operating an electric machine as claimed in claim 7, being a converter.
11. An electric machine comprising: a rotor and a stator, a current-measuring device for measuring an electric stator current, and a control loop for receiving the measured stator current as an input value and for subsequently calculating a torque, wherein the torque is compared to a specified setpoint torque, wherein based on a comparison difference a control value is calculated which ascertains a stator current to be impressed on the stator of the electric machine, wherein the control loop comprises a functional block for generating a PI control value or a PID control value based on the comparison difference, a functional block for ascertaining a pulse-width modulation value based on the PI control value PID control value and a functional block for ascertaining the stator current based on the pulse-width modulation value.
12. The electric machine as claimed in claim 11, wherein the control loop comprises a functional block for calculating the torque and a functional block for comparing the torque to the specified setpoint torque.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) The invention is explained in more detail below with reference to exemplary embodiments, in which:
(2)
(3)
(4)
(5)
(6)
(7)
DETAILED DESCRIPTION
(8) The present invention is based on the finding that the aforementioned object may be achieved in that the magnetic flux .sub.m(I) caused by the stator current, rather than the stator current I, is used as a controlled variable for controlling the torque as a target value.
(9) As a result, the effects of the magnetic saturation due to the stator current may be efficiently taken into consideration, thus making the corresponding calculations and the corresponding control more efficient, so that as a result, the electric machine may be operated more efficiently.
(10) In one embodiment, the electric machine is a three-phase machine, for example, an asynchronous machine or a synchronous machine, a direct-current machine, or a linear motor.
(11) In one embodiment, it is provided that the calculation of the magnetic flux comprises approximating the magnetic flux by means of a polynomial which is of at least the second degree.
(12) As a result, for example, the technical advantage is achieved that the magnetic flux may be efficiently calculated.
(13) According to one embodiment, the degree of the polynomial is equal to three.
(14) In another embodiment, it is provided that the polynomial is of the third degree and the absolute term of the polynomial is equal to zero.
(15) As a result, for example, the technical advantage is achieved that the magnetic flux may be efficiently calculated.
(16) In one embodiment, the polynomial is a Taylor polynomial.
(17) According to another embodiment, it is provided that a differential inductance resulting from the corresponding stator current is measured in each case at different stator currents, wherein the coefficients of the polynomial are ascertained based on the measured differential inductances.
(18) As a result, for example, the technical advantage is achieved that the coefficients of the polynomial may be efficiently calculated.
(19) It is thus in particular provided to measure the corresponding differential inductances as a function of the stator current, via a variation of the operating point (different stator currents), wherein the coefficients are calculated from the differential inductances.
(20) In one embodiment, it is provided that for ascertaining the coefficients, it is specified that the linear term of the polynomial is equal to 1 or is equal to the torque constant of the electric machine.
(21) As a result, for example, the technical advantage is achieved that the coefficients of the polynomial may be efficiently calculated.
(22) The ascertainment of the coefficients will be described below in greater detail by way of example, using a third-order polynomial (degree=3) and under the stipulation that the absolute term of the polynomial is equal to zero:
(23) The polynomial is:
.sub.m(I)b.sub.1.Math.I+b.sub.2.Math.I.Math.|I|+b.sub.3.Math.I.sup.3
where .sub.m is the magnetic flux, I is the stator current, and b.sub.1, b.sub.2, b.sub.3 are the coefficients of the polynomial.
(24) For reasons of symmetry, the function .sub.m(I) is an odd (point-symmetric) function. I.sup.2 has therefore been replaced by I.Math.|I|, and the coefficient b.sub.0 (absolute term) is not required (=0). According to one embodiment, an approximation using a polynomial of a higher order is optionally provided.
(25) If a characteristic curve having many experimentally ascertained data points is available, an approximation to a third-order polynomial may efficiently compensate for random measuring errors, and results in a smooth profile similar to a linear regression. In practice, an approximation to a second-order polynomial is sufficiently accurate, and is therefore provided according to one embodiment.
(26) The coefficients are, for example, ascertained via curve fitting. An approximation of the curve (the polynomial) according to the method of least squares is particularly suitable and is therefore also thus provided according to one embodiment.
(27) Advantageously and according to one embodiment, a simplification according to Horner's method (factoring out) is thus also provided:
.sub.m(I)I.Math.(b.sub.1+|I|.Math.(b.sub.2+b.sub.3.Math.|I|+ . . . ))
(28) As a result, the number of required calculation operations (multiplication, addition) decreases.
(29) One advantage of the polynomial is that the algorithm for calculation is less complex than an algorithm for linear interpolation of tabular data points. This applies both to a process-based implementation (for example, in the C programming language) and to a VHDL (very high-speed integrated circuit hardware description language) implementation, or rather a FPGA (field-programmable gate array) implementation, or an ASIC (application-specific integrated circuit) implementation.
(30) In one embodiment, the apparatus has a function for creating, or rather depicting, a Bode diagram (Bode plot) of the control loops. Preferably, a frequency-dependent transfer function of the current control loop of the electric machine is measured open-loop and closed-loop in a small-signal range.
(31) Preferably, the transfer function is measured for multiple, i.e., in particular for a series of, frequencies f or , and, for example, is graphically depicted as an amplitude response (magnitude) and a phase response in the form of a Bode plot.
(32) The open-loop amplitude response |G.sub.0(j)| thus measured comprises the frequency-dependent transfer function of the controller G.sub.controller(j) and of the controlled system G.sub.controlled system(j).
(33) Alternatively or in addition, according to one embodiment, it is provided that an open-loop transfer function G.sub.0(I) is measured at a constant frequency .sub.1 (this is the frequency which is also used below in the formula
(34)
and as a function of a current I (operating point) causing saturation.
(35) Since the frequency remains unchanged, the transfer function of the controller does not change. As a result, the effect of the saturation on the controlled system may be measured in a simple manner.
(36) The excitation for measurement takes place, for example, as a small signal, for example, sinusoidally at, for example, 500 Hz, about an operating point, defined by the current I.
(37) For example,
(38) The following hold true: I.sub.0: Open-loop transfer function measured at I.sub.0=0 A, and I.sub.n: Open-loop transfer function measured at a point 607.
(39) An additional advantage is that a uniform approach may be used for the parameterization. A small-signal gain K.sub.diff of the characteristic curve may be calculated closed, according to the processing rules for polynomials, because an increase in the tangent at a point of the polynomial is proportional to the differential inductance:
(40)
(41) For reasons of symmetry, the function K.sub.diff(I) must be an even (axis-symmetric) function. I.sup.1 has therefore been replaced by |I|.
(42) Dividing by b.sub.1 results in:
(43)
(44) It is also possible, and according to one embodiment, it is also thus provided, to determine these coefficients b.sub.1, b.sub.2, b.sub.3 within the scope of an automated parameter ascertainment in a comparatively simple manner, internally in the apparatus, or rather in the electric machine, without a test bench.
(45) One possibility for determining the coefficients is the measurement of an open-loop current controller Bode plot with subsequent curve fitting:
(46)
G.sub.0 is calculated from G.sub.controller multiplied by G.sub.controlled system:
(47)
(48) By varying the operating point (I), the differential inductance L.sub.diff(I) may be measured as a function of the stator current I.
(49) With the specification b.sub.1=1, the coefficients b.sub.2 and b.sub.3 may be calculated in a particularly simple manner from the measured L.sub.diff(I) data.
(50) Alternatively, b.sub.1=K.sub.T (K.sub.T is the torque constant of the electric machine without considering magnetic saturation) may also be specified before b.sub.2 and b.sub.3 are calculated.
(51) In one embodiment, it is provided that the ascertained stator current is divided into a torque-generating current and into a field-generating current by means of the Park transformation, wherein for calculating the (torque-generating) magnetic flux, only the torque-generating current is used.
(52) As a result, for example, the technical advantage is achieved that the magnetic flux may be efficiently calculated.
(53) For calculating the magnetic flux, it is also in particular provided to ignore the field-generating current.
(54) Experimental tests have surprisingly shown that the constant permanently excited flux (d-axis) which is perpendicular to the torque-generating flux (q-axis) may be ignored.
(55) Magnetic leakage fields have an effect on the winding inductance, but not on the torque. Here as well, tests have surprisingly shown that these leakage fields may be ignored.
(56) Thus, according to one embodiment, it is provided that in the embodiments made above and below, the stator current I may be replaced by the torque-generating current I.sub.q.
(57) In an additional embodiment, it is provided that a torque corresponding to the calculated magnetic flux is calculated, which, being an actual torque, is compared to a specified setpoint torque, wherein the control of the torque provided by means of the electric machine is carried out based on the comparison.
(58) As a result, for example, the technical advantage is achieved that the control may be efficiently carried out.
(59) In one embodiment, it is provided that the method for operating an electric machine is implemented or carried out by means of the apparatus for operating an electric machine.
(60) In one embodiment, it is provided that the apparatus for operating an electric machine is configured to implement or carry out the method for operating an electric machine.
(61) Technical functionalities of the method for operating an electric machine result analogously from corresponding technical functionalities of the apparatus for operating an electric machine, and vice-versa.
(62) Apparatus features result analogously from corresponding method features, and vice-versa.
(63) The formulation or rather comprises in particular the formulation and/or.
(64)
(65) The method includes the following steps: ascertaining 101 an electric stator current of the electric machine, calculating 103 a magnetic flux which is generated based on the ascertained electric stator current, and controlling 105 a torque generated by means of the electric machine, as a function of the calculated magnetic flux.
(66) In one embodiment, the ascertainment 101 comprises measuring the electric stator current of the electric machine. The ascertained stator current is based, for example, on the measured electric stator current. The ascertained stator current corresponds, for example, to the measured electric stator current.
(67) The ascertained stator current corresponds, for example, to a current to be impressed on the stator.
(68) In one embodiment, it is provided that the ascertained stator current is divided into a torque-generating current and into a field-generating current by means of the Park transformation. In one embodiment, it is provided that only the torque-generating current is used for calculating the magnetic flux. The field-generating current is, for example, ignored.
(69)
(70) The apparatus 201 comprises: an ascertainment device 203 for ascertaining an electric stator current of the electric machine, a processor 205 for calculating a magnetic flux which is generated based on the ascertained electric stator current, and a control device 207 for controlling a torque generated by means of the electric machine, as a function of the calculated magnetic flux.
(71) The ascertainment device 203 comprises, for example, a current-measuring device for measuring the electric stator current. The measuring device comprises, for example, a current sensor.
(72) The apparatus 201 comprises, for example, an inductance-measuring device for measuring an inductance, or rather a differential inductance, which results from the electric stator current.
(73) The processor 205 is, for example, configured to calculate the coefficients of a polynomial which approximates the magnetic flux which is generated based on the ascertained electric stator current, based on (for example, differential) inductances measured at different stator currents.
(74) According to one embodiment, the control device 207 comprises the functional blocks 417 and 421 according to
(75) In the context of the description, a combination made up of an electric machine and a converter may be referred to as a drive. For example, a combination made up of a three-phase machine and a converter may be referred to as a three-phase drive.
(76)
(77) The electric machine 301 comprises the apparatus 201 of
(78) Electric machines as such are known to those skilled in the art. Therefore, details, for example, an armature and a rotor of the electric machine 301, are not depicted.
(79) An electric machine in the context of the description thus comprises in particular a rotor and a stator.
(80) The stator comprises in particular one or more coils.
(81) The rotor comprises in particular one or more coils and/or one or more permanent magnets.
(82)
(83) In a feedback path 403 of the control loop 401, according to a functional block 407, an ascertained electric stator current 405, which may also be referred to by I in the context of the description, is converted into a magnetic flux .sub.m, or rather into a torque M. For this purpose, for example, a polynomial is used to approximate the magnetic flux. For example, the polynomial is of the third degree.
(84) For example, it is provided that the electric stator current is divided into a torque-generating current, which may also be referred to by I.sub.q in the context of the description, and into a field-generating current, which may also be referred to by I.sub.d in the context of the description, using the Park transformation. In the calculations above and below, in this respect, according to one embodiment, the stator current I is replaced by the torque-generating current I.sub.q.
(85) In the context of this description, it is, for example, assumed that in a first approximation, the magnetic flux is proportional to the torque M which is provided by means of the electric machine. Thus, for example, the following relationship holds true:
M(I)=C.sub.m.Math..sub.m(I),
where C.sub.m is a torque constant of the electric machine.
(86) The stator current is thus initially converted in the feedback path 403 into the associated torque-proportional magnetic flux, or rather into a torque, so that the nonlinear formation of a desired variable
I.sub.setpoint=f(M.sub.setpoint)
(87) may be replaced by a linear calculation:
(88)
(89) In the apparatus, or rather in the electric machine, for example, a respective characteristic curve is stored for .sub.m(I) and/or M(I).
(90) A particular advantage of the concept according to the present invention is in particular that via this .sub.m(I), or rather M(I) characteristic curve, the current-dependent saturation of the inductance is also simultaneously considered. Exactly as the differential inductance decreases with increasing current, a small-signal gain K.sub.diff of the .sub.m(I) characteristic curve also decreases accordingly:
(91)
(92) An additional current-dependent adaptation of a current regulator Kr for controlling an electric stator current is no longer required. In one exemplary field-oriented control of the multiphase motor of the electric machine, the torque-generating current I.sub.q is accordingly converted into a magnetic flux, or rather into a torque.
(93) Using a third-order polynomial, it will be described in greater detail how these characteristic curves are calculated or ascertained according to one embodiment for .sub.m(I) and/or M(I), wherein
(94) In one embodiment, it is provided that the function .sub.m(I) is approximated by means of a third-order polynomial (degree=3), wherein it is further specified that the absolute term b.sub.0 of the polynomial is equal to zero.
.sub.m(I)b.sub.1.Math.I+b.sub.2.Math.I.Math.|I|+b.sub.3.Math.I.sup.3
where .sub.m is the magnetic flux, I is the stator current, and b.sub.1, b.sub.2, b.sub.3 are the coefficients of the polynomial.
(95) For reasons of symmetry, the function .sub.m(I) is an odd (point-symmetric) function. I.sup.2 has therefore been replaced by I.Math.|I|, and the coefficient b.sub.0 (absolute term) is not required (=0). According to one embodiment, an approximation using a polynomial of a higher order is optionally provided.
(96) If a characteristic curve having many experimentally ascertained data points is available, an approximation to a third-order polynomial may efficiently compensate for random measuring errors, and results in a smooth profile similar to a linear regression. In one embodiment, a second-order polynomial is used for the approximation.
(97) The coefficients of the polynomial are, for example, ascertained via curve fitting. An approximation of the curve (the polynomial) according to the method of least squares is particularly suitable and is therefore also thus provided according to one embodiment.
(98) Advantageously and according to one embodiment, a simplification according to Horner's method (factoring out) is thus also provided:
.sub.m(I)I.Math.(b.sub.1+|I|.Math.(b.sub.2+b.sub.3.Math.|I|+ . . . ))
(99) As a result, the number of required calculation operations (multiplication, addition) decreases in an efficient manner.
(100) One advantage of the polynomial is that the algorithm for the calculation is less complex than an algorithm for the linear interpolation of tabular data points. This applies both to a process-based implementation (for example, in the C programming language) and to a VHDL (very high-speed integrated circuit hardware description language) implementation, or rather a FPGA (field-programmable gate array) implementation, or an ASIC (application-specific integrated circuit) implementation.
(101) An additional advantage is that a uniform approach may advantageously be used for the parameterization. The small-signal gain K.sub.diff of the characteristic curve may be calculated closed, according to the processing rules for polynomials:
(102)
(103) For reasons of symmetry, the function K.sub.diff(I) must be an even (axis-symmetrical) function. I.sup.1 has therefore been replaced by |I|.
(104) Dividing by b.sub.1 results in:
(105)
(106) Within the scope of an automated parameter ascertainment, the determination of these coefficients b.sub.1, b.sub.2, b.sub.3 is also possible internally in the apparatus, or rather in the electric machine, in a comparatively simple manner, without a test bench.
(107) One possibility for determining the coefficients is the measurement of an open-loop current controller Bode plot with subsequent curve fitting:
(108)
(109) By varying the operating point (I), the differential inductance L.sub.diff(I) may be measured as a function of the stator current I.
(110) With the specification b.sub.1=1, the coefficients b.sub.2 and b.sub.3 may be calculated in a particularly simple manner from the measured L.sub.diff(I) data.
(111) Alternatively, b.sub.1=K.sub.T (K.sub.T is a torque constant of the electric machine) may also be specified before b.sub.2 and b.sub.3 are calculated.
(112) Thus, in the case of different stator currents, it is provided to measure the resulting differential inductances, in order to accordingly calculate the coefficients of the polynomial. Subsequently, the accordingly ascertained characteristic curves are available, in order, with knowledge of the electric stator current of the electric machine, to ascertain the associated torque, as is carried out in functional block 407.
(113) The functional block 407 thus receives the ascertained stator current 405 as an input value and, as described above, subsequently calculates the associated torque, which is indicated symbolically in
(114) According to a functional block 413, the torque {circumflex over (M)} is compared to a specified setpoint torque M*, denoted in
(115) The PI control value is thus, for example, created via a PI controller. For example, other controllers, for example, a PID controller (proportional integral derivative controller), are provided.
(116) The control device 207 according to
(117) The PI control value 419 is supplied to a functional block 421, which ascertains a pulse-width modulation value 423 which is based on the PI control value 419 and which may also be denoted by U.
(118) The PWM value 423 is supplied to a functional block 425 which ascertains a stator current I to be impressed on the stator of the electric machine, based on the PWM value 423, said stator current being denoted in
(119) The ascertained stator current according to the box 427 is, for example, impressed on the stator of the electric machine.
(120) For ascertaining the stator current according to the box 427, the following formula is used in the functional block 425:
(121)
(122) The term
(123)
represents, for example, the transfer function of a motor winding in the Laplace domain.
(124) The voltage of the PWM is provided as an input value. The current is provided as an output value. The induced voltage is not depicted here. The Laplace operated is denoted by s.
(125) The term
(126)
is thus a modeling of the motor winding of an electric machine in the Laplace domain (also referred to as the image domain). The s is used due to the depiction in the Laplace domain. The motor winding constitutes a so-called PT1 element.
(127) In control technology, the following depiction of transmission elements is generally used:
(128)
(129) The prior art is in particular interpolated using a table:
(130)
where R is the electrical resistance of the stator, i.e., the coil of the stator; L.sub.diff is the differential inductance, and s is the Laplace operator.
(131) The current control takes place via an adjustable voltage. In order to achieve proper efficiency, the voltage is generated via PWM (pulse width modulation). The small-signal behavior about the operating point is important for the stability of the control loop. Thus, the inductance Liar measured at the operating point is applied.
(132) The core of the present invention thus comprises in particular the idea that for controlling the torque M as a target value, the magnetic flux .sub.m(I) generated by the torque-generating current I.sub.q, rather than the stator current I, is used as a controlled variable. The control loop 401 may thus in particular be referred to as a flux control loop, or rather a torque control loop, which advantageously takes saturation effects into consideration. The controlled variable of the control loop 401 is the magnetic flux, or rather the torque.
(133) The concept according to the present invention has in particular the advantage that only one approximation, for example, by means of the polynomial, must be carried out in order to efficiently take the saturation effects into consideration. In contrast, two different interpolations are required in the prior art.
(134) An additional advantage may thus in particular be seen in the fact that, for considering the saturation effects, only one characteristic curve must be ascertained. In contrast, in the known prior art, two different characteristic curves must be ascertained.
(135) A resulting advantage is in particular to be seen in a reduction in the processor load. In particular, memory may advantageously be saved by storing only one characteristic curve, compared to storing two characteristic curves.
(136)
(137) The abscissa 503 indicates the electric stator current in amperes. The ordinate 505 indicates the torque in Nm.
(138) The experimentally ascertained torques at different stator currents are indicated as squares having the reference numeral 507.
(139) The experimentally ascertained torques are approached or approximated by a curve which has the reference numeral 509. The curve 509 is a third-order polynomial and has the following coefficients, wherein b.sub.1=1.0 is fixedly specified:
b.sub.0=0; b.sub.1=1.0; b.sub.2=0.07 and b.sub.3=0.003.
(140) The approximation of the experimentally ascertained torques is, for example, carried out using the method of least squares. For this purpose, the resulting differential inductances were measured at different stator currents (cf.
(141)
(142) The abscissa 603 indicates the electric stator current in amperes. The ordinate 605 indicates the relative small-signal gain in relation to the small-signal gain at I=0 A (K.sub.diff (I=0 A)=1).
(143) The experimentally ascertained values for the small-signal gain at different stator currents are indicated as squares having the reference numeral 607. As described above, it holds true that
(144)
so that by measuring the differential inductance at different operating points (I) and under the (exemplary) stipulation that b.sub.1=1.0, the small-signal gain may be experimentally ascertained at different stator currents.
(145) The experimentally ascertained values for the small-signal gain are approached or approximated via a curve which has the reference numeral 609. The curve 609 corresponds to a second-order polynomial which is the derivative with respect to the stator current of the third-order polynomial according to
(146)
(147) The coefficients of this polynomial are thus the following coefficients:
b.sub.0=0; b.sub.1=1.0; b.sub.2=0.07 and b.sub.3=0.003.
(148) Thus, it is advantageously made possible that the coefficients may be efficiently determined or ascertained internally in the electric machine, without a test bench, by varying the operating point and measuring the corresponding differential inductances.
(149) In summary, the present invention provides an efficient concept which considers, addresses, and compensates for the two effects of the saturation together: the nonlinear relationship between the stator current and the torque, and the reduction in the winding inductance due to saturation, so that corresponding parameterizations and algorithms may be simply and efficiently designed, compared to the known prior art, which separately considers and separately compensates for the two effects of the saturation.