PatentDe  


Dokumentenidentifikation EP0924852 29.03.2007
EP-Veröffentlichungsnummer 0000924852
Titel Automatische Kalibrierung der Parameter eines feldorientierten Aufzugmotorantriebs unter Verwendung von Messungen im Motorstillstand
Anmelder Otis Elevator Co., Farmington, Conn., US
Erfinder Colby, Roy Stephen, Raleigh, North Carolina 27615, US;
Vecchiotti, Alberto, Middletown, Connecticut 06457, US;
Lamontagne, Leslie M., Prospect, Connecticut 06712, US
Vertreter Klunker, Schmitt-Nilson, Hirsch, 80797 München
DE-Aktenzeichen 69837073
Vertragsstaaten DE, FR, GB
Sprache des Dokument EN
EP-Anmeldetag 18.12.1998
EP-Aktenzeichen 983104100
EP-Offenlegungsdatum 23.06.1999
EP date of grant 14.02.2007
Veröffentlichungstag im Patentblatt 29.03.2007
IPC-Hauptklasse H02P 21/00(2006.01)A, F, I, 20070116, B, H, EP

Beschreibung[en]

Co-pending EP Patent Applications, (EP 98310433.2, EP 98310406.8, EP 98310409.2, EP 98310470.4), filed contemporaneously herewith, contain subject matter related to that disclosed herein.

This invention relates to calibrating a field-oriented motor drive and more particularly to automatic calibration of field-oriented (or vector-controlled) drive parameters for an elevator motor drive.

It is known in the art of field-oriented (or vector-controlled) motor drives and motor speed controls that such drives and controls require knowledge of the motor parameters such as the rotor time constant (&tgr;R), torque constant (KT*), and rated magnetizing current IdRATED.

One technique used to determine these motor parameters is to analyze the motor in an engineering laboratory using a dynamometer and expensive test equipment and performing time consuming, highly technical tasks by a skilled technician or engineer. However, in modernization or retrofit applications, where a new drive replaces an older drive in an existing elevator system, it is not convenient or cost effective to remove the motor from the elevator system for evaluation of the motor.

Also, it is desirable to determine the motor parameters based on measurements made on the motor at standstill. While techniques exist to determine the motor parameters while the motor is running (at no-load and under-load), it is not always practical to run such tests in an elevator application. In particular, a no-load test is not practical because this would require unroping the elevator or disconnecting the motor from the gearbox. Further, performing a test under load, i.e., with the motor connected to an elevator, is not practical because it is necessary to have approximately correct motor parameters to start the elevator moving to obtain measurements of the motor running under-load. Also, it is desirable that the technique for determining such motor parameters be wholly contained within the drive control itself, so that field commissioning of retrofit or modernization drives may be performed by installers and service personnel without requiring specific motor/drive tuning skills.

Objects of the invention include provision of automatic, on-site, calibration of motor parameters for field-oriented drives and/or controls for elevators, which does not require removal or uncoupling of the motor from the elevator system, and which uses only standstill measurements of the motor.

According to the present invention, a method for calculating at least one parameter of an elevator motor operated by field-oriented control, using standstill measurements of the motor, the motor having a motor impedance ZM, a rotor impedance ZR, and a transient inductance L&sgr;, includes: a) providing a sinusoidal torque current reference signal at a high frequency (FHIGH) high enough such that the motor impedance (ZM) is dominated by the transient inductance (L&sgr;); b) measuring a feedback torque current (Iq) and a feedback torque voltage (Vq); c) calculating the transient inductance (L&sgr;) at the high frequency (FHIGH), by the equation: L&sgr; = Imag (ZM) @ FHIGH /(2&pgr; FHIGH), where ZM = Vq/Iq; d) providing a sinusoidal torque current reference signal having a variable input frequency; e) measuring the feedback torque current (Iq) and the feedback torque voltage (Vq); f) calculating an imaginary part of the rotor impedance Imag(ZR) as follows: Imag(ZR) = Imag(ZM) - &ohgr;L&sgr;, where &ohgr; is the input frequency, and ZM = Vq/Iq; g) varying the input frequency and performing steps (d) - (f) to obtain the frequency (FPEAK) at which the maximum value of Imag(ZR) occurs; and h) calculating a rotor time constant (&tgr;R) based on FPEAK.

According further to the present invention, &tgr;R is calculated as follows: &tgr;R = 1/(2&pgr;FPEAK). According still further to the present invention, an additional step of calculating a magnetizing inductance Lϕ is performed as follows: L&phgr; = 2Imag(ZR)/&ohgr; @ &ohgr; = 1/&tgr;R.

According still further to the present invention, the following additional steps are performed: i) calculating a motor torque constant (KT*) using the equation: KT* = (3/2)(P/2) L&phgr; Id*, where: P= number of motor poles, Id* = Vph_RATED/(&ohgr;R_RATED x L&phgr;), and Vph_RATED = rated motor voltage, &ohgr;R_RATED = rated motor speed; j) calculating a motor voltage (Vm*) using the equation: Vm* = (Vd*2 + Vq*2)S, where: Vd* = &ohgr;EL&sgr;Iq*, Vq* = &ohgr;ELsId*, Iq* = T_RATED/KT*, T_RATED is the rated motor torque, &ohgr;E* = &ohgr;R_RATED + &ohgr;S*, and &ohgr;S*= (1/&tgr;R)(Iq*/Id*); k) calculating a ratio of Vph_RATED to Vm*; and 1) varying Id* and performing steps (h) - (k) until the ratio is within a predetermined tolerance of 1.

The present invention represents a significant improvement over the prior art by providing automatic calibration of a field-oriented (or vector-controlled) induction motor controller for an elevator system based on standstill measurements of the induction motor. The invention provides motor parameters such as the rotor time constant (&tgr;R), torque constant (KT*), and rated magnetizing current IdRATED, without disconnecting the motor from the elevator system or from the gearbox. The invention also computes transient inductance L&sgr;, magnetizing inductance L&phgr;, stator resistance Rs, and rated torque current IqRATED, as needed by a given application. Further, the invention does not require a specially trained engineer with special test equipment to tune the motor/drive system. Thus, the invention greatly reduces cost associated with tuning the motor drive when new motors drives are retrofit into job sites. Accordingly, automatic calibration of motor parameters at the field site saves both time and money. As a result, the present invention makes it more attractive for building owners to upgrade their elevator systems to modern controls, which are currently economically impractical due to the high cost of determining parameters of older motors found in modernization job sites.

