WO2021038670A1 - 推定装置および交流電動機の駆動装置 - Google Patents

推定装置および交流電動機の駆動装置 Download PDF

Info

Publication number
WO2021038670A1
WO2021038670A1 PCT/JP2019/033182 JP2019033182W WO2021038670A1 WO 2021038670 A1 WO2021038670 A1 WO 2021038670A1 JP 2019033182 W JP2019033182 W JP 2019033182W WO 2021038670 A1 WO2021038670 A1 WO 2021038670A1
Authority
WO
WIPO (PCT)
Prior art keywords
time
current
phase
voltage
unit
Prior art date
Application number
PCT/JP2019/033182
Other languages
English (en)
French (fr)
Inventor
健治 ▲高▼橋
慎也 豊留
鹿嶋 美津夫
知宏 沓木
Original Assignee
三菱電機株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 三菱電機株式会社 filed Critical 三菱電機株式会社
Priority to US17/625,096 priority Critical patent/US11705841B2/en
Priority to JP2021541799A priority patent/JP7130143B2/ja
Priority to CN201980099465.8A priority patent/CN114270695B/zh
Priority to DE112019007652.3T priority patent/DE112019007652T5/de
Priority to PCT/JP2019/033182 priority patent/WO2021038670A1/ja
Publication of WO2021038670A1 publication Critical patent/WO2021038670A1/ja

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P23/00Arrangements or methods for the control of AC motors characterised by a control method other than vector control
    • H02P23/14Estimation or adaptation of motor parameters, e.g. rotor time constant, flux, speed, current or voltage
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • H02P21/18Estimation of position or speed
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/24Vector control not involving the use of rotor position or rotor speed sensors
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/14Electronic commutators
    • H02P6/16Circuit arrangements for detecting position
    • H02P6/18Circuit arrangements for detecting position without separate position detecting elements
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/34Modelling or simulation for control purposes

