Method and apparatus for estimating torque

Abstract

A method of estimating the torque of an axle of a running motor, comprising receiving data indicative of the total electrical active power, P, supplied to the motor, determining the electromechanical power, P.sub.em, supplied to the axle using the data indicative of the total electrical active power, P, receiving data indicative of a rotor speed, n.sub.r, or rotor frequency, f.sub.r, of the axle, determining an angular rotor frequency, .sub.r, of the axle using the data indicative of rotor speed, n.sub.r, or rotor frequency, f.sub.r, of the axle, and determining electromechanical axle torque, T.sub.em, using the determined electromechanical power, P.sub.em, and the determined angular rotor frequency, .sub.r.

Claims

1. A processor-implemented method of estimating and communicating the torque of an axle of a running motor, comprising: obtaining, by a control unit comprised of one or more processors, data indicative of a total electrical active power, P, supplied to the motor, by measuring the total electrical active power, P, supplied to the motor; decomposing, by the control unit, the total electrical active power, P, into a fundamental frequency component, P.sub.1, and harmonic frequency components; decomposing, by the control unit, the fundamental frequency component, P.sub.1, into symmetrical components including at least a positive sequence power component, P.sup.+, and a negative sequence power component, P.sup.; determining, by the control unit, an electromechanical power, P.sub.em, supplied to the axle using the data indicative of the total electrical active power, P, by subtracting a sum of the negative sequence power component, P.sup., and harmonic frequency components that provide a negative contribution to mechanical power at the axle from a sum of the positive sequence power component, P.sup.+, and harmonic frequency components that provide a positive contribution to mechanical power at the axle; determining, by the control unit, a rotor speed, n.sub.r, or rotor frequency, f.sub.r, of the axle based on data indicative of a load dependent slip, s, and a synchronous speed, n.sub.s, or synchronous frequency, f.sub.s, of the motor; determining, by the control unit, an angular rotor frequency, .sub.r, of the axle using the data indicative of the rotor speed, n.sub.r, or rotor frequency, f.sub.r, of the axle; determining electromechanical axle torque, T.sub.em, using the determined electromechanical power, P.sub.em, and the determined angular rotor frequency, .sub.r; and outputting, by the control unit, a signal indicative of the determined electromechanical axle torque, T.sub.em.

2. The method according to claim 1, wherein measuring the total electrical active power, P, supplied to the motor includes measuring one or both of an instantaneous motor voltage and an instantaneous motor current on each connected phase of the motor.

3. The method of claim 1, wherein the fundamental and harmonic frequency components of the total electrical active power, P, are obtained by Fourier analysis of total electrical active power, P.

4. The method of claim 1, wherein determining the rotor speed, n.sub.r, or rotor frequency, f.sub.r, of the axle includes receiving data indicative of one or more of the fundamental electrical frequency, f.sub.1, of the current supplied to the motor, the load dependent slip, s, of the motor, the frequency shift, f, due to rotor slip, and the number of poles, p, of the motor.

5. The method of claim 4, wherein the rotor speed, n.sub.r, or rotor frequency, f.sub.r, of the axle is determined where f.sub.r(Hz)=(1-s).Math.f.sub.s, =(1-s).Math.2.Math.f.sub.1/p =n.sub.r(rpm)/60, and f.sub.s is the synchronous frequency of the motor.

6. The method of claim 4, wherein the method further includes determining the number of poles, p, of the motor, where p =INTEGER((120.Math.f.sub.nom)/n.sub.nom), where n.sub.nom is the nominal speed of the motor operating at its nominal frequency, f.sub.nom.

7. The method of claim 6, wherein the method further includes determining the load dependent slip, s, by analysis of a frequency spectrum of the current supplied to the motor.

8. The method of claim 7, wherein the frequency shift, f, due to rotor slip, and the fundamental electrical frequency, f.sub.1, of the current supplied to the motor are determined from the frequency spectrum of the current supplied to the motor, and the load dependent slip, s, is determined as s=2.Math.f/f.sub.1.

9. The method of claim 7, wherein the method further comprises creating a frequency spectrum of the current supplied to the motor.