Various embodiments of the present invention will now be described by way of example, and with reference to the accompanying drawings in which:

  • Fig. 1 is a schematic block diagram of a motor controller with auto-calibration logic, in accordance with the present invention.
  • Fig. 2 is a schematic block diagram of a current regulator/motor drive, in accordance with the present invention.
  • Fig. 3 is a schematic drawing of an equivalent circuit model of an induction motor controlled by field orientation, in accordance with the present invention.
  • Fig. 4 is a simplified schematic diagram of the equivalent circuit of Fig. 3, in accordance with the present invention.
  • Fig. 5 is a logic flow diagram of the auto-calibration logic of Fig. 1, in accordance with the present invention.
  • Fig. 6 is a graph of the imaginary part of the rotor impedance and of the motor impedance versus frequency, in accordance with the present invention.
  • Fig. 7 is a logic flow diagram of a portion of the auto-calibration logic of Fig. 5, in accordance with the present invention.

Referring to Fig. 1, a portion of an elevator motor controller 7 shown to the left of the line 9, includes a field-oriented (or vector-based) motor control that has two control loops each corresponding to a different control axis, a d-axis relating to motor magnetization, and a q-axis relating to torque. The d-axis loop has a d-axis current reference input signal IdREF provided on a line 14. IdREF is set to a predetermined constant value so as to provide appropriate magnetic flux in the motor based on motor magnetization curves, e.g., IdRATED or INO-LOAD, discussed more hereinafter. The IdREF signal is fed to a field-oriented current regulator/motor drive circuit 20, described more hereinafter with Fig. 2.

The q-axis current loop has a first q-axis current reference input signal IqREF1 on a line 15 is fed to one input of a switch 19. IqREF1 is provided by other logic (not shown), such as speed loop compensation logic (not shown), which closes a motor speed control loop which provides the q-axis current reference signal to the controller when it is not in auto-calibration.

The other input to the switch 19 is a second q-axis current reference input signal IqREF2 on a line 17. The output of the switch 19 is the q-axis current loop reference signal I qREF on a line 18, which is set equal to I qREF1 or I qREF2 based on the state of the MODE signal provided to the switch 19 on the line 13. The IqREF signal is fed to the field-oriented current regulator/motor drive circuit 20, described more hereinafter with Fig. 2.

Two examples of three phase AC induction motors used with the present invention are, Model LUGA-225LB-04A, by Loher, having a rated power of 45KW, rated voltage of 355 volts, rated speed of 1480, and rated frequency of 50 Hz, in a geared configuration; and Model 156MST, by Tatung (of Taiwan), having a rated power of 40KW, rated voltage of 500 volts, rated speed of 251, and rated frequency of 16.7 Hz, in a gearless configuration. Other motors having other rated parameters may be used if desired.

The circuit 20 provides three phase voltage signals Vx,Vy,Vz on lines 22 to a motor 24, e.g., a three phase induction motor. The motor 24 is connected by a mechanical linkage 26, e.g, a shaft and/or a gearbox, to a sheave 28. A rope or cable 30 is wrapped around the sheave 28 and has one end connected to an elevator car 32 and the other end connected to a counterweight 34. The weight of the counterweight is typically equal to the weight of an empty car plus 40-50% of the max load in the car.

Other elevator system configurations, and with or without a counterweight, with or without a gearbox, may be used if desired to convert the output torque of the motor 24 to movement of the elevator cab 32, such as dual lift (where two elevator cars are connected to a single rope, the cars move in opposite directions and each car provides a counterweight for the other car), drum machine (where the rope is wrapped around a drum driven by a motor), etc.

A brake 37, e.g., an electro-magnetic actuated disk brake, is disposed on the shaft 26, and is driven by an electrical brake command signal BRKCMD on a line 38 from the circuit 20. The brake 37, when activated or "dropped", clamps onto the shaft 26 and prevents the motor shaft 26 from turning, i.e., locks the rotor, and thus prevents the sheave 28 from moving.

Referring to Fig. 2, it is known in the art of field-oriented motor control that such control uses current (Id,Iq) and voltage (Vd,Vq) parameters corresponding to the d and q axes, respectively. Using field orientation, the motor magnetic field (or flux) will be controlled by Id and the motor torque will be controlled by Iq, as is known. In particular, the field-oriented current regulator/motor drive 20 of Fig. 1 comprises two current control loops, one for the d-axis current Id and one for q-axis current Iq. The Id loop receives the IdREF signal on the line 14 to a positive input to a summer 102. A measured or feedback d-axis current signal Id on a line 104 is fed to a negative input to the summer 102. The output of the summer 102 is an error signal IdERR on a line 106 which is fed to control compensation logic 108, such as proportional-plus-integral (P-I) current loop control. Other current loop control compensation may be used if desired. The logic 108 provides a d-axis voltage command signal VdCMD on a line 110.

For the q-axis, the Iq loop receives the IqREF signal on the line 18 to a positive input to a summer 114. A measured or feedback q-axis current signal Iq on a line 116 is fed to a negative input to the summer 114. The output of the summer 114 is an error signal IqERR on a line 118 which is fed to control compensation logic 120, e.g., proportional-plus-integral (P-I) logic similar to the logic 108. The output of the logic 120 is a q-axis voltage command signal VqCMD on a line 122. Other control compensation, e.g., proportional, lead-lag, etc., may be used for the logics 108,120. The form of compensation used is not critical to the present invention.

The voltage commands VdCMD and VqCMD are fed to known field-oriented to three-phase conversion logic 124 which converts the d-axis and q-axis voltage commands to three phase voltage commands VXCMD, VYCMD, VZCMD on lines 126. The phase voltage commands VXCMD, VYCMD, VZCMD are fed to a known three phase drive circuit (or inverter) 128 which provides three phase voltages VX,VY,VZ on lines 130, 132, 134, respectively (collectively, the lines 22), to drive the motor 24.