Definitions

  • the present invention relates to an estimation device that estimates at least one of the position and speed of an AC motor and a drive device for the AC motor.
  • Information indicating the position of the rotor is used to drive AC motors such as induction machines and synchronous machines.
  • Using a position sensor or speed sensor to obtain information indicating the position of the rotor has problems such as increased manufacturing costs. Therefore, with respect to the drive device of the AC motor, many studies have been made on position sensorless control without using a position sensor or a speed sensor.
  • the 1-shunt current detection method is a method of measuring the phase current flowing through an AC motor by using a current sensor provided on the DC bus of the inverter.
  • the method of measuring the phase current using the current sensor provided on the DC bus is called the 1-shunt current detection method because the shunt resistance is often used as the current sensor, but the shunt resistance is used as the current sensor provided on the DC bus.
  • it is generally called a one-shunt current detection method.
  • the one-shunt current detection method is also known to use a current sensor different from the shunt resistor, such as a current transformer called CT (Current Transformer).
  • CT Current Transformer
  • Patent Document 1 discloses a technique for driving an AC motor by using a position sensorless control and a one-shunt current detection method in combination.
  • the current of each phase cannot be detected at the apex of the carrier signal at the same time. Therefore, in the technique described in Patent Document 1, the phase current obtained from the generatrix current is interpolated between the latter half of the first carrier cycle and the first half of the second carrier cycle following the first carrier cycle.
  • the phase current at the apex of the carrier signal which is the timing of the boundary between the carrier period of 1 and the carrier period of the second, is calculated.
  • a voltage command is generated based on the phase current at the apex of the carrier signal, and information indicating the rotational state of the rotor of the AC motor is calculated based on the generated voltage command.
  • the information indicating the rotational state of the rotor is at least one of the position and speed of the rotor.
  • the technique described in Patent Document 1 uses the bus currents detected in the latter half of the first carrier cycle and the first half of the second carrier cycle following the first carrier cycle, respectively.
  • the detection cycle of the phase current used for calculating the rotational state of the rotor of the AC motor is twice or more the carrier cycle.
  • the switching loss of the inverter increases as the carrier frequency rises, there is generally an upper limit on the carrier frequency due to the cooling performance of the drive device of the AC motor or the convenience of power efficiency.
  • the control cycle becomes longer with respect to the rotation frequency of the AC motor, and it may be difficult to secure the estimation accuracy of the rotation state of the rotor included in the AC motor.
  • the present invention has been made in view of the above, and an object of the present invention is to obtain an estimation device capable of improving the estimation accuracy of at least one of the position and speed of an AC motor.
  • the estimation device of the present invention includes a bus current detection unit, a phase current determination unit, a time difference calculation unit, and an estimation unit.
  • the bus current detection unit detects the value of the bus current, which is the current flowing through the DC bus of the voltage inverter driven by a plurality of gate pulse signals generated based on the voltage command.
  • the phase current determination unit determines the values of a plurality of phase currents supplied from the voltage inverter to the AC motor based on the value of the bus current detected by the bus current detection unit and the states of the plurality of gate pulse signals. Judgment processing is performed.
  • the time difference calculation unit is the first current detection time specified as the detection time of the bus current used in the previous judgment process by the phase current determination unit and the bus current used in the current determination process by the phase current determination unit.
  • the time difference from the second current detection time specified as the detection time of is calculated.
  • the estimation unit estimates at least one of the position and speed of the AC motor based on the values of the plurality of phase currents determined by the phase current determination unit and the time difference calculated by the time difference calculation unit.
  • the figure which shows an example of the change of the voltage command of an intermediate phase when the voltage command of three phases which concerns on Embodiment 1 is a sine wave.
  • the figure which shows the image of the integral operation by the rectangle approximation The figure which shows the image of the integral operation by the rectangle approximation
  • the figure which shows the image of the integral operation when the fluctuation of the time difference is ignored.
  • Diagram to explain the error of differential operation The figure which shows an example of the configuration example of the electric motor system including the estimation device which concerns on Embodiment 2 of this invention.
  • the figure which shows an example of the result of speed estimation when the time difference between the current detection times concerning Embodiment 3 is a fixed value.
  • FIG. 1 is a diagram showing a configuration example of an electric motor system including an estimation device according to a first embodiment of the present invention.
  • the electric motor system 100 according to the first embodiment includes an AC electric motor 1 and a drive device 2 for driving the AC electric motor 1.
  • the AC electric motor 1 is a permanent magnet synchronous motor in which a permanent magnet is provided on the rotor, but it may be a winding field type synchronous motor in which a field winding is wound around the rotor, and the rotor bumps. It may be a relaxation type synchronous motor that obtains rotational torque by utilizing the polarity. Further, the arrangement of the permanent magnets in the AC motor 1 may be an embedded type arrangement or a surface type arrangement. Further, although an example in which the AC motor 1 is a three-phase AC motor will be described here, the AC motor 1 may be an AC motor other than the three-phase one. For example, the AC motor 1 may be a two-phase AC motor or a five-phase AC motor.
  • the drive device 2 converts the DC voltage supplied from the DC power supply 3 into an AC voltage and outputs the converted AC voltage to the AC electric motor 1, and controls the voltage type inverter 4 to control the AC electric motor 1.
  • a control device 5 for driving is provided.
  • the voltage type inverter 4 includes a main circuit 6 and a gate driver 7.
  • the main circuit 6 includes a plurality of switching elements Q1, Q2, Q3, Q4, Q5, Q6.
  • one end of the switching elements Q1, Q3 and Q5 is connected to the DC bus 61 on the high potential side, and one end of the switching elements Q2, Q4 and Q6 is connected to the DC bus 61 on the low potential side.
  • the other ends of the switching element Q1 and the switching element Q2 are connected to each other to form a U-phase leg.
  • the other ends of the switching element Q3 and the switching element Q4 are connected to each other to form a V-phase leg.
  • the other ends of the switching element Q5 and the switching element Q6 are connected to each other to form a W-phase leg.
  • the voltage type inverter 4 includes a three-phase bridge circuit including a U-phase leg, a V-phase leg, and a W-phase leg.
  • the voltage type inverter 4 can output an AC voltage having an arbitrary amplitude and an arbitrary frequency by switching the plurality of switching elements Q1, Q2, Q3, Q4, Q5, and Q6 on and off.
  • switching elements Q1, Q2, Q3, Q4, Q5, and Q6 when each of the switching elements Q1, Q2, Q3, Q4, Q5, and Q6 is shown without distinction, it may be described as the switching element Q.
  • Each switching element Q is an IGBT (Insulated Gate Bipolar Transistor) with a built-in anti-parallel diode, but may be a MOSFET (Metal-Oxide-Semiconductor Field-Effect Transistor) with a built-in anti-parallel diode.
  • MOSFET Metal-Oxide-Semiconductor Field-Effect Transistor
  • the gate driver 7, the gate pulse signal G u outputted from the control unit 5, G v, a gate pulse signal by inverting amplifying G w G un, G vn, generates G wn, generated G un, G vn and G wn are output to the gates of the switching elements Q2, Q4 and Q6.
  • the gate pulse signals G u , G v , G w , G up , Gun , G vp , G vn , G wp , G wn are on when the high potential level is in the state, and are on at the low potential level. If it is in the state, it is in the off state.
  • the gate driver 7 has a function of insulating the low-voltage system control device 5 and the high-voltage system main circuit 6, and has a role of preventing a failure of the control device 5 when the main circuit 6 is abnormal. ..
  • the gate pulse signal G when each of the gate pulse signals G u , G v , and G w is shown without distinction, it may be described as the gate pulse signal G.
  • the control device 5 includes a zero vector modulation unit 34 that performs zero vector modulation, a gate pulse generation unit 35 that generates gate pulse signals Gu , G v , and G w , a magnetic pole position ⁇ e of the rotor of the AC motor 1, and It is provided with an estimation device 9 for estimating the rotation speed ⁇ e.
  • the magnetic pole position ⁇ e is the electric angle of the rotor included in the AC motor 1, and is an example of the position of the AC motor 1.
  • the rotation speed ⁇ e is the electric angular velocity of the rotor included in the AC motor 1, and is an example of the speed of the AC motor 1.
  • the estimation of the magnetic pole position ⁇ e may be described as position estimation, and the estimation of the rotational speed ⁇ e may be described as velocity estimation.
  • the zero vector modulation unit 34 irregularly changes the output ratio of the two types of zero voltage vectors output from the voltage type inverter 4.
  • the two types of zero voltage vectors are a first zero voltage vector and a second zero voltage vector.
  • the first zero voltage vector is output from the voltage inverter 4 with the switching elements Q1, Q3, and Q5, which are the upper arms, all on.
  • the second zero voltage vector is output from the voltage inverter 4 in a state where the switching elements Q1, Q3, and Q5, which are the upper arms, are all off.
  • the carrier noise is noise generated by vibration of the AC electric motor 1 or the voltage inverter 4 or the like according to the carrier frequency which is the frequency of the carrier wave Sc described later.
  • the output ratio of the first zero voltage vector and the second zero voltage vector changes depending on the modulation method by the gate pulse generation unit 35. Even if the modulation method changes, the output ratio between the first zero voltage vector and the second zero voltage vector can be changed by adding the same value to the three-phase voltage commands v u * , v v * , and v w *. Can be made to.
  • the voltage command v u * is a u-phase voltage command
  • the voltage command v v * is a v-phase voltage command
  • the voltage command v w * is a w-phase voltage command.
  • the time for all gate pulse signals Gu , G v , and G w to turn on is extended, and the voltage type inverter The time from which the first zero voltage vector is output from the fourth is extended.
  • the time during which the gate pulse signals Gu , G v , and G w are all turned off is extended, and the voltage is increased. The time for which the second zero voltage vector is output from the type inverter 4 is extended.
  • the zero vector modulator 34 adds a random value to the three-phase voltage commands v u * , v v * , v w * to make the output ratio between the first zero voltage vector and the second zero voltage vector inconsistent. Performs zero vector modulation that changes to a rule. If the zero vector modulation unit 34 does not satisfy the preset conditions, the input three-phase voltage commands v u * , v v * , and v w * may be output to the gate pulse generation unit 35 as they are. it can.
  • the preset condition is, for example, that the AC motor 1 is rotating at a speed equal to or lower than a preset speed. Further, when the zero vector modulation unit 34 is set not to perform zero vector modulation, the input three-phase voltage commands v u * , v v * , v w * are output to the gate pulse generation unit 35 as they are. can do.
  • the gate pulse generation unit 35 generates gate pulse signals Gu , G v , G w based on the voltage commands v u * , v v * , v w * .
  • the gate pulse generation unit 35 outputs the generated gate pulse signals Gu , G v , and G w to the gate driver 7 of the voltage type inverter 4.
  • the gate pulse generation unit 35 is a comparison unit 21 that compares the carrier wave Sc, which is a high-frequency periodic signal, with the three-phase voltage commands v u * , v v * , v w *, and the gate output from the comparison unit 21.
  • a pulse shift processing unit 22 that shifts the pulse signals Gu , G v , and G w is provided.
  • the carrier wave Sc is a triangular wave signal
  • the carrier wave Sc may be a signal having another waveform such as a sawtooth wave.
  • the carrier wave Sc is also called a carrier signal.
  • Comparing section 21 when the instantaneous value of the voltage command v u * is equal to or less than the instantaneous value of the carrier wave Sc, and turns on the gate pulse signal G u, the instantaneous value of the voltage command v u * is than the instantaneous value of the carrier wave Sc If is also large, turn off the gate pulse signal Gu.
  • the comparison unit 21 turns on the gate pulse signal G v, and the instantaneous value of the voltage command v v * is from the instantaneous value of the carrier wave Sc. If is also large, the gate pulse signal Gv is turned off.
  • the comparison unit 21 turns on the gate pulse signal G w, and the instantaneous value of the voltage command v w * is from the instantaneous value of the carrier wave Sc. If is also large, the gate pulse signal G w is turned off.
  • the gate pulse generation unit 35 generates the gate pulse signals Gu , G v , and G w by using the carrier comparison modulation method, but generates the gate pulse signals G u , G v , and G w .
  • the method is not limited to the carrier comparison modulation method.
  • the gate pulse generation unit 35 may generate gate pulse signals Gu , G v , and G w by using another modulation method such as a space vector modulation method instead of the carrier comparison modulation method.
  • the gate pulse generator 35, the gate pulse signal G u, G v, instead of G w, the gate pulse signal G up, G un, G vp , G vn, G wp, may generate G wn ..
  • the gate driver 7 the gate pulse signal G up, G un, G vp , G vn, G wp, to amplify the G wn.
  • the gate driver 7, amplified gate pulse signal G up, G un, G vp , G vn, G wp, G wn switching elements Q1, Q2, Q3, Q4, Q5, and outputs Q6 to the gate of the.
  • the pulse shift processing unit 22 performs pulse shift processing for shifting the switching timing of the gate pulse signals Gu , G v , and G w from on and off.
  • the pulse shift processing of the pulse shift processing unit 22 will be described in detail later.
  • the estimation device 9 determines the magnetic pole position ⁇ e and the rotation speed ⁇ e of the AC motor 1 based on the gate pulse signals Gu , G v , G w and the voltage commands v u * , v v * , v w *. presume.
  • the estimation device 9 can also estimate only one of the magnetic pole position ⁇ e and the rotation speed ⁇ e of the AC motor 1.
  • Estimating unit 9 the bus current detection unit 11 for detecting a value of the bus current i bus, phase currents i u, i v, and determines phase current determination unit 12 the value of i w, the gate pulse signal G u, G
  • An estimation unit 15 for estimating the magnetic pole position ⁇ e and the rotation speed ⁇ e based on v v * and v w * is provided.
  • the value of the bus current i bus detected by the bus current detection unit 11 is the instantaneous value of the bus current is a current flowing through the DC bus 61 between the DC power supply 3 and the voltage source inverter 4.
  • the bus current detection unit 11 may be a current sensor of a type using a current transformer called CT, or a current sensor of a type using a shunt resistor. In the following, sometimes simply described as the bus current i bus the value of the bus current i bus.
  • the bus current detection unit 11 is provided on the DC bus 61 on the low potential side and detects the value of the current flowing through the DC bus 61 on the low potential side, but on the DC bus 61 on the high potential side. It may be provided.
  • the bus current detection unit 11 is a type of current sensor that uses a shunt resistor, the component cost of the insulation circuit in the bus current detection unit 11 can be reduced by providing the bus current detection unit 11 on the DC bus 61 on the low potential side. It can be suppressed.
  • Phase current determination unit 12 determines, based on the value of the bus current i bus detected by the bus current detection unit 11, the phase currents i u, i v, the value of i w.
  • the values of the phase currents i u , iv , and i w determined by the phase current determination unit 12 are the u-phase, v-phase, and w-phase currents flowing between the voltage-type inverter 4 and the AC motor 1. It is the instantaneous value of one phase current.
  • phase current i u described with u-phase current i u
  • phase current i v is described as v-phase current i v
  • phase currents i u, i v, i simply phase currents i u values of w, i v, i w.
  • FIGS. 2 to 4 are diagrams for explaining a method of determining the value of the phase current by the phase current determining unit according to the first embodiment.
  • a Y-connected three-phase resistor load 1a is connected to the voltage inverter 4.
  • the current flowing from the connection points of the switching elements Q1 and Q2 to the three-phase resistance load 1a is the u-phase current iu
  • the current flowing from the connection points of the switching elements Q3 and Q4 to the three-phase resistance load 1a is the v-phase current i. v
  • the current flowing from the connection point of the switching elements Q5, Q6 to the three-phase resistor loads 1a is w-phase current i w.
  • the direction of the phase current flowing from the voltage type inverter 4 to the three-phase resistance load 1a is the positive direction
  • the direction of the phase current flowing from the three-phase resistance load 1a to the voltage type inverter 4 is the negative direction.
  • the direction of the u-phase current i u is the positive direction.
  • the direction of the V-phase current i v flows in the direction of the arrow shown in FIG. 2
  • the direction of the V-phase current i v is positive
  • the w-phase current i w flows in the direction of the arrow shown in FIG. 2
  • the direction of the w-phase current i w is the positive direction.
  • Phase current determination unit 12 compares the value of the bus current i bus detected by the bus current detection unit 11, a plurality of gate pulse signal G u, G v, based on the state of the G w, the phase currents i u, i Determine the values of v and i w.
  • the phase current determining unit 12 based on the six values of the bus current i bus the combination patterns of the switching elements Q on and off is detected by the bus current detection unit 11 at the timing is a specific combination pattern, Determine the values of the phase currents i u , iv , and i w.
  • a u-phase current i u flows from the voltage inverter 4 to the three-phase resistance load 1a, and currents having the same magnitude as the u-phase current i u are the v-phase current iv and the w-phase current. It is divided into i w and flows from the three-phase resistance load 1a to the voltage inverter 4.
  • Current flowing through the DC bus 61 is the same magnitude of current and the U-phase current i u, the phase current determination unit 12, in the state shown in FIG. 3, the bus current i bus detected by the bus current detection unit 11
  • the value of the u-phase current iu can be determined from the value of.
  • the bus current detection unit 11 the current flowing through the DC bus 61 is arranged to the DC bus 61 as in the case of orientation shown in FIG. 3, and outputs a positive bus current i bus.
  • FIG. 3 there is shown an example of determining the value of the positive direction of the U-phase current i u, the value of the negative direction of the U-phase current i u, the positive direction and value of the negative direction of the v-phase current i v, and the value of the positive and negative directions of the w-phase current i w are also determined similarly.
  • the switching elements Q1, Q3 and Q6 are on and the switching elements Q2, Q4 and Q5 are off.
  • the phase current determination unit 12 can determine the value of the negative direction of the w-phase current i w by the value of the bus current i bus detected by the bus current detection unit 11.
  • phase current determining unit 12 determines by calculation the value of the remaining one phase current from the values of the two determined phase currents. In this way, the phase current determination unit 12 determines the values of the three phase currents i u , iv , and i w, with the two bus current detections as one set.
  • the voltage type inverter 4 is a two-level inverter
  • Each of the six combination patterns out of the eight combination patterns is the specific pattern described above. If six combinations of patterns of on and off of the switching element Q is a specific pattern, the bus current i bus of the phase current from the values i u is detected by the bus current detection unit 11, i v, either i w One value can be determined.
  • the remaining two combination pattern since the output of the voltage-source inverter 4 becomes zero voltage vector as described above, the bus current i values phase currents from i u of bus detected by the bus current detection unit 11, i v, i w It is difficult to determine the value of.
  • FIG. 5 is a diagram showing an example of the relationship between the detection timing of the bus current used for the determination of the phase current by the phase current determination unit according to the first embodiment, the carrier wave, the voltage command, and the gate pulse signal. ..
  • Tc is the carrier period which is the period of the carrier wave Sc
  • fc is the carrier frequency which is the frequency of the carrier wave Sc.
  • the instantaneous value of the voltage command v u * is the largest, and the instantaneous value of the voltage command v v * is next. It is large and the instantaneous value of the voltage command v w * is the smallest.
  • the phase of the voltage command whose absolute value is in the middle among the three phases of the voltage commands v u * , v v * , and v w * at a certain time is referred to as an intermediate phase.
  • the intermediate phase is the V phase.
  • the period from time t1 to t7 is the first downlink half cycle of the carrier wave Sc.
  • the downlink half cycle is a half cycle of the carrier wave Sc, and is a period in which the value of the carrier wave Sc gradually decreases.
  • the switching elements Q1, Q4 and Q6 are on and the switching element Q2.
  • Q3 and Q5 are off.
  • the value of the bus current i bus detected by the bus current detection unit 11 is the same as the value of u-phase current i u.
  • Phase current determination unit 12 determines the value of the bus current i bus detected by the bus current detection unit 11 at time t3 as the value of u-phase current i u.
  • the switching elements Q1, Q3 and Q6 are on and the switching element Q2. , Q4 and Q5 are off.
  • the value of the bus current i bus detected by the bus current detection unit 11 is the same as the value of w-phase current i w.
  • Phase current determination unit 12 determines the value of the bus current i bus detected by the bus current detection unit 11 at time t5 as the value of w-phase current i w.
  • Phase current determination unit 12 the u-phase current i u obtained from the value of the detected bus current i bus at time t3 values and obtained from the value of the detected bus current i bus at time t5 w-phase current i and a value of w, based on the Kirchhoff's law, to calculate the value of the v-phase current i v.
  • the phase current determining unit 12 the bus current i values phase currents from i u of bus detected by the bus current detection unit 11 at a plurality of current detection timing in the downlink half cycle of the carrier wave Sc, i v, i Determine the value of w.
  • phase current determining unit 12 the value of the bus current i bus detected by the bus current detection unit 11 a time t13, t15 of the period from the time t11 which is the next downstream half period t16 as the current detection timing It is used to determine the values of the phase currents i u , iv , i w.
  • the phase current determination unit 12 the downlink half cycle of the carrier wave Sc, 2 times of current detection timing in the bus current i bus of the phase current from the values i u is detected by the bus current detection unit 11, i v , it determines the value of i w.
  • the timing of these two current detections is not limited to the downlink half cycle of the carrier wave Sc.
  • the phase current determination unit 12, instead of the downstream half period of the carrier wave Sc, the upstream half-period of the carrier wave Sc, the bus current i bus detected by the bus current detection unit 11 with two current detection timing It is also possible to determine the values of the phase currents i u , iv , and i w from the values.
  • the upstream half cycle is a half cycle of the carrier wave Sc, and is a period in which the value of the carrier wave Sc gradually increases.
  • the phase current determination unit 12 determines waits a predetermined time as the time until the ringing subsides, the value of the phase current with the value of the bus current i bus detected by the bus current detection unit 11 ..
  • the phase current determination unit 12 is detected by the bus current detection unit 11 at time t3 is the gate pulse signal of an intermediate phase gate pulse signal G v is immediately before the time t4 from OFF to ON
  • the value of the bus current i bus is determined as the value of the phase current i u .
  • the phase current determining unit 12 the value of the bus current i bus gate pulse signal G v is detected by the bus current detection unit 11 at time t5 after a preset time TA from the time t4 which is turned from OFF Determined as the value of the phase current i w.
  • the phase current determination unit 12 determines the state switching timing, which is the timing at which the gate pulse signal of the intermediate phase is switched on and off among the gate pulse signals Gu , G v , and G w. Then, the phase current determination unit 12 determines the timing immediately before the state switching timing and the timing after the time TA from the state switching timing as the bus current detection timing, which is the timing for detecting the value of the bus current ibus, respectively. Phase current determination unit 12, these in each of the two bus current detection timing from the value of the bus current i bus detected by the bus current detection unit 11, the phase currents i u, i v, any two of i w Determine the value of the phase current.
  • the phase current determination unit 12 calculates the values of the remaining phase currents based on the determined values of the two phase currents. Since the determination timings of the two phase currents can be brought close to each other, the phase current determination unit 12 can accurately determine the values of the remaining phase currents based on the detected values of the two phase currents.
  • the phase current determination unit 12 can store the value of the bus current ibus that is repeatedly detected by the bus current detection unit 11. Phase current determination unit 12, among the values of the stored bus current i bus, it is possible to extract the value of the bus current i bus detected by the bus current detection unit 11 just before the gate pulse signal of an intermediate phase is switched. The phase current determination unit 12 can determine the value of the phase current based on the value of the extracted bus current ibus.
  • phase current determination unit 12 can also use a value of the detected bus current i bus on and off and is switched time of the gate pulse signal of the intermediate phase by the bus current detection unit 11.
  • Phase current determining unit 12 in addition to the value of the bus current i bus to the on and off gate pulse signal of an intermediate phase is detected in the switching time, since switching the gate pulse signal on and off of the intermediate phase using the value of the bus current i bus detected time after TC.
  • the time TC is the time obtained by adding the delay time Td to the time TA.
  • the phase current determination unit 12 determines the two times of the bus current detection timing based on the timing of the gate pulse signal of an intermediate phase is changed among a plurality of gate pulse signal G u, G v, G w .
  • These two bus current detection timings are timings for detecting the value of the bus current ibus used for determining the value of the two phase currents of the phase currents i u , iv , and i w.
  • Phase current determination unit 12 based on the values of a plurality of bus current i bus detected by the bus current detection unit 11 for each detection period including two bus current detection timing determined phase currents i u, i v , it determines the value of i w.
  • the above-mentioned eight combination patterns are switched by the change of the voltage command v u * , v v * , v w * , and the switching interval of the combination pattern is also changed by the change of the voltage command v u * , v v * , v w *. ..
  • Combination pattern switching of the gate pulse signal G u of three-phase, G v performed by the switched and any of the on and off of G w, ringing appears on the bus current i bus immediately thereafter.
  • Ringing generally converges on the order of several microseconds, but if any of the three-phase gate pulse signals Gu , G v , and G w is switched on and off again before the ringing converges, the state before switching is achieved. It is difficult to obtain the value of the phase current in. That is, when the state change of the gate pulse signal group composed of the three-phase gate pulse signals Gu , G v , and G w occurs twice in the range of several microseconds, it is detected by the bus current detection unit 11. It is difficult to determine the value of the phase current from the value of the bus current ibus. Thus, switching intervals of the combination patterns is shorter than the convergence time of ringing, it is difficult to determine the value of the phase current from the value of the bus current i bus detected by the bus current detection unit 11.
  • the pulse shift processing unit 22 of the control device 5 turns on at least one of the three-phase gate pulse signals Gu , G v , and G w so that the switching interval of the combination pattern becomes equal to or longer than the preset time TA. Performs processing to shift the timing of switching between and off.
  • the time TA is set to a value equal to or longer than the time from the occurrence of ringing to the convergence of ringing, that is, the convergence time of ringing.
  • the pulse shift processing unit 22 can shift the switching interval of the combination pattern by the ringing convergence time or more, until the ringing converges regardless of the values of the voltage commands v u * , v v * , and v w *. Waiting time can be secured.
  • the pulse shift processing unit 22 performs pulse shift processing for shifting the timing at which at least one of the three-phase gate pulse signals Gu , G v , and G w is switched on and off within one carrier cycle. For example, the pulse shift processing unit 22 shifts the gate pulse signal G in the period from the ascending half cycle to the descending half cycle of the carrier wave Sc, and the gate pulse signal G in the ascending half cycle and the descending half cycle of the carrier wave Sc, respectively. Process to change the duty ratio of.
  • the pulse shift processing unit 22 can suppress a change in the output voltage of the three phases of the voltage inverter 4 in units of one carrier cycle.
  • the phase current determination unit 12 determines the phase current i from the value of the bus current i bus detected by the bus current detection unit 11 regardless of the values of the voltage commands v u * , v v * , and v w * by the pulse shift process. it is possible to determine u, i v, the value of i w.
  • the phase current determination unit 12 determines the values of the three phase currents i u , iv , and i w in the downlink half cycle of the carrier wave Sc, but the three phase currents i u , i v, timing determine the value of i w is not limited to the downstream half period of the carrier wave Sc.
  • the phase current determination unit 12 instead of the downstream half period of the carrier wave Sc, the upstream half-period of the carrier wave Sc, 2 times of the bus current i bus detected by the bus current detection unit 11 by the bus current detection timing
  • the values of the phase currents i u , iv , and i w can be determined based on the values of.
  • the phase current determination unit 12 performs determination processing for determining the values of the phase currents i u , iv , and i w in the downlink half cycle of the carrier wave Sc, and the phase current i u in the uplink half cycle of the carrier wave Sc. , i v, and a determination processing for determining the value of i w can be performed alternately.
  • the phase current determination unit 12 since the pulse shift processing unit 22 suppresses the change in the output voltage of the voltage inverter 4 in units of one carrier cycle, the phase current determination unit 12 has a cycle of 1.5 times or more the cycle of the carrier wave Sc. Once every time, the values of the phase currents i u , iv , and i w are determined.
  • the phase current determination unit 12 determines the values of the three phase currents i u , iv , and i w for each cycle of the carrier wave Sc, and instead of determining the values of the carrier wave Sc for every two or more cycles. It is also possible to determine the values of the three phase currents i u , iv , and i w.
  • the frequency of the carrier wave Sc has an upper limit
  • the three-phase current determination cycle which is the cycle of the determination process for determining the values of the three phase currents i u , iv , and i w, becomes long, the AC motor 1 is rotating at high speed.
  • the frequencies of the phase currents i u , iv , i w approach the frequency of the carrier wave Sc.
  • the time resolution for one cycle of the phase currents i u , iv , i w determined by the phase current determination unit 12 becomes low.
  • the estimation accuracy of the estimation device 9 and the control performance by the control device 5 decrease. Therefore, it is preferable that the three-phase current determination cycle is not long at the time of high-speed rotation. Therefore, when there is an upper limit to the frequency of the carrier wave Sc, it may be desirable that the three-phase current determination cycle is one cycle of the carrier wave Sc.
  • the state switching timing of the gate pulse signal of the intermediate phase described above changes according to the voltage command of the intermediate phase, and the bus current detection timing which is the detection timing of the value of the bus current ibus depends on the state switching timing of the gate pulse signal of the intermediate phase. Changes. Further, the bus current detection timing is also changed by the pulse shift processing by the pulse shift processing unit 22. The change in the bus current detection timing due to the pulse shift processing is small, but when the voltage command of the intermediate phase fluctuates greatly, the bus current detection timing also fluctuates greatly.
  • the value of the voltage command v v * of the v phase which is the voltage command of the intermediate phase, changes significantly at the time t11 when the carrier wave Sc reaches the maximum value. Therefore, the timing of a gate pulse signal of an intermediate phase gate pulse signal G v is switched also be around greater. Along with this, the time from the time t1 when the carrier wave Sc reaches the maximum value to the time t3 and t5 which is the bus current detection timing, and the time t13 and t15 which are the bus current detection timings from the time t11 when the carrier wave Sc reaches the maximum value. The difference from the time until is large.
  • the bus current detection timing changes for each carrier cycle. Therefore, when the estimation unit 15 estimates the magnetic pole position ⁇ e and the rotation speed ⁇ e using the values of the phase currents i u , iv , and i w determined by the phase current determination unit 12, the current detection time interval The time difference Tb of is used.
  • the current detection time is a time specified as a detection time by the bus current detection unit 11 of the bus current ibus used for the determination process by the phase current determination unit 12. For example, as for the current detection time, the detection time of the bus current ibus used by the phase current determination unit 12 to determine the current value of the first phase and the phase current determination unit 12 determine the current value of the second phase.
  • the phase current determination unit 12 detects the detection time of the bus current ibus used to determine the current value of the first phase and the bus current i used to determine the current value of the second phase. The time obtained by adding the bus detection time and dividing by 2 is obtained as the current detection time. Detection time of the bus current i bus is the time at which bus current i bus is detected by bus current detection unit 11.
  • the current detection time may be the detection time by the bus current detection unit 11 of the bus current ibus used to determine the current value of the first phase, and the current value of the second phase is determined. It may be the detection time by the bus current detection unit 11 of the bus current ibus used for this purpose.
  • the bus current i bus used to determine the value of u-phase current i u is detected, then in order to determine the value of u-phase current i u bus current i bus is the time tu2 the time t13 is detected to be used in.
  • the time t15 when the current ibus is detected be the time tw2.
  • the time TA which is the waiting time for ringing convergence, is known.
  • the time difference between the time tu1 and the time tw1 is known, and the time difference between the time tu2 and the time tw2 is also known.
  • the time difference between the time tu1 and the time tu2 or the time difference between the time tw1 and the time tw2 is not known. That is, the time difference Tb between the current detection times is unknown.
  • the time difference calculation unit 13 calculates the time difference Tb between the current detection times based on the gate pulse signals Gu , G v , and G w. Specifically, the time difference calculating section 13, the gate pulse signal G u, G v, based on the timing of the gate pulse signal of an intermediate phase is switched between on and off of the G w, time tu1, tw1, Determine tu2 and tw2. The time difference calculation unit 13 calculates time tabg1, which is an intermediate time between time tu1 and time tw1, based on the determined time tu1, tw1, tu2, and tw2, and at a time intermediate between time tu2 and time tw2. Calculate tavg2 at a certain time.
  • the time difference calculation unit 13 calculates the time difference between the time tabg1 and the time tabg2 as the time difference Tb. As described above, the time difference calculation unit 13 is used in the current detection time specified as the detection time of the bus current ibus used in the determination process by the phase current determination unit 12 and in the next determination process by the phase current determination unit 12. calculating a time difference Tb between the current detection time specified as the detection time of the bus current i bus that is.
  • the time difference between the time tu1 and the time tu2 and the time difference between the time tw1 and the time tw2 are the time differences between the time tabg1 and the time tabg2, respectively.
  • the time difference calculation unit 13 can calculate the time difference between the time tu1 and the time tu2 as the time difference Tb, and can also calculate the time difference between the time tw1 and the time tw2 as the time difference Tb.
  • FIG. 5 shows an example in which the time difference between the time tw1 and the time tw2 is calculated as the time difference Tb.
  • FIG. 6 is a diagram showing an example of the relationship between the three-phase voltage command and the carrier wave when the AC motor according to the first embodiment is rotating at a low speed.
  • FIG. 7 is a diagram showing an example of the relationship between the three-phase voltage command and the carrier wave when the AC motor according to the first embodiment is rotating at high speed.
  • the vertical axis represents the modulation factor and the horizontal axis represents the time.
  • the three-phase voltage commands v u * , v v * , and v w * shown in FIGS. 6 and 7 are not subjected to zero vector modulation by the zero vector modulation unit 34.
  • the frequencies of the three-phase voltage commands v u * , v v * , and v w * are sufficiently low with respect to the frequency of the carrier wave Sc.
  • the speed electromotive force of the AC motor 1 increases in proportion to the speed of the AC motor 1
  • the amplitude of * is relatively small. Therefore, the three-phase voltage commands v u * , v v * , and v w * have a gentle slope, and the intermediate phase voltage commands have a gentle slope.
  • the intermediate phase is the v phase from the first half to a part of the latter half, and the u phase is the rest of the latter half.
  • the three-phase voltage commands v u * , v v * , v w * are compared with the case where the AC motor 1 is rotating at a low speed.
  • the frequency is high.
  • PWM Pulse Width Modulation
  • the ratio of the frequencies of the three-phase voltage commands v u * , v v * , and v w * to the frequency of the carrier wave Sc is 10% to 15%. It is not uncommon to have a degree.
  • the AC motor 1 is rotating at a high speed, as shown in FIG.
  • the intermediate phases are v-phase, u-phase, w-phase, v-phase, and so on, as compared with the case where the AC motor 1 is rotating at a low speed.
  • ⁇ ⁇ It changes at high speed in the order of, and the speed electromotive force also increases. Therefore, the amplitudes of the three-phase voltage commands v u * , v v * , and v w * are also larger than those shown in FIG. 6, and the slopes of the intermediate phase voltage commands are also considerably larger. Therefore, the fluctuation of the time difference Tb between the current detection times is large.
  • the phase current determining unit 12 the bus current i the phase current i u of three-phase based on the value bus is detected every cycle of the carrier wave Sc by the bus current detection unit 11, i v, the process of determining the value of i w do.
  • the fluctuation of the time difference Tb between the current detection times is the largest with respect to the average time difference between the current detection times, once per cycle of the carrier wave Sc, and the three phase currents i u , iv , This is a case of determining the value of i w.
  • the time difference Tb between the current detection times changes in the range of 0.5 ⁇ Tbave ⁇ Tb ⁇ 1.5 ⁇ Tbave, where the average time difference between the current detection times is “Tbave”.
  • Tbave is the average time difference between the current detection times as described above. The longer the average time difference between the current detection times, the less the fluctuation of the time difference Tb between the current detection times, but the lower the time resolution for the waveforms of the phase currents i u , iv , and i w. Therefore, when the AC motor 1 is rotating at a high speed, it may be difficult to increase the estimation accuracy of the rotation speed ⁇ e.
  • FIG. 8 is a diagram showing an example of a change in the voltage command of the intermediate phase when the voltage command of the three phases according to the first embodiment is a sine wave.
  • FIG. 9 is a diagram showing an example of a change in the voltage command of the intermediate phase among the voltage commands of the three phases deformed by using the third harmonic superimposition according to the first embodiment.
  • the waveform of the intermediate phase voltage command is substantially triangular. Further, as shown in FIG.
  • the waveform of the voltage command of the intermediate phase is , Almost a triangular wave.
  • the intermediate phase waveform may not be a triangular wave as shown in FIGS. 8 and 9, but the degree of freedom in selecting the zero voltage vector is reduced near the modulation factor of 1. Therefore, the waveform of the intermediate phase is quite close to the triangular wave.
  • the frequency of these triangular waves is three times the three-phase voltage commands v u * , v v * , and v w *. Therefore, the waveform of the intermediate phase contains a frequency component of the 3nth order with respect to the fundamental frequency of the three-phase voltage commands v u * , v v * , and v w *. n is an odd number of 1 or more.
  • FIG. 10 is a diagram showing an example of a three-phase voltage command in which zero vector modulation is performed by the zero vector modulation unit according to the first embodiment.
  • the vertical axis represents the modulation factor and the horizontal axis represents the time.
  • random noise is generated with respect to the waveforms of the voltage commands v u * , v v * , v w * on which the third harmonic shown in FIG. 9 is superimposed. It is a waveform as if it was added.
  • the spectral peaks of the carrier noise are dispersed by the zero vector modulation, and the carrier noise can be reduced.
  • the waveform of the intermediate phase may change abruptly even when the AC motor 1 is rotating at a low speed.
  • the zero vector modulation unit 34 increases the random value to be added to the voltage commands v u * , v v * , and v w * as the rotation of the AC motor 1 is slower.
  • the waveform of the voltage command of the intermediate phase changes abruptly even when the AC motor 1 is rotating at a low speed, and the change in the time difference Tb between the current detection times becomes large. In this case, if the fluctuation of the time difference Tb between the current detection times is ignored, good results may not be obtained even if the speed estimation or the position estimation is performed.
  • 11 and 12 are diagrams showing an image of an integral calculation by rectangular approximation.
  • FIG. 13 is a diagram showing an image of an integral calculation when the fluctuation of the time difference is ignored.
  • the vertical axis represents the value of an arbitrary curve x (t) at time t, and the horizontal axis represents the time.
  • the integration result can be obtained by laying out a plurality of rectangles each having a width of a minute time T in the width direction and adding up the areas of each rectangle. In this case, each rectangle is laid out so as not to overlap with the other rectangles and to have no gap between adjacent rectangles.
  • FIG. 11 shows an example of integration by rectangular approximation when the minute time T is constant, and an error of integration calculation occurs, but in many cases, there is no problem in integration by rectangular approximation. If you want to perform more accurate integration calculation, you can handle it by performing integration by trapezoidal approximation.
  • the integration operation When the integration operation is performed without using ⁇ T (t), as shown in FIG. 13, it may overlap with other rectangles or a gap may occur between adjacent rectangles. Regarding the total area, the overlapped portion and the gap portion may cancel each other out to form a large error, but a large error occurs in the partial area. Therefore, if the integration operation is performed with the time difference Tb between the current detection times as a fixed value, an error occurs in the result of the integration operation.
  • FIG. 14 is a diagram for explaining an error in the differential operation.
  • Tb T + ⁇ T (t)
  • ⁇ T (t) is a fluctuation amount of the time difference Tb, and it is assumed that it changes at each calculation timing.
  • an error occurs in the result of the differential operation by the amount of ⁇ T (t).
  • the estimation unit 15 estimates the magnetic pole position ⁇ e and the rotation speed ⁇ e of the AC motor 1 based on the time difference Tb calculated by the time difference calculation unit 13.
  • the estimation unit 15 estimates the magnetic pole position ⁇ e and the rotation speed ⁇ e by using the arc tangent method.
  • the method of estimating the magnetic pole position ⁇ e and the rotation speed ⁇ e is not limited to the arc tangent method.
  • a method of estimating the magnetic pole position ⁇ e and the rotation speed ⁇ e there is a method using an adaptive observer. The adaptive observer will be described in the second embodiment.
  • the arctangent method is the most primitive position estimation method and is widely known.
  • the following equation (1) is a voltage equation of the surface magnet type synchronous AC motor on the stator coordinates.
  • "p” is a differential operator
  • "R a” is the armature resistance
  • "L a” is the armature inductance
  • "v alpha, v beta” the voltage on the stator coordinates
  • the stator coordinates may be described as static coordinates.
  • the second term on the right side of the above equation (1) represents the velocity electromotive force.
  • the term of velocity electromotive force can be expressed as the following equation (2).
  • equation (2) "e ⁇ " is the ⁇ -axis velocity electromotive force
  • e ⁇ " is the ⁇ -axis velocity electromotive force
  • p is the differential operator
  • ⁇ ⁇ r is the rotor ⁇ -axis magnetic flux
  • is the rotor ⁇ -axis magnetic flux
  • ⁇ a is the number of armature interlinkage magnetic fluxes.
  • the velocity electromotive force includes the information of the magnetic pole position ⁇ e which is the magnetic pole position information
  • the above equation (1) is arranged and the magnetic pole position ⁇ e is calculated.
  • the following equation (3) is obtained by summarizing the rotor magnetic flux term on the left side and the other terms on the right side. Since the differential calculation amplifies noise, both sides of the following equation (3) are integrated to obtain the following equation (4).
  • the voltage sensor has a DC offset
  • the integrated value will diverge if pure integration is used. Therefore, when calculating the following equation (4), use an approximate integral using a low-pass filter and do not integrate the DC component. It is customary to do.
  • the symbol " ⁇ " represents an estimated value.
  • the rotor magnetic flux is calculated by calculating the above equation (4), and the arc tangent shown in the following equation (5) is calculated using the obtained rotor magnetic flux to estimate the magnetic pole position ⁇ e of the rotor. be able to. Since the rotation speed ⁇ e can be calculated by using the estimated magnetic pole position ⁇ e of the rotor, the estimated speed ⁇ e ⁇ can be calculated by the following equation (6). However, in order to avoid the influence of differential noise, it is common to apply a low-pass filter when this estimated speed ⁇ e ⁇ is used for control.
  • the estimated speed ⁇ e ⁇ can also be calculated by estimating the speed electromotive force and dividing the amplitude of the speed electromotive force by the armature interlinkage magnetic flux number ⁇ a as shown in the following equation (7). ..
  • the velocity electromotive force can be calculated by calculating the right side of the above equation (3).
  • the error by the speed estimation method by the following equation (6) is smaller than the error by the method shown by the following equation (7).
  • the arc tangent method has the advantage that the amount of calculation is small. In the first embodiment, the case where the speed is estimated by the following equation (6) will be described.
  • the following equation (8) can be said to be a discretized integral of the first term on the right side of the above equation (4).
  • the low-pass filter is discretized by the backward difference, but the low-pass filter may be discretized by using the forward difference or the bilinear transform.
  • the rotor magnetic flux of the stationary coordinate system can be discretized as in the following equation (9). Therefore, the estimated position ⁇ e ⁇ (n) at the sampling point this time can be calculated by the following equation (10). Further, using the estimated position ⁇ e ⁇ (n) at the estimated position ⁇ e ⁇ (n) and the preceding sampling point in the present sampling point, the estimated velocity omega e ⁇ is between two sampling points, the following equation (11 ) Can be calculated. “ ⁇ e ⁇ ” is an estimated value of “ ⁇ e ”, and “ ⁇ e ⁇ ” is an estimated value of “ ⁇ e”.
  • the time difference Tb between the current detection times fluctuates from moment to moment. Therefore, if the calculation is not performed using the time difference Tb between the current detection times, an error will occur in the estimated position ⁇ e ⁇ and the estimated speed ⁇ e ⁇ .
  • the estimation unit 15 performs an calculation by the arctangent method using the time difference Tb between the current detection times.
  • Tb between the current detection time expressed as T S + [Delta] T (n)
  • operation of the low-pass filter represented by the above formula (8) can be rewritten as the following equation (12).
  • the estimation unit 15 coordinates-converts the phase currents i u , iv , i w determined by the phase current determination unit 12 for each sampling point into the currents i ⁇ , i ⁇ of the static coordinate system. Further, the estimation unit 15 coordinates the voltage commands v u * , v v * , v w * into the voltages v ⁇ and v ⁇ of the static coordinate system. The estimation unit 15 performs the calculation of the above equation (12) based on the currents i ⁇ , i ⁇ and the voltages v ⁇ , v ⁇ of the stationary coordinate system and the time difference Tb calculated by the time difference calculation unit 13. ..
  • the estimation unit 15 substitutes the value calculated by the above equation (12) and the currents i ⁇ and i ⁇ of the rest coordinate system into the above equation (9) to calculate the rotor magnetic flux of the rest coordinate system.
  • the estimation unit 15 substitutes the rotor magnetic flux of the stationary coordinate system calculated by the above equation (9) into the above equation (10) to calculate the estimated position ⁇ e ⁇ .
  • the estimation unit 15 substitutes the estimated position ⁇ e ⁇ calculated by the calculation of the above equation (10) and the time difference Tb calculated by the time difference calculation unit 13 into the above equation (13) to estimate the speed. Calculate ⁇ e ⁇ .
  • the carrier frequency which is the frequency of the carrier wave Sc
  • the phase current determination unit 12 determines the values of the three phase currents i u , iv , and i w once in one cycle of the carrier wave Sc. This will be described by taking as an example.
  • Control period T S is the inverse of the carrier frequency is 200 [mu] s.
  • the magnetic pole position ⁇ e of the AC motor 1 advances 0.5 rad during 200 ⁇ s.
  • the time difference Tb between the current detection times varies in the range of 100 ⁇ s to 300 ⁇ s. Therefore, for example, when the time difference Tb is 300 ⁇ s, the magnetic pole position ⁇ e advances 0.75 rad during 300 ⁇ s.
  • the estimated velocity omega e ⁇ will become 3750Rad, it becomes 1.5 times the true value of the rotational speed omega e It ends up. This estimation error is so large that it cannot be ignored in controlling the AC motor 1.
  • the estimation unit 15 Since the estimation unit 15 performs the calculation of the above equation (13), 0.75 rad is divided by the time difference Tb. Therefore, the estimated speed ⁇ e ⁇ is 2500 rad, which matches the true value of the rotation speed ⁇ e. In this way, by performing the calculation using the above equation (13), the estimation unit 15 can significantly reduce the error of the speed estimation as compared with the case where the calculation using the above equation (11) is performed. .. Further, by performing the calculation of the above formula (12), the estimation unit 15 can significantly reduce the error of the position estimation as compared with the case of performing the calculation of the above formula (8).
  • the estimation device 9 includes a bus current detection unit 11, a phase current determination unit 12, a time difference calculation unit 13, and an estimation unit 15.
  • the bus current detection unit 11 is the DC bus of the voltage inverter 4 driven by a plurality of gate pulse signals Gu , G v , G w generated based on the voltage commands v u * , v v * , v w *.
  • the value of the bus current ibus which is the current flowing through 61, is detected.
  • Phase current determination unit 12 the bus current i bus values and a plurality of gate pulse signal G u detected by the bus current detection unit 11, G v, based on the state of the G w, the AC motor from the voltage-source inverter 4
  • the values of the phase currents i u , iv , and i w supplied to 1 are determined.
  • Time difference calculating section 13 was used in this determination processing by the current detection time and the phase current determining unit 12, which is identified as the detection time of the bus current i bus used in the previous determination process by the phase current determining section 12 calculating a time difference Tb between the current detection time specified as the detection time of the bus current i bus.
  • the estimation unit 15 determines the magnetic pole position ⁇ e of the AC motor 1 based on the voltage commands v u * , v v * , v w * , the values of the phase currents i u , i v , i w, and the time difference Tb. And at least one of the rotational speeds ⁇ e is estimated.
  • the estimation device 9 performs the estimation process using the time difference Tb between the current detection times of the phase current determination unit 12. Therefore, the estimation device 9 is compared with the case of determining the value of the bus current i value three-phase current by interpolation from the bus detected by each of the two half cycles of continuous in the carrier wave Sc, 3 single phase current
  • the determination cycle which is the cycle of processing for determining the values of i u , iv , and i w, can be shortened. Further, since the estimation device 9 uses the time difference Tb between the current detection times to estimate the magnetic pole position ⁇ e and the rotation speed ⁇ e , even when the time difference Tb between the current detection times fluctuates greatly, the estimation device 9 uses the time difference Tb.
  • the estimation accuracy of the magnetic pole position ⁇ e and the rotation speed ⁇ e can be improved. Therefore, when the AC motor 1 is rotating at high speed or when zero vector modulation is performed, the estimation accuracy of the magnetic pole position ⁇ e and the rotation speed ⁇ e can be improved.
  • the gate pulse signals G u , G v , G w are generated based on the comparison between the voltage commands v u * , v v * , v w * and the carrier wave Sc.
  • the time difference Tb between the current detection times is longer than the half-cycle period of the carrier wave Sc and shorter than 1.5 times the period of the carrier wave Sc.
  • the drive device 2 including the estimation device 9 includes a pulse shift processing unit 22 that temporally shifts at least one of a plurality of gate pulse signals Gu , G v , and G w input to the voltage inverter 4. ..
  • the pulse shift processing unit 22 shifts at least one of the two gate pulse signals Gu , G v , and G w of the gate pulse signals G u, G v, and G w whose change interval is less than the preset time TA by time.
  • the interval between changes of the two gate pulse signals is set to be equal to or longer than the preset time TA. As a result, it is possible to secure a waiting time until the ringing converges regardless of the values of the voltage commands v u * , v v * , and v w *.
  • the drive device 2 irregularly changes the output ratio of the two types of zero voltage vectors output from the voltage type inverter 4 by adding a random value to the voltage commands v u * , v v * , and v w *.
  • a zero vector modulation unit 34 is provided.
  • the spectral peaks of the carrier noise are dispersed by the zero vector modulation, and the carrier noise can be reduced.
  • the waveform of the voltage command of the intermediate phase changes abruptly even when the AC motor 1 is rotating at a low speed, and the change in the time difference Tb between the current detection times becomes large. Even in such a case, the estimation device 9 can improve the estimation accuracy of the magnetic pole position ⁇ e and the rotation speed ⁇ e.
  • Embodiment 2 The estimation unit of the estimation device according to the second embodiment is different from the estimation unit of the estimation device according to the first embodiment in that the position and speed of the AC motor are estimated by using the adaptive observer.
  • the components having the same functions as those in the first embodiment will be designated by the same reference numerals and the description thereof will be omitted, and the differences from the first embodiment will be mainly described.
  • FIG. 15 is a diagram showing an example of a configuration example of an electric motor system including the estimation device according to the second embodiment of the present invention.
  • the electric motor system 100A according to the second embodiment includes a drive device 2A having a control device 5A instead of the drive device 2 having the control device 5.
  • the control device 5A includes an estimation device 9A having an estimation unit 15A instead of the estimation device 9 having the estimation unit 15.
  • FIG. 16 is a diagram showing a configuration example of the estimation unit according to the second embodiment.
  • the estimation unit 15A is composed of an adaptive observer, and estimates the magnetic pole position ⁇ e and the rotation speed ⁇ e of the AC motor 1.
  • the adaptive observer is also called an adaptive magnetic flux observer because it is defined by an equation of state in which the stator magnetic flux and the rotor magnetic flux of the AC motor 1 are state quantities.
  • An adaptive observer can also be configured by adopting an extended induced voltage or current as a state quantity.
  • the estimation unit 15A shown in FIG. 16 calculates the estimated speed ⁇ e ⁇ and the estimated position ⁇ e ⁇ by using the voltage vector, the current vector, and the time difference Tb between the current detection times, and the calculated estimation.
  • the velocity ⁇ e ⁇ and the estimated position ⁇ e ⁇ are output.
  • the voltage vector is a vector of the voltage output from the voltage type inverter 4 to the AC motor 1
  • the current vector is a vector of the current output from the voltage type inverter 4 to the AC motor 1.
  • the estimation unit 15A includes a model deviation calculation unit 51, an angular velocity estimator 52, a primary angular frequency calculation unit 53, and an integrator 54.
  • the model deviation calculation unit 51 calculates the model deviation ⁇ based on the voltage vector, the current vector, the primary angular frequency ⁇ 1 , the estimated speed ⁇ e ⁇ , and the time difference Tb between the current detection times.
  • the angular velocity estimator 52 calculates the estimated velocity ⁇ e ⁇ based on the model deviation ⁇ and the time difference Tb between the current detection times.
  • the primary angular frequency calculator 53 calculates the primary angular frequency ⁇ 1 using the estimated magnetic flux vector, the estimated current vector, and the estimated velocity ⁇ e ⁇ .
  • the integrator 54 integrates the primary angular frequency ⁇ 1 based on the time difference Tb between the current detection times and outputs the estimated position ⁇ e ⁇ .
  • the model deviation calculation unit 51 includes a current estimator 71, a subtractor 72, and a deviation calculation device 73.
  • the current estimator 71 calculates the estimated magnetic flux vector and the estimated current vector based on the voltage vector, the current vector, the primary angular frequency ⁇ 1 , the estimated speed ⁇ e ⁇ , and the time difference Tb between the current detection times. Then, the calculated estimated magnetic flux vector and the estimated current vector are output.
  • the subtractor 72 calculates the current deviation vector by subtracting the current vector from the estimated current vector, and outputs the calculated current deviation vector.
  • the deviation calculator 73 takes the current deviation vector from the subtractor 72 as an input, extracts the orthogonal component of the estimated magnetic flux vector as a scalar amount, and outputs the extracted scalar amount as a model deviation ⁇ .
  • a method of extracting the orthogonal component of the estimated magnetic flux vector as a scalar quantity a method of converting the coordinates of the current deviation vector onto a rotating Cartesian coordinate system and a method of calculating the magnitude of the outer product value of the current deviation vector and the estimated magnetic flux vector are calculated. The method is known.
  • the current estimator 71 estimates the current and the magnetic flux from the equation of state of the AC motor 1.
  • the AC motor 1 is a general embedded permanent magnet synchronous motor, but if the equation can be formulated, it is an induction motor, a surface type permanent magnet synchronous motor, a wound field magnetic synchronous motor, or a relaxation. It may be another kind of electric motor such as a type synchronous motor. That is, the current estimator 71 can estimate the current for other types of electric motors in the same manner.
  • Ra is the armature resistance
  • ⁇ 1 is the primary angular frequency
  • v d is the d-axis voltage
  • v q is the q-axis voltage
  • h 11 is the q-axis voltage
  • h 11 is the q-axis voltage
  • h 11 is the q-axis voltage
  • h 11 to “h 32 ” are observer gains.
  • the primary angular frequency ⁇ 1 is given by the following equation (16).
  • H 41 ” and “h 42 ” represent observer gains.
  • the above equations (14) and (15) are equations based on the usual induced voltage, but the same calculation can be performed even if the above equations (14) and (15) are modified and expressed in the form of extended induced voltage. Can be done.
  • the above equations (14) and (15) are mathematical expressions in dq rotating coordinates, but in other coordinate systems such as stationary ⁇ coordinates or three-phase coordinates by modifying the above equations (14) and (15). The same calculation can be done even if it is expressed. Since the above equation (14) includes the estimated speed ⁇ e ⁇ , if the estimated speed ⁇ e ⁇ and the actual rotation speed ⁇ e do not match, an error occurs in the current estimation.
  • the model deviation ⁇ is defined by the following equation (17), and the estimation unit 15A adjusts the value of the estimated velocity ⁇ e ⁇ by using the angular velocity estimator 52 so that the model deviation ⁇ becomes zero.
  • the angular velocity estimator 52 is configured by, for example, connecting an integrator in series to a proportional integration controller.
  • the primary angular frequency calculator 53 calculates the primary angular frequency ⁇ 1 from the estimated magnetic flux vector, the estimated current vector, and the estimated speed ⁇ e ⁇ based on the above equation (16).
  • the integrator 54 estimates the magnetic pole position ⁇ e by integrating the primary angular frequency ⁇ 1.
  • the adaptive observer is a high-performance velocity estimation method because it is robust to fluctuations in the number of interlinkage magnetic fluxes and is excellent in that a steady velocity estimation error does not occur.
  • equation (21) can be obtained from the following equation (18).
  • T S is the control cycle
  • n is a value representing the current sampling point
  • n-1 is a value representing the previous sampling point
  • t avg the carrier wave Sc from the current detection sampling point Represents the time to the valley of.
  • the valley of the carrier wave Sc indicates, for example, the position where the carrier wave Sc becomes the minimum value, and appears at the times t7 and t16 shown in FIG. 5, respectively.
  • the current detection sample point is, for example, the above-mentioned current detection time.
  • the primary angular frequency ⁇ 1 is given by the following equation (22).
  • the model deviation calculation unit 51 of the estimation unit 15A performs coordinate conversion using the estimated position ⁇ e ⁇ . Using the estimated position ⁇ e ⁇ obtained by the integrator 54, the model deviation calculation unit 51 coordinates the three phase currents i u , iv , i w to the currents in the dq coordinate system, which is a rotating Cartesian coordinate system. By converting, the d-axis current id and the q-axis current i q are obtained. Further, the model deviation calculation unit 51 uses the estimated position ⁇ e ⁇ obtained by the integrator 54 to rotate the three-phase voltage commands v u * , v v * , v w * , which is a Cartesian coordinate system dq.
  • the d-axis voltage v d and the q-axis voltage v q are obtained by converting the coordinates to the voltage of the coordinate system.
  • the estimated position ⁇ e ⁇ used for coordinate conversion to obtain the d-axis voltage v d and the q-axis voltage v q may be, for example, the one at the center time between the current detection times.
  • the equations of state of the above equations (19) and (20) use the stator magnetic flux and the rotor magnetic flux as the state quantities, the current or the extended induced voltage may be used as the state quantities.
  • the model deviation ⁇ is defined as the following equation (23), and the estimation unit 15A adjusts the value of the estimated velocity ⁇ e ⁇ by using the angular velocity estimator 52 so that the model deviation ⁇ becomes zero. ..
  • the angular velocity estimator 52 calculates the estimated velocity ⁇ e ⁇ by the calculation of the following equation (24) to equation (26).
  • "x a " is the state quantity of proportional integration control
  • "n + 1" is the next sampling point.
  • the estimation unit 15A can accurately calculate the estimation speed ⁇ e ⁇ even when the time difference Tb between the current detection times fluctuates.
  • the primary angular frequency calculator 53 calculates the primary angular frequency ⁇ 1 by substituting the estimated magnetic flux vector, the estimated current vector, and the estimated velocity ⁇ e ⁇ into the above equation (22).
  • the estimation unit 15A includes a current estimator 71, a subtractor 72, a deviation calculator 73, and an angular velocity estimator 52.
  • the current estimator 71 is an example of a state quantity estimation unit
  • the subtractor 72 and the deviation calculator 73 are an example of an error calculation unit
  • the angular velocity estimator 52 is an example of a velocity estimation unit.
  • the current estimator 71 estimates the magnetic flux vector and the current vector of the AC motor 1, and outputs the estimated magnetic flux vector and the estimated current vector which are the estimation results.
  • the magnetic flux vector is an example of a state quantity.
  • the deviation calculator 73 calculates the model deviation ⁇ .
  • the model deviation ⁇ is an example of the error of the state quantity.
  • the angular velocity estimator 52 calculates the estimated velocity ⁇ e ⁇ based on the model deviation ⁇ and the time difference Tb. Further, the current estimator 71 uses an estimated magnetic flux vector and an estimated current vector based on the voltage vector, the current vector, the primary angular frequency ⁇ 1 , the estimated speed ⁇ e ⁇ , and the time difference Tb between the current detection times. Is calculated. Since the estimation unit 15A does not generate a steady speed estimation error, the speed estimation can be performed with high accuracy.
  • Embodiment 3 is different from the drive device according to the first and second embodiments in that the position and speed of the AC motor are estimated by using the voltage applied to the AC motor from the voltage inverter during the current detection time.
  • components having the same functions as those of the first embodiment will be designated by the same reference numerals and description thereof will be omitted, and the differences from the control devices 5 and 5A of the first and second embodiments will be mainly described.
  • FIG. 17 is a diagram showing an example of a configuration example of an electric motor system including the estimation device according to the third embodiment of the present invention.
  • the electric motor system 100B according to the third embodiment includes a drive device 2B having a control device 5B in place of the drive device 2 having the control device 5 or the drive device 2A having the control device 5A.
  • the control device 5B includes an estimation device 9B in place of the estimation device 9 or the estimation device 9A.
  • the estimation device 9B includes a bus current detection unit 11, a phase current determination unit 12, a time difference calculation unit 13, a detection time interval voltage calculation unit 14, and an estimation unit 15B.
  • Detection time voltage calculation unit 14 the gate pulse signal G u, G v, and G w, based on the time difference Tb between the current detection time, current detection time between the three phase voltages v ub, v vb, Calculate the value of v wb.
  • the values of the phase voltages v ub , v v b , and v wb calculated by the detection time interval voltage calculation unit 14 are the instantaneous values of the three phase voltages which are the u-phase, v-phase, and w-phase voltages during the current detection time.
  • the phase voltage v ub is the u-phase voltage
  • the phase voltage v vb is the v-phase voltage
  • the phase voltage v w b is the w-phase voltage.
  • phase voltage v ub described as u-phase voltage v ub the phase voltage v vb described as v-phase voltage v vb
  • phase voltage v wb to as w-phase voltage v wb
  • phase voltage v ub, v vb, v simply phase voltage values of wb v ub, v vb, may be referred to as v wb.
  • Detection time voltage calculation unit 14 the current respective gate pulse signal G u between detection time, G v, based on the duty ratio of G w, 3 single phase voltage v ub between current detection time, v vb, v wb Calculate the value of. Specifically, the detection time interval voltage calculation unit 14 determines the u-phase voltage v u during the current detection time based on the ratio of the on time when the gate pulse signal Gu is on during the current detection time. Calculate the value. Similarly, the detection time interval voltage calculation unit 14 is based on the on-time, which is the time during which the gate pulse signal Gv is on during the current detection time, and the value of the v-phase voltage v v during the current detection time. Is calculated. Further, the detection time interval voltage calculation unit 14 sets the value of the w-phase voltage v w during the current detection time based on the on time, which is the time during which the gate pulse signal G w is on during the current detection time. calculate.
  • FIG. 18 is a diagram showing a configuration example of the detection time interval voltage calculation unit according to the third embodiment.
  • the detection time interval voltage calculation unit 14 includes an on-time calculation unit 41 and a voltage calculation unit 42.
  • the on-time calculation unit 41 calculates the on-time, which is the time during which each of the plurality of gate pulse signals Gu , G v , and G w is on during the current detection time.
  • the voltage calculation unit 42 has three phase voltages v ub based on the value of the bus voltage V dc , the on time calculated by the on time calculation unit 41, and the time difference Tb calculated by the time difference calculation unit 13. , V vb , v wb values are calculated.
  • Bus voltage V dc is the voltage bus voltage is between the DC bus 61 between the DC bus 61 on the low potential side high potential side
  • the value of the bus voltage V dc is the instantaneous value of the bus voltage V dc ..
  • FIG. 19 is a diagram for explaining a method of calculating three phase voltages during the current detection time by the detection time interval voltage calculation unit according to the third embodiment.
  • time t5 is the detection time of the bus current i bus used to determine the value of the second phase of the currents in the first phase current determination cycle is identified as the first current detection time .
  • the gate pulse signals Gu , G v , and G w are not subjected to pulse shift processing by the pulse shift processing unit 22.
  • the detection time interval voltage calculation unit 14 can calculate the value of the u-phase voltage v ub during the current detection time by the calculation of the following equation (28).
  • “tuon1” is the on-time of the gate pulse signal G u u-phase in the first half of the current between the detection time
  • “tuon2” gate pulse signal G of the u-phase in the second half of the current between the detection time It is the on-time of u.
  • V dc is a bus voltage, which is detected by a bus voltage detecting means (not shown). When the bus voltage is constant, it is not necessary to use the bus voltage detecting means.
  • the detection time interval voltage calculation unit 14 calculates the value of the v-phase voltage v vb during the current detection time and the w-phase voltage v wb between the current detection times by the calculation of the following equations (29) and (30). The value can be calculated.
  • “Tvon1” is the time that is gate pulse signal G v of the first half in the v-phase between the current detection time is turned on
  • “tvon2” is a gate pulse signal G v v-phase in the second half of the current between the detection time on It is the time that has become.
  • the time difference Tb between the current detection times does not coincide with the period of the carrier wave Sc in most cases. Therefore, the value of the u-phase voltage v ub during the current detection time does not match the voltage command v u * in most cases. For the same reason, the value of the v-phase voltage v vb also, in most cases, does not match the voltage command v v *, the value of the w-phase voltage v wb also, in most cases, does not match the voltage command v w *. Therefore, when the voltage commands v u * , v v * , and v w * are used for speed estimation and position estimation, an error occurs in speed estimation and position estimation due to the voltage error.
  • the time difference Tb between the current detection times is set to the carrier cycle Tc.
  • the carrier cycle Tc varies in the range of 0.5 ⁇ Tc ⁇ Tb ⁇ 1.5 ⁇ Tc. In this case, the fluctuation of the voltage error is large, and the fluctuation of the error that occurs in the speed estimation and the position estimation is also large.
  • Estimating unit 15B the speed estimation and the position estimated by using the three phase voltages v ub between current detection time calculated by the detection time between the voltage calculating portion 14, v vb, v the value of wb. Therefore, the estimation unit 15B can perform speed estimation and position estimation with high estimation accuracy.
  • Estimating unit 15B three phase voltages v ub between current detection time calculated by the detection time between the voltage calculating portion 14, v vb, v voltage of stationary coordinate system wb v alpha, v coordinate transformation into beta. That is, the estimation unit 15B * voltage command v u, v v *, v w * in place of the current detection time between the three phase voltages v ub, v vb, stationary coordinate system obtained the v wb by coordinate transformation It differs from the estimation units 15 and 15A in that the voltages v ⁇ and v ⁇ of are used.
  • the estimation unit 15B calculates the difference between the estimated value and the true value of the state quantity and measures the error.
  • the subtractor 72 calculates the difference between the estimated current vector and the current vector. In the process of calculating the difference between the estimated current vector and the current vector, if there is a time lag between the estimated current vector and the current vector, the calculation result of the difference between the estimated current vector and the current vector may not be accurate. ..
  • the current estimator 71 of the estimation unit 15B accurately obtains the estimated current vector of the current detection time by using the values of the three phase voltages vub , vvb , and vwb during the current detection time. be able to.
  • the three phase voltages v ub between current detection time by the detection time between the voltage calculating portion 14, v vb specifically described an example of a method of calculating the value of v wb.
  • a method for detecting the on-time of each gate pulse signal Gu , G v , and G w during the current detection time will be described.
  • the on-time calculation unit 41 of the detection time interval voltage calculation unit 14 divides the period between the current detection times into three or more sections.
  • the on-time calculation unit 41 sets the on-time, which is the time during which the gate pulse signals Gu , G v , and G w are on in each of the three or more divided sections, as the gate pulse signals G u , G v , and G. Obtained for each w.
  • the on-time calculation unit 41 adds up the on-time of each of the three or more divided sections for each gate pulse signal G u , G v , G w , and sets the on-time between the current detection times to the gate pulse signal G. Calculated for each u , G v , and G w. As a result, the on-time calculation unit 41 can accurately calculate the on-time of each gate pulse signal Gu , G v , and G w during the current detection time.
  • On-time computing unit 41 of the detection time between the voltage calculating portion 14 calculates each of the on-time tuon1 and on-time Tuon2, by summing the ON time tuon1 and on time Tuon2, on-time of the gate pulse signal G u Ask for.
  • the relative time from the time of the apex of the carrier wave Sc to the current detection time is defined as Tdet1 and Tdet2.
  • the apex of the carrier wave Sc is the valley of the carrier wave Sc.
  • the time of the valley of the carrier wave Sc is time t7 and t16.
  • the relative time Tdet1 is the time until the time t5, which is the current detection time, and the time t7, when the carrier wave Sc reaches the apex.
  • the relative time Tdet2 is the time between the time t15, which is the current detection time, and the time t16, when the carrier wave Sc reaches the apex.
  • the phase current determination unit 12 can calculate the relative times Tdet1 and Tdet2 indicating the relative time difference between the timing when the carrier wave Sc reaches the apex and the current detection time.
  • the phase current determination unit 12 notifies the detection time interval voltage calculation unit 14 of the calculated relative times Tdet1 and Tdet2.
  • the carrier cycle from time t1 to time t11 is defined as the first carrier cycle
  • the carrier cycle from time t11 to time t17 is defined as the second carrier cycle.
  • the on-time of the downlink half cycle is the same as the relative time Tdate1.
  • the on-time calculation unit 41 determines that the on-time of the downlink half cycle of the on-time tuon1 is the same as the relative time Tdate1.
  • the ON time of the downlink half-cycle of the on-time tuon1 can be obtained from the gate pulse signal G u u-phase.
  • the on-time calculation unit 41 the time from the time that the gate pulse signal G u u-phase is switched from off to on until the time t7 of the vertices of the carrier wave Sc, as the on-time of the downlink half-cycle of the on-time tuon1 Can be calculated.
  • On-time computing unit 41 the on-time of the uplink semi-period of the ON time Tuon1, calculated from the gate pulse signal G u u-phase.
  • On-time calculation unit 41 calculates the time from time t7 vertex of the carrier wave Sc to the time that the gate pulse signal G u u-phase is switched from on to off, as the on time of the uplink semi-period of the ON time tuon1 be able to.
  • On-time computing unit 41 calculates the ON time of the gate pulse signal G u u-phase in the downlink half-life of the second carrier period, by subtracting the relative time Tdet2 from the calculated time to obtain the on-time tuon2 be able to. For example, the on-time calculation unit 41, the second carrier period, and calculates the time from the time t12 to the gate pulse signal G u u-phase is switched from off to on until the time t16 vertex of the carrier wave Sc. The on-time calculation unit 41 obtains the on-time tuon2 by subtracting the relative time Tdate2 from the time from the time t12 to the time t16. However, since the minimum value of the on-time tuon2 is zero, the on-time calculation unit 41 forcibly sets the on-time tuon2 to zero when the calculation result becomes negative.
  • the phase current determination unit 12 can also calculate the time difference Tb between the current detection times based on the relative times Tdet1 and Tdet2.
  • the variation ⁇ T (t) of the time difference Tb of the current detection time is the difference between the relative time Tdet1 and the relative time Tdet2.
  • the time difference Tb of the current detection time is obtained by adding the variation ⁇ T (t) to the carrier period Tc.
  • the phase current determination unit 12 can calculate the time difference Tb between the current detection times by the calculation of the following equation (31).
  • control device 5B does not detect the absolute time of the current detection time, but based on the relative times Tdate1 and Tdate2, the time difference Tb between the current detection times and the current detection time.
  • the values of the three phase voltages v ub , v vb , and v wb can be calculated.
  • the on-time calculation unit 41 replaces the gate pulse signals G u , G v , and G w with the gate pulse signals G u , based on the three-phase voltage commands v u * , v v * , v w *. It is also possible to calculate the on-time of G v and G w. In this case, the on-time calculation unit 41 calculates the time when the gate pulse signals Gu , G v , and G w are switched from off to on from the magnitudes of the three-phase voltage commands v u * , v v * , and v w *. To do. Further, the time from when the gate pulse signal G of the intermediate phase is switched from off to on until the time TA elapses is calculated.
  • the on-time calculation unit 41 determines the time when the gate pulse signals G u , G v , and G w are switched from off to on, and the time TA elapses after the intermediate phase gate pulse signal G is switched from off to on.
  • the on-time of the gate pulse signals Gu , G v , and G w can be calculated based on the time of.
  • the on-time calculation unit 41 of the detection time interval voltage calculation unit 14 divides the period between the current detection times into three sections, and obtains the on time of the gate pulse signal G in each of the three divided sections. ..
  • the on-time calculation unit 41 calculates the on-time of the gate pulse signal G between the current detection times by adding up the on-time of each of the three divided sections.
  • the three sections are a downlink half cycle in which current detection is performed first, an uplink half cycle in which current detection is not performed, and a downlink half cycle in which current detection is performed next.
  • the on-time calculation unit 41 of the detection time-to-detection voltage calculation unit 14 divides the on-time of the gate pulse signal G between the current detection times by the time difference Tb between the current detection times to obtain a gate pulse signal between the current detection times. Calculate the duty ratio of G.
  • the voltage calculation unit 42 of the detection time interval voltage calculation unit 14 calculates the voltage during the current detection time from the duty ratio of the gate pulse signal G during the current detection time.
  • the detection time interval voltage calculation unit 14 sets the current detection time interval to four or more sections even when the period in which the phase current is detected by the phase current determination unit 12 is 1.5 times or more the period of the carrier wave Sc.
  • the voltage between the current detection times can be calculated by dividing and calculating by case.
  • the detection time interval voltage calculation unit 14 calculates the period during which the gate pulse signals Gu , G v , and G w are on, but the gate pulse signals G u , G v , and G w are calculated. You can also calculate how long is off. The length of the down half cycle and the up half cycle of the carrier wave Sc is known. Therefore, by calculating the period during which the gate pulse signals G u , G v , G w are off during the current detection time, indirectly, the gate pulse signals G u , G v , G during the current detection time. You can find the time when w is on.
  • the time difference Tb between the current detection times does not coincide with the carrier period Tc in most cases, and the higher the rotation speed of the AC motor 1, the larger the fluctuation.
  • the waveform of the intermediate phase includes a frequency component of the 3nth order with respect to the fundamental frequency of the three-phase voltage commands v u * , v v * , v w *.
  • n is an integer of 1 or more. Therefore, when the AC motor 1 is rotating at high speed, the deviation between the time difference Tb between the current detection times and the carrier cycle Tc becomes remarkable in the cycle of the electric angle 3f.
  • the electric angle 3f is a frequency three times the electric angular frequency of the AC motor 1.
  • FIG. 20 is a diagram showing an example of the result of speed estimation when the time difference between the current detection times according to the third embodiment is set as a fixed value.
  • FIG. 21 is a diagram showing the result of FFT analysis of the result of speed estimation shown in FIG. 20.
  • the frequency of the pulsating load is equal to the rotation frequency of the rotation shaft of the AC motor 1, and the number of pole pairs of the AC motor 1 is 3.
  • the vertical axis represents the rotation speed ⁇ e of the AC motor 1
  • the horizontal axis represents the time.
  • the vertical axis represents the rotation speed ⁇ e of the AC motor 1
  • the horizontal axis represents the frequency.
  • the waveform of the estimated speed has a large high-frequency vibration component as compared with the true speed value.
  • the main component of the vibration is a frequency component three times the electrical frequency fe of the AC motor 1.
  • the current detection time changes at a frequency three times the electrical frequency fe , but since the fluctuation of the time difference Tb between the current detection times is ignored, FIG. 20 shows. As shown, there is a large error in the estimated speed. Further, as shown in FIG.
  • the mechanical frequency is f is vibrating at 1 times the frequency of the m, the waveform of the estimated speed machine so as to cover the vibration a peak 1.4 times the frequency of the specific frequency f m. Therefore, 1 times the pulsation of mechanical frequency f m has become almost invisible.
  • the number of pole pairs of the AC motor 1 is 3, but the peaks of the mechanical 7.6f component and the mechanical 10.4f component stand as sidebands of the electric 3f component and the mechanical 9f component.
  • Electrical 3f component is three times the frequency component of the electrical frequency f e
  • mechanical 9f component is 9 times the frequency components of the mechanical frequency f m.
  • the machine 7.6f component is 7.6 times the frequency components of the mechanical frequency f m
  • the machine 10.4f component is 10.4 times the frequency components of the mechanical frequency f m.
  • the true speed value is due to the change in the current detection time by the one-shunt current detection method and the voltage error accompanying the change in the current detection time. There is an error in estimating the speed of the components not included in.
  • FIG. 22 is a diagram showing an example of the result of speed estimation by the estimation unit according to the third embodiment.
  • FIG. 23 is a diagram showing the result of FFT analysis of the result of speed estimation shown in FIG. 22.
  • the frequency of the pulsating load is equal to the mechanical frequency f m
  • the pole pair number of the AC motor 1 is 3.
  • the vertical axis is the rotation speed ⁇ e of the AC motor 1
  • the horizontal axis is the time.
  • the vertical axis represents the rotation speed ⁇ e of the AC motor 1
  • the horizontal axis represents the frequency.
  • the estimation unit 15B according to the third embodiment has a time difference Tb between the current detection times specified by the phase current determination unit 12 and a phase voltage vub between the current detection times calculated by the detection time voltage calculation unit 14. , V vb , v wb values are used to estimate the speed. As a result, the estimation unit 15B can significantly reduce the error in speed estimation, as shown in FIGS. 22 and 23.
  • the estimation unit 15B according to the third embodiment has a phase between the time difference Tb between the current detection times detected by the phase current determination unit 12 and the current detection time specified by the detection time voltage calculation unit 14. Position estimation is performed based on the values of the voltages v ub , v vb , and v wb. As a result, the estimation unit 15B can significantly reduce the error in position estimation.
  • the estimating apparatus 9B the current detection time between in the voltage source inverter 4 phase voltages applied to the AC motor 1 from v ub, v vb, v value voltage command of wb v u
  • a detection time interval voltage calculation unit 14 calculated based on * , v v * , and v w * is provided.
  • the detection time interval voltage calculation unit 14 is an example of a phase voltage calculation unit.
  • Estimating unit 15B based on the phase currents i u, i v, and the value of i w, and the time difference Tb, the detection time between the voltage phase voltage v ub calculated by the computation unit 14, v vb, the value of v wb Then, the estimated position ⁇ e ⁇ and the estimated speed ⁇ e ⁇ are calculated.
  • the detection time interval voltage calculation unit 14 includes an on-time calculation unit 41 and a voltage calculation unit 42.
  • the on-time calculation unit 41 calculates the on-time, which is the time during which each of the plurality of gate pulse signals Gu , G v , and G w is on during the current detection time.
  • the voltage calculation unit 42 has a phase voltage v ub , v based on the value of the bus voltage V dc , the on time calculated by the on time calculation unit 41, and the time difference Tb calculated by the time difference calculation unit 13. Calculate the values of vb and v wb. Accordingly, the detection time between the voltage calculation unit 14 can calculate the gate pulse signal G u, G v, the phase voltage v ub between the current detecting time from the G w, v vb, v the value of wb.
  • the on-time calculation unit 41 divides the period between the current detection times into three or more sections, and each of the plurality of gate pulse signals Gu , G v , and G w is turned on in each of the three or more sections. Find the on-time, which is the time that is set to.
  • the on-time calculation unit 41 calculates the on-time of each gate pulse signal G during the current detection time by adding up the on-time of each of the three or more sections for each gate pulse signal G.
  • the detection time interval voltage calculation unit 14 can accurately calculate the on-time of each gate pulse signal Gu , G v , and G w during the current detection time.
  • Embodiment 4 the configuration of the portion of the drive device of the AC motor excluding the estimation device and the voltage type inverter described above will be specifically described.
  • components having the same functions as those in the third embodiment will be designated by the same reference numerals and description thereof will be omitted, and the differences from the electric motor system 100B of the third embodiment will be mainly described.
  • FIG. 24 is a diagram showing a configuration example of the electric motor system according to the fourth embodiment of the present invention.
  • the electric motor system 100C includes an AC electric motor 1 and a drive device 2C.
  • the drive device 2C includes a voltage type inverter 4 and a control device 5C.
  • the control device 5C includes an estimation device 9B, a speed control unit 30, coordinate conversion units 31, 33, a dq axis current control unit 32, a zero vector modulation unit 34, and a gate pulse generation unit 35.
  • the control device 5C may be configured to include the estimation device 9 or the estimation device 9A instead of the estimation device 9B. It is assumed that the control device 5C is sensorless vector control based on the dq coordinate system, but other control methods or coordinate systems may be used.
  • the speed control unit 30 generates a d-axis current command id * and a q-axis current command i q * so that the estimated speed ⁇ e ⁇ calculated by the estimation device 9B matches the speed command ⁇ *.
  • the speed control unit 30 is composed of a feedback controller such as a PID (Proportional Integral Differential) controller, but may be composed of a controller other than the feedback controller. Further, the speed control unit 30 may have a feedforward controller in parallel with the PID controller, whereby the control response can be accelerated. Further, the speed control unit 30 may be configured to enhance the disturbance response of a specific frequency by using a vibration suppression controller such as repetitive control in combination with the PID controller.
  • the q-axis current command i q * is determined so that the deviation between the estimated speed ⁇ e ⁇ and the speed command ⁇ * converges to zero, but the d-axis current command i d * is determined by the type of AC motor 1.
  • the speed control unit 30 can also perform "Maximum Torque Per Ampere control".
  • Maximum torque / current control is a control system that gives a negative direction of the d-axis current command i d *, it is possible to reduce the copper loss.
  • the speed control unit 30 can also perform weak magnetic flux control, which is a control method for giving a larger negative direction d-axis current command id *. By performing weakening magnetic flux control, when the voltage is saturated, the mechanical output is increased in exchange for the increase in copper loss regardless of whether the AC motor 1 is a surface type permanent magnet synchronous motor or an embedded type permanent magnet synchronous motor. Can be made to.
  • Coordinate conversion unit 31 based on the calculated estimated position theta e ⁇ by estimator 9B, the phase currents i u by the coordinate transformation, i v, i is the current of the dq coordinate system w d-axis current i d and the q Convert to axial current i q.
  • dq-axis current control unit 32 adjusts the d-axis voltage command v d * as the d-axis current i d is coincident with the d-axis current command i d *, q-axis current i q is q-axis current command i q * Adjust the q-axis voltage command v q * so that it matches.
  • dq-axis current control unit 32 for example, the d-axis current i and not shown in the PID (Proportional Integral Differential) controller for feedback control of the d, q-axis current i q (not shown) of the PID controller for feedback control of the, dq It is provided with a non-interference control (not shown) that feedforward-compensates the interference component of the shaft.
  • PID Proportional Integral Differential
  • dq axis current control unit 32 the dq coordinate system, the d-axis current i d and the q-axis current i q are individually controlled, thereby realizes a good current control performance .
  • the speed control unit 30 and the dq-axis current control unit 32 include the calculation of the calculus
  • the dq-axis current control unit 32 uses the time difference Tb between the current detection times to prevent the calculation error of the calculus. There is.
  • the coordinate conversion unit 33 issues the d-axis voltage command v d * and the q-axis voltage command v q * to the voltage commands v u * , v v * , v w based on the estimated position ⁇ e ⁇ calculated by the estimation device 9B. Convert coordinates to *.
  • the zero vector modulation unit 34 performs zero vector modulation on the voltage commands v u * , v v * , and v w * coordinate-converted by the coordinate conversion unit 33.
  • the gate pulse generation unit 35 generates gate pulse signals Gu , G v , G w based on the voltage commands v u * , v v * , v w * , and the generated gate pulse signals G u , G v , G w is output to the voltage type inverter 4. Since the gate pulse generation unit 35 has the pulse shift processing unit 22 as described above, the three-phase current is determined by the phase current determination unit 12 regardless of the values of the voltage commands v u * , v v * , and v w *. The value of can be determined. Therefore, the estimation device 9B can further improve the accuracy of position estimation and speed estimation, and can further improve the control performance.
  • the estimation unit 15B of the estimation device 9B is based on the time difference Tb of the current detection time by the one-shunt current detection method and the voltage between the current detection times that changes with the time difference Tb of the current detection time. Speed estimation and position estimation are performed. That is, the estimation unit 15B estimates the speed based on the time difference Tb between the current detection times detected by the phase current determination unit 12 and the voltage between the current detection times calculated by the detection time voltage calculation unit 14. And position estimation. As a result, the estimation device 9B can significantly reduce the error of the speed estimation and the position estimation.
  • FIG. 25 is a diagram showing an example of the hardware configuration of the control device according to the fourth embodiment.
  • the control device 5C according to the fourth embodiment includes a bus current detection unit 11 and a processing circuit 80.
  • the processing circuit 80 includes a processor 81, a memory 82, an input circuit 83, an AD (Analog to Digital) conversion circuit 84, and a DA (Digital to Analog) conversion circuit 85.
  • the processor 81 includes one or more of a CPU (Central Processing Unit), a DSP (Digital Signal Processor), and a system LSI (Large Scale Integration).
  • the memory 82 is one or more of RAM (Random Access Memory), ROM (Read Only Memory), flash memory, EPROM (Erasable Programmable Read Only Memory), and EEPROM (registered trademark) (Electrically Erasable Programmable Read Only Memory). Including.
  • the memory 82 also includes a recording medium on which a computer-readable program is recorded. Such recording media include one or more of non-volatile or volatile semiconductor memories, magnetic disks, flexible memories, optical disks, compact disks, and DVDs (Digital Versatile Discs).
  • the processor 81 executes the functions of the phase current determination unit 12, the time difference calculation unit 13, the detection time interval voltage calculation unit 14, and the estimation unit 15B by reading and executing the program stored in the recording medium of the memory 82. To do.
  • the input circuit 83 notifies the processor 81 of an external speed command ⁇ * or a three-phase voltage command v u * , v v * , v w *.
  • the processor 81 performs a process of estimating the magnetic pole position ⁇ e and the rotation speed ⁇ e based on the speed command ⁇ * from the outside or the three-phase voltage commands v u * , v v * , v w *.
  • AD conversion circuit 84 converts the value of the bus current i bus detected by the bus current detection unit 11 from an analog signal to a digital signal.
  • the processor 81 acquires the digitized value of the bus current ibus from the AD conversion circuit 84, and performs processing based on the acquired value of the bus current ibus. Further, the processor 81 modulates the three-phase voltage commands v u * , v v * , v w * to generate gate pulse signals Gu , G v , G w, and the generated gate pulse signals Gu , G. Outputs v and G w to the DA conversion circuit 85.
  • the DA conversion circuit 85 converts the gate pulse signals Gu , G v , and G w from digital signals to analog signals, and outputs the converted gate pulse signals Gu , G v , and G w to the gate driver 7. .
  • the processor 81 can transmit the three-phase voltage commands v u * , v v * , and v w * to the gate driver 7.
  • the three-phase voltage commands v u * , v v * , v w * may be digital values or analog values.
  • Data is transmitted and received between the processor 81, the memory 82, the input circuit 83, the AD conversion circuit 84, and the DA conversion circuit 85.
  • the data is transmitted and received as a digital signal, but may be transmitted and received as an analog signal. Further, the data transmission / reception method may be parallel communication or serial communication. Further, the analog signal and the digital signal may be appropriately converted by a converter (not shown).
  • estimation devices 9 and 9A of the first and second embodiments can also be realized by the hardware configuration shown in FIG. 25, similarly to the estimation device 9B of the fourth embodiment.
  • the processor 81 can execute the functions of the phase current determination unit 12, the time difference calculation unit 13, and the estimation unit 15 by reading and executing the program stored in the recording medium in the memory 82. Further, the processor 81 can execute the functions of the phase current determination unit 12, the time difference calculation unit 13, and the estimation unit 15A by reading and executing the program stored in the recording medium in the memory 82.
  • the processor 81 includes a speed control unit 30, a coordinate conversion unit 31, 33, a dq axis current control unit 32, a zero vector modulation unit 34, and a gate pulse generation unit 35. You can also perform some or all of the functions. In this case, the processor 81 reads and executes the program stored in the memory 82, thereby causing the speed control unit 30, the coordinate conversion units 31, 33, the dq axis current control unit 32, the zero vector modulation unit 34, and the gate pulse. The function of the generation unit 35 is executed.
  • the estimation device 9B having high estimation accuracy of the magnetic pole position ⁇ e and the rotation speed ⁇ e is used. Therefore, in the drive device 2C, the control performance for the AC motor 1 is improved when the AC motor 1 is rotating at high speed or when zero vector modulation is performed.
  • the configuration shown in the above-described embodiment shows an example of the content of the present invention, can be combined with another known technique, and is one of the configurations without departing from the gist of the present invention. It is also possible to omit or change the part.

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Ac Motors In General (AREA)
  • Inverter Devices (AREA)

