WO2016121751A1 - インバータ制御装置及びモータ駆動システム - Google Patents

インバータ制御装置及びモータ駆動システム Download PDF

Info

Publication number
WO2016121751A1
WO2016121751A1 PCT/JP2016/052161 JP2016052161W WO2016121751A1 WO 2016121751 A1 WO2016121751 A1 WO 2016121751A1 JP 2016052161 W JP2016052161 W JP 2016052161W WO 2016121751 A1 WO2016121751 A1 WO 2016121751A1
Authority
WO
WIPO (PCT)
Prior art keywords
command value
current command
current
calculation unit
value
Prior art date
Application number
PCT/JP2016/052161
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 EP16743354.9A priority Critical patent/EP3252941B1/en
Priority to CN201680003392.4A priority patent/CN107078674B/zh
Priority to JP2016535264A priority patent/JP6367332B2/ja
Publication of WO2016121751A1 publication Critical patent/WO2016121751A1/ja
Priority to US15/636,198 priority patent/US10158305B2/en

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00—Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/14—Electronic commutators
    • H02P6/16—Circuit arrangements for detecting position
    • H02P6/18—Circuit arrangements for detecting position without separate position detecting elements
    • H02P6/182—Circuit arrangements for detecting position without separate position detecting elements using back-emf in windings
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/22—Current control, e.g. using a current control loop
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/24—Vector control not involving the use of rotor position or rotor speed sensors
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00—Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/14—Electronic commutators
    • H02P6/16—Circuit arrangements for detecting position
    • H02P6/18—Circuit arrangements for detecting position without separate position detecting elements
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P25/00—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
    • H02P25/08—Reluctance motors