Within the drive circuit 128 (details not shown), each of the voltage commands VXCMD, VYCMD, VZCMD on lines 126 are converted to percent duty cycle commands indicative of the corresponding input voltage level. The percent duty cycle is converted into a pulse-width-modulated drive signal which drives power transistors to provide the pulse-width-modulated, variable frequency, three phase voltages VX,VY,YZ on lines 130, 132, 134, respectively. The conversions within the drive 128 are performed using electronic components and/or software well known in the art of motor drive circuits. Any other type of drive circuit that receives input voltage commands and provides output phase voltages may be used, and the phase voltages need not be pulse-width modulated.

Phase currents IX, IY, IZ associated with the voltages VX,VY,VZ, respectively, are measured by known current sensors 136, 138, 140, e.g., closed-loop Hall-effect current sensors (such as LEMS), respectively, and are provided on lines 132, 134, 136, respectively. The phase currents IX, IY, IZ are fed to known three phase to field oriented conversion logic 142, which provides a known conversion from phase currents to d-axis and q-axis currents on the lines 104, 116 which are fed to the summers 102,114, respectively.

The converters 124,150 provide known conversions between vector (d and q axis) parameters and per-phase parameters, such as that described in D. Novotny, et al, "Vector Control and Dynamics of AC Drives", Oxford University Press, 1996, Ch 5, pp 203-251. The converters 124,15 may likely implement such conversions in software using a microprocessor or the like.

It is known in the art of field oriented drives that the value of the rotor time constant &tgr;R of the motor being controlled is required to perform the conversion to and from the field oriented d and q axes. In particular, &tgr;R is used to establish the correct slip frequency &ohgr;S to achieve field orientation. The value of the rotor time constant &tgr;R is provided to the two converters 124, 150 on a line 144.

The motor drive logic 111 also includes a brake drive circuit 145 which receives an input signal BRK on a line 146 and provides a BRKCMD signal on the line 38.

Referring to Fig. 1, the present invention comprises auto-calibration logic 48 which automatically computes the motor parameters &tgr;R, KT*, IdRATED, and provides &tgr;R on the line 144 to the circuit 20, provides IdRATED as IdREF on the line 14 to the circuit 20, and provides KT* on a line 160 to speed loop compensation logic (not shown). The logic 48 also computes other motor parameters such as transient inductance L&sgr;, magnetizing inductance L&phgr;, stator resistance Rs (or R1) and rated torque current IqRATED. The logic 48 receives Vq and Iq from the circuit 20. The logic 48 also provides the current reference signal IqREF2 to the switch 19, and provides IdREF to the circuit 20 on the line 14.

The logic 48 comprises known electronic components, which may include a microprocessor, interface circuitry, memory, software, and/or firmware, capable of performing the functions described herein.

The logic 48 also provides the MODE1 signal on the line 13 to the switch 19. The MODE1 flag causes the current reference signal IqREF2 from the calibration logic 48 to be fed to the logic 20. The logic 48 also provides a break request signal BRK on the line 146 to the circuit 20.

The calculation logic 48 also communicates with a service tool 80 over a serial link 82. The service tool 80 includes a display 84 and a keypad (or keyboard) 86 for entering data into the service tool 80 and over the link 82 to the controller 7. In particular, the logic 48 receives a start command over the link 82 from the service tool 80, which controls when auto-calibration is started. The logic 66 also provides a DONE signal and a FAULT signal to the service tool 80 over the link 82. The DONE signal indicates when auto-calibration has completed without faults and the FAULT signal indicates when a fault has been detected during auto-calibration.

Referring to Fig. 3, a known equivalent circuit 90 of an induction motor is similar to that described in "Vector Control and Dynamics of AC Drives", Novotny and Lipo, Oxford 1996, Chapter 5. Fig. 3 is a per-phase equivalent circuit for AC steady state operation where the current I 1 and voltage V 1 are phasor quantities. The circuit 90 comprises a resistor Rs in series with an equivalent "transient" inductor L&sgr; in series with a rotor impedance ZR which comprises a "magnetizing" inductance L&phgr; in parallel with an equivalent resistance R2/S. Where:

  • Rs (or R1) = stator winding resistance
  • Ls = stator winding inductance
  • Lr = rotor winding inductance
  • Lm = mutual inductance
  • Rr = the rotor winding resistance
  • L&sgr; = Ls - Lm2/Lr = transient inductance
  • L&phgr; = Lm2/Lr = magnetizing inductance
  • &ohgr;E = electrical frequency of the input current I1
  • &ohgr;R = motor output rotational speed in radians per second referred to an electrical reference frame
  • S = slip = (&ohgr;E - &ohgr;R)/&ohgr;E
  • &ohgr;S= slip frequency =&ohgr;E - &ohgr;R= (1/&tgr;R)(Iq/Id)
where &tgr;R= rotor time constant Iq= q-axis (or torque) current, and Id = d-axis (or magnetizing) current R 2 = ( Lm 2 / Lr 2 ) * Rr

Also, the rotor time constant &tgr;R and motor torque constant KT* are related to the parameters of the circuit 90 as follows: &tgr; R = Lr / Rr = L&phgr; / R 2 K T * = 3 / 2 P / 2 L&phgr; Id = torque / current

where P = number of poles.

Referring to Fig. 4, the circuit 92 is an equivalent to the circuit 90 of Fig. 3 with the rotor impedance ZR transformed into an equivalent series circuit impedance having a real part Real(ZR) and an imaginary part Imag(ZR) equal to &ohgr;Lx. The equivalent circuit 92 with the transformation of ZR is useful for determining the rotor time constant &tgr;R (discussed more hereinafter).

Referring now to Fig. 5, a top level flow chart for the logic 48 begins at a step 200 which determines whether a start command has been received from the service tool 80. If it has not, the logic exits. If a start command has been received, a step 202 requests and receives motor parameters from the service tool 80 over the link 82 (Fig. 1), which are entered by service personnel. The motor parameters received are: the rated motor shaft power in watts (PWR_RATED); the rated motor speed in rpm (RPM_RATED); the rated rms line-to-line voltage in volts (VLL_RATED); the rated frequency in hertz (HZ_RATED); and number of poles (POLES), all of which may be obtained from the motor nameplate data.