Abstract

推定装置(9)は、相電流判定部(12)と、時刻差算出部(13)と、推定部(15)と、を備える。相電流判定部(12)は、母線電流検出部(11)で検出される母線電流の値と、複数のゲートパルス信号の状態とに基づいて、複数の相電流の値を判定する判定処理を行う。時刻差算出部(13)は、相電流判定部(12)による前回の判定処理で用いられた母線電流の検出時刻として特定される第1の電流検出時刻と相電流判定部(12)による今回の判定処理で用いられた母線電流の検出時刻として特定される第2の電流検出時刻との時刻差を算出する。推定部(15)は、相電流判定部(12)で判定された複数の相電流の値と、時刻差算出部(13)で算出された時刻差とに基づいて、交流電動機(1)の位置および速度のうち少なくとも一方を推定する。

Description

推定装置および交流電動機の駆動装置
 本発明は、交流電動機の位置および速度のうち少なくとも一方を推定する推定装置および交流電動機の駆動装置に関する。
 誘導機および同期機といった交流電動機の駆動には、回転子の位置を示す情報が用いられる。回転子の位置を示す情報を取得するために位置センサまたは速度センサを使用すると、製造コストの増加などの問題がある。そのため、交流電動機の駆動装置に関し、位置センサまたは速度センサなどを用いない位置センサレス制御に関する多くの検討がなされている。
 また、交流電動機の駆動装置に関し、製造コストを低減するため、電流センサの個数削減についても多くの検討がなされており、安価な電流検出方式として1シャント電流検出方式が広く用いられている。1シャント電流検出方式は、インバータの直流母線に設けられる電流センサを用いて、交流電動機に流れる相電流を測定する方式である。直流母線に設けられる電流センサを用いて相電流を測定する方式は、電流センサとしてシャント抵抗が用いられることが多いため、1シャント電流検出方式と呼ばれるが、直流母線に設けられる電流センサとしてシャント抵抗以外を用いる場合も、一般に1シャント電流検出方式と呼ばれる。例えば、1シャント電流検出方式は、CT(Current Transformer)と呼ばれる変流器などのようにシャント抵抗とは異なる電流センサを用いた構成も知られている。
 特許文献1には、位置センサレス制御と1シャント電流検出方式とを併用して交流電動機を駆動する技術が開示されている。1シャント電流検出方式では、各相の電流をキャリア信号の頂点で同時に検出できない。そこで、特許文献1に記載の技術では、第1のキャリア周期の後半と第1のキャリア周期に続く第2のキャリア周期の前半とで母線電流から得られる相電流を補間処理することで、第1のキャリア周期と第2のキャリア周期との境界のタイミングであるキャリア信号の頂点での相電流を算出する。特許文献1に記載の技術では、キャリア信号の頂点での相電流に基づいて電圧指令を生成し、生成した電圧指令に基づいて交流電動機の回転子の回転状態を示す情報を算出する。回転子の回転状態を示す情報は、回転子の位置および速度のうち少なくとも1つの情報である。