Definitions

  • Embodiments of the present invention relate to an inverter control device and a motor drive system.
  • a rotational phase angle sensorless control device for a permanent magnet synchronous motor (PMSM) or a synchronous reluctance motor (SynRM) a rotational phase angle estimation method using an induced voltage is used in a high speed range.
  • PMSM permanent magnet synchronous motor
  • SynRM and PMSM with a small magnet magnetic flux have a problem that the estimation accuracy of the rotational phase angle is deteriorated because the induced voltage is small in a low load state even in a high speed region.
  • an inverter control device and a motor drive system that can accurately estimate the rotational phase angle of a motor in a high-speed and low-load state.
  • An inverter control device includes an inverter main circuit, a current detector, a current command value calculation unit, a voltage command value calculation unit, and an estimation unit.
  • the inverter main circuit can be electrically connected to a predetermined rotational drive target.
  • the current detector detects a current value output from the inverter main circuit.
  • the current command value calculation unit calculates a current command value at which the output voltage output from the inverter main circuit is equal to or higher than a predetermined target value.
  • the voltage command value calculation unit calculates a voltage command value at which the current value is equal to the current command value.
  • the estimation unit calculates an estimated rotation phase angle of the rotation drive target based on the voltage command value and the current value.
  • the figure which shows the structure of the electric current command value calculation part of FIG. The figure explaining the calculation method of an electric current command value.
  • FIG. The figure explaining the calculation method of the electric current command value in the modification 2.
  • FIG. The figure which shows the structure of the modification 3 of an electric current command value calculation part.
  • FIG. 1 is a diagram illustrating a configuration of a system according to the present embodiment. As shown in FIG. 1, the system according to the present embodiment includes a motor 1 and an inverter control device 2.
  • the motor 1 is a rotational drive target of the control device 2 and is connected to the control device 2.
  • the SynRM 1 includes a stator and a rotor.
  • the stator has three excitation phases (U phase, V phase, and W phase).
  • the stator generates a magnetic field by a three-phase alternating current flowing in each excitation phase.
  • the rotor does not have a permanent magnet. The rotor rotates due to magnetic interaction with the magnetic field generated by the stator.
  • the inverter control device 2 (hereinafter referred to as “control device 2”) controls the rotational phase angle of the SynRM 1 without a sensor.
  • the control device 2 includes an inverter 21, a current detector 22, a coordinate conversion unit 23, a current command value calculation unit 24, a voltage command value calculation unit 25, and a coordinate A conversion unit 26, a PWM modulation unit 27, and a speed / rotation phase angle estimation unit 28 are provided.
  • the inverter 21 is a circuit including a switching element (transistor).
  • the inverter 21 converts electric power from a power source (not shown) into alternating current by switching ON / OFF of the switching element, and supplies the alternating current to the SynRM 1.
  • the inverter 21 receives a control signal for controlling ON / OFF of each switching element from the PWM modulator 22.
  • the current detector 22 detects a two-phase or three-phase current value among the three-phase alternating current flowing through the stator of the SynRM 1.
  • FIG. 1 shows a configuration for detecting current values iu and iw of two phases (U phase and W phase).
  • the coordinate conversion unit 23 converts the current values iu and iw detected by the current detector 22 from the three-phase fixed coordinate system to the dcqc axis rotation coordinate system.
  • the three-phase fixed coordinate system and the dcqc axis rotation coordinate system will be described with reference to FIG.
  • the three-phase fixed coordinate system is a fixed coordinate system composed of an ⁇ axis and a ⁇ axis.
  • the ⁇ axis is set in the U-phase direction
  • the ⁇ axis is set in a direction perpendicular to the ⁇ axis.
  • the current values iu and iw detected by the current detector 22 are expressed on such three-phase fixed coordinates.
  • the dcqc axis rotational coordinate system is a rotational coordinate system composed of a dc axis and a qc axis.
  • the dc axis is set in the direction estimated by the control device 2 as the d-axis direction (the direction in which the rotor inductance is minimum), and the qc axis is the q-axis direction in the q-axis direction (the direction in which the rotor inductance is maximum). Is set in the estimated direction.
  • the inductance ellipse in FIG. 2 indicates the inductance of the rotor.
  • the dcqc axis and the dq axis do not always coincide with each other.
  • the actual rotational phase angle ⁇ of the rotor is represented by an angle from the ⁇ axis to the d axis.
  • the estimated rotational phase angle ⁇ est of the rotor estimated by the control device 2 is represented by an angle from the ⁇ axis to the dc axis. The closer the angle between the rotational phase angle ⁇ and the estimated rotational phase angle ⁇ est, the higher the estimation accuracy of the rotational phase angle.
  • the coordinate conversion unit 23 can convert the three-phase fixed coordinate system into the dcqc axis rotation coordinate system by using the estimated rotation phase angle ⁇ est output from the speed rotation phase angle estimation unit 28.
  • the current values iu and iw coordinate-converted by the coordinate conversion unit 23 are referred to as current values idc and iqc.
  • the current value idc is a dc axis component of the current flowing through the stator
  • the current value iqc is a qc axis component of the current flowing through the stator.
  • the current command value calculation unit 24 calculates current command values idc * and iqc * based on the torque command value T * and the estimated speed ⁇ est.
  • the torque command value T * is a torque value generated in the rotor. In the present embodiment, it is assumed that the torque command value T * is input from an external device.
  • the estimated speed ⁇ est is the angular speed of the rotor estimated by the control device 2.
  • the current command value idc * is a dc axis component of a current that flows through SynRM1.
  • the current command value iqc * is the qc-axis component of the current that flows through SynRM1. Details of the current command value calculation unit 24 will be described later.
  • the voltage command value calculation unit 25 calculates voltage command values vdc * and vqc * at which the current values idc and iqc of SynRM1 are equal to the current command values idc * and iqc * .
  • the voltage command value vdc * is a dc axis component of the voltage applied to the stator of SynRM1.
  • the voltage command value vqc * is a qc-axis component of the voltage applied to the stator of SynRM1.
  • the coordinate conversion unit 26 converts the voltage command values vdc * and vqc * output from the voltage command calculation unit 25 from a dcqc axis rotation coordinate system to a three-phase fixed coordinate system. Similar to the coordinate conversion unit 23, the coordinate conversion unit 26 can convert the dcqc axis rotation coordinate system into a three-phase fixed coordinate system by using the estimated rotation phase angle ⁇ est.
  • the voltage command values vdc * and vqc * coordinate-converted by the coordinate conversion unit 26 are referred to as voltage command values vu * , vv * and vw * .
  • the voltage command value vu * is a voltage applied to the U phase of the stator
  • the voltage command value vv * is a voltage applied to the V phase of the stator
  • the voltage command value vw * is a W phase of the stator. Is the voltage applied to.
  • the PWM modulation unit 27 modulates the voltage command values vu * , vv * , vw * by PWM (Pulse-Width Modulation) using a triangular wave, and outputs a binary value corresponding to ON or OFF of each switching element of the inverter 21. Generate a control signal.
  • the PWM modulation unit 27 inputs the generated control signal to the inverter 21.
  • the speed / rotation phase angle estimator 28 determines the rotor of the SynRM 1 based on the torque command value T * , the voltage command values vdc * and vqc * , and the current values idc and iqc.
  • the speed ⁇ and the rotational phase angle ⁇ are estimated, and the estimated speed ⁇ est and the estimated rotational phase angle ⁇ est are calculated.
  • FIG. 3 is a diagram illustrating a configuration of the estimation unit 28.
  • the estimation unit 28 includes a phase difference ⁇ setting unit 31, a ⁇ voltage calculation unit 32, a ⁇ voltage estimation unit 33, a subtractor 34, a PI controller 35, an integrator 36, Is provided.
  • the phase difference ⁇ setting unit 31 outputs a phase difference ⁇ corresponding to the torque command value T * from among a plurality of phase differences ⁇ stored in advance.
  • the phase difference ⁇ is a phase difference value or range in which the voltage change is the largest due to the influence of the error ⁇ between the rotational phase angle ⁇ and the estimated rotational phase angle ⁇ est.
  • the plurality of phase differences ⁇ are calculated in advance analytically or experimentally for each torque value and stored in the phase difference ⁇ setting unit 31.
  • the ⁇ voltage calculation unit 32 calculates the voltage value v ⁇ of the ⁇ voltage based on the voltage command values vdc * and vqc * and the phase difference ⁇ set (output) by the phase difference ⁇ setting unit 31.
  • the ⁇ voltage is a feature amount that changes according to the error ⁇ .
  • the voltage value v ⁇ is calculated by the following equation, for example.
  • the voltage value v ⁇ calculated by the ⁇ voltage calculation unit 32 is input to the subtractor 34.
  • the ⁇ voltage estimation unit 33 calculates an estimated voltage value v ⁇ est of the ⁇ voltage based on the voltage values vdcest, vqcest and the phase difference ⁇ set (output) by the phase difference ⁇ setting unit 31.
  • the ⁇ voltage estimation unit 33 calculates voltage values vdcest and vqcest based on the current values idc and iqc and the estimated speed ⁇ est.
  • the voltage value vdcest is an estimated value of the dc axis component of the voltage applied to the stator of SynRM1.
  • the voltage value vqcest is an estimated value of the qc-axis component of the voltage applied to the stator of SynRM1.
  • the voltage values vdcest and vqcest are calculated by the following equations.
  • Equation (2) Rm is a winding resistance of the stator, Ld is an inductance in the d-axis direction, Lq is an inductance in the q-axis direction, and p is a differential operator (d / dt).
  • the ⁇ voltage estimation unit 33 stores these values in advance.
  • the ⁇ voltage estimation unit 33 calculates an estimated voltage value v ⁇ est of the ⁇ voltage based on the voltage values vdest, vqcest and the phase difference ⁇ .
  • the estimated voltage value v ⁇ est is calculated by the following equation, for example.
  • the estimated voltage value v ⁇ est calculated by the ⁇ voltage estimation unit 33 is input to the subtractor 34.
  • the subtracter 34 subtracts the voltage value v ⁇ from the estimated voltage value v ⁇ est to calculate an error ⁇ v ⁇ of the ⁇ voltage. Since the ⁇ voltage changes according to the error ⁇ , the error ⁇ v ⁇ is proportional to the error ⁇ .
  • the error ⁇ v ⁇ calculated by the subtractor 34 is input to the PI controller 35.
  • the PI controller 35 performs the PI control so that the error ⁇ v ⁇ becomes 0, thereby estimating the rotor speed ⁇ and calculating the estimated speed ⁇ est.
  • the estimated speed ⁇ est calculated by the PI controller 35 is sequentially fed back to the ⁇ voltage estimator 33 and input to the integrator 36.
  • the integrator 36 integrates the estimated speed ⁇ est calculated by the PI controller 35 to calculate the estimated rotational phase angle ⁇ est.
  • the estimation unit 28 can calculate the estimated speed ⁇ est and the estimated rotational phase angle ⁇ est.
  • the estimated speed ⁇ est calculated by the estimation unit 28 is input to the current command value calculation unit 24.
  • the estimated rotational phase angle ⁇ est is input to the coordinate conversion units 23 and 26 and used for coordinate conversion.
  • the estimation method of the speed ⁇ and the rotational phase angle ⁇ by the estimation unit 28 is not limited to this, and can be arbitrarily selected from known estimation methods.
  • the estimation unit 28 may estimate the rotational phase angle ⁇ by another method using a voltage generated by the flux linkage, or may estimate the rotation phase angle ⁇ by using the flux linkage itself. Then, the rotational phase angle ⁇ may be estimated by performing PI control so that the deviation of the q-axis component of the current value becomes zero.
  • the current command value calculation unit 24 calculates the current command value so that the output voltage of the SynRM 1 is equal to or higher than a predetermined target value V SET .
  • the target value V SET is a voltage value obtained experimentally or analytically in advance as an output voltage capable of accurately estimating the speed ⁇ and the rotational phase angle ⁇ .
  • FIG. 4 is a diagram showing a configuration of the current command value calculation unit 24.
  • the current command value calculation unit 24 includes a first calculation unit 41, a second calculation unit 42, and a selection unit 43.
  • the first calculator 41 generates first current command values id1 * and iq1 * based on the torque command value T * .
  • the first current command value id1 * is a dc-axis component of a current that flows through SynRM1.
  • the first current command value iq1 * is a qc-axis component of a current that flows through SynRM1.
  • the first calculator 41 calculates the first current command values id1 * and iq1 * so that the torque of SynRM1 becomes the torque command value T * .
  • FIG. 5 is a diagram for explaining a method of calculating a current command value.
  • the horizontal axis is idc
  • the vertical axis is iqc
  • the arrow 51 is the first current vector
  • the arrow 52 is the second current vector
  • the curve 53 is the equal torque curve of the torque command value T *
  • the curve 62 is the target value V. It is an isovoltage ellipse of SET .
  • the first current vector 51 is a vector corresponding to the first current command values id1 * and iq1 * .
  • the first current command value is represented by two values (id1 * , iq1 * ).
  • the current magnitude (id1 * 2 + iq1 * 2 ) 1/2 , the current phase angle ⁇ 1 It is also possible to express by. These representations can be converted to each other. Selecting an arbitrary point from the plane of FIG. 5 corresponds to calculating (selecting) the first current command value.
  • the first calculator 41 selects an arbitrary point on the equal torque curve 53 of the torque command value T * as the first current command value.
  • the first calculator 41 can select the first current command value from the equal torque curve 53 by an arbitrary method according to the control desired to be realized by the control device 2.
  • the first calculation unit 41 selects the first current command value so that the magnitude of the current flowing through the stator (the magnitude of the first current vector) is minimized.
  • the first calculation unit 41 may select the first current command value at which the current phase angle ⁇ 1 is 135 degrees ignoring the magnetic saturation, or the current phase angle ⁇ 1 in consideration of the magnetic saturation. May be selected as the first current command value at an angle more than 135 degrees.
  • the first calculation unit 41 may select the first current command value so that the efficiency and power factor of the SynRM 1 are maximized.
  • the first calculation unit 41 may select the first current command value as described above with reference to a table in which the first current command value for each torque value is stored, or may be obtained by calculation.
  • the first current command value calculated by the first calculation unit 41 is input to the selection unit 43.
  • the second calculator 42 calculates the second current command values id2 * and iq2 * based on the torque command value T * and the estimated speed ⁇ est.
  • the second current command value id2 * is a dc-axis component of a current that flows through SynRM1.
  • the current command value iq2 * is the qc-axis component of the current that flows through SynRM1.
  • the second calculation unit 42 as the torque of the SynRM1 is the torque command value T *, and, as the output voltage of SynRM1 becomes the target value V SET, the second current command value id2 *, calculates the iq2 * .
  • the second current vector 52 in FIG. 5 is a vector corresponding to the second current command values id2 * and iq2 * .
  • the second current command value is represented by two values (id2 * , iq2 * ), but the current magnitude (id2 * 2 + iq2 * 2 ) 1/2 , the current phase angle ⁇ 2, It is also possible to express by. These representations can be converted to each other. Selecting an arbitrary point from the plane of FIG. 5 corresponds to calculating (selecting) the second current command value.
  • the second calculation unit 42 uses any one of intersections of the equal torque curve 53 of the torque command value T * and the equal voltage ellipse 54 of the target value V SET as the second current command value. Select one. When there are two intersections A and B, the second calculation unit 42 preferably selects the intersection A as shown in FIG.
  • Intersection A is an intersection on the iqc axis side (q axis side) of the two intersections.
  • the equal voltage ellipse 54 is wide in the d-axis direction and narrow in the q-axis direction.
  • the magnitude (the magnitude of the second current vector) becomes small. For this reason, when the 2nd calculation part 42 selects the intersection A, the power of the control apparatus 2 can be saved. In the following, it is assumed that the second current command value is the intersection A.
  • the second calculation unit 42 may select the second current command value as described above with reference to a table in which the second current command value for each torque value is stored, or solve the following equation: You may ask for it.
  • the second current command value calculated by the second calculation unit 42 is input to the selection unit 43.
  • the selector 43 outputs the first current command value id1 * , iq1 * or the second current command value id2 * , iq2 * as the current command value idc * , iqc * .
  • the selection unit 43, the first current command value id1 *, iq1 * is output voltages V 1 corresponding to determines whether less than the target value V SET.
  • the output voltage V 1, is that of SynRM1 the output voltage when the output of the first current command value as a current command value. The determination method will be described later.
  • the selection unit 43 selects one of the first current command values id1 * and iq1 * and the second current command values id2 * and iq2 * according to the determination result, and outputs the selected current command values idc * and iqc *. To do.
  • the selection unit 43 If the output voltage V 1 is lower than the target value V SET (V 1 ⁇ V SET ), the selection unit 43, the second current command value id2 *, iq2 * of the current command value idc *, selected as iqc *.
  • the case where the output voltage V 1 is lower than the target value V SET is a case where the first current command value is included in the equal voltage ellipse 54 as shown in FIG.
  • the second current command value is closer to the iqc axis (q axis) than the first current command value. That is, the current phase angle ⁇ 2 of the second current command value is on the iqc axis (q axis) side from the current phase angle ⁇ 1 of the first current command value.
  • the q-axis component iq2 * of the second current command value is larger than the q-axis component iq1 * of the first current command value, and the d-axis component id2 * of the second current command value is the d-axis of the first current command value. Less than component id1 * .
  • the selection unit 43 selects the first current command values id1 * and iq1 * as the current command values idc * and iqc * .
  • the case where the output voltage V 1 is equal to or higher than the target value V SET is a case where the second current command value is on the equal voltage ellipse 54 or outside the equal voltage ellipse 54.
  • the first current command value is closer to the iqc axis (q axis) than the second current command value. That is, the current phase angle ⁇ 1 of the first current command value is on the iqc axis (q axis) side from the current phase angle ⁇ 2 of the second current command value. Further, q-axis component of the first current command value iq1 * is, q-axis component iq2 * greater second current command value, d-axis component of the first current command value id1 * is, d-axis of the second current command value Less than component id2 * .
  • the selection unit 43 outputs the selected first current command value or second current command value as a current command value.
  • the current command value output by the selection unit 43 is input to the voltage command calculation unit 25.
  • Selection unit 43 acquires the output voltage V 1, by comparing the target value V SET, a determination.
  • the selection unit 43 may select the output voltage V 1 with reference to a table in which the output voltage for each first current command value is stored, or may be calculated by the following equation.
  • the selection unit 43 may perform the determination by comparing the current phase angle ⁇ 1 of the first current command value with the current phase angle ⁇ 2 of the second current command value. Selection unit 43, when the current phase angle ⁇ 2 of the q-axis side from the current phase angle .beta.1, determines the output voltage V 1 is lower than the target value V SET.
  • the selection unit 43 a q-axis component iq1 * of the first current command value, the q-axis component iq2 * of the second current command value, by comparing may perform determination.
  • Selection unit 43 when the q-axis component iq1 * is the q-axis component iq2 * smaller, determines the output voltage V 1 is lower than the target value V SET.
  • the selection unit 43 may perform the determination by comparing the d-axis component id1 * of the first current command value with the d-axis component id2 * of the second current command value. Selection unit 43 determines a d-axis component id1 * if d-axis component id2 * is greater than the output voltage V 1 is lower than the target value V SET.
  • the current command value calculation unit 24 calculates the current command values idc * and iqc * in this way, so that the output voltage of the SynRM1 is set as a target for an arbitrary torque command value T * as shown by a solid line in FIG. It can be greater than or equal to the value V SET .
  • control device 2 since the control device 2 according to the present embodiment can set the output voltage to be equal to or higher than the target value V SET even in the case of a low load, based on the output voltage of SynRM1, the speed ⁇ and rotation of the SynRM1 The phase angle ⁇ can be estimated with high accuracy.
  • control device 2 controls the operation of the SynRM 1
  • the control device 2 can also be used as a PMSM control device with a small magnet magnetic flux.
  • PMSM having a small magnet magnetic flux has a low induced voltage at low load, similar to SynRM1.
  • FIG. 7 is a diagram illustrating a configuration of the first modification of the current command value calculation unit 24.
  • the second calculation unit 43 includes a q-axis component calculation unit 44 and a d-axis component calculation unit 45.
  • Other configurations of the current command value calculation unit 24 are the same as those in FIG.
  • the q-axis component calculation unit 44 calculates the q-axis component iq2 * of the second current command value based on the estimated speed ⁇ est. Specifically, the q-axis component calculation unit 44, the current value iq2 * to calculate the q-axis component iq2 * to be inversely proportional with the estimated speed .omega.est.
  • the q-axis component iq2 * is calculated by the following equation, for example.
  • the d-axis component calculation unit 45 calculates the d-axis component id2 * of the second current command value based on the torque command value T * and the q-axis component iq2 * .
  • the d-axis component id2 * is calculated by the following equation, for example.
  • control device 2 can set the output voltage at a low load to be equal to or higher than the target value V SET .
  • current loss in the high speed region can be suppressed.
  • FIG. 8 is a diagram illustrating a configuration of a second modification of the current command value calculation unit 24. As shown in FIG. 8, the current command value calculation unit 24 calculates current command values iqc * and idc * based on the torque command value T * and the estimated speed ⁇ est.
  • the current command value calculation unit 24 calculates the q-axis component iqc * such that the q-axis component iqc * of the current command value is equal to or greater than a predetermined target value Iq SET , as shown in FIG. . Further, as shown in FIG. 10, the current command value calculation unit 24 generates the q-axis component iqc * so that the target value Iq SET decreases as the estimated speed ⁇ est increases.
  • the motor output voltage at the time of low load can be set to the target value V SET or more.
  • current loss in the high speed region can be suppressed.
  • Modification 3 of the current command value calculation unit 24 according to the present embodiment will be described with reference to FIGS. 11 and 12.
  • the current command value calculation unit 24 described above calculates a current command value at which the output voltage is equal to or higher than the target value V SET by selecting one of the two current command values.
  • the current command value at which the output voltage is equal to or higher than the target value V SET is calculated by correcting one current command value.
  • FIG. 11 is a diagram illustrating a configuration of a third modification of the current command value calculation unit 24.
  • the current command value calculation unit 24 includes a q-axis component calculation unit 46, a correction unit 47, an adder 48, and a d-axis component calculation unit 49.
  • the q-axis component calculation unit 46 calculates a current command value iq3 * based on the torque command value T * .
  • the current command value iq3 * is a value before correction of the q-axis component of the current command value.
  • the current command value calculation unit 24 in the present modification calculates the q-axis component of the current command value by correcting the current command value iq3 * .
  • the current command value iq3 * is calculated by an arbitrary method according to the control to be realized by the control device so that the torque of the SynRM1 becomes the torque command value T * .
  • the q-axis component calculation unit 46 calculates the current command value iq3 * so that the magnitude of the current flowing through the stator is minimized.
  • the q-axis component calculation unit 46 may calculate the current command value iq3 * at which the current phase angle ⁇ is 135 degrees, ignoring the magnetic saturation, or the current phase angle in consideration of the magnetic saturation.
  • a current command value iq3 * at which ⁇ is an angle advanced from 135 degrees may be calculated.
  • the q-axis component calculation unit 46 may calculate the current command value iq3 * so that the efficiency and power factor of SynRM1 are maximized.
  • the current command value iq3 * calculated by the q-axis component calculation unit 46 is input to the adder 48.
  • the correction unit 47 receives the voltage command values vdc * and vqc * from the voltage command value calculation unit 25. Correction unit 47 calculates correction value ⁇ iq * based on voltage command values vdc * and vqc * . Correction value? Iq * corrects the current command value IQ3 *, a current command value for calculating a q-axis component of the current command value iqc *. Details of the correction unit 47 will be described later.
  • the correction value ⁇ iq * calculated by the correction unit 47 is input to the adder 48.
  • the adder 48 adds the correction value? Iq * to the current command value IQ3 *, calculates the current command value iqc *.
  • the current command value iqc * calculated by the adder 48 is input to the d-axis component calculation unit 49.
  • the d-axis component calculation unit 49 calculates the d-axis component idc * of the current command value based on the torque command value T * and the current command value iqc * so that the torque of the SynRM1 becomes the torque command value T *. calculate.
  • the d-axis component idc * is calculated by the following equation, for example.
  • FIG. 12 is a diagram illustrating a configuration of the correction unit 47.
  • the correction unit 47 includes an output voltage calculation unit 61, a subtracter 62, a limiter 63, and a PI controller 64.
  • the output voltage calculation unit 61 receives the voltage command values vdc * and vqc * calculated by the voltage command value calculation unit 25. Output voltage calculation unit 61, voltage command value vdc *, based on vqc *, calculates the output voltage V 1 of the SynRM1. Output voltages V 1, for example, is calculated by the following equation.
  • the output voltage V 1 calculated by the output voltage calculation unit 61 is input to the subtractor 62.
  • the subtractor 62 subtracts the output voltage V 1 from the output voltage target value V SET to calculate an output voltage error ⁇ V.
  • the error ⁇ V calculated by the subtractor 62 is input to the limiter 63.
  • the limiter 63 limits the error ⁇ V to 0 or more. That is, only an error ⁇ V of 0 or more is output.
  • the error ⁇ V output from the limiter 63 is input to the PI controller 64.
  • the PI controller 64 performs PI control based on the error ⁇ V limited to 0 or more, and calculates a correction value ⁇ iq * at which the output voltage V 1 of the SynRM 1 is equal to or higher than the target value V SET .
  • the correction value ⁇ iq * calculated by the PI controller 64 is input to the adder 48.
  • Modification 4 of the current command value calculation unit 24 will be described with reference to FIG.
  • a current command value at which the output voltage is equal to or higher than the target value V SET is calculated.
  • a correction value ⁇ iq * is generated by the correction unit 47 and added to the current command value iq3 * by the adder 48 to generate a current command value iqc * .
  • the lower limiter 91 by limiting the lower limit of the current command value IQ3 *, generates a current command value iqc *.
  • iq SET is calculated from the minimum speed ⁇ MIN used in the present embodiment and the output voltage target value V SET by the following equation.
  • the output voltage V 1 can be made equal to or higher than the target value V SET even when the load is low.
  • the output voltage V 1 can be made equal to or higher than the target value V SET with a simple configuration in which only the lower limiter 91 is added to the first embodiment and Modifications 1 to 3.
  • FIG. 14 is a diagram illustrating a configuration of a system according to the present embodiment. As shown in FIG. 14, the control device 2 according to the present embodiment further includes a torque command value calculation unit 29. Other configurations are the same as those of the first embodiment.
  • the torque command value calculation unit (speed controller) 29 calculates a torque command value T * based on the speed command value ⁇ * and the estimated speed ⁇ est.
  • the speed command value ⁇ * is an angular speed for rotating the rotor of SynRM1.
  • the torque command value T * calculated by the torque command value calculation unit 29 is input to the current command value calculation unit 24.
  • the speed command value ⁇ * is input from an external device.
  • FIG. 15 is a diagram illustrating a configuration of the torque command value calculation unit 29. As shown in FIG. 15, the torque command value calculation unit 29 includes a subtracter 71 and a PI controller 72.
  • the subtracter 71 receives the speed command value ⁇ * and the estimated speed ⁇ est.
  • the subtracter 71 subtracts the estimated speed ⁇ est from the speed command value ⁇ * to calculate an angular speed error ⁇ of the rotor.
  • the error ⁇ calculated by the subtracter 71 is input to the PI controller 72.
  • the PI controller 72 performs PI control based on the error ⁇ , and calculates a torque command value T * at which the error ⁇ becomes zero.
  • the torque command value T * calculated by the PI controller 72 is input to the current command value calculation unit 24.
  • the current command value calculation unit 24 may have the configuration described in the first embodiment, but may have another configuration.
  • FIG. 16 is a diagram illustrating another configuration of the current command value calculation unit 24.
  • the current command value calculation unit 24 calculates current command values idc * and idc * based on the torque command value T * and the phase angle command value ⁇ * .
  • the current command value calculation unit 24 includes an output voltage calculation unit 81, a subtracter 82, a limiter 83, a PI controller 84, a phase angle command value calculation unit 85, and a subtractor 86. And a current command value acquisition unit 87.
  • the output voltage calculation unit 81 receives the voltage command values vdc * and vqc * calculated by the voltage command value calculation unit 25. Output voltage calculation unit 81, voltage command value vdc *, based on vqc *, calculates the output voltage V 1 of the SynRM1. The output voltage V 1 calculated by the output voltage calculation unit 81 is input to the subtractor 82.
  • the subtracter 82 subtracts the output voltage V 1 from the output voltage target value V SET to calculate an output voltage error ⁇ V.
  • the error ⁇ V calculated by the subtractor 82 is input to the limiter 83.
  • the limiter 83 limits the error ⁇ V to 0 or more. That is, only an error ⁇ V of 0 or more is output.
  • the error ⁇ V output from the limiter 83 is input to the PI controller 84.
  • the PI controller 84 performs PI control based on the error ⁇ V limited to 0 or more, and calculates a correction value ⁇ * at which the output voltage V 1 of the SynRM 1 is equal to or higher than the target value V SET .
  • Correction value [Delta] [beta] * is the phase angle command value to calculate the phase angle command value beta * to correct the phase angle command value .beta.1 *.
  • the correction value ⁇ * calculated by the PI controller 84 is input to the subtractor 86.
  • the phase angle command value calculation unit 85 calculates the phase angle command value ⁇ 1 * based on the torque command value T * .
  • the phase angle command value ⁇ 1 * is a value before correction of the phase angle command value ⁇ * .
  • the phase angle command value ⁇ 1 * is calculated by an arbitrary method according to the control desired to be realized by the control device so that the torque of the SynRM1 becomes the torque command value T * .
  • the phase angle command value calculation unit 85 calculates the phase angle command value ⁇ 1 * so that the magnitude of the current flowing through the stator is minimized.
  • the phase angle command value calculation unit 85 may ignore the magnetic saturation and set the phase angle command value ⁇ 1 * to 135 degrees, or consider the magnetic saturation and set the phase angle command value ⁇ 1 * from 135 degrees. It may be an advanced angle. Further, the phase angle command value calculation unit 85 may calculate the phase angle command value ⁇ 1 * so that the efficiency and the power factor of SynRM1 are maximized.
  • phase angle command value ⁇ 1 * calculated by the phase angle command value calculation unit 85 is input to the subtractor 86.
  • the subtracter 86 subtracts the correction value ⁇ * from the phase angle command value ⁇ 1 * to calculate the phase angle command value ⁇ * .
  • the phase angle command value ⁇ * calculated by the subtractor 86 is input to the current command value acquisition unit 87.
  • the current command value acquisition unit 87 calculates current command values idc * and iqc * based on the torque command value T * and the phase angle command value ⁇ * . Specifically, the current command value acquisition unit 87 may acquire the coordinates of the point where the current phase angle becomes the phase angle command value ⁇ * on the equal torque curve of the torque command value T * . The current command values idc * and iqc * acquired by the current command value acquisition unit 87 are input to the voltage command value calculation unit 25.
  • the current command value calculation unit 24 in FIG. 16 corrects the current phase angle so that the output voltage V 1 is equal to or higher than the target value V SET . Specifically, the current phase angle is corrected in the delay direction so that the current command value moves to the q-axis side. Thus, it is possible to output voltage V 1 of the low load and the target value V SET more.
  • control device 2 according to the present embodiment has lower required accuracy of torque control than the first embodiment. This is because the control device 2 according to the present embodiment performs control so that the speed ⁇ of the SynRM 1 becomes the speed command value ⁇ * .
  • the current command value calculation unit 24 can calculate the current command value without using a table (such as a table in which the output voltage for each current command value is stored) for accurately controlling the torque of the SynRM1. .
  • the present invention is not limited to the above-described embodiments as they are, and can be embodied by modifying the constituent elements without departing from the scope of the invention in the implementation stage.
  • various inventions can be formed by appropriately combining a plurality of constituent elements disclosed in the above embodiments. Further, for example, a configuration in which some components are deleted from all the components shown in each embodiment is also conceivable. Furthermore, you may combine suitably the component described in different embodiment.
  • 1 synchronous reluctance motor (SynRM), 2: inverter control device, 21: inverter, 22: current detector, 23: coordinate conversion unit, 24: current command value calculation unit, 25: voltage command value calculation unit, 26: Coordinate conversion unit, 27: PWM modulator, 28: Speed / rotation phase angle estimation unit, 29: Torque command value calculation unit, 31: Phase difference ⁇ setting unit, 32: ⁇ voltage calculation unit, 33: ⁇ voltage estimation unit, 34: subtractor, 35: PI controller, 36: integrator, 41: first calculator, 42: second calculator, 43: selector, 44: q-axis component calculator, 45: d-axis component calculator 46: q-axis component calculation unit, 47: correction unit, 48: adder, 49d-axis component calculation unit, 51: first current vector, 52: second current vector, 53: isotorque curve, 54: isovoltage ellipse 61: output voltage calculation unit, 62: subtractor, 63 Limiter, 64: PI controller, 71

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Ac Motors In General (AREA)
  • Control Of Motors That Do Not Use Commutators (AREA)
  • Control Of Electric Motors In General (AREA)