Then, a step 203 sets MODE1 = 1, BRK = 1 to cause the brake 37 (Fig.1) to lock the rotor, and IdREF2 = 0 amps. For each of the tests described herein, the rotor remains locked (rotor speed &ohgr;R = 0) and IdREF2 = 0 amps. When &ohgr;R = 0 and IdREF2 = 0, the slip S=1, and the motor current I1 is equal to the q-axis current Iq and the motor voltage V1 is equal to the q-axis voltage Vq. When Iq=0, the motor is operated in single phase operation, in accordance with the circuits of Figs. 3,4.

Next, a step 204 measures the transient inductance L&sgr;, by providing a sinusoidal current signal into the q-axis of the reference current IqREF2 on the line 17 (Fig. 1) at a frequency FHIGH high enough such that the motor impedance will be dominated by the transient inductance L&sgr;, e.g., 31.25 Hertz. Other frequencies may be used if desired, e.g., greater than 30 Hz. The sinewave input signal is generated digitally by a signal processor, such as a digital signal processor, e.g., a Motorola DSP 56002 processor, with an update (or sample) rate of 5KHz. Other hardware and/or software techniques or update rates may be used to generate the sinusoidal input signals.

The step 204 reads the q-axis feedback current Iq and the q-axis output voltage Vq (equal to the motor current I1 and motor voltage V1, respectively, as discussed hereinbefore). Next, the step 204 uses the aforementioned digital signal processor to perform a Discrete Fourier Transform (DFT) of Iq and Vq to determine the first harmonic Fourier coefficients. The fundamental or first harmonic component of a measured signal from a DFT is A sin(&ohgr;t) + B cos(&ohgr;t), where &ohgr;=2&pgr;f is the input frequency (in rad/sec). The first harmonic is used to calculate the impedance primarily so that non-linearities in the system do not distort the calculation.

To compute a DFT, as is known, standard sine and cosine waves of unit amplitude at the test frequency are generated within the logic 48. The measured signal (Iq,Vq) is multiplied by the standard sinewave and the product is integrated over one period of the excitation to yield the Fourier series coefficient A of the signal. Multiplying the signal by the standard cosine and integrating yields the B coefficient. We have found that integrating over 15 periods of the input signal it is sufficient to filter out any transients in the system response. Other numbers of periods may be used if desired. Also, for any DFT discussed herein, other types of Fourier transforms may be used if desired, e.g., a Fast Fourier Transform (FFT), etc., provided the first harmonic of the desired signal is obtained. Further, instead of a Fourier transform, any other filtering or spectrum analysis technique may be used for determining the first harmonic of the desired signals.

Then, the step 204 computes the motor impedance ZM by calculating the ratio of voltage to current (V1/I1 = Vq/Iq) using the first harmonic components of voltage and current computed above. The step 204 then computes the real and imaginary parts of ZM from the Fourier coefficients. The imaginary part of the motor impedance ZM at FHIGH Hertz is dominated by the transient inductance term &ohgr;L&sgr;. Thus, the transient inductance L&sgr; is the transient reactance (or imaginary part of ZM) with the input frequency equal to FHIGH Hz, divided by the frequency &ohgr; in radians/sec (2&pgr; FHIGH), or : L&sgr; = Imag Z M @ F HIGH Hz / 2 &pgr; F HIGH

Next, an optional step 206 measures the total resistance of the circuit impedance (RTOT = RS+R2), i.e., the sum of the stator and rotor resistances, as the real part of the motor impedance ZM determined in step 204. Thus: R TOT = Real Zm @ F HIGH In particular, at the relatively high frequency FHIGH used in step 204, the inductance Lϕ in the circuit 90 is large and the real part of ZM will be equal to RTOT. The value of RTOT is saved for later use to calculate Rs (see step 212).

Next, a step 208 measures the rotor time constant &tgr;R as follows. The step 208 produces a progression of low frequency sinusoidal input q-axis reference currents IqREF2 from 0.1 to 8.0 Hertz in increments defined by a search algorithm, discussed hereinafter. The sinewave input signal is generated digitally as discussed hereinbefore with step 204. At each frequency, the motor current Iq and voltage Vq (equal to the motor current I1 and voltage V1, respectively, as discussed hereinbefore) are measured and a DFT of the current I1 and the motor voltage V1 signals are computed separately. The fundamental or first harmonic Fourier coefficients are obtained as discussed hereinbefore with step 204.

Then, the step 208 computes the motor impedance ZM at each frequency by calculating the ratio of voltage to current (V1/I1). The step 208 then calculates the real and imaginary parts of the ZM from the Fourier coefficients. Then, the step 208 calculates the imaginary part of the rotor impedance Imag(ZR) = &ohgr;Lx by subtracting the transient reactance (&ohgr;L&sgr;) from the imaginary part of the motor impedance ZM, where L&sgr; was previously calculated in step 204 and &ohgr; is the input frequency, as follows: Imag Z R = &ohgr;Lx = Imag Z M - &ohgr;L&sgr;

Referring to Fig. 6, a curve 250 shows the imaginary part of the motor impedance Imag(ZM) = &ohgr; (L&sgr; + Lx) and a curve 252 shows the imaginary part of the rotor impedance Imag(ZR) = &ohgr;Lx. The frequency &ohgr; (radians/second) at which the maximum 254 of the curve 252 occurs is the inverse of the rotor time constant, i.e., &ohgr; = 1/&tgr;R. A known search algorithm, e.g., a "golden section line search" algorithm, varies the input frequency and determines the frequency Fpeak at which the maximum value &ohgr;Lx occurs. The type of search algorithm used is not critical to the present invention, and any search algorithm that varies an input parameter and determines the maximum value of an output parameter may be used. The rotor time constant &tgr;R is then calculated as follows: &tgr; R = 1 / &ohgr;paek = 1 / 2 &pgr;Fpeak