特開2015-139359号公報
 しかしながら、特許文献1に記載の技術は、第1のキャリア周期の後半と第1のキャリア周期に続く第2のキャリア周期の前半とで各々検出される母線電流を用いる。インバータを構成する複数のスイッチング素子のオンおよびオフの切り替わりタイミングによっては、キャリア周期における前半と後半との両方で母線電流から3相の電流を得ることが難しい場合がある。そのため、特許文献1に記載の技術では、交流電動機の回転子の回転状態の算出に用いられる相電流の検出周期がキャリア周期の2倍以上になる。一般に、インバータのスイッチング損失はキャリア周波数の上昇とともに増加することから、交流電動機の駆動装置の冷却性能または電力効率の都合などによりキャリア周波数には上限が存在するのが一般的である。交流電動機の回転周波数がキャリア周波数に近づくと交流電動機の回転周波数に対して制御周期が長くなり、交流電動機が備える回転子の回転状態の推定精度を確保することが難しい場合がある。
 本発明は、上記に鑑みてなされたものであって、交流電動機の位置および速度のうち少なくとも一方の推定精度を向上させることができる推定装置を得ることを目的とする。
 上述した課題を解決し、目的を達成するために、本発明の推定装置は、母線電流検出部と、相電流判定部と、時刻差算出部と、推定部と、を備える。母線電流検出部は、電圧指令に基づいて生成される複数のゲートパルス信号によって駆動される電圧形インバータの直流母線に流れる電流である母線電流の値を検出する。相電流判定部は、母線電流検出部によって検出される母線電流の値と、複数のゲートパルス信号の状態とに基づいて、電圧形インバータから交流電動機へ供給される複数の相電流の値を判定する判定処理を行う。時刻差算出部は、相電流判定部による前回の判定処理で用いられた母線電流の検出時刻として特定される第1の電流検出時刻と相電流判定部による今回の判定処理で用いられた母線電流の検出時刻として特定される第2の電流検出時刻との時刻差を算出する。推定部は、相電流判定部によって判定された複数の相電流の値と、時刻差算出部によって算出された時刻差とに基づいて、交流電動機の位置および速度のうち少なくとも一方を推定する。
 本発明によれば、交流電動機の位置および速度のうち少なくとも一方の推定精度を向上させることができる、という効果を奏する。
本発明の実施の形態1にかかる推定装置を含む電動機システムの構成例を示す図 実施の形態1にかかる相電流判定部による相電流の値の判定方法を説明するための図 実施の形態1にかかる相電流判定部による相電流の値の判定方法を説明するための図 実施の形態1にかかる相電流判定部による相電流の値の判定方法を説明するための図 実施の形態1にかかる相電流判定部による相電流の判定に用いられる母線電流の検出タイミングと、キャリア波と、電圧指令と、ゲートパルス信号との関係の一例を示す図 実施の形態1にかかる交流電動機が低速で回転している場合の3相の電圧指令とキャリア波との関係の一例を示す図 実施の形態1にかかる交流電動機が高速で回転している場合の3相の電圧指令とキャリア波との関係の一例を示す図 実施の形態1にかかる3相の電圧指令が正弦波である場合の中間相の電圧指令の変化の一例を示す図 実施の形態1にかかる三次高調波重畳を用いて変形させた3相の電圧指令のうち中間相の電圧指令の変化の一例を示す図 実施の形態1にかかる零ベクトル変調部によって零ベクトル変調が行われた3相の電圧指令の一例を示す図 長方形近似による積分演算のイメージを示す図 長方形近似による積分演算のイメージを示す図 時刻差の変動を無視した場合の積分演算のイメージを示す図 微分演算の誤差を説明するための図 本発明の実施の形態2にかかる推定装置を含む電動機システムの構成例の一例を示す図 実施の形態2にかかる推定部の構成例を示す図 本発明の実施の形態3にかかる推定装置を含む電動機システムの構成例の一例を示す図 実施の形態3にかかる検出時刻間電圧演算部の構成例を示す図 実施の形態3にかかる検出時刻間電圧演算部による電流検出時刻間の3つの相電圧の算出方法を説明するための図 実施の形態3にかかる電流検出時刻間の時刻差を固定値とした場合の速度推定の結果の一例を示す図 図20に示す速度推定の結果をFFT(Fast Fourier Transform)解析した結果を示す図 実施の形態3にかかる推定部による速度推定の結果の一例を示す図 図22に示す速度推定の結果をFFT解析した結果を示す図 本発明の実施の形態4にかかる電動機システムの構成例を示す図 実施の形態4にかかる制御装置のハードウェア構成の一例を示す図
 以下に、本発明の実施の形態にかかる推定装置および交流電動機の駆動装置を図面に基づいて詳細に説明する。なお、この実施の形態によりこの発明が限定されるものではない。