Abstract

【課題】　高速かつ低負荷の状態におけるモータの回転位相角を精度よく推定できるインバータ制御装置及び駆動システムを提供する。 【解決手段】　一実施形態に係るインバータ制御装置は、インバータ主回路と、電流検出器と、電流指令値算出部と、電圧指令値算出部と、推定部と、を備える。インバータ主回路は、所定の回転駆動対象と電気的に接続可能である。電流検出器は、インバータ主回路から出力される電流値を検出する。電流指令値算出部は、インバータ主回路から出力される出力電圧が所定の目標値以上となる電流指令値を算出する。電圧指令値算出部は、電流値が、電流指令値と等しくなる電圧指令値を算出する。推定部は、電圧指令値及び電流値に基づいて、回転駆動対象の推定回転位相角を算出する。

Description

インバータ制御装置及びモータ駆動システム
　本発明の実施形態は、インバータ制御装置及びモータ駆動システムに関する。
　従来、永久磁石同期モータ（ＰＭＳＭ）やシンクロナスリラクタンスモータ（ＳｙｎＲＭ）の回転位相角センサレス制御装置において、高速域では誘起電圧を利用した回転位相角の推定方法が利用されている。しかしながら、ＳｙｎＲＭや、磁石磁束が小さいＰＭＳＭでは、高速域であっても、低負荷の状態では誘起電圧が小さいため、回転位相角の推定精度が悪化するという問題があった。
特開２００３－２５０２９３号公報 特開２００６－７４９０２号公報
　高速かつ低負荷の状態におけるモータの回転位相角を精度よく推定できるインバータ制御装置及びモータ駆動システムを提供する。
　一実施形態に係るインバータ制御装置は、インバータ主回路と、電流検出器と、電流指令値算出部と、電圧指令値算出部と、推定部と、を備える。インバータ主回路は、所定の回転駆動対象と電気的に接続可能である。電流検出器は、インバータ主回路から出力される電流値を検出する。電流指令値算出部は、インバータ主回路から出力される出力電圧が所定の目標値以上となる電流指令値を算出する。電圧指令値算出部は、電流値が、電流指令値と等しくなる電圧指令値を算出する。推定部は、電圧指令値及び電流値に基づいて、回転駆動対象の推定回転位相角を算出する。
第１実施形態に係るモータ駆動システムの構成を示す図。 三相固定座標系及びｄｃｑｃ軸回転座標系を説明する図。 図１の速度・回転位相角推定部の構成を示す図。 図１の電流指令値算出部の構成を示す図。 電流指令値の算出方法を説明する図。 図１の制御装置のトルク指令値Ｔ＊と出力電圧との関係を示す図。 電流指令値算出部の変形例１の構成を示す図。 電流指令値算出部の変形例２の構成を示す図。 変形例２における電流指令値の算出方法を説明する図。 変形例２における電流指令値の算出方法を説明する図。 電流指令値算出部の変形例３の構成を示す図。 補正部の構成を示す図。 電流指令値算出部の変形例４の構成を示す図。 第２実施形態に係るモータ駆動システムの構成を示す図。 図１４のトルク指令値算出部の構成を示す図。 図１４の電流指令値算出部の構成を示す図。
　以下、本発明の実施形態について図面を参照して説明する。
（第１実施形態）
　第１実施形態に係るモータ駆動システム（以下、「システム」という）について、図１～図１３を参照して説明する。図１は、本実施形態に係るシステムの構成を示す図である。図１に示すように、本実施形態に係るシステムは、モータ１と、インバータ制御装置２と、を備える。
　モータ１は、制御装置２の回転駆動対象であり、制御装置２に接続される。以下では、モータ１がシンクロナスリラクタンスモータ（以下、「ＳｙｎＲＭ１」という）である場合を例として説明する。ＳｙｎＲＭ１は、固定子と、回転子と、を備える。固定子は、３つの励磁相（Ｕ相、Ｖ相、及びＷ相）を有する。固定子は、各励磁相に流れる３相交流電流によって磁界を発生させる。回転子は、永久磁石を有さない。回転子は、固定子が発生させた磁界との磁気的相互作用により回転する。
　インバータ制御装置２（以下、「制御装置２」という）は、ＳｙｎＲＭ１の回転位相角を、センサレスで制御する。図１に示すように、本実施形態に係る制御装置２は、インバータ２１と、電流検出器２２と、座標変換部２３と、電流指令値算出部２４と、電圧指令値算出部２５と、座標変換部２６と、ＰＷＭ変調部２７と、速度・回転位相角推定部２８と、を備える。
　インバータ２１は、スイッチング素子（トランジスタ）を備える回路である。インバータ２１は、スイッチング素子のＯＮ／ＯＦＦを切替えることにより、電源（図示省略）からの電力を交流に変換して、ＳｙｎＲＭ１に供給する。インバータ２１は、各スイッチング素子のＯＮ／ＯＦＦを制御する制御信号をＰＷＭ変調部２２から入力される。
　電流検出器２２は、ＳｙｎＲＭ１の固定子に流れる３相交流電流のうち、２相又は３相の電流値を検出する。図１は、２相（Ｕ相及びＷ相）の電流値ｉｕ，ｉｗを検出する構成を示している。
　座標変換部２３は、電流検出器２２が検出した電流値ｉｕ，ｉｗを、三相固定座標系からｄｃｑｃ軸回転座標系に座標変換する。ここで、三相固定座標系及びｄｃｑｃ軸回転座標系について、図２を参照して説明する。
　図２に示すように、三相固定座標系は、α軸とβ軸とからなる固定座標系である。図２において、α軸は、Ｕ相方向に設定され、β軸は、α軸と垂直な方向に設定されている。電流検出器２２により検出された電流値ｉｕ，ｉｗは、このような三相固定座標上で表される。
　これに対して、ｄｃｑｃ軸回転座標系は、ｄｃ軸とｑｃ軸とからなる回転座標系である。ｄｃ軸は、制御装置２がｄ軸方向（回転子のインダクタンスが最小の方向）と推定した方向に設定され、ｑｃ軸は、制御装置２がｑ軸方向（回転子のインダクタンスが最大の方向）と推定した方向に設定される。図２のインダクタンス楕円は、回転子のインダクタンスを示している。
　図２に示すように、ｄｃｑｃ軸と、ｄｑ軸と、は必ずしも一致するとは限らない。回転子の実際の回転位相角θは、α軸からｄ軸までの角度で表される。また、制御装置２が推定した回転子の推定回転位相角θｅｓｔは、α軸からｄｃ軸までの角度で表される。回転位相角θと推定回転位相角θｅｓｔとの角度が近いほど、回転位相角の推定精度が高いことを意味する。
　座標変換部２３は、速度回転位相角推定部２８が出力した推定回転位相角θｅｓｔを用いることにより、三相固定座標系をｄｃｑｃ軸回転座標系に変換することができる。以下では、座標変換部２３が座標変換した電流値ｉｕ，ｉｗを、電流値ｉｄｃ，ｉｑｃという。電流値ｉｄｃは、固定子に流れる電流のｄｃ軸成分であり、電流値ｉｑｃは、固定子に流れる電流のｑｃ軸成分である。
　電流指令値算出部２４は、トルク指令値Ｔ＊及び推定速度ωｅｓｔに基づいて、電流指令値ｉｄｃ＊，ｉｑｃ＊を算出する。トルク指令値Ｔ＊とは、回転子に発生させるトルク値である。本実施形態では、トルク指令値Ｔ＊は、外部装置から入力されるものとする。推定速度ωｅｓｔとは、制御装置２が推定した回転子の角速度のことである。電流指令値ｉｄｃ＊とは、ＳｙｎＲＭ１に流す電流のｄｃ軸成分である。電流指令値ｉｑｃ＊とは、ＳｙｎＲＭ１に流す電流のｑｃ軸成分である。電流指令値算出部２４の詳細については後述する。
　電圧指令値算出部２５（電流制御部）は、ＳｙｎＲＭ１の電流値ｉｄｃ，ｉｑｃが、電流指令値ｉｄｃ＊，ｉｑｃ＊と等しくなる電圧指令値ｖｄｃ＊，ｖｑｃ＊を算出する。電圧指令値ｖｄｃ＊は、ＳｙｎＲＭ１の固定子に印加する電圧のｄｃ軸成分である。電圧指令値ｖｑｃ＊は、ＳｙｎＲＭ１の固定子に印加する電圧のｑｃ軸成分である。
　座標変換部２６は、電圧指令算出部２５が出力した電圧指令値ｖｄｃ＊，ｖｑｃ＊を、ｄｃｑｃ軸回転座標系から三相固定座標系に座標変換する。座標変換部２６は、座標変換部２３と同様、推定回転位相角θｅｓｔを用いることにより、ｄｃｑｃ軸回転座標系を三相固定座標系に変換することができる。以下では、座標変換部２６が座標変換した電圧指令値ｖｄｃ＊，ｖｑｃ＊を、電圧指令値ｖｕ＊，ｖｖ＊，ｖｗ＊という。電圧指令値ｖｕ＊は、固定子のＵ相に印加する電圧であり、電圧指令値ｖｖ＊は、固定子のＶ相に印加する電圧であり、電圧指令値ｖｗ＊は、固定子のＷ相に印加する電圧である。
　ＰＷＭ変調部２７は、電圧指令値ｖｕ＊，ｖｖ＊，ｖｗ＊を、三角波を用いたＰＷＭ（Pulse-Width Modulation）によって変調し、インバータ２１の各スイッチング素子のＯＮ又はＯＦＦに対応した２値の制御信号を生成する。ＰＷＭ変調部２７は、生成した制御信号をインバータ２１に入力する。
　速度・回転位相角推定部２８（以下、「推定部２８」という）は、トルク指令値Ｔ＊、電圧指令値ｖｄｃ＊，ｖｑｃ＊、及び電流値ｉｄｃ，ｉｑｃに基づいて、ＳｙｎＲＭ１の回転子の速度ω及び回転位相角θを推定し、推定速度ωｅｓｔ及び推定回転位相角θｅｓｔを算出する。
　図３は、推定部２８の構成を示す図である。図３に示すように、推定部２８は、位相差δ設定部３１と、γ電圧算出部３２と、γ電圧推定部３３と、減算器３４と、ＰＩ制御器３５と、積分器３６と、を備える。
　位相差δ設定部３１は、予め記憶した複数の位相差δの中から、トルク指令値Ｔ＊に応じた位相差δを出力する。位相差δとは、回転位相角θと推定回転位相角θｅｓｔとの誤差Δθの影響による、電圧の変化が一番大きくなる位相差の値又は範囲のことである。複数の位相差δは、トルク値毎に解析的又は実験的に予め算出され、位相差δ設定部３１に記憶される。
　γ電圧算出部３２は、電圧指令値ｖｄｃ＊，ｖｑｃ＊と、位相差δ設定部３１が設定（出力）した位相差δと、に基づいて、γ電圧の電圧値ｖγを算出する。γ電圧とは、誤差Δθに応じて変化する特徴量である。電圧値ｖγは、例えば、以下の式により算出される。
Figure JPOXMLDOC01-appb-M000001
　γ電圧算出部３２が算出した電圧値ｖγは、減算器３４に入力される。
　γ電圧推定部３３は、電圧値ｖｄｃｅｓｔ，ｖｑｃｅｓｔと、位相差δ設定部３１が設定（出力）した位相差δと、に基づいて、γ電圧の推定電圧値ｖγｅｓｔを算出する。
　まず、γ電圧推定部３３は、電流値ｉｄｃ，ｉｑｃと、推定速度ωｅｓｔと、に基づいて、電圧値ｖｄｃｅｓｔ，ｖｑｃｅｓｔを算出する。電圧値ｖｄｃｅｓｔは、ＳｙｎＲＭ１の固定子に印加された電圧のｄｃ軸成分の推定値である。電圧値ｖｑｃｅｓｔは、ＳｙｎＲＭ１の固定子に印加された電圧のｑｃ軸成分の推定値である。電圧値ｖｄｃｅｓｔ，ｖｑｃｅｓｔは、以下の式により算出される。
Figure JPOXMLDOC01-appb-M000002
　式（２）において、Ｒｍは固定子の巻線抵抗、Ｌｄはｄ軸方向のインダクタンス、Ｌｑはｑ軸方向のインダクタンス、ｐは微分演算子（ｄ／ｄｔ）である。γ電圧推定部３３は、これらの値を予め記憶する。
　次に、γ電圧推定部３３は、電圧値ｖｄｃｅｓｔ，ｖｑｃｅｓｔと、位相差δと、に基づいて、γ電圧の推定電圧値ｖγｅｓｔを算出する。推定電圧値ｖγｅｓｔは、例えば、以下の式により算出される。
Figure JPOXMLDOC01-appb-M000003
　γ電圧推定部３３が算出した推定電圧値ｖγｅｓｔは、減算器３４に入力される。
　減算器３４は、推定電圧値ｖγｅｓｔから電圧値ｖγを減算し、γ電圧の誤差Δｖγを算出する。γ電圧は、誤差Δθに応じて変化するため、誤差Δｖγは、誤差Δθに比例する。減算器３４が算出した誤差Δｖγは、ＰＩ制御器３５に入力される。
　なお、γ電圧を式（１）で算出することにより、誤差Δθに対する誤差Δｖγの線形性を向上させることができる。すなわち、誤差Δｖγが誤差Δθに比例する誤差Δθの範囲を広げることができる。
　ＰＩ制御器３５は、誤差Δｖγが０となるようにＰＩ制御を行うことで、回転子の速度ωを推定し、推定速度ωｅｓｔを算出する。ＰＩ制御器３５が算出した推定速度ωｅｓｔは、γ電圧推定部３３に逐次フィードバックされると共に、積分器３６に入力される。
　積分器３６は、ＰＩ制御器３５が算出した推定速度ωｅｓｔを積分し、推定回転位相角θｅｓｔを算出する。
　以上のような構成により、推定部２８は、推定速度ωｅｓｔ及び推定回転位相角θｅｓｔを算出することができる。推定部２８が算出した推定速度ωｅｓｔは、電流指令値算出部２４に入力される。また、推定回転位相角θｅｓｔは、座標変換部２３，２６に入力され、座標変換に利用される。
　なお、推定部２８による速度ω及び回転位相角θの推定方法はこれに限られず、既知の推定方法から任意に選択可能である。例えば、推定部２８は、鎖交磁束により生じる電圧を利用する他の方法により回転位相角θを推定してもよいし、鎖交磁束自体を利用して回転位相角θを推定してもよいし、電流値のｑ軸成分の偏差が０になるようにＰＩ制御することにより回転位相角θを推定してもよい。
　ここで、電流指令値算出部２４について詳細に説明する。本実施形態に係る電流指令値算出部２４は、ＳｙｎＲＭ１の出力電圧が所定の目標値ＶＳＥＴ以上となるように、電流指令値を算出する。目標値ＶＳＥＴは、速度ω及び回転位相角θを精度よく推定可能な出力電圧として、実験的又は解析的に予め求められた電圧値である。
　ここで、図４は、電流指令値算出部２４の構成を示す図である。図４に示すように、電流指令値算出部２４は、第１算出部４１と、第２算出部４２と、選択部４３と、を備える。
　第１算出部４１は、トルク指令値Ｔ＊に基づいて、第１電流指令値ｉｄ１＊，ｉｑ１＊を生成する。第１電流指令値ｉｄ１＊とは、ＳｙｎＲＭ１に流す電流のｄｃ軸成分である。第１電流指令値ｉｑ１＊とは、ＳｙｎＲＭ１に流す電流のｑｃ軸成分である。第１算出部４１は、ＳｙｎＲＭ１のトルクがトルク指令値Ｔ＊となるように、第１電流指令値ｉｄ１＊，ｉｑ１＊を算出する。
　図５は、電流指令値の算出方法を説明する図である。図５において、横軸はｉｄｃ、縦軸はｉｑｃ、矢印５１は第１電流ベクトル、矢印５２は第２電流ベクトル、曲線５３はトルク指令値Ｔ＊の等トルク曲線、及び曲線６２は目標値ＶＳＥＴの等電圧楕円である。
　第１電流ベクトル５１は、第１電流指令値ｉｄ１＊，ｉｑ１＊に対応するベクトルである。以上の説明では、第１電流指令値は、（ｉｄ１＊，ｉｑ１＊）の２値で表されたが、電流の大きさ（ｉｄ１＊２＋ｉｑ１＊２）１／２と、電流位相角β１と、で表すことも可能である。これらの表現は、相互に変換可能である。図５の平面上から任意の点を選択することは、第１電流指令値を算出（選択）することと対応する。
　図５に示すように、第１算出部４１は、第１電流指令値として、トルク指令値Ｔ＊の等トルク曲線５３上の任意の点を選択する。第１算出部４１は、制御装置２により実現したい制御に応じた任意の方法で、等トルク曲線５３上から第１電流指令値を選択することができる。
　例えば、第１算出部４１は、固定子に流す電流の大きさ（第１電流ベクトルの大きさ）が最小となるように、第１電流指令値を選択する。この場合、第１算出部４１は、磁気飽和を無視して、電流位相角β１が１３５度となる第１電流指令値を選択してもよいし、磁気飽和を考慮して、電流位相角β１が１３５度より進んだ角度となる第１電流指令値を選択してもよい。また、第１算出部４１は、ＳｙｎＲＭ１の効率や力率が最大となるように第１電流指令値を選択してもよい。
　第１算出部４１は、上記のような第１電流指令値を、トルク値毎の第１電流指令値が記憶されたテーブルを参照して選択してもよいし、計算により求めてもよい。第１算出部４１により算出された第１電流指令値は、選択部４３に入力される。
　第２算出部４２は、トルク指令値Ｔ＊及び推定速度ωｅｓｔに基づいて、第２電流指令値ｉｄ２＊，ｉｑ２＊を算出する。第２電流指令値ｉｄ２＊とは、ＳｙｎＲＭ１に流す電流のｄｃ軸成分である。電流指令値ｉｑ２＊とは、ＳｙｎＲＭ１に流す電流のｑｃ軸成分である。第２算出部４２は、ＳｙｎＲＭ１のトルクがトルク指令値Ｔ＊となるように、かつ、ＳｙｎＲＭ１の出力電圧が目標値ＶＳＥＴとなるように、第２電流指令値ｉｄ２＊，ｉｑ２＊を算出する。
　図５の第２電流ベクトル５２は、第２電流指令値ｉｄ２＊，ｉｑ２＊に対応するベクトルである。以上の説明では、第２電流指令値は、（ｉｄ２＊，ｉｑ２＊）の２値で表されたが、電流の大きさ（ｉｄ２＊２＋ｉｑ２＊２）１／２と、電流位相角β２と、で表すことも可能である。これらの表現は、相互に変換可能である。図５の平面上から任意の点を選択することは、第２電流指令値を算出（選択）することと対応する。
　図５に示すように、第２算出部４２は、第２電流指令値として、トルク指令値Ｔ＊の等トルク曲線５３と、目標値ＶＳＥＴの等電圧楕円５４と、の交点のいずれか１つを選択する。第２算出部４２は、２つの交点Ａ，Ｂが存在する場合、図５に示すように、交点Ａを選択するのが好ましい。
　交点Ａとは、２つの交点のうち、ｉｑｃ軸側（ｑ軸側）の交点のことである。一般に、等電圧楕円５４は、ｄ軸方向に広く、ｑ軸方向に狭くなるため、ｑ軸側の交点Ａを選択すると、ｄ軸側の交点Ｂを選択した場合より、固定子に流す電流の大きさ（第２電流ベクトルの大きさ）が小さくなる。このため、第２算出部４２が交点Ａを選択することにより、制御装置２を省電力することができる。以下では、第２電流指令値は交点Ａであるものとする。
　なお、第２算出部４２は、上記のような第２電流指令値を、トルク値毎の第２電流指令値が記憶されたテーブルを参照して選択してもよいし、以下の式を解くことにより求めてもよい。
Figure JPOXMLDOC01-appb-M000004
　第２算出部４２により算出された第２電流指令値は、選択部４３に入力される。
　選択部４３は、第１電流指令値ｉｄ１＊，ｉｑ１＊又は第２電流指令値ｉｄ２＊，ｉｑ２＊を電流指令値ｉｄｃ＊，ｉｑｃ＊として出力する。まず、選択部４３は、第１電流指令値ｉｄ１＊，ｉｑ１＊に応じた出力電圧Ｖ１が、目標値ＶＳＥＴより低いか否か判定する。出力電圧Ｖ１とは、電流指令値として第１電流指令値を出力した場合のＳｙｎＲＭ１の出力電圧のことである。判定方法については後述する。
　次に、選択部４３は、判定結果に応じて第１電流指令値ｉｄ１＊，ｉｑ１＊及び第２電流指令値ｉｄ２＊，ｉｑ２＊の一方を選択し、電流指令値ｉｄｃ＊，ｉｑｃ＊として出力する。
　出力電圧Ｖ１が目標値ＶＳＥＴより低い場合（Ｖ１＜ＶＳＥＴ）、選択部４３は、第２電流指令値ｉｄ２＊，ｉｑ２＊を電流指令値ｉｄｃ＊，ｉｑｃ＊として選択する。出力電圧Ｖ１が目標値ＶＳＥＴより低い場合とは、図５に示すように、第１電流指令値が、等電圧楕円５４に含まれる場合である。
　このとき、第２電流指令値は、第１電流指令値よりｉｑｃ軸（ｑ軸）に近い。すなわち、第２電流指令値の電流位相角β２は、第１電流指令値の電流位相角β１よりｉｑｃ軸（ｑ軸）側である。また、第２電流指令値のｑ軸成分ｉｑ２＊は、第１電流指令値のｑ軸成分ｉｑ１＊より大きく、第２電流指令値のｄ軸成分ｉｄ２＊は、第１電流指令値のｄ軸成分ｉｄ１＊より小さい。
　一方、出力電圧Ｖ１が目標値ＶＳＥＴ以上の場合（Ｖ１≧ＶＳＥＴ）、選択部４３は、第１電流指令値ｉｄ１＊，ｉｑ１＊を電流指令値ｉｄｃ＊，ｉｑｃ＊として選択する。出力電圧Ｖ１が目標値ＶＳＥＴ以上の場合とは、第２電流指令値が、等電圧楕円５４上又は等電圧楕円５４の外側にある場合である。
　このとき、第１電流指令値は、第２電流指令値よりｉｑｃ軸（ｑ軸）に近い。すなわち、第１電流指令値の電流位相角β１は、第２電流指令値の電流位相角β２よりｉｑｃ軸（ｑ軸）側である。また、第１電流指令値のｑ軸成分ｉｑ１＊は、第２電流指令値のｑ軸成分ｉｑ２＊より大きく、第１電流指令値のｄ軸成分ｉｄ１＊は、第２電流指令値のｄ軸成分ｉｄ２＊より小さい。
　選択部４３は、選択した第１電流指令値又は第２電流指令値を電流指令値として出力する。選択部４３により出力された電流指令値は、電圧指令算出部２５に入力される。
　次に、選択部４３による、出力電圧Ｖ１が目標値ＶＳＥＴより低いか否かの判定方法について説明する。
　選択部４３は、例えば、出力電圧Ｖ１を取得し、目標値ＶＳＥＴと比較することにより、判定を行う。選択部４３は、第１電流指令値毎の出力電圧が記憶されたテーブルを参照して出力電圧Ｖ１を選択してもよいし、以下の式により算出してもよい。
Figure JPOXMLDOC01-appb-M000005
　また、選択部４３は、第１電流指令値の電流位相角β１と、第２電流指令値の電流位相角β２と、を比較することにより、判定を行なってもよい。選択部４３は、電流位相角β２が電流位相角β１よりｑ軸側の場合、出力電圧Ｖ１が目標値ＶＳＥＴより低いと判定する。
　さらに、選択部４３は、第１電流指令値のｑ軸成分ｉｑ１＊と、第２電流指令値のｑ軸成分ｉｑ２＊と、を比較することにより、判定を行なってもよい。選択部４３は、ｑ軸成分ｉｑ１＊がｑ軸成分ｉｑ２＊より小さい場合、出力電圧Ｖ１が目標値ＶＳＥＴより低いと判定する。
　またさらに、選択部４３は、第１電流指令値のｄ軸成分ｉｄ１＊と、第２電流指令値のｄ軸成分ｉｄ２＊と、を比較することにより、判定を行なってもよい。選択部４３は、ｄ軸成分ｉｄ１＊がｄ軸成分ｉｄ２＊より大きい場合、出力電圧Ｖ１が目標値ＶＳＥＴより低いと判定する。
　電流指令値算出部２４がこのように電流指令値ｉｄｃ＊，ｉｑｃ＊を算出することにより、図６に実線で示すように、任意のトルク指令値Ｔ＊に対して、ＳｙｎＲＭ１の出力電圧を目標値ＶＳＥＴ以上とすることができる。
　従来の制御装置では、ＳｙｎＲＭ１が高速域で動作しても、トルク指令値Ｔ＊が小さい低負荷の状態では、図６に破線で示すように、ＳｙｎＲＭ１の出力電圧（誘起電圧）が小さくなった。このため、誘起電圧を利用して精度よく速度ωや回転位相角θを推定することが困難であった。
　しかしながら、本実施形態に係る制御装置２は、低負荷の場合であっても、出力電圧を目標値ＶＳＥＴ以上とすることができるため、ＳｙｎＲＭ１の出力電圧に基づいて、ＳｙｎＲＭ１の速度ωや回転位相角θを精度よく推定することができる。
　なお、以上の説明では、制御装置２がＳｙｎＲＭ１の動作を制御する場合について説明したが、この制御装置２は、磁石磁束が小さいＰＭＳＭの制御装置として利用することもできる。磁石磁束が小さいＰＭＳＭは、ＳｙｎＲＭ１と同様に、低負荷時の誘起電圧が小さい。このようなＰＭＳＭに本実施形態に係る制御装置２を適用することにより、ＰＭＳＭの速度ωや回転位相角θの低負荷時における推定精度を向上させることができる。
（第１実施形態の変形例１）
　本実施形態に係る電流指令値算出部２４の変形例１について、図７を参照して説明する。図７は、電流指令値算出部２４の変形例１の構成を示す図である。図７に示すように、第２算出部４３は、ｑ軸成分算出部４４と、ｄ軸成分算出部４５と、を備える。電流指令値算出部２４の他の構成は図４と同様である。
　ｑ軸成分算出部４４は、推定速度ωｅｓｔに基づいて第２電流指令値のｑ軸成分ｉｑ２＊を算出する。具体的には、ｑ軸成分算出部４４は、電流値ｉｑ２＊が推定速度ωｅｓｔと反比例するようにｑ軸成分ｉｑ２＊を算出する。ｑ軸成分ｉｑ２＊は、例えば、以下の式により算出される。
Figure JPOXMLDOC01-appb-M000006
　ｄ軸成分算出部４５は、トルク指令値Ｔ＊及びｑ軸成分ｉｑ２＊に基づいて、第２電流指令値のｄ軸成分ｉｄ２＊を算出する。ｄ軸成分ｉｄ２＊は、例えば、以下の式により算出される。
Figure JPOXMLDOC01-appb-M000007
　このような構成により、制御装置２は、低負荷時の出力電圧を目標値ＶＳＥＴ以上とすることができる。また、高速域における電流損失を抑制することができる。
（第１実施形態の変形例２）
　本実施形態に係る電流指令値算出部２４の変形例２について、図８～図１０を参照して説明する。図８は、電流指令値算出部２４の変形例２の構成を示す図である。図８に示すように、この電流指令値算出部２４は、トルク指令値Ｔ＊及び推定速度ωｅｓｔに基づいて、電流指令値ｉｑｃ＊，ｉｄｃ＊を算出する。
　本変形例において、電流指令値算出部２４は、図９に示すように、電流指令値のｑ軸成分ｉｑｃ＊が所定の目標値ＩｑＳＥＴ以上となるように、ｑ軸成分ｉｑｃ＊を算出する。また、電流指令値算出部２４は、図１０に示すように、推定速度ωｅｓｔが大きくなるほど、目標値ＩｑＳＥＴが小さくなるように、ｑ軸成分ｉｑｃ＊を生成する。
　このような構成により、低負荷時のモータ出力電圧を目標値ＶＳＥＴ以上とすることができる。また、高速域における電流損失を抑制することができる。
（第１実施形態の変形例３）
　本実施形態に係る電流指令値算出部２４の変形例３について、図１１及び図１２を参照して説明する。以上で説明した電流指令値算出部２４は、２つの電流指令値から一方を選択することにより、出力電圧が目標値ＶＳＥＴ以上となる電流指令値を算出した。これに対して、本変形例では、１つの電流指令値を補正することにより、出力電圧が目標値ＶＳＥＴ以上となる電流指令値を算出する。
　図１１は、電流指令値算出部２４の変形例３の構成を示す図である。図１１に示すように、電流指令値算出部２４は、ｑ軸成分算出部４６と、補正部４７と、加算器４８と、ｄ軸成分算出部４９と、を備える。
　ｑ軸成分算出部４６は、トルク指令値Ｔ＊に基づいて、電流指令値ｉｑ３＊を算出する。電流指令値ｉｑ３＊は、電流指令値のｑ軸成分の補正前の値である。本変形例における電流指令値算出部２４は、電流指令値ｉｑ３＊を補正することにより、電流指令値のｑ軸成分を算出する。電流指令値ｉｑ３＊は、ＳｙｎＲＭ１のトルクがトルク指令値Ｔ＊となるように、制御装置で実現したい制御に応じた任意の方法で算出される。
　例えば、ｑ軸成分算出部４６は、固定子に流す電流の大きさが最小となるように、電流指令値ｉｑ３＊を算出する。この場合、ｑ軸成分算出部４６は、磁気飽和を無視して、電流位相角βが１３５度となる電流指令値ｉｑ３＊を算出してもよいし、磁気飽和を考慮して、電流位相角βが１３５度より進んだ角度となる電流指令値ｉｑ３＊を算出してもよい。また、ｑ軸成分算出部４６は、ＳｙｎＲＭ１の効率や力率が最大となるように電流指令値ｉｑ３＊を算出してもよい。
　ｑ軸成分算出部４６が算出した電流指令値ｉｑ３＊は、加算器４８に入力される。
　補正部４７は、電圧指令値算出部２５から、電圧指令値ｖｄｃ＊，ｖｑｃ＊を入力される。補正部４７は、電圧指令値ｖｄｃ＊，ｖｑｃ＊に基づいて、補正値Δｉｑ＊を算出する。補正値Δｉｑ＊は、電流指令値ｉｑ３＊を補正して、電流指令値のｑ軸成分ｉｑｃ＊を算出するための電流指令値である。補正部４７の詳細については後述する。補正部４７が算出した補正値Δｉｑ＊は、加算器４８に入力される。
　加算器４８は、電流指令値ｉｑ３＊に補正値Δｉｑ＊を加算し、電流指令値ｉｑｃ＊を算出する。加算器４８が算出した電流指令値ｉｑｃ＊は、ｄ軸成分算出部４９に入力される。
　ｄ軸成分算出部４９は、トルク指令値Ｔ＊と、電流指令値ｉｑｃ＊と、に基づいて、ＳｙｎＲＭ１のトルクがトルク指令値Ｔ＊となるように、電流指令値のｄ軸成分ｉｄｃ＊を算出する。ｄ軸成分ｉｄｃ＊は、例えば、以下の式により算出される。
Figure JPOXMLDOC01-appb-M000008
　ここで、補正部４７について、図１２を参照して説明する。図１２は、補正部４７の構成を示す図である。図１２に示すように、補正部４７は、出力電圧算出部６１と、減算器６２と、リミッタ６３と、ＰＩ制御器６４と、を備える。
　出力電圧算出部６１は、電圧指令値算出部２５が算出した電圧指令値ｖｄｃ＊，ｖｑｃ＊を入力される。出力電圧算出部６１は、電圧指令値ｖｄｃ＊，ｖｑｃ＊に基づいて、ＳｙｎＲＭ１の出力電圧Ｖ１を算出する。出力電圧Ｖ１は、例えば、以下の式により算出される。
Figure JPOXMLDOC01-appb-M000009
　出力電圧算出部６１が算出した出力電圧Ｖ１は、減算器６２に入力される。
　減算器６２は、出力電圧の目標値ＶＳＥＴから出力電圧Ｖ１を減算し、出力電圧の誤差ΔＶを算出する。減算器６２が算出した誤差ΔＶは、リミッタ６３に入力される。
　リミッタ６３は、誤差ΔＶを０以上にリミットする。すなわち、０以上の誤差ΔＶだけを出力する。リミッタ６３が出力した誤差ΔＶは、ＰＩ制御器６４に入力される。
　ＰＩ制御器６４は、０以上にリミットされた誤差ΔＶに基づいてＰＩ制御を行い、ＳｙｎＲＭ１の出力電圧Ｖ１が目標値ＶＳＥＴ以上となる補正値Δｉｑ＊を算出する。ＰＩ制御器６４が算出した補正値Δｉｑ＊は、加算器４８に入力される。
　このような構成により、電流指令値に応じた出力電圧Ｖ１が目標値ＶＳＥＴより低い場合すなわち、ΔＶが０以上の場合、出力電圧Ｖ１が目標値ＶＳＥＴ以上となるように、電流指令値に補正値Δｉｑ＊がオフセット電流として追加される。これにより、低負荷時の出力電圧Ｖ１を目標値ＶＳＥＴ以上とすることができる。
（第１実施形態の変形例４）
　本実施形態に係る電流指令値算出部２４の変形例４について、図１６を参照して説明する。本変形例では、変形例３と同様に、１つの電流指令値を補正することにより、出力電圧が目標値ＶＳＥＴ以上となる電流指令値を算出する。変形例３では、補正部４７によって補正値Δｉｑ＊を生成し、これを加算器４８で電流指令値ｉｑ３＊に加算することで、電流指令値ｉｑｃ＊を生成した。これに対して、本変形例では、下限リミッタ９１によって、電流指令値ｉｑ３＊の下限をリミットすることで、電流指令値ｉｑｃ＊を生成する。
　下限リミッタ９１は、図９のように、ｉｑｃ＊がｉｑＳＥＴ以上になるように、ｉｑ３＊を補正（リミット）する。ｉｑＳＥＴは、例えば本実施例が使用される最小速度ωＭＩＮと、出力電圧の目標値ＶＳＥＴと、から次式で計算される。
Figure JPOXMLDOC01-appb-M000010
　このような構成により、低負荷時であっても出力電圧Ｖ１を目標値ＶＳＥＴ以上とすることができる。本変形例では、第１実施形態および変形例１～３に対して、下限リミッタ９１を追加するだけという簡単な構成で、出力電圧Ｖ１を目標値ＶＳＥＴ以上とすることができる。
（第２実施形態）
　次に、第２実施形態に係る制御装置２について、図１４～図１６を参照して説明する。図１４は、本実施形態に係るシステムの構成を示す図である。図１４に示すように、本実施形態に係る制御装置２は、トルク指令値算出部２９を更に備える。他の構成は、第１実施形態と同様である。
　トルク指令値算出部（速度制御器）２９は、速度指令値ω＊及び推定速度ωｅｓｔに基づいて、トルク指令値Ｔ＊を算出する。速度指令値ω＊とは、ＳｙｎＲＭ１の回転子を回転させる角速度である。本実施形態において、電流指令値算出部２４には、トルク指令値算出部２９が算出したトルク指令値Ｔ＊が入力される。なお、本実施形態では、速度指令値ω＊は、外部装置から入力されるものとする。
　図１５は、トルク指令値算出部２９の構成を示す図である。図１５に示すように、トルク指令値算出部２９は、減算器７１と、ＰＩ制御器７２と、を備える。
　減算器７１は、速度指令値ω＊と推定速度ωｅｓｔとを入力される。減算器７１は、速度指令値ω＊から推定速度ωｅｓｔを減算し、回転子の角速度の誤差Δωを算出する。減算器７１が算出した誤差Δωは、ＰＩ制御器７２に入力される。
　ＰＩ制御器７２は、誤差Δωに基づいてＰＩ制御を行い、誤差Δωが０となるトルク指令値Ｔ＊を算出する。ＰＩ制御器７２が算出したトルク指令値Ｔ＊は、電流指令値算出部２４に入力される。
　上述の通り、電流指令値算出部２４は、第１実施形態で説明した構成としてもよいが、他の構成とすることもできる。図１６は、電流指令値算出部２４の他の構成を示す図である。この電流指令値算出部２４は、トルク指令値Ｔ＊及び位相角指令値β＊に基づいて、電流指令値ｉｄｃ＊，ｉｄｃ＊を算出する。
　図１６に示すように、この電流指令値算出部２４は、出力電圧算出部８１と、減算器８２と、リミッタ８３と、ＰＩ制御器８４と、位相角指令値算出部８５と、減算器８６と、電流指令値取得部８７と、を備える。
　出力電圧算出部８１は、電圧指令値算出部２５が算出した電圧指令値ｖｄｃ＊，ｖｑｃ＊を入力される。出力電圧算出部８１は、電圧指令値ｖｄｃ＊，ｖｑｃ＊に基づいて、ＳｙｎＲＭ１の出力電圧Ｖ１を算出する。出力電圧算出部８１が算出した出力電圧Ｖ１は、減算器８２に入力される。
　減算器８２は、出力電圧の目標値ＶＳＥＴから出力電圧Ｖ１を減算し、出力電圧の誤差ΔＶを算出する。減算器８２が算出した誤差ΔＶは、リミッタ８３に入力される。
　リミッタ８３は、誤差ΔＶを０以上にリミットする。すなわち、０以上の誤差ΔＶだけを出力する。リミッタ８３が出力した誤差ΔＶは、ＰＩ制御器８４に入力される。
　ＰＩ制御器８４は、０以上にリミットされた誤差ΔＶに基づいてＰＩ制御を行い、ＳｙｎＲＭ１の出力電圧Ｖ１が目標値ＶＳＥＴ以上となる補正値Δβ＊を算出する。補正値Δβ＊は、位相角指令値β１＊を補正して位相角指令値β＊を算出するための位相角指令値である。ＰＩ制御器８４が算出した補正値Δβ＊は、減算器８６に入力される。
　位相角指令値算出部８５は、トルク指令値Ｔ＊に基づいて、位相角指令値β１＊を算出する。位相角指令値β１＊は、位相角指令値β＊の補正前の値である。位相角指令値β１＊は、ＳｙｎＲＭ１のトルクがトルク指令値Ｔ＊となるように、制御装置で実現したい制御に応じた任意の方法で算出される。
　例えば、位相角指令値算出部８５は、固定子に流す電流の大きさが最小となるように、位相角指令値β１＊を算出する。この場合、位相角指令値算出部８５は、磁気飽和を無視して、位相角指令値β１＊を１３５度としてもよいし、磁気飽和を考慮して、位相角指令値β１＊を１３５度より進んだ角度としてもよい。また、位相角指令値算出部８５は、ＳｙｎＲＭ１の効率や力率が最大となるように位相角指令値β１＊を算出してもよい。
　位相角指令値算出部８５が算出した位相角指令値β１＊は、減算器８６に入力される。
　減算器８６は、位相角指令値β１＊から補正値Δβ＊を減算して、位相角指令値β＊を算出する。減算器８６が算出した位相角指令値β＊は、電流指令値取得部８７に入力される。
　電流指令値取得部８７は、トルク指令値Ｔ＊及び位相角指令値β＊に基づいて、電流指令値ｉｄｃ＊，ｉｑｃ＊を算出する。具体的には、電流指令値取得部８７は、トルク指令値Ｔ＊の等トルク曲線上の、電流位相角が位相角指令値β＊となる点の座標を取得すればよい。電流指令値取得部８７により取得された電流指令値ｉｄｃ＊，ｉｑｃ＊は、電圧指令値算出部２５に入力される。
　以上説明した通り、図１６の電流指令値算出部２４は、出力電圧Ｖ１が目標値ＶＳＥＴ以上となるように、電流位相角を補正する。具体的には、電流指令値がｑ軸側に移動するように、電流位相角を遅れ方向に補正する。これにより、低負荷時の出力電圧Ｖ１を目標値ＶＳＥＴ以上とすることができる。
　また、本実施形態に係る制御装置２は、要求されるトルク制御の精度が第１実施形態に比べて低い。これは、本実施形態に係る制御装置２は、ＳｙｎＲＭ１の速度ωが速度指令値ω＊となるように制御するためである。
　このため、電流指令値算出部２４は、ＳｙｎＲＭ１のトルクを正確に制御するためのテーブル（電流指令値毎の出力電圧が記憶されたテーブルなど）を用いずに電流指令値を算出することができる。
　なお、本発明は上記各実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記各実施形態に開示されている複数の構成要素を適宜組み合わせることによって種々の発明を形成できる。また例えば、各実施形態に示される全構成要素からいくつかの構成要素を削除した構成も考えられる。さらに、異なる実施形態に記載した構成要素を適宜組み合わせてもよい。
１：シンクロナスリラクタンスモータ（ＳｙｎＲＭ）、２：インバータ制御装置、２１：インバータ、２２：電流検出器、２３：座標変換部、２４：電流指令値算出部、２５：電圧指令値算出部、２６：座標変換部、２７：ＰＷＭ変調器、２８：速度・回転位相角推定部、２９：トルク指令値算出部、３１：位相差δ設定部、３２：γ電圧算出部、３３：γ電圧推定部、３４：減算器、３５：ＰＩ制御器、３６：積分器、４１：第１算出部、４２：第２算出部、４３：選択部、４４：ｑ軸成分算出部、４５：ｄ軸成分算出部、４６：ｑ軸成分算出部、４７：補正部、４８：加算器、４９ｄ軸成分算出部、５１：第１電流ベクトル、５２：第２電流ベクトル、５３：等トルク曲線、５４：等電圧楕円、６１：出力電圧算出部、６２：減算器、６３：リミッタ、６４：ＰＩ制御器、７１：減算器、７２：ＰＩ制御器、８１：出力電圧算出部、８２減算器、８３：リミッタ、８４：ＰＩ制御器、８５：位相角指令値算出部、８６：減算器、８７：電流指令値取得部