Next a step 210 calculates the magnetizing inductance L&phgr;. In particular, at the frequency of rotor time constant (&ohgr; = 1/&tgr;R), which is also the break frequency of the motor transfer function, the real and imaginary parts of the rotor impedance ZR are equal to each other, i.e., &ohgr;LX = RX. Also, at this same frequency, it can be shown (below) that the &ohgr;LX is also equal to S&ohgr;L&phgr; (the magnetizing reactance). In particular, the rotor impedance ZR is equal to j&ohgr;L&phgr; in parallel with R2, as shown below: Z R = j&ohgr;L &phgr; R 2 / R 2 + j&ohgr;L &phgr; Multiplying the numerator and denominator by the complex conjugate of the denominator (R2 - j&ohgr;L&phgr;), gives Z R = &ohgr; 2 L &phgr; 2 R 2 / R 2 2 + &ohgr; 2 L &phgr; 2 + j&ohgr;L &phgr; R 2 2 / R 2 2 + &ohgr; 2 L &phgr; 2 Which has the form of a series combination of impedances, or a real and imaginary part as indicated below: Z R = Rx + j&ohgr;Lx Z R = Real + j Imaginary At the peak 254 of the curve 252 of Imag(ZR), the real and imaginary parts are equal, which gives: &ohgr; 2 L &phgr; 2 R 2 / R 2 2 + &ohgr; 2 L &phgr; 2 = &ohgr;L &phgr; R 2 2 / R 2 2 + &ohgr; 2 L &phgr; 2 Simplifying Eq. 2, gives: &ohgr;L &phgr; = R 2 Substituting R2 = &ohgr;L&phgr; into the Imaginary part of ZR, and setting equal to &ohgr;Lx, gives: Imag Z R = &ohgr;L &phgr; &ohgr; 2 L &phgr; 2 / &ohgr; 2 L &phgr; 2 + &ohgr; 2 L &phgr; 2 = &ohgr;Lx Simplifying Eq. 3, gives: &ohgr;Lx = &ohgr;L &phgr; / 2 Thus, the magnetizing inductance Lϕ is calculated as follows: L&phgr; = 2 Imag Z R / &ohgr; @ &ohgr; = 1 / &tgr; R

Next, an optional step 212 calculates the stator resistance RS by first calculating the value of R2. It can be shown (below) that the real part of the rotor impedance Real(ZR) at &ohgr; = 1/&tgr;R is equal to R2/2. In particular, the real part of Eq. 1 is: Real Z R = Rx = L&phgr; 2 R 2 / R 2 2 + &ohgr; 2 L&phgr; 2 Substituting R2=&ohgr;L&phgr;, and simplifying, gives: Rx = R 2 / 2 Thus, R 2 = 2 Real Z R @ &ohgr; = 1 / &tgr; R Alternatively, R2 may be calculated using the equation: R 2 = L&phgr; / &tgr; R

where L&phgr; and &tgr;R were previously calculated in steps 204,208, respectively. In either case, the stator resistance RS is then determined by subtracting R2 from the total resistance (RTOT= RS + R2) calculated in step 206. Thus, R S = R TOT - R 2 If the value for RS for the motor is already known, e.g., from the data sheet, it may be provided to the control over the link 82 and then RS may also be range checked in the step 212 to ensure it is within a predetermined percentage of the expected value. If RS is not within the desired range, the step 212 sets a fault flag FAULT 1 = 1. Alternatively, the value of RS may be calculated and provided to the service tool to help service personnel determine the type of motor installed in the system.

Next, a step 214 uses L&phgr;, &tgr;R, and the input parameters PWR_RATED, RPM_RATED, VLL_RATED, HZ_RATED, and POLES, obtained in the step 202, to simulate motor parameters and to iterate and calculate the rated magnetizing current IdRATED, and the torque constant KT*, as shown in Fig. 7.

Referring to Fig. 7, the simulated motor parameters are indicated by an asterisk (*) to avoid confusion with actual measured motor parameters discussed hereinbefore. In particular, a step 300 calculates the rated rotational speed of the motor referred to the electrical reference frame &ohgr;R_RATED. Next, a step 302 converts rated line-to-line voltage (VLL_RATED) to rated line-to-neutral voltage (or per-phase voltage) Vph_RATED. Next, a step 303 calculates the rated torque T_RATED based on rated power and rated RPM. Then, a step 304 calculates the stator inductance Ls as the sum of the transient inductance L&sgr; and the magnetizing L&phgr;. Next, a step 306 calculates an initial value for simulated d-axis current Id* based on a first order approximation of Id using rated voltage and speed. Next, a step 308 sets a variable COUNT equal to zero.

Next, a series of steps 310 - 322 calculates KT* and a simulated motor voltage VM* using various simulated motor parameters based on the value of L&phgr; calculated in step 210 (Fig. 5), the parameters calculated in steps 300-308 above, and using known relationships for a field-oriented motor controller, some of which are discussed hereinbefore. In particular, a step 310 calculates the torque constant KT* based on L&phgr; calculated in step 210 (Fig. 5) and the current value of the magnetizing current Id*. Next, a step 312 calculates the torque current Iq*. Next, a step 314 calculates a simulated slip frequency &ohgr;S* which is used in a next step 316 to calculate a simulated electrical current frequency &ohgr;E* which is equal to the rotational frequency (or speed) of the motor (emulated as the rated speed) &ohgr;R_RATED plus the slip frequency &ohgr;S*.

Next, a step 318 calculates a simulated q-axis output voltage Vq* based on the magnetizing current Id* and a step 320 calculates a simulated d-axis output voltage Vd* based on the torque current Iq*. Then, a step 322 calculates a simulated vector sum total motor voltage Vm* equal to the square root of the sum of the squares of the d-axis and q-axis output voltages Vd*,Vq*, respectively.

Next, a step 324 calculates a Ratio parameter equal to the ratio of the rated phase voltage Vph_RATED to the simulated per-phase motor voltage Vm*. The logic iterates until the Ratio goes to 1 within the desired tolerance, e.g., 0.001. When the ratio equals 1 the value of Id* will produce the rated voltage at the rated RPM and rated torque.