実施の形態1.
 図1は、本発明の実施の形態1にかかる推定装置を含む電動機システムの構成例を示す図である。図1に示すように、実施の形態1にかかる電動機システム100は、交流電動機1と、交流電動機1を駆動する駆動装置2とを備える。
 交流電動機1は、回転子に永久磁石が設けられる永久磁石同期電動機であるが、回転子に界磁巻線が巻かれている巻線界磁式同期電動機であってもよく、回転子の突極性を利用して回転トルクを得るリラクタンス式同期電動機であってもよい。また、交流電動機1における永久磁石の配置は、埋め込み型の配置であってもよく表面型の配置であってもよい。また、ここでは、交流電動機1が3相の交流電動機である例を説明するが、交流電動機1は、3相以外の交流電動機であってもよい。例えば、交流電動機1は、2相の交流電動機であってもよく、5相の交流電動機であってもよい。
 駆動装置2は、直流電源3から供給される直流電圧を交流電圧へ変換し、変換した交流電圧を交流電動機1へ出力する電圧形インバータ4と、電圧形インバータ4を制御して交流電動機1を駆動する制御装置5とを備える。電圧形インバータ4は、図1に示すように、主回路6と、ゲートドライバ7とを備える。
 主回路6は、複数のスイッチング素子Q1,Q2,Q3,Q4,Q5,Q6を備える。主回路6において、スイッチング素子Q1,Q3,Q5は、高電位側の直流母線61に一端が接続され、スイッチング素子Q2,Q4,Q6は、低電位側の直流母線61に一端が接続される。スイッチング素子Q1とスイッチング素子Q2とは、互いの他端同士が接続され、U相のレグを構成する。スイッチング素子Q3とスイッチング素子Q4とは、互いの他端同士が接続され、V相のレグを構成する。スイッチング素子Q5とスイッチング素子Q6とは、互いの他端同士が接続され、W相のレグを構成する。
 このように、電圧形インバータ4は、U相のレグと、V相のレグと、W相のレグとを含む3相ブリッジ回路を備える。電圧形インバータ4は、複数のスイッチング素子Q1,Q2,Q3,Q4,Q5,Q6のオンとオフとの切り替えによって、任意の振幅および任意の周波数の交流電圧を出力することができる。以下において、スイッチング素子Q1,Q2,Q3,Q4,Q5,Q6の各々を区別せずに示す場合、スイッチング素子Qと記載する場合がある。
 各スイッチング素子Qは、逆並列ダイオードが内蔵されたIGBT(Insulated Gate Bipolar Transistor)であるが、逆並列ダイオードが内蔵されたMOSFET(Metal-Oxide-Semiconductor Field-Effect Transistor)であってもよい。また、ここでは、電圧形インバータ4が、出力する電圧のレベルが2つである2レベルインバータである例を説明するが、電圧形インバータ4は、出力する電圧のレベルが3つ以上であるマルチレベルインバータであってもよい。
 ゲートドライバ7は、制御装置5から出力されるゲートパルス信号G,G,Gを増幅し、増幅したゲートパルス信号G,G,Gをゲートパルス信号Gup,Gvp,Gwpとしてスイッチング素子Q1,Q3,Q5のゲートへ出力する。また、ゲートドライバ7は、制御装置5から出力されるゲートパルス信号G,G,Gを反転増幅してゲートパルス信号Gun,Gvn,Gwnを生成し、生成したGun,Gvn,Gwnをスイッチング素子Q2,Q4,Q6のゲートへ出力する。
 ゲートパルス信号Gup,Gunは、一方がオンである場合、他方はオフである。これにより、スイッチング素子Q1がオンになると、スイッチング素子Q2がオフになり、スイッチング素子Q1がオフになると、スイッチング素子Q2はオンになる。このように、ゲートパルス信号Gup,Gunによって、スイッチング素子Q1,Q2は相補的に動作する。同様に、ゲートパルス信号Gvp,Gvnによって、スイッチング素子Q3,Q4は相補的に動作し、ゲートパルス信号Gwp,Gwnによって、スイッチング素子Q5,Q6は相補的に動作する。ゲートパルス信号G,G,G,Gup,Gun,Gvp,Gvn,Gwp,Gwnは、高電位レベルの状態である場合にオンの状態であり、低電位レベルの状態である場合にオフの状態である。
 また、ゲートドライバ7は、低圧系の制御装置5と高圧系の主回路6とを絶縁する機能を有しており、主回路6の異常時に制御装置5の故障を防ぐ役割を有している。以下において、ゲートパルス信号G,G,Gの各々を区別せずに示す場合、ゲートパルス信号Gと記載する場合がある。
 制御装置5は、零ベクトル変調を行う零ベクトル変調部34と、ゲートパルス信号G,G,Gを生成するゲートパルス生成部35と、交流電動機1の回転子の磁極位置θおよび回転速度ωを推定する推定装置9とを備える。磁極位置θは、交流電動機1が備える回転子の電気角であり、交流電動機1の位置の一例である。また、回転速度ωは、交流電動機1が備える回転子の電気角速度であり、交流電動機1の速度の一例である。以下において、磁極位置θの推定を位置推定と記載する場合があり、回転速度ωの推定を速度推定と記載する場合がある。
 零ベクトル変調部34は、電圧形インバータ4から出力される2種類の零電圧ベクトルの出力比率を不規則に変化させる。2種類の零電圧ベクトルは、第1の零電圧ベクトルと第2の零電圧ベクトルである。第1の零電圧ベクトルは、上側アームであるスイッチング素子Q1,Q3,Q5がすべてオンである状態で電圧形インバータ4から出力される。第2の零電圧ベクトルは、上側アームであるスイッチング素子Q1,Q3,Q5がすべてオフである状態で電圧形インバータ4から出力される。
 第1の零電圧ベクトルと第2の零電圧ベクトルとの出力比率を不規則に変化させることで、キャリア騒音のスペクトルピークが分散され、キャリア騒音を低減することができる。キャリア騒音は、後述するキャリア波Scの周波数であるキャリア周波数に応じて交流電動機1または電圧形インバータ4などが振動して発する騒音である。第1の零電圧ベクトルと第2の零電圧ベクトルとの出力比率は、ゲートパルス生成部35による変調方式によって変わる。変調方式が変わる場合でも、3相の電圧指令v ,v ,v に同じ値を加えることによって、第1の零電圧ベクトルと第2の零電圧ベクトルとの出力比率を変化させることができる。電圧指令v は、u相の電圧指令であり、電圧指令v は、v相の電圧指令であり、電圧指令v は、w相の電圧指令である。
 3相の電圧指令v ,v ,v のすべてに同じ正の値を加えることによってゲートパルス信号G,G,Gがすべてオンになる時間が延び、電圧形インバータ4から第1の零電圧ベクトルが出力される時間が延びる。また、3相の電圧指令v ,v ,v のすべてに同じ負の値を加えることによってゲートパルス信号G,G,Gがすべてオフになる時間が延び、電圧形インバータ4から第2の零電圧ベクトルが出力される時間が延びる。
 零ベクトル変調部34は、3相の電圧指令v ,v ,v に乱数値を加えることによって、第1の零電圧ベクトルと第2の零電圧ベクトルとの出力比率を不規則に変化させる零ベクトル変調を行う。なお、零ベクトル変調部34は、予め設定された条件を満たさない場合、入力される3相の電圧指令v ,v ,v をそのままゲートパルス生成部35へ出力することができる。予め設定された条件は、例えば、交流電動機1が予め設定された速度以下で回転していることなどである。また、零ベクトル変調部34は、零ベクトル変調を行わない設定がされている場合、入力される3相の電圧指令v ,v ,v をそのままゲートパルス生成部35へ出力することができる。
 ゲートパルス生成部35は、電圧指令v ,v ,v に基づいて、ゲートパルス信号G,G,Gを生成する。ゲートパルス生成部35は、生成したゲートパルス信号G,G,Gを電圧形インバータ4のゲートドライバ7へ出力する。
 ゲートパルス生成部35は、高周波の周期信号であるキャリア波Scと3相の電圧指令v ,v ,v とを比較する比較部21と、比較部21から出力されるゲートパルス信号G,G,Gのシフト処理を行うパルスシフト処理部22とを備える。ここでは、キャリア波Scが三角波の信号である例を説明するが、キャリア波Scは、鋸波などの他の波形の信号であってもよい。キャリア波Scは、キャリア信号とも呼ばれる。
 比較部21は、電圧指令v の瞬時値がキャリア波Scの瞬時値以下である場合、ゲートパルス信号Gをオンにし、電圧指令v の瞬時値がキャリア波Scの瞬時値よりも大きい場合、ゲートパルス信号Gをオフにする。比較部21は、電圧指令v の瞬時値がキャリア波Scの瞬時値以下である場合、ゲートパルス信号Gをオンにし、電圧指令v の瞬時値がキャリア波Scの瞬時値よりも大きい場合、ゲートパルス信号Gをオフにする。比較部21は、電圧指令v の瞬時値がキャリア波Scの瞬時値以下である場合、ゲートパルス信号Gをオンにし、電圧指令v の瞬時値がキャリア波Scの瞬時値よりも大きい場合、ゲートパルス信号Gをオフにする。
 なお、上述した例では、ゲートパルス生成部35は、キャリア比較変調方式を用いてゲートパルス信号G,G,Gを生成するが、ゲートパルス信号G,G,Gを生成する方法はキャリア比較変調方式に限定されない。例えば、ゲートパルス生成部35は、キャリア比較変調方式に代えて、空間ベクトル変調方式などの他変調方式を用いて、ゲートパルス信号G,G,Gを生成してもよい。また、ゲートパルス生成部35は、ゲートパルス信号G,G,Gに代えて、ゲートパルス信号Gup,Gun,Gvp,Gvn,Gwp,Gwnを生成してもよい。この場合、ゲートドライバ7は、ゲートパルス信号Gup,Gun,Gvp,Gvn,Gwp,Gwnを増幅する。ゲートドライバ7は、増幅したゲートパルス信号Gup,Gun,Gvp,Gvn,Gwp,Gwnをスイッチング素子Q1,Q2,Q3,Q4,Q5,Q6のゲートへ出力する。
 パルスシフト処理部22は、ゲートパルス信号G,G,Gのオンとオフとの切り替わりタイミングをずらすパルスシフト処理を行う。パルスシフト処理部22のパルスシフト処理については、後で詳述する。
 推定装置9は、ゲートパルス信号G,G,Gと、電圧指令v ,v ,v とに基づいて、交流電動機1の磁極位置θおよび回転速度ωを推定する。推定装置9は、交流電動機1の磁極位置θおよび回転速度ωのうち一方のみを推定することもできる。
 推定装置9は、母線電流ibusの値を検出する母線電流検出部11と、相電流i,i,iの値を判定する相電流判定部12と、ゲートパルス信号G,G,Gに基づいて、後述する時刻差Tbを算出する時刻差算出部13とを備える。また、推定装置9は、時刻差算出部13によって算出された時刻差Tbと、相電流判定部12によって判定された相電流i,i,iの値と、電圧指令v ,v ,v とに基づいて、磁極位置θおよび回転速度ωを推定する推定部15を備える。
 母線電流検出部11によって検出される母線電流ibusの値は、直流電源3と電圧形インバータ4との間の直流母線61に流れる電流である母線電流の瞬時値である。母線電流検出部11は、CTと呼ばれる変流器を用いたタイプの電流センサであってもよく、シャント抵抗を用いるタイプの電流センサであってもよい。以下において、母線電流ibusの値を単に母線電流ibusと記載する場合がある。
 図1に示す例では、母線電流検出部11は、低電位側の直流母線61に設けられ、低電位側の直流母線61に流れる電流の値を検出するが、高電位側の直流母線61に設けられてもよい。なお、母線電流検出部11がシャント抵抗を用いるタイプの電流センサである場合、低電位側の直流母線61に母線電流検出部11を設けることで、母線電流検出部11における絶縁回路の部品コストを抑えることができる。
 相電流判定部12は、母線電流検出部11によって検出される母線電流ibusの値に基づいて、相電流i,i,iの値を判定する。相電流判定部12によって判定される相電流i,i,iの値は、電圧形インバータ4と交流電動機1との間に流れるu相、v相、およびw相の電流である3つの相電流の瞬時値である。以下、相電流iをu相電流iと記載し、相電流iをv相電流iと記載し、相電流iをw相電流iと記載する場合がある。また、相電流i,i,iの値を単に相電流i,i,iと記載する場合がある。
 図2から図4は、実施の形態1にかかる相電流判定部による相電流の値の判定方法を説明するための図である。図2から図4に示す例では、電圧形インバータ4にY結線の3相抵抗負荷1aが接続されている。また、スイッチング素子Q1,Q2の接続点から3相抵抗負荷1aへ流れる電流がu相電流iであり、スイッチング素子Q3,Q4の接続点から3相抵抗負荷1aへ流れる電流がv相電流iである。また、スイッチング素子Q5,Q6の接続点から3相抵抗負荷1aへ流れる電流がw相電流iである。
 また、電圧形インバータ4から3相抵抗負荷1aへ流れる相電流の向きは、正方向であり、3相抵抗負荷1aから電圧形インバータ4へ流れる相電流の向きは、負方向である。例えば、u相電流iが図2に示す矢印の方向に流れる場合、u相電流iの向きは正方向である。同様に、v相電流iが図2に示す矢印の方向に流れる場合、v相電流iの向きは正方向であり、w相電流iが図2に示す矢印の方向に流れる場合、w相電流iの向きは正方向である。
 相電流判定部12は、母線電流検出部11によって検出される母線電流ibusの値と、複数のゲートパルス信号G,G,Gの状態とに基づいて、相電流i,i,iの値を判定する。例えば、相電流判定部12は、6つのスイッチング素子Qのオンとオフとの組み合わせパターンが特定の組み合わせパターンであるタイミングで母線電流検出部11によって検出される母線電流ibusの値に基づいて、相電流i,i,iの値を判定する。
 例えば、スイッチング素子Q1,Q4,Q6がオンであり、スイッチング素子Q2,Q3,Q5がオフであるとする。この場合、図3に示すように、電圧形インバータ4から3相抵抗負荷1aへu相電流iが流れ、u相電流iと同じ大きさの電流がv相電流iおよびw相電流iに分割されて3相抵抗負荷1aから電圧形インバータ4へ流れる。直流母線61に流れる電流は、u相電流iと同じ大きさの電流であるため、相電流判定部12は、図3に示す状態において、母線電流検出部11によって検出される母線電流ibusの値によってu相電流iの値を判定することができる。
 図3に示すように、u相電流iの向きに対して直流母線61に流れる電流の向きが逆方向になる。そのため、母線電流検出部11は、直流母線61に流れる電流が図3に示す向きである場合に、正の母線電流ibusを出力するように直流母線61に配置される。なお、母線電流検出部11は、直流母線61に流れる電流が図3に示す向きである場合に、負の母線電流ibusを出力するように直流母線61に配置されてもよい。この場合、母線電流ibusの値は、相電流判定部12によって正負が反転される。
 図3では、正方向のu相電流iの値を判定する例を示しているが、負方向のu相電流iの値、正方向および負方向のv相電流iの値、および正方向および負方向のw相電流iの値も同様に判定される。例えば、図4に示すように、スイッチング素子Q1,Q3,Q6がオンであり、スイッチング素子Q2,Q4,Q5がオフであるとする。この場合、相電流判定部12は、母線電流検出部11によって検出される母線電流ibusの値によって負方向のw相電流iの値を判定することができる。
 3つの相電流i,i,iのうち2つの相電流の値がわかれば、残りの1つの相電流の値は、キルヒホッフの電流則によって算出することができる。そのため、相電流判定部12は、母線電流検出部11によって異なる2つのタイミングで検出される母線電流ibusの値から3つの相電流i,i,iのうちいずれか2つの相電流の値を判定する。そして、相電流判定部12は、判定した2つの相電流の値から残りの1つの相電流の値を算出により判定する。このように、相電流判定部12は、2回の母線電流検出を1セットとして、3つの相電流i,i,iの値を判定する。
 電圧形インバータ4が2レベルインバータである場合、6つのスイッチング素子Qにおけるオンとオフとの組み合わせのパターンは、8つある。8つの組み合わせパターンのうち6つの組み合わせパターンの各々は、上述した特定パターンである。6つのスイッチング素子Qにおけるオンとオフとの組み合わせのパターンが特定パターンである場合、母線電流検出部11によって検出される母線電流ibusの値から相電流i,i,iのいずれか1つの値を判定することができる。残りの2つの組み合わせパターンでは、電圧形インバータ4の出力は上述した零電圧ベクトルになるため、母線電流検出部11によって検出される母線電流ibusの値から相電流i,i,iの値を判定することが難しい。
 図5は、実施の形態1にかかる相電流判定部による相電流の判定に用いられる母線電流の検出タイミングと、キャリア波と、電圧指令と、ゲートパルス信号との関係の一例を示す図である。図5において、Tcは、キャリア波Scの周期であるキャリア周期であり、fcは、キャリア波Scの周波数であるキャリア周波数である。
 図5に示す例では、3相の電圧指令v ,v ,v のうち、電圧指令v の瞬時値が最も大きく、電圧指令v の瞬時値がその次に大きく、電圧指令v の瞬時値が最も小さい。本実施の形態では、ある時刻における3相の電圧指令v ,v ,v のうち絶対値が中間にある電圧指令の相を中間相と呼ぶ。図5に示す例では、電圧指令v が中間の大きさであるため、中間相は、V相である。
 図5において、時刻t1からt7の期間は、キャリア波Scの最初の下り半周期である。下り半周期は、キャリア波Scの半周期であってキャリア波Scの値が漸次小さくなる期間である。キャリア波Scの最初の下り半周期のうち、ゲートパルス信号Gがオンでゲートパルス信号G,Gがオフである時刻t3では、スイッチング素子Q1,Q4,Q6がオンで、スイッチング素子Q2,Q3,Q5がオフである。したがって、時刻t3において、母線電流検出部11で検出される母線電流ibusの値は、u相電流iの値と同じである。相電流判定部12は、時刻t3で母線電流検出部11によって検出される母線電流ibusの値をu相電流iの値として判定する。
 キャリア波Scの最初の下り半周期のうち、ゲートパルス信号G,Gがオンでゲートパルス信号Gがオフである時刻t5では、スイッチング素子Q1,Q3,Q6がオンで、スイッチング素子Q2,Q4,Q5がオフである。したがって、時刻t5において、母線電流検出部11で検出される母線電流ibusの値はw相電流iの値と同じである。相電流判定部12は、時刻t5で母線電流検出部11によって検出される母線電流ibusの値をw相電流iの値として判定する。
 相電流判定部12は、時刻t3で検出された母線電流ibusの値から得られるu相電流iの値と、時刻t5で検出された母線電流ibusの値から得られるw相電流iの値とから、キルヒホッフの法則に基づき、v相電流iの値を算出する。このように、相電流判定部12は、キャリア波Scの下り半周期において複数の電流検出タイミングで母線電流検出部11によって検出される母線電流ibusの値から相電流i,i,iの値を判定する。同様に、相電流判定部12は、次の下り半周期である時刻t11からt16の期間のうち時刻t13,t15を電流検出タイミングとして母線電流検出部11によって検出される母線電流ibusの値を用いて相電流i,i,iの値を判定する。
 このように、相電流判定部12は、キャリア波Scの下り半周期において、2回の電流検出タイミングで母線電流検出部11によって検出される母線電流ibusの値から相電流i,i,iの値を判定する。なお、これらの2回の電流検出タイミングは、キャリア波Scの下り半周期に限定されない。例えば、相電流判定部12は、キャリア波Scの下り半周期に代えて、キャリア波Scの上り半周期において、2回の電流検出タイミングで母線電流検出部11によって検出される母線電流ibusの値から相電流i,i,iの値を判定することもできる。上り半周期は、キャリア波Scの半周期であってキャリア波Scの値が漸次大きくなる期間である。
 スイッチング素子Qのオンとオフが切り替わった直後は、スイッチング素子Qのオンとオフとの切り替えによって母線電流ibusにリンギングが現われる。リンギングが生じている母線電流ibusの値から相電流i,i,iの値を精度よく判定することが難しい。そのため、相電流判定部12は、リンギングが収まるまでの時間として予め定められた時間待ってから、母線電流検出部11によって検出される母線電流ibusの値を用いて相電流の値を判定する。
 図5に示す例では、相電流判定部12は、中間相のゲートパルス信号であるゲートパルス信号Gがオフからオンになる時刻t4の直前である時刻t3で母線電流検出部11によって検出される母線電流ibusの値を相電流iの値として判定する。また、相電流判定部12は、ゲートパルス信号Gがオフからオンになる時刻t4から予め設定された時間TA後の時刻t5で母線電流検出部11によって検出される母線電流ibusの値を相電流iの値として判定する。
 このように、相電流判定部12は、ゲートパルス信号G,G,Gのうち中間相のゲートパルス信号のオンとオフとの切り替わりであるタイミングである状態切り替わりタイミングを判定する。そして、相電流判定部12は、状態切り替わりタイミングの直前のタイミングと状態切り替わりタイミングから時間TA後のタイミングとを各々母線電流ibusの値を検出するタイミングである母線電流検出タイミングとして決定する。相電流判定部12は、これら2回の母線電流検出タイミングの各々で母線電流検出部11によって検出される母線電流ibusの値から、相電流i,i,iのいずれか2つの相電流の値を判定する。相電流判定部12は、判定した2つの相電流の値に基づいて残りの相電流の値を算出する。2つの相電流の判定タイミングを近づけることができることから、相電流判定部12は、検出した2つの相電流の値に基づいて残りの相電流の値を精度よく判定することができる。
 相電流判定部12は、母線電流検出部11によって繰り返し検出される母線電流ibusの値を記憶することができる。相電流判定部12は、記憶した母線電流ibusの値のうち、中間相のゲートパルス信号が切り替わる直前に母線電流検出部11によって検出された母線電流ibusの値を抽出することができる。相電流判定部12は、抽出した母線電流ibusの値に基づいて、相電流の値を判定することができる。
 また、中間相のゲートパルス信号がオフからオンになってから中間相のレグの上側アームがオフからオンになるまでに遅延時間Tdがあるとする。この場合、相電流判定部12は、中間相のゲートパルス信号のオンとオフとが切り替わった時刻に母線電流検出部11によって検出された母線電流ibusの値を用いることもできる。相電流判定部12は、中間相のゲートパルス信号のオンとオフとが切り替わった時刻に検出される母線電流ibusの値に加え、中間相のゲートパルス信号のオンとオフとが切り替わってから時間TC後に検出される母線電流ibusの値を用いる。時間TCは、時間TAに遅延時間Tdを加算した時間である。なお、上述した例では、中間相がv相である例を示したが、相電流判定部12は、中間相がu相である場合および中間相がw相である場合も同様の処理を行うことができる。
 このように、相電流判定部12は、複数のゲートパルス信号G,G,Gのうち中間相のゲートパルス信号が変化するタイミングに基づき2回の母線電流検出タイミングを決定する。これら2回の母線電流検出タイミングは、相電流i,i,iのうち2つの相電流の値を判定するために用いられる母線電流ibusの値を検出するタイミングである。相電流判定部12は、決定した2回の母線電流検出タイミングを含む検出期間毎に母線電流検出部11によって検出される複数の母線電流ibusの値に基づいて、相電流i,i,iの値を判定する。
 上述した8つの組み合わせパターンは、電圧指令v ,v ,v の変化によって切り替わり、組み合わせパターンの切り替わり間隔も、電圧指令v ,v ,v の変化によって変わる。組み合わせパターンの切り替わりは、3相のゲートパルス信号G,G,Gのいずれかのオンとオフとが切り替わることによって行われ、その直後に母線電流ibusにリンギングが現われる。
 リンギングは一般に数マイクロ秒オーダで収束するが、リンギングが収束する前に再度、3相のゲートパルス信号G,G,Gのいずれかのオンとオフとが切り替わると、切り替わり前の状態における相電流の値を得ることが難しい。すなわち、数マイクロ秒の範囲において、3相のゲートパルス信号G,G,Gで構成されるゲートパルス信号群の状態変化が2回発生する場合、母線電流検出部11によって検出される母線電流ibusの値から相電流の値を判定することが難しい。このように、組み合わせパターンの切り替わり間隔がリンギングの収束時間よりも短い場合、母線電流検出部11によって検出される母線電流ibusの値から相電流の値を判定することが難しい。
 そこで、制御装置5のパルスシフト処理部22は、組み合わせパターンの切り替わり間隔が予め設定された時間TA以上になるように、3相のゲートパルス信号G,G,Gの少なくとも1つのオンとオフとが切り替わるタイミングをずらす処理を行う。時間TAは、リンギングが生じてからリンギングが収束するまでの時間、すなわち、リンギングの収束時間以上の値に設定される。これにより、パルスシフト処理部22は、組み合わせパターンの切り替わり間隔をリンギングの収束時間以上ずらすことができ、電圧指令v ,v ,v の値にかかわらず、リンギングが収束するまでの待ち時間を確保することができる。
 パルスシフト処理部22は、キャリア1周期内で、3相のゲートパルス信号G,G,Gの少なくとも1つのオンとオフとが切り替わるタイミングをずらすパルスシフト処理を行う。例えば、パルスシフト処理部22は、キャリア波Scの上り半周期から下り半周期までの期間において、ゲートパルス信号Gをずらし、キャリア波Scの上り半周期と下り半周期とで各々ゲートパルス信号Gのデューティ比を変える処理を行う。
 パルスシフト処理によって、キャリア波Scの半周期では電圧形インバータ4から出力される3相の電圧が変化するが、キャリア1周期では、オンとオフとが切り替わるタイミングをずらす前と後とでゲートパルス信号Gのデューティ比は同じである。これにより、電圧形インバータ4から出力される3相の電圧のキャリア1周期毎の平均値が電圧指令v ,v ,v と一致するようにスイッチング素子Q1,Q2,Q3,Q4,Q5,Q6のオンおよびオフが行われる。したがって、パルスシフト処理部22は、キャリア1周期単位で電圧形インバータ4の3相の出力電圧の変化を抑制することができる。
 相電流判定部12は、パルスシフト処理によって、電圧指令v ,v ,v の値とは無関係に母線電流検出部11によって検出される母線電流ibusの値から相電流i,i,iの値を判定することができる。
 図5に示す例では、相電流判定部12は、キャリア波Scの下り半周期において、3つの相電流i,i,iの値を判定するが、3つの相電流i,i,iの値を判定するタイミングは、キャリア波Scの下り半周期に限定されない。例えば、相電流判定部12は、キャリア波Scの下り半周期に代えて、キャリア波Scの上り半周期において、2回の母線電流検出タイミングで母線電流検出部11によって検出される母線電流ibusの値に基づき相電流i,i,iの値を判定できる。
 また、相電流判定部12は、キャリア波Scの下り半周期での相電流i,i,iの値を判定する判定処理と、キャリア波Scの上り半周期での相電流i,i,iの値を判定する判定処理とを交互に行うこともできる。この場合、相電流判定部12は、パルスシフト処理部22によってキャリア1周期単位で電圧形インバータ4の出力電圧の変化が抑制されるため、キャリア波Scの周期の1.5倍以上の周期毎に1度、相電流i,i,iの値を判定する。
 また、相電流判定部12は、キャリア波Scの1周期毎に3つの相電流i,i,iの値を判定することに代えて、キャリア波Scの2周期以上の期間毎に3つの相電流i,i,iの値を判定することもできる。キャリア波Scの周波数に上限がある場合、3つの相電流i,i,iの値を判定する判定処理の周期である3相電流判定周期が長くなると、交流電動機1の高速回転時において、相電流i,i,iの周波数がキャリア波Scの周波数へ近づく。この場合、相電流判定部12によって判定される相電流i,i,iの1周期分の波形に対する時間分解能が低くなる。時間分解能が低くなるほど、推定装置9の推定精度および制御装置5による制御性能などが下がることから、高速回転時には3相電流判定周期は長くない方がよい。したがって、キャリア波Scの周波数に上限がある場合、3相電流判定周期は、キャリア波Scの1周期であることが望ましい場合がある。
 上述した中間相のゲートパルス信号の状態切り替わりタイミングは、中間相の電圧指令によって変化し、中間相のゲートパルス信号の状態切り替わりタイミングによって、母線電流ibusの値の検出タイミングである母線電流検出タイミングが変化する。また、母線電流検出タイミングは、パルスシフト処理部22によるパルスシフト処理によっても変化する。パルスシフト処理による母線電流検出タイミングの変化は小さいが、中間相の電圧指令が大きく変動する場合、母線電流検出タイミングも大きく変動する。
 図5に示す例では、キャリア波Scが極大値になる時刻t11を境に中間相の電圧指令であるv相の電圧指令vv*の値が大きく変化している。そのため、中間相のゲートパルス信号であるゲートパルス信号Gが切り替わるタイミングも大きく前後する。それに伴い、キャリア波Scが極大値になる時刻t1から母線電流検出タイミングである時刻t3,t5までの時間と、キャリア波Scが極大値になる時刻t11から母線電流検出タイミングである時刻t13,t15までの時間との差が大きくなる。
 このように、母線電流検出タイミングは、キャリア周期毎に変化する。そのため、推定部15は、相電流判定部12によって判定される相電流i,i,iの値を用いて磁極位置θおよび回転速度ωを推定する際に、電流検出時刻間の時刻差Tbを用いている。電流検出時刻は、相電流判定部12による判定処理に用いられる母線電流ibusの母線電流検出部11による検出時刻として特定される時刻である。例えば、電流検出時刻は、相電流判定部12が1相目の電流の値を判定するために用いた母線電流ibusの検出時刻と相電流判定部12が2相目の電流の値を判定するために用いた母線電流ibusの検出時刻との平均の時刻である。この場合、相電流判定部12は、1相目の電流の値を判定するために用いた母線電流ibusの検出時刻と、2相目の電流の値を判定するために用いた母線電流ibusの検出時刻とを加算して2で除して得られる時刻を電流検出時刻として求める。母線電流ibusの検出時刻は、母線電流ibusが母線電流検出部11によって検出された時刻である。
 なお、電流検出時刻は、1相目の電流の値を判定するために用いられた母線電流ibusの母線電流検出部11による検出時刻であってもよく、2相目の電流の値を判定するために用いられた母線電流ibusの母線電流検出部11による検出時刻であってもよい。
 例えば、図5において、最初にu相電流iの値を判定するために用いられる母線電流ibusが検出される時刻t3を時刻tu1とし、次にu相電流iの値を判定するために用いられる母線電流ibusが検出される時刻t13を時刻tu2とする。また、最初にw相電流iの値を判定するために用いられる母線電流ibusが検出される時刻t5を時刻tw1とし、次にw相電流iの値を判定するために用いられる母線電流ibusが検出される時刻t15を時刻tw2とする。リンギングの収束のための待ち時間である時間TAは既知である。したがって、時刻tu1と時刻tw1との時刻差は既知であり、時刻tu2と時刻tw2との時刻差も既知である。しかしながら、時刻tu1と時刻tu2との時刻差または時刻tw1と時刻tw2との時刻差は、既知ではない。すなわち、電流検出時刻間の時刻差Tbは既知ではない。
 そこで、時刻差算出部13は、ゲートパルス信号G,G,Gに基づいて、電流検出時刻間の時刻差Tbを算出する。具体的には、時刻差算出部13は、ゲートパルス信号G,G,Gのうち中間相のゲートパルス信号がオンとオフとの間で切り替わるタイミングに基づいて、時刻tu1,tw1,tu2,tw2を判定する。時刻差算出部13は、判定した時刻tu1,tw1,tu2,tw2に基づいて、時刻tu1と時刻tw1との中間の時刻である時刻tavg1を算出し、時刻tu2と時刻tw2との中間の時刻である時刻tavg2を算出する。時刻差算出部13は、時刻tavg1と時刻tavg2との時刻差を時刻差Tbとして算出する。このように、時刻差算出部13は、相電流判定部12による判定処理で用いられた母線電流ibusの検出時刻として特定される電流検出時刻と相電流判定部12による次の判定処理で用いられた母線電流ibusの検出時刻として特定される電流検出時刻との時刻差Tbを算出する。
 リンギングの収束のための待ち時間である時間TAが一定である場合、時刻tu1と時刻tu2との時刻差および時刻tw1と時刻tw2との時刻差の各々は、時刻tavg1と時刻tavg2との時刻差と同じである。この場合、時刻差算出部13は、時刻tu1と時刻tu2との時刻差を時刻差Tbとして算出することもでき、時刻tw1と時刻tw2との時刻差を時刻差Tbとして算出することもできる。図5では、時刻tw1と時刻tw2との時刻差が時刻差Tbとして算出された例を示している。
 時刻差Tbの変動についてさらに具体的に説明する。図6は、実施の形態1にかかる交流電動機が低速で回転している場合の3相の電圧指令とキャリア波との関係の一例を示す図である。図7は、実施の形態1にかかる交流電動機が高速で回転している場合の3相の電圧指令とキャリア波との関係の一例を示す図である。図6および図7において、縦軸は変調率を示し、横軸は時刻を示す。また、図6および図7に示す3相の電圧指令v ,v ,v は、零ベクトル変調部34による零ベクトル変調は行われていない。
 図6に示すように、交流電動機1が低速で回転している場合、キャリア波Scの周波数に対して3相の電圧指令v ,v ,v の周波数が十分に低い。また、交流電動機1の速度起電力は交流電動機1の速度に比例して大きくなるため、交流電動機1が低速で回転している場合、3相の電圧指令v ,v ,v の振幅は比較的小さい。したがって、3相の電圧指令v ,v ,v は、傾きが緩やかであり、中間相の電圧指令は、傾きが緩やかである。なお、図6に示す例では、中間相は、前半から後半の一部までがv相で、後半の残りがu相である。
 図7に示すように、交流電動機1が高速で回転している場合、交流電動機1が低速で回転している場合に比べ、3相の電圧指令v ,v ,v の周波数が高い。例えば、非同期PWM(Pulse Width Modulation)で交流電動機1を回転させる場合、キャリア波Scの周波数に対する3相の電圧指令v ,v ,v の周波数の比率が10%~15%程度になることは珍しくない。交流電動機1が高速で回転している場合、図7に示すように、交流電動機1が低速で回転している場合に比べて、中間相がv相、u相、w相、v相、・・・の順で高速で変化し、速度起電力も大きくなっている。そのため、3相の電圧指令v ,v ,v の振幅も、図6に示す場合に比べて大きく、中間相の電圧指令の傾きもかなり大きくなっている。したがって、電流検出時刻間の時刻差Tbの変動は大きい。
 交流電動機1を高速で回転させる場合、相電流i,i,iの波形に対する時間分解能は高いほどよい。そのため、図5に示す例では、相電流判定部12は、キャリア波Scの1周期毎に母線電流検出部11によって検出される母線電流ibusの値に基づいて3相の相電流i,i,iの値を判定する処理を行う。電流検出時刻間の平均時刻差に対して、電流検出時刻間の時刻差Tbの変動が最も大きくなるのは、キャリア波Scの1周期毎に1度、3つの相電流i,i,iの値を判定する場合である。この場合、電流検出時刻間の時刻差Tbは、電流検出時刻間の平均時刻差を「Tbave」とすると、0.5×Tbave<Tb<1.5×Tbaveの範囲で変化する。
 キャリア波Scの2周期毎に1度、3つの相電流i,i,iの値を判定する場合、電流検出時刻間の時刻差Tbの変動は、0.75×Tbave<Tb<1.25×Tbaveの範囲に緩和される。「Tbave」は、上述したように電流検出時刻間の平均時刻差である。電流検出時刻間の平均時刻差を延ばすほど、電流検出時刻間の時刻差Tbの変動は緩和されるが、相電流i,i,iの波形に対する時間分解能が低くなる。そのため、交流電動機1が高速で回転している場合に回転速度ωの推定精度を高くすることが難しくなる場合がある。
 図8は、実施の形態1にかかる3相の電圧指令が正弦波である場合の中間相の電圧指令の変化の一例を示す図である。図9は、実施の形態1にかかる三次高調波重畳を用いて変形させた3相の電圧指令のうち中間相の電圧指令の変化の一例を示す図である。図8に示すように、3相の電圧指令v ,v ,v が正弦波である場合、中間相の電圧指令の波形は、ほぼ三角波である。また、図9に示すように、3相の電圧指令v ,v ,v が正弦波に三次高調波を重畳した電圧指令である場合も、中間相の電圧指令の波形は、ほぼ三角波である。また、他の変調方式を用いた場合、中間相波形が図8および図9に示す程度の三角波にならないこともあるが、変調率が1の付近では零電圧ベクトルの選択の自由度が減少することから、中間相の波形は三角波にかなり近づく。
 これらの三角波の周波数は、図8および図9に示すように、3相の電圧指令v ,v ,v の3倍である。そのため、中間相の波形は3相の電圧指令v ,v ,v の基本周波数に対して3n次の周波数成分を含んでいる。nは1以上の奇数である。
 電流検出時刻間の時刻差Tbの変動は、交流電動機1を高速で回転させる場合だけではなく、交流電動機1を低速で回転させる場合においても大きくなる場合がある。図10は、実施の形態1にかかる零ベクトル変調部によって零ベクトル変調が行われた3相の電圧指令の一例を示す図である。図10において、縦軸は変調率を示し、横軸は時刻を示す。図10に示すv ,v ,v は、図9に示す三次高調波が重畳された電圧指令v ,v ,v の波形に対して、ランダムノイズが加えられたような波形である。
 上述したように、零ベクトル変調によりキャリア騒音のスペクトルピークが分散され、キャリア騒音を低減することができる。零ベクトル変調を行った場合、交流電動機1が低速で回転している場合でも中間相の波形が急激に変化する場合がある。交流電動機1が低速で回転している場合、電圧形インバータ4から出力される電圧が小さく、電圧形インバータ4から出力される電圧が小さいときほど、零電圧ベクトルの選択の自由度が増す。そのため、交流電動機1が低速で回転している場合、零ベクトル変調で電圧指令v ,v ,v に加算できる乱数値は電圧指令v ,v ,v の変調率が小さいときほど大きい。零ベクトル変調部34は、交流電動機1の回転が低速であるほど、電圧指令v ,v ,v に加算する乱数値を大きくする。
 したがって、零ベクトル変調を行った場合、交流電動機1が低速で回転しているときでも中間相の電圧指令の波形が急激に変化し、電流検出時刻間の時刻差Tbの変化が大きくなる。この場合、電流検出時刻間の時刻差Tbの変動を無視すると速度推定または位置推定を行っても良好な結果が得られない場合がある。
 ここで、電流検出時刻間の時刻差Tbの変化に伴う微積分演算の誤差について説明する。まず、電流検出時刻間の時刻差Tbの変化に伴う積分演算の誤差について説明する。図11および図12は、長方形近似による積分演算のイメージを示す図である。図13は、時刻差の変動を無視した場合の積分演算のイメージを示す図である。図11から図13において、縦軸は、任意の曲線x(t)の時刻tにおける値を示し、横軸は、時刻を示す。
 任意の曲線x(t)を離散的に積分する場合、幅が各々微小時間Tである複数の長方形を幅方向に敷き詰め、各長方形の面積を合算することによって積分結果が得られる。この場合、各長方形は、他の長方形と重ならず且つ隣接する長方形との間で隙間がないように敷き詰められる。図11は、微小時間Tが一定の場合における長方形近似による積分の例を示しており、積分演算の誤差が生じるが、多くの場合は長方形近似による積分でも問題がない。より高精度な積分演算を行いたい場合、台形近似による積分を行うことで対応できる。
 微小時間Tが一定でない場合も同様に、図12に示すように、長方形と重ならず且つ隣接する長方形との間で隙間がないように各長方形を敷き詰めることで、適切な積分演算を行うことができる。図12に示す例では、T=T+ΔT(t)であり、微小時間Tは一定であり、ΔT(t)が時刻の経過に伴って変動する時間である。
 ΔT(t)を用いずに積分演算を行った場合、図13に示すように、他の長方形と重なったり、隣接する長方形との間で隙間が生じたりする。全体の面積については、重なった分と隙間の分が相殺しあって大きな誤差にならない場合があるが、部分的な面積には大きな誤差が発生する。したがって、電流検出時刻間の時刻差Tbを固定値として積分演算を行うと積分演算の結果に誤差が生じる。
 次に、電流検出時刻間の時刻差Tbの変化に伴う微分演算の誤差について説明する。 図14は、微分演算の誤差を説明するための図である。図14に示す例では、任意の曲線x2(t)の傾きを求めるのに使用できる値は、時刻t1におけるx2(t1)=aと、時刻t2におけるx2(t2)=bの2点であり、時刻t1と時刻t2との時刻差Tbは、Tb=T+ΔT(t)であるとする。また、ΔT(t)は時刻差Tbの変動量であり、計算するタイミング毎に変化するものとする。この場合、ΔT(t)を用いずに微分演算を行った場合、ΔT(t)の分だけ微分演算の結果に誤差が生じる。
 このように、電流検出時刻間の時刻差Tbが変動する場合において、電流検出時刻間の時刻差Tbの変動を無視して微積分演算を行う場合、微積分演算の誤差が生じる。そこで、推定部15は、時刻差算出部13によって算出された時刻差Tbに基づいて、交流電動機1の磁極位置θおよび回転速度ωを推定する。
 磁極位置θおよび回転速度ωを推定する方法として様々な方法が検討されている。交流電動機1の回転速度の全域のうち中高速の領域では交流電動機1の速度起電力の情報を利用して磁極位置θを求めるのが一般的である。実施の形態1にかかる推定部15は、アークタンジェント法を用いて磁極位置θおよび回転速度ωを推定する。なお、磁極位置θおよび回転速度ωの推定方法は、アークタンジェント法に限定されない。例えば、磁極位置θおよび回転速度ωの推定方法には、適応オブザーバを用いる方法がある。適応オブザーバについては、実施の形態2において説明する。
 アークタンジェント法は最もプリミティブな位置推定法であり、広く知られている。下記式(1)は固定子座標上における表面磁石型同期交流電動機の電圧方程式である。下記式(1)において、「p」は微分演算子、「R」は電機子抵抗、「L」は電機子インダクタンス、「vα,vβ」は固定子座標上の電圧、「Φ」は電機子鎖交磁束数、「iα,iβ」は固定子座標上の電流である。以下において、固定子座標を静止座標と記載する場合がある。