10. The method of claim 4, wherein the angular rotor frequency, .sub.r, of the axle is determined by calculating .sub.r (rad/s)=2.Math..Math.f.sub.r=4.Math..Math.(1s).Math.f.sub.1/p.

11. The method of claim 1, wherein the electromechanical axle torque, T.sub.em, is determined by calculating T.sub.em=P.sub.em/.sub.r.

12. The method of claim 1, wherein the motor is an asynchronous induction motor.

13. The method of claim 1, wherein the positive and negative sequence power components, P.sup.+, P.sup., of the fundamental frequency component, P.sub.1, are obtained by symmetrical decomposition of the fundamental frequency component, P.sub.1.

14. An apparatus configured to estimate and communicate the torque of an axle of a running motor, the apparatus comprising one or more processing components configured to: obtain data indicative of a total electrical active power, P, supplied to the motor, wherein the total electrical active power, P, supplied to the motor is derived from measurement of voltage and current supplied to the motor; decompose the total electrical active power, P, into a fundamental frequency component, P.sub.1, and harmonic frequency components; decompose the fundamental frequency component, P.sub.1, into symmetrical components including at least a positive sequence power component, P.sup.+, and a negative sequence power component, P.sup.; determine an electromechanical power, P.sub.em, supplied to the axle using the data indicative of the total electrical active power, P, by subtracting a sum of the negative sequence power component, P.sup.and harmonic frequency components that provide a negative contribution to mechanical power at the axle from a sum of the positive sequence power component, P.sup.+, and harmonic frequency components that provide a positive contribution to mechanical power at the axle; determine a rotor speed, n.sub.r, or rotor frequency, f.sub.r, of the axle based on data indicative of a load dependent slip, s, and a synchronous speed, n.sub.s, or frequency, f.sub.s, of the motor; determine an angular rotor frequency, .sub.r, of the axle using the data indicative of the rotor speed, n.sub.r, or rotor frequency, f.sub.r, of the axle; and determine electromechanical axle torque, T.sub.em, using the determined electromechanical power, P.sub.em, and the determined angular rotor frequency, .sub.r, wherein the apparatus further comprises one or more components configured to display and/or transmit the determined electromechanical axle torque, T.sub.em.

15. The apparatus of claim 14, wherein the fundamental and harmonic frequency components of the total electrical active power, P, are obtained by Fourier analysis of total electrical active power, P.

16. The apparatus of claim 14, wherein the positive and negative sequence power components, P.sup.+, P.sup., of the fundamental frequency component, P.sub.1, are obtained by symmetrical decomposition of the fundamental frequency component, P.sub.1.

17. The apparatus of claim 14, wherein the electromechanical axle torque, T.sub.em, is determined by calculating T.sub.em=P.sub.em/.sub.r.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) Embodiments of the disclosure are further described hereinafter with reference to the accompanying drawings, in which:

(2) FIG. 1 shows a method according to an embodiment of the present disclosure;

(3) FIG. 2 shows a schematic diagram of an apparatus according to an embodiment of the present disclosure; and

(4) FIG. 3 shows a frequency spectrum associated with the motor current.

DETAILED DESCRIPTION

(5) Embodiments of the present disclosure provide a method and apparatus for estimating the torque of an axle of a running motor. Advantageously, the method and apparatus may be used to estimate torque without interrupting operation of the motor.

(6) FIG. 1 illustrates a method 10 of determining the electromechanical torque, T.sub.em, of an axle of a running motor.

(7) The method 10 may be performed by an apparatus 40 as illustrated in FIG. 2. The apparatus 40 may be positioned between a power source 30 and a motor 50.

(8) Returning to FIG. 1, the method 10 comprises a step 12 of receiving data that is indicative of the total electrical active power, P, supplied to the motor 50 by the power source 30. Additionally at step 12, data is received that is indicative of a rotor speed, n.sub.r, or rotor frequency, f.sub.r, of the axle. At step 14, the electromechanical power, P.sub.em, supplied to an axle of the motor 50 is determined using the data indicative of the total electrical active power, P. Additionally, at step 16, the angular rotor frequency, .sub.r, of the axle is determined using the data indicative of rotor speed, n.sub.r, or rotor frequency, f.sub.r, of the axle. Using the determined electromechanical power, P.sub.em, and the determined angular rotor frequency, .sub.r, the electromechanical axle torque, T.sub.em, is determined at step 20.