Next, a step 326 calculates a next value for Id* equal to the value of the Ratio times the current value of Id*. Next, a step 328 checks whether Ratio is within a predetermined tolerance of 1, e.g., 0.001. If it is not within the desired tolerance, a step 330 checks whether COUNT is greater than or equal to ten (i.e., whether the loop has iterated at least ten times). If the loop has iterated at least ten times, a FAULT flag is set equal to 1 at a step 332 and output to the service tool 80 over the link 82 (Fig. 1) and the logic is exited. If it has iterated less than ten times, a step 334 increments COUNT by one and the logic 214 goes the step 310 to iterate again.

If Ratio is within the desired tolerance in step 328, the logic is deemed to have converged and at convergence the values of Id* and Iq* are equal to the rated d-axis current IdRATED and the rated q-axis current IqRATED, respectively. Accordingly, a step 340 sets the d-axis current reference IdREF equal to Id* which is equal to IdRATED and a step 344 sets IqRATED equal to Iq*. Then the logic 214 exits and returns to the logic 48 of Fig. 5.

Referring to Fig. 5, next, a step 216 determines whether an error has been detected in any of the above steps 202 - 214 (i.e., if FAULT 1 = 1). If a fault has been detected, a step 218 sets FAULT = 1 which is sent to the service tool 80 (Fig. 1) over the serial link 82 and a step 220 sets MODE1 = 0, BRK=0, and the logic exits. If a fault has not occurred, a step 222 sets the DONE flag equal to 1 which is transmitted via the serial link 82 to the service tool 80. Next, some or all of the motor parameters &tgr;R, KT*, IdRATED, L&sgr;, L&phgr;, Rs, and IqRATED are transmitted via the serial link 82 to the service tool 80 in a step 224. The service tool 80 displays the parameters for use by the service personnel. Next, the step 220 sets MODE1=0, BRK=0, and the logic 48 exits.

Although the invention has been described and illustrated with respect to exemplary embodiments thereof, it should be understood by those skilled in the art that the foregoing, and various other changes, omissions and additions may be made without departing from the scope of the present invention as defined in the claims.


Anspruch[de]
Verfahren zum Berechnen mindestens eines Parameters eines mittels feldorientierter Regelung betriebenen Aufzugmotors unter Verwendung von Stillstandsmessungen des Motors, welcher eine Motorimpedanz ZM, eine Rotorimpedanz ZR und eine Einschwinginduktivität L&sgr; aufweist, umfassend folgende Schritte: a) Bereitstellen eines sinusförmigen Drehmoment-Stromreferenzsignals bei einer hohen Frequenz FHIGH, hoch genug, damit die Motorimpedanz ZM durch die Einschwinginduktivität L&sgr; beherrscht wird; b) Messen eines Rückkopplungsdrehmomentstroms Iq und einer Rückkopplungsdrehmomentspannung Vq; c) Berechnen der Überschwingerinduktivität L&sgr; bei der erwähnten hohen Frequenz FHIGH nach folgender Gleichung: L&sgr; = Imag Z M @ F HIGH / 2 &pgr; F HIGH

wobei ZM = Vq/Iq;
d) Bereitstellen eines sinusförmigen Drehmoment-Stromreferenzsignals mit einer veränderlichen Eingangsfrequenz; e) Messen des Rückkopplungdrehmomentstroms Iq und der Rückkopplungsdrehmomentspannung Vq; f) Berechnen eines Imaginärteils der Rotorimpedanz Imag (ZR) wie folgt: Imag Z R = Imag Z M - &ohgr;L&sgr;

wobei &ohgr; die Eingangsfrequenz ist und ZM = Vq/Iq;
g) Variieren der Eingangsfrequenz und Ausführen der Schritte (d) - (f), um die Frequenz FPEAK zu erhalten, bei der der Maximalwert von Imag (ZR) auftritt; und h) Berechnen einer Rotorzeitkonstanten &tgr;R anhand von FPEAK.
Verfahren nach Anspruch 1, bei dem &tgr;R folgendermaßen berechnet wird: &tgr; R = 1 / 2 &pgr;F PEAK . Verfahren nach Anspruch 1, bei dem der Schritt des Berechnens von Zm das Berechnen einer ersten Harmonischen von Vq und einer ersten Harmonischen Iq beinhaltet. Verfahren nach Anspruch 3, bei dem der Schritt des Berechnens der ersten Harmonischen von Iq und Vq die Berechnung von Fouriertransformationen von Iq und Vq beinhaltet. Verfahren nach Anspruch 2, weiterhin umfassend einen Schritt des Berechnens einer Magnetisierungsinduktivität L&PHgr; wie folgt: L&PHgr; = 2 Imag Z R / &ohgr; @ &ohgr; = 1 / &tgr; R . Verfahren nach Anspruch 5, weiterhin umfassend die Schritte: i) Berechnen einer Motordrehmomentkonstanten KT* mithilfe der Gleichung: K T * = 3 / 2 P / 2 L&phgr; Id *

wobei: P = Anzahl von Motorpolen; Id* = Vph_RATED/(&ohgr;R_RATED x L&PHgr;) und Vph_RATED = Motornennspannung; &ohgr;R_RATED = Motornenndrehzahl;
j) Berechnen einer Motorspannung Vm* unter Verwendung der Gleichung: Vm * = Vd * 2 + Vq * 2

wobei: Vd* = &ohgr;&egr;L&sgr;Iq*; Vq* = &ohgr;&egr;Lsld*; Iq* = T_RATED/KT*; T_RATED das Motor-Nenndrehmoment ist; &ohgr;E* = &ohgr;R_RATED + &ohgr;S*; und &ohgr;S* = (1/&tgr;R) (Iq*/Id*);
k) Berechnen eines Verhältnisses von Vph_RATED zu Vm*; und l) Variieren von Id* und Ausführen der Schritte (h) - (k) so lange, bis jenes Verhältnis innerhalb einer vorbestimmten Toleranz von 1 liegt.
Verfahren nach Anspruch 6, bei dem der Schritt des Variierens Id* das Berechnen des nächsten Werts von Id* als Id* × Verhältnis beinhaltet. Verfahren nach Anspruch 5, weiterhin umfassend die Schritte: m) Berechnen eines Gesamtwiderstands RTOT bei jener hohen Frequenz mittels der Gleichung: R TOT = Real Zm @ F HIGH n) Berechnen eines Statorwiderstands Rs wie folgt: R s = R TOT - R 2