Figure JPOXMLDOC01-appb-M000001
 上記式(1)の右辺第2項は速度起電力を表している。なお、速度起電力の項は下記式(2)のように表現可能である。下記式(2)において、「eα」はα軸速度起電力、「eβ」はβ軸速度起電力、「p」は微分演算子、「φαr」は回転子α軸磁束、「φβr」は回転子β軸磁束、「Φ」は電機子鎖交磁束数である。
Figure JPOXMLDOC01-appb-M000002
 上記式(1)を見て分かるとおり、速度起電力に磁極位置情報である磁極位置θの情報が含まれるため、上記式(1)を整理して磁極位置θを演算する。まず、回転子磁束項を左辺に、それ以外の項を右辺にまとめると、下記式(3)が得られる。微分計算はノイズを増幅させるため、下記式(3)の両辺を積分して下記式(4)を得る。ここで、電圧センサなどに直流オフセットがある場合、純粋積分を使うと積分値が発散するため、下記式(4)の計算を行う際はローパスフィルタによる近似積分を使い、直流分は積分しないようにするのが慣例である。
Figure JPOXMLDOC01-appb-M000003
Figure JPOXMLDOC01-appb-M000004
 上記式(4)において、記号「^」は推定値を表す。上記式(4)を計算して、回転子磁束を求め、求めた回転子磁束を用いて下記式(5)に示すアークタンジェントの計算を行うことで、回転子の磁極位置θを推定することができる。推定された回転子の磁極位置θを用いれば、回転速度ωを計算可能なため、下記式(6)で推定速度ω^を計算することができる。但し、微分ノイズの影響を避けるため、この推定速度ω^を制御に使用する場合にはローパスフィルタをかけるのが一般的である。また、速度起電力の推定を行い、下記式(7)のように、速度起電力の振幅を電機子鎖交磁束数Φで除算することにより、推定速度ω^を計算することもできる。速度起電力は上記式(3)の右辺を計算することによって算出することができる。しかしながら、永久磁石の磁束は温度変化により変動するため、下記式(7)の計算方法では温度変化によって定常的な速度推定誤差が生じる。そのため、下記式(6)による速度推定方法による誤差は、下記式(7)に示す方法による誤差よりも少ない。また、アークタンジェント法は計算量が少ないというメリットがある。実施の形態1では下記式(6)による速度推定を行った場合について説明する。
Figure JPOXMLDOC01-appb-M000005
Figure JPOXMLDOC01-appb-M000006
Figure JPOXMLDOC01-appb-M000007
 まず、上述した電流検出時刻間の時刻差Tbが変動しないものとして、上記式(4)から式(6)で表されるアークタンジェント法の演算を離散化する場合について説明する。交流電動機1に印加する静止座標系の電圧から交流電動機1の電機子抵抗による電圧降下を差し引き、ローパスフィルタによって下記式(8)のように近似積分の演算を行うことができる。下記式(8)において、「T」は制御周期、「TLPF」は近似積分のローパスフィルタの時定数、「n」は今回のサンプリング点を表す値、「n-1」は前回のサンプリング点を表す値である。今回のサンプリング点は、現在のサンプリング点ということもできる。下記式(8)は、上記式(4)の右辺第1項の積分を離散化したものということができる。ここでは、ローパスフィルタの離散化は後退差分によって行っているが、前進差分または双一次変換を用いてローパスフィルタの離散化を行ってもよい。
Figure JPOXMLDOC01-appb-M000008
 また、静止座標系の回転子磁束は、下記式(9)のように離散化することができる。したがって、今回のサンプリング点における推定位置θ^(n)は、下記式(10)によって算出することができる。また、今回のサンプリング点における推定位置θ^(n)と前回のサンプリング点における推定位置θ^(n)を用いて、2つのサンプリング点間における推定速度ω^は、下記式(11)によって算出することができる。「θ^」は、「θ」の推定値であり、「ω^」は、「ω」の推定値である。
Figure JPOXMLDOC01-appb-M000009
Figure JPOXMLDOC01-appb-M000010
Figure JPOXMLDOC01-appb-M000011
 しかしながら、上述したように、電流検出時刻間の時刻差Tbは時々刻々と変動する。したがって、電流検出時刻間の時刻差Tbを用いて演算しなければ、推定位置θ^および推定速度ω^に誤差が生じる。
 そこで、推定部15は、電流検出時刻間の時刻差Tbを用いて、アークタンジェント法による演算を行う。電流検出時刻間の時刻差TbをT+ΔT(n)と表したとき、上記式(8)で示されるローパスフィルタの演算は、下記式(12)のように書き換えることができる。
Figure JPOXMLDOC01-appb-M000012
 また、静止座標系の回転子磁束は、上記式(12)で求めた値を上記式(9)に代入することによって算出することができる。したがって、今回のサンプリング点における推定位置θ^(n)は、上記式(10)によって算出することができる。また、2つのサンプリング点間における推定速度ω^は、電流検出時刻間の時刻差Tb=T+ΔT(n)を用いた下記式(13)により算出することができる。
Figure JPOXMLDOC01-appb-M000013
 推定部15は、サンプリング点毎に相電流判定部12によって判定された相電流i,i,iを静止座標系の電流iα,iβへ座標変換する。また、推定部15は、電圧指令v ,v ,v を静止座標系の電圧vα,vβへ座標変換する。推定部15は、静止座標系の電流iα,iβおよび電圧vα,vβと、時刻差算出部13によって算出された時刻差Tbとに基づいて、上記式(12)の演算を行う。推定部15は、上記式(12)によって算出された値と、静止座標系の電流iα,iβとを上記式(9)へ代入して、静止座標系の回転子磁束を算出する。推定部15は、上記式(9)によって算出された静止座標系の回転子磁束を上記式(10)に代入して、推定位置θ^を算出する。また、推定部15は、上記式(10)の演算によって算出した推定位置θ^と、時刻差算出部13によって算出された時刻差Tbとを上記式(13)に代入して、推定速度ω^を算出する。
 上記式(12),(13)で示されるΔT(n)の定義域は、図5からも分かるように、-0.5×T<ΔT(n)<0.5×Tである。したがって、上記式(12),(13)を用いた演算と、上記式(8),(11)を用いた演算とで誤差に大きな差が生じる。すなわち、ΔTを用いて推定処理を行うか否かによって推定誤差に大きな差が生じる。
 ここで、上記式(11)を用いた演算による速度推定について説明する。ここでは、キャリア波Scの周波数であるキャリア周波数が5kHzであり、相電流判定部12は、キャリア波Scの1周期に1回、3つの相電流i,i,iの値を判定する場合を例に挙げて説明する。制御周期Tは、キャリア周波数の逆数であり、200μsである。
 交流電動機1の回転速度ωが2500rad/sである場合、交流電動機1の磁極位置θは、200μsの間に0.5rad進む。しかしながら、電流検出時刻間の時刻差Tbは、100μsから300μsの範囲で変動する。したがって、例えば、時刻差Tbが300μsであった場合、磁極位置θは、300μsの間に0.75rad進む。この場合、上記式(11)のように0.75radを制御周期Tで除算すると、推定速度ω^は、3750radになり、回転速度ωの真値に対して1.5倍になってしまう。この推定誤差は、交流電動機1を制御する上で無視できないほど大きい。
 ΔT(n)の平均値は零であることから、推定速度ω^にローパスフィルタをかけることによって、ΔT(n)による誤差を除去することができるが、制御ループ内にローパスフィルタを挿入すると、制御応答が下がってしまう。そのため、上記式(11)の演算による推定処理では高速応答を実現することが難しい。
 推定部15は、上記式(13)の演算を行うため、0.75radは時刻差Tbで除算される。そのため、推定速度ω^は、2500radであり、回転速度ωの真値と一致する。このように、推定部15は、上記式(13)を用いた演算を行うことによって、上記式(11)を用いた演算を行う場合に比べ、速度推定の誤差を大幅に低減することができる。また、推定部15は、上記式(12)の演算を行うことで、上記式(8)の演算を行う場合に比べ、位置推定の誤差を大幅に低減することができる。
 以上のように、実施の形態1にかかる推定装置9は、母線電流検出部11と、相電流判定部12と、時刻差算出部13と、推定部15とを備える。母線電流検出部11は、電圧指令v ,v ,v に基づいて生成される複数のゲートパルス信号G,G,Gによって駆動される電圧形インバータ4の直流母線61に流れる電流である母線電流ibusの値を検出する。相電流判定部12は、母線電流検出部11によって検出される母線電流ibusの値と複数のゲートパルス信号G,G,Gの状態とに基づいて、電圧形インバータ4から交流電動機1へ供給される相電流i,i,iの値を判定する。時刻差算出部13は、相電流判定部12による前回の判定処理で用いられた母線電流ibusの検出時刻として特定される電流検出時刻と相電流判定部12による今回の判定処理で用いられた母線電流ibusの検出時刻として特定される電流検出時刻との時刻差Tbを算出する。相電流判定部12による前回の判定処理で用いられた母線電流ibusの検出時刻として特定される電流検出時刻は、第1の電流検出時刻の一例である。相電流判定部12による今回の判定処理で用いられた母線電流ibusの検出時刻として特定される電流検出時刻は、第2の電流検出時刻の一例である。推定部15は、電圧指令v ,v ,v と、相電流i,i,iの値と、時刻差Tbとに基づいて、交流電動機1の磁極位置θおよび回転速度ωのうち少なくとも一方を推定する。このように、推定装置9は、相電流判定部12の電流検出時刻間の時刻差Tbを用いて推定処理を行う。そのため、推定装置9は、キャリア波Scにおける連続する2つの半周期の各々で検出される母線電流ibusの値から補間処理により3相の電流の値を判定する場合に比べ、3つの相電流i,i,iの値を判定する処理の周期である判定周期を短くできる。また、推定装置9は、電流検出時刻間の時刻差Tbを磁極位置θおよび回転速度ωの推定に用いることから、電流検出時刻間の時刻差Tbの変動が大きい場合であっても、磁極位置θおよび回転速度ωの推定精度を高めることができる。したがって、交流電動機1が高速に回転している場合または零ベクトル変調が行われている場合において、磁極位置θおよび回転速度ωの推定精度を高めることができる。
 また、ゲートパルス信号G,G,Gは、電圧指令v ,v ,v とキャリア波Scとの比較に基づいて生成される。そして、電流検出時刻間の時刻差Tbは、キャリア波Scの半周期の期間よりも長く、キャリア波Scの1.5倍の周期よりも短い。これにより、推定装置9は、交流電動機1が高速に回転している場合または零ベクトル変調が行われている場合において、磁極位置θおよび回転速度ωの推定精度を高めることができる。
 また、推定装置9を備える駆動装置2は、電圧形インバータ4へ入力される複数のゲートパルス信号G,G,Gのうち少なくとも1つを時間的にずらすパルスシフト処理部22を備える。パルスシフト処理部22は、ゲートパルス信号G,G,Gのうち変化の間隔が予め設定された時間TA未満である2つのゲートパルス信号のうち少なくとも1つを時間的にずらして2つのゲートパルス信号の変化の間隔を予め設定された時間TA以上にする。これにより、電圧指令v ,v ,v の値にかかわらず、リンギングが収束するまでの待ち時間を確保することができる。
 また、駆動装置2は、電圧指令v ,v ,v に乱数値を加えることによって、電圧形インバータ4から出力される2種類の零電圧ベクトルの出力比率を不規則に変化させる零ベクトル変調部34を備える。これにより、零ベクトル変調によりキャリア騒音のスペクトルピークが分散され、キャリア騒音を低減することができる。零ベクトル変調を行った場合、交流電動機1が低速で回転しているときでも中間相の電圧指令の波形が急激に変化し、電流検出時刻間の時刻差Tbの変化が大きくなる。このような場合であっても、推定装置9は、磁極位置θおよび回転速度ωの推定精度を高めることができる。