Claims (15)

  1. 　所定の回転駆動対象と電気的に接続可能なインバータ主回路と、
    　前記インバータ主回路から出力される電流値を検出する電流検出器と、
    　前記インバータ主回路から出力される出力電圧が所定の目標値以上となる電流指令値を算出する電流指令値算出部と、
    　前記電流値が、前記電流指令値と等しくなる電圧指令値を算出する電圧指令値算出部と、
    　前記電圧指令値及び前記電流値に基づいて、前記回転駆動対象の推定回転位相角を算出する推定部と、
    を備えるインバータ制御装置。
  2. 　前記電流指令値算出部は、トルク指令値に基づいて第１電流指令値を算出し、前記トルク指令値及び推定速度に基づいて第２電流指令値を算出し、前記第１電流指令値に応じた前記出力電圧が前記目標値より低いか判定し、前記目標値以上のとき、前記第１電流指令値を前記電流指令値として出力し、前記目標値より低いとき、前記第２電流指令値を前記電流指令値として出力する
    請求項１に記載のインバータ制御装置。
  3. 　前記第２電流指令値に応じた前記出力電圧は前記目標値である
    請求項２に記載のインバータ制御装置。
  4. 　前記電流指令値算出部は、前記第１電流指令値及び前記第２電流指令値のうち、電流位相角がｑ軸に近い方を前記電流指令値として出力する
    請求項２又は請求項３に記載のインバータ制御装置。
  5. 　前記電流指令値算出部は、前記第１電流指令値及び前記第２電流指令値のうち、ｑ軸成分が大きい方を前記電流指令値として出力する
    請求項２乃至請求項４のいずれか１項に記載のインバータ制御装置。
  6. 　前記電流指令値算出部は、トルク指令値、推定速度、及び前記電圧指令値の少なくとも１つに基づいて、前記電流指令値を算出する
    請求項１乃至請求項５のいずれか１項に記載のインバータ制御装置。
  7. 　前記電流指令値算出部は、前記第１電流指令値に基づいて算出した前記出力電圧の電圧値と、前記目標値と、を比較することにより前記判定を行う
    請求項２乃至請求項６のいずれか１項に記載のインバータ制御装置。
  8. 　前記電流指令値算出部は、前記第１電流指令値のｑ軸成分と、前記第２電流指令値のｑ軸成分と、を比較することにより前記判定を行う
    請求項２乃至請求項６のいずれか１項に記載のインバータ制御装置。
  9. 　前記電流指令値算出部は、前記第１電流指令値のｄ軸成分と、前記第２電流指令値のｄ軸成分と、を比較することにより前記判定を行う
    請求項２乃至請求項６のいずれか１項に記載のインバータ制御装置。
  10. 　前記電流指令値算出部は、前記第１電流指令値の電流位相角と、前記第２電流指令値の電流位相角と、を比較することにより前記判定を行う
    請求項２乃至請求項６のいずれか１項に記載のインバータ制御装置。
  11. 　前記電流指令値算出部は、前記出力電圧が前記目標値より低いとき、前記出力電圧と前記目標値との誤差に基づいて補正値を算出し、前記電流指令値に前記補正値を加算する
    請求項１に記載のインバータ制御装置。
  12. 　速度指令値及び推定速度に基づいて、前記トルク指令値を算出するトルク指令値算出部を更に備える
    請求項１乃至請求項１１のいずれか１項に記載のインバータ制御装置。
  13. 前記電流指令値算出部は、前記電流指令値の下限をリミットする下限リミッタを備える
    請求項１乃至請求項１２のいずれか１項に記載のインバータ制御装置。
  14. 　前記下限リミッタは、最低速度における前記出力電圧が前記目標値となるように設定される
    請求項１３に記載のインバータ制御装置。
  15. 　モータと、
    　前記モータに接続されたインバータ主回路と、
    　前記インバータ主回路から出力される電流値を検出する電流検出器と、
    　前記インバータ主回路から出力される出力電圧が所定の目標値以上となる電流指令値を算出する電流指令値算出部と、
    　前記電流値が、前記電流指令値と等しくなる電圧指令値を算出する電圧指令値算出部と、
    　前記電圧指令値及び前記電流値に基づいて、前記モータの推定回転位相角を算出する推定部と、
    を備えるモータ駆動システム。