(9) Therefore, the electromechanical axle torque, T.sub.em, is estimated using only data relating to the total electrical active power, P, supplied to the motor 50 by the power source and the rotor speed, n.sub.r, or rotor frequency, f.sub.r, of the axle. Advantageously, this data may be received without interrupting the operation of the motor 50.

(10) In at least one embodiment, data indicative of the total electrical active power, P, supplied to the motor 50 is obtained by measuring the total electrical active power, P, supplied to the motor 50. For example, the total electrical active power, P, supplied to the motor 50 may be determined by measurement of the instantaneous motor voltage and instantaneous motor current on each connected phase of the motor 50. In particular, the total electrical active power, P, supplied to the motor 50 may be derived using the measured voltages and currents. For example, the voltage and/or current may be measured in accordance with IEC 61000-4-30 which is incorporated herein by reference, and/or the total electrical active power, P, may be calculated in accordance with IEEE Std 1459 which is also incorporated herein by reference.

(11) The mechanical power at the axle of a motor arises due to the electrical active power, P, supplied to the motor. However, not all electrical active power, P, is converted into mechanical power. The total electrical active power may therefore be decomposed into three categories, namely (i) electrical active power components positively contributing to mechanical power, (ii) electrical active power negatively contributing to mechanical power, and (iii) electrical active power components not contributing to mechanical power.

(12) The total electrical active power, P, may be decomposed into its individual frequency components. For example, this decomposition may be achieved by Fourier analysis. In a particularly preferable embodiment, this decomposition may be achieved by a fast Fourier transform (FFT) in accordance with the standard IEC 61000-4-7 which is incorporated herein by reference. Once decomposed, the total electrical active power, P, may be expressed as a series of harmonic components:
P=P.sub.1+P.sub.2+P.sub.3+P.sub.4. . . +P.sub.N,
where P.sub.1 is the power supplied by the fundamental frequency of the electrical voltage connected to the motor and P.sub.N is the power supplied by the Nth harmonic component of this fundamental frequency. The fundamental frequency is typically the mains frequency (e.g. 50 Hz or 60 Hz), however, in the case of electronic variable speed drive circuits, the fundamental frequency may be a varying frequency.

(13) The fundamental component P.sub.1 may be further decomposed into its symmetrical components using the established method symmetrical decomposition (originally proposed by C. L. Fortescue). In particular, the fundamental component P1 may be decomposed as follows:
P.sub.1=P.sup.++P.sup.+P.sup.0,
where P.sup.+ is the positive sequence power, P.sup. is the negative sequence power and P.sup.0 is the zero sequence power as described by Fortescue.

(14) Substituting this expression for P1 into the decomposed series for total active electrical power, P, and terminating the series at the 25.sup.th harmonic component gives:
P=P.sup.++P.sup.+P.sup.0+(P.sub.2+P.sub.3+P.sub.4 . . . +P.sub.25).

(15) As noted above, it is understood that certain terms of this expression provide a positive contribution to the mechanical power at the axle, whilst certain terms provide a negative contribution or do not contribute at all. For example, the fundamental wave revolves synchronously at the synchronous speed, n.sub.s, whereas the 3.sup.rd, 5.sup.th, 7.sup.th . . . harmonics may rotate in the forward or backward direction at speeds of n.sub.s/3, n.sub.s/5, n.sub.s/7, . . . , respectively. Even harmonics do not contribute to torque. The following table demonstrates how certain harmonics contribute positively or negatively (or not at all) to the torque (in respect of odd harmonics up to the 25.sup.th).