実施の形態2.
 実施の形態2にかかる推定装置の推定部は、適応オブザーバを用いて、交流電動機の位置および速度を推定する点で、実施の形態1にかかる推定装置の推定部と異なる。以下においては、実施の形態1と同様の機能を有する構成要素については同一符号を付して説明を省略し、実施の形態1と異なる点を中心に説明する。
 図15は、本発明の実施の形態2にかかる推定装置を含む電動機システムの構成例の一例を示す図である。図15に示すように、実施の形態2にかかる電動機システム100Aは、制御装置5を有する駆動装置2に代えて、制御装置5Aを有する駆動装置2Aを備える。制御装置5Aは、推定部15を有する推定装置9に代えて、推定部15Aを有する推定装置9Aを備える。
 図16は、実施の形態2にかかる推定部の構成例を示す図である。推定部15Aは、適応オブザーバによって構成されており、交流電動機1の磁極位置θおよび回転速度ωを推定する。適応オブザーバは、交流電動機1の固定子磁束および回転子磁束を状態量とする状態方程式で規定されているため、適応磁束オブザーバとも呼ばれる。なお、状態量として拡張誘起電圧または電流などを採用して適応オブザーバを構成することもできる。
 図16に示す推定部15Aは、電圧ベクトルと、電流ベクトルと、電流検出時刻間の時刻差Tbとを用いて、推定速度ω^と、推定位置θ^とを算出し、算出した推定速度ω^と推定位置θ^とを出力する。電圧ベクトルは、電圧形インバータ4から交流電動機1へ出力される電圧のベクトルであり、電流ベクトルは、電圧形インバータ4から交流電動機1へ出力される電流のベクトルである。
 推定部15Aは、モデル偏差演算部51と、角速度推定器52と、一次角周波数演算器53と、積分器54とを備える。モデル偏差演算部51は、電圧ベクトルと、電流ベクトルと、一次角周波数ωと、推定速度ω^と、電流検出時刻間の時刻差Tbとに基づきモデル偏差εを演算する。角速度推定器52は、モデル偏差εと電流検出時刻間の時刻差Tbとに基づき推定速度ω^を演算する。一次角周波数演算器53は、推定磁束ベクトルと、推定電流ベクトルと、推定速度ω^とを用いて、一次角周波数ωを演算する。積分器54は、電流検出時刻間の時刻差Tbに基づいて一次角周波数ωを積分して推定位置θ^を出力する。
 モデル偏差演算部51は、電流推定器71と、減算器72と、偏差演算器73とを備える。電流推定器71は、電圧ベクトルと、電流ベクトルと、一次角周波数ωと、推定速度ω^と、電流検出時刻間の時刻差Tbとに基づき、推定磁束ベクトルと推定電流ベクトルとを算出し、算出した推定磁束ベクトルと推定電流ベクトルとを出力する。減算器72は、推定電流ベクトルから電流ベクトルを減算することによって電流偏差ベクトルを算出し、算出した電流偏差ベクトルを出力する。
 偏差演算器73は、減算器72からの電流偏差ベクトルを入力とし、推定磁束ベクトルの直交成分をスカラ量として抽出し、抽出したスカラ量をモデル偏差εとして出力する。推定磁束ベクトルの直交成分をスカラ量として抽出する手法としては、電流偏差ベクトルを回転する直交座標系上に座標変換する手法と、電流偏差ベクトルと推定磁束ベクトルとの外積値の大きさを演算する手法とが公知である。
 電流推定器71は、交流電動機1の状態方程式から電流と磁束を推定する。ここでは、交流電動機1は一般的な埋込型永久磁石同期電動機であると仮定するが、方程式が立式できれば、誘導電動機、表面型永久磁石同期電動機、巻線界磁式同期電動機、またはリラクタンス式同期電動機などのように他種の電動機であってもよい。すなわち、電流推定器71は、他種の電動機についても、同様の方法で電流推定を行うことができる。
 交流電動機1が永久磁石埋込型同期電動機の場合、状態方程式は下記式(14)および式(15)のように表される。但し、「L」はd軸のインダクタンス、「L」はq軸のインダクタンス、「i」はd軸電流、「i」はq軸電流、「φds」はd軸固定子磁束、「φqs」はq軸固定子磁束、「φdr」はd軸回転子磁束、記号「^」は推定値を表す。また、「R」は電機子抵抗、「ω」は一次角周波数、「v」はd軸電圧、「v」はq軸電圧、「h11」から「h32」はオブザーバゲインを表す。また、一次角周波数ωは下記式(16)のように与えられる。「h41」,「h42」はオブザーバゲインを表す。
Figure JPOXMLDOC01-appb-M000014
Figure JPOXMLDOC01-appb-M000015
Figure JPOXMLDOC01-appb-M000016
 上記式(14)および式(15)は通常の誘起電圧に基づく式であるが、上記式(14)および式(15)に変形を加えて拡張誘起電圧の形式で表現しても同様の計算ができる。なお、上記式(14)および式(15)はdq回転座標における数式であるが、上記式(14)および式(15)に変形を加えて静止αβ座標または3相座標といった他の座標系で表現しても同様の計算ができる。上記式(14)には推定速度ω^が含まれるため、推定速度ω^と実際の回転速度ωとが一致していない場合、電流推定に誤差が生じる。ここではモデル偏差εを下記式(17)のように定義し、推定部15Aはモデル偏差εが零になるように、角速度推定器52を用いて推定速度ω^の値を調整する。角速度推定器52は、例えば、比例積分制御器に積分器を直列接続して構成される。
Figure JPOXMLDOC01-appb-M000017
 一次角周波数演算器53は、上記式(16)を基に、推定磁束ベクトルと、推定電流ベクトルと、推定速度ω^とから一次角周波数ωを演算する。積分器54は、一次角周波数ωを積分することにより磁極位置θを推定する。適応オブザーバは、鎖交磁束数の変動にロバストであり、定常的な速度推定誤差が発生しない点で優れており、高性能な速度推定法である。
 上記式(14)および式(15)の状態方程式は、電流検出時刻間の時刻差Tbが変動することから、離散化する際には注意を要する。オイラー法による離散化を行うと、下記式(18)から式(21)が得られる。但し、「T」は制御周期、「n」は今回のサンプリング点を表す値、「n-1」は前回のサンプリング点を表す値、「tavg」は、電流検出サンプル点からキャリア波Scの谷までの時間を表す。キャリア波Scの谷は、例えば、キャリア波Scが極小値になる位置を示し、図5に示す時刻t7,t16で各々現われる。電流検出サンプル点は、例えば、上述した電流検出時刻である。また、一次角周波数ωは下記式(22)のように与えられる。
Figure JPOXMLDOC01-appb-M000018
Figure JPOXMLDOC01-appb-M000019
Figure JPOXMLDOC01-appb-M000020
Figure JPOXMLDOC01-appb-M000021
Figure JPOXMLDOC01-appb-M000022
 推定部15Aのモデル偏差演算部51は、推定位置θ^を用いて座標変換を行う。モデル偏差演算部51は、積分器54によって得られる推定位置θ^を用いて、3つの相電流i,i,iを、回転する直交座標系であるdq座標系の電流へ座標変換することで、d軸電流iおよびq軸電流iを求める。また、モデル偏差演算部51は、積分器54によって得られる推定位置θ^を用いて、3相の電圧指令v ,v ,v を、回転する直交座標系であるdq座標系の電圧へ座標変換することで、d軸電圧vおよびq軸電圧vを求める。なお、d軸電圧vおよびq軸電圧vを得るために座標変換に用いる推定位置θ^は、例えば、電流検出時刻間の中央の時刻のものを用いるとよい。なお、上記式(19)および式(20)の状態方程式は、固定子磁束と回転子磁束を状態量として用いているが、電流または拡張誘起電圧を状態量としてもよい。
 上記式(19)には、推定速度ω^が含まれているため、推定速度ω^と実際の回転速度ωとが一致していない場合、電流推定に誤差が生じる。ここでは、モデル偏差εを下記式(23)のように定義し、推定部15Aは、モデル偏差εが零になるように、角速度推定器52を用いて推定速度ω^の値を調整する。角速度推定器52は、下記式(24)から式(26)の演算によって、推定速度ω^を算出する。下記式(25)において、「x」は、比例積分制御の状態量であり、「n+1」は、次回のサンプリング点である。
Figure JPOXMLDOC01-appb-M000023
Figure JPOXMLDOC01-appb-M000024
Figure JPOXMLDOC01-appb-M000025
Figure JPOXMLDOC01-appb-M000026
 上記式(25)では、ΔT(n)が用いられ、上記式(26)では、ΔT(n-1)が用いられている。そのため、推定部15Aは、電流検出時刻間の時刻差Tbが変動する場合においても、精度よく推定速度ω^を算出することができる。
 一次角周波数演算器53は、推定磁束ベクトルと、推定電流ベクトルと、推定速度ω^とを上記式(22)に代入することによって、一次角周波数ωを算出する。積分器54は、下記式(27)を用いて、一次角周波数ωを積分することによって推定位置θ^を算出する。なお、Tb=T+ΔT(n)である。下記式(27)では、ΔT(n)が用いられているため、推定部15Aは、電流検出時刻間の時刻差Tbの変動がある場合においても、精度よく推定位置θ^を算出することができる。
Figure JPOXMLDOC01-appb-M000027
 以上のように、実施の形態2にかかる推定部15Aは、電流推定器71と、減算器72と、偏差演算器73と、角速度推定器52とを備える。電流推定器71は、状態量推定部の一例であり、減算器72および偏差演算器73は、誤差算出部の一例であり、角速度推定器52は、速度推定部の一例である。電流推定器71は、交流電動機1の磁束ベクトルと電流ベクトルとを推定し、推定した結果である推定磁束ベクトルと推定電流ベクトルとを出力する。磁束ベクトルは、状態量の一例である。偏差演算器73は、モデル偏差εを算出する。モデル偏差εは、状態量の誤差の一例である。角速度推定器52は、モデル偏差εと時刻差Tbとに基づいて、推定速度ω^を算出する。また、電流推定器71は、電圧ベクトルと、電流ベクトルと、一次角周波数ωと、推定速度ω^と、電流検出時刻間の時刻差Tbとに基づいて、推定磁束ベクトルと推定電流ベクトルを算出する。推定部15Aでは、定常的な速度推定誤差が発生しないことから、速度推定を精度よく行うことができる。
実施の形態3.
 実施の形態3では、電流検出時刻間で電圧形インバータから交流電動機に印加される電圧を用いて交流電動機の位置および速度を推定する点で、実施の形態1,2にかかる駆動装置と異なる。以下においては、実施の形態1と同様の機能を有する構成要素については同一符号を付して説明を省略し、実施の形態1,2の制御装置5,5Aと異なる点を中心に説明する。
 図17は、本発明の実施の形態3にかかる推定装置を含む電動機システムの構成例の一例を示す図である。図17に示すように、実施の形態3にかかる電動機システム100Bは、制御装置5を有する駆動装置2または制御装置5Aを有する駆動装置2Aに代えて、制御装置5Bを有する駆動装置2Bを備える。制御装置5Bは、推定装置9または推定装置9Aに代えて、推定装置9Bを備える。
 推定装置9Bは、母線電流検出部11と、相電流判定部12と、時刻差算出部13と、検出時刻間電圧演算部14と、推定部15Bとを備える。検出時刻間電圧演算部14は、ゲートパルス信号G,G,Gと、電流検出時刻間の時刻差Tbとに基づいて、電流検出時刻間の3つの相電圧vub,vvb,vwbの値を算出する。検出時刻間電圧演算部14によって算出される相電圧vub,vvb,vwbの値は、電流検出時刻間におけるu相、v相、およびw相の電圧である3つの相電圧の瞬時値である。推定部15Bは、電圧指令v ,v ,v に代えて、電流検出時刻間の3つの相電圧vub,vvb,vwbの値を用いて、推定速度ω^および推定位置θ^を算出する点で、推定部15または推定部15Aと異なる。相電圧vubは、u相の電圧であり、相電圧vvbは、v相の電圧であり、相電圧vwbは、w相の電圧である。以下において、相電圧vubをu相電圧vubと記載し、相電圧vvbをv相電圧vvbと記載し、相電圧vwbをw相電圧vwbと記載する場合がある。また、相電圧vub,vvb,vwbの値を単に相電圧vub,vvb,vwbと記載する場合がある。
 検出時刻間電圧演算部14は、電流検出時刻間における各ゲートパルス信号G,G,Gのデューティ比に基づいて、電流検出時刻間の3つの相電圧vub,vvb,vwbの値を算出する。具体的には、検出時刻間電圧演算部14は、電流検出時刻間においてゲートパルス信号Gがオンになっているオン時間の比に基づいて、電流検出時刻間のu相の電圧vの値を算出する。同様に、検出時刻間電圧演算部14は、電流検出時刻間においてゲートパルス信号Gがオンになっている時間であるオン時間に基づいて、電流検出時刻間のv相の電圧vの値を算出する。また、検出時刻間電圧演算部14は、電流検出時刻間においてゲートパルス信号Gがオンになっている時間であるオン時間に基づいて、電流検出時刻間のw相の電圧vの値を算出する。
 図18は、実施の形態3にかかる検出時刻間電圧演算部の構成例を示す図である。図18に示すように、検出時刻間電圧演算部14は、オン時間算出部41と、電圧算出部42とを備える。オン時間算出部41は、電流検出時刻間において複数のゲートパルス信号G,G,Gの各々がオンになっている時間であるオン時間を算出する。電圧算出部42は、母線電圧Vdcの値と、オン時間算出部41によって算出されたオン時間と、時刻差算出部13によって算出された時刻差Tbとに基づいて、3つの相電圧vub,vvb,vwbの値を算出する。母線電圧Vdcは、低電位側の直流母線61と高電位側の直流母線61との間の電圧である母線電圧であり、母線電圧Vdcの値は、母線電圧Vdcの瞬時値である。
 図19は、実施の形態3にかかる検出時刻間電圧演算部による電流検出時刻間の3つの相電圧の算出方法を説明するための図である。図19に示す例では、最初の相電流判定周期で2相目の電流の値を判定するために用いられた母線電流ibusの検出時刻である時刻t5が最初の電流検出時刻として特定される。また、次の相電流判定周期で同じ2相目の電流の値を判定するために用いられた母線電流ibusの検出時刻である時刻t15が次の電流検出時刻として特定される。したがって、図19に示す例では、電流検出時刻間の時刻差Tbは、Tb=t15-t5で求められる。なお、図19では、各ゲートパルス信号G,G,Gは、パルスシフト処理部22によるパルスシフト処理が行われていない。
 検出時刻間電圧演算部14は、下記式(28)の演算によって、電流検出時刻間のu相電圧vubの値を算出することができる。下記式(28)において、「tuon1」は電流検出時刻間の前半でu相のゲートパルス信号Gのオン時間であり、「tuon2」は電流検出時刻間の後半でu相のゲートパルス信号Gのオン時間である。また、「Vdc」は母線電圧であり、不図示の母線電圧検出手段により検出される。なお、母線電圧が一定である場合、母線電圧検出手段を用いなくてもよい。
Figure JPOXMLDOC01-appb-M000028
 同様に、検出時刻間電圧演算部14は、下記式(29)および式(30)の演算によって、電流検出時刻間のv相電圧vvbの値および電流検出時刻間のw相電圧vwbの値を算出することができる。「tvon1」は電流検出時刻間の前半でv相のゲートパルス信号Gがオンとなっている時間であり、「tvon2」は電流検出時刻間の後半でv相のゲートパルス信号Gがオンとなっている時間である。また、「twon1」は電流検出時刻間の前半でw相のゲートパルス信号Gがオンとなっている時間であり、「twon2」は電流検出時刻間の後半でw相のゲートパルス信号Gがオンとなっている時間である。
Figure JPOXMLDOC01-appb-M000029
Figure JPOXMLDOC01-appb-M000030
 相電流判定部12で用いられる母線電流ibusの検出時刻は変動することから、電流検出時刻間の時刻差Tbは、ほとんどの場合、キャリア波Scの周期と一致しない。そのため、電流検出時刻間のu相電圧vubの値は、ほとんどの場合、電圧指令v と一致しない。同様の理由で、v相電圧vvbの値も、ほとんどの場合、電圧指令v と一致せず、w相電圧vwbの値も、ほとんどの場合、電圧指令v と一致しない。したがって、電圧指令v ,v ,v を速度推定および位置推定に用いた場合、電圧誤差によって速度推定および位置推定に誤差が生じる。
 相電流判定部12による3つの相電流i,i,iの判定処理がキャリア波Scの1周期毎に1度行われる場合、電流検出時刻間の時刻差Tbは、キャリア周期Tcに対し、0.5×Tc<Tb<1.5×Tcの範囲で変動する。この場合、電圧誤差の変動が大きく、速度推定および位置推定に生じる誤差の変動も大きくなる。
 推定部15Bは、検出時刻間電圧演算部14によって算出された電流検出時刻間の3つの相電圧vub,vvb,vwbの値を用いて速度推定および位置推定を行う。そのため、推定部15Bは、推定精度が高い速度推定および位置推定を行うことができる。
 推定部15Bは、検出時刻間電圧演算部14によって算出された電流検出時刻間の3つの相電圧vub,vvb,vwbを静止座標系の電圧vα,vβへ座標変換する。すなわち、推定部15Bは、電圧指令v ,v ,v に代えて電流検出時刻間の3つの相電圧vub,vvb,vwbを座標変換して得られる静止座標系の電圧vα,vβを用いる点で、推定部15,15Aと異なる。
 推定部15Bが適応オブザーバのように状態量の推定誤差から速度を算出するアルゴリズムの場合、推定部15Bは、状態量の推定値と真値の差分を算出して誤差を測定する。図16に示す適応オブザーバでは、減算器72は推定電流ベクトルと電流ベクトルとの差分を算出する。推定電流ベクトルと電流ベクトルとの差分を算出する処理において、推定電流ベクトルと電流ベクトルとに時刻のずれがあると、推定電流ベクトルと電流ベクトルとの差分の算出結果を精度よく行えない場合がある。推定部15Bの電流推定器71は、上述したように、電流検出時刻間の3つの相電圧vub,vvb,vwbの値を用いることで、電流検出時刻の推定電流ベクトルを精度よく求めることができる。
 ここで、検出時刻間電圧演算部14による電流検出時刻間の3つの相電圧vub,vvb,vwbの値の算出方法の一例について具体的に説明する。まず、電流検出時刻間における各ゲートパルス信号G,G,Gのオン時間の検出方法の一例について説明する。
 検出時刻間電圧演算部14のオン時間算出部41は、電流検出時刻間の期間を3つ以上の区間に分割する。オン時間算出部41は、分割した3つ以上の区間の各々においてゲートパルス信号G,G,Gがオンになっている時間であるオン時間をゲートパルス信号G,G,G毎に求める。オン時間算出部41は、分割した3つ以上の区間の各々のオン時間をゲートパルス信号G,G,G毎に合算することで、電流検出時刻間におけるオン時間をゲートパルス信号G,G,G毎に算出する。これにより、オン時間算出部41は、電流検出時刻間における各ゲートパルス信号G,G,Gのオン時間を精度よく算出することができる。
 以下、電流検出時刻間におけるゲートパルス信号G,G,Gのオン時間のうちゲートパルス信号Gのオン時間を例に挙げて説明する。検出時刻間電圧演算部14のオン時間算出部41は、オン時間tuon1とオン時間tuon2とを各々算出し、オン時間tuon1とオン時間tuon2とを合算することで、ゲートパルス信号Gのオン時間を求める。
 まず、オン時間tuon1の算出方法について説明する。ここでは、キャリア波Scの頂点の時刻から電流検出時刻への相対時間をTdet1,Tdet2と定義する。図19に示す例では、キャリア波Scの頂点は、キャリア波Scの谷である。キャリア波Scの谷の時刻は、時刻t7,t16である。
 相対時間Tdet1は、電流検出時刻である時刻t5とキャリア波Scが頂点になる時刻t7までの時間である。相対時間Tdet2は、電流検出時刻である時刻t15とキャリア波Scが頂点になる時刻t16までの時間である。相電流判定部12は、キャリア波Scが頂点になるタイミングと電流検出時刻との相対的な時刻差を示す相対時間Tdet1,Tdet2を算出することができる。相電流判定部12は、算出した相対時間Tdet1,Tdet2を検出時刻間電圧演算部14へ通知する。
 図19では、上述したパルスシフト処理を行っていない例を示しているが、ゲートパルス信号Gに対するパルスシフト処理が行われた場合、ゲートパルス信号Gのオンとオフとの切り替わりのタイミングが変わる。オン時間tuon1は、キャリア波Scの下り半周期と上り半周期に跨っているため、パルスシフト処理が行われた場合、オン時間tuon1のうちキャリア波Scの下り半周期の時間幅と上り半周期の時間幅とが変わる。そのため、オン時間算出部41は、オン時間tuon1を下り半周期と上り半周期に分けて算出する。
 ここで、時刻t1から時刻t11までのキャリア周期を第1のキャリア周期とし、時刻t11から時刻t17でのキャリア周期を第2のキャリア周期とする。第1のキャリア周期の下り半周期において、2相目の電流であるw相の電流iの値を判定するために用いられる母線電流ibusの値が検出される時刻t5の時点で、u相のゲートパルス信号Gはすでにオフからオンに切り替わっている。したがって、オン時間tuon1のうち下り半周期のオン時間は相対時間Tdet1と同じである。この場合、オン時間算出部41は、オン時間tuon1のうち下り半周期のオン時間が相対時間Tdet1と同じ時間であると判定する。
 仮に、時刻t5においてu相のゲートパルス信号Gがオフからオンへ切り替わっていない場合、オン時間tuon1のうち下り半周期のオン時間は、u相のゲートパルス信号Gから求めることができる。例えば、オン時間算出部41は、u相のゲートパルス信号Gがオフからオンへ切り替わる時刻からキャリア波Scの頂点の時刻t7までの時間を、オン時間tuon1のうち下り半周期のオン時間として算出することができる。
 オン時間算出部41は、オン時間tuon1のうち上り半周期のオン時間を、u相のゲートパルス信号Gから算出する。オン時間算出部41は、キャリア波Scの頂点の時刻t7からu相のゲートパルス信号Gがオンからオフへ切り替わる時刻までの時間を、オン時間tuon1のうち上り半周期のオン時間として算出することができる。
 次に、オン時間tuon2の算出方法について説明する。オン時間算出部41は、第2のキャリア周期の下り半周期におけるu相のゲートパルス信号Gのオン時間を算出し、算出した時間から相対時間Tdet2を減算することによって、オン時間tuon2を得ることができる。例えば、オン時間算出部41は、第2のキャリア周期において、u相のゲートパルス信号Gがオフからオンへ切り替わる時刻t12からキャリア波Scの頂点の時刻t16までの時間を算出する。オン時間算出部41は、時刻t12から時刻t16までの時間から相対時間Tdet2を減算することによって、オン時間tuon2を得る。但し、オン時間tuon2の最小値は零であるため、オン時間算出部41は、計算結果が負になる場合は強制的にオン時間tuon2を零にする。
 相電流判定部12は、相対時間Tdet1,Tdet2に基づいて、電流検出時刻間の時刻差Tbを算出することもできる。電流検出時刻の時刻差Tbの変動分ΔT(t)は、相対時間Tdet1と相対時間Tdet2との差である。また、電流検出時刻の時刻差Tbは、キャリア周期Tcに変動分ΔT(t)を加えたものである。相電流判定部12は、下記式(31)の演算によって、電流検出時刻間の時刻差Tbを算出することができる。