wobei R2 = L&PHgr;/&tgr;R; und
o) Ausführen des Schritts (m) zwischen den Schritten (b) und (d).
Verfahren nach Anspruch 1, bei dem die Schritte (a) - (h) automatisch bei Empfang eines Befehls von einem Servicewerkzeug durchgeführt werden. Verfahren nach Anspruch 6, bei dem die Schritte (a) - (I) automatisch bei Empfang eines Befehls von einem Servicewerkzeug durchgeführt werden. Aufzugregelung (7) umfassend eine Logik (48), die automatisch eine Motorzeitkonstante &tgr;R, eine Drehmomentkonstante KT* und einen Magnetisierungsstrom Id für ein feldorientiertes Motor-/Treibersystem berechnet, mit einer Einrichtung zum Berechnen einer Einschwinginduktivität L&sgr;, einer Einrichtung zum Injizieren einer eine variable Eingangsfrequenz aufweisenden Sinuswelle in den q-Achsen-Referenzstrom IqREF und eine Einrichtung zum Variieren der Eingangsfrequenz, während ein Imaginärteil einer Rotorimpedanz Imag (ZR) berechnet wird, um eine Frequenz FPEAK zu erhalten, bei der der Maximalwert von Imag (ZR) auftritt, zusätzlich mit einer Einrichtung zum Berechnen von &tgr;R aus FPEAK, wobei KT*, Id, und einer Gesamtmotorspannung VM berechnet werden, und einer Einrichtung zum Variieren von Id so lange, bis ein Verhältnis der Nennmotorspannung Vph_RATED zu VM in einer vorbestimmten Toleranz von 1 liegt, und mit einer Einrichtung zum Durchführen der Prozedur bei verriegeltem Rotor.
Anspruch[en]
A method of calculating at least one parameter of an elevator motor operated by field-oriented control, using standstill measurements of the motor, the motor having a motor impedance ZM, a rotor impedance ZR, and a transient inductance L&sgr;, comprising the steps of: a) providing a sinusoidal torque current reference signal at a high frequency FHIGH high enough such that the motor impedance ZM is dominated by the transient inductance L&sgr;; b) measuring a feedback torque current Iq and a feedback torque voltage Vq; c) calculating said transient inductance (L&sgr;) at said high frequency FHIGH, by the equation: Lo = Imag ( Z M ) @ F HIGH / 2 &pgr; F HIGH

where ZM = Vq/Iq;
d) providing a sinusoidal torque current reference signal having a variable input frequency; e) measuring the feedback torque current Iq and the feedback torque voltage Vq; f) calculating an imaginary part of the rotor impedance Imag(ZR) as follows: Imag Z R = Imag Z M - &ohgr;L&sgr;

where &ohgr; is said input frequency, and ZM = Vq/Iq;
g) varying said input frequency and performing steps (d) - (f) to obtain the frequency FPEAK at which the maximum value of Imag(ZR) occurs; and h) calculating a rotor time constant &tgr;R based on FPEAK.
The method of claim 1, wherein said &tgr;R is calculated as follows: &tgr; R = 1 / 2 &pgr; F PEAK . The method of claim 1, wherein said step of calculating Zm comprises computing a first harmonic of Vq and a first harmonic of Iq. The method of claim 3, wherein said step of computing the first harmonic of Iq and Vq, comprises computing Fourier transforms of Iq and of Vq. The method of claim 2, further comprising a step of calculating a magnetizing inductance Lϕ as follows: L&phgr; = 2 Imag Z R / &ohgr; @ &ohgr; = 1 / &tgr; R . The method of claim 5, further comprising the steps of: i) calculating a motor torque constant KT* using the equation: K T * = 3 / 2 P / 2 L&phgr; Id *

where: P= number of motor poles; Id* = Vph_RATED/(&ohgr;R_RATED x L&phgr;); and Vph_RATED = rated motor voltage; &ohgr;R_RATED = rated motor speed;
j) calculating a motor voltage Vm* using the equation: Vm * = Vd * 2 + Vq * 2 S

where : Vd* = &ohgr;EL&sgr;Iq*; Vq* = &ohgr;ELsId*; Iq* = T_RATED/KT*; T_RATED is the rated motor torque; &ohgr;E* = &ohgr;R_RATED + &ohgr;S*; and &ohgr;S* = (1/&tgr;R) (Iq*/Id*);
k) calculating a ratio of Vph_RATED to Vm*; and l) varying Id* and performing steps (h) - (k) until said ratio is within a predetermined tolerance of 1.
The method of claim 6, wherein said step of varying Id* comprises calculating the next value of Id* as Id* x Ratio. The method of claim 5, further comprising the steps of: m) calculating a total resistance RTOT at said high frequency, by the equation: R TOT = Real Zm @ F HIGH n) calculating a stator resistance Rs as follows : R S = R TOT - R 2