(16) TABLE-US-00001 Total phase Nearest whole Effective Harmonic shift (degrees) cycle (degrees) phase shift Torque 1 120 0 +120 Pos 3 360 360 0 n/a 5 600 720 120 Neg 7 840 720 +120 Pos 9 1080 1080 0 n/a 11 1320 1440 120 Neg 13 1560 1440 +120 Pos 15 1800 1800 0 n/a 17 2040 2160 120 Neg 19 2280 2160 +120 Pos 21 2520 2520 0 n/a 23 2760 2880 120 Neg 25 3000 2880 +120 Pos

(17) Considering the above, the following relationships apply (considering a limited series):
Positive contribution: P.sub.pos=P.sup.++P.sub.7+P.sub.13+P.sub.19;
Negative contribution: P.sub.neg=P.sup.+P.sub.5+P.sub.11+P.sub.17+P.sub.23;

(18) No contribution: all remaining power components.

(19) The power components not contributing to the mechanical power at the axle may, for example, generate heat, noise or mechanical malfunction, but contribute neither positively nor negatively to the mechanical power supplied to the axle.

(20) The electromechanical power, P.sub.em, may be expressed as:
P.sub.em=P.sub.posP.sub.neg=P.sup.++P.sub.7+P.sub.13+P.sub.19P.sup.P.sub.5P.sub.11P.sub.17P.sub.23.

(21) Therefore, one method for determining the electromechanical power, P.sub.em, is to measure or otherwise obtain data indicative of the total electrical active power, P, supplied to the axle (e.g., by measuring the instantaneous voltage and current as described above), and then decomposing the total active power, P, to obtain the components required to determine P.sub.em. Any method or apparatus in which data that is indicative of the total electrical active power, P, is received and then used to determine electromechanical power, P.sub.em, may be used within the scope of the present disclosure.

(22) As described above, step 12 also includes receiving data that is indicative of a rotor speed, n.sub.r, or rotor frequency, f.sub.r, of the axle (n.sub.r (rpm)=60.Math.f.sub.r (Hz)). The rotor speed, n.sub.r, may be expressed as:
n.sub.r=(1s).Math.n.sub.s,
where s is the load dependent slip of the motor and n.sub.s is the synchronous speed (i.e., the speed that the motor would run at without slip). Similarly, the rotor frequency, f.sub.r, may be expressed as:
f.sub.r=(1s).Math.f.sub.s,
where f.sub.s is the synchronous frequency, and n.sub.s (rpm)=60.Math.f.sub.s (Hz). Therefore, the rotor speed, n.sub.r, or rotor frequency, f.sub.r, of the axle may be ascertained by receiving data indicative of the slip, s, and the synchronous speed, n.sub.s, or frequency, f.sub.s. The synchronous speed, n.sub.s, may be expressed as:
n.sub.s=(120.Math.f.sub.1)/p,
where f.sub.1 is the fundamental electrical frequency, and p is the number of motor poles. Similarly, the synchronous frequency, f.sub.s, may be expressed as:
f.sub.s=2.Math.f.sub.1/p.
Therefore, the rotor speed, n.sub.r, or rotor frequency, f.sub.r, of the axle may be ascertained by receiving data indicative of the slip, s, the fundamental electrical frequency, f.sub.1, and the number of poles of the motor, p.

(23) The number of motor poles, p, can be deduced from the nominal speed n.sub.nom of the motor at the nominal frequency f.sub.nom, (which is usually 50 Hz or 60 Hz) by the expression:
p=INTEGER((120.Math.f.sub.nom)/n.sub.nom),
where the INTEGER function truncates the result.

(24) The nominal speed n.sub.nom and nominal frequency f.sub.nom, are quantities defined on the motor nameplate (also known as the motor type shield or motor rating plate) found on most, if not all, motors (typically listed simply as speed and frequency).

(25) If the load dependent slip, s, is unknown, it may be calculated using the expression:
s=2.Math.f/f.sub.1,
where f is the frequency shift due to rotor slip. f (and f.sub.1, if unknown) can be obtained by analysis of the frequency spectrum of the motor current. FIG. 3 shows an example of a motor current frequency spectrum where the fundamental frequency, f.sub.1, is 60 Hz. On either side of the fundamental frequency, f.sub.1, there is a notable sideband frequency that is shifted relative to the fundamental frequency, f.sub.1, by f. Thus, analysis of the frequency spectrum can provide the value of f. In order to successfully derive f from a measured frequency spectrum, the resolution of the frequency spectrum should be sufficiently great. In certain preferable embodiments, a resolution of 0.1 to 0.5 Hz is sufficient. Also, given that the sideband signals may be significantly lower than the signal at the fundamental frequency, f.sub.1, the frequency spectrum is preferably obtained with a sufficiently great dynamic range to permit the sidebands to be accurately identified. In the example shown in FIG. 3, the sidebands are approximately 60 dB lower than the fundamental component, f.sub.1.