PCT/JP2016/052161 2015-01-28 2016-01-26 インバータ制御装置及びモータ駆動システム WO2016121751A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
EP16743354.9A EP3252941B1 (en) 2015-01-28 2016-01-26 Inverter control apparatus and motor driving system
CN201680003392.4A CN107078674B (zh) 2015-01-28 2016-01-26 逆变器控制装置以及电机驱动系统
JP2016535264A JP6367332B2 (ja) 2015-01-28 2016-01-26 インバータ制御装置及びモータ駆動システム
US15/636,198 US10158305B2 (en) 2015-01-28 2017-06-28 Inverter controller and motor driving system

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2015014722 2015-01-28
JP2015-014722 2015-01-28

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US15/636,198 Continuation US10158305B2 (en) 2015-01-28 2017-06-28 Inverter controller and motor driving system

Publications (1)

Publication Number Publication Date
WO2016121751A1 true WO2016121751A1 (ja) 2016-08-04

Family

ID=56543362

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2016/052161 WO2016121751A1 (ja) 2015-01-28 2016-01-26 インバータ制御装置及びモータ駆動システム

Country Status (5)

Country Link
US (1) US10158305B2 (ja)
EP (1) EP3252941B1 (ja)
JP (1) JP6367332B2 (ja)
CN (1) CN107078674B (ja)
WO (1) WO2016121751A1 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018038111A1 (ja) * 2016-08-22 2018-03-01 株式会社 東芝 インバータ制御装置およびドライブシステム
WO2018123426A1 (ja) * 2016-12-28 2018-07-05 アルプス電気株式会社 直流整流子電動機の回転に関する情報を取得する装置及び方法
KR20190032556A (ko) * 2016-09-05 2019-03-27 가부시끼가이샤 도시바 인버터 제어 장치 및 모터 구동 시스템
EP3509210A4 (en) * 2016-09-05 2020-04-22 Toshiba Infrastructure Systems & Solutions Corporation CONVERTER CONTROL DEVICE AND ELECTRIC MOTOR DRIVE SYSTEM
US11296633B2 (en) 2018-01-12 2022-04-05 Mitsubishi Electric Corporation Rotary machine control device

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10763740B2 (en) 2016-04-15 2020-09-01 Emerson Climate Technologies, Inc. Switch off time control systems and methods
US10277115B2 (en) 2016-04-15 2019-04-30 Emerson Climate Technologies, Inc. Filtering systems and methods for voltage control
US10312798B2 (en) 2016-04-15 2019-06-04 Emerson Electric Co. Power factor correction circuits and methods including partial power factor correction operation for boost and buck power converters
US10305373B2 (en) 2016-04-15 2019-05-28 Emerson Climate Technologies, Inc. Input reference signal generation systems and methods
US10284132B2 (en) 2016-04-15 2019-05-07 Emerson Climate Technologies, Inc. Driver for high-frequency switching voltage converters
US9933842B2 (en) 2016-04-15 2018-04-03 Emerson Climate Technologies, Inc. Microcontroller architecture for power factor correction converter
US10656026B2 (en) 2016-04-15 2020-05-19 Emerson Climate Technologies, Inc. Temperature sensing circuit for transmitting data across isolation barrier
DK3522362T3 (da) * 2018-02-01 2024-03-04 Siemens Gamesa Renewable Energy As Styring af en elektrisk permanentmagnetmaskine med multiviklingssæt
JP7032250B2 (ja) * 2018-06-28 2022-03-08 株式会社日立産機システム 電力変換装置
CN112039401A (zh) * 2019-05-17 2020-12-04 杭州三花研究院有限公司 一种控制方法和控制装置及其系统
CN113839605B (zh) * 2020-06-23 2024-03-15 苏州宝时得电动工具有限公司 电机转速控制方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000037099A (ja) * 1998-07-16 2000-02-02 Yaskawa Electric Corp 交流電動機の可変速制御装置
JP2006121768A (ja) * 2004-10-19 2006-05-11 Fuji Electric Fa Components & Systems Co Ltd 誘導電動機の可変速制御装置
JP2009261103A (ja) * 2008-04-15 2009-11-05 Jtekt Corp モータ制御装置
JP2010130751A (ja) * 2008-11-26 2010-06-10 Honda Motor Co Ltd 電動機の相電流推定装置および電動機の磁極位置推定装置
JP2010136586A (ja) * 2008-12-08 2010-06-17 Honda Motor Co Ltd 電動機の磁極位置推定装置

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100354775B1 (ko) * 2000-03-25 2002-11-04 엘지전자 주식회사 동기 릴럭턴스 모터의 속도 제어장치
KR100421373B1 (ko) * 2001-06-20 2004-03-06 엘지전자 주식회사 동기 릴럭턴스 모터의 회전 속도 제어장치
US6850030B2 (en) * 2001-07-04 2005-02-01 Kabushiki Kaisha Yaskawa Denki Method and device for controlling currents of synchronous motor
JP3692085B2 (ja) 2002-02-21 2005-09-07 株式会社東芝 モータ制御方法及び装置
US6822417B2 (en) * 2002-03-22 2004-11-23 Matsushita Electric Industrial Co., Ltd. Synchronous reluctance motor control device
JP4579627B2 (ja) 2004-09-02 2010-11-10 三菱電機株式会社 回転機の制御装置
JP2010119284A (ja) * 2008-10-16 2010-05-27 Panasonic Corp モータ駆動装置
EP2360831B1 (en) * 2008-12-15 2019-05-29 Mitsubishi Electric Corporation Power converter for driving motor
JP5971707B2 (ja) * 2011-08-29 2016-08-17 株式会社東芝 同期電動機のセンサレス制御装置ならびにインバータ装置
US9438153B2 (en) * 2013-03-28 2016-09-06 Aisin Aw Co., Ltd. Rotary electric machine control device
JP6193006B2 (ja) * 2013-06-20 2017-09-06 株式会社東芝 電気車制御装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000037099A (ja) * 1998-07-16 2000-02-02 Yaskawa Electric Corp 交流電動機の可変速制御装置
JP2006121768A (ja) * 2004-10-19 2006-05-11 Fuji Electric Fa Components & Systems Co Ltd 誘導電動機の可変速制御装置
JP2009261103A (ja) * 2008-04-15 2009-11-05 Jtekt Corp モータ制御装置
JP2010130751A (ja) * 2008-11-26 2010-06-10 Honda Motor Co Ltd 電動機の相電流推定装置および電動機の磁極位置推定装置
JP2010136586A (ja) * 2008-12-08 2010-06-17 Honda Motor Co Ltd 電動機の磁極位置推定装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP3252941A4 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018038111A1 (ja) * 2016-08-22 2018-03-01 株式会社 東芝 インバータ制御装置およびドライブシステム
TWI668953B (zh) * 2016-08-22 2019-08-11 日商東芝股份有限公司 Inverter control device and drive system
US10637381B2 (en) 2016-08-22 2020-04-28 Kabushiki Kaisha Toshiba Inverter control device and drive system
KR20190032556A (ko) * 2016-09-05 2019-03-27 가부시끼가이샤 도시바 인버터 제어 장치 및 모터 구동 시스템
EP3509211A4 (en) * 2016-09-05 2020-04-15 Kabushiki Kaisha Toshiba CONVERTER CONTROL DEVICE AND DRIVE SYSTEM
EP3509210A4 (en) * 2016-09-05 2020-04-22 Toshiba Infrastructure Systems & Solutions Corporation CONVERTER CONTROL DEVICE AND ELECTRIC MOTOR DRIVE SYSTEM
US10742151B2 (en) 2016-09-05 2020-08-11 Toshiba Infrastructure Systems & Solutions Corporation Inverter control device and motor drive system
US10833613B2 (en) 2016-09-05 2020-11-10 Kabushiki Kaisha Toshiba Inverter control apparatus and motor drive system
KR102285041B1 (ko) * 2016-09-05 2021-08-04 가부시끼가이샤 도시바 인버터 제어 장치 및 모터 구동 시스템
WO2018123426A1 (ja) * 2016-12-28 2018-07-05 アルプス電気株式会社 直流整流子電動機の回転に関する情報を取得する装置及び方法
US11296633B2 (en) 2018-01-12 2022-04-05 Mitsubishi Electric Corporation Rotary machine control device

