WO2019138691A1 - 回転機の制御装置 - Google Patents

回転機の制御装置 Download PDF

Info

Publication number
WO2019138691A1
WO2019138691A1 PCT/JP2018/042568 JP2018042568W WO2019138691A1 WO 2019138691 A1 WO2019138691 A1 WO 2019138691A1 JP 2018042568 W JP2018042568 W JP 2018042568W WO 2019138691 A1 WO2019138691 A1 WO 2019138691A1
Authority
WO
WIPO (PCT)
Prior art keywords
rotating machine
rotor position
component
estimated
inductance
Prior art date
Application number
PCT/JP2018/042568
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 DE112018006843.9T priority Critical patent/DE112018006843T5/de
Priority to JP2019564320A priority patent/JP6847268B2/ja
Priority to US16/959,446 priority patent/US11296633B2/en
Publication of WO2019138691A1 publication Critical patent/WO2019138691A1/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
    • 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
    • H02P6/185Circuit arrangements for detecting position without separate position detecting elements using inductance sensing, e.g. pulse excitation
    • 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/13Observer control, e.g. using Luenberger observers or Kalman filters
    • 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
    • H02P25/00Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02Arrangements 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/08Reluctance motors
    • H02P25/086Commutation
    • H02P25/089Sensorless control

Definitions

  • the present invention relates to a control device for a rotating machine that obtains and controls rotor position information without using a position sensor that detects a rotor position, and a rotating machine having an inductance alternating current component whose inductance changes according to the rotor position.
  • the position sensorless control method of the rotating machine includes a method of estimating a rotor position by applying a high frequency voltage to the rotating machine, a rotor position from an induced voltage of the rotating machine, a linkage flux, etc. without applying a high frequency voltage. There is a method to estimate
  • Patent Document 1 discloses a method of estimating a rotor position by applying a high frequency voltage to a rotating machine.
  • the rotating machine current when a high frequency voltage is applied to the rotating machine is detected, and a high frequency current having the same frequency component as the high frequency voltage is extracted.
  • the rotor position is estimated using the inductance of the rotating machine, that is, the fact that the amplitude of the high frequency current changes at a frequency twice that of the rotor position electrical angle.
  • torque ripples or noise is generated by the superimposed high frequency voltage.
  • the maximum torque of the rotating machine or the speed of the rotating machine is increased by the amount of high frequency voltage superimposed on the voltage applied to the winding of the rotating machine and the current flowing in the winding of the rotating machine.
  • the rotational speed decreases or both decrease.
  • Patent Documents 2, 3 and 4 disclose methods of estimating a rotor position without applying a high frequency voltage.
  • the product of the inductance of the rotor q-axis or d-axis and the motor current is subtracted from the linkage flux of the rotary machine,
  • the components that rotate in synchronization with the rotor position are extracted.
  • the components include ActiveFlux based on d axis and ActiveFlux based on q axis.
  • the direction of the rotor in which the inductance is maximized is taken as d axis
  • the direction in which the inductance is minimized is taken as q axis.
  • the d axis is an axis called a magnetic flux axis
  • the q axis is an axis called a torque axis.
  • a rotor position is estimated by estimating ActiveFlux with an observer.
  • a rotor position is estimated using an induced voltage generated by ActiveFlux.
  • the present invention has been made in view of the above, and it is an object of the present invention to provide a control device of a rotating machine capable of estimating a rotor position with a simple configuration regardless of the current flow direction of the rotating machine.
  • a control device for a rotating machine is a control device for a rotating machine that controls a rotating machine having an inductance AC component whose inductance changes according to a rotor position. And a current detector for detecting a rotating machine current flowing to the rotating machine.
  • a controller for a rotating machine includes a position estimator for estimating a rotor position, a controller for outputting a rotating machine voltage command for driving the rotating machine based on the rotating machine current and the rotor position, and a rotating machine voltage And a voltage applying device for applying a voltage to the rotating machine based on the command.
  • the position estimator is characterized by estimating the rotor position from the flux linkage inductance alternating current generated by the inductance alternating current component and the rotating machine current.
  • a diagram showing a configuration of a control device for a rotating machine according to Embodiment 1 of the present invention The figure which shows the structure of the position estimator shown in FIG.
  • the figure which shows the structure of the position estimator shown in FIG. The figure which shows the structure of the control apparatus of the rotary machine which concerns on Embodiment 3 of this invention.
  • the figure which shows the structure of the position estimator shown in FIG. The figure which shows the structure of the control coordinate angle calculator shown in FIG.
  • FIG. 1 A diagram showing the configuration of the adaptive estimator shown in FIG.
  • the figure which shows the structure of the position estimator shown in FIG. A diagram showing a first hardware configuration example of a control device for a rotary electric machine according to Embodiments 1 to 5 of the present invention
  • FIG. 1 is a view showing the configuration of a control device for a rotary machine according to Embodiment 1 of the present invention.
  • Control device 100 for a rotating machine according to the first embodiment applies a voltage application device to apply an AC voltage to rotating machine 1 according to rotating machine voltage commands v su * , v sv * , v sw * for driving the rotating machine.
  • a current detector 2 that detects an alternating current supplied from the voltage application unit 5 to the rotating machine 1 and outputs the detected alternating current as a rotating machine current isu , isv , isw .
  • control device 100 of the rotating machine uses the rotating machine voltage commands v su * , v sv * , v sw * and the rotating machine current i su , isv , isw to estimate an estimated rotor position ⁇ which is an electrical angle.
  • a position estimator 3 for calculating ⁇ r and a controller 4 are provided.
  • the estimated rotor position ⁇ ⁇ r is an estimated value of a rotor position which is a rotational position of a rotor (not shown) constituting the rotary machine 1.
  • the estimated rotor position ⁇ ⁇ r is expressed by an electrical angle.
  • the rotating machine 1 controlled by the controller 100 of the rotating machine is a synchronous reluctance motor.
  • a synchronous reluctance motor is a motor having an inductance AC component whose inductance changes depending on the rotor position.
  • the inductance is a value obtained by dividing the flux linkage of the stator winding by the stator current value.
  • the rotating machine 1 may be simply referred to as a "rotating machine”.
  • the voltage applicator 5 is connected to the rotating machine via three wires of U, V and W.
  • the current detector 2 detects an alternating current flowing through the wiring.
  • the direction of the rotor in which the inductance is maximized is taken as d axis
  • the direction in which the inductance is minimized is taken as q axis
  • the rotor position is based on the d axis of the rotor.
  • the controller 4 includes a current command calculator 401, a three-phase current command calculator 402, and a three-phase current controller 403, and a rotary machine voltage command v for the rotary machine to generate an output corresponding to the torque command value T *. Calculate su * , v sv * , v sw * .
  • the current command calculator 401 calculates the current effective value command I ph * and the energization angle command ⁇ i * necessary for the rotary machine to generate an output corresponding to the torque command value T * .
  • the energization angle command ⁇ i * is a command of the energization angle indicating the angle formed by the actual rotor position and the rotating machine current.
  • the current effective value command I ph * and the energization angle command ⁇ i * are selected such that the current effective value with respect to torque, that is, the copper loss of the rotating machine is minimized.
  • the three-phase current command computing unit 402 uses the current effective value command I ph * , the energization angle command ⁇ i *, and the estimated rotor position ⁇ ⁇ r to generate a three-phase current command as shown in the following equation (1) Calculate i su * , i sv * , i sw * .
  • the three-phase current controller 403 calculates rotating machine voltage commands v su * , v sv * , v sw * .
  • the rotating machine voltage commands v su * , v sv * , v sw * are the rotating machine currents i su , isv , i sw on the three-phase coordinates, and the three-phase current commands i su * , isv * , isw *
  • the rotation machine current isu , isv , isw can be obtained, for example, by controlling Proportional Integral (PI).
  • FIG. 2 is a diagram showing the configuration of the position estimator shown in FIG.
  • the position estimator 3 includes a three-phase to two-phase converter 301, a three-phase to two-phase converter 302, a flux linkage alternating current calculator 303, a linkage flux inductance alternating current component estimator 304, A rotor position estimation error calculator 305 and a rotor position calculator 306 are provided.
  • the linkage flux inductance AC component calculator 303 is referred to as "AC component calculator 303”
  • the linkage flux inductance AC component estimator 304 is referred to as “AC component estimator 304”
  • the rotor position estimation error calculator 305 may be referred to as “estimated error calculator 305”.
  • Three-phase to two-phase converter 301 converts rotating machine voltage commands v su * , v sv * , v sw * on three-phase coordinates into rotating machine voltage commands v s ⁇ * , v s ⁇ * on two-phase coordinates .
  • v s ⁇ *, v s ⁇ * a transformation matrix C 32 of the following equation (2) is used for three-phase to two-phase transformation.
  • three-phase - two-phase converter 302 the rotating machine current i su on the three-phase coordinate, i sv, the rotating machine current i sa on a two-phase coordinate i sw, converted to i S [beta.
  • i s ⁇ , i s ⁇ the rotating machine current i su on the three-phase coordinate
  • i sa the rotating machine current i sa on a two-phase coordinate i sw
  • the flux linkage ⁇ s ⁇ of the rotating machine can be obtained by the following equation (3).
  • the cutoff frequency is sufficiently low with respect to the fundamental frequency component of the flux linkage s s ⁇ of the rotating machine.
  • the transfer function of the HPF is expressed by the following equation (5), where the cutoff frequency is ⁇ hpf .
  • L sdc in the equation (8) is an inductance direct current component which does not change depending on the rotor position
  • L mac is an inductance alternating current component which changes depending on the rotor position
  • L sac ⁇ represents an inductance AC component.
  • the change in inductance is generally represented by a sine or cosine function of 2 ⁇ r , where ⁇ r is a rotor position represented by an electrical angle.
  • the flux linkage generated by the rotating machine current i s .alpha..beta inductance AC component, and flux linkage inductance AC component.
  • the AC component calculator 303 calculates the flux linkage inductance AC division ⁇ sac, calc ⁇ according to the following equation (9) using the flux linkage ⁇ shpf ⁇ calculated according to the above equation (7).
  • the linkage flux inductance alternating current ratio, sac, calc ⁇ calculated by the above equation (9) is a linkage flux inductance alternating current calculation value calculated by the control device 100 of the rotary machine according to the first embodiment.
  • AC component estimator 304 uses a rotary motor current i s .alpha..beta and the estimated rotor position theta ⁇ r, estimate the flux linkage inductance AC component.
  • the ⁇ ⁇ sac ⁇ calculated by the above equation (10) is the flux linkage inductance AC component estimated value calculated by the control device 100 of the rotary machine according to the first embodiment.
  • “L ⁇ sac ⁇ ” in the above equation (10) is an “inductance alternating current estimated value”.
  • the estimation error computing unit 305 computes an estimation error of the rotor position using the flux linkage inductance alternating component calculation value and the linkage flux inductance alternating component estimation value.
  • the outer product of the Flux Inductance AC component calculation value and the Flux Inductance AC component estimated value is the value obtained by computing the second term of the above equation (9), that is, the above equation (8), and the above equation (10) Using, it is represented by the following (11) formula.
  • the estimated error of the rotor position can be calculated by the following equation (12).
  • the estimation error calculator 305 outputs the calculated estimation error of the rotor position to the rotor position calculator 306.
  • the rotor position calculator 306 PI-controls the estimation error of the rotor position and then integrates it to calculate an estimated rotor position ⁇ ⁇ r .
  • I ph is a current effective value
  • ⁇ i is a conduction angle, that is, an angle difference between the rotor position and the rotating machine current.
  • the dq axis currents i sd and isq are obtained by rotational coordinate conversion of the ⁇ - ⁇ axis current using the rotor position ⁇ r as shown in the following equation (16).
  • the dq axis currents i sd and i sq may be referred to as “rotary machine current i s dq ”.
  • the transform matrix C dq of the above equation (16) is used for the rotational coordinate transformation of the first embodiment. Since the direction of ActiveFlux in the above equation (15) is the rotor d-axis direction, the rotor position can be calculated by estimating this using an observer or the like. In the technique disclosed in Patent Document 3 described above, the rotor position is estimated using an induced voltage generated by ActiveFlux based on d-axis. Further, in Patent Document 2 mentioned above, as shown in the following equation (17), also regarding ActiveFlux based on q axis, which is obtained by subtracting the product of d axis inductance of the rotor and rotor current from the linkage flux of the rotor. Have been described.
  • ActiveFlux based on the q-axis is in the rotor q-axis direction
  • the rotor position can be estimated by estimating this using an observer or the like.
  • the position sensorless control method disclosed in the above-mentioned Patent Document 4 uses the d-axis component of the rotating machine current or flux linkage and the q-axis in ActiveFlux both in the d-axis reference and the q-axis reference. This problem is solved by weighting based on the components and using it to estimate the rotor position.
  • weighting based on the d-axis component and the q-axis component of the rotating machine current or linkage flux complicates not only control design of rotor position estimation but also control processing.
  • the control device 100 of the rotary machine according to the first embodiment estimates the rotor position using the flux linkage alternating current component, that is, the second term of the equation (8) and the equation (10).
  • the flux linkage AC component has components of both two phases, and by using this flux linkage AC component, it is possible to estimate the rotor position without depending on the current flow direction. Further, in the control device 100 for a rotating machine according to the first embodiment, it is not necessary to perform weighting based on the direction of energization. Therefore, the control device 100 of the rotating machine according to the first embodiment has a remarkable effect which can not estimate the rotor position with a simple configuration regardless of the current flow direction of the rotating machine.
  • FIG. 3 is a diagram showing the configuration of a control device for a rotary machine according to Embodiment 2 of the present invention.
  • the change in inductance of the rotating machine is represented by a sine function or cosine function of 2 [Theta] r as rotor position theta r.
  • the actual inductance change is precisely not the sine function or cosine function of 2 [Theta] r, an error occurs in the estimated rotor position.
  • control device 100 of the rotating machine according to the first embodiment requires a sine function and a cosine function of 2 ⁇ r to estimate the flux linkage alternating current component of equation (10) above, and the control operation amount is large.
  • Control device 100A for a rotary machine according to Embodiment 2 is configured to be able to accurately estimate the rotor position with a small amount of control calculation.
  • Control device 100A for a rotary machine according to Embodiment 2 includes position estimator 6 and controller 7 instead of position estimator 3 and controller 4 shown in FIG.
  • the other configuration is the same as or equivalent to the configuration of the first embodiment, and the same or equivalent components are denoted by the same reference numerals and redundant description will be omitted.
  • the controller 7 includes a current command calculator 701, a dq current controller 702, a rotational coordinate inverse converter 703, a two-phase to three-phase converter 704, a three-phase to two-phase converter 705, and a rotational coordinate converter 706. provided, the rotating machine voltage command v for the rotating machine to generate an output corresponding to the torque command value T * su *, v sv * , v computes a sw *.
  • the current command calculator 701 calculates current commands i sd * , i sq * on the rotational coordinates necessary for the rotating machine to generate an output corresponding to the torque command value T * .
  • the current commands i sd * and i sq * on the rotational coordinates are selected so as to minimize the current effective value with respect to torque, that is, the copper loss of the rotating machine.
  • the dq current controller 702 performs control such that the dq axis currents i sd and i sq subjected to rotational coordinate conversion by the rotational coordinate converter 706 become the current commands i sd * and i sq * , and the rotational coordinate is The upper rotary machine voltage command v sd * , v sq * is calculated.
  • proportional integral control or the like is used for this current control.
  • the rotational coordinate inverse converter 703 uses the estimated rotor position ⁇ ⁇ r calculated by the position estimator 6 as shown in the following equation (19), to set the rotating machine voltage command v sd * , v on the rotational coordinate. sq * the rotating machine voltage command v on the two-phase coordinate sa *, to rotational coordinate reverse conversion into v S [beta *.
  • the transformation matrix C dq ⁇ 1 ( ⁇ ⁇ r ) of the following equation (19) is used for rotational coordinate inverse transformation.
  • the two-phase to three-phase converter 704 converts a rotating machine voltage command v s ⁇ * , v s ⁇ * on two phase coordinates into a rotating machine voltage command v su * , on three phase coordinates Convert to v sv * , v sw * .
  • a transformation matrix C 23 of the following equation (20) is used for two-phase to three-phase transformation.
  • Three-phase - two-phase converter 705 uses the estimated rotor position ⁇ ⁇ r to perform rotational coordinate conversion from the two-phase coordinate rotating machine current i s ⁇ , i s ⁇ to the dq axis current i sd on the rotational coordinate, i sq Do.
  • FIG. 4 is a diagram showing the configuration of the position estimator shown in FIG.
  • the position estimator 6 includes a three-phase to two-phase converter 601, a three-phase to two-phase converter 602, a rotational coordinate converter 603, an interlinkage flux inductance alternating current component calculator 604, an interlinkage flux inductance An AC component estimator 605, a rotor position estimation error calculator 606, and a rotor position calculator 607 are provided.
  • the operations of the three-phase to two-phase converter 602, the linkage flux inductance alternating current component calculator 604, and the linkage flux inductance alternating current component estimator 605 are different from those of the first embodiment.
  • the linkage flux inductance AC component calculator 604 is referred to as an "AC component calculator 604"
  • the linkage flux inductance AC component estimator 605 is referred to as an “AC component estimator 605".
  • the 606 may be referred to as an “estimated error computing unit 606”.
  • the three-phase to two-phase converter 601 has two-phase rotating machine voltage commands v su * , v sv * , v sw * on three-phase coordinates. It converts into rotary machine voltage command v s ⁇ * , v s ⁇ * on the coordinates. Similar to the three-phase to two-phase converter 302 of the first embodiment, the three-phase to two-phase converter 602 rotates the rotating machine current isu , isv , isw on three-phase coordinates on two-phase coordinates. Convert to machine current i s ⁇ , i s ⁇ .
  • the flux linkage ⁇ s ⁇ of the rotating machine is expressed by the above equation (8) on the two-phase coordinates, but the flux linkage s s dq of the rotating machine shown by the following equation (21) is the control coordinate angle ⁇ s Is used to carry out rotational coordinate transformation of the flux linkage ⁇ s ⁇ of the rotating machine.
  • the second term of the following equation (21) is a flux linkage alternating current component generated by the inductance alternating current component and the rotating machine current.
  • the AC component calculator 604 calculates the flux linkage ⁇ shpf ⁇ in the two-phase coordinates using the equation (7) as in the first embodiment. Further, the AC component calculator 604 performs rotational coordinate conversion of the flux linkage ⁇ shpf ⁇ to be a linkage flux ⁇ shpf dq .
  • the linkage flux inductance alternating current ratio, sac, calc dq on the rotational coordinate is calculated by the following equation (22) according to the equation (21).
  • the linkage flux inductance alternating current ratio ⁇ sac, calc dq calculated by the following equation (22) is hereinafter referred to as a linkage flux inductance alternating current calculation value.
  • the rotational coordinate converter 603 uses the estimated rotor position ⁇ ⁇ r to convert the rotating machine current i s ⁇ , i s ⁇ on two-phase coordinates into the dq axis current i sd , i sq on the rotating coordinate. Do.
  • the AC component estimator 605 estimates the flux linkage inductance AC component using the rotating machine current i s dq , as shown in the following equation (23).
  • [Psi ⁇ sac dq above (24) is a flux linkage inductance AC component estimation value calculated by the control unit 100A of the rotating machine according to the second embodiment.
  • the estimation error computing unit 606 uses the calculated value of the linkage flux inductance alternating current component and the estimated value of the linkage flux inductance alternating current component to estimate the estimation error of the rotor position.
  • the rotor position calculator 607 performs PI control on the estimation error of the rotor position and then integrates the calculated error to calculate an estimated rotor position ⁇ ⁇ r .
  • the control device 100A of the rotating machine according to the second embodiment calculates the flux linkage alternating current component as shown in the above equation (22) on the rotational coordinates synchronized with the rotor position, and further calculates the above equation (24) As shown, the flux linkage alternating current component is estimated.
  • the flux linkage AC component has components of both two phases, and by using this flux linkage AC component, it is possible to estimate the rotor position without depending on the current flow direction. Further, in control apparatus 100A for a rotating machine according to Embodiment 2, it is not necessary to perform weighting according to the current-carrying direction.
  • FIG. 5 is a diagram showing the configuration of a control device for a rotary machine according to Embodiment 3 of the present invention.
  • the control device 100B for a rotary machine according to the third embodiment is configured to be able to estimate the rotational speed and rotational position of the rotor with higher response and stability.
  • the control device 100B for the rotating machine includes a position estimator 8 and a controller 10 instead of the position estimator 6 and the controller 7 according to the second embodiment shown in FIG.
  • the control device 100B of the rotating machine includes a control coordinate angle calculator 9.
  • the control coordinate angle calculator 9 may be referred to as an "angle calculator 9".
  • the other configuration is the same as or equivalent to the configuration of the second embodiment, and the same or equivalent components are denoted by the same reference numerals and redundant description will be omitted.
  • the controller 10 includes a speed calculator 1001, a current command calculator 701, a dq current controller 702, a rotational coordinate inverse converter 1002, a two-phase to three-phase converter 704, a three-phase to two-phase converter 705, and a rotation.
  • a coordinate converter 1003 is provided, and a rotating machine voltage command v su * , v sv * , v sw * is calculated for generating a rotational speed corresponding to the rotational speed command value ⁇ r * .
  • the speed calculator 1001 performs, for example, proportional integral control such that the estimated rotation speed ⁇ ⁇ r becomes the rotation speed command value ⁇ r * , and calculates the torque command value T * .
  • the rotary coordinate inverse converter 1002 uses the control coordinate angle ⁇ s calculated by the angle calculator 9 to rotate the rotary machine voltage command v sd * , v sq * on the rotary coordinate, and the rotary machine voltage command on the two-phase coordinate Inverse rotation transformation of v s ⁇ * and v s ⁇ * .
  • the rotational coordinate converter 1003 uses the control coordinate angle ⁇ s calculated by the angle calculator 9 to convert the rotating machine current i s ⁇ , i s ⁇ on the two-phase coordinate into the dq axis current i sd on the rotational coordinate. Convert rotational coordinates to i sq .
  • the position estimator 8 calculates the estimated rotational speed ⁇ ⁇ r and the estimated rotor position ⁇ ⁇ r .
  • FIG. 6 is a diagram showing the configuration of the position estimator shown in FIG.
  • the position estimator 8 includes a three-phase to two-phase converter 601, a three-phase to two-phase converter 602, a rotational coordinate converter 801, an interlinkage flux inductance alternating current component calculator 802, an interlinkage flux inductance An AC component estimator 803, an estimation error calculator 606, and a rotor position calculator 804 are provided.
  • the linkage flux inductance AC component calculator 802 may be referred to as an “AC component calculator 802”
  • the linkage flux inductance AC component estimator 803 may be referred to as an “AC component estimator 803”.
  • the AC component calculator 802 calculates the flux linkage inductance alternating current ratio ⁇ sac, calc dq on the rotational coordinates, that is, the flux linkage component AC component computation value.
  • the control coordinate angle ⁇ s is used instead of the estimated rotor position ⁇ ⁇ r .
  • the rotary coordinate converter 801 uses the control coordinate angle ⁇ s to perform rotary coordinate conversion of the rotary machine currents i s ⁇ and i s ⁇ on two-phase coordinates into dq axis currents i sd and i sq on rotary coordinates.
  • the AC component estimator 803 estimates the flux linkage inductance AC component using the rotating machine current i s dq as shown in the above equation (23). Assuming that the position estimator 8 is operating such that the estimated rotor position ⁇ ⁇ r calculated by the position estimator 8 is approximately equal to the control coordinate angle ⁇ s calculated by the angle calculator 9,
  • the equation (23), which is an arithmetic equation for the flux linkage inductance alternating current, is simplified as the following equation (25).
  • [Psi ⁇ sac dq above (25) is a flux linkage inductance AC component estimation value calculated by the control unit 100B of the rotating machine according to the third embodiment.
  • the rotor position calculator 804 PI-controls the estimation error of the rotor position to calculate the estimated rotational speed ⁇ ⁇ r . Further, the rotor position calculator 804 integrates the estimated rotational speed ⁇ ⁇ r to calculate an estimated rotor position ⁇ ⁇ r .
  • FIG. 7 is a diagram showing the configuration of the control coordinate angle calculator 9 shown in FIG. Angle calculator 9, the error between the estimated rotor position theta ⁇ r and the control coordinate angle theta s, by integrating after PI control calculates a control coordinate angle theta s, estimated rotor position control coordinate angle theta s Let ⁇ ⁇ r follow.
  • the tracking response is arbitrarily adjustable by the gain of proportional integral control.
  • the control device 100B of the rotary machine according to the third embodiment does not directly use the estimated rotor position ⁇ ⁇ r as an angle used for rotational coordinate conversion of each of the position estimator 8 and the controller 10.
  • the control coordinate angle ⁇ s to follow in response to is used.
  • the response of the estimated rotor position ⁇ ⁇ r may be oscillatory.
  • the control of the rotating machine including the estimation of the rotor position can be stably operated. Therefore, the control device 100B for a rotating machine according to the third embodiment has a remarkable effect which is unprecedented in that the rotational speed and the rotational position of the rotor can be estimated with higher response and stability.
  • FIG. 8 is a diagram showing the configuration of a control device for a rotary machine according to Embodiment 4 of the present invention.
  • Control devices 100, 100A, 100B of the rotating machine according to Embodiments 1, 2, 3 estimate the rotor position using the calculated value and the estimated value of the flux linkage inductance AC component.
  • the control device 100C for a rotating machine according to the fourth embodiment is configured to estimate a rotor position using an adaptive observer.
  • An adaptive observer shows the whole observer in the case of estimating the variable (rotor position) contained in Formula (34) which is an arithmetic expression of the coefficient H of Formula (33) mentioned later by adaptive estimation.
  • the control device 100C for the rotary machine includes a position estimator 11 instead of the position estimator 6 according to the second embodiment shown in FIG.
  • the other configuration is the same as or equivalent to the configuration of the second embodiment, and the same or equivalent components are denoted by the same reference numerals and redundant description will be omitted.
  • FIG. 9 is a diagram showing the configuration of the position estimator shown in FIG.
  • the position estimator 11 includes a three-phase to two-phase converter 1101, a three-phase to two-phase converter 1102, a flux linkage calculator 1103, a rotational coordinate converter 1104, an adaptive observer 1105, and an adaptive estimator 1106.
  • the three-phase to two-phase converter 1101 coordinates-converts the rotating machine voltage command on the three-phase coordinates into the rotating machine voltage command v s ⁇ * on the two-phase coordinates.
  • Three-phase - two-phase converter 1102 coordinate transformation the rotating machine current on the three-phase coordinates to the rotating machine current i s .alpha..beta on the two-phase coordinate.
  • the flux linkage calculator 1103 computes the flux linkage ⁇ shpf ⁇ on two-phase coordinates, and further uses the estimated rotor position ⁇ ⁇ r to generate a chain on two-phase coordinates.
  • Interchange flux ⁇ shpf ⁇ is rotationally transformed into interlinkage flux s s dq on rotational coordinates.
  • rotary coordinate converter 1104 converts rotary machine voltage command v s ⁇ * on two-phase coordinates into rotary machine voltage command v s dq * on rotary coordinates. Then, the adaptive observer 1105 calculates the estimated rotor position ⁇ ⁇ r using the rotating machine voltage command v s dq * on the rotational coordinates and the linkage flux ⁇ s dq .
  • Equation (31) is an observer that estimates the flux linkage ⁇ s dq .
  • the observer is configured as the following equation (33).
  • the observer gain H by appropriately designing the observer gain H, the estimated interlinking flux ⁇ dq of the observer converges with the response ⁇ cobs that is a true value.
  • the design of the observer gain H is not limited to the example of the following equation (34).
  • the adaptive estimator 1106 calculates the rotor position from the estimated error of the flux linkage with respect to the flux linkage AC component caused by the estimated error of the rotor position.
  • FIG. 10 is a diagram showing the configuration of the adaptive estimator shown in FIG. As shown in FIG. 10, the adaptive estimator 1106 calculates the estimated rotor position ⁇ ⁇ r of the rotor using the adaptation law of the following equation (35).
  • k ap of the above equation (35) is a value calculated by the following equation (36)
  • ⁇ ai of the above equation (35) is a value calculated by the following equation (37).
  • the response is ⁇ ca.
  • k ap represents an adaptive estimated proportional gain
  • ⁇ ai represents an adaptive estimated integral response.
  • the adaptive estimator 1106 extracts an estimation error of the rotor position using the flux linkage alternating current component according to e 1 shown in the following equation (38), that is, the value included in the equation (35). There is.
  • e 1 represents “adaptive estimation error 1”.
  • Control device 100C for a rotary machine according to Embodiment 4 uses the estimation error of the flux linkage relating to the flux linkage alternating current component corresponding to the term relating to ⁇ r shown in the above equation (27) to the above equation (32) To estimate the rotor position.
  • the flux linkage AC component has components of both two phases, and by using this flux linkage AC component, it is possible to estimate the rotor position without depending on the current flow direction. Further, in control apparatus 100C for a rotating machine according to Embodiment 4, it is not necessary to perform weighting based on the direction of energization. Therefore, the control device 100C for a rotating machine according to the fourth embodiment produces a remarkable effect not hitherto found, such as being able to estimate the rotor position with a simple configuration regardless of the current flow direction of the rotating machine.
  • the controller 100C of the rotating machine according to the fourth embodiment calculates the flux linkage on the rotational coordinates synchronized with the rotor position, and further estimates the flux linkage. Therefore, the remarkable effect which has not been achieved conventionally, such as being able to estimate the rotor position accurately with a small amount of control operation, is achieved.
  • estimation of the rotor position by the control device 100C of the rotating machine according to the fourth embodiment can be carried out on two-phase coordinates based on the equations (26) to (28). It is.
  • the control device 100C for the rotary machine calculates the rotor position using the above equation (35) using the calculated value and the estimated value of the flux linkage.
  • the flux linkage inductance alternating current component estimated value ⁇ ⁇ sac dq can be calculated by the following equation (39).
  • the controller 100C for the rotary machine according to the fourth embodiment is the same as in the first, second and third embodiments, and the flux linkage inductance AC component estimation value calculated by the above equation (39) and flux linkage flux AC It goes without saying that the rotor position can be estimated using the minute operation value. Further, as in the third embodiment, control device 100C for a rotary machine according to the fourth embodiment estimates the estimated rotor position ⁇ ⁇ r as an angle used for rotational coordinate conversion of each of position estimator 11 and controller 7. It is needless to say that the control coordinate angle ⁇ s that follows this can be used, not directly.
  • FIG. 11 is a diagram showing the configuration of a control device for a rotary machine according to Embodiment 5 of the present invention.
  • Control devices 100, 100A, and 100B of the rotary machine according to Embodiments 1, 2, and 3 estimate the rotor position using the calculated value and the estimated value of the flux linkage inductance AC component, and relates to Embodiment 4.
  • the control device 100C of the rotating machine estimates the rotor position using an adaptive observer.
  • the control device 100D for the rotary machine according to the fifth embodiment is configured to estimate the rotor position using the calculated value of the flux linkage inductance AC, the energization angle, and the rotor position.
  • the control device 100D for the rotary machine includes a position estimator 12 instead of the position estimator 3 shown in FIG.
  • the other configuration is the same as or equivalent to the configuration of the first embodiment, and the same or equivalent components are denoted by the same reference numerals and redundant description will be omitted.
  • FIG. 12 is a diagram showing the configuration of the position estimator shown in FIG.
  • the position estimator 12 includes a three-phase to two-phase converter 301, a three-phase to two-phase converter 302, an AC component calculator 303, an interlinkage flux inductance AC component vector estimator 1201, a rotor position An estimation error calculator 1202 and a rotor position calculator 306 are provided.
  • the position estimator 12 differs from that of the first embodiment in the operations of the flux linkage alternating current component vector estimator 1201 and the rotor position estimation error calculator 1202.
  • the flux linkage alternating current component rotates at a phase obtained by subtracting the energization angle of the rotating machine current from the rotor position.
  • “e ⁇ sac ⁇ ” of the above equation (40) which is a vector in the same direction as the flux linkage flux AC component, is referred to as “ link flux inductance AC component vector”.
  • the linkage flux inductance AC component vector estimator 1201 uses the conduction angle command ⁇ i * and the estimated rotor position ⁇ ⁇ r to generate a linkage flux inductance AC component vector estimated value Calculate e ⁇ ⁇ sac ⁇ .
  • the rotor position estimation error calculator 1202 calculates an estimation error of the rotor position using the linkage flux inductance alternating current calculation value ⁇ sac, calc ⁇ and the linkage flux inductance alternating current vector estimation value e ⁇ ⁇ sac ⁇ Do. Assuming that the conduction angle ⁇ i and the conduction angle command ⁇ i * are equal, the outer product of the flux linkage inductance alternating current calculation value ⁇ sac, calc ⁇ ⁇ and the linkage flux inductance alternating current component vector estimated value e ⁇ ⁇ sac ⁇ ⁇ It can be expressed by the equation (42).
  • the interlinkage flux inductance alternating current calculation value ⁇ sac, calc ⁇ is obtained by calculating the equation (40) according to the equation (9) using the rotating machine voltage command and the rotating machine current.
  • the estimated error of the rotor position can be calculated by the following equation (43).
  • the estimated rotor position is calculated from this rotor position estimation error in the same manner as in the first embodiment.
  • Control device 100D for a rotary machine uses the flux linkage inductance AC component and a vector in the same direction, that is, using the above equations (40) and (41) to determine the rotor position.
  • the flux linkage alternating current component vector which is a vector in the same direction, is calculated using the energization angle command and the estimated rotor position.
  • the flux linkage AC component and flux linkage AC component vector have both two-phase components, and by using the flux linkage AC component and flux inductance AC component vector, it depends on the direction of current flow. It is possible to estimate the rotor position without.
  • control device 100D for a rotary machine according to Embodiment 5 it is not necessary to perform weighting based on the direction of energization. Therefore, the control device 100D for a rotating machine according to the fifth embodiment produces a remarkable effect not heretofore found, such as being able to estimate the rotor position with a simple configuration regardless of the current flow direction of the rotating machine.
  • control device 100, 100A, 100B, 100C, 100D of a rotating machine according to the first to fifth embodiments can be realized using a processing circuit.
  • the respective functions are a current detector 2, a voltage applicator 5, position estimators 3, 6, 8, 11 and 12, controllers 4, 7 and 10, and an angle calculator 9.
  • FIG. 13 is a diagram showing a first hardware configuration example of a control device for a rotary electric machine according to Embodiments 1 to 5 of the present invention.
  • FIG. 13 shows an example in which the above processing circuit is realized by dedicated hardware such as the dedicated processing circuit 13.
  • the dedicated processing circuit 13 may be a single circuit, a composite circuit, a programmed processor, a parallel programmed processor, an application specific integrated circuit (ASIC), an FPGA (field) -Programmable Gate Array) or a combination thereof.
  • ASIC application specific integrated circuit
  • FPGA field-Programmable Gate Array
  • FIG. 14 is a diagram showing a second hardware configuration example of the control device for a rotary electric machine according to Embodiments 1 to 5.
  • FIG. 14 shows an example in which the processor 14 and the storage unit 15 realize the above processing circuit.
  • each of the functions described above is realized by software, firmware or a combination thereof.
  • the software or firmware is described as a program and stored in the storage device 15.
  • the processor 14 reads and executes the program stored in the storage device 15. It can also be said that these programs cause a computer to execute the procedures and methods performed by each of the above functions.
  • the storage device 15 is a semiconductor memory such as a random access memory (RAM), a read only memory (ROM), a flash memory, an erasable programmable read only memory (EPROM), or an EEPROM (registered trademark) (electrically erasable programmable read only memory). Do.
  • the semiconductor memory may be non-volatile memory or volatile memory.
  • the storage device 15 corresponds to a magnetic disk, a flexible disk, an optical disk, a compact disk, a mini disk, or a DVD (Digital Versatile Disc).
  • each of the functions described above may be partially realized by hardware and partially by software or firmware.
  • the functions of the current detector 2 and the voltage applicator 5 are realized using dedicated hardware, and the position estimators 3, 6, 8, 11, 12; the controllers 4, 7, 10;
  • the function of the computing unit 9 is realized using the processor 14 and the storage unit 15.
  • the current command for the torque of the rotating machine is selected such that the effective value, ie, the copper loss of the rotating machine is minimized, but the flux linkage, ie, the induced voltage of the rotating machine is reduced to rotate
  • the machine may be set to operate at higher speeds or may be set to minimize motor loss.
  • the phase current flowing to the rotating machine is the current incorporated in the inverter that is the voltage application unit 5 It may be detected by a sensor, a shunt resistor or the like.
  • the configuration shown in the above embodiment shows an example of the contents of the present invention, and can be combined with another known technique, and one of the configurations is possible within the scope of the present invention. Parts can be omitted or changed.

Landscapes

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

Abstract

回転機の制御装置(100)は、回転機(1)に流れる回転機電流を検出する電流検出器(2)と、回転子位置を推定する位置推定器(3)と、回転機電流と回転子位置とに基づいて回転機(1)を駆動するための回転機電圧指令を出力する制御器(4)と、回転機電圧指令に基づいて回転機(1)に電圧を印加する電圧印加器(5)とを備える。位置推定器(3)は、インダクタンス交流分と回転機電流とによって生成される鎖交磁束インダクタンス交流分から回転子位置を推定する。

Description

回転機の制御装置
 本発明は、インダクタンスが回転子位置によって変化するインダクタンス交流分を有する回転機を、回転子位置を検出する位置センサを用いることなく回転子位置情報を得て制御する、回転機の制御装置に関する。
 回転機の性能を十分に引き出して回転機を駆動するためには、回転子の位置情報が必要である。そのため、従来の回転機の制御装置は、回転機に取付けられた位置センサで検出される位置情報を用いている。しかしながら、回転機の製造コストのより一層の低減、回転機の小型化、及び回転機の信頼性の向上といった観点から、位置センサレスで回転機を駆動する技術が開発されてきた。回転機の位置センサレス制御方法には、高周波電圧を回転機に印加することによって回転子位置を推定する方法と、高周波電圧を印加せずに回転機の誘起電圧、鎖交磁束などから回転子位置を推定する方法とがある。
 特許文献1には、高周波電圧を回転機に印加することによって回転子位置を推定する方法が開示されている。特許文献1に開示される位置センサレス制御方法では、高周波電圧を回転機に印加したときの回転機電流が検出され、高周波電圧と同じ周波数成分の高周波電流が抽出される。そして、回転機のインダクタンス、すなわち、高周波電流の振幅が回転子位置電気角に対して2倍の周波数で変化することを利用して、回転子位置が推定される。このような高周波電圧を利用する方式では、回転機が零速又は零速に近い低速でも、正確に回転子位置を推定できる一方で、重畳される高周波電圧によってトルク脈動又は騒音が発生する。また、高周波電圧を利用する方式では、回転機の巻線に印加される電圧と、回転機の巻線に流れる電流とに高周波電圧が重畳される分だけ、回転機の最大トルク又は回転機の回転速度が小さくなり、もしくはその両方が小さくなる。
 特許文献2,3,4には、高周波電圧を印加せずに回転子位置を推定する方法が開示されている。高周波電圧を印加せずに回転子位置を推定する方法では、回転機の鎖交磁束から回転子のq軸又はd軸のインダクタンスと回転機電流との積を減算することで、鎖交磁束の内、回転子位置に同期して回転する成分を抽出する。当該成分には、d軸基準のActiveFluxと、q軸基準のActiveFluxとがそれぞれ含まれる。ここで、インダクタンスが最大になる回転子の方向をd軸とし、インダクタンスが最小になる方向をq軸とする。d軸は磁束軸と称される軸であり、q軸はトルク軸と称される軸である。特許文献2では、ActiveFluxをオブザーバで推定することによって回転子位置が推定される。特許文献3では、ActiveFluxによって発生する誘起電圧を用いて回転子位置が推定される。
 特許文献2,3に開示される位置センサレス制御方法では、回転機電流のd軸成分又はq軸成分が零の場合、ActiveFluxが零となり、回転子位置を推定できない。すなわち、特許文献2,3に開示される位置センサレス制御方法では、回転機の通電方向によっては、回転子位置を推定できない。特許文献4に開示される位置センサレス制御方法では、d軸基準とq軸基準の両方のActiveFluxに、回転機電流又は鎖交磁束のd軸成分とq軸成分に基づいた重み付けを行い、回転子位置の推定に利用することによって、この問題を解決している。
特許第5069306号公報 特許第4644010号公報 特許第3571698号公報 欧州特許出願公開第2493067号明細書
 しかしながら、特許文献4に開示される、回転機電流又は鎖交磁束のd軸成分とq軸成分に基づいた重み付けは、回転子位置推定の制御設計だけでなく制御処理を複雑にするという課題がある。そのため、特許文献4に開示される技術では、回転子位置を推定するための構成が複雑になるという課題がある。
 本発明は、上記に鑑みてなされたものであって、回転機の通電方向に依らず簡単な構成で回転子位置を推定することができる回転機の制御装置を得ることを目的とする。
 上述した課題を解決し、目的を達成するために、本発明の回転機の制御装置は、インダクタンスが回転子位置によって変化するインダクタンス交流分を有する回転機を制御する回転機の制御装置であって、回転機に流れる回転機電流を検出する電流検出器を備える。回転機の制御装置は、回転子位置を推定する位置推定器と、回転機電流と回転子位置とに基づいて回転機を駆動するための回転機電圧指令を出力する制御器と、回転機電圧指令に基づいて回転機に電圧を印加する電圧印加器とを備える。位置推定器は、インダクタンス交流分と回転機電流とによって生成される鎖交磁束インダクタンス交流分から回転子位置を推定することを特徴とする。
 本発明によれば、回転機の通電方向に依らず簡単な構成で回転子位置を推定することができる、という効果を奏する。
本発明の実施の形態1に係る回転機の制御装置の構成を示す図 図1に示す位置推定器の構成を示す図 本発明の実施の形態2に係る回転機の制御装置の構成を示す図 図3に示す位置推定器の構成を示す図 本発明の実施の形態3に係る回転機の制御装置の構成を示す図 図5に示す位置推定器の構成を示す図 図5に示す制御座標角度演算器の構成を示す図 本発明の実施の形態4に係る回転機の制御装置の構成を示す図 図8に示す位置推定器の構成を示す図 図9に示す適応推定器の構成を示す図 本発明の実施の形態5に係る回転機の制御装置の構成を示す図 図11に示す位置推定器の構成を示す図 本発明の実施の形態1から5に係る回転電機の制御装置の第1のハードウェア構成例を示す図 実施の形態1から5に係る回転電機の制御装置の第2のハードウェア構成例を示す図
 以下に、本発明の実施の形態に係る回転機の制御装置を図面に基づいて詳細に説明する。なお、この実施の形態によりこの発明が限定されるものではない。
実施の形態1.
 図1は本発明の実施の形態1に係る回転機の制御装置の構成を示す図である。実施の形態1に係る回転機の制御装置100は、回転機を駆動するための回転機電圧指令vsu ,vsv ,vsw に従って、交流電圧を回転機1に印加する電圧印加器5と、電圧印加器5から回転機1に供給される交流電流を検出し、検出した交流電流を回転機電流isu,isv,iswとして出力する電流検出器2とを備える。また回転機の制御装置100は、回転機電圧指令vsu ,vsv ,vsw と回転機電流isu,isv,iswとを用いて、電気角である推定回転子位置θ を演算する位置推定器3と、制御器4とを備える。推定回転子位置θ は、回転機1を構成する図示しない回転子の回転位置である回転子位置の推定値である。推定回転子位置θ は、電気角で表わされる。回転機の制御装置100によって制御される回転機1は同期リラクタンスモータである。同期リラクタンスモータは、インダクタンスが回転子位置によって変化するインダクタンス交流分を有するモータである。インダクタンスとは、固定子巻線の鎖交磁束を固定子電流値で除した値である。以下では回転機1を単に「回転機」と称する場合がある。電圧印加器5は、U,V,Wの3本の配線を介して、回転機に接続される。電流検出器2では、当該配線に流れる交流電流が検出される。以下では、インダクタンスが最大になる回転子の方向をd軸とし、インダクタンスが最小になる方向をq軸とし、回転子位置は回転子のd軸を基準とする。
 制御器4は、電流指令演算器401、三相電流指令演算器402及び三相電流制御器403を備え、回転機がトルク指令値Tに対応した出力を発生するための回転機電圧指令vsu ,vsv ,vsw を演算する。
 電流指令演算器401は、回転機がトルク指令値Tに対応した出力を発生するために必要な電流実効値指令Iph 及び通電角度指令φ を演算する。通電角度指令φ は、実際の回転子位置と回転機電流との成す角度を示す通電角の指令である。ここでは、電流実効値指令Iph と通電角度指令φ とは、トルクに対する電流実効値、すなわち回転機の銅損が最小になるように選定される。
 三相電流指令演算器402は、下記(1)式に示すように、電流実効値指令Iph と通電角度指令φ と推定回転子位置θ とを用いて、三相電流指令isu ,isv ,isw を演算する。
Figure JPOXMLDOC01-appb-M000002
 三相電流制御器403は、回転機電圧指令vsu ,vsv ,vsw を演算する。回転機電圧指令vsu ,vsv ,vsw は、三相座標上の回転機電流isu,isv,iswが三相電流指令isu ,isv ,isw となるように、回転機電流isu,isv,iswを例えば比例積分(Proportional Integral:PI)制御することで得られる。
 図2は図1に示す位置推定器の構成を示す図である。図2に示すように位置推定器3は、三相-二相変換器301、三相-二相変換器302、鎖交磁束インダクタンス交流分演算器303、鎖交磁束インダクタンス交流分推定器304、回転子位置推定誤差演算器305及び回転子位置演算器306を備える。以下では、鎖交磁束インダクタンス交流分演算器303を「交流分演算器303」と称し、鎖交磁束インダクタンス交流分推定器304を「交流分推定器304」と称し、回転子位置推定誤差演算器305を「推定誤差演算器305」と称する場合がある。
 三相-二相変換器301は、三相座標上の回転機電圧指令vsu ,vsv ,vsw を二相座標上の回転機電圧指令vsα ,vsβ に変換する。なお、以下では「vsα ,vsβ 」を「v αβ*」と表記する場合がある。実施の形態1では、三相-二相変換に下記(2)式の変換行列C32が利用される。
Figure JPOXMLDOC01-appb-M000003
 同様に、三相-二相変換器302は、三相座標上の回転機電流isu,isv,iswを二相座標上の回転機電流isα,isβに変換する。なお、以下では「isα,isβ」を「i αβ*」と表記する場合がある。
 次に、交流分演算器303の動作について説明する。まず、回転機の鎖交磁束Ψ αβは、下記(3)式により求められる。
Figure JPOXMLDOC01-appb-M000004
 上記(3)式のv αβ*は二相座標での回転機電圧指令であり、回転機電流i αβは二相座標上の回転機電流であり、Rは巻線抵抗である。また上記(3)式の積分は、ラプラス変換におけるs領域において下記(4)式で表される。
Figure JPOXMLDOC01-appb-M000005
 回転機の鎖交磁束Ψ αβを積分で演算する場合、通常は初期値が不明であるため、回転機の鎖交磁束Ψ αβの基本波周波数成分に対して十分にカットオフ周波数の低いハイパスフィルタ(High Pass Filter:HPF)を利用する。そのHPFの伝達関数は、カットオフ周波数をωhpfとして下記(5)式で表される。
Figure JPOXMLDOC01-appb-M000006
 このHPFを上記(4)式に適用すると、鎖交磁束Ψshpf αβは下記(6)式で計算される。
Figure JPOXMLDOC01-appb-M000007
 上記(6)式を変形すると下記(7)式となる。
Figure JPOXMLDOC01-appb-M000008
 一方で、回転機の鎖交磁束Ψ αβは下記(8)式で表される。
Figure JPOXMLDOC01-appb-M000009
 上記(8)式のLsdcは、回転子位置によって変化しないインダクタンス直流分であり、Lmacは回転子位置によって変化するインダクタンス交流分である。Lsac αβはインダクタンス交流分を表す。また、インダクタンスの変化は一般的に、電気角で表される回転子位置をθとすれば、2θの正弦関数又は余弦関数で表される。実施の形態1では、インダクタンス交流分と回転機電流i αβとによって生成される鎖交磁束を、鎖交磁束インダクタンス交流分とする。交流分演算器303は、上記(7)式で計算した鎖交磁束Ψshpf αβを用いて、下記(9)式により、鎖交磁束インダクタンス交流分Ψsac,calc αβを演算する。
Figure JPOXMLDOC01-appb-M000010
 上記(9)式で演算される鎖交磁束インダクタンス交流分Ψsac,calc αβは、実施の形態1に係る回転機の制御装置100で演算される鎖交磁束インダクタンス交流分演算値である。
 交流分推定器304は、下記(10)式に示すように、回転機電流i αβと推定回転子位置θ とを用いて、鎖交磁束インダクタンス交流分を推定する。
Figure JPOXMLDOC01-appb-M000011
 上記(10)式で演算されるΨ sac αβは、実施の形態1に係る回転機の制御装置100で演算される鎖交磁束インダクタンス交流分推定値である。上記(10)式の「L sac αβ」は、「インダクタンス交流分推定値」である。
 推定誤差演算器305は、鎖交磁束インダクタンス交流分演算値と鎖交磁束インダクタンス交流分推定値とを用いて、回転子位置の推定誤差を演算する。鎖交磁束インダクタンス交流分演算値と鎖交磁束インダクタンス交流分推定値との外積は、上記(9)式すなわち上記(8)式の第2項を演算した値と、上記(10)式とを用いて、下記(11)式で表される。
Figure JPOXMLDOC01-appb-M000012
 ここで、回転子位置の推定値と真値がおよそ等しい、すなわちθ ≒θとすると、回転子位置の推定誤差は下記(12)式で演算できる。推定誤差演算器305は、演算した回転子位置の推定誤差を回転子位置演算器306へ出力する。
Figure JPOXMLDOC01-appb-M000013
 回転子位置演算器306は、回転子位置の推定誤差をPI制御した後に積分して、推定回転子位置θ を演算する。
 前述した特許文献2に開示される技術では、下記(13)式に示すように、回転機の鎖交磁束Ψ αβから、回転子のq軸インダクタンスLsqと回転機電流i αβとの積を減算することで、鎖交磁束のうち回転子位置に同期して回転する成分であるd軸基準のActiveFlux(ΨAFd αβ)が抽出される。
Figure JPOXMLDOC01-appb-M000014
 上記(13)式に示す回転子のq軸インダクタンスLsqは、Lsq=Lsdc-Lmacで表すことができる。また、回転機電流i αβは下記(14)式で表される。
Figure JPOXMLDOC01-appb-M000015
 上記(14)式のIphは電流実効値であり、φは通電角度、すなわち回転子位置と回転機電流との角度差である。
 上記(8)式と上記(14)式とを上記(13)式の右辺へ代入すると、d軸基準のActiveFlux(ΨAFd αβ)は下記(15)式で表される。
Figure JPOXMLDOC01-appb-M000016
 上記(15)式のLsdは、Lsd=Lsdc+Lmacで表すことができる。また、d-q軸電流isd,isqは、下記(16)式に示すように、α-β軸電流を回転子位置θを用いて回転座標変換して得られる。なお、以下ではd-q軸電流isd,isqを「回転機電流i dq」と称する場合がある。
Figure JPOXMLDOC01-appb-M000017
 実施の形態1の回転座標変換には、上記(16)式の変換行列Cdqが利用される。上記(15)式のActiveFluxの方向は回転子d軸方向であるので、これをオブザーバ等を用いて推定することで、回転子位置を演算できる。前述した特許文献3に開示される技術では、d軸基準のActiveFluxによって発生する誘起電圧を用いて回転子位置の推定が行われている。また前述した特許文献2には、下記(17)式に示すように、回転機の鎖交磁束から回転子のd軸インダクタンスと回転機電流との積を減算した、q軸基準のActiveFluxについても記載されている。
Figure JPOXMLDOC01-appb-M000018
 d軸基準の場合と同様に、上記(17)式の右辺に上記(8)式と上記(14)式を代入すると、q軸基準のActiveFlux(ΨAFq αβ)は下記(18)式で表される。
Figure JPOXMLDOC01-appb-M000019
 q軸基準のActiveFluxは回転子q軸方向であるため、これをオブザーバ等を用いて推定することで、回転子位置を推定できる。しかしながら、d軸基準のActiveFluxは、通電方向が回転子q軸方向のみ場合、すなわちisd=0の場合には零となり、これを用いた回転子位置の推定ができない。同様に、q軸基準のActiveFluxは、通電方向が回転子d軸方向のみ場合、すなわちisq=0の場合には零となり、これを用いた回転子位置の推定ができない。このような問題に対して、前述した特許文献4に開示される位置センサレス制御方法は、d軸基準とq軸基準の両方のActiveFluxに、回転機電流又は鎖交磁束のd軸成分とq軸成分に基づいた重み付けを行い、回転子位置の推定に利用することによって、この問題を解決している。しかしながら、回転機電流又は鎖交磁束のd軸成分とq軸成分に基づいた重み付けは、回転子位置推定の制御設計だけでなく制御処理を複雑にする。
 実施の形態1に係る回転機の制御装置100は、鎖交磁束インダクタンス交流分、すなわち上記(8)式の第2項と上記(10)式とを用いて回転子位置を推定する。この鎖交磁束インダクタンス交流分は二相両方の成分を持っており、この鎖交磁束インダクタンス交流分を用いることによって、電流の通電方向に依存せずに回転子位置の推定が可能である。また、実施の形態1に係る回転機の制御装置100では、通電方向によって重み付けを行う必要もない。従って実施の形態1に係る回転機の制御装置100は、回転機の通電方向に依らず簡単な構成で回転子位置を推定できる、という従来にない顕著な効果を奏する。
実施の形態2.
 図3は本発明の実施の形態2に係る回転機の制御装置の構成を示す図である。実施の形態1に係る回転機の制御装置100では、回転機のインダクタンスの変化が、回転子位置θとして2θの正弦関数又は余弦関数で表される。しかしながら実際のインダクタンスの変化は、正確には2θの正弦関数又は余弦関数とはならず、推定回転子位置に誤差が発生する。また実施の形態1に係る回転機の制御装置100は、上記(10)式の鎖交磁束インダクタンス交流分を推定するために2θの正弦関数と余弦関数が必要であり、制御演算量が多くなる。実施の形態2に係る回転機の制御装置100Aは、少ない制御演算量で正確に回転子位置を推定できるように構成されている。実施の形態2に係る回転機の制御装置100Aは、図1に示す位置推定器3及び制御器4の代わりに、位置推定器6及び制御器7を備える。その他の構成については、実施の形態1の構成と同一又は同等であり、同一又は同等の構成部には同一の符号を付して、重複する説明は省略する。
 制御器7は、電流指令演算器701、d-q電流制御器702、回転座標逆変換器703、二相-三相変換器704、三相-二相変換器705及び回転座標変換器706を備え、回転機がトルク指令値Tに対応した出力を発生するための回転機電圧指令vsu ,vsv ,vsw を演算する。
 電流指令演算器701は、回転機がトルク指令値Tに対応した出力を発生するために必要な回転座標上の電流指令isd ,isq を演算する。ここでは、回転座標上の電流指令isd ,isq は、トルクに対する電流実効値、すなわち回転機の銅損が最小になるように選定される。
 d-q電流制御器702は、回転座標変換器706で回転座標変換されるd-q軸電流isd,isqが電流指令isd ,isq となるように制御を行い、回転座標上の回転機電圧指令vsd ,vsq を演算する。この電流制御には、例えば比例積分制御などを利用する。
 回転座標逆変換器703は、下記(19)式に示すように、位置推定器6で演算された推定回転子位置θ を用いて、回転座標上の回転機電圧指令vsd ,vsq を二相座標上の回転機電圧指令vsα ,vsβ へ回転座標逆変換する。実施の形態2では、回転座標逆変換に下記(19)式の変換行列Cdq -1(θ )が利用される。
Figure JPOXMLDOC01-appb-M000020
 二相-三相変換器704は、下記(20)式に示すように、二相座標上の回転機電圧指令vsα ,vsβ を三相座標上の回転機電圧指令vsu ,vsv ,vsw に変換する。実施の形態2では、二相-三相変換に下記(20)式の変換行列C23が利用される。
Figure JPOXMLDOC01-appb-M000021
 三相-二相変換器705は、三相座標上の回転機電流isu,isv,iswを二相座標上の回転機電流isα,isβへ三相-二相変換する。回転座標変換器706は、推定回転子位置θ を用いて、二相座標上の回転機電流isα,isβを回転座標上のd-q軸電流isd,isqへ回転座標変換する。
 図4は図3に示す位置推定器の構成を示す図である。図4に示すように位置推定器6は、三相-二相変換器601、三相-二相変換器602、回転座標変換器603、鎖交磁束インダクタンス交流分演算器604、鎖交磁束インダクタンス交流分推定器605、回転子位置推定誤差演算器606及び回転子位置演算器607を備える。位置推定器6では、三相-二相変換器602、鎖交磁束インダクタンス交流分演算器604及び鎖交磁束インダクタンス交流分推定器605の動作が実施の形態1と異なる。以下では、鎖交磁束インダクタンス交流分演算器604を「交流分演算器604」と称し、鎖交磁束インダクタンス交流分推定器605を「交流分推定器605」と称し、回転子位置推定誤差演算器606を「推定誤差演算器606」と称する場合がある。
 三相-二相変換器601は、実施の形態1の三相-二相変換器301と同様に、三相座標上の回転機電圧指令vsu ,vsv ,vsw を二相座標上の回転機電圧指令vsα ,vsβ に変換する。三相-二相変換器602は、実施の形態1の三相-二相変換器302と同様に、三相座標上の回転機電流isu,isv,iswを二相座標上の回転機電流isα,isβに変換する。
 回転機の鎖交磁束Ψ αβは二相座標上において上記(8)式で表されるが、下記(21)式で示される回転機の鎖交磁束Ψ dqは、制御座標角度θを用いて、回転機の鎖交磁束Ψ αβを回転座標変換したものである。下記(21)式の第2項は、インダクタンス交流分と回転機電流とによって生成される鎖交磁束インダクタンス交流分である。
Figure JPOXMLDOC01-appb-M000022
 交流分演算器604は、まず実施の形態1と同様に上記(7)式を用いて、二相座標での鎖交磁束Ψshpf αβを演算する。さらに、交流分演算器604は、鎖交磁束Ψshpf αβを回転座標変換して、鎖交磁束Ψshpf dqとする。回転座標上での鎖交磁束インダクタンス交流分Ψsac,calc dqは、上記(21)式に従って、下記(22)式で演算される。下記(22)式で演算される鎖交磁束インダクタンス交流分Ψsac,calc dqは、以下では鎖交磁束インダクタンス交流分演算値と称する。
Figure JPOXMLDOC01-appb-M000023
 回転座標変換器603は、推定回転子位置θ を用いて、二相座標上の回転機電流isα,isβを回転座標上のd-q軸電流isd,isqへ回転座標変換する。
 交流分推定器605は、下記(23)式に示すように、回転機電流i dqを用いて、鎖交磁束インダクタンス交流分を推定する。
Figure JPOXMLDOC01-appb-M000024
 ここで、上記(23)式では、制御座標角度θとして推定回転子位置θ を用いているので、上記(23)式は下記(24)式に示すように簡略化される。
Figure JPOXMLDOC01-appb-M000025
 上記(24)式のΨ sac dqは、実施の形態2に係る回転機の制御装置100Aで演算される鎖交磁束インダクタンス交流分推定値である。
 推定誤差演算器606は、実施の形態1の推定誤差演算器305と同様に、鎖交磁束インダクタンス交流分演算値と鎖交磁束インダクタンス交流分推定値とを用いて、回転子位置の推定誤差を演算する。回転子位置演算器607は、実施の形態1の回転子位置演算器306と同様に、回転子位置の推定誤差をPI制御した後に積分して、推定回転子位置θ を演算する。
 実施の形態2に係る回転機の制御装置100Aは、回転子位置に同期した回転座標上で上記(22)式に示すように鎖交磁束インダクタンス交流分を演算し、さらに上記(24)式に示すように鎖交磁束インダクタンス交流分を推定している。この鎖交磁束インダクタンス交流分は二相両方の成分を持っており、この鎖交磁束インダクタンス交流分を用いることによって、電流の通電方向に依存せずに回転子位置の推定が可能である。また、実施の形態2に係る回転機の制御装置100Aでは、通通電方向によって重み付けを行う必要もない。さらに、実施の形態2では、回転座標上での鎖交磁束インダクタンス交流分の演算式である上記(24)式において、回転座標上のd軸におけるインダクタンス値と、回転座標上のq軸におけるインダクタンス値とを使用しているため、回転子位置をθとしたときに、インダクタンスを2θの正弦関数又は余弦関数で表す必要がない。従って、実際の回転機のインダクタンスは正確には2θの正弦関数又は余弦関数とならないことによる回転子推定誤差が発生しない。また、実施の形態2では、鎖交磁束インダクタンス交流分の推定に、2θの正弦関数と余弦関数を必要としない。従って、実施の形態2に係る回転機の制御装置100Aは、少ない制御演算量で正確に回転子位置を推定できる、といった従来にない顕著な効果を奏する。
実施の形態3. 
 図5は本発明の実施の形態3に係る回転機の制御装置の構成を示す図である。実施の形態3に係る回転機の制御装置100Bは、より高応答かつ安定に回転子の回転速度と回転位置を推定できるように構成されている。回転機の制御装置100Bは、図3に示す実施の形態2に係る位置推定器6及び制御器7の代わりに位置推定器8及び制御器10を備える。さらに回転機の制御装置100Bは、制御座標角度演算器9を備える。以下では、制御座標角度演算器9を「角度演算器9」と称する場合がある。その他の構成については、実施の形態2の構成と同一又は同等であり、同一又は同等の構成部には同一の符号を付して、重複する説明は省略する。
 制御器10は、速度演算器1001、電流指令演算器701、d-q電流制御器702、回転座標逆変換器1002、二相-三相変換器704、三相-二相変換器705及び回転座標変換器1003を備え、回転機が回転速度指令値ω に対応した回転速度を発生するための回転機電圧指令vsu ,vsv ,vsw を演算する。
 速度演算器1001は、推定回転速度ω が回転速度指令値ω になるように、例えば比例積分制御を行って、トルク指令値Tを演算する。
 回転座標逆変換器1002は、角度演算器9で演算された制御座標角度θを用いて、回転座標上の回転機電圧指令vsd ,vsq を二相座標上の回転機電圧指令vsα ,vsβ へ回転座標逆変換する。
 回転座標変換器1003は、角度演算器9で演算された制御座標角度θを用いて、二相座標上の回転機電流isα,isβを回転座標上のd-q軸電流isd,isqへ回転座標変換する。
 位置推定器8は、推定回転速度ω と推定回転子位置θ を演算する。図6は図5に示す位置推定器の構成を示す図である。図6に示すように位置推定器8は、三相-二相変換器601、三相-二相変換器602、回転座標変換器801、鎖交磁束インダクタンス交流分演算器802、鎖交磁束インダクタンス交流分推定器803、推定誤差演算器606及び回転子位置演算器804を備える。以下では鎖交磁束インダクタンス交流分演算器802を「交流分演算器802」と称し、鎖交磁束インダクタンス交流分推定器803を「交流分推定器803」と称する場合がある。
 交流分演算器802は、実施の形態2と同様に、回転座標上での鎖交磁束インダクタンス交流分Ψsac,calc dq、すなわち鎖交磁束インダクタンス交流分演算値を計算する。ただし回転座標変換には、推定回転子位置θ の代わりに制御座標角度θが用いられる。
 回転座標変換器801は、制御座標角度θを用いて、二相座標上の回転機電流isα,isβを回転座標上のd-q軸電流isd,isqへ回転座標変換する。
 交流分推定器803は、上記(23)式に示すように、回転機電流i dqを用いて、鎖交磁束インダクタンス交流分を推定する。位置推定器8で演算される推定回転子位置θ が角度演算器9で演算される制御座標角度θとおよそ等しくなるように、位置推定器8が動作していると仮定した場合、鎖交磁束インダクタンス交流分の演算式である上記(23)式は、下記(25)式のように簡略化される。
Figure JPOXMLDOC01-appb-M000026
 上記(25)式のΨ sac dqは、実施の形態3に係る回転機の制御装置100Bで演算される鎖交磁束インダクタンス交流分推定値である。
 回転子位置演算器804は、回転子位置の推定誤差をPI制御して推定回転速度ω を演算する。さらに回転子位置演算器804は、推定回転速度ω を積分して、推定回転子位置θ を演算する。
 図7は図5に示す制御座標角度演算器9の構成を示す図である。角度演算器9は、推定回転子位置θ と制御座標角度θとの誤差を、PI制御した後に積分して制御座標角度θを演算し、制御座標角度θを推定回転子位置θ に追従させる。その追従応答は、比例積分制御のゲインによって任意に調整可能である。
 実施の形態3に係る回転機の制御装置100Bは、位置推定器8及び制御器10のそれぞれの回転座標変換に用いる角度として、推定回転子位置θ を直接用いるのではなく、これに任意の応答で追従する制御座標角度θを利用している。回転子位置及び回転速度をより高応答に推定する場合には、推定回転子位置θ の応答が振動的になる場合がある。この場合も、制御座標角度θの追従をその振動を除去できるように位置推定応答より遅い応答に設定すれば、回転子位置の推定を含む回転機の制御を安定に動作させることができる。従って、実施の形態3に係る回転機の制御装置100Bは、より高応答かつ安定に回転子の回転速度と回転位置を推定できる、といった従来にない顕著な効果を奏する。
実施の形態4.
 図8は本発明の実施の形態4に係る回転機の制御装置の構成を示す図である。実施の形態1,2,3に係る回転機の制御装置100,100A,100Bは、鎖交磁束インダクタンス交流分の演算値と推定値とを用いて回転子位置を推定する。これに対して実施の形態4に係る回転機の制御装置100Cでは、適応オブザーバを用いて回転子位置を推定するように構成されている。適応オブザーバは、後述する(33)式の係数Hの演算式である(34)式に含まれる変数(回転子位置)を適応推定によって推定する場合のオブザーバ全体を示す。回転機の制御装置100Cは、図3に示す実施の形態2に係る位置推定器6の代わりに位置推定器11を備える。その他の構成については、実施の形態2の構成と同一又は同等であり、同一又は同等の構成部には同一の符号を付して、重複する説明は省略する。
 図9は図8に示す位置推定器の構成を示す図である。位置推定器11は、三相―二相変換器1101、三相―二相変換器1102、鎖交磁束インダクタンス演算器1103、回転座標変換器1104、適応オブザーバ1105及び適応推定器1106を備える。
 三相―二相変換器1101は、三相座標上の回転機電圧指令を二相座標上の回転機電圧指令v αβ*へ座標変換する。三相―二相変換器1102は、三相座標上の回転機電流を二相座標上の回転機電流i αβへ座標変換する。鎖交磁束インダクタンス演算器1103は、実施の形態2と同様に、二相座標上の鎖交磁束Ψshpf αβを演算し、さらに推定回転子位置θ を用いて、二相座標上の鎖交磁束Ψshpf αβを回転座標上の鎖交磁束Ψ dqへ回転座標変換する。
 一方、回転座標変換器1104は、二相座標上の回転機電圧指令v αβ*を回転座標上での回転機電圧指令v dq*に変換する。そして、適応オブザーバ1105は、回転座標上の回転機電圧指令v dq*と鎖交磁束Ψ dqとを用いて、推定回転子位置θ を演算する。
 適応オブザーバ1105の動作について説明する。まず、回転機のモデルは二相座標上において下記(26)式及び(27)式で表される。
Figure JPOXMLDOC01-appb-M000027
Figure JPOXMLDOC01-appb-M000028
 鎖交磁束を変数とすれば、上記(26)式及び(27)式のモデルは下記(28)式で表される。
Figure JPOXMLDOC01-appb-M000029
 上記(26)式及び(27)式のモデルを、回転座標上では下記(29)式及び(30)式で表される。
Figure JPOXMLDOC01-appb-M000030
Figure JPOXMLDOC01-appb-M000031
 鎖交磁束を変数とすれば、上記(29)式及び(30)式のモデルは下記(31)式で表される。下記(31)式は、鎖交磁束Ψ dqを推定するオブザーバである。
Figure JPOXMLDOC01-appb-M000032
 さらに、上記(31)式のモデルは、回転子位置θと制御座標角度θがおよそ等しい、すなわちθ≒θとすれば下記(32)式で表される。
Figure JPOXMLDOC01-appb-M000033
 回転子位置を推定するため、まずオブザーバを下記(33)式のように構成する。
Figure JPOXMLDOC01-appb-M000034
 ここで、下記(34)式に示すように、オブザーバゲインHを適切に設計することで、オブザーバの推定鎖交磁束Ψ^dqは真値である応答ωcobsで収束する。なお、オブザーバゲインHの設計が下記(34)式の例に限定されるものではない。
Figure JPOXMLDOC01-appb-M000035
 オブザーバにおいて、推定回転子位置θ に誤差がある場合、推定鎖交磁束Ψ^dqに誤差が表れる。適応推定器1106は、回転子位置の推定誤差によって生じる鎖交磁束インダクタンス交流分に関する鎖交磁束の推定誤差から、回転子位置を演算する。
 図10は図9に示す適応推定器の構成を示す図である。図10に示すように適応推定器1106は、下記(35)式の適応則を用いて回転子の推定回転子位置θ を演算する。
Figure JPOXMLDOC01-appb-M000036
 ここで、上記(35)式のkapを下記(36)式で演算される値とし、上記(35)式のωaiを下記(37)式で演算される値とし、回転子位置の推定応答はωcaとなる。上記(35)式のkapは適応推定比例ゲインを表し、上記(35)式のωaiは適応推定積分応答を表す。
Figure JPOXMLDOC01-appb-M000037
Figure JPOXMLDOC01-appb-M000038
 また、適応推定器1106は、下記(38)式に示すe、すなわち上記(35)式に含まれる値によって、鎖交磁束インダクタンス交流分を利用して回転子位置の推定誤差を抽出している。下記(38)式のeは「適応推定誤差1」を表す。
Figure JPOXMLDOC01-appb-M000039
 なお、回転子位置を推定する適応則の選び方は上記の方法に限定されず、例えば、以下の参考文献1に従って決めることができる。
 (参考文献1)ヨアン・D・ランダウ、富塚誠義「適応制御システムの理論と実際」オーム社、1981
 実施の形態4に係る回転機の制御装置100Cは、上記(27)式から上記(32)式までに示すθに関する項に対応する鎖交磁束インダクタンス交流分に関する鎖交磁束の推定誤差を用いて、回転子位置を推定する。この鎖交磁束インダクタンス交流分は二相両方の成分を持っており、この鎖交磁束インダクタンス交流分を用いることによって、電流の通電方向に依存せずに回転子位置の推定が可能である。また、実施の形態4に係る回転機の制御装置100Cでは、通電方向によって重み付けを行う必要もない。従って実施の形態4に係る回転機の制御装置100Cは、回転機の通電方向に依らず簡単な構成で回転子位置を推定できる、といった従来にない顕著な効果を奏する。
 実施の形態4に係る回転機の制御装置100Cは、実施の形態2と同様に、回転子位置に同期した回転座標上で鎖交磁束を演算し、さらに鎖交磁束を推定している。従って、少ない制御演算量で正確に回転子位置を推定できる、といった従来にない顕著な効果を奏する。当然ながら、実施の形態4に係る回転機の制御装置100Cによる回転子位置を推定は、上記(26)式から上記(28)式に基づいて、二相座標上で実施可能なことは自明なことである。
 実施の形態4に係る回転機の制御装置100Cでは、鎖交磁束の演算値と推定値とを用いて、上記(35)式を用いて回転子位置を演算している。ここで、鎖交磁束インダクタンス交流分推定値Ψ sac dqは、下記(39)式により演算できる。
Figure JPOXMLDOC01-appb-M000040
 実施の形態4に係る回転機の制御装置100Cは、実施の形態1,2,3と同様にして、上記(39)式により演算される鎖交磁束インダクタンス交流分推定値と鎖交磁束インダクタンス交流分演算値とを用いて、回転子位置を推定できることは言うまでもない。また実施の形態4に係る回転機の制御装置100Cは、実施の形態3と同様に、位置推定器11及び制御器7のそれぞれの回転座標変換に用いる角度として、推定回転子位置θ を直接用いるのではなく、これに追従する制御座標角度θを利用できることは言うまでもない。
実施の形態5.
 図11は本発明の実施の形態5に係る回転機の制御装置の構成を示す図である。実施の形態1,2,3に係る回転機の制御装置100,100A,100Bは、鎖交磁束インダクタンス交流分の演算値と推定値を用いて回転子位置を推定し、実施の形態4に係る回転機の制御装置100Cは、適応オブザーバを用いて回転子位置を推定する。実施の形態5に係る回転機の制御装置100Dは、鎖交磁束インダクタンス交流分の演算値と通電角度と回転子位置とを用いて、回転子位置の推定を行うように構成されている。回転機の制御装置100Dは、図1に示す位置推定器3の代わりに、位置推定器12を備える。その他の構成については、実施の形態1の構成と同一又は同等であり、同一又は同等の構成部には同一の符号を付して、重複する説明は省略する。
 図12は図11に示す位置推定器の構成を示す図である。図12に示すように位置推定器12は、三相-二相変換器301、三相-二相変換器302、交流分演算器303、鎖交磁束インダクタンス交流分ベクトル推定器1201、回転子位置推定誤差演算器1202及び回転子位置演算器306を備える。位置推定器12では、鎖交磁束インダクタンス交流分ベクトル推定器1201及び回転子位置推定誤差演算器1202の動作が実施の形態1と異なる。
 まず、上記(8)式の第2項に示す鎖交磁束インダクタンス交流分を、上記(14)式の電流実効値指令Iphを用いて演算すると、下記(40)式に示すΨsac αβを演算できる。下記(40)式のΨsac αβは、二相座標上での鎖交磁束インダクタンス交流分を表す。
Figure JPOXMLDOC01-appb-M000041
 上記(40)式より、鎖交磁束インダクタンス交流分は、回転子位置から回転機電流の通電角度を減じた位相で回転している。また、実施の形態5では、鎖交磁束インダクタンス交流分と同方向のベクトルである上記(40)式の「eΨsac αβ」を「鎖交磁束インダクタンス交流分ベクトル」と称する。
 鎖交磁束インダクタンス交流分ベクトル推定器1201は、下記(41)式に示すように、通電角度指令φ と推定回転子位置θ とを用いて、鎖交磁束インダクタンス交流分ベクトル推定値e Ψsac αβを演算する。
Figure JPOXMLDOC01-appb-M000042
 回転子位置推定誤差演算器1202は、鎖交磁束インダクタンス交流分演算値Ψsac,calc αβと鎖交磁束インダクタンス交流分ベクトル推定値e Ψsac αβとを用いて、回転子位置の推定誤差を演算する。鎖交磁束インダクタンス交流分演算値Ψsac,calc αβと鎖交磁束インダクタンス交流分ベクトル推定値e Ψsac αβとの外積は、通電角度φと通電角度指令φ とが等しいとすると、下記(42)式で表すことができる。なお、鎖交磁束インダクタンス交流分演算値Ψsac,calc αβは、上記(40)式を、回転機電圧指令と回転機電流とを用いて上記(9)式により演算したものである。
Figure JPOXMLDOC01-appb-M000043
 ここで、回転子位置の推定値と真値がおよそ等しい、すなわちθ ≒θとすると、回転子位置の推定誤差は下記(43)式で演算できる。この回転子位置の推定誤差から、実施の形態1と同様にして、推定回転子位置が演算される。
Figure JPOXMLDOC01-appb-M000044
 実施の形態5に係る回転機の制御装置100Dは、鎖交磁束インダクタンス交流分とこれと同方向のベクトルとを用いて、すなわち上記(40)式及び(41)式を用いて回転子位置を推定する。なお、この同方向のベクトルである鎖交磁束インダクタンス交流分ベクトルは、通電角度指令と推定回転子位置とを用いて演算される。これらの鎖交磁束インダクタンス交流分と磁束インダクタンス交流分ベクトルは、二相両方の成分を持っており、鎖交磁束インダクタンス交流分と磁束インダクタンス交流分ベクトルを用いることによって、電流の通電方向に依存せずに回転子位置の推定が可能である。また、実施の形態5に係る回転機の制御装置100Dでは、通電方向によって重み付けを行う必要もない。従って実施の形態5に係る回転機の制御装置100Dは、回転機の通電方向に依らず簡単な構成で回転子位置を推定できる、といった従来にない顕著な効果を奏する。
 実施の形態1から5に係る回転機の制御装置100,100A,100B,100C,100Dが備える各機能は処理回路を用いて実現することができる。各機能とは、電流検出器2、電圧印加器5、位置推定器3,6,8,11,12、制御器4,7,10及び角度演算器9である。
 図13は本発明の実施の形態1から5に係る回転電機の制御装置の第1のハードウェア構成例を示す図である。図13には専用処理回路13のような専用のハードウェアにより上記の処理回路を実現する例が示される。図13に示すように専用のハードウェアを利用する場合、専用処理回路13は、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ASIC(Application Specific Integrated Circuit)、FPGA(Field-Programmable Gate Array)、又はこれらを組み合わせたものが該当する。上記の各機能のそれぞれを、処理回路で実現してもよいし、まとめて処理回路で実現してもよい。
 図14は実施の形態1から5に係る回転電機の制御装置の第2のハードウェア構成例を示す図である。図14にはプロセッサ14及び記憶装置15により上記の処理回路を実現する例が示される。図14に示すようにプロセッサ14及び記憶装置15を利用する場合、上記の各機能のそれぞれは、ソフトウェア、ファームウェア又はこれらの組合せにより実現される。ソフトウェア又はファームウェアはプログラムとして記述され、記憶装置15に記憶される。プロセッサ14は記憶装置15に記憶されたプログラムを読み出して実行する。またこれらのプログラムは、上記の各機能のそれぞれが実行する手順及び方法をコンピュータに実行させるものであるとも言える。記憶装置15は、RAM(Random Access Memory)、ROM(Read Only Memory)、フラッシュメモリー、EPROM(Erasable Programmable Read Only Memory)、又はEEPROM(登録商標)(Electrically Erasable Programmable Read Only Memory)といった半導体メモリが該当する。半導体メモリは不揮発性メモリでもよいし揮発性メモリでもよい。また記憶装置15は、半導体メモリ以外にも、磁気ディスク、フレキシブルディスク、光ディスク、コンパクトディスク、ミニディスク又はDVD(Digital Versatile Disc)が該当する。
 また上記の各機能のそれぞれは、一部をハードウェアで実現し、一部をソフトウェア又はファームウェアで実現してもよい。具体例としては、電流検出器2及び電圧印加器5は専用のハードウェアを用いてその機能が実現され、位置推定器3,6,8,11,12、制御器4,7,10及び角度演算器9はプロセッサ14及び記憶装置15を用いてその機能が実現される。
 実施の形態1から5では、回転機のトルクに対する電流指令は、実効値すなわち回転機の銅損が最小になるように選定されたが、鎖交磁束すなわち回転機の誘起電圧を小さくして回転機をより高速まで運転できるように設定してもよいし、モータ損失を最小化するように設定してもよいことは言うまでもない。また実施の形態1から5では、電流検出器2が回転機に流れる相電流を検出する構成例を説明したが、回転機に流れる相電流は、電圧印加器5であるインバータに内蔵された電流センサ、シャント抵抗等により検出してもよい。また、実施の形態1から5では、説明を分かりやすくするため、制御器における回転機の電圧及び電流の座標変換と、位置推定器における回転機の電圧及び電流の座標変換とが重複して行われているが、これらの座標変換は制御器及び位置推定器に共通の機能で実施してもよい。
 以上の実施の形態に示した構成は、本発明の内容の一例を示すものであり、別の公知の技術と組み合わせることも可能であるし、本発明の要旨を逸脱しない範囲で、構成の一部を省略、変更することも可能である。
 1 回転機、2 電流検出器、3,6,8,11,12 位置推定器、4,7,10 制御器、5 電圧印加器、9 制御座標角度演算器、13 専用処理回路、14 プロセッサ、15 記憶装置、100,100A,100B,100C,100D 制御装置、301,302,601,602,705,1101,1102 三相―二相変換器、303,304,604,802 鎖交磁束インダクタンス交流分演算器、305,606,1202 回転子位置推定誤差演算器、306,607,804 回転子位置演算器、401,701 電流指令演算器、402 三相電流指令演算器、403 三相電流制御器、603 回転座標変換器、706,801,1003,1104 回転座標変換器、605,803 鎖交磁束インダクタンス交流分推定器、702 d-q電流制御器、703,1002 回転座標逆変換器、704 二相-三相変換器、1001 速度演算器、1103 鎖交磁束インダクタンス演算器、1105 適応オブザーバ、1106 適応推定器、1201 鎖交磁束インダクタンス交流分ベクトル推定器。

Claims (13)

  1.  インダクタンスが回転子位置によって変化するインダクタンス交流分を有する回転機を制御する回転機の制御装置であって、
     前記回転機に流れる回転機電流を検出する電流検出器と、
     前記回転子位置を推定する位置推定器と、
     前記回転機電流と前記回転子位置とに基づいて前記回転機を駆動するための回転機電圧指令を出力する制御器と、
     前記回転機電圧指令に基づいて回転機に電圧を印加する電圧印加器とを備え、
     前記位置推定器は、
     前記インダクタンス交流分と前記回転機電流とによって生成される鎖交磁束インダクタンス交流分から前記回転子位置を推定することを特徴とする回転機の制御装置。
  2.  前記回転機のインダクタンスは、前記回転子位置によって変化しない第1の成分と、前記回転子位置に対して2倍の周波数で変化する第2の成分とから構成され、
     前記インダクタンス交流分は、前記第2の成分であることを特徴とする請求項1に記載の回転機の制御装置。
  3.  前記位置推定器は、
     前記回転機電圧指令及び前記回転機電流から演算された前記鎖交磁束インダクタンス交流分である鎖交磁束インダクタンス交流分演算値と、前記回転機電流、前記インダクタンス交流分及び前記回転子位置から推定された前記鎖交磁束インダクタンス交流分である鎖交磁束インダクタンス交流分推定値とを用いて、前記回転子位置を推定することを特徴とする請求項1又は2に記載の回転機の制御装置。
  4.  前記位置推定器は、
     前記鎖交磁束インダクタンス交流分演算値と前記鎖交磁束インダクタンス交流分推定値の外積から前記回転子位置の推定誤差を演算することを特徴とする請求項3に記載の回転機の制御装置。
  5.  前記位置推定器は、
     前記鎖交磁束インダクタンス交流分を前記回転機の回転に同期した回転座標上で演算又は推定することを特徴とする請求項1から4の何れか一項に記載の回転機の制御装置。
  6.  前記位置推定器は、
     前記回転機電圧指令及び前記回転機電流から演算した鎖交磁束インダクタンス交流分演算値と、前記回転機電流と前記回転子位置との角度差である通電角度と、前記回転子位置とを用いて、前記回転子位置の推定を行うことを特徴とする請求項1又は2に記載の回転機の制御装置。
  7.  前記位置推定器は、
     前記通電角度及び前記回転子位置から前記鎖交磁束インダクタンス交流分と同方向のベクトルである鎖交磁束インダクタンス交流分ベクトルを推定し、
     前記鎖交磁束インダクタンス交流分演算値と、推定された前記鎖交磁束インダクタンス交流分ベクトルの推定値との外積から前記回転子位置の推定誤差を演算することを特徴とする請求項6に記載の回転機の制御装置。
  8.  前記位置推定器は、
     オブザーバを用いて鎖交磁束を推定し、推定された当該鎖交磁束インダクタンス交流分の推定誤差を用いて前記回転子位置を推定することを特徴とする請求項1から4の何れか一項に記載の回転機の制御装置。
  9.  前記位置推定器は、
     前記鎖交磁束を前記回転機の回転に同期した回転座標上で演算又は推定することを特徴とする請求項8に記載の回転機の制御装置。
  10.  適応推定誤差をeとし、回転座標上での前記鎖交磁束をΨ dqとし、回転座標上での前記鎖交磁束の推定値をΨ dqとしたとき、
     前記位置推定器は、回転座標上での前記鎖交磁束と、回転座標上での前記鎖交磁束の推定値とを含む下記(1)式によって、前記回転子位置を演算することを特徴とする請求項9に記載の回転機の制御装置。
    Figure JPOXMLDOC01-appb-M000001
  11.  推定した前記回転子位置に制御座標角度を追従させる制御座標角度演算器を備え、
     前記制御器は、
     前記回転機電流を前記制御座標角度で回転座標変換して、前記回転機電圧指令を演算することを特徴とする請求項1から10の何れか一項に記載の回転機の制御装置。
  12.  推定した前記回転子位置に制御座標角度を追従させる制御座標角度演算器を備え、
     前記位置推定器は、
     前記鎖交磁束インダクタンス交流分を前記制御座標角度で変換した回転座標上で演算し、また前記鎖交磁束インダクタンス交流分を前記制御座標角度で変換した回転座標上で推定することを特徴とする請求項5、9又は11に記載の回転機の制御装置。
  13.  前記制御座標角度の追従応答は、前記回転子位置の推定応答よりも遅いことを特徴とする請求項11又は12に記載の回転機の制御装置。
PCT/JP2018/042568 2018-01-12 2018-11-16 回転機の制御装置 WO2019138691A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
DE112018006843.9T DE112018006843T5 (de) 2018-01-12 2018-11-16 Rotationsmaschinen-steuervorrichtung
JP2019564320A JP6847268B2 (ja) 2018-01-12 2018-11-16 回転機の制御装置
US16/959,446 US11296633B2 (en) 2018-01-12 2018-11-16 Rotary machine control device

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2018-003786 2018-01-12
JP2018003786 2018-01-12

Publications (1)

Publication Number Publication Date
WO2019138691A1 true WO2019138691A1 (ja) 2019-07-18

Family

ID=67219506

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2018/042568 WO2019138691A1 (ja) 2018-01-12 2018-11-16 回転機の制御装置

Country Status (4)

Country Link
US (1) US11296633B2 (ja)
JP (1) JP6847268B2 (ja)
DE (1) DE112018006843T5 (ja)
WO (1) WO2019138691A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110572102A (zh) * 2019-10-11 2019-12-13 杭州兆鼎科技实业有限公司 一种电机的软件容错控制方法及其控制系统

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI774315B (zh) * 2021-04-08 2022-08-11 台達電子工業股份有限公司 馬達控制裝置及馬達控制方法
CN115208262A (zh) 2021-04-08 2022-10-18 台达电子工业股份有限公司 马达控制装置及马达控制方法
TWI814596B (zh) * 2022-09-30 2023-09-01 台達電子工業股份有限公司 獲取同步馬達參數的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003500999A (ja) * 1999-05-19 2003-01-07 エービービー インダストリー オサケ ユキチュア 同期機における開ループベクトル制御の始動手順
JP2003189669A (ja) * 2001-12-18 2003-07-04 Switched Reluctance Drives Ltd スイッチトリラクタンス駆動装置の回転子位置検出
JP2004364498A (ja) * 2003-06-04 2004-12-24 Switched Reluctance Drives Ltd スイッチトリラクタンス駆動装置の回転子位置検出

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5144564A (en) * 1991-01-08 1992-09-01 University Of Tennessee Research Corp. Rotor position estimation of a permanent magnet synchronous-machine for high performance drive
KR100421373B1 (ko) 2001-06-20 2004-03-06 엘지전자 주식회사 동기 릴럭턴스 모터의 회전 속도 제어장치
EP1492223A4 (en) 2002-03-22 2005-12-14 Matsushita Electric Ind Co Ltd DEVICE FOR CONTROLLING A SYNCHRONOUS RELUCTANCE ENGINE
US6756753B1 (en) * 2002-12-11 2004-06-29 Emerson Electric Co. Sensorless control system and method for a permanent magnet rotating machine
JP4644010B2 (ja) 2005-02-14 2011-03-02 日本電産株式会社 同期リラクタンス電動機のベクトル制御方法及び同装置
JP5167631B2 (ja) 2006-11-30 2013-03-21 株式会社デンソー モータの制御方法及びそれを利用するモータ制御装置
JP5002335B2 (ja) * 2007-05-29 2012-08-15 株式会社東芝 モータ制御装置,洗濯機及びモータ制御方法
MX2010003444A (es) 2007-09-27 2010-04-21 Mitsubishi Electric Corp Controlador de maquina electrica rotatoria.
KR101313386B1 (ko) 2009-03-25 2013-10-14 미쓰비시덴키 가부시키가이샤 회전 전기 기기의 제어 장치 및 제어 방법
EP2493067B1 (en) 2011-02-22 2018-04-04 ABB Oy Method and apparatus for estimating rotor angle of synchronous reluctance motor
WO2016121751A1 (ja) 2015-01-28 2016-08-04 株式会社 東芝 インバータ制御装置及びモータ駆動システム

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003500999A (ja) * 1999-05-19 2003-01-07 エービービー インダストリー オサケ ユキチュア 同期機における開ループベクトル制御の始動手順
JP2003189669A (ja) * 2001-12-18 2003-07-04 Switched Reluctance Drives Ltd スイッチトリラクタンス駆動装置の回転子位置検出
JP2004364498A (ja) * 2003-06-04 2004-12-24 Switched Reluctance Drives Ltd スイッチトリラクタンス駆動装置の回転子位置検出

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110572102A (zh) * 2019-10-11 2019-12-13 杭州兆鼎科技实业有限公司 一种电机的软件容错控制方法及其控制系统

Also Published As

Publication number Publication date
US11296633B2 (en) 2022-04-05
DE112018006843T5 (de) 2020-10-01
US20200373864A1 (en) 2020-11-26
JPWO2019138691A1 (ja) 2020-04-23
JP6847268B2 (ja) 2021-03-24

Similar Documents

Publication Publication Date Title
WO2019138691A1 (ja) 回転機の制御装置
JP6279151B2 (ja) 交流回転機の制御装置および電動パワーステアリング装置
JP6367332B2 (ja) インバータ制御装置及びモータ駆動システム
JP6375757B2 (ja) 電動機制御装置、電動機の磁束推定装置および電動機の磁束推定方法
JP5929863B2 (ja) 制御装置
TWI545885B (zh) 交流旋轉機的控制裝置
JP6027176B2 (ja) 同期電動機のセンサレスベクトル制御のための回転角推定装置
JP6166601B2 (ja) モータ制御装置及び発電機制御装置
US11515821B2 (en) Control device for rotating electrical machine
JP2010035352A (ja) 同期電動機のロータ位置推定装置
WO2015186406A1 (ja) 周期外乱自動抑制装置
JP6847269B2 (ja) 回転機の制御装置
JP4359546B2 (ja) 交流モータの制御装置
JP6766398B2 (ja) 磁石温度推定方法、及び、磁石温度推定装置
JP5660191B2 (ja) 電動機制御装置
JP7199605B1 (ja) 回転機制御装置
JP7362003B2 (ja) 回転機の制御装置
WO2023238219A1 (ja) 回転機の制御装置
JP7357112B2 (ja) 回転機制御装置
WO2023026352A1 (ja) 回転機の制御装置
JP7166443B2 (ja) 交流回転機の制御装置及び交流回転機の制御方法
JP2897373B2 (ja) 直流ブラシレスモータの制御装置
JP2022131542A (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: 18899378

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2019564320

Country of ref document: JP

Kind code of ref document: A

122 Ep: pct application non-entry in european phase

Ref document number: 18899378

Country of ref document: EP

Kind code of ref document: A1