(26) Considering the above expression for the rotor frequency, f.sub.r, and substituting the expression for the synchronous frequency, f.sub.s, provides:
f.sub.r=(1s).Math.2.Math.f.sub.1/p.

(27) Given that the angular rotor frequency, .sub.r, is defined as .sub.r=2.Math..Math.f.sub.r, it may be expressed (in rad/s) as:
.sub.r=2.Math..Math.(1s).Math.2.Math.f.sub.1/p.

(28) Therefore, simply by analyzing the motor current frequency spectrum and obtaining data indicative of the nominal speed n.sub.nom and nominal frequency f.sub.nom of the motor (which may be readily obtained from the motor nameplate), the angular rotor frequency, .sub.r, may be determined.

(29) Once the electromechanical power, P.sub.em, and the angular rotor frequency, .sub.r, is determined, the electromechanical axle torque, T.sub.em, of the motor can be calculated using the expression:
T.sub.em=P.sub.em/.sub.r.

(30) Certain methods according to embodiments of the present disclosure are particularly suited to calculating electromechanical axle torque in asynchronous motors. The skilled person will appreciate that for synchronous motors, the effects of slip will not apply.

(31) Certain embodiments of the present disclosure comprise apparatus that is configured to carry out the method. Returning to FIG. 2, the apparatus 40 may be any suitable apparatus for obtaining the required data and determining electromechanical axle torque, T.sub.em, in accordance with methods according to embodiments of the present disclosure. The apparatus 40 may include any suitable means, sensor, or system for measuring the instantaneous voltage and current on each connected phase of the motor 50. The apparatus 40 may include a control means in the form of a control unit that comprises one or more processors for executing computer software instructions that may be arranged, when executed, to perform a method according to an embodiment of the present disclosure. The software instructions may be stored in a non-transitory computer-readable medium, such as a memory of the control unit. The control unit may be communicably coupled to the means, sensor or system (if present) for measuring the instantaneous voltage and current on each connected phase of the motor 50. The apparatus 40 may include a display for displaying information, which may include, for example, the calculated electromechanical axle torque, T.sub.em. The apparatus may include one or more transmitters and/or receivers for transmitting and/or receiving data. For example, the apparatus 40 may be configured to transmit data indicative of the electromechanical axle torque, T.sub.em.

(32) Throughout the description and claims of this specification, the words comprise and contain and variations of them mean including but not limited to, and they are not intended to (and do not) exclude other moieties, additives, components, integers or steps. Throughout the description and claims of this specification, the singular encompasses the plural unless the context otherwise requires. In particular, where the indefinite article is used, the specification is to be understood as contemplating plurality as well as singularity, unless the context requires otherwise.

(33) Features, integers, characteristics, compounds, chemical moieties or groups described in conjunction with a particular aspect, embodiment or example disclosed herein are to be understood to be applicable to any other aspect, embodiment or example described herein unless incompatible therewith. All of the features disclosed in this specification (including any accompanying claims, abstract and drawings), and/or all of the steps of any method or process so disclosed, may be combined in any combination, except combinations where at least some of such features and/or steps are mutually exclusive. The invention is not restricted to the details of any foregoing embodiments. The invention extends to any novel one, or any novel combination, of the features disclosed in this specification (including any accompanying claims, abstract and drawings), or to any novel one, or any novel combination, of the steps of any method or process so disclosed.

(34) The reader's attention is directed to all papers and documents which are filed concurrently with or previous to this specification in connection with this application and which are open to public inspection with this specification, and the contents of all such papers and documents are incorporated herein by reference.