Also Published As

Publication number Publication date
CN107078674A (zh) 2017-08-18
EP3252941B1 (en) 2021-06-30
EP3252941A1 (en) 2017-12-06
US20170317623A1 (en) 2017-11-02
JP6367332B2 (ja) 2018-08-01
CN107078674B (zh) 2019-09-17
US10158305B2 (en) 2018-12-18
EP3252941A4 (en) 2018-10-03
JPWO2016121751A1 (ja) 2017-04-27

Similar Documents

Publication Publication Date Title
JP6367332B2 (ja) インバータ制御装置及びモータ駆動システム
JP5494760B2 (ja) 電動機制御装置
JP4989075B2 (ja) 電動機駆動制御装置及び電動機駆動システム
WO2016121237A1 (ja) インバータ制御装置及びモータ駆動システム
TWI654827B (zh) 換流器控制裝置及馬達驅動系統
JP6776066B2 (ja) インバータ制御装置および電動機駆動システム
EP2424105A2 (en) Vector control apparatus and motor control system
JP2002095300A (ja) 永久磁石同期電動機の制御方法
JP2009095135A (ja) 同期電動機の制御装置
JP2009060688A (ja) 同期電動機の制御装置
JP5499965B2 (ja) 交流回転機の制御装置
JP2006197712A (ja) 同期電動機の駆動システム及び同期電動機の駆動方法
JP5648310B2 (ja) 同期モータの制御装置、及び同期モータの制御方法
JP2000037098A (ja) 速度センサレスベクトル制御を用いた電力変換装置
JP5660191B2 (ja) 電動機制御装置
JP2008148437A (ja) 永久磁石型同期モータの制御装置
JP3692085B2 (ja) モータ制御方法及び装置
KR102409792B1 (ko) 영구 자석 동기 전동기의 제어 장치, 마이크로 컴퓨터, 전동기 시스템 및 영구 자석 동기 전동기의 운전 방법
JP7251424B2 (ja) インバータ装置及びインバータ装置の制御方法
JP5390970B2 (ja) モータ制御装置
JP2020036513A (ja) モータ制御装置
JP7099226B2 (ja) モータ制御装置
JP7226211B2 (ja) インバータ装置及びインバータ装置の制御方法
JP5456873B1 (ja) 同期機制御装置
JP7012931B2 (ja) インバータ制御装置および電動機駆動システム

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref document number: 2016535264

Country of ref document: JP

Kind code of ref document: A

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 16743354

Country of ref document: EP

Kind code of ref document: A1

REEP Request for entry into the european phase

Ref document number: 2016743354

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: DE