WO2019093064A1 - 誘導電動機の駆動装置および駆動方法 - Google Patents
誘導電動機の駆動装置および駆動方法 Download PDFInfo
- Publication number
- WO2019093064A1 WO2019093064A1 PCT/JP2018/037936 JP2018037936W WO2019093064A1 WO 2019093064 A1 WO2019093064 A1 WO 2019093064A1 JP 2018037936 W JP2018037936 W JP 2018037936W WO 2019093064 A1 WO2019093064 A1 WO 2019093064A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- induction motor
- value
- resistance
- frequency
- constant
- Prior art date
Links
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
- H02P21/16—Estimation of constants, e.g. the rotor time constant
Definitions
- the present invention relates to a drive device and a drive method of an induction motor, and is particularly suitable as a drive device and a drive method of an induction motor mounted on a railway vehicle.
- vector control capable of realizing high response and high accuracy of torque control is widely spread.
- magnetic flux and current are defined in an orthogonal coordinate system that rotates in synchronization with the rotating magnetic field generated in the stator of the induction motor, and based on the equivalent circuit constant of the induction motor, AC voltage has a predetermined amplitude, phase and frequency
- the torque of the induction motor is freely controlled by applying at The resistance value, which is one of the equivalent circuit constants, fluctuates with the temperature change of the stator winding and the rotor conductor due to the change of the outside air temperature and the heat generation of the induction motor.
- the influence of the resistance constant error on the torque accuracy is determined by the proportion of the resistance component in the impedance of the induction motor.
- the impedance consists of a combination of reactance and resistance.
- the reactance is the product of the inductance and the angular frequency, and increases as the frequency increases.
- resistance has frequency dependency due to the effects of skin effect and proximity effect, its amount of change is smaller than reactance. That is, as the frequency approaches a low speed region, the proportion of the resistance component in the impedance increases. Therefore, when the induction motor is driven at a low speed range in a state where the resistance constant error is large, torque pulsation and deviation occur, and a problem occurs that the torque according to the command is not output.
- a technique for estimating the resistance constant has been proposed.
- Patent Document 1 a motor constant error at inverter restart is realized by interpolating the motor constant correction value during inverter stop and setting the complemented motor constant correction value as an initial value of the motor constant correction value at inverter restart.
- Patent Document 2 discloses a technique of newly providing a primary resistance estimating means for identifying a primary resistance, and using the primary resistance identified by the primary resistance estimating means for an initial value estimator that outputs an initial velocity and an initial secondary magnetic flux. It is done. This is a technology to identify the primary resistance based on the result of estimating the deviation of the d-axis secondary magnetic flux by the magnetic flux observer.
- Patent Document 1 in the resistance constant correction at the time of restart, since the correction is performed using the previous value (speed information at the time of stopping the inverter), the accuracy of the correction when the time from the previous stop time to the start is long. May decrease. In particular, after the induction motor is driven all day, at the first start of the next day, a temperature difference of about 100 ° C. with respect to the previous value may occur, which may cause an error in the correction value. Further, in the correction of the resistance constant, the change in the outside air temperature is not taken into consideration, and there is a problem in the estimation accuracy of the resistance constant.
- Patent Document 2 there is a possibility that sufficient accuracy can not be obtained in resistance resistance estimation of a low slippage high efficiency induction motor in recent years.
- the secondary resistance decreases and the secondary time constant (second flux time constant) increases, making it difficult for the magnetic flux to rise.
- the rising of the magnetic flux becomes particularly difficult under the condition where the primary frequency and the rotor frequency largely diverge (state in which the slip frequency is large).
- the method disclosed in Patent Document 2 does not take into consideration the influence of interference to the d axis due to the q axis secondary magnetic flux, and when the q axis secondary magnetic flux is generated, the error in the resistance estimated value is generated. May occur. For this reason, at the time of restart, there is a problem that sufficient estimation accuracy can not be obtained with the method based on the magnetic flux deviation amount.
- An object of the present invention is to drive an induction motor by applying an AC voltage at a primary frequency that intentionally increases the slip frequency, and based on the current detection value at that time and the setting value of the equivalent circuit constant By correcting the set value used in the control device to improve the torque accuracy when restarting from the low speed free run state.
- the drive device of the conductive motor detects a voltage output device that outputs an arbitrary AC voltage waveform according to the voltage command value and the primary frequency command value and applies the waveform to the induction motor, and detects the current flowing in the induction motor.
- the control device includes a current detector and a control device that controls the voltage output device, and the control device performs calculation based on the equivalent circuit of the induction motor using the current detection value detected by the current detector, and performs a resistance constant of the induction motor To estimate.
- the resistance constant can be estimated with high accuracy even at the time of restart even for a low-slip, high-efficiency induction motor, and the resistance constant used in the control device is corrected using the estimated value. It is possible to improve the torque accuracy when restarting from the run state.
- FIG. 5 is a diagram showing an outline of a flow of resistance constant estimation in the first embodiment.
- FIG. 7 is a diagram showing an example of a resistance constant estimation operation in the resistance estimation unit of the first embodiment.
- FIG. 8 is a view showing the effect of improving the torque accuracy according to the first embodiment. It is a figure which shows the functional block of the drive device of the induction motor which concerns on Example 2 of this invention. It is a figure which shows the characteristic of the relationship of the slip frequency and q axis current at the time of resistance change.
- FIG. 8 is a diagram showing an outline of a flow of resistance constant estimation in the second embodiment. It is a figure which shows the outline
- FIG. 17 is a diagram showing a schematic configuration of a portion of a railway vehicle equipped with the drive device for an induction motor according to any one of the first to fourth embodiments as a fifth embodiment of the present invention.
- FIG. 1 is a diagram showing a functional block of a drive device for an induction motor according to a first embodiment of the present invention.
- the configuration at the time of estimating the resistance constant and restarting by the speed sensorless vector control is shown as a functional block.
- FIG. 1 shows only the functional block of the minimum configuration necessary to calculate resistance constant estimation and initial velocity estimation, and power devices such as driving transistors such as IGBTs (Insulated Gate Bipolar Transistors) and diodes are illustrated.
- the power converter (main circuit) and other control configurations that are configured from the above are not shown.
- the drive device 1 drives the induction motor 3 by switching the control mode at the time of estimating the resistance constant, at the time of initial speed estimation and at the time of normal operation by the control mode switching device 20.
- the search signal command generator 31 for resistance resistance estimation is operated, and then the initial speed estimation command generator 8 is operated.
- the control mode is switched by the control mode switch 20, and the command value of the vector control unit 30 responsible for control during power running or regenerative operation is output.
- the integrator 21 integrates the phase angle ⁇ .
- the coordinate converter 5 performs three-phase conversion from dq coordinates, which are rotational coordinates, using the d-axis voltage command V d * , the q-axis voltage command V q * and the phase angle ⁇ to obtain a command value of a three-phase AC voltage waveform. Generate V u * , V v * and V w * .
- the voltage output device 2 converts the command values V u * , V v * and V w * of the three-phase AC voltage waveform into PWM (Pulse Width Modulation) signals, and then the drive circuit and main circuit (voltage output device (not shown) V.sub.u , V.sub.v and V.sub.w are applied to the induction motor 3 as a three-phase alternating current voltage.
- PWM Pulse Width Modulation
- the induction motor 3 generates rotational torque by the interaction of the rotating magnetic field generated on the stator side by the application of the three-phase AC voltage and the induced current flowing through the rotor conductor when the rotor is slipped.
- the current detector 4 is composed of a Hall CT (Current Transformer) or the like, and detects U-phase, V-phase and W-phase three-phase currents I u , I v and I w flowing through the induction motor 3 together with their waveforms. However, it is not necessary to detect currents of all three phases by the current detector 4. Even if one of the two phases is detected and the remaining one phase is calculated by calculation assuming that the three phase currents are in equilibrium, Good.
- Hall CT Current Transformer
- the coordinate converter 6 converts the three-phase currents I u , I v and I w of the induction motor detected by the current detector 4 into dq coordinates of the rotational coordinate system using the phase angle ⁇ , and calculates I d and I q. It outputs to the initial velocity estimation calculating unit 7 as
- the initial speed estimation calculation unit 7 estimates the rotational frequency of the rotor by an observer method or the like based on the current detection values of I d and I q , and outputs a rotor frequency estimated value f r ⁇ .
- the functional block has been shown above as the configuration of the driving device 1, but the configuration in which the voltage output device 2 and the current detector 4 are removed from the driving device 1 constitutes a control device, and the second to fifth embodiments described below. The same is true for
- Resistance estimation unit 25 calculates rotor frequency estimated value f r ⁇ , voltage command values V d * and V q * , primary frequency command value f 1 * , current detection value I d and equivalent circuit constant setting values r 1 * , R 2 * , l 1 * , l 2 * , and L M * are input, and a resistance variation ratio k r ⁇ indicating a ratio normalized by the resistance value of the reference temperature set in the control device is output.
- * represents a command value
- ⁇ represents an estimated value.
- FIG. 2 is a diagram showing a T-type equivalent circuit of the induction motor.
- Phase voltage vector V 1 Equation (1) the phase current vector I 1 is defined as in Equation (2), the synthesized impedance of the primary side and the secondary side, the combined impedance Z of the induction motor 3, the formula ( It is represented by 6).
- r 1 primary resistance
- r 2 secondary resistance
- l 1 primary leakage inductance
- l 2 secondary leakage inductance
- Z 1 primary impedance
- Z 2 S secondary impedance
- Y 0 excitation Admittance
- L M Excitation reactance
- slip Slip.
- slip is defined by equation (7) with f 1 as the primary frequency and f r as the rotor frequency.
- the drive range of the rotor frequency (electrical angle) is assumed to be -10 Hz to 200 Hz, but this is merely an example. However, even if the number of poles and the drive range change, the principle is the same, and the same effect can be obtained.
- the impedance of the induction motor is determined by the combination of reactance and resistance as shown in equations (3) to (6).
- the reactance is the product of the inductance and the angular frequency, and increases as the frequency increases.
- resistance has frequency dependency due to the effects of skin effect and proximity effect, its amount of change is smaller than reactance. That is, as the frequency approaches a low speed region, the proportion of the resistance component in the impedance increases. For this reason, if the resistance constant error occurs in a low speed range of about -10 Hz to 10 Hz, the torque accuracy will be degraded.
- FIG. 3 is a graph showing the characteristics of the relationship between the rotor frequency fr and the d-axis current at the time of resistance fluctuation when an AC voltage is applied at a primary frequency of 100 Hz. That explores signal command generator 31 from the primary frequency f 1 100 Hz, and outputs the constant amplitude of the alternating voltage, the rotor frequency at which the resistance variation ratio k r is changed from 0.7, 1.0 and 1.3 It is the calculation result of the d-axis current when f r is ⁇ 10 Hz to 200 Hz.
- the slip frequency f s which is a difference frequency between the primary frequency f 1 and the rotor frequency f r is described in the parentheses shown on the horizontal axis.
- the rotor frequency f r is for a range of ⁇ 10 Hz It can be seen that the amount of current deviation due to resistance change is small.
- FIG. 4 is a diagram showing the characteristics of the relationship between the rotor frequency and the d-axis current at the time of resistance fluctuation when an AC voltage is applied at a primary frequency of -20 Hz. That explores signal command generator 31 from the primary frequency f 1 -20Hz, outputs a constant amplitude of the AC voltage, the rotor at the time of the resistance variation ratio k r is changed from 0.7, 1.0 and 1.3 It is the calculation result of the d-axis current at the frequency f r of ⁇ 10 Hz to 200 Hz.
- the value of the voltage is determined so that
- the primary and secondary leakage inductances l 1 , l 2 and the excitation inductance L M of the induction motor coincide with the values set in the control device, and the voltage output device 2 receives the voltage V 1 as instructed.
- the primary frequency f 1 is output and the error in current detection is sufficiently small, at the time before the initial velocity estimation in the configuration shown in FIG.
- the child frequency f r and the resistance variation ratio k r are two.
- the resistance constant setting used in each control system is set using the resistance constant estimation result at the time of starting the rise of torque. It is possible to correct the value.
- FIG. 5 is a diagram schematically showing the flow of resistance constant estimation.
- step1 before initial rate estimate, an AC voltage of negative frequency opposite applied against the forward direction at a predetermined amplitude, detects a d-axis motor current I d at that time.
- Step 2 initial velocity estimation is performed to search for the rotor frequency f r .
- the rotor frequency estimate f r ⁇ is obtained, it can be one of the resistance variation ratio k r unknowns in equation (8) from equation (1).
- Step 3 the resistance constant used in the control device is corrected using the estimated resistance change ratio k r ⁇ , and the induction motor 3 is restarted by raising the torque using the corrected resistance constant.
- FIG. 6 is a diagram showing an example of the resistance constant estimation operation in the resistance estimation unit 25.
- the resistance fluctuation ratio k r ⁇ is obtained by comparing the d-axis current value for the rotor frequency virtually calculated based on the equivalent circuit constant in the resistance estimation unit 25 with the d-axis current detection value I d actually flowing through the motor. Calculate It is possible to calculate k r ⁇ from the d-axis current I d of Step 1 and the speed estimated value f r ⁇ of Step 2 for three curves having different k r calculated virtually.
- the resistance constant is corrected only when it is judged as the low speed region for simplification of processing. It may be configured to
- the lower limit of the driving range assuming-10 Hz, although the primary frequency f 1 of the probing signal is a -20Hz, if equal to or less than the lower limit value of the drive frequency, frequency is set to a number It does not matter. If a probe signal for constant estimation is applied within the drive range -10 Hz to 200 Hz and the rotor frequency and probe signal frequency become close (if the slip frequency is small), dq axis The magnetic flux rises up and affects the initial velocity estimation to be processed later. As apparent from the equivalent circuit shown in FIG. 2, the excitation current I 0 contributing to the magnetic flux increases as the slip slip decreases, and decreases as the slip slip increases.
- the primary frequency f 1 which explores signal command generator 31 for resistance constant estimation is output to the non-driving range-10 Hz ⁇ 200 Hz for the rotor frequency f r, and, significantly current deviation due to resistance change from both purposes to appear, it is desirable lower than the lower limit value of the driving range of the rotor frequency f r (in the first embodiment assumed-10 Hz).
- the set value of the resistance constant set in the control device can be corrected using estimated value k r ⁇ at the start of torque rise after initial speed estimation, and resistance constant error It is possible to stabilize the restart from a slow free run condition that is susceptible to
- FIG. 7 is a view showing the effect of improving the torque accuracy according to the first embodiment, and shows an example of the calculation result at the time of restarting with or without correction by resistance resistance estimation.
- the rotor frequency f r is in a free run state of 0.5 Hz.
- k r is a resistance change ratio of the induction motor
- k r * is a resistance change ratio set by the control device.
- Example 1 is a method of estimating a resistance constant without using a magnetic flux, and a second-order time constant T 2 (second-order self-inductance L 2 / first-order converted second-order resistance r It is also possible to cope with low slip induction motors where 2 ') is large and magnetic flux does not easily rise.
- the resistance constant estimation method using the d-axis current detection value has been described, but this may be a q-axis current detection value.
- the characteristics of the current depend on the AC voltage to be applied and the equivalent circuit constant of the induction motor, and are determined based on the equations (1) to (8). Therefore, an axis with high sensitivity to current change at resistance fluctuation is appropriately selected Shall be used.
- the same speed can be obtained by using the detected speed information f r instead of the speed estimated value f r ⁇ of Step 2 shown in FIG.
- the resistance constant can be estimated by the principle.
- the first embodiment is a resistance estimation method mainly focusing on the change of the term of r 2 / slip shown in FIG. 2, and it is assumed that the primary resistance also fluctuates at a similar ratio to the secondary resistance. . However, since the temperature rise and the temperature coefficient are different between the primary resistance of the stator winding and the secondary resistance of the rotor conductor, the resistance variation ratio between the two is strictly different.
- the compensation ratio of the primary resistance may be changed based on the fluctuation ratio of the secondary resistance estimated based on the temperature rise ratio and the temperature coefficient ratio of the two.
- FIG. 8 is a diagram showing a functional block of a drive device for an induction motor according to a second embodiment of the present invention. In the stop or free run state of the induction motor, while estimating the resistance constant, the configuration at the time of restarting by the speed sensorless vector control is shown as a functional block.
- the second embodiment aims to shorten the time (resistance constant estimation and initial speed estimation) required to start the torque rise of the restart more than the first embodiment. Therefore, in the second embodiment, the q-axis current detection value I q is input to the resistance constant estimator 25 without using the search signal generator 31 of the first embodiment, and the resistance variation ratio k r ⁇ is output. It differs from the first embodiment.
- the resistance fluctuation ratio k r ⁇ is estimated.
- the rotor frequency f r needs to be an existing value.
- the initial value of the rotor frequency used for the initial velocity estimation is assumed to be approximately the maximum rotation frequency (assumed to be around 200 Hz in the second embodiment), and the target velocity range for resistance constant estimation is limited to the low speed region (only the scope of the present embodiment 2 in-10 Hz ⁇ 10 Hz assumed), at the time of initial velocity estimation starts, the slip frequency becomes sufficiently large state, to reduce the sensitivity of the rotor frequency f r is unknown, the initial It is possible to estimate the resistance constant in parallel processing (initial speed estimation processing and resistance constant estimation processing) during the speed estimation period. The principle is described below.
- the q-axis current I q of the induction motor is defined by the following equation.
- L ⁇ is a linear converted inductance
- L 2 is a secondary self inductance.
- the primary and secondary leakage inductances l 1 and l 2 and the excitation inductance L M are constant, and the voltage output device 2 outputs the voltage V 1 and the primary frequency f 1 as instructed. and the error of the current detection is sufficiently small, and assuming the formula (9) to the rotor frequency f r and the secondary resistance r 2 is not included, that variations of the primary resistance r 1 appears as a q-axis current error Recognize.
- FIG. 9 is a diagram showing the characteristic of the relationship between the slip frequency and the q-axis current at the time of resistance change.
- k r resistance variation ratio
- k r resistance variation ratio
- 1.3 show the results rotor frequency f r is calculated characteristics of the q-axis current I q at -20Hz ⁇ 190Hz.
- the slip frequency f s is sufficiently large (around 190 Hz)
- the slip frequency f s fluctuates by about ⁇ 10 Hz
- the influence on the q-axis current I q is small.
- the primary frequency f 1 is set to a high value so that the q-axis current I q becomes asymptotically slip frequency f s If the voltage is set to 0, it means that the influence of the resistance fluctuation appears as a current deviation without being influenced by the magnitude of the rotor frequency f r .
- FIG. 10 is a diagram showing an outline of a flow of resistance constant estimation in the second embodiment.
- the resistance fluctuation is determined from the q-axis current detection value I q , the d-axis voltage command value V d * , the q-axis voltage command value V q * , the primary frequency command value f 1 * and the equivalent circuit constant Calculate the ratio k r ⁇ .
- the resistance constant estimation unit 25 outputs k r ⁇ by the calculation, but the low speed determiner 32 determines the low speed range of the rotor frequency estimated value f r ⁇ by the set threshold F x ( For example, the resistance constant is updated only when the frequency becomes 10 Hz or less.
- the estimation calculation of the resistance constant is performed only under the condition that the primary frequency f 1 is large (for example, 150 to 200 Hz).
- the primary frequency f 1 for example, 150 to 200 Hz.
- the estimated value of the resistance constant is not affected. It is because it assumes that.
- the second embodiment is a resistance estimating method focusing on the change of the term of r 1 shown in FIG. 2 as described in the equation (9), and the ratio of the secondary resistance r 2 to the primary resistance r 1 is also the same. It is assumed that it fluctuates in However, the primary resistance r 1 of the stator winding, the secondary resistance r 2 of the rotor conductor, also differ each temperature increase and temperature coefficient, both the resistance variation ratio of strictly different. Based on the temperature rise ratio and temperature coefficient ratio of the two, based on the estimated variation ratio of the primary resistance r 1, it may be configured to change the compensation ratio of the secondary resistance r 2.
- the resistance change ratio can be estimated by parallel processing during the initial speed estimation period, and restart can be performed in a shorter time than the first embodiment.
- FIG. 11 is a diagram showing an outline of a flow of resistance constant estimation in the third embodiment of the present invention.
- the third embodiment aims to stabilize the estimation of the resistance change ratio more than the second embodiment. Therefore, although the configuration of the drive device for an induction motor according to the third embodiment is the same as that of the second embodiment, in order to achieve the above object, an integrator for the q-axis current I q serving as the input of the resistance estimator 25 is provided. This embodiment differs from the second embodiment in that 26 is provided and integration processing is performed.
- the integration process can reduce the above-mentioned influence, and resistance estimation can be made more stably than in the second embodiment.
- the resistance estimation unit 25 integration processing is similarly performed on the current value virtually calculated using Expressions (1) to (8).
- the range for resistance estimation is set to the range of the high speed range (for example, 150 Hz to 200 Hz) as in the second embodiment.
- the third embodiment can estimate the resistance constant more stably than the second embodiment.
- FIG. 12 is a diagram showing an outline of a flow of resistance constant estimation in the fourth embodiment of the present invention.
- the fourth embodiment aims to improve the accuracy of resistance constant estimation more than the third embodiment. Therefore, the drive motor configuration of the induction motor according to the fourth embodiment is the same as that of the second embodiment, but in order to achieve the above object, the induction motor of the third embodiment provided in the third embodiment.
- the second embodiment differs from the second embodiment and the third embodiment in that a low pass filter 27 set to the time constant of the motor electrical time constant is provided with respect to the current detection value (q axis current I q in FIG. 12).
- the fourth embodiment can improve the accuracy of the resistance constant estimation and can further improve the torque accuracy more than the third embodiment.
- FIG. 13 is a view showing a schematic configuration of a portion of a railway vehicle equipped with the drive device for an induction motor according to any one of the first to fourth embodiments as a fifth embodiment according to the present invention.
- drive device 1 mounted on a railway vehicle receives supply of power from overhead wire 101 via a current collector, outputs AC power and supplies it to induction motor 3 so that electric energy becomes mechanical torque. It is converted.
- the induction motor 3 is connected to the axle of the railway vehicle via a reduction gear, and the railway vehicle travels by the tangential force generated between the wheel 103 and the rail 102 connected to the axle.
- the drive device 1 for an induction motor according to any one of the first to fourth embodiments described above is applied as the drive device 1 mounted on a railway vehicle.
- the torque shock of the induction motor 3 can be reduced at the time of restart from the low speed free run state.
- the body vibration of the railway vehicle can be reduced, and the passenger's ride quality can be further improved.
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
インバータ等の電力変換器により駆動する誘導電動機を低速フリーラン状態から再起動する際に、誘導電動機の温度変動によって生じる抵抗定数の誤差によるトルクの脈動や偏差を低減しトルク精度を向上させるために、誘導電動機の駆動装置は、電圧指令値および一次周波数指令値に応じて任意の交流電圧の波形を出力して誘導電動機に印加する電圧出力装置と、誘導電動機に流れる電流を検出する電流検出器と、電圧出力装置を制御する制御装置とを備え、誘導電動機の抵抗定数を補正するために、制御装置は、電流検出器が検出した電流検出値を用いて誘導電動機の等価回路に基づく演算を行い当該誘導電動機の抵抗定数を推定する。
Description
本発明は、誘導電動機の駆動装置および駆動方法に関し、特に、鉄道車両に搭載する誘導電動機の駆動装置および駆動方法として好適である。
鉄道車両用の誘導電動機の駆動装置では、トルク制御の高応答化、高精度化が実現可能なベクトル制御が広く普及している。誘導電動機のベクトル制御では、誘導電動機の固定子に生じる回転磁界と同期回転する直行座標系で磁束と電流を定義し、誘導電動機の等価回路定数に基づき、交流電圧を所定の振幅、位相および周波数で印加することで、誘導電動機のトルクを自在に制御している。
この等価回路定数の1つである抵抗値は、外気温度の変化や誘導電動機の発熱による固定子巻線や回転子導体の温度変化に伴って変動する。
この等価回路定数の1つである抵抗値は、外気温度の変化や誘導電動機の発熱による固定子巻線や回転子導体の温度変化に伴って変動する。
ベクトル制御を実装した制御器にセットした抵抗定数の設定値と誘導電動機の実際の抵抗定数とに差異が生じると、抵抗定数誤差となり、制御の不安定化やトルク精度が低下する問題が発生する。
抵抗定数誤差がトルク精度に与える影響は、誘導電動機のインピーダンスにおける抵抗成分の占める割合によって決まる。インピーダンスは、リアクタンスと抵抗の合成から成る。リアクタンスは、インダクタンスと角周波数の積であるため、周波数が高くなる程大きくなる。一方、抵抗は、表皮効果や近接効果の影響で周波数の依存性はあるものの、リアクタンスに比べるとその変化量は小さい。つまり、周波数が低い低速域に近付く程、インピーダンスにおける抵抗成分の割合が大きくなる。
そのため、抵抗定数誤差が大きい状態で、低速域で誘導電動機を駆動すると、トルクの脈動や偏差が生じ、指令通りのトルクが出力されない問題が発生する。この抵抗定数誤差に起因したトルク精度の低下を解決すべく、抵抗定数の推定技術が提案されている。
そのため、抵抗定数誤差が大きい状態で、低速域で誘導電動機を駆動すると、トルクの脈動や偏差が生じ、指令通りのトルクが出力されない問題が発生する。この抵抗定数誤差に起因したトルク精度の低下を解決すべく、抵抗定数の推定技術が提案されている。
特許文献1では、インバータ停止中のモータ定数補正値を補間し、補完したモータ定数補正値をインバータ再起動時のモータ定数補正値の初期値として設定することにより、インバータ再起動時のモータ定数誤差を抑制する技術が開示されている。
特許文献2では、一次抵抗を同定する一次抵抗推定手段を新たに設け、初期速度と初期二次磁束を出力する初期値推定器に対して一次抵抗推定手段で同定した一次抵抗を用いる技術が開示されている。磁束オブザーバによりd軸二次磁束の偏差量を推定した結果を基に一次抵抗を同定する技術である。
本願発明者が、低速域における抵抗定数誤差の影響を小さくし、トルク精度の改善ついて鋭意検討した結果、次の知見を得るに至った。
特許文献1では、再起動時の抵抗定数補正においては、前回値(インバータ停止時の速度情報)を用いて補正するため、前回の停止時間から起動までの時間が長い場合には、補正の精度が低下する可能性がある。特に、誘導電動機を終日駆動させた後、翌日の初回の起動では、前回値に対して約100℃近くの温度差が生じることもあり、補正値に誤差が生じる可能性がある。また、抵抗定数の補正には、外気温度の変化も考慮されておらず、抵抗定数の推定精度に課題がある。
特許文献1では、再起動時の抵抗定数補正においては、前回値(インバータ停止時の速度情報)を用いて補正するため、前回の停止時間から起動までの時間が長い場合には、補正の精度が低下する可能性がある。特に、誘導電動機を終日駆動させた後、翌日の初回の起動では、前回値に対して約100℃近くの温度差が生じることもあり、補正値に誤差が生じる可能性がある。また、抵抗定数の補正には、外気温度の変化も考慮されておらず、抵抗定数の推定精度に課題がある。
特許文献2では、近年の低すべり化された高効率誘導電動機の抵抗定数推定では、十分な精度が得られない可能性がある。低すべりの誘導電動機では、二次抵抗が小さくなることで二次時定数(二次磁束の時定数)が大きくなり、磁束が立ち上がり難くなっている。特に、速度情報が不定である初期速度推定期間で、一次周波数と回転子周波数が大きく乖離した条件(すべり周波数が大きい状態)では、磁束の立ち上げは特に難しくなる。また、特許文献2で開示されている方式は、q軸二次磁束によるd軸への干渉の影響が考慮されておらず、q軸二次磁束が発生した場合には、抵抗推定値に誤差が発生することが考えられる。このような理由から、再起動時においては、磁束偏差量に基づいた方式では十分な推定精度が得られない課題がある。
本発明の目的は、誘導電動機の駆動装置において、意図的にすべり周波数が大きくなるような一次周波数で交流電圧を印加し、その際の電流検出値と等価回路定数の設定値に基づいて抵抗定数を推定し、制御装置で用いる設定値を補正することで低速フリーラン状態から再起動する時のトルク精度を向上させることである。
本発明に係る導電動機の駆動装置は、電圧指令値および一次周波数指令値に応じて任意の交流電圧の波形を出力して誘導電動機に印加する電圧出力装置と、誘導電動機に流れる電流を検出する電流検出器と、電圧出力装置を制御する制御装置とを備え、該制御装置は、電流検出器が検出した電流検出値を用いて誘導電動機の等価回路に基づく演算を行い当該誘導電動機の抵抗定数を推定することを特徴とする。
本発明によれば、低すべり化された高効率誘導電動機に対しても、再起動時に抵抗定数を高精度に推定でき、推定値を用いて制御装置で用いる抵抗定数を補正することで低速フリーラン状態からの再起動時にトルク精度を向上させることが可能になる。
以下、本発明を実施するための形態として、実施例1~5について、それぞれ図面に従い詳細に説明する。各実施例において、参照番号が同一のものは、同一の構成要件または類似の機能を備えた構成要件を示している。なお、以下に説明する構成は、あくまでも実施例に過ぎず、本発明に係る実施態様は、以下の実施例に限定されるものではない。
図1は、本発明の実施例1に係る誘導電動機の駆動装置の機能ブロックを示す図である。誘導電動機の停止またはフリーラン状態において、抵抗定数を推定し、速度センサレスのベクトル制御によって再起動する際の構成を、機能ブロックとして示したものである。
図1に示す駆動装置1には、負荷として接続する誘導電動機3を速度センサレスでベクトル制御するための制御プログラムが実装されている。ただし、図1は、抵抗定数推定と初期速度推定を演算するために必要な最小構成の機能ブロックのみを示したもので、IGBT(Insulated Gate Bipora Transistor)等の駆動用トランジスタやダイオード等のパワーデバイスから構成される電力変換器(主回路)および他の制御構成については、図示を省略している。
駆動装置1は、制御モード切替器20によって、抵抗定数推定時、初期速度推定時および通常運転時とで制御モードを切り替えて誘導電動機3を駆動する。誘導電動機3の回転子が停止またはフリーラン状態から通電を開始して再起動する場合、抵抗定数推定用の探り信号指令発生器31を動作させた後に、初期速度推定用指令発生器8を動作させて、ベクトル制御における回転座標上の励磁軸(磁束軸)であるd軸の電圧指令Vd
*、d軸と直交するq軸の電圧指令Vq
*および一次周波数指令f1
*を出力する。抵抗定数推定と初期速度推定の完了後、制御モード切換器20によって制御モードを切り替え、力行または回生運転時の制御を担うベクトル制御部30の指令値を出力させる。
制御モード切換器20を介して出力される3つの指令(d軸電圧指令Vd
*、q軸電圧指令Vq
*および一次周波数指令f1
*)の内、一次周波数指令f1
*に関しては、積分器21により位相角θに積分演算する。ただし、位相角θの演算は、実際には一次角周波数指令ω1
*(=2π・f1
*)であるが、簡略的にf1
*で表記している(以下も同様)。
座標変換器5は、d軸電圧指令Vd
*、q軸電圧指令Vq
*および位相角θを用いて、回転座標であるdq座標から3相変換を行い、3相交流電圧波形の指令値Vu
*、Vv
*およびVw
*を生成する。
電圧出力装置2は、3相交流電圧波形の指令値Vu
*、Vv
*およびVw
*をPWM(Pulse Width Modulation)信号に変換した後、駆動回路および主回路(図示しないが電圧出力装置2に含まれる)を介して、Vu、VvおよびVwを3相交流電圧として誘導電動機3に印加する。
誘導電動機3は、3相交流電圧の印加により固定子側に発生した回転磁界と、回転子にすべりが生じたときに回転子導体に流れる誘導電流との相互作用により、回転トルクを発生する。
電流検出器4は、ホールCT(Current Transformer)等から構成され、誘導電動機3に流れるU相、V相およびW相の3相電流Iu、IvおよびIwをその波形と共に検出する。ただし、電流検出器4によって必ずしも3相全ての電流を検出する必要はなく、いずれかの2相を検出し、残る1相は3相電流が平衡状態であると仮定して演算により求める構成でもよい。
座標変換器6は、電流検出器4で検出した誘導電動機の3相電流Iu、IvおよびIwを、位相角θを用いて回転座標系のdq座標に変換し、IdおよびIqとして初期速度推定演算部7に出力する。
初期速度推定演算部7は、IdおよびIqの電流検出値に基づき、オブザーバ方式等で回転子の回転周波数を推定し、回転子周波数推定値fr^を出力する。
以上、駆動装置1の構成として機能ブロックを示したが、該駆動装置1から電圧出力装置2および電流検出器4を除いた構成が制御装置を構成するものであり、以下の実施例2~5においても同様である。
次に、本実施例1の要点となる、抵抗推定部25で抵抗定数を推定する方法について説明する。
抵抗推定部25は、回転子周波数推定値fr^、電圧指令値Vd *とVq *、一次周波数指令値f1 *、電流検出値Idおよび等価回路定数の各設定値r1 *、r2 *、l1 *、l2 *、LM *を入力とし、制御装置において設定した基準温度の抵抗値で正規化した比率を示す抵抗変動比kr^を出力する。
ここにおいて、抵抗変動比kr^の推定原理について説明する。ただし、本実施例1において、*は指令値、^は推定値を表すものとする。
抵抗推定部25は、回転子周波数推定値fr^、電圧指令値Vd *とVq *、一次周波数指令値f1 *、電流検出値Idおよび等価回路定数の各設定値r1 *、r2 *、l1 *、l2 *、LM *を入力とし、制御装置において設定した基準温度の抵抗値で正規化した比率を示す抵抗変動比kr^を出力する。
ここにおいて、抵抗変動比kr^の推定原理について説明する。ただし、本実施例1において、*は指令値、^は推定値を表すものとする。
図2は、誘導電動機のT型等価回路を示す図である。相電圧ベクトルV1を式(1)、相電流ベクトルI1を式(2)のように定義し、一次側および二次側のインピーダンスを合成すると、誘導電動機3の合成インピーダンスZは、式(6)で表わされる。
ここで、r1:一次抵抗、r2:二次抵抗、l1:一次漏れインダクタンス、l2:二次漏れインダクタンス、Z1:一次側インピーダンス、Z2S:二次側インピーダンス、Y0:励磁アドミタンス、LM:励磁リアクタンス、slip:すべり、とする。
本実施例1では、4極の誘導電動機を搭載した鉄道車両において、回転子周波数(電気角)の駆動範囲を-10Hz~200Hzと仮定して説明するが、これはあくまでも一例に過ぎない。しかし、極数と駆動範囲が変わった場合でも原理は同一であるため、同様の効果が得られる。
誘導電動機のインピーダンスは、式(3)から式(6)に示されるように、リアクタンスと抵抗の合成によって決まる。リアクタンスは、インダクタンスと角周波数の積であるため、周波数が高くなる程大きくなる。一方、抵抗は、表皮効果や近接効果の影響で周波数の依存性はあるものの、リアクタンスに比べるとその変化量は小さい。つまり、周波数が低い低速域に近付く程、インピーダンスにおける抵抗成分の割合が大きくなる。このような理由から、-10Hz~10Hz程度の低速域で抵抗定数誤差が生じると、トルク精度が低下してしまうことになる。
図3は、一次周波数100Hzで交流電圧を印加した際の抵抗変動時における回転子周波数frとd軸電流との関係の特性を示す図である。すなわち、探り信号指令発生器31から一次周波数f1100Hz、振幅一定の交流電圧を出力し、抵抗変動比krを0.7、1.0および1.3と変化させた際の回転子周波数frが-10Hz~200Hz時におけるd軸電流の算出結果である。なお、横軸に示す丸括弧内には、一次周波数f1と回転子周波数frとの差分周波数となるすべり周波数fsを記載している。
一次周波数f1と回転子周波数frが近い領域(すべり周波数fsが±30Hz程度以内)では、抵抗変動による電流値の変化が大きいものの、回転子周波数frが±10Hzの範囲に対しては、抵抗変動による電流偏差量が小さいことがわかる。
発明者らは、回転子の正転方向に対して逆向き(負の周波数)で、回転子周波数frの下限値より低い周波数(-10Hz以下、絶対値でいえば下限値よりも大きい周波数となる)を印加することで、低速域における電流検出値の変化量を大きくできることを明らかにした。
図4は、一次周波数-20Hzで交流電圧を印加した際の抵抗変動時における回転子周波数とd軸電流との関係の特性を示す図である。すなわち、探り信号指令発生器31から一次周波数f1-20Hz、振幅一定の交流電圧を出力し、抵抗変動比krを0.7、1.0および1.3と変化させた際の回転子周波数frが-10Hz~200Hz時におけるd軸電流の算出結果である。図4では、図3の条件に対し、|V/F|(電圧と一次周波数の絶対値の比)が同一となるように電圧の値を決めており、縦軸のスケールも同一としている。
回転子周波数frが-10Hz~10Hz程度の領域の温度変動に対する検出電流の変化量が図3よりも顕著に表れていることから、本発明に係る手法を用いることで、抵抗定数変動の影響を電流検出値から推定することが可能となる。
ここで、誘導電動機の一次、二次の漏れインダクタンスl1、l2および励磁インダクタンスLMが制御装置に設定している値と一致しており、電圧出力装置2が指令通りの電圧V1と一次周波数f1を出力し、電流検出の誤差が十分に小さい、と仮定すると、図1に示す構成での初期速度推定前の時点では、式(1)から式(8)における未知数は、回転子周波数frと抵抗変動比krの2つとなる。
すなわち、初期速度推定が完了した時点で回転子周波数推定値fr^が得られれば、トルクの立上げを開始する時点で、抵抗定数推定結果を用いて、各制御系に用いられる抵抗定数設定値を補正することが可能となる。
図5は、抵抗定数推定のフローの概略を示す図である。
Step1では、初期速度推定前に、所定の振幅で正転方向に対して逆向きの負の周波数の交流電圧を印加し、その際のd軸モータ電流Idを検出する。
Step1では、初期速度推定前に、所定の振幅で正転方向に対して逆向きの負の周波数の交流電圧を印加し、その際のd軸モータ電流Idを検出する。
Step2では、初期速度推定を実施し回転子周波数frを探索する。初期速度推定を完了した時点で、回転子周波数推定値fr^が得られることから、式(1)から式(8)における未知数を抵抗変動比krの1つにすることができる。
Step3では、推定した抵抗変動比kr^を用いて制御装置で用いる抵抗定数を補正し、補正後の抵抗定数を用いて、トルクを立ち上げて誘導電動機3を再起動する。
図6は、抵抗推定部25における抵抗定数推定演算の一例を示す図である。
探り信号指令発生器31で発生させた交流電圧と、抵抗変動比krを変化させた抵抗変動時のd軸電流を、式(1)から式(8)に基づき、抵抗推定部25で仮想的に算出する(図に示すkrが異なる3本の曲線)。抵抗推定部25で、等価回路定数に基づいて仮想的に算出したロータ周波数に対するd軸電流値と、実際にモータに流れたd軸電流検出値Idを比較することで抵抗変動比kr^を算定する。仮想的に算出したkrが異なる3本の曲線に対し、Step1のd軸電流IdおよびStep2の速度推定値fr^により、kr^が算出可能である。
探り信号指令発生器31で発生させた交流電圧と、抵抗変動比krを変化させた抵抗変動時のd軸電流を、式(1)から式(8)に基づき、抵抗推定部25で仮想的に算出する(図に示すkrが異なる3本の曲線)。抵抗推定部25で、等価回路定数に基づいて仮想的に算出したロータ周波数に対するd軸電流値と、実際にモータに流れたd軸電流検出値Idを比較することで抵抗変動比kr^を算定する。仮想的に算出したkrが異なる3本の曲線に対し、Step1のd軸電流IdおよびStep2の速度推定値fr^により、kr^が算出可能である。
なお、本実施例1では、kr=0.7、1.0および1.3に対して、仮想的なd軸電流を算出しているが、式(1)から式(8)の数式を用いて制御ソフト内で演算させる方法や、計算結果や実測結果を事前にマップ化して実装する方法でもよい。
また、前述したように、中高速域では抵抗定数誤差の影響が生じても、トルク精度に与える影響が小さいため、処理の簡素化のために、低速域と判断した場合のみ、抵抗定数を補正する構成にしてもよい。
更に、本実施例1では、駆動範囲の下限値を-10Hzと仮定し、探り信号の一次周波数f1を-20Hzとしたが、駆動周波数の下限値以下となれば、周波数はいくつに設定しても構わない。仮に、駆動範囲内-10Hzから200Hzの範囲で、定数推定用の探り信号を印加し、回転子周波数と探り信号周波数が近くなってしまった場合(すべり周波数が小さい場合)には、dq軸の磁束が立ち上がり、その後の処理となる初期速度推定に影響を与えてしまう。図2に示す等価回路からも明らかであるように、磁束に寄与する励磁電流I0は、すべりslipが小さくなると増加し、すべりslipが大きいと減少する。
このような理由から、抵抗定数推定用の探り信号指令発生器31が出力する一次周波数f1は、回転子周波数frの駆動範囲-10Hz~200Hz以外とし、かつ、抵抗変動による電流偏差が顕著に現れるようにする両目的から、回転子周波数frの駆動範囲の下限値(本実施例1では-10Hzと仮定)よりも低いことが望ましい。
以上の構成とすることで、初期速度推定後のトルク立上げ開始時点で、制御装置に設定している抵抗定数の設定値を推定値kr^を用いて補正することができ、抵抗定数誤差の影響を受け易い低速フリーラン状態からの再起動を安定化できる。
図7は、本実施例1によるトルク精度向上の効果を示す図であり、抵抗定数推定による補正の有無における再起動時の計算結果の一例を示す。回転子周波数frは0.5Hzのフリーラン状態としている。また、krは誘導電動機の抵抗変動比で、kr
*は制御装置で設定している抵抗変動比である。
図7(a)に示す抵抗定数推定なしの場合では、トルク脈動が発生し、さらに指令値に対し定常偏差が発生し、トルクが小さくなっている。
一方、図7(b)に示す抵抗定数推定ありの場合では、Step3のトルク立上げに移行した時点で、kr *がkrとほぼ一致しており、抵抗定数推定した結果を用いて制御装置で用いる設定値を更新することで、実際のモータの抵抗値とほぼ同等に設定できる。その結果、トルクの脈動と定常偏差を低減でき、トルク精度を向上させることができる。
一方、図7(b)に示す抵抗定数推定ありの場合では、Step3のトルク立上げに移行した時点で、kr *がkrとほぼ一致しており、抵抗定数推定した結果を用いて制御装置で用いる設定値を更新することで、実際のモータの抵抗値とほぼ同等に設定できる。その結果、トルクの脈動と定常偏差を低減でき、トルク精度を向上させることができる。
また、Step1および2の期間において、探り信号を印加しているため、d軸二次磁束Φ2dが僅かに振動しているものの非常に小さい振動であり、直流量は発生していないことがわかる。本結果からも明らかであるように、本実施例1は、磁束を用いずに抵抗定数を推定する方式であり、二次時定数T2(二次自己インダクタンスL2/一次換算二次抵抗r2’)が大きく、磁束が立ち上がり難い低すべり誘導電動機にも対応可能である。
なお、本実施例1では、d軸電流検出値を用いた抵抗定数推定手法を説明したが、これはq軸電流検出値としても構わない。電流の特性は、印加する交流電圧や誘導電動機の等価回路定数に依存し、式(1)から式(8)に基づいて決まることから、抵抗変動時における電流変化の感度が高い軸を適宜選択し用いるものとする。
また、速度センサレス制御の構成を用いて説明したが、速度センサ付きの場合は、図5に示すStep2の速度推定値fr^の代わりに、検出した速度情報frを用いることにより、同一の原理で抵抗定数を推定可能である。
本実施例1は、主に、図2に示すr2/slipの項の変化に着目した抵抗推定手法であり、二次抵抗に対し一次抵抗も同様の比率で変動するという仮定を置いている。しかし、固定子巻線の一次抵抗と回転子導体の二次抵抗とでは、それぞれ温度上昇や温度係数も異なるため、両者の抵抗変動比は厳密には異なる。両者の温度上昇比や温度係数比を踏まえ、推定した二次抵抗の変動比を基準に、一次抵抗の補償比を変える構成にしてもよい。
以上の構成により、初期速度推定完了時点で、誘導電動機の抵抗値を推定可能であるため、トルク立上げ時のトルクショックを低減し、トルク精度を向上させることができる。
図8は、本発明の実施例2に係る誘導電動機の駆動装置の機能ブロックを示す図である。誘導電動機の停止またはフリーラン状態において、抵抗定数を推定しつつ、速度センサレスのベクトル制御によって再起動する際の構成を、機能ブロックとして示したものである。
本実施例2は、実施例1よりも再起動のトルク立上げ開始までにかかる時間(抵抗定数推定と初期速度推定)を短縮することを目的とする。そのため、本実施例2では、実施例1の探り信号発生器31を用いずに、q軸電流検出値Iqを抵抗定数推定器25に入力し、抵抗変動比kr^を出力する点で実施例1と異なる。
抵抗変動比kr^の推定には、実施例1に記載したように式(1)から式(8)における未知数が2つ(抵抗変動比kr、回転子周波数fr)あるため、本来は回転子周波数frが既値である必要がある。
本実施例2では、初期速度推定に用いる回転子周波数の初期値を略最高回転周波数とし(本実施例2では200Hz付近と仮定)、抵抗定数推定する対象速度範囲を低速域に限定することで(本実施例2では-10Hz~10Hzの範囲のみと仮定)、初期速度推定開始時においては、すべり周波数が十分に大きい状態となり、未知数である回転子周波数frの感度を小さくして、初期速度推定期間中に並列処理(初期速度推定処理と抵抗定数推定処理)で抵抗定数を推定することが可能である。以下に、その原理について説明する。
すべり周波数を無限大に近付けると、誘導電動機のq軸電流Iqは次式にて定義される。
ただし、Lσは一次換算したインダクタンス、L2は二次自己インダクタンスとする。ここで、実施例1と同様に、一次、二次の漏れインダクタンスl1、l2および励磁インダクタンスLMが一定であり、電圧出力装置2が指令通りの電圧V1と一次周波数f1を出力し、電流検出の誤差が十分に小さい、と仮定すると、式(9)に回転子周波数frや二次抵抗r2は含まれず、一次抵抗r1の変動がq軸電流誤差として現れることがわかる。
図9は、抵抗変動時におけるすべり周波数とq軸電流との関係の特性を示す図である。抵抗変動比krを0.7、1.0および1.3と変化させた際の、回転子周波数frが-20Hz~190Hz時におけるq軸電流Iqの特性を算出した結果を示す。すべり周波数fsが十分に大きい領域(190Hz付近)では、すべり周波数fsが±10Hz程度変動しても、q軸電流Iqに与える影響は小さいことがわかる。言い換えれば、実際の回転子周波数frが低速(-10Hz~10Hz)で回転していたとして、q軸電流Iqが漸近する程のすべり周波数fsになるように一次周波数f1を高い値に設定し電圧を印加すれば、回転子周波数frの大きさに影響されることなく、抵抗変動の影響が電流偏差として現れることを意味する。
誘導電動機3の二次抵抗r2が小さくなる程、q軸電流検出値は漸近し易い傾向にあり、1時間定格出力におけるすべりが1.0%以下となるような低すべり化された高効率誘導電動機に対して、適した方式である。
図10は、本実施例2における抵抗定数推定のフローの概要を示す図である。Step1で示す初期速度推定期間に、q軸電流検出値Iq、d軸電圧指令値Vd
*、q軸電圧指令値Vq
*、一次周波数指令値f1
*および等価回路定数から、抵抗変動比kr^を算出する。
初期速度推定演算を実施する都度、抵抗定数推定部25による演算によりkr^を出力するが、低速判定器32によって、回転子周波数推定値fr^が低速域を判断する設定閾値Fx(例えば、10Hz)以下になった場合のみ、抵抗定数を更新するものとする。
また、抵抗定数の推定演算は一次周波数f1が大きい条件のみ(例えば150~200Hz)実施するものとする。これは、本実施例2では、図9に示したように、すべり周波数fsが十分に大きく、回転子周波数frが±10Hz程度変動しても、抵抗定数の推定値へ影響を与えないことを前提としているためである。
本実施例2は、式(9)で説明したように、図2に示すr1の項の変化に着目した抵抗推定手法であり、一次抵抗r1に対し二次抵抗r2も同様の比率で変動するという仮定を置いている。しかし、固定子巻線の一次抵抗r1と、回転子導体の二次抵抗r2では、それぞれ温度上昇や温度係数も異なるため、両者の抵抗変動比は厳密には異なる。両者の温度上昇比や温度係数比を踏まえ、推定した一次抵抗r1の変動比を基準に、二次抵抗r2の補償比を変える構成にしてもよい。
以上の構成により、本実施例2は、初期速度推定期間中に並列処理で抵抗変動比を推定でき、実施例1よりも短時間での再起動が可能となる。
図11は、本発明の実施例3における抵抗定数推定のフローの概要を示す図である。
本実施例3は、実施例2よりも抵抗変動比の推定を安定化することを目的とする。したがって、本実施例3に係る誘導電動機の駆動装置の構成は、実施例2と同様であるが、上記目的を達成するために、抵抗推定器25の入力となるq軸電流Iqに積分器26を設け積分処理する点で実施例2と異なる。
本実施例3は、実施例2よりも抵抗変動比の推定を安定化することを目的とする。したがって、本実施例3に係る誘導電動機の駆動装置の構成は、実施例2と同様であるが、上記目的を達成するために、抵抗推定器25の入力となるq軸電流Iqに積分器26を設け積分処理する点で実施例2と異なる。
電流検出値の瞬時値を用いた抵抗定数推定の場合、PWMによる電流高調波の影響やノイズの影響を拾い、抵抗定数推定の精度を低下させる可能性がある。これに対し、積分処理を入れることで、上記の影響を小さくでき、実施例2よりも安定した抵抗推定が可能となる。
なお、抵抗推定部25の内部で、式(1)から式(8)を用いて仮想的に算出する電流値に対しても、同様に積分処理を行うものとする。
なお、抵抗推定部25の内部で、式(1)から式(8)を用いて仮想的に算出する電流値に対しても、同様に積分処理を行うものとする。
また、積分演算を追加しても、抵抗推定する範囲は、実施例2と同様に高速域(例えば150Hz~200Hz)の範囲とする。
以上の構成により、本実施例3は、実施例2よりも安定した抵抗定数推定が可能となる。
以上の構成により、本実施例3は、実施例2よりも安定した抵抗定数推定が可能となる。
図12は、本発明の実施例4における抵抗定数推定のフローの概要を示す図である。
本実施例4は、実施例3よりも抵抗定数推定の精度を向上することを目的とする。したがって、本実施例4に係る誘導電動機の駆動装置の構成も、実施例2と同様であるが、上記目的を達成するために、実施例3において設けた積分器26の前段に、誘導電動機の電流検出値(図12では、q軸電流Iq)に対し、モータ電気時定数の時定数に設定したローパスフィルター27を設ける点で、実施例2および実施例3と異なる。
本実施例4は、実施例3よりも抵抗定数推定の精度を向上することを目的とする。したがって、本実施例4に係る誘導電動機の駆動装置の構成も、実施例2と同様であるが、上記目的を達成するために、実施例3において設けた積分器26の前段に、誘導電動機の電流検出値(図12では、q軸電流Iq)に対し、モータ電気時定数の時定数に設定したローパスフィルター27を設ける点で、実施例2および実施例3と異なる。
等価回路演算は、定常的な値を演算するものであり、過渡的な値を演算する場合には、時間軸に誘導電動機の電気時定数Tσ(=一次換算インダクタンスLσ/一次換算抵抗Rσ)の差が生じる。これに対し、電流検出値に対して、誘導電動機の電気時定数Tσによるカットオフを呈するローパスフィルター27を追加することで、等価回路演算と、実際のモータ電流検出値の時間軸の基準を同一にすることができる。
以上の構成により、本実施例4は、実施例3よりも、抵抗定数推定の精度を向上させ、トルク精度を更に向上させることができる。
図13は、本発明に係る実施例5として、先の実施例1~4に係るいずれかの誘導電動機の駆動装置を搭載する鉄道車両の一部の概略構成を示す図である。
図13において、鉄道車両に搭載する駆動装置1は、架線101から集電装置を介して電力の供給を受け、交流電力を出力して誘導電動機3に供給することで、電気エネルギーが機械トルクに変換される。誘導電動機3は、減速ギアを介して鉄道車両の車軸と連結されており、車軸に接続された車輪103とレール102間に生じる接線力により鉄道車両は走行する。
図13において、鉄道車両に搭載する駆動装置1は、架線101から集電装置を介して電力の供給を受け、交流電力を出力して誘導電動機3に供給することで、電気エネルギーが機械トルクに変換される。誘導電動機3は、減速ギアを介して鉄道車両の車軸と連結されており、車軸に接続された車輪103とレール102間に生じる接線力により鉄道車両は走行する。
本実施例5では、鉄道車両に搭載する駆動装置1として、先に説明した実施例1~4に係るいずれかの誘導電動機の駆動装置1を適用する。それにより、低速フリーラン状態からの再起動時において、誘導電動機3のトルクショックを低減できる。また、鉄道車両の低速惰行状態からの再起動時においても、鉄道車両の車体振動を低減し、乗客の乗り心地を更に改善することができる。
1…駆動装置、2…電圧出力装置、3…誘導電動機、4…電流検出器、5…座標変換器、6…座標変換器、7…初期速度推定演算部、8…初期速度推定指令発生装置、
20…制御モード切換器、21…積分器、25…抵抗推定部、26…積分器、
27…ローパスフィルター、30…ベクトル制御演算部、31…探り信号指令発生器、
32…低速判定器、101…架線、102…レール、103…車輪
20…制御モード切換器、21…積分器、25…抵抗推定部、26…積分器、
27…ローパスフィルター、30…ベクトル制御演算部、31…探り信号指令発生器、
32…低速判定器、101…架線、102…レール、103…車輪
Claims (13)
- 電圧指令値および一次周波数指令値に応じて任意の交流電圧の波形を出力して誘導電動機に印加する電圧出力装置と、
前記誘導電動機に流れる電流を検出する電流検出器と、
前記電圧出力装置を制御する制御装置と
を備える誘導電動機の駆動装置であって、
前記制御装置は、前記電流検出器が検出した電流検出値を用いて前記誘導電動機の等価回路に基づく演算を行い当該誘導電動機の抵抗定数を推定する
ことを特徴とする誘導電動機の駆動装置。 - 請求項1に記載の誘導電動機の駆動装置であって、
前記制御装置は、前記電流検出値と、前記誘導電動機の回転子速度情報と、前記電圧出力装置に対する電圧指令値および一次周波数指令値と、前記誘導電動機の等価回路定数とを用いた演算を行い前記抵抗定数を推定する
ことを特徴とする誘導電動機の駆動装置。 - 請求項1または請求項2に記載の誘導電動機の駆動装置であって、
前記制御装置が前記抵抗定数を推定する際に、前記誘導電動機に印加する前記交流電圧の一次周波数を、前記誘導電動機の回転子の正転方向に対し逆向きの負の周波数とし、かつ、該周波数の絶対値を前記回転子の周波数の駆動範囲の下限値よりも大きくする
ことを特徴とする誘導電動機の駆動装置。 - 請求項2を引用する請求項3に記載の誘導電動機の駆動装置であって、
前記駆動装置は、速度センサレスで前記誘導電動機を駆動し、
前記回転子速度情報は、前記誘導電動機の停止またはフリーランの状態から前記回転子の初期速度を推定した際の推定値であり、
前記電流検出値は、前記初期速度を推定する前に前記交流電圧を前記誘導電動機に印加して前記電流検出器が検出した値である
ことを特徴とする誘導電動機の駆動装置。 - 請求項1に記載の誘導電動機の駆動装置であって、
前記駆動装置は、速度センサレスで前記誘導電動機を駆動し、
前記制御装置は、前記誘導電動機の停止またはフリーランの状態から前記誘導電動機の回転子の初期速度を推定する演算を行うと共に、当該初期速度を推定する演算を実行する時に検出した前記電流検出値と、前記電圧出力装置に対する電圧指令値および一次周波数指令値と、前記誘導電動機の等価回路定数とを用いた演算を行い前記抵抗定数を推定する
ことを特徴とする誘導電動機の駆動装置。 - 請求項5に記載の誘導電動機の駆動装置であって、
前記初期速度を推定する際の前記誘導電動機の回転子周波数の初期値を、当該回転子周波数の駆動範囲の略最大値とすることで、前記抵抗定数を推定する対象速度範囲を低速域に限定する
ことを特徴とする誘導電動機の駆動装置。 - 請求項5または請求項6に記載の誘導電動機の駆動装置であって、
前記抵抗定数を推定する際に、前記電流検出値を積分演算した値を用いる
ことを特徴とする誘導電動機の駆動装置。 - 請求項7に記載の誘導電動機の駆動装置であって、
前記電流検出値を積分演算する前段に、前記誘導電動機の電気時定数によるカットオフを呈するローパスフィルターを備える
ことを特徴とする誘導電動機の駆動装置。 - 自らに搭載する誘導電動機を、請求項1から請求項8のいずれか1項に記載の誘導電動機の駆動装置を用いて駆動することを特徴とする鉄道車両。
- 誘導電動機に対して所定の振幅で正転方向に対して逆向きの負の周波数の交流電圧を印加し、当該印加の際に前記誘導電動機に流れる電流を検出する第1のステップと、
前記誘導電動機の初期速度推定を実行して前記誘導電動機の回転子周波数を推定すると共に、推定した当該回転子周波数と、前記電流の検出値と、前記誘導電動機を駆動制御するための電圧指令値および一次周波数指令値と、前記誘導電動機の等価回路定数とから前記誘導電動機の抵抗変動比を推定する第2のステップと、
推定した前記抵抗変動比により補正した前記誘導電動機の抵抗定数を用いて前記誘導電動機を再起動させる第3のステップと
を有する誘導電動機の駆動方法。 - 誘導電動機の初期速度推定を実行して前記誘導電動機の回転子周波数を推定すると共に、前記誘導電動機に流れる電流の検出値と、前記誘導電動機を駆動制御するための電圧指令値および一次周波数指令値と、前記誘導電動機の等価回路定数とから前記誘導電動機の抵抗変動比を推定する第1のステップと、
推定した前記回転子周波数が所定の低速域以下になった場合に、推定した前記抵抗変動比により補正した前記誘導電動機の抵抗定数を用いて前記誘導電動機を再起動させる第2のステップと
を有する誘導電動機の駆動方法。 - 請求項11に記載の誘導電動機の駆動方法であって、
前記誘導電動機に流れる電流の検出値に替えて、該検出値を積分演算した値を用いる誘導電動機の駆動方法。 - 請求項11に記載の誘導電動機の駆動方法であって、
前記誘導電動機に流れる電流の検出値に替えて、該検出値を前記誘導電動機の電気時定数によるカットオフを呈するローパスフィルターを介した後に積分演算した値を用いる誘導電動機の駆動方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017215283A JP2021040353A (ja) | 2017-11-08 | 2017-11-08 | 誘導電動機の駆動装置および駆動方法 |
JP2017-215283 | 2017-11-08 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2019093064A1 true WO2019093064A1 (ja) | 2019-05-16 |
Family
ID=66439165
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2018/037936 WO2019093064A1 (ja) | 2017-11-08 | 2018-10-11 | 誘導電動機の駆動装置および駆動方法 |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP2021040353A (ja) |
WO (1) | WO2019093064A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20230353075A1 (en) * | 2021-05-20 | 2023-11-02 | Mitsubishi Electric Corporation | Rotary machine control apparatus, machine learning apparatus, and inference apparatus |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007049843A (ja) * | 2005-08-11 | 2007-02-22 | Hitachi Ltd | 永久磁石同期モータのベクトル制御装置 |
JP2011024281A (ja) * | 2009-07-13 | 2011-02-03 | Fuji Electric Systems Co Ltd | 永久磁石形同期電動機の制御装置 |
WO2011064846A1 (ja) * | 2009-11-25 | 2011-06-03 | 三菱電機株式会社 | 電気車の電力変換装置 |
-
2017
- 2017-11-08 JP JP2017215283A patent/JP2021040353A/ja active Pending
-
2018
- 2018-10-11 WO PCT/JP2018/037936 patent/WO2019093064A1/ja active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007049843A (ja) * | 2005-08-11 | 2007-02-22 | Hitachi Ltd | 永久磁石同期モータのベクトル制御装置 |
JP2011024281A (ja) * | 2009-07-13 | 2011-02-03 | Fuji Electric Systems Co Ltd | 永久磁石形同期電動機の制御装置 |
WO2011064846A1 (ja) * | 2009-11-25 | 2011-06-03 | 三菱電機株式会社 | 電気車の電力変換装置 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20230353075A1 (en) * | 2021-05-20 | 2023-11-02 | Mitsubishi Electric Corporation | Rotary machine control apparatus, machine learning apparatus, and inference apparatus |
US11929694B2 (en) * | 2021-05-20 | 2024-03-12 | Mitsubishi Electric Corporation | Rotary machine control apparatus, machine learning apparatus, and inference apparatus |
Also Published As
Publication number | Publication date |
---|---|
JP2021040353A (ja) | 2021-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6978602B2 (ja) | 永久磁石同期電動機の駆動装置、駆動システムおよび駆動方法 | |
JP4082444B1 (ja) | 永久磁石同期電動機のベクトル制御装置 | |
JP4205157B1 (ja) | 電動機の制御装置 | |
US9054630B2 (en) | Synchronous machine controller | |
JP6868772B2 (ja) | モータ制御装置およびモータ制御方法 | |
JP4578700B2 (ja) | ブラシレスdcモータの制御装置 | |
WO2017033508A1 (ja) | ドライブシステムおよびインバータ装置 | |
US9935568B2 (en) | Control apparatus of rotary electric machine | |
US20190305713A1 (en) | Method of controlling motor and device of controlling motor | |
JP6762825B2 (ja) | 電動機制御装置およびドライブシステム | |
JP6536473B2 (ja) | 回転電機の制御装置 | |
JP5109790B2 (ja) | 永久磁石形同期電動機の制御装置 | |
JP7361924B2 (ja) | モータ制御装置、モータ制御方法 | |
JP6287715B2 (ja) | 回転機の制御装置 | |
WO2019093064A1 (ja) | 誘導電動機の駆動装置および駆動方法 | |
JP2013141345A (ja) | モータ制御装置及び空気調和機 | |
WO2018179822A1 (ja) | 誘導電動機の駆動装置 | |
JP2023048833A (ja) | モータユニットの状態推定方法及び状態推定装置 | |
JP2012114974A (ja) | 電気車制御装置 | |
JP2017205017A (ja) | 空気調和機のモータ制御装置及び空気調和機 | |
JP7245739B2 (ja) | 誘導電動機の駆動装置、駆動方法、および電気車 | |
JP5996485B2 (ja) | モータの駆動制御装置 | |
WO2019008838A1 (ja) | 誘導電動機の駆動装置及び駆動方法 | |
JP2021002949A (ja) | 永久磁石同期電動機の駆動装置、駆動方法および電気車 | |
US11545924B2 (en) | Power conversion device, power conversion method and program |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 18877003 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 18877003 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: JP |