where R2 = L&phgr;/&tgr;R ; and
o) performing step (m) between the steps (b) and (d).
The method of claim 1, wherein said steps (a) - (h) are performed automatically upon receiving a command from a service tool. The method of claim 6, wherein said steps (a) - (l) are performed automatically upon receiving a command from a service tool. An elevator controller (7) comprising logic (48) which automatically calculates a motor time constant &tgr;R, a torque constant KT*, and a magnetizing current Id for a field-oriented motor/drive system with means for calculating a transient inductance L&sgr;, means for injecting a variable input frequency sinewave into a q-axis reference current IqREF and means for varying the input frequency while calculating an imaginary part of a rotor impedance Imag(ZR) to obtain a frequency FPEAK at which the maximum value of Imag(ZR) occurs, additionally having means for calculating &tgr;R from FPEAK, wherein KT*, Id, and a total motor voltage VM are calculated and means for varying Id is until a ratio of rated motor voltage Vph_RATED to VM is within a predetermined tolerance of one, means to perform said procedure with the rotor locked.
Anspruch[fr]
Méthode de calcul d'au moins un paramètre d'un moteur d'ascenseur entraîné par une commande à champ orienté, utilisant des mesures à l'arrêt du moteur, le moteur ayant une impédance de moteur ZM, une impédance de rotor ZR, et une inductance transitoire L&sgr;, la méthode comprenant les étapes de : a) fournir un signal de référence de courant de couple sinusoïdal à une fréquence élevée FHIGH, suffisamment élevée pour que l'impédance du moteur ZM soit dominée par l'inductance transitoire L&sgr; ; b) mesurer un courant de couple de rétroaction Iq et une tension de couple de rétroaction Vq ; c) calculer ladite inductance transitoire L&sgr; à ladite fréquence élevée FHIGH, en utilisant l'équation : L&sgr; = Imag ( ZM ) @ F HIGH / 2 &pgr; F HIGH où ZM = Vq/Iq ; d) fournir un signal de référence de courant de couple sinusoïdal ayant une fréquence d'entrée variable ; e) mesurer le courant de couple de rétroaction Iq et la tension de couple de rétroaction Vq ; f) calculer une partie imaginaire de l'impédance du rotor Imag(ZR) de la façon suivante : Imag ( Z R ) = Imag ( Z M ) - &ohgr;L&sgr; où &ohgr; est ladite fréquence d'entrée, et ZM= Vq/Iq ; g) faire varier ladite fréquence d'entrée et accomplir les étapes (d) à (f) de façon à obtenir la fréquence FPEAK à laquelle la valeur maximum de Imag(ZR) est obtenue ; et h) calculer une constante de temps du rotor &tgr;R sur la base de FPEAK. Méthode selon la revendication 1, dans laquelle &tgr;R est calculée de la façon suivante : &tgr; R = 1 / 2 &pgr;F PEAK Méthode selon la revendication 1, dans laquelle ladite étape de calcul de Zm comprend l'étape de calcul d'une première harmonique de Vq et d'une première harmonique de Iq. Méthode selon la revendication 3, dans laquelle ladite étape de calcul de la première harmonique de Iq et de Vq, comprend l'étape de calcul des transformées de Fourier de Iq et de Vq. Méthode selon la revendication 2, comprenant en outre une étape de calcul d'une inductance de magnétisation Lϕ de la façon suivante : L&phgr; = 2 Imag Z R / &ohgr; @ &ohgr; = 1 / &tgr; R . Méthode selon la revendication 5, comprenant en outre les étapes de : i) calculer une constante de couple moteur KT* en utilisant l'équation : K T * = 3 / 2 P / 2 Lϕ Id *

où : P = le nombre de pôles du moteur ; Id* = Vph_RATED (&ohgr;R_RATED × Lϕ) ; et Vph_RATED = la tension nominale du moteur ; &ohgr;R_RATED = la vitesse nominale du moteur ;
j) calculer une tension moteur Vm* en utilisant l'équation : Vm * = Vd * 2 + Vq * 2 1 / 2

où : Vd* = &ohgr;EL&sgr;Iq* ; Vq* = &ohgr;ELSId* ; Iq* = T_RATED/K&tgr;* ; T_RATED est le couple nominal du moteur ; &ohgr;E* = &ohgr;R_RATED + &ohgr;S* ; et &ohgr;s* = (1/&tgr;R) (Iq*/Id*) ;
k) calculer un rapport de Vph_RATED par rapport à Vm* ; et 1) faire varier Id* et accomplir les étapes (h) à (k) jusqu'à ce que ledit rapport se situe à l'intérieur d'une tolérance prédéterminée de 1.
Méthode selon la revendication 6, dans laquelle ladite étape consistant à faire varier Id* comprend l'étape de calcul de la valeur suivante de Id* en tant que Id* × Ratio (le rapport). Méthode selon la revendication 5, comprenant en outre les étapes de : m) calculer une résistance totale RTOT à ladite fréquence élevée, en utilisant l'équation : R TOT = Real Zm @ F HIGH n) calculer une résistance du stator Rs de la façon suivante : R s = R TOT - R 2

où R2 = Lϕ/&tgr;R ; et
o) accomplir l'étape (m) entre les étapes (b) et (d).
Méthode selon la revendication 1, dans laquelle lesdites étapes (a) à (h) sont exécutées automatiquement à réception d'une commande envoyée par un outil d'entretien. Méthode selon la revendication 6, dans laquelle lesdites étapes (a) à (1) sont exécutées automatiquement à réception d'une commande envoyée par un outil d'entretien. Contrôleur d'ascenseur (7) comprenant un circuit logique (48) qui calcule automatiquement une constante de temps d'un moteur &tgr;R, une constante de couple KT*, et un courant de magnétisation Id pour un moteur/système d'entraînement à champ orienté en employant des moyens pour calculer une inductance transitoire L&sgr;, des moyens pour injecter une onde sinusoïdale ayant une fréquence d'entrée variable dans un courant de référence d'axe q IqREF et des moyens pour faire varier la fréquence d'entrée tout en calculant une partie imaginaire de l'impédance d'un rotor Imag(ZR) de façon à obtenir une fréquence FPEAK à laquelle la valeur maximum de Imag(ZR) est obtenue, et comprenant en outre des moyens pour calculer &tgr;R à partir de FPEAK, dans lequel les valeurs de KT*, Id, et une tension de moteur totale VM sont calculées, et des moyens pour faire varier Id jusqu'à ce qu'un rapport de la tension nominale du moteur Vph_RATED à VM se situe à l'intérieur d'une tolérance prédéterminée de 1, et des moyens pour accomplir ladite procédure avec le rotor bloqué.






IPC
A Täglicher Lebensbedarf
B Arbeitsverfahren; Transportieren
C Chemie; Hüttenwesen
D Textilien; Papier
E Bauwesen; Erdbohren; Bergbau
F Maschinenbau; Beleuchtung; Heizung; Waffen; Sprengen
G Physik
H Elektrotechnik

Anmelder
Datum

Patentrecherche

Patent Zeichnungen (PDF)

Copyright © 2008 Patent-De Alle Rechte vorbehalten. eMail: info@patent-de.com