Figure JPOXMLDOC01-appb-M000031
 このように、実施の形態3にかかる制御装置5Bは、電流検出時刻の絶対時間を検出することなく、相対時間Tdet1,Tdet2に基づいて、電流検出時刻間の時刻差Tbと電流検出時刻間の3つの相電圧vub,vvb,vwbの値とを算出することができる。
 なお、オン時間算出部41は、ゲートパルス信号G,G,Gに代えて、3相の電圧指令v ,v ,v に基づいて、ゲートパルス信号G,G,Gのオン時間を算出することもできる。この場合、オン時間算出部41は、3相の電圧指令v ,v ,v の大きさからゲートパルス信号G,G,Gがオフからオンへ切り替わる時刻を算出する。また、中間相のゲートパルス信号Gがオフからオンへ切り替わってから時間TAが経過するまでの時間を算出する。そして、オン時間算出部41は、ゲートパルス信号G,G,Gがオフからオンへ切り替わる時刻と、中間相のゲートパルス信号Gがオフからオンへ切り替わってから時間TAが経過するまでの時間とに基づいて、ゲートパルス信号G,G,Gのオン時間を算出することができる。
 このように、検出時刻間電圧演算部14のオン時間算出部41は、電流検出時刻間の期間を3つの区間に分割し、分割した3つの区間の各々においてゲートパルス信号Gのオン時間を求める。オン時間算出部41は、分割した3つの区間の各々のオン時間を合算することで電流検出時刻間におけるゲートパルス信号Gのオン時間を算出する。上述した例では、3つの区間は、最初に電流検出が行われる下り半周期と、電流検出が行われない上り半周期と、次に電流検出が行われる下り半周期である。検出時刻間電圧演算部14のオン時間算出部41は、電流検出時刻間におけるゲートパルス信号Gのオン時間を電流検出時刻間の時刻差Tbで除算することで、電流検出時刻間におけるゲートパルス信号Gのデューティ比を算出する。検出時刻間電圧演算部14の電圧算出部42は、電流検出時刻間におけるゲートパルス信号Gのデューティ比から電流検出時刻間における電圧を算出する。
 検出時刻間電圧演算部14は、相電流判定部12による相電流の検出が行われる周期がキャリア波Scの周期の1.5倍以上の場合でも、電流検出時刻間を4つ以上の区間に分割し、場合分けして計算することで、電流検出時刻間の電圧を算出することができる。
 また、上述した例では、検出時刻間電圧演算部14は、ゲートパルス信号G,G,Gがオンになっている期間を算出するが、ゲートパルス信号G,G,Gがオフになっている期間を算出することもできる。キャリア波Scの下り半周期および上り半周期の長さは既知である。したがって、電流検出時刻間においてゲートパルス信号G,G,Gがオフになっている期間を算出することで、間接的に、電流検出時刻間においてゲートパルス信号G,G,Gがオンになっている時間を求めることができる。
 上述したように、電流検出時刻間の時刻差Tbは、ほとんどの場合、キャリア周期Tcと一致せず、また、交流電動機1の回転が高速であるほど、変動が大きい。また、図8および図9を用いて説明したように、中間相の波形は3相の電圧指令v ,v ,v の基本周波数に対して3n次の周波数成分を含んでいる。nは1以上の整数である。そのため、交流電動機1が高速で回転している場合、電流検出時刻間の時刻差Tbとキャリア周期Tcとのズレは、電気角3fの周期で顕著になる。電気角3fは、交流電動機1の電気角周波数の3倍の周波数である。
 図20は、実施の形態3にかかる電流検出時刻間の時刻差を固定値とした場合の速度推定の結果の一例を示す図である。図21は、図20に示す速度推定の結果をFFT解析した結果を示す図である。図20および図21では、脈動負荷の周波数は交流電動機1の回転軸の回転周波数と等しく、交流電動機1の極対数は3である。図20において、縦軸は交流電動機1の回転速度ωであり、横軸は時刻である。また、図21において、縦軸は交流電動機1の回転速度ωであり、横軸は周波数である。以下、交流電動機1の相電圧または相電流の周波数である電気的な周波数を「f」と記載し、交流電動機1の回転軸の回転周波数である機械的な周波数を「f」と記載する。なお、ここで、回転速度ωと電気的な周波数fとの関係は、ω=2π・fで表される。πは円周率である。また、交流電動機1の極対数を「P」と表したとき、電気的な周波数fと機械的な周波数fとの関係は、f=P・fで表される。
 図20に示すように、電流検出時刻間の時刻差Tbを固定値とした場合、速度真値に比べて推定速度の波形は、高周波の振動成分が大きい。図21に示すように、図20に示す推定速度をFFT解析すると、振動の主要な成分は交流電動機1の電気的な周波数fの3倍の周波数成分である。1シャント電流検出を用いた場合、電気的な周波数fの3倍の周波数で電流検出時刻が変化するが、電流検出時刻間の時刻差Tbの変動を無視しているために、図20に示すように、推定速度には、大きな誤差が生じている。また、図20に示すように、脈動負荷の影響によって速度真値は、機械的な周波数fの1倍の周波数で振動しているが、推定速度の波形ではこの振動を覆い隠すように機械的な周波数fの1.4倍の周波数にピークがある。そのため、機械的な周波数fの1倍の脈動がほとんど見えなくなっている。
 上述したように交流電動機1の極対数は3であるが、電気3f成分および機械9f成分の側波帯として、機械7.6f成分と機械10.4f成分のピークが立っている。電気3f成分は、電気的な周波数fの3倍の周波数成分であり、機械9f成分は、機械的な周波数fの9倍の周波数成分である。また、機械7.6f成分は、機械的な周波数fの7.6倍の周波数成分であり、機械10.4f成分は、機械的な周波数fの10.4倍の周波数成分である。
 このように、電流検出時刻間の時刻差Tbを固定値とした場合の速度推定では、1シャント電流検出方式による電流検出時刻の変化と電流検出時刻の変化に付随する電圧誤差とによって速度真値に含まれない成分の速度推定誤差が生じる。
 図22は、実施の形態3にかかる推定部による速度推定の結果の一例を示す図である。図23は、図22に示す速度推定の結果をFFT解析した結果を示す図である。図22および図23では、脈動負荷の周波数は機械的な周波数fと等しく、交流電動機1の極対数は3である。図22において、縦軸は交流電動機1の回転速度ωであり、横軸は時刻である。また、図23において、縦軸は交流電動機1の回転速度ωであり、横軸は周波数である。
 実施の形態3にかかる推定部15Bは、相電流判定部12によって特定される電流検出時刻間の時刻差Tbと、検出時刻間電圧演算部14によって算出される電流検出時刻間の相電圧vub,vvb,vwbの値とに基づいて、速度推定を行う。これにより、推定部15Bは、図22および図23に示すように、速度推定の誤差を大幅に低減することができる。同様に、実施の形態3にかかる推定部15Bは、相電流判定部12によって検出される電流検出時刻間の時刻差Tbと、検出時刻間電圧演算部14によって特定される電流検出時刻間の相電圧vub,vvb,vwbの値とに基づいて、位置推定を行う。これにより、推定部15Bは、位置推定の誤差を大幅に低減することができる。
 以上のように、実施の形態3にかかる推定装置9Bは、電流検出時刻間で電圧形インバータ4から交流電動機1に印加される相電圧vub,vvb,vwbの値を電圧指令v ,v ,v に基づいて算出する検出時刻間電圧演算部14を備える。検出時刻間電圧演算部14は、相電圧算出部の一例である。推定部15Bは、相電流i,i,iの値と、時刻差Tbと、検出時刻間電圧演算部14によって算出される相電圧vub,vvb,vwbの値とに基づいて、推定位置θ^と推定速度ω^とを算出する。
 また、検出時刻間電圧演算部14は、オン時間算出部41と、電圧算出部42とを備える。オン時間算出部41は、電流検出時刻間において複数のゲートパルス信号G,G,Gの各々がオンになっている時間であるオン時間を算出する。電圧算出部42は、母線電圧Vdcの値と、オン時間算出部41によって算出されたオン時間と、時刻差算出部13によって算出された時刻差Tbとに基づいて、相電圧vub,vvb,vwbの値を算出する。これにより、検出時刻間電圧演算部14は、ゲートパルス信号G,G,Gから電流検出時刻間における相電圧vub,vvb,vwbの値を算出することができる。
 また、オン時間算出部41は、電流検出時刻間の期間を3つ以上の区間に分割し、3つ以上の区間の各々において複数のゲートパルス信号G,G,Gの各々がオンになっている時間であるオン時間を求める。オン時間算出部41は、3つ以上の区間の各々のオン時間をゲートパルス信号G毎に合算することで、電流検出時刻間における各ゲートパルス信号Gのオン時間を算出する。これにより、検出時刻間電圧演算部14は、電流検出時刻間における各ゲートパルス信号G,G,Gのオン時間を精度よく算出することができる。
実施の形態4.
 実施の形態4では、交流電動機の駆動装置のうち上述した推定装置および電圧形インバータを除く部分の構成について具体的に説明する。以下においては、実施の形態3と同様の機能を有する構成要素については同一符号を付して説明を省略し、実施の形態3の電動機システム100Bと異なる点を中心に説明する。
 図24は、本発明の実施の形態4にかかる電動機システムの構成例を示す図である。図24に示すように、電動機システム100Cは、交流電動機1と、駆動装置2Cとを備える。駆動装置2Cは、電圧形インバータ4と、制御装置5Cを備える。制御装置5Cは、推定装置9Bと、速度制御部30と、座標変換部31,33と、dq軸電流制御部32と、零ベクトル変調部34と、ゲートパルス生成部35とを備える。なお、制御装置5Cは、推定装置9Bに代えて、推定装置9または推定装置9Aを備える構成であってもよい。なお、制御装置5Cは、dq座標系に基づくセンサレスベクトル制御であるとするが、他の制御法や座標系を用いても構わない。
 速度制御部30は、推定装置9Bによって算出された推定速度ω^が速度指令ωと一致するように、d軸電流指令i およびq軸電流指令i を生成する。速度制御部30は、PID(Proportional Integral Differential)制御器などのフィードバック制御器で構成されるが、フィードバック制御器以外の制御器で構成されてもよい。また、速度制御部30は、PID制御器と並列にフィードフォワード制御器を有する構成であってもよく、これにより、制御応答を早めることができる。また、速度制御部30は、繰り返し制御などの振動抑制制御器をPID制御器と併用して特定周波数の外乱応答を高める構成であってもよい。
 q軸電流指令i は、推定速度ω^と速度指令ωとの偏差がゼロに収束するよう決定されるが、d軸電流指令i の決定方法は、交流電動機1の種類によって種々の変更が可能である。例えば、交流電動機1が表面型永久磁石同期電動機である場合、速度制御部30は、電圧が飽和していない状況下ではd軸電流指令i をゼロにする制御方式である「i=0制御」を行うことができる。「i=0制御」によって、銅損を低減することができる。また、交流電動機1が埋込型永久磁石同期電動機である場合、リラクタンストルクを併用できるため、速度制御部30は、「最大トルク/電流制御(Maximum Torque Per Ampere control)」を行うこともできる。「最大トルク/電流制御」は、マイナス方向のd軸電流指令i を与える制御方式であり、銅損を低減することができる。また、速度制御部30は、より大きなマイナス方向のd軸電流指令i を与える制御方式である弱め磁束制御を行うこともできる。弱め磁束制御を行うことによって、電圧飽和時には、交流電動機1が表面型永久磁石同期電動機と埋込型永久磁石同期電動機のどちらの場合であっても、銅損の増加と引き換えに機械出力を増大させることができる。
 座標変換部31は、推定装置9Bによって算出された推定位置θ^に基づいて、座標変換によって相電流i,i,iをdq座標系の電流であるd軸電流iおよびq軸電流iへ変換する。dq軸電流制御部32は、d軸電流iがd軸電流指令i と一致するようにd軸電圧指令v を調整し、q軸電流iがq軸電流指令i と一致するようにq軸電圧指令v を調整する。dq軸電流制御部32は、例えば、d軸電流iをフィードバック制御する不図示のPID(Proportional Integral Differential)制御器と、q軸電流iをフィードバック制御する不図示のPID制御器と、dq軸の干渉成分をフィードフォワード補償する不図示の非干渉化制御とを備える。
 このように、dq軸電流制御部32は、dq座標系において、d軸電流iとq軸電流iとを個別に制御しており、これにより、良好な電流制御性能を実現している。なお、速度制御部30およびdq軸電流制御部32には微積分の演算が含まれるため、dq軸電流制御部32は、電流検出時刻間の時刻差Tbを使って、微積分の演算誤差を防いでいる。
 座標変換部33は、推定装置9Bによって算出された推定位置θ^に基づいて、d軸電圧指令v およびq軸電圧指令v を電圧指令v ,v ,v へ座標変換する。零ベクトル変調部34は、座標変換部33によって座標変換された電圧指令v ,v ,v に対して零ベクトル変調を行う。ゲートパルス生成部35は、電圧指令v ,v ,v に基づいて、ゲートパルス信号G,G,Gを生成し、生成したゲートパルス信号G,G,Gを電圧形インバータ4へ出力する。ゲートパルス生成部35は、上述したようにパルスシフト処理部22を有しているため、電圧指令v ,v ,v の値にかかわらず相電流判定部12によって3相電流の値を判定することができる。そのため、推定装置9Bは、位置推定および速度推定の精度をさらに向上させることができ、制御性能を一層高めることができる。
 推定装置9Bの推定部15Bは、上述したように、1シャント電流検出方式による電流検出時刻の時刻差Tbと電流検出時刻の時刻差Tbに付随して変化する電流検出時刻間の電圧とに基づいて速度推定および位置推定を行う。すなわち、推定部15Bは、相電流判定部12によって検出される電流検出時刻間の時刻差Tbと、検出時刻間電圧演算部14によって算出される電流検出時刻間の電圧とに基づいて、速度推定および位置推定を行う。これにより、推定装置9Bは、速度推定および位置推定の誤差を大幅に低減することができる。
 図25は、実施の形態4にかかる制御装置のハードウェア構成の一例を示す図である。図25に示すように、実施の形態4にかかる制御装置5Cは、母線電流検出部11と、処理回路80とを備える。処理回路80は、プロセッサ81と、メモリ82と、入力回路83と、AD(Analog to Digital)変換回路84と、DA(Digital to Analog)変換回路85とを備える。プロセッサ81は、CPU(Central Processing Unit)、DSP(Digital Signal Processer)、およびシステムLSI(Large Scale Integration)のうち1つ以上を含む。メモリ82は、RAM(Random Access Memory)、ROM(Read Only Memory)、フラッシュメモリ、EPROM(Erasable Programmable Read Only Memory)、およびEEPROM(登録商標)(Electrically Erasable Programmable Read Only Memory)のうち1つ以上を含む。また、メモリ82は、コンピュータが読み取り可能なプログラムが記録された記録媒体を含む。かかる記録媒体は、不揮発性または揮発性の半導体メモリ、磁気ディスク、フレキシブルメモリ、光ディスク、コンパクトディスク、およびDVD(Digital Versatile Disc)のうち1つ以上を含む。
 プロセッサ81は、メモリ82の記録媒体に記憶されたプログラムを読み出して実行することによって、相電流判定部12、時刻差算出部13、検出時刻間電圧演算部14、および推定部15Bの機能を実行する。入力回路83は、外部からの速度指令ωまたは3相の電圧指令v ,v ,v をプロセッサ81へ通知する。プロセッサ81は、外部からの速度指令ωまたは3相の電圧指令v ,v ,v に基づいて磁極位置θおよび回転速度ωを推定する処理を行う。AD変換回路84は、母線電流検出部11で検出される母線電流ibusの値をアナログ信号からデジタル信号へ変換する。プロセッサ81は、AD変換回路84からデジタル化された母線電流ibusの値を取得し、取得した母線電流ibusの値に基づく処理を行う。また、プロセッサ81は、3相の電圧指令v ,v ,v を変調してゲートパルス信号G,G,Gを生成し、生成したゲートパルス信号G,G,GをDA変換回路85へ出力する。DA変換回路85は、ゲートパルス信号G,G,Gをデジタル信号からアナログ信号へ変換し、アナログ信号へ変換したゲートパルス信号G,G,Gをゲートドライバ7へ出力する。なお、ゲートドライバ7に変調手段が設けられている場合、プロセッサ81は、3相の電圧指令v ,v ,v をゲートドライバ7へ送信することができる。この場合、3相の電圧指令v ,v ,v は、デジタル値であってもアナログ値であってもよい。
 プロセッサ81、メモリ82、入力回路83、AD変換回路84、およびDA変換回路85の間でデータの送受信が行われる。データは、デジタル信号で送受信されるが、アナログ信号で送受信されてもよい。また、データの送受信方法は、パラレル通信でもシリアル通信でもよい。また、アナログ信号とデジタル信号は不図示の変換器によって適宜変換しても構わない。
 なお、実施の形態1,2の推定装置9,9Aも、実施の形態4の推定装置9Bと同様に、図25に示すハードウェア構成によって実現できる。例えば、プロセッサ81は、メモリ82のうち記録媒体に記憶されたプログラムを読み出して実行することによって、相電流判定部12、時刻差算出部13、および推定部15の機能を実行することができる。また、プロセッサ81は、メモリ82のうち記録媒体に記憶されたプログラムを読み出して実行することによって、相電流判定部12、時刻差算出部13、および推定部15Aの機能を実行することができる。
 また、プロセッサ81は、推定装置9,9A,9Bの機能に加え、速度制御部30、座標変換部31,33、dq軸電流制御部32、零ベクトル変調部34、およびゲートパルス生成部35の機能の一部または全部を実行することもできる。この場合、プロセッサ81は、メモリ82に記憶されたプログラムを読み出して実行することによって、速度制御部30、座標変換部31,33、dq軸電流制御部32、零ベクトル変調部34、およびゲートパルス生成部35の機能を実行する。
 以上のように、実施の形態4にかかる駆動装置2Cでは、磁極位置θおよび回転速度ωの推定精度が高い推定装置9Bを用いている。そのため、駆動装置2Cでは、交流電動機1が高速に回転している場合または零ベクトル変調が行われている場合において、交流電動機1に対する制御の性能が向上する。
 以上の実施の形態に示した構成は、本発明の内容の一例を示すものであり、別の公知の技術と組み合わせることも可能であるし、本発明の要旨を逸脱しない範囲で、構成の一部を省略、変更することも可能である。
 1 交流電動機、2,2A,2B,2C 駆動装置、3 直流電源、4 電圧形インバータ、5,5A,5B,5C 制御装置、6 主回路、7 ゲートドライバ、9,9A,9B 推定装置、11 母線電流検出部、12 相電流判定部、13 時刻差算出部、14 検出時刻間電圧演算部、15,15A,15B 推定部、21 比較部、22 パルスシフト処理部、30 速度制御部、31,33 座標変換部、32 dq軸電流制御部、34 零ベクトル変調部、35 ゲートパルス生成部、41 オン時間算出部、42 電圧算出部、51 モデル偏差演算部、52 角速度推定器、53 一次角周波数演算器、54 積分器、61 直流母線、71 電流推定器、72 減算器、73 偏差演算器、100,100A,100B,100C 電動機システム、G,G,G,G,Gup,Gun,Gvp,Gvn,Gwp,Gwn ゲートパルス信号、ibus 母線電流、i,i,i 相電流、Sc キャリア波、Tb 時刻差、tuon1,tuon2 オン時間、vub,vvb,vwb 相電圧、v ,v ,v  電圧指令、θ 磁極位置、θ^ 推定位置、ω 回転速度、ω^ 推定速度。

Claims (9)

  1.  電圧指令に基づいて生成される複数のゲートパルス信号によって駆動される電圧形インバータの直流母線に流れる電流である母線電流の値を検出する母線電流検出部と、
     前記母線電流検出部によって検出される母線電流の値と、前記複数のゲートパルス信号の状態とに基づいて、前記電圧形インバータから交流電動機へ供給される複数の相電流の値を判定する判定処理を行う相電流判定部と、
     前記相電流判定部による前回の判定処理で用いられた前記母線電流の検出時刻として特定される第1の電流検出時刻と前記相電流判定部による今回の判定処理で用いられた前記母線電流の検出時刻として特定される第2の電流検出時刻との時刻差を算出する時刻差算出部と、
     前記相電流判定部によって判定された前記複数の相電流の値と、前記時刻差算出部によって算出された前記時刻差とに基づいて、前記交流電動機の位置および速度のうち少なくとも一方を推定する推定部と、を備える
     ことを特徴とする推定装置。
  2.  前記第1の電流検出時刻と前記第2の電流検出時刻との間に前記電圧形インバータから前記交流電動機に印加される相電圧の値を前記電圧指令に基づいて算出する相電圧算出部を備え、
     前記推定部は、
     前記複数の相電流の値と、前記時刻差と、前記相電圧算出部によって算出された前記相電圧の値とに基づいて、前記交流電動機の位置および速度のうち少なくとも一方を推定する
     ことを特徴とする請求項1に記載の推定装置。
  3.  前記相電圧算出部は、
     前記第1の電流検出時刻と前記第2の電流検出時刻との間において前記複数のゲートパルス信号の各々がオンになっている時間であるオン時間を算出するオン時間算出部と、
     前記直流母線の電圧の値と、前記オン時間算出部によって算出された前記オン時間と、前記時刻差算出部によって算出された前記時刻差とに基づいて、前記第1の電流検出時刻と前記第2の電流検出時刻との間の前記相電圧の値を算出する電圧算出部と、を備える
     ことを特徴とする請求項2に記載の推定装置。
  4.  前記オン時間算出部は、
     前記第1の電流検出時刻と前記第2の電流検出時刻との間の期間を3つ以上の区間に分割し、前記3つ以上の区間の各々において前記複数のゲートパルス信号の各々がオンになっている時間であるオン時間を求め、前記3つ以上の区間の各々の前記オン時間を合算することで前記第1の電流検出時刻と前記第2の電流検出時刻との間における前記オン時間を算出する
     ことを特徴とする請求項3に記載の推定装置。
  5.  前記ゲートパルス信号は、
     前記電圧指令とキャリア波との比較に基づいて生成され、
     前記時刻差は、前記キャリア波の半周期の期間よりも長く、前記キャリア波の1.5倍の周期よりも短い
     ことを特徴とする請求項1から4のいずれか1つに記載の推定装置。
  6.  前記推定部は、
     前記交流電動機の状態量を推定する状態量推定部と、
     前記状態量推定部によって推定された前記状態量の誤差を算出する誤差算出部と、
     前記誤差算出部によって算出された前記状態量の誤差に基づいて、前記交流電動機の速度を推定する速度推定部と、を備え、
     前記状態量推定部および前記速度推定部のうち少なくとも一方は、
     前記時刻差に基づく演算を行う
     ことを特徴とする請求項1から5のいずれか1つに記載の推定装置。
  7.  請求項1から6のいずれか1つに記載の推定装置を備える
     ことを特徴とする交流電動機の駆動装置。
  8.  前記電圧形インバータへ入力される前記複数のゲートパルス信号のうち少なくとも1つを時間的にずらすパルスシフト処理部を備え、
     前記パルスシフト処理部は、
     前記複数のゲートパルス信号のうち2つのゲートパルス信号の変化の間隔が予め設定された時間未満である場合に、前記2つのゲートパルス信号のうち少なくとも1つを時間的にずらして前記2つのゲートパルス信号の変化の間隔を前記予め設定された時間以上にする
     ことを特徴とする請求項7に記載の交流電動機の駆動装置。
  9.  前記電圧指令に乱数値を加えることによって、前記電圧形インバータから出力される2種類の零電圧ベクトルの出力比率を不規則に変化させる零ベクトル変調部を備える
     ことを特徴とする請求項7または8に記載の交流電動機の駆動装置。
PCT/JP2019/033182 2019-08-23 2019-08-23 推定装置および交流電動機の駆動装置 WO2021038670A1 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
US17/625,096 US11705841B2 (en) 2019-08-23 2019-08-23 Estimation device and AC motor drive device
JP2021541799A JP7130143B2 (ja) 2019-08-23 2019-08-23 推定装置および交流電動機の駆動装置
CN201980099465.8A CN114270695B (zh) 2019-08-23 2019-08-23 推测装置以及交流电动机的驱动装置
DE112019007652.3T DE112019007652T5 (de) 2019-08-23 2019-08-23 Schätzeinrichtung und wechselstrommotor-ansteuereinrichtung
PCT/JP2019/033182 WO2021038670A1 (ja) 2019-08-23 2019-08-23 推定装置および交流電動機の駆動装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2019/033182 WO2021038670A1 (ja) 2019-08-23 2019-08-23 推定装置および交流電動機の駆動装置

Publications (1)

Publication Number Publication Date
WO2021038670A1 true WO2021038670A1 (ja) 2021-03-04

Family

ID=74683919

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2019/033182 WO2021038670A1 (ja) 2019-08-23 2019-08-23 推定装置および交流電動機の駆動装置

Country Status (5)

Country Link
US (1) US11705841B2 (ja)
JP (1) JP7130143B2 (ja)
CN (1) CN114270695B (ja)
DE (1) DE112019007652T5 (ja)
WO (1) WO2021038670A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7280170B2 (ja) * 2019-11-22 2023-05-23 日立Astemo株式会社 モータ制御装置、モータ制御方法、ハイブリッドシステム、昇圧コンバータシステム、電動パワーステアリングシステム

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008263667A (ja) * 2007-04-10 2008-10-30 Mitsubishi Electric Corp インバータ制御装置並びに冷凍空調装置
JP2012070532A (ja) * 2010-09-24 2012-04-05 Mitsubishi Electric Corp 回転機の制御装置
JP2015139359A (ja) * 2014-01-24 2015-07-30 株式会社富士通ゼネラル モータ制御装置
JP2015198546A (ja) * 2014-04-03 2015-11-09 日本電産株式会社 モータ制御方法およびモータ制御装置

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3668870B2 (ja) 2001-08-09 2005-07-06 株式会社日立製作所 同期電動機駆動システム
JP4269686B2 (ja) * 2003-01-09 2009-05-27 株式会社明電舎 同期電動機の位置・速度・電圧振幅のセンサレス計測方法およびセンサレス可変速装置
JP5428325B2 (ja) * 2008-08-25 2014-02-26 株式会社ジェイテクト モータ制御装置および電動パワーステアリング装置
JP4746667B2 (ja) * 2008-11-26 2011-08-10 本田技研工業株式会社 電動機の相電流推定装置および電動機の磁極位置推定装置
JP6102768B2 (ja) * 2014-01-24 2017-03-29 株式会社富士通ゼネラル モータ制御装置
JP6091446B2 (ja) * 2014-02-10 2017-03-08 三菱電機株式会社 電動機制御装置
JP2016032342A (ja) * 2014-07-29 2016-03-07 株式会社明電舎 同期電動機の速度・位相推定方法
JP2016103886A (ja) * 2014-11-27 2016-06-02 株式会社富士通ゼネラル モータ制御装置
CN107681929B (zh) * 2017-10-30 2023-11-28 广东美的制冷设备有限公司 电机控制系统、变频空调器

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008263667A (ja) * 2007-04-10 2008-10-30 Mitsubishi Electric Corp インバータ制御装置並びに冷凍空調装置
JP2012070532A (ja) * 2010-09-24 2012-04-05 Mitsubishi Electric Corp 回転機の制御装置
JP2015139359A (ja) * 2014-01-24 2015-07-30 株式会社富士通ゼネラル モータ制御装置
JP2015198546A (ja) * 2014-04-03 2015-11-09 日本電産株式会社 モータ制御方法およびモータ制御装置

Also Published As

Publication number Publication date
CN114270695B (zh) 2024-05-14
US11705841B2 (en) 2023-07-18
JPWO2021038670A1 (ja) 2021-12-02
CN114270695A (zh) 2022-04-01
JP7130143B2 (ja) 2022-09-02
DE112019007652T5 (de) 2022-06-15
US20220294377A1 (en) 2022-09-15

Similar Documents

Publication Publication Date Title
JP5069306B2 (ja) 回転電機の制御装置
JP4631672B2 (ja) 磁極位置推定方法、モータ速度推定方法及びモータ制御装置
WO2017141513A1 (ja) 電力変換装置
JP5333256B2 (ja) 交流回転機の制御装置
KR100960043B1 (ko) 2상 영구자석 동기 전동기의 공간전압벡터 제어 장치
JP5637155B2 (ja) モータ制御装置及びモータ制御方法
WO2012066800A1 (ja) 電流検出装置及びモータ制御装置
JP6685452B1 (ja) 回転電機の制御装置
WO2021038670A1 (ja) 推定装置および交流電動機の駆動装置
JP2008206330A (ja) 同期電動機の磁極位置推定装置および磁極位置推定方法
JP7053955B2 (ja) 回転機の制御装置
JP2019057981A (ja) モータ制御用集積回路
JP5574925B2 (ja) 位置・速度センサレス制御装置
WO2017030055A1 (ja) 回転機の制御装置および制御方法
JP3939481B2 (ja) 交流モータの制御装置
JP2013066304A (ja) 電力変換装置
JP6753326B2 (ja) モータ制御装置
JP4566695B2 (ja) 巻線界磁式同期機の磁極位置検出装置
JP2021069187A (ja) 電力変換装置
JP2020141553A (ja) 埋込磁石同期電動機を用いたトルク加振制御システムおよびトルク加振制御方法
JP7362003B2 (ja) 回転機の制御装置
JP5590106B2 (ja) マトリクスコンバータ
WO2023026352A1 (ja) 回転機の制御装置
JP6848680B2 (ja) 同期電動機の制御装置
JP6756293B2 (ja) 回転機の制御装置および制御方法

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: 19943193

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2021541799

Country of ref document: JP

Kind code of ref document: A

122 Ep: pct application non-entry in european phase

Ref document number: 19943193

Country of ref document: EP

Kind code of ref document: A1