WO2020148900A1 - 回転電機の制御装置 - Google Patents

回転電機の制御装置 Download PDF

Info

Publication number
WO2020148900A1
WO2020148900A1 PCT/JP2019/001516 JP2019001516W WO2020148900A1 WO 2020148900 A1 WO2020148900 A1 WO 2020148900A1 JP 2019001516 W JP2019001516 W JP 2019001516W WO 2020148900 A1 WO2020148900 A1 WO 2020148900A1
Authority
WO
WIPO (PCT)
Prior art keywords
estimated
electric machine
rotor position
speed
rotating electric
Prior art date
Application number
PCT/JP2019/001516
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 PCT/JP2019/001516 priority Critical patent/WO2020148900A1/ja
Priority to DE112019006688.9T priority patent/DE112019006688T5/de
Priority to JP2020566081A priority patent/JP7019077B2/ja
Publication of WO2020148900A1 publication Critical patent/WO2020148900A1/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
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/24Vector control not involving the use of rotor position or rotor speed sensors
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/14Electronic commutators
    • H02P6/16Circuit arrangements for detecting position
    • H02P6/18Circuit arrangements for detecting position without separate position detecting elements
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/14Electronic commutators
    • H02P6/16Circuit arrangements for detecting position
    • H02P6/18Circuit arrangements for detecting position without separate position detecting elements
    • H02P6/181Circuit arrangements for detecting position without separate position detecting elements using different methods depending on the speed

Definitions

  • the present invention relates to a rotary electric machine control device that controls a rotary electric machine by obtaining rotor position information without using a position sensor that detects the rotational position of the rotor.
  • a position sensor for acquiring the rotor position is used in the control device for the rotating electric machine.
  • the use of the position sensor causes problems such as an increase in size of the system, an increase in cost, and deterioration of environmental resistance. Therefore, it is required for the control device for the rotating electric machine to apply the position sensorless control for driving the rotating electric machine without using the position sensor.
  • the position sensorless control is roughly divided into a method of estimating the rotor position by utilizing the salient pole of the rotor and a method of estimating the rotor position by using the interlinkage magnetic flux calculated by the induced voltage generated in the rotating electric machine.
  • the former is referred to as a "salient pole method” and the latter is referred to as an "induced voltage and interlinkage magnetic flux method”.
  • the saliency of the rotor is magnetic anisotropy of the inductance of the rotor, and is a characteristic that the inductance changes depending on the rotor position. Therefore, it can be said that the salient pole method is a method that utilizes the angle dependence of the inductance.
  • the salient pole method information about salient polarity is excited by superimposing a position estimation voltage or position estimation current on a rotating electric machine, and the rotor position is estimated based on the excited information.
  • the salient pole method is used in a low speed region where the induced voltage required for position estimation is not sufficiently obtained, and the induced voltage and interlinkage magnetic flux method is used in a high speed region where the induced voltage is sufficiently obtained.
  • the position sensorless control in the torque speed range is realized by switching between the two position sensorless control methods according to the torque speed range.
  • Patent Document 1 discloses a method of estimating the rotor position using the anisotropy of magnetic saturation. This method is called “magnetic saturation method".
  • the position estimation voltage is superimposed on each of the dq axes of the rotation coordinate, and the multiplication value of the d axis amplitude of the position estimation current generated by the superposition of the position estimation voltage and the q axis amplitude is added. The position estimation is performed based on this.
  • Patent Document 1 has a salient-pole position estimator and a magnetic saturation position estimator.
  • the position estimation error correlation amount ⁇ 1 is calculated by the salient pole method
  • the position estimation error correlation amount ⁇ 2 is calculated by the magnetic saturation method.
  • the rotor position is estimated using the position estimation error correlation amount ⁇ 1
  • the rotor position is estimated using the position estimation error correlation amount ⁇ 2.
  • the contributions of the position estimation error correlation amount ⁇ 1 and the position estimation error correlation amount ⁇ 2 are weighted averaged according to the degree of magnetic saturation, and the rotational position is estimated based on the weighted average contribution.
  • Patent Document 1 does not include a position estimator for driving in a high-speed area. For this reason, with the technique of Patent Document 1, it is difficult to perform position sensorless control in the high speed region.
  • the present invention has been made in view of the above, and an object of the present invention is to obtain a controller for a rotary electric machine that can perform position sensorless control covering an expanded torque speed range.
  • the present invention is a controller for a rotating electric machine that controls the driving of a multi-phase rotating electric machine, wherein the current detecting means detects a rotating electric machine current flowing through the rotating electric machine.
  • a position estimating means for estimating a rotor position which is a rotational position of a rotor of the rotating electric machine based on the rotating electric machine current; and a rotation based on the rotating electric machine current and an estimated rotor position which is an estimated value of the rotor position.
  • a controller is provided that calculates a rotating electric machine drive voltage command for driving the electric machine and a position estimation voltage command for each phase for estimating the rotor position.
  • the position estimating means includes a first estimated rotor position based on an AC component of a current amplitude for position estimation which is an amplitude of a current for position estimation generated by applying a voltage for position estimation based on a voltage command for position estimation, and for position estimation.
  • a second estimated rotor position based on the DC component of the current amplitude, a rotating electrical machine current, and a third estimated rotor position based on the interlinkage magnetic flux calculated by the induced voltage caused by the rotor saliency.
  • the position estimation means obtains an estimated rotor position obtained by selecting any one of the first estimated rotor position, the second estimated rotor position, and the third estimated rotor position, or An estimated rotor position obtained by combining using at least two of the first estimated rotor position, the second estimated rotor position, and the third estimated rotor position is output.
  • control device for a rotating electric machine of the present invention it is possible to perform position sensorless control covering the expanded torque speed range.
  • the figure which shows the structural example of the detail of the position estimator shown in FIG. FIG. 1 is a first diagram for explaining the operation of the second estimated position calculator shown in FIG.
  • the figure which shows the structural example of the detail of the 3rd estimated position calculator shown in FIG. The figure which shows the structural example of the control apparatus of the rotary electric machine which concerns on Embodiment 2.
  • FIG. 8 is a diagram for explaining a drive region in the estimated position switch shown in FIGS. 8 and 12.
  • the figure which shows the structural example of the detail of the position estimator shown in FIG. The figure which shows the example of a detailed structure of the estimated position switching device shown in FIG. FIG.
  • FIG. 15 is a diagram for explaining a drive region in the estimated position switch shown in FIGS. 15 and 16.
  • the figure which shows the structural example of the detail of the position estimator shown in FIG. The figure which shows the structural example of the detail of the estimated position switching device shown in FIG. FIG. 19 is a diagram provided for explaining a drive region in the estimated position switch shown in FIGS. 19 and 20.
  • the figure which shows the structural example of the detail of the estimated position switching device shown in FIG. The figure which shows the 1st hardware structural example of the control apparatus of the rotary electric machine which concerns on Embodiment 1-5.
  • control device for a rotating electric machine will be described in detail with reference to the drawings.
  • the present invention is not limited to the embodiments described below.
  • control device for the rotating electric machine may be simply referred to as the “control device”.
  • the control device for a rotary electric machine is a control device that controls driving of a multi-phase rotary electric machine.
  • the multi-phase rotating electric machine is a rotating electric machine configured to be able to apply an alternating voltage of three phases or four phases or more.
  • FIG. 1 is a diagram showing a configuration example of a control device 100 for a rotary electric machine according to the first embodiment.
  • the control device 100 according to the first embodiment includes a voltage applicator 3 that is a voltage applying unit, a current detector 2 that is a current detecting unit, and a position estimator 4 that is a position estimating unit. , And a controller 5 as a control means.
  • the voltage applicator 3 supplies AC power to the rotary electric machine 1 based on the rotary electric machine voltage commands v u * , v v * , v w * for driving the rotary electric machine 1.
  • An example of the rotating electric machine 1 is a three-phase synchronous reluctance motor that generates torque by utilizing the salient pole of the rotor 1a.
  • the synchronous reluctance motor is a motor in which the magnetic resistance of the rotor 1a changes depending on the rotor position.
  • the current detector 2 detects rotating electric machine currents i u , iv , and i w that are currents flowing through the rotating electric machine 1.
  • the rotating electric machine currents i u , iv , and i w are alternating currents supplied from the voltage applying device 3 to each phase of the rotating electric machine 1.
  • the current detector 2 outputs the detected rotating electrical machine currents iu , iv , and iw to the position estimator 4 and the controller 5.
  • the position estimator 4 calculates an estimated rotor position ⁇ ⁇ based on the rotary electric machine currents i u , i v , i w .
  • the controller 5 controls the rotating electric machine 1 so that the output torque of the rotating electric machine 1 becomes a value instructed by the torque command value ⁇ * based on the rotating electric machine currents i u , i v , i w and the estimated rotor position ⁇ ⁇ .
  • the voltage commands v u * , v v * , v w * are calculated.
  • the estimated rotor position ⁇ ⁇ is an estimated value of the rotor position, which is the rotational position of the rotor 1a included in the rotary electric machine 1.
  • the estimated rotor position ⁇ ⁇ is represented by an electrical angle.
  • the controller 5 is divided into a drive voltage command calculation unit 5a which is a first calculation unit and a position estimation voltage calculation unit 5b which is a second calculation unit.
  • the drive voltage command calculator 5a includes a current command calculator 6, a current controller 7, a rotation coordinate inverse converter 8, a two-phase/three-phase converter 9, a drive current extractor 11, a three-phase/two-phase converter 12, and a rotation coordinate.
  • the converter 13 and the adder 14 are provided.
  • the division of the drive voltage command calculation unit 5a and the position estimation voltage calculation unit 5b shown in FIG. 1 is an example, and the constituent elements of the controller 5 may be divided in any manner.
  • the position estimating voltage calculator 5b calculates position estimating voltage commands v uf * , v vf * , v wf * for estimating the rotor position.
  • Driving voltage calculation unit 5a the rotary electric machine currents i u, i v, i w and rotation for driving the rotary electric machine 1 on the basis of the estimated rotor position theta ⁇ machine driving voltage command v uf *, v vf *, Calculate v wf * .
  • the rotating electric machine driving voltage command v uf *, v vf *, v wf * is, the adder 14 in the position estimation voltage instruction v uh *, v vh *, v wh * is added.
  • the output of the adder 14 is output to the voltage applicator 3 as the rotating electrical machine voltage commands v u * , v v * , v w * .
  • the torque command value ⁇ * which is a command value of the output torque of the rotary electric machine 1, is input to the current command calculator 6.
  • the current command calculator 6 uses the torque command value ⁇ * to calculate the rotary electric machine drive current commands i df * , i qf * on the rotary two-phase coordinates.
  • the rotary electric machine drive current commands i df * , i qf * are current commands on the rotary two-phase coordinates necessary for the rotary electric machine 1 to generate an output corresponding to the torque command value ⁇ * .
  • the current command calculator 6 of the first embodiment calculates a current command that minimizes the effective current value for torque, in other words, minimizes copper loss for torque.
  • the rotary electric machine drive current command i df * is the d-axis direction armature current component that minimizes the magnetic resistance of the rotor 1 a. It is a command value of the d-axis drive current shown.
  • the rotary electric machine drive current command i qf * is a command value of the q-axis drive current indicating the armature current component in the q-axis direction, which is the direction orthogonal to the d-axis.
  • the motor constant of the rotary electric machine 1 is used to calculate the rotary electric machine drive current commands i df * , i qf * on the two-phase rotary coordinates.
  • the motor constant include the mutual inductance of the rotary electric machine 1 and the number of poles of the rotary electric machine 1.
  • a relational expression or a table between the current command and the torque obtained in advance may be used.
  • the drive current extractor 11 uses the rotary electric machine currents i u , i v , i w on the three-phase coordinates detected by the current detector 2 to calculate the rotary electric machine drive currents i uf , i vf , i wf on the three-phase coordinates. To extract.
  • the rotary electric machine drive currents i uf , i vf , i wf on the three-phase coordinates are generated by the rotary electric machine drive voltage commands v uf * , v vf * , v wf * for driving the rotary electric machine 1 on the three-phase coordinates. This is the rotating electric machine drive current on the three-phase coordinates.
  • the rotary electric machine drive voltage commands v uf * , v vf * , v wf * are drive voltage commands output from the two-phase/three-phase converter 9 and input to the adder 14.
  • the adder 14 includes rotating electric machine drive voltage commands v uf * , v vf * , v wf * and position estimating voltage commands v uh * , v that are voltage commands for estimating the rotor position of the rotating electric machine 1.
  • vh * and vwh * are input.
  • the adder 14 outputs the rotary electric machine drive voltage commands v uf * , v vf * , v wf * on the three-phase coordinates and the position estimation voltage commands v uh * , v vh * , v wh * on the three-phase coordinates.
  • the rotating electric machine voltage commands v u * , v v * , v w * on the three-phase coordinates are generated and added to the voltage applying device 3.
  • the position estimating voltage commands v uh * , v vh * , v wh * are calculated by the position estimating voltage calculator 5b.
  • FIG. 2 is a diagram showing a position estimation voltage command output from the position estimation voltage calculation unit 5b shown in FIG.
  • the position estimation voltage commands v uh * , v vh * , and v wh * on the three-phase coordinates square wave-shaped voltages having a phase difference of 120° from each other are applied to the u-phase from the upper side. , V phase and w phase are shown in this order.
  • FIG. 2 shows a case where the position estimation voltage commands v uh * , v vh * , v wh * on the three-phase coordinates are square wave voltages, the present invention is not limited to this.
  • a sine wave voltage may be used instead of the square wave voltage.
  • the position estimation voltage commands v uh * , v vh are calculated from the rotary electric machine currents i u , i v , i w on the three-phase coordinates by a notch filter, for example.
  • a notch filter for example.
  • the rotary electric machine currents i u of the three-phase coordinate, i v, rotary machine driving current on the three-phase coordinates from i w i uf, i vf, a method of extracting i wf is not limited to the notch filter, low pass A filter or a high pass filter may be used.
  • the three-phase/two-phase converter 12 converts the rotary electric machine drive currents i uf , i vf , i wf on the three-phase coordinates extracted by the drive current extractor 11 into the rotary electric machine drive currents i ⁇ f , i on the stationary two-phase coordinates. Convert to ⁇ f .
  • the rotary coordinate converter 13 performs coordinate conversion using the estimated rotor position ⁇ ⁇ estimated by the position estimator 4 to convert the rotary electric machine drive currents i ⁇ f and i ⁇ f on the stationary two-phase coordinates to the rotary two-phase coordinates. Convert to the above rotating electric machine drive currents i df and i qf .
  • the rotary electric machine drive currents i df and i qf on the rotary two-phase coordinates converted by the rotary coordinate converter 13 are calculated by the current command calculator 6, and the rotary electric machine drive current commands i df * , I qf * is controlled so that the rotating electric machine drive voltage commands v df * , v qf * on the rotating two-phase coordinates are calculated.
  • An example of the current control in the current controller 7 is proportional integral (PI) control.
  • the rotary coordinate inverse converter 8 uses the estimated rotor position ⁇ ⁇ to convert the rotary electric machine drive voltage commands v df * , v qf * on the rotary two-phase coordinates calculated by the current controller 7 to the stationary two-phase coordinates.
  • the above rotary electric machine drive voltage command is converted into v ⁇ f * and v ⁇ f * .
  • the two-phase/three-phase converter 9 converts the rotary electric machine drive voltage commands v ⁇ f * , v ⁇ f * on the stationary two-phase coordinates to the rotary electric machine drive voltage commands v uf * , v vf * , v on the three-phase coordinates described above. Convert to wf * .
  • FIG. 3 is a diagram showing a detailed configuration example of the position estimator 4 shown in FIG.
  • the position estimator 4 includes a signal processor 41, a position estimation current amplitude calculator 403, a signal processor 42, an estimated position calculator 406 which is a first estimated position calculator, and a second estimated position calculator 406.
  • An estimated position calculator 407 that is an estimated position calculator, an estimated position calculator 408 that is a third estimated position calculator, and an estimated position switch 409 are provided.
  • the signal processor 41 includes a position estimation current extractor 401 and a drive current extractor 402.
  • the position estimation current extractor 401 extracts the position estimation currents i uh , i vh , i wh from the rotary electric machine currents i u , i v , i w on the three-phase coordinates.
  • the drive current extractor 402 extracts the rotary electric machine drive currents i uf , i vf , i wf on the three-phase coordinates from the rotary electric machine currents i u , i v , i w on the three-phase coordinates.
  • the position estimation current amplitude calculator 403 calculates the position estimation current amplitudes I uh , I vh , I wh on the three-phase coordinates based on the position estimation currents i uh , i vh , i wh on the three-phase coordinates. Calculate
  • the signal processor 42 includes an AC component extractor 404 and a DC component extractor 405.
  • the AC component extractor 404 extracts the three-phase AC components I uhac , I vhac , and I whac of the current amplitude for position estimation from the current amplitudes I uh , I vh , and I wh for position estimation.
  • the estimated position calculator 406 estimates the first estimated rotor position ⁇ ⁇ 1 based on the three-phase AC components I uhac , I vhac , and I whac of the position estimation current amplitude.
  • the rotor position estimation method used in the estimated position calculator 406 is a salient pole method. That is, the first estimated rotor position ⁇ ⁇ 1 is an estimated rotor position calculated based on the AC component of the position estimation current amplitude, which is the amplitude of the position estimation current generated by the application of the position estimation voltage. ..
  • the first estimated rotor position ⁇ ⁇ 1 is one of the candidate values of the estimated rotor position ⁇ ⁇ that is finally output from the position estimator 4.
  • the estimated position calculator 407 calculates the second estimated rotor position ⁇ ⁇ 2 based on the DC component I hdc of the position estimation current amplitude.
  • the rotor position estimation method used in the estimated position calculator 407 is a magnetic saturation method.
  • the second estimated rotor position ⁇ ⁇ 2 is one of the candidate values of the estimated rotor position ⁇ ⁇ that is finally output from the position estimator 4.
  • Estimated position calculator 408 the rotary electric machine drive current i uf on the three-phase coordinate, i vf, i wf and a flux linkage that is calculated by the induced voltage generated by the saliency of the rotor 1a, the third on the basis of the The estimated rotor position ⁇ ⁇ 3 of is calculated.
  • the estimation method of the rotor position used in the estimated position calculator 408 is an induced voltage and flux linkage method.
  • the third estimated rotor position ⁇ ⁇ 3 is one of the candidate values of the estimated rotor position ⁇ ⁇ finally output from the position estimator 4.
  • the estimated position switch 409 selects any one of the first estimated rotor position ⁇ ⁇ 1 , the second estimated rotor position ⁇ ⁇ 2, and the third estimated rotor position ⁇ ⁇ 3.
  • the position information obtained in step 1 or the position information obtained by combining using at least two estimated rotor positions is output as the estimated rotor position ⁇ ⁇ .
  • FIG. 4 is a first diagram for explaining the operation of the estimated position calculator 407 shown in FIG.
  • FIG. 5 is a second diagram for explaining the operation of the estimated position calculator 407 shown in FIG.
  • FIG. 6 is a diagram showing a detailed configuration example of the estimated position calculator 408 shown in FIG.
  • Drive current extractor 402 uses, for example, a notch filter, the rotary electric machine currents i u of the three-phase coordinate, i v, from i w, position estimation voltage command v uh *, v vh *, v wh * application of By removing the position estimation current amplitudes I uh , I vh , and I wh on the three-phase coordinates, the rotary electric machine drive currents i uf , i vf , and i wf on the three-phase coordinates are extracted.
  • Location estimation current extractor 401 for example, the rotary electric machine currents i u of the three-phase coordinate, i v, from i w, rotation of the three-phase coordinate calculated in the drive current extractor 402 machine driving current i uf, i vf, by subtracting the i wf, location estimation current amplitude I uh in a three-phase coordinates, I vh, calculates the I wh.
  • the method for extracting or calculating the rotary electric machine drive currents i uf , i vf , i wf on the three-phase coordinates and the position estimation current amplitudes I uh , I vh , I wh on the three-phase coordinates is not limited to this.
  • a bandpass filter, a bandstop filter, a lowpass filter or a highpass filter may be used.
  • the position estimation current amplitudes I uh , I vh , and I wh on the three-phase coordinates are calculated by using the above-described DC component I hdc of the position estimation current amplitude and the AC component I hac of the position estimation current amplitude. , Can be expressed as in the following equation (1).
  • the position estimation current amplitude I uh in a three-phase coordinates, I vh are those utilizing the relative relationship of the AC component of the I wh computes the magnetic pole position, the position estimation current amplitude I uh , I vh , I wh need not be calculated. Therefore, the following equation (2) can be used to calculate the position estimation current amplitudes I uh , I vh , and I wh on the three-phase coordinates.
  • the DC component extractor 405 extracts the DC component I hdc of the position estimation current amplitude from the position estimation current amplitudes I uh , I vh , and I wh . Specifically, the DC component extractor 405 calculates the position estimation current amplitudes I uh , I vh on the three-phase coordinates calculated by the position estimation current amplitude calculator 403 as shown in the following equation (3). By calculating the average value of Iwh , the DC component Ihdc of the current amplitude for position estimation is extracted.
  • the AC component extractor 404 extracts the three-phase AC components I uhac , I vhac , and I whac of the position estimation current amplitude from the position estimation current amplitudes I uh , I vh , and I wh .
  • the AC component extractor 404 extracts a DC component I hdc of the position estimation current amplitude extracted by the DC component extractor 405 from each of the input position estimation current amplitudes I uh , I vh , and I wh.
  • the three-phase AC components I uhac , I vhac , and I whac of the position estimation current amplitude are calculated.
  • the three-phase AC components I uhac , I vhac , and I whac of the current amplitude for position estimation can be expressed by the following expression (4).
  • DC components I hdc and position estimation current amplitudes of the three-phase AC component I Uhac position estimation current amplitude, I Vhac, a method of calculating the I Whac is not limited thereto, using a low-pass filter or high-pass filter Good.
  • the estimated position calculator 406 utilizes that the three-phase AC components I uhac , I vhac , and I whac of the position estimation current amplitude are functions of the rotor position ⁇ , that is, sin(2 ⁇ ) or cos(2 ⁇ ). Then, the first estimated rotor position ⁇ ⁇ 1 is calculated. Specifically, the first estimated rotor position ⁇ ⁇ 1 is obtained by inverting the signal of any one of the three-phase AC components I uhac , I vhac , and I whac of the current amplitude for position estimation shown in the above equation (4). It can be obtained by calculating the cosine.
  • the first estimated rotor position ⁇ ⁇ 1 is reversed after three-phase two-phase conversion of the three-phase AC components I uhac , I vhac , and I whac of the position estimation current amplitude expressed on the three-phase coordinates. It can be obtained by calculating the tangent.
  • the first estimated rotor position ⁇ ⁇ 1 has six sections for each 60° electrical angle about the zero-cross points of the three-phase AC components I uhac , I vhac , and I whac of the position estimation current amplitude. It is also possible to perform the calculation by linearly approximating the zero-crossing one of the three-phase AC components I uhac , I vhac , and I whac of the position estimation current amplitude in each section.
  • the estimated position calculator 407 calculates the second estimated rotor position ⁇ ⁇ 2 based on the DC component I hdc of the position estimation current amplitude.
  • FIG. 4 the behavior of the DC component I hdc of the current amplitude for position estimation in maximum torque (Maximum Torque Per Ampere: MTPA) control is shown in relation to the current phase ⁇ and the magnitude
  • the MTPA control is control that minimizes the magnitude
  • of the drive current vector is d when the rotary electric machine drive currents i uf , i vf , and i wf on the three-phase coordinates are subjected to rotary coordinate conversion and divided into a d-axis component and a q-axis component. It corresponds to the square root of the sum of squares of the axis component and the q-axis component.
  • is given in 20% current increments between 20% and 100% of the rated current.
  • the current phase ⁇ is defined as a lead phase with respect to the d axis.
  • the control range of current phase ⁇ is assumed to be 45° to 55°.
  • of the drive current vector is 60% or more of the rated current and the current phase ⁇ is in the range of 45° to 55°
  • the DC component I hdc of the current amplitude for position estimation is , And monotonically decreases with respect to the current phase ⁇ .
  • the estimated position calculator 407 in Embodiment 1 utilizes this characteristic to calculate the second estimated rotor position ⁇ ⁇ 2 .
  • the characteristic required for estimation of rotor position is unique with the DC component I hdc and the driving current vector position estimation current amplitude, position estimation current amplitude DC component I hdc and the current phase ⁇
  • the relationship is not limited to monotonically decreasing.
  • FIG. 5 shows, the magnitude of the drive current vector shown in FIG. 4
  • the DC component Ihdc of the working current amplitude is shown.
  • the current phase command value is ⁇ *
  • the DC component I hdc of the position estimation current amplitude when driven with the current phase command value ⁇ * is expressed as I hdc ( ⁇ * ).
  • the position estimation error of the rotor position is ⁇ 2
  • the DC component I hdc of the position estimation current amplitude when the current phase command value ⁇ * is driven by the actual current phase ⁇ deviated by the position estimation error ⁇ 2 is calculated. It is defined as I hdc ( ⁇ ).
  • the position estimation error ⁇ 2 is proportional to the difference between I hdc ( ⁇ * ) and I hdc ( ⁇ ). In this specification, the difference between I hdc ( ⁇ * ) and I hdc ( ⁇ ) is referred to as “Ihdc error”.
  • the estimated position calculator 407 When the estimated position calculator 407 is composed of a phase locked loop (PLL), the estimated position calculator 407 operates the PLL so that the Ihdc error becomes zero, and thus the second estimated rotor position is calculated. Calculate ⁇ ⁇ 2 .
  • the PLL may have any configuration as long as the Ihdc error becomes zero, and a proportional integrator or a proportional integral integrator is exemplified.
  • the proportional-plus-integrator has a configuration in which an integrator is further provided after the proportional-plus-integrator.
  • the relationship between the current phase ⁇ and the DC component I hdc of the position estimation current amplitude is stored in advance in association with the assumed torque or the range of the rotating electric machine current. To do.
  • the estimated position calculator 408 includes a position estimation error calculator 4080 and a PLL 4081, as shown in FIG. As described above, the estimated position calculator 408 is a calculator that estimates the third estimated rotor position ⁇ ⁇ 3 based on the interlinkage magnetic flux calculated by the induced voltage generated by the saliency of the rotor 1a. ..
  • the position estimation error calculator 4080 outputs the rotary electric machine drive voltage commands v uf * , v vf * , v wf * on the three-phase coordinates and the rotary electric machine drive currents i uf , i vf , i wf on the three-phase coordinates. Based on this, the position estimation error “ ⁇ ( ⁇ ⁇ 3 ⁇ )” of the rotor position is calculated.
  • the PLL 4081 calculates the third estimated rotor position ⁇ ⁇ 3 from the position estimation error “ ⁇ ( ⁇ ⁇ 3 ⁇ )” of the rotor position.
  • the position estimation error calculator 4080 is a three-phase/two-phase converter 40800, 40801, a rotating coordinate converter 40802, an interlinkage magnetic flux inductance AC component calculator 40803, an interlinkage magnetic flux inductance AC component estimator 40804, and a rotor position estimation error calculator. 40805 is provided.
  • the three-phase/two-phase converter 40800 converts the rotary electric machine drive voltage commands v uf * , v vf * , v wf * on the three-phase coordinates to the rotary electric machine drive voltage commands v ⁇ f * , v ⁇ f * on the stationary two-phase coordinates . Convert to.
  • Three-phase two-phase converter 40 801 the rotary electric machine drive current i uf on the three-phase coordinate, i vf, the i wf, rotary machine driving current i .alpha.f on the stationary two-phase coordinates converted to i .beta.f.
  • the rotary coordinate converter 40802 performs coordinate conversion using the estimated rotor position ⁇ ⁇ 3 , and outputs the rotary electric machine drive currents i ⁇ f and i ⁇ f on the stationary two-phase coordinates to the rotary electric machine drive currents i df on the rotary two-phase coordinates. , I qf .
  • the PLL 4081 calculates the third estimated rotor position ⁇ ⁇ 3 by performing the PLL operation so that the position estimation error “ ⁇ ( ⁇ ⁇ 3 ⁇ )” of the rotor position becomes zero.
  • the PLL 4081 may be configured so that the position estimation error “ ⁇ ( ⁇ ⁇ 3 ⁇ )” becomes zero, and a proportional integrator or a proportional integral differentiator may be used.
  • v dq is a vector composed of the d-axis rotary electric machine voltage v d and the q-axis rotary electric machine voltage v q
  • i dq is the d-axis rotary electric machine current i d and the q-axis rotary electric machine current i q.
  • R s is the winding resistance of the rotary electric machine 1
  • ⁇ s is the rotational angular velocity of the coordinates representing the model
  • ⁇ dq in the above equation (5) is the interlinkage magnetic flux.
  • J in the above equation (6) is a transformation matrix.
  • the interlinkage magnetic flux ⁇ dq in the above equation (5) can be expressed as in the above equation (7).
  • L dq in the above equation (7) can be expressed by a matrix as in the above equation (8) by using L sdc , L mac, and the electrical angle ⁇ of the rotor position.
  • L sdc in the above equation (8) is an inductance DC component that does not change depending on the rotor position
  • L mac is an inductance AC component that changes depending on the rotor position.
  • the change in inductance is generally represented by a sine function or cosine function of 2 ⁇ with respect to the electrical angle ⁇ of the rotor position.
  • the first term of the above equation (9) is the term of the inductance DC component L sdc that does not change depending on the rotor position.
  • the second term of the equation (9) is a term based on the inductance AC component L mac that changes depending on the rotor position, and this term is the “linkage flux inductance AC component”. That is, the interlinkage magnetic flux inductance AC component is the interlinkage magnetic flux generated by the inductance AC component and the rotating electric machine current.
  • the interlinkage magnetic flux inductance AC component calculator 40803 performs the following computations in order to compute the interlinkage magnetic flux inductance AC component. First, the interlinkage magnetic flux inductance AC component calculator 40803 calculates the interlinkage magnetic flux ⁇ dq of the rotating electric machine by using the following equation (10).
  • v ⁇ * is a vector composed of an ⁇ -axis rotating electric machine voltage command v ⁇ * and a ⁇ -axis rotating electric machine voltage command v ⁇ * .
  • the transfer function of the high-pass filter is represented by the following equation (12), where the cutoff frequency is ⁇ hpf .
  • the interlinkage magnetic flux inductance AC component calculator 40803 uses the third estimated rotor position ⁇ ⁇ 3 to calculate the interlinkage flux ⁇ ⁇ hpf ⁇ on the stationary two-phase coordinate to the chain on the rotating two-phase coordinate. Coordinate conversion is performed into the magnetic flux ⁇ ⁇ hpfdq .
  • the interlinkage magnetic flux inductance AC components ⁇ ⁇ acdq , calc on the rotating coordinates are calculated by the following equation (15) according to the above equation (9).
  • the interlinkage magnetic flux inductance alternating current component ⁇ ⁇ addq , calc calculated by the above equation (15) is hereinafter referred to as a “linkage magnetic flux inductance alternating current component calculation value”.
  • the interlinkage magnetic flux inductance AC component estimator 40804 is the second term of the above equation (9) using the estimated rotor position ⁇ ⁇ 3 and the rotating electrical machine current i dq as shown in the following equation (16).
  • the flux linkage inductance AC component is estimated.
  • ⁇ ⁇ acdq in the above equation (17) is an estimated value of the interlinkage magnetic flux inductance AC component calculated by the interlinkage magnetic flux inductance AC component estimator 40804.
  • this estimated value is referred to as an “interlinkage magnetic flux inductance AC component estimated value”.
  • the position estimation error calculator 4080 uses the interlinkage magnetic flux inductance AC component calculated values ⁇ ⁇ acdq , calc and the interlinkage magnetic flux inductance AC component estimated value ⁇ ⁇ acdq to estimate the rotor position position error "-( ⁇ ⁇ 3 ⁇ ⁇ )” is calculated.
  • the cross product of the interlinkage magnetic flux inductance AC component calculation values ⁇ ⁇ acdq , calc and the interlinkage magnetic flux inductance AC component estimated value ⁇ ⁇ acdq is expressed by the above equation (15), that is, the second term of the above equation (9).
  • the calculated value and the above equation (16) it is represented by the following equation (18).
  • the estimation error of the rotor position can be calculated by the following equation (19).
  • the above is the calculation processing by the estimated position calculator 408.
  • the rotary electric machine voltage command and the rotary electric machine current used for position estimation are the rotary electric machine drive voltage command and the rotary electric machine drive current, respectively.
  • the estimated position switch 409 selects one of the first estimated rotor position ⁇ ⁇ 1 , the second estimated rotor position ⁇ ⁇ 2, and the third estimated rotor position ⁇ ⁇ 3. Select one, that is, switch and output.
  • the estimated position switcher 409 uses at least two estimated rotors among the first estimated rotor position ⁇ ⁇ 1 , the second estimated rotor position ⁇ ⁇ 2, and the third estimated rotor position ⁇ ⁇ 3.
  • the position information is selected, and a value obtained by weighting and averaging these at a preset ratio is output as the estimated rotor position ⁇ ⁇ .
  • the estimated position switch 409 selects or switches the estimated rotor position information and outputs it.
  • the position estimator 4 uses the first estimated rotor position ⁇ ⁇ 1 , the second estimated rotor position ⁇ ⁇ 2, and the third estimated rotor position ⁇ ⁇ . selects and outputs any one of the three.
  • the position estimator 4 uses at least two estimated rotor positions of the first estimated rotor position ⁇ ⁇ 1 , the second estimated rotor position ⁇ ⁇ 2, and the third estimated rotor position ⁇ ⁇ 3.
  • the estimated rotor position ⁇ ⁇ synthesized using is output.
  • Embodiment 2 In the first embodiment, by selecting any one of the first estimated rotor position ⁇ ⁇ 1 , the second estimated rotor position ⁇ ⁇ 2, and the third estimated rotor position ⁇ ⁇ 3. In the embodiment, the obtained position information or the position information obtained by combining at least two estimated rotor positions is output as the estimated rotor position ⁇ ⁇ .
  • the first estimated rotor position ⁇ ⁇ 1 is position information estimated by the salient pole method
  • the second estimated rotor position ⁇ ⁇ 2 is position information estimated by the magnetic saturation method.
  • the third estimated rotor position ⁇ ⁇ 3 is position information estimated by the induced voltage and the flux linkage method.
  • each rotor position estimation method has its own characteristics.
  • the salient pole method has a high estimation accuracy in the low-speed range and the drive range in which the degree of magnetic saturation is small
  • the magnetic saturation method has a high estimation accuracy in the low-speed range and the drive range in which the degree of magnetic saturation is high and the induced voltage.
  • the interlinkage magnetic flux method is characterized by high estimation accuracy in the high speed range.
  • the relationship between the estimated position information and the driving area was not particularly mentioned. Therefore, in the second embodiment, the rotation speed of the rotor 1a is estimated in addition to the rotor position, and the estimated rotor position information is switched based on the estimated speed information and the magnetic saturation information. Will be described.
  • the magnetic saturation information is defined as information having a correlation with the degree of magnetic saturation of the rotary electric machine 1.
  • FIG. 7 is a diagram showing a configuration example of a control device 100A for a rotary electric machine according to the second embodiment.
  • the position estimator 4 is replaced with the position estimator 4A in the configuration of the control device 100 according to the first embodiment shown in FIG.
  • the other configurations are the same as or equivalent to those of the first embodiment, and the same or similar components are designated by the same reference numerals, and the duplicate description will be omitted.
  • FIG. 8 is a diagram showing a detailed configuration example of the position estimator 4A shown in FIG.
  • the position estimator 4A according to the second embodiment shown in FIG. 8 is the same as the position estimator 4 according to the first embodiment shown in FIG. 3, except that the estimated position calculators 406, 407, 408 and the estimated position switch 409 are respectively It is replaced by the estimated position/speed calculators 406A, 407A, 408A and the estimated position switch 409A.
  • the estimated position/velocity calculators 406A, 407A, and 408A are distinguished without a code, each of them is referred to as a “first estimated position/velocity calculator”, a “second estimated position/velocity calculator”, and a “third”. Of the estimated position/velocity calculator”.
  • the estimated position/velocity calculator 406A calculates the first estimated rotor position ⁇ ⁇ 1 and the first estimated speed ⁇ ⁇ 1 based on the three-phase AC components I uhac , I vhac , and I whac of the current amplitude for position estimation. And are calculated. Estimated position/velocity calculator 406A uses the salient pole method as in the first embodiment. That is, the first estimated speed ⁇ ⁇ 1 is an estimated speed calculated based on the saliency of the rotor 1a detected from the position estimation current.
  • Estimated position speed calculator 407A based on the DC component I hdc position estimation current amplitude, and calculates the second estimated rotor position theta ⁇ 2, and a second estimated speed omega ⁇ 2.
  • the estimated position/velocity calculator 407A uses the magnetic saturation method as in the first embodiment. That is, the second estimated speed ⁇ ⁇ 2 is an estimated speed calculated based on the DC component I hdc of the position estimation current amplitude.
  • Estimated position speed calculator 408A is the rotary electric machine drive current i uf on the three-phase coordinate, i vf, i wf and a flux linkage that is calculated by the induced voltage generated by the saliency of the rotor, the third on the basis of the The estimated rotor position ⁇ ⁇ 3 and the third estimated speed ⁇ ⁇ 3 are calculated.
  • the estimated position/velocity calculator 408A uses the induced voltage and the interlinkage magnetic flux method as in the first embodiment.
  • the third estimated speed ⁇ ⁇ 3 is the rotary electric machine drive currents i uf , i vf , and i wf on the three-phase coordinates, and the interlinkage magnetic flux calculated by the induced voltage generated by the saliency of the rotor 1a. It is an estimated speed calculated based on.
  • the estimated position switch 409A, ⁇ *, ⁇ f * , ⁇ ⁇ m, ⁇ ⁇ 1, i f *, i ff *, i f is input.
  • ⁇ * is the torque command value described above, and if is the drive current.
  • ⁇ f * is a value obtained by filtering the torque command value ⁇ * by taking into account the delay of the control system.
  • ⁇ ⁇ m is an estimated torque calculated from the drive current if and the rotary electric machine parameter based on the mathematical model of the rotary electric machine.
  • ⁇ ⁇ 1 is an estimated torque obtained by a look-up table using the drive current if as an argument. if * is a rotary electric machine drive current command.
  • Iff * is a value obtained through a filter that considers the delay of the control system in the rotating electric machine drive current command if * .
  • the frequency of the drive current i f is utilized.
  • the estimated position switch 409A at least one of these pieces of input information is used as magnetic saturation information.
  • FIG. 9 is a diagram showing a detailed configuration example of the estimated position/speed calculator 406A shown in FIG.
  • the estimated position/speed calculator 406A includes an estimated position calculator 406 and an estimated speed calculator 4060, as shown in FIG.
  • the estimated position calculator 406 shown in FIG. 9 is a component equivalent to the estimated position calculator 406 shown in FIG.
  • the estimated speed calculator 4060 calculates the first estimated speed ⁇ ⁇ 1 by pseudo-differentiating the first estimated rotor position ⁇ ⁇ 1 calculated by the estimated position calculator 406.
  • the estimated position/velocity calculator 406A outputs the calculated first estimated rotor position ⁇ ⁇ 1 and the first estimated speed ⁇ ⁇ 1 .
  • Pseudo-derivative is obtained by subjecting an input value to differential operation and filter processing.
  • the pseudo-differential processor can be realized by a differentiator and a low-pass filter.
  • FIG. 10 is a diagram showing a detailed configuration example of the estimated position/velocity calculator 407A shown in FIG.
  • the estimated position/velocity calculator 407A is a type 2 equipped with a lookup table (Look Up Table: LUT) 4070, a LUT 4071, a subtractor 4072, a divider 4073, a proportional integrator 4074, and an integrator 4075. It is configured as an arithmetic unit of.
  • LUT Look Up Table
  • the LUT 4070 stores the torque command value ⁇ * and I hdc * which is the DC component I hdc of the position estimation current amplitude obtained when driving with the torque command value ⁇ * .
  • the estimated position/velocity calculator 407A refers to the table value I hdc * stored in the LUT 4070.
  • the subtractor 4072 calculates a deviation ⁇ I h between the DC component I hdc of the position estimation current amplitude and the table value I hdc * .
  • the LUT 4071 stores a coefficient K Ih ⁇ 2 for converting the position estimation error ⁇ 2 from the torque command value ⁇ * .
  • the estimated position/velocity calculator 407A refers to the coefficient K Ih ⁇ 2 stored in the LUT 4071 .
  • the divider 4073 obtains the position estimation error ⁇ 2 by dividing the deviation ⁇ I h by K Ih ⁇ 2 .
  • the proportional integrator 4074 calculates the second estimated velocity ⁇ ⁇ 2 based on the position estimation error ⁇ 2 .
  • the integrator 4075 calculates the second estimated rotor position ⁇ ⁇ 2 based on the second estimated speed ⁇ ⁇ 2 .
  • FIG. 10 has a configuration in which the second estimated speed ⁇ ⁇ 2 is calculated first, and the second estimated rotor position ⁇ ⁇ 2 is calculated based on the second estimated speed ⁇ ⁇ 2. Not limited to. On the contrary, the second estimated rotor position ⁇ ⁇ 2 is calculated first, and the second estimated speed ⁇ ⁇ 2 is calculated by pseudo-differentiating the second estimated rotor position ⁇ ⁇ 2. May be
  • FIG. 10 shows an example in which the argument of the LUT 4070 and the LUT 4071 is the torque command value ⁇ * , but the argument is not limited to this.
  • the rotary electric machine drive current command if * may be used as an argument of the LUT 4070 and the LUT 4071.
  • FIG. 11 is a diagram showing a detailed configuration example of the estimated position/velocity calculator 408A shown in FIG.
  • the estimated position/velocity calculator 408A includes a position estimation error calculator 4080 and a PLL4081A, as shown in FIG.
  • the position estimation error calculator 4080 shown in FIG. 11 is a component equivalent to the position estimation error calculator 4080 shown in FIG.
  • the PLL4081A has an arithmetic unit configured such that the position estimation error “ ⁇ ( ⁇ 3- ⁇ )” becomes zero.
  • the PLL 4081A is configured as a type 2 computing unit including a proportional integrator 4081Aa and an integrator 4081Ab, like the estimated position/velocity computing unit 407A shown in FIG.
  • the proportional integrator 4081Aa calculates the third estimated velocity ⁇ ⁇ 3 based on the position estimation error "-( ⁇ ⁇ 3 - ⁇ )".
  • the integrator 4081Ab calculates the third estimated rotor position ⁇ ⁇ 3 based on the third estimated speed ⁇ ⁇ 3 .
  • the PLL 4081A in FIG. 11 has a configuration in which the third estimated speed ⁇ ⁇ 3 is calculated first, and then the third estimated rotor position ⁇ ⁇ 3 is calculated based on the third estimated speed ⁇ ⁇ 3. , But is not limited to this. On the contrary, the third estimated rotor position ⁇ ⁇ 3 is calculated first, and the third estimated speed ⁇ ⁇ 3 is calculated by pseudo differentiating the third estimated rotor position ⁇ ⁇ 3.
  • the PLL4081A may be configured in the above.
  • FIG. 12 is a diagram showing a detailed configuration example of the estimated position switch 409A shown in FIG.
  • the estimated position switch 409A includes an estimated speed selector 411A and an estimated position selector 410A.
  • the estimated speed selector 411A selects any one of the first estimated speed ⁇ ⁇ 1 , the second estimated speed ⁇ ⁇ 2, and the third estimated speed ⁇ ⁇ 3 , and the selected estimated speed ⁇ ⁇ Is output to the estimated position selector 410A as speed information.
  • the estimated position selector 410A selects or switches the estimated rotor position information based on the speed information output from the estimated speed selector 411A and the magnetic saturation information described above, and outputs the selected rotor position estimated information.
  • the speed information used for selecting or switching the estimated rotor position ⁇ ⁇ includes the first estimated speed ⁇ ⁇ 1 , the second estimated speed ⁇ ⁇ 2, and the third estimated speed ⁇ ⁇ 3 . Any one of them. That is, the estimated speed selector 411A selects one of the first estimated speed ⁇ ⁇ 1 , the second estimated speed ⁇ ⁇ 2, and the third estimated speed ⁇ ⁇ 3 to estimate the estimated position. Output to the selector 410A.
  • FIG. 13 is a diagram for explaining the drive region in the estimated position switch 409A shown in FIGS. 8 and 12.
  • the horizontal axis represents the rotation speed of the rotating electric machine
  • the vertical axis represents the torque of the rotating electric machine.
  • the rotation speed is an example of speed information.
  • the torque is information that correlates with the magnetic saturation information described above.
  • (1) is a drive area driven by a salient pole method
  • (2) is a drive area driven by a magnetic saturation method
  • (3) is a drive area driven by an induced voltage and interlinkage magnetic flux method.
  • the broken line in FIG. 13 means a boundary line with each of the above-described rotor position estimation methods, and serves as a boundary determination threshold value.
  • This boundary line is defined by magnetic saturation information and velocity information. That is, the magnetic saturation information and the velocity information are used to determine the region, and the rotor position estimation method defined in the region is selected.
  • the estimated position switch 409A selects and outputs the first estimated rotor position ⁇ ⁇ 1 when the rotation speed is smaller than the threshold and the torque is smaller than the threshold. Further, when the rotation speed is lower than the threshold value and the torque is higher than the threshold value, the estimated position switch 409A selects and outputs the second estimated rotor position ⁇ ⁇ 2 . Then, the estimated position switch 409A selects and outputs the third estimated rotor position ⁇ ⁇ 3 when the rotation speed is higher than the threshold value.
  • the boundary line for switching is shown by a straight line, but it does not have to be a straight line and may be a curved line.
  • an appropriate rotor position estimation method can be switched based on the rotor magnetic saturation information and the rotation speed. This enables desired position sensorless control in the torque speed range of the operating range.
  • the estimated position switch 409A selects or switches the estimated rotor position information based on the magnetic saturation information and the rotation speed of the rotor 1a.
  • the estimated information of the rotor position is selected or switched, the estimated information may change discontinuously before and after the switching, and a shock at the time of switching becomes a problem. Therefore, in the third embodiment, an embodiment in which the shock at the time of switching can be reduced will be described.
  • FIG. 14 is a diagram showing a configuration example of a control device 100B for a rotary electric machine according to the third embodiment.
  • the position estimator 4A is replaced with a position estimator 4B in the configuration of the control device 100A according to the second embodiment shown in FIG.
  • the other configurations are the same as or equivalent to those of the second embodiment, and the same or equivalent components are designated by the same reference numerals, and duplicate description will be omitted.
  • FIG. 15 is a diagram showing a detailed configuration example of the position estimator 4B shown in FIG.
  • the position estimator 4B in the third embodiment shown in FIG. 15 has the estimated position switch 409A replaced with the estimated position switch 409B in the configuration of the position estimator 4A in the second embodiment shown in FIG.
  • the other configurations are the same as or equivalent to those of FIG. 8, and the same or equivalent components are designated by the same reference numerals, and overlapping description will be omitted.
  • FIG. 16 is a diagram showing a detailed configuration example of the estimated position switch 409B shown in FIG.
  • the estimated position switch 409B in the third embodiment shown in FIG. 16 has the estimated position selector 410A replaced with the estimated position combiner 410B in the configuration of the estimated position switch 409A in the second embodiment shown in FIG. ..
  • the estimated position combiner 410B uses the magnetic saturation information and the speed information to estimate the first estimated rotor position ⁇ ⁇ 1 , the second estimated rotor position ⁇ ⁇ 2, and the third estimated position. and it outputs the combined position information using at least two estimated rotor position of the rotor position theta ⁇ 3 as the estimated rotor position theta ⁇ .
  • the magnetic saturation information used for combining the estimated rotor positions in the estimated position switch 409B is the same as the magnetic saturation information used in the second embodiment.
  • the speed information used to switch the estimated rotor position is the same as the speed information used in the second embodiment.
  • the estimated position combiner 410B based on the magnetic saturation information and the velocity information, estimates the first estimated rotor position ⁇ ⁇ 1 , the second estimated rotor position ⁇ ⁇ 2, and the third estimated rotor position ⁇ ⁇ 3. Are synthesized by a weighted average.
  • FIG. 17 is a diagram for explaining the drive region in the estimated position switch 409B shown in FIGS. 15 and 16.
  • the horizontal axis represents the rotation speed of the rotating electric machine
  • the vertical axis represents the torque of the rotating electric machine.
  • the torque is an example of information correlated with magnetic saturation information, as in the second embodiment.
  • (1) is a drive area that is driven only by the salient pole method
  • (2) is a drive area that is driven only by the magnetic saturation method
  • (3) is a drive area that is driven only by the induced voltage and the interlinkage magnetic flux method. is there.
  • the regions represented by a plurality of numbers are drive regions that use the estimation method corresponding to those numbers.
  • (1) and (2) are drive regions driven by the salient pole method and the magnetic saturation method.
  • the estimated rotor position to be output is determined by the first estimated rotor position ⁇ ⁇ 1 estimated by the salient pole method and the second estimated rotor position ⁇ ⁇ 2 estimated by the magnetic saturation method.
  • the position ⁇ ⁇ is synthesized.
  • the boundary of the drive region is shown by a broken line as in FIG. Further, in FIG. 17, the boundary of the drive region is shown by a straight line as in FIG. 13, but it is not necessary to be a straight line and may be a curved line.
  • W spd is a weight based on velocity information
  • W ms is a weight based on magnetic saturation information.
  • the values of the weight W spd and the weight W ms change in the range of 0 to 1 according to the velocity information and the magnetic saturation information, respectively.
  • the estimated position combiner 410B combines the estimated rotor position ⁇ ⁇ in the drive region defined in FIG. 17 using the following formulas (20) and (21).
  • Theta ⁇ 4 in the above formula (21) is a "weighted average by the weighting W ms of salient pole type magnetic saturation method.”
  • this weighted average value ⁇ ⁇ 4 is referred to as a “fourth estimated rotor position”.
  • ⁇ ⁇ 5 in the above equation (20) is defined as “the fourth estimated rotor position ⁇ ⁇ 4 and the third estimated rotor position ⁇ ⁇ 3 which is a calculated value only by the induced voltage and the interlinkage magnetic flux method. Weighted average value by the weight W spd of ". Hereinafter, a combination of the weighted averages will be referred to as a “fifth estimated rotor position”.
  • W spd 1 case
  • the estimated rotor position output from the estimated position synthesizer 410B theta ⁇ Becomes
  • the fourth estimated rotor position ⁇ ⁇ 4 which is a weighted average value according to the magnetic saturation method, becomes the estimated rotor position ⁇ ⁇ output from the estimated position combiner 410B.
  • the weighted average of at least two methods of the salient pole method, the magnetic saturation method, the induced voltage and the interlinkage magnetic flux method is used for the combination according to the value of the weight W ms.
  • the calculated fifth estimated rotor position ⁇ ⁇ 5 is output from the estimated position combiner 410B.
  • the order of the operations of the expressions (20) and (21) may be exchanged.
  • the fifth estimated rotor position may be calculated without calculating the fourth estimated rotor position by using an arithmetic expression obtained by substituting the equation (21) into the equation (20). Further, in the drive range defined in FIG.
  • the estimated rotor position ⁇ is synthesized by combining the defined first estimated rotor position, second estimated rotor position, and third estimated rotor position. If the calculation formula is, the calculation formula of the estimated rotor position ⁇ does not have to be the formulas (20) and (21).
  • the position estimation voltage calculator 5b in the controller 5 in FIG. 14 and the estimated position/speed calculators 406A, 407A, 408A in the position estimator 4B in FIG. 15 do not need to be constantly operated and are defined in FIG. It may be stopped depending on the drive area. For example, it is possible to operate only the components of the estimation method corresponding to the numbers in FIG. 17 and stop the operations of the components of the estimation method corresponding to the numbers not described. Further, in the area (3) in FIG. 17, the operation of the position estimation voltage calculation unit 5b can be stopped and its output can be set to zero. As a result, it is possible to reduce the power consumption associated with the arithmetic processing.
  • the estimated position/velocity calculator 406A, 407A, 408A starts the operation from the state in which the operation is stopped, as the initial value of the position/speed estimator for starting the operation, the other It goes without saying that the estimated rotor position of the position/speed estimator, the estimated speed, the output ⁇ ⁇ of the integrator output or estimated position combiner 410B, and the output ⁇ ⁇ of the estimated speed selector 411A can be given.
  • the estimated information of the three rotor positions obtained in the second embodiment is combined and output by the weighted average.
  • the second embodiment it is possible to reduce shock when switching estimated rotor position information.
  • the torque command value ⁇ * for the rotary electric machine 1 and the filter output value ⁇ f * of the torque command value ⁇ * are used as the magnetic saturation information used when switching the rotor position estimation method.
  • Estimated torque ⁇ ⁇ m obtained by the mathematical model of the rotating electric machine, estimated torque ⁇ ⁇ 1 obtained by a lookup table using the drive current if as arguments, rotating electric machine drive current command if * for the rotating electric machine 1, rotating electric machine drive current command i f * of the filter output value i ff *, and using at least one of the frequency of the drive current i f.
  • the magnetic saturation information used in the second and third embodiments cannot accurately grasp the degree of magnetic saturation. Is assumed. Therefore, in the fourth embodiment, the DC component Ihdc of the position estimation current amplitude described above is used as the magnetic saturation information when switching the rotor position estimation method.
  • the DC component Ihdc of the position estimation current amplitude is a component of the position estimation current amplitudes Iuh , Ivh , and Iwh that does not change depending on the rotor position. That is, the fourth embodiment is an embodiment in which the DC component Ihdc of the position estimation current amplitude that does not change depending on the rotor position is used as the magnetic saturation information.
  • FIG. 18 is a diagram showing a configuration example of a control device 100C for a rotary electric machine according to the fourth embodiment.
  • a control device 100C according to the fourth embodiment shown in FIG. 18 has the position estimator 4A replaced with a position estimator 4C in the configuration of the control device 100A according to the second embodiment shown in FIG.
  • the other configurations are the same as or equivalent to those of the second embodiment, and the same or equivalent components are designated by the same reference numerals, and duplicate description will be omitted.
  • FIG. 19 is a diagram showing a detailed configuration example of the position estimator 4C shown in FIG.
  • the position estimator 4C in the fourth embodiment shown in FIG. 18 has the estimated position switch 409A replaced with the estimated position switch 409C in the configuration of the position estimator 4A in the second embodiment shown in FIG.
  • a configuration in which the DC component I hdc position estimation current amplitudes extracted by the DC component extractor 405 is input to the estimated position switch 409C as the magnetic saturation information is input to the estimated position switch 409C as the magnetic saturation information.
  • the other configurations are the same as or equivalent to those of FIG. 8, and the same or equivalent components are designated by the same reference numerals, and overlapping description will be omitted.
  • FIG. 20 is a diagram showing a detailed configuration example of the estimated position switch 409C shown in FIG.
  • Estimated position selector 409C according to the fourth embodiment shown in FIG. 20 has estimated position selector 410A replaced with estimated position combiner 410C in the configuration of estimated position selector 409A according to the second embodiment shown in FIG. ..
  • the DC component Ihdc of the position estimation current amplitude is input to the estimated position combiner 410C as the magnetic saturation information.
  • the estimated position combiner 410C based on the magnetic saturation information and the velocity information, estimates the first estimated rotor position ⁇ ⁇ 1 , the second estimated rotor position ⁇ ⁇ 2, and the third estimated rotor position ⁇ ⁇ 3. Are synthesized by a weighted average.
  • FIG. 21 is a diagram for explaining the drive region in the estimated position switch 409C shown in FIGS. 19 and 20.
  • the definitions of the horizontal axis, the vertical axis, and the drive area are the same as those in the second embodiment.
  • It is the rotational speed of the rotating electric machine
  • the vertical axis is the torque of the rotating electric machine.
  • the torque on the vertical axis is an example of information correlated with the DC component I hdc of the position estimation current amplitude, as in the second embodiment.
  • W ms (I hdc ) that divides the vertical axis is a weight defined by the DC component I hdc of the position estimation current amplitude and velocity information, that is, magnetic saturation information and velocity information.
  • the boundary of the drive region is shown by a straight line as in FIGS. 13 and 17, but the boundary does not have to be a straight line and may be a curved line.
  • the value of the weight W spd changes in the range of 0 to 1 according to the speed information. Further, the value of the weight W ms (I hdc ) changes in the range of 0 to 1 according to the magnetic saturation information and the velocity information.
  • the estimated position combiner 410C combines the estimated rotor position ⁇ ⁇ in the drive region defined in FIG. 21, using the following expressions (22) and (23).
  • ⁇ ⁇ 4 in the above equation (23) is a “fourth estimated rotor” which is a “weighted average value of salient pole system and magnetic saturation system”, except that ⁇ ms 4 is weighted by W ms (I hdc ). It means "position”.
  • the above formula (22) is the same formula as the above formula (20), and the fourth estimated rotor position ⁇ ⁇ 4 weighted by the weight W ms (I hdc ), the induced voltage and the interlinkage magnetic flux are used. It is a fifth estimated rotor position ⁇ ⁇ 5 which is a weighted average value by the weight W spd with the third estimated rotor position ⁇ ⁇ 3 which is a calculated value only by the method. It should be noted that the order of operations of the expressions (22) and (23) may be exchanged. Further, the fifth estimated rotor position may be calculated without calculating the fourth estimated rotor position by using an arithmetic expression obtained by substituting the equation (23) into the equation (22).
  • the estimated rotor position ⁇ is synthesized by combining the defined first estimated rotor position, second estimated rotor position, and third estimated rotor position. If the calculation formula is, the calculation formula of the estimated rotor position ⁇ does not need to be the formulas (22) and (23).
  • the boundary of the drive region and the weight W ms (I hdc ) are defined based on the DC component I hdc of the current amplitude for position estimation, but the present invention is not limited to this.
  • the magnetic saturation information used in the second and third embodiments may be jointly specified.
  • the drive region in which a plurality of estimation methods are used together is defined by the weighted average, but the present invention is not limited to this.
  • the first estimated rotor position ⁇ ⁇ 1 , the second estimated rotor position ⁇ ⁇ 2 and the third estimated rotor position ⁇ ⁇ 3 are calculated .
  • the drive area may be defined so that any one of them is selected.
  • the DC component Ihdc of the current amplitude for position estimation can also be used for the abnormality determination of the computing unit that performs the processing of each estimation method and the availability determination when switching the estimation method.
  • the estimated position/velocity calculator 407A stores the relationship between the drive current vector and the DC component I hdc of the current amplitude for position estimation in advance as shown in FIG. 4 on the estimation principle. That is, since the estimated position/velocity calculator 407A grasps the DC component I hdc of the position estimation current amplitude obtained according to the driving condition, it stores the stored DC component I hdc and the actual DC component I hdc .
  • An abnormality can be determined by comparing Further, since it is possible to grasp whether or not the drive region can be driven by the magnetic saturation method from the relationship between the drive current vector and the direct current component I hdc , it is possible to switch the estimation method by using the actual direct current component I hdc. It becomes possible to determine whether or not.
  • the position estimation voltage calculator 5b in the controller 5 in FIG. 18 and the estimated position/speed calculators 406A, 407A, 408A in the position estimator 4C in FIG. 19 do not need to be constantly operated and are defined in FIG. It may be stopped depending on the drive area. For example, it is possible to operate only the components of the estimation method corresponding to the numbers in FIG. 21 and stop the operations of the components of the estimation method corresponding to the numbers not described. Further, in the region (3) in FIG. 21, the operation of the position estimation voltage calculation unit 5b can be stopped and its output can be set to zero. As a result, it is possible to reduce the power consumption associated with the arithmetic processing.
  • the estimated position/velocity calculator 406A, 407A, 408A starts the operation from the state in which the operation is stopped, as the initial value of the position/speed estimator for starting the operation, the other It goes without saying that the estimated rotor position of the position/speed estimator, the estimated speed, the output ⁇ ⁇ of the integrator output or estimated position combiner 410C, and the output ⁇ ⁇ of the estimated speed selector 411A can be given.
  • the degree of magnetic saturation can be accurately grasped regardless of the torque control error, the current control error, and the position estimation error. Therefore, in addition to the effect of the third embodiment, it is possible to improve the estimation accuracy of the rotor position.
  • Embodiment 5 the first estimated speed ⁇ ⁇ 1 , the second estimated speed ⁇ ⁇ 2, and the third estimated speed ⁇ are used as the speed information used to switch the estimated information of the rotor position. It was an embodiment using any one of ⁇ 3 .
  • the estimation accuracy of the first estimated speed ⁇ ⁇ 1 , the second estimated speed ⁇ ⁇ 2, and the third estimated speed ⁇ ⁇ 3 differs depending on the driving condition.
  • the first estimated speed ⁇ ⁇ 1 has high estimation accuracy in a low speed region and a drive region in which the degree of magnetic saturation is small
  • the second estimated speed ⁇ ⁇ 2 is a drive region in a low speed region and a high degree of magnetic saturation.
  • the estimation accuracy is high, and the third estimation speed ⁇ ⁇ 3 has a high estimation accuracy in the high speed range.
  • the relationship between the estimated speed information and the drive area is not particularly mentioned. Therefore, in the fifth embodiment, an embodiment will be described in which the accuracy of the speed information used to switch the estimated rotor position information is increased.
  • FIG. 22 is a diagram showing a configuration example of a control device 100D for a rotary electric machine according to the fifth embodiment.
  • the control device 100D according to the fifth embodiment shown in FIG. 22 has the position estimator 4C replaced with a position estimator 4D in the configuration of the control device 100C according to the fourth embodiment shown in FIG.
  • Other configurations are the same as or equivalent to those of the fourth embodiment, and the same or equivalent components are designated by the same reference numerals, and duplicated description will be omitted.
  • FIG. 23 is a diagram showing a detailed configuration example of the position estimator 4D shown in FIG.
  • the position estimator 4D according to the fifth embodiment shown in FIG. 23 has the configuration of the position estimator 4C according to the fourth embodiment shown in FIG. 19, in which the estimated position switch 409C is replaced with the estimated position switch 409D.
  • the other configurations are the same as or equivalent to those of FIG. 8, and the same or equivalent components are designated by the same reference numerals, and overlapping description will be omitted.
  • FIG. 24 is a diagram showing a detailed configuration example of the estimated position switch 409D shown in FIG.
  • the estimated position switch 409D according to the fifth embodiment shown in FIG. 24 is the same as the estimated position switch 409C according to the fourth embodiment shown in FIG. 20 except that the estimated position combiner 410C is replaced with the estimated position combiner 410D.
  • the speed selector 411A is replaced with the estimated speed synthesizer 411D.
  • the DC component I hdc position estimation current amplitude is input to the estimated position synthesizer 410D as the magnetic saturation information, the output of the estimated velocity combiner 411D is fed back to itself.
  • the estimated speed synthesizer 411D selects one of the first estimated speed ⁇ ⁇ 1 , the second estimated speed ⁇ ⁇ 2, and the third estimated speed ⁇ ⁇ 3 based on the magnetic saturation information and the speed information. Using the speed information obtained by selecting, or at least two estimated speeds, the speed information combined by the weighted average is output as the estimated speed ⁇ ⁇ . As the speed information for selecting or combining the estimated speed ⁇ ⁇ , the output of the estimated speed combiner 411D, that is, the output of itself is used. Instead of the estimated speed ⁇ ⁇ , any one of the first estimated speed ⁇ ⁇ 1 , the second estimated speed ⁇ ⁇ 2, and the third estimated speed ⁇ ⁇ 3 may be used.
  • the magnetic saturation information in the fifth embodiment is the DC component Ihdc of the position estimation current amplitudes I uh , I vh , and I wh that does not change depending on the rotor position.
  • the magnetic saturation information described in the second and third embodiments may be used. Therefore, the estimated position combiner 410D shown in FIG. 24 is replaced with any of the estimated position selector 410A shown in FIG. 12, the estimated position combiner 410B shown in FIG. 16, and the estimated position combiner 410C shown in FIG. can do.
  • the numerical values shown in FIG. 21 can be regarded as the estimated speed calculated in each region.
  • (1) is an area for calculating the first estimated speed ⁇ ⁇ 1
  • (2) is an area for calculating the second estimated speed ⁇ ⁇ 2
  • (3) is a third estimated speed ⁇ 2. It is an area for calculating ⁇ 3 .
  • the regions represented by a plurality of numbers as in (1) and (2) are regions in which a plurality of estimated velocities corresponding to those numbers are combined based on the magnetic saturation information and the velocity information described above.
  • the weighted average of the first estimated speed ⁇ ⁇ 1 and the second estimated speed ⁇ ⁇ 2 is calculated.
  • the broken line in the figure is a boundary when switching which estimated speed is used, and this boundary is defined by the magnetic saturation information and the speed information described above. Although this boundary is shown as a straight line, it does not have to be a straight line and may be a curved line.
  • the value of the weight W spd changes in the range of 0 to 1 according to the speed information. Further, the value of the weight W ms (I hdc ) changes in the range of 0 to 1 according to the magnetic saturation information and the velocity information.
  • the estimated speed combiner 411D combines the estimated speed ⁇ ⁇ in the region defined in FIG. 21 using the following expressions (24) and (25).
  • ⁇ ⁇ 4 in the above equation (25) is a “weighted average value by weights W ms (I hdc ) of the salient pole method and the magnetic saturation method”.
  • this weighted average value ⁇ ⁇ 4 is referred to as a “fourth estimated speed”.
  • ⁇ ⁇ 5 in the above equation (24) is the weight W spd of “the fourth estimated speed ⁇ ⁇ 4 and the third estimated speed ⁇ ⁇ 3 which is a calculated value only by the induced voltage and the flux linkage method. Is the weighted average.
  • a combination of the weighted tree average values will be referred to as a “fifth estimated speed”.
  • the estimated speed ⁇ ⁇ synthesized by the weighted average is output from the estimated speed synthesizer 411D.
  • the area in which the estimated speed is combined by the weighted average is defined, but the present invention is not limited to this. Any one of the first estimated speed ⁇ ⁇ 1 , the second estimated speed ⁇ ⁇ 2, and the third estimated speed ⁇ ⁇ 3 is based on the DC component I hdc of the position estimation current amplitude and the velocity information.
  • the area may be defined to be selected. It should be noted that the order of the operations of the expressions (24) and (25) may be exchanged.
  • the fifth estimated speed may be calculated without calculating the fourth estimated speed by using an arithmetic expression obtained by substituting the expression (25) into the expression (24). Further, in the driving area defined in FIG. 21, if the calculation formula is such that the estimated speed ⁇ is combined by the combination of the defined first estimated speed, second estimated speed, and third estimated speed, The equation for calculating the speed ⁇ need not be equations (24) and (25).
  • the fifth embodiment it is possible to prevent a sharp change in the estimated speed and to accurately grasp the speed information when switching the estimated position.
  • the degree of magnetic saturation can be accurately grasped regardless of the torque control error, the current control error, and the position estimation error. Therefore, the estimation accuracy of the rotor position can be improved as compared with the fourth embodiment.
  • FIG. 25 is a diagram showing a first hardware configuration example for realizing each function or part of the control device for a rotating electric machine according to the first to fifth embodiments.
  • FIG. 26 is a diagram showing a second hardware configuration example for realizing each function or part of the control device for a rotating electric machine according to the first to fifth embodiments.
  • FIG. 25 shows an example in which the above processing circuit is realized by dedicated hardware such as the dedicated processing circuit 1000.
  • FIG. 26 shows an example in which the above processing circuit is realized by the processor 1001 and the storage device 1002.
  • the functions of the current detector 2 and the voltage applicator 3 in FIGS. 1, 7, 14, 18, and 22 are realized using dedicated hardware, and the position estimator is used.
  • 4, 4A, 4B, 4C, 4D and the controller 5 are realized by the dedicated processing circuit 1000.
  • the dedicated processing circuit 1000 may be a single circuit, a composite circuit, a programmed processor, a parallel programmed processor, an ASIC, an FPGA, or a combination thereof. Applicable
  • Each of the above functions may be realized by a processing circuit or may be collectively realized by a processing circuit.
  • the functions of the current detector 2 and the voltage applicator 3 in FIGS. 1, 7, 14, 18, and 22 are realized using dedicated hardware, and the position estimator is used.
  • 4, 4A, 4B, 4C, 4D and the controller 5 are realized by the processor 1001 that executes the program recorded in the storage device 1002.
  • the plurality of processors 1001 and the plurality of storage devices 1002 may cooperate with each other to realize the above-described functions.
  • each function 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 1002.
  • the processor 1001 reads and executes the program stored in the storage device 1002. It can also be said that these programs cause a computer to execute procedures and methods for executing each function.
  • the memory device 1002 includes a RAM (Random Access Memory), a ROM (Read Only Memory), a flash memory, an EPROM (Erasable Programmable Read Only Memory), or an EEPROM (registered trademark) (Electrical Anywhere Memory). To do.
  • the semiconductor memory may be a non-volatile memory or a volatile memory.
  • the storage device 1002 corresponds to a magnetic disk, a flexible disk, an optical disk, a compact disk, a mini disk, or a DVD (Digital Versatile Disc).
  • the rotary electric machine 1 is described as a synchronous reluctance motor, but the invention is not limited to this.
  • the control devices 100, 100A, 100B, 100C, 100D according to the first to fifth embodiments are also applicable to a motor having a salient polarity such as an embedded magnet type synchronous motor or a surface magnet type synchronous motor.
  • the rotor position is estimated from the direct current component of the position estimation current amplitude.
  • the position estimation voltage command is used to improve the position estimation accuracy or noise.
  • the size of may be variable.
  • the ratio of the DC component and the position estimation voltage command for example, the DC component of the inductance can be used.
  • the controller 5 of the control devices 100, 100A, 100B, 100C, 100D has been described as controlling the torque, but the present invention is not limited to this.
  • the controller 5 can also be configured to control the rotation speed.
  • the current command of the rotary electric machine 1 with respect to the torque is selected so as to minimize the current effective value, that is, the copper loss, but the present invention is not limited to this.
  • the flux linkage or the loss of the rotary electric machine 1 may be set to be minimum.
  • the present invention is not limited to this.
  • the control devices 100, 100A, 100B, 100C, and 100D only need to be able to detect the phase current, and detect the phase current of the rotary electric machine 1 by a current sensor built in an inverter (not shown) that constitutes the voltage applicator 3. It may be configured.

Landscapes

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

Abstract

回転電機の制御装置(100)は、回転電機電流を検出する電流検出器(2)、回転電機電流に基づいて回転電機(1)の回転子位置を推定する位置推定器(4)、及び回転電機(1)を駆動するための回転電機駆動電圧指令と、回転子位置を推定するための位置推定用電圧指令とを演算する制御器(5)を備える。位置推定器(4)は、位置推定用電流振幅の交流成分に基づく第1の推定回転子位置と、位置推定用電流振幅の直流成分に基づく第2の推定回転子位置と、回転電機電流と鎖交磁束とに基づく第3の推定回転子位置とを演算し、第1の推定回転子位置、第2の推定回転子位置及び第3の推定回転子位置のうちの何れか1つを選択することで得られる推定回転子位置、又は、第1の推定回転子位置、第2の推定回転子位置及び第3の推定回転子位置のうちの少なくとも2つを用いて合成することで得られる推定回転子位置を出力する。

Description

回転電機の制御装置
 本発明は、回転子の回転位置を検出する位置センサを用いることなく回転子位置情報を得て回転電機を制御する、回転電機の制御装置に関する。
 回転電機の駆動には、回転子の位置情報が必要である。そのため、一般に、回転電機の制御装置には、回転子位置取得のための位置センサが用いられる。しかしながら、位置センサの使用により、システムの大型化、高コスト化、耐環境性の低下といった問題が生じる。このため、回転電機の制御装置には、位置センサを用いずに回転電機を駆動する位置センサレス制御の適用が求められる。
 位置センサレス制御は、回転子の突極性を利用して回転子位置を推定する方式と、回転電機に生じる誘起電圧によって演算される鎖交磁束を利用して回転子位置を推定する方式とに概ね大別される。以下、前者を「突極方式」と呼び、後者を「誘起電圧及び鎖交磁束方式」と呼ぶ。ここで、回転子の突極性とは、回転子のインダクタンスの磁気的異方性であり、インダクタンスが回転子位置によって変化する特性である。このため、突極方式は、インダクタンスの角度依存性を利用する方式であるとも言える。
 突極方式では、回転電機に位置推定用電圧もしくは位置推定電流を重畳することにより突極性に関する情報を励起させ、励起された情報に基づいて回転子位置が推定される。一般に、位置推定に必要な誘起電圧が十分に得られない低速領域では突極方式が使用され、誘起電圧が十分に得られる高速領域では誘起電圧及び鎖交磁束方式が使用される。このように、従来の位置センサレス制御では、トルク速度領域に応じて2つの位置センサレス制御方式を切り替えることで、トルク速度領域での位置センサレス制御を実現している。
 上述した技術的な背景の下、近年、回転電機の高出力密度化の要求が高まり、磁気飽和領域を積極的に利用する磁気的設計がなされた回転電機が出現している。以下、このような磁気的設計がなされた回転電機を「高出力密度回転電機」と呼ぶ。また、高出力密度回転電機の出現によって、トルク速度領域は拡大される。以下、高出力密度回転電機の出現後のトルク速度領域を「拡大されたトルク速度領域」と呼ぶ。
 高出力密度回転電機に対して、突極方式と、誘起電圧及び鎖交磁束方式とを適用した場合、拡大されたトルク速度領域の全ての範囲を網羅的に駆動することは困難である。特に、低速領域、且つ、磁気飽和の程度が大きい低速高トルク領域においては、推定位置に磁気飽和の程度に応じた補正処理を行ったとしても、十分な位置推定精度が得られないか、或いは、位置推定誤差が増大して、回転電機の動作が不安定化する現象が発生する。
 この課題に対して、下記特許文献1には、磁気飽和の異方性を用いて回転子位置を推定する方式が開示されている。この方式は「磁気飽和方式」と呼ばれる。磁気飽和方式では、回転座標のdq軸のそれぞれに位置推定用電圧が重畳され、位置推定用電圧の重畳によって発生する位置推定用電流のd軸の振幅と、q軸の振幅との乗算値に基づいて位置推定が行われる。
 特許文献1は、突極方式の位置推定器と、磁気飽和方式の位置推定器とを有している。特許文献1では、突極方式によって位置推定誤差相関量Δθ1が演算され、磁気飽和方式によって位置推定誤差相関量Δθ2が演算される。そして、低速低トルク領域では、位置推定誤差相関量Δθ1を用いて回転子位置が推定され、低速高トルク領域では位置推定誤差相関量Δθ2を用いて回転子位置が推定される。そして、位置推定誤差相関量Δθ1及び位置推定誤差相関量Δθ2の寄与度を磁気飽和の程度に応じて加重平均し、加重平均した寄与度に基づいて、回転位置を推定している。
特許第5145850号公報
 しかしながら、特許文献1は、高速領域駆動用の位置推定器を備えていない。このため、特許文献1の技術では、高速領域における位置センサレス制御は困難である。
 また、特許文献1の技術である磁気飽和方式を、上述した突極方式、並びに誘起電圧及び鎖交磁束方式と組み合わせることも考えられる。しかしながら、特許文献1を含む従来技術では、突極方式、誘起電圧及び鎖交磁束方式、並びに磁気飽和方式という3つの推定方式をどのように切り替えてその推定値を出力するのか、或いはどのように組み合わせてそれらの合成値を出力するのかの十分な検討は為されていない。このため、特許文献1を含む従来技術では、高出力密度回転電機の制御において、拡大されたトルク速度領域を網羅した位置センサレス制御を実施することは困難である。
 本発明は、上記に鑑みてなされたものであって、拡大されたトルク速度領域を網羅した位置センサレス制御を実施できる回転電機の制御装置を得ることを目的とする。
 上述した課題を解決し、目的を達成するために、本発明は、多相の回転電機の駆動制御を行う回転電機の制御装置であって、回転電機に流れる回転電機電流を検出する電流検出手段、回転電機電流に基づいて回転電機の回転子の回転位置である回転子位置を推定する位置推定手段、及び回転電機電流と、回転子位置の推定値である推定回転子位置とに基づいて回転電機を駆動するための回転電機駆動電圧指令と、回転子位置を推定するための各相の位置推定用電圧指令とを演算する制御器を備える。位置推定手段は、位置推定用電圧指令に基づく位置推定用電圧の印加により生じる位置推定用電流の振幅である位置推定用電流振幅の交流成分に基づく第1の推定回転子位置と、位置推定用電流振幅の直流成分に基づく第2の推定回転子位置と、回転電機電流と、回転子の突極性により生じる誘起電圧により演算される鎖交磁束とに基づく第3の推定回転子位置とを演算する。位置推定手段は、第1の推定回転子位置、第2の推定回転子位置及び第3の推定回転子位置のうちの何れか1つを選択することで得られる推定回転子位置、又は、第1の推定回転子位置、第2の推定回転子位置及び第3の推定回転子位置のうちの少なくとも2つを用いて合成することで得られる推定回転子位置を出力する。
 本発明に係る回転電機の制御装置によれば、拡大されたトルク速度領域を網羅した位置センサレス制御を実施できるという効果を奏する。
実施の形態1に係る回転電機の制御装置の構成例を示す図 図1に示す位置推定用電圧演算部から出力される位置推定用電圧指令を示す図 図1に示す位置推定器の細部の構成例を示す図 図3に示す第2の推定位置演算器の動作説明に供する第1の図 図3に示す第2の推定位置演算器の動作説明に供する第2の図 図3に示す第3の推定位置演算器の細部の構成例を示す図 実施の形態2に係る回転電機の制御装置の構成例を示す図 図7に示す位置推定器の細部の構成例を示す図 図8に示す第1の推定位置速度演算器の細部の構成例を示す図 図8に示す第2の推定位置速度演算器の細部の構成例を示す図 図8に示す第3の推定位置速度演算器の細部の構成例を示す図 図8に示す推定位置切替器の細部の構成例を示す図 図8及び図12に示す推定位置切替器における駆動領域の説明に供する図 実施の形態3に係る回転電機の制御装置の構成例を示す図 図14に示す位置推定器の細部の構成例を示す図 図15に示す推定位置切替器の細部の構成例を示す図 図15及び図16に示す推定位置切替器における駆動領域の説明に供する図 実施の形態4に係る回転電機の制御装置の構成例を示す図 図18に示す位置推定器の細部の構成例を示す図 図19に示す推定位置切替器の細部の構成例を示す図 図19及び図20に示す推定位置切替器における駆動領域の説明に供する図 実施の形態5に係る回転電機の制御装置の構成例を示す図 図22に示す位置推定器の細部の構成例を示す図 図23に示す推定位置切替器の細部の構成例を示す図 実施の形態1から5に係る回転電機の制御装置の第1のハードウェア構成例を示す図 実施の形態1から5に係る回転電機の制御装置の第2のハードウェア構成例を示す図
 以下に、本発明の実施の形態に係る回転電機の制御装置を図面に基づいて詳細に説明する。なお、以下の実施の形態により、本発明が限定されるものではない。また、以下では「回転電機の制御装置」を単に「制御装置」と称する場合がある。
実施の形態1.
 実施の形態1に係る回転電機の制御装置は、多相の回転電機の駆動制御を行う制御装置である。多相の回転電機とは、三相もしくは四相以上の交流電圧が印加可能に構成されている回転電機である。
 図1は、実施の形態1に係る回転電機の制御装置100の構成例を示す図である。実施の形態1に係る制御装置100は、図1に示すように、電圧印加手段である電圧印加器3と、電流検出手段である電流検出器2と、位置推定手段である位置推定器4と、制御手段である制御器5とを備える。
 電圧印加器3は、回転電機1を駆動するための回転電機電圧指令v ,v ,v に基づいて、回転電機1に交流電力を供給する。
 回転電機1の一例は、回転子1aの突極性を利用してトルクを発生する三相の同期リラクタンスモータである。同期リラクタンスモータは、回転子1aの磁気抵抗が回転子位置によって変化するモータである。
 電流検出器2は、回転電機1に流れる電流である回転電機電流i,i,iを検出する。回転電機電流i,i,iは、電圧印加器3から回転電機1の各相に供給される交流電流である。電流検出器2は、検出した回転電機電流i,i,iを位置推定器4及び制御器5に出力する。
 位置推定器4は、回転電機電流i,i,iに基づいて推定回転子位置θを演算する。制御器5は、回転電機電流i,i,i及び推定回転子位置θに基づいて、回転電機1の出力トルクがトルク指令値τで指示される値となるように回転電機電圧指令v ,v ,v を演算する。推定回転子位置θは、回転電機1を構成する回転子1aの回転位置である回転子位置の推定値である。推定回転子位置θは、電気角で表される。
 制御器5は、第1の演算部である駆動電圧指令演算部5aと、第2の演算部である位置推定用電圧演算部5bとに区分される。駆動電圧指令演算部5aは、電流指令演算器6、電流制御器7、回転座標逆変換器8、二相三相変換器9、駆動電流抽出器11、三相二相変換器12、回転座標変換器13及び加算器14を備える。なお、図1に示す駆動電圧指令演算部5a及び位置推定用電圧演算部5bの区分は一例であり、制御器5の構成要素をどのように区分してもよい。
 位置推定用電圧演算部5bは、回転子位置を推定するための各相の位置推定用電圧指令vuf ,vvf ,vwf を演算する。駆動電圧指令演算部5aは、回転電機電流i,i,i及び推定回転子位置θに基づいて回転電機1を駆動するための回転電機駆動電圧指令vuf ,vvf ,vwf を演算する。なお、回転電機駆動電圧指令vuf ,vvf ,vwf は、加算器14で位置推定用電圧指令vuh ,vvh ,vwh が加算される。加算器14の出力は、回転電機電圧指令v ,v ,v として、電圧印加器3に出力される。
 駆動電圧指令演算部5aにおいて、電流指令演算器6には、回転電機1の出力トルクの指令値であるトルク指令値τが入力される。
 電流指令演算器6は、トルク指令値τを用いて、回転二相座標上の回転電機駆動電流指令idf ,iqf を演算する。回転電機駆動電流指令idf ,iqf は、回転電機1がトルク指令値τに対応した出力を発生するために必要な回転二相座標上の電流指令である。実施の形態1の電流指令演算器6では、トルクに対する電流実効値が最小、別言すると、トルクに対する銅損が最小になるような電流指令が演算される。
 回転二相座標上の回転電機駆動電流指令idf ,iqf のうち、回転電機駆動電流指令idf は、回転子1aの磁気抵抗が最も小さくなるd軸方向の電機子電流成分を示すd軸駆動電流の指令値である。回転電機駆動電流指令iqf は、d軸に直交する方向となるq軸方向の電機子電流成分を示すq軸駆動電流の指令値である。回転二相座標上の回転電機駆動電流指令idf ,iqf の演算には、トルク指令値τ以外にも、回転電機1のモータ定数が用いられる。モータ定数としては、回転電機1の相互インダクタンスと、回転電機1の極数とが例示される。なお、モータ定数に代えて、予め求めた電流指令とトルクとの間の関係式又はテーブルを用いてもよい。
 駆動電流抽出器11は、電流検出器2で検出された三相座標上の回転電機電流i,i,iから、三相座標上の回転電機駆動電流iuf,ivf,iwfを抽出する。三相座標上の回転電機駆動電流iuf,ivf,iwfは、回転電機1を駆動するための三相座標上の回転電機駆動電圧指令vuf ,vvf ,vwf により発生した三相座標上の回転電機駆動電流である。回転電機駆動電圧指令vuf ,vvf ,vwf は、二相三相変換器9から出力されて加算器14に入力される駆動電圧指令である。加算器14には、回転電機駆動電圧指令vuf ,vvf ,vwf と、回転電機1の回転子位置を推定するための電圧指令である位置推定用電圧指令vuh ,vvh ,vwh とが入力される。
 加算器14は、三相座標上の回転電機駆動電圧指令vuf ,vvf ,vwf と、三相座標上の位置推定用電圧指令vuh ,vvh ,vwh とを加算した三相座標上の回転電機電圧指令v ,v ,v を生成して、電圧印加器3へ出力する。位置推定用電圧指令vuh ,vvh ,vwh は、位置推定用電圧演算部5bによって演算される。
 図2は、図1に示す位置推定用電圧演算部5bから出力される位置推定用電圧指令を示す図である。図2には、三相座標上の位置推定用電圧指令vuh ,vvh ,vwh の例として、互いに120°の位相差を有する方形波状の電圧が、上段部側からu相、v相及びw相の順で示されている。なお、図2では、三相座標上の位置推定用電圧指令vuh ,vvh ,vwh が方形波電圧である場合を示しているが、これに限定されない。方形波電圧に代えて、正弦波電圧を用いてもよい。
 図1に戻り、実施の形態1の駆動電流抽出器11では、例えばノッチフィルタにより三相座標上の回転電機電流i,i,iから、位置推定用電圧指令vuh ,vvh ,vwh の印加により発生する三相座標上の位置推定用電流振幅Iuh,Ivh,Iwhを除去することにより、三相座標上の回転電機駆動電流iuf,ivf,iwfを抽出する。なお、三相座標上の回転電機電流i,i,iから三相座標上の回転電機駆動電流iuf,ivf,iwfを抽出する方法は、ノッチフィルタに限定されず、ローパスフィルタ又はハイパスフィルタを用いてもよい。
 三相二相変換器12は、駆動電流抽出器11で抽出された三相座標上の回転電機駆動電流iuf,ivf,iwfを静止二相座標上の回転電機駆動電流iαf,iβfに変換する。回転座標変換器13は、位置推定器4によって推定された推定回転子位置θを用いて座標変換を行い、静止二相座標上の回転電機駆動電流iαf,iβfを、回転二相座標上の回転電機駆動電流idf,iqfに変換する。
 電流制御器7は、回転座標変換器13で変換された、回転二相座標上の回転電機駆動電流idf,iqfが、電流指令演算器6で演算された回転電機駆動電流指令idf ,iqf となるように電流制御を行い、回転二相座標上の回転電機駆動電圧指令vdf ,vqf を演算する。電流制御器7における電流制御としては、比例積分(Proportional Integral:PI)制御を例示できる。
 回転座標逆変換器8は、推定回転子位置θを用いて、電流制御器7で演算された回転二相座標上の回転電機駆動電圧指令vdf ,vqf を、静止二相座標上の回転電機駆動電圧指令vαf ,vβf に変換する。二相三相変換器9は、静止二相座標上の回転電機駆動電圧指令vαf ,vβf を、前述した三相座標上の回転電機駆動電圧指令vuf ,vvf ,vwf に変換する。
 次に、位置推定器4の細部について説明する。図3は、図1に示す位置推定器4の細部の構成例を示す図である。位置推定器4は、図3に示すように、信号処理器41、位置推定用電流振幅演算器403、信号処理器42、第1の推定位置演算器である推定位置演算器406、第2の推定位置演算器である推定位置演算器407、第3の推定位置演算器である推定位置演算器408、及び推定位置切替器409を備える。
 信号処理器41は、位置推定用電流抽出器401及び駆動電流抽出器402を備える。位置推定用電流抽出器401は、三相座標上の回転電機電流i,i,iから位置推定用電流iuh,ivh,iwhを抽出する。駆動電流抽出器402は、三相座標上の回転電機電流i,i,iから三相座標上の回転電機駆動電流iuf,ivf,iwfを抽出する。
 位置推定用電流振幅演算器403は、三相座標上の位置推定用電流iuh,ivh,iwhに基づいて、三相座標上の位置推定用電流振幅Iuh,Ivh,Iwhを演算する。
 信号処理器42は、交流成分抽出器404及び直流成分抽出器405を備える。交流成分抽出器404は、位置推定用電流振幅Iuh,Ivh,Iwhから位置推定用電流振幅の三相交流成分Iuhac,Ivhac,Iwhacを抽出する。直流成分抽出器405は、位置推定用電流振幅Iuh,Ivh,Iwhから位置推定用電流振幅の直流成分Ihdcを抽出する。
 推定位置演算器406は、位置推定用電流振幅の三相交流成分Iuhac,Ivhac,Iwhacに基づいて、第1の推定回転子位置θ を推定する。推定位置演算器406において用いられる回転子位置の推定方式は、突極方式である。即ち、第1の推定回転子位置θ は、位置推定用電圧の印加により生じる位置推定用電流の振幅である位置推定用電流振幅の交流成分に基づいて演算される推定回転子位置である。第1の推定回転子位置θ は、位置推定器4から最終的に出力される推定回転子位置θの候補値の1つである。
 推定位置演算器407は、位置推定用電流振幅の直流成分Ihdcに基づいて、第2の推定回転子位置θ を演算する。推定位置演算器407において用いられる回転子位置の推定方式は、磁気飽和方式である。第2の推定回転子位置θ は、位置推定器4から最終的に出力される推定回転子位置θの候補値の1つである。
 推定位置演算器408は、三相座標上の回転電機駆動電流iuf,ivf,iwfと、回転子1aの突極性により生じる誘起電圧により演算される鎖交磁束と、に基づいて第3の推定回転子位置θ を演算する。推定位置演算器408において用いられる回転子位置の推定方式は、誘起電圧及び鎖交磁束方式である。第3の推定回転子位置θ は、位置推定器4から最終的に出力される推定回転子位置θの候補値の1つである。
 推定位置切替器409は、第1の推定回転子位置θ 、第2の推定回転子位置θ 及び第3の推定回転子位置θ のうちの何れか1つを選択することで得られる位置情報、又は、少なくとも2つの推定回転子位置を用いて合成することで得られる位置情報を推定回転子位置θとして出力する。
 次に、位置推定器4の更に詳細な機能について、図3に加え、更に図4から図6の図面を参照して説明する。図4は、図3に示す推定位置演算器407の動作説明に供する第1の図である。図5は、図3に示す推定位置演算器407の動作説明に供する第2の図である。図6は、図3に示す推定位置演算器408の細部の構成例を示す図である。
 駆動電流抽出器402は、例えばノッチフィルタを用いて、三相座標上の回転電機電流i,i,iから、位置推定用電圧指令vuh ,vvh ,vwh の印加により発生する三相座標上の位置推定用電流振幅Iuh,Ivh,Iwhを除去することにより、三相座標上の回転電機駆動電流iuf,ivf,iwfを抽出する。位置推定用電流抽出器401は、例えば三相座標上の回転電機電流i,i,iから、駆動電流抽出器402で演算された三相座標上の回転電機駆動電流iuf,ivf,iwfを減算することにより、三相座標上の位置推定用電流振幅Iuh,Ivh,Iwhを演算する。なお、三相座標上の回転電機駆動電流iuf,ivf,iwf及び三相座標上の位置推定用電流振幅Iuh,Ivh,Iwhの抽出方法又は演算方法は、これに限定されず、バンドパスフィルタ、バンドストップフィルタ、ローパスフィルタ又はハイパスフィルタを用いてもよい。
 前述の通り、位置推定用電流振幅演算器403は、三相座標上の位置推定用電流振幅Iuh,Ivh,Iwhに基づいて、三相座標上の位置推定用電流振幅Iuh,Ivh,Iwhを演算する。なお、三相座標上の位置推定用電流振幅Iuh,Ivh,Iwhは、前述した位置推定用電流振幅の直流成分Ihdcと、位置推定用電流振幅の交流成分Ihacとを用いて、下記(1)式のように表すことができる。
Figure JPOXMLDOC01-appb-M000001
 なお、位置推定用電流振幅演算器403による演算方法は公知であり、ここでの詳細な説明は割愛する。演算方法の詳細は、例えば特許第5324646号公報の明細書の段落[0034]から[0055]に記載されており、当該記載内容を参照されたい。当該記載内容は、本明細書に取り込まれて本明細書の一部を構成する。
 また、本手法は、三相座標上の位置推定用電流振幅Iuh,Ivh,Iwhの交流成分の相対関係を利用して磁極位置を演算するものであり、位置推定用電流振幅Iuh,Ivh,Iwhの絶対値を演算する必要がない。このため、三相座標上の位置推定用電流振幅Iuh,Ivh,Iwhの演算には、下記(2)式を用いることができる。
Figure JPOXMLDOC01-appb-M000002
 上記(2)式において、“t”は時間であり、“T”は交流の1周期である。上記(2)式では、自己相関により絶対値を求めるときに積分記号の前に付される“√(2/T)”の係数が省略されている。即ち、位置推定用電流振幅演算器403による演算では、乗算及び平方根演算は不要である。従って、実施の形態1の位置推定用電流振幅演算器403によれば、演算処理の高速化及び演算時間の短縮化が可能である。
 前述の通り、直流成分抽出器405は、位置推定用電流振幅Iuh,Ivh,Iwhから位置推定用電流振幅の直流成分Ihdcを抽出する。具体的に、直流成分抽出器405は、下記(3)式に示されるように、位置推定用電流振幅演算器403で演算された三相座標上の位置推定用電流振幅Iuh,Ivh,Iwhの平均値を演算することで、位置推定用電流振幅の直流成分Ihdcを抽出する。
Figure JPOXMLDOC01-appb-M000003
 また、前述の通り、交流成分抽出器404は、位置推定用電流振幅Iuh,Ivh,Iwhから位置推定用電流振幅の三相交流成分Iuhac,Ivhac,Iwhacを抽出する。具体的に、交流成分抽出器404は、入力された位置推定用電流振幅Iuh,Ivh,Iwhのそれぞれから、直流成分抽出器405によって抽出された位置推定用電流振幅の直流成分Ihdcを減算することにより、位置推定用電流振幅の三相交流成分Iuhac,Ivhac,Iwhacを演算する。上記(1)及び(3)式を用いれば、位置推定用電流振幅の三相交流成分Iuhac,Ivhac,Iwhacは、下記(4)式で表すことができる。なお、位置推定用電流振幅の直流成分Ihdc及び位置推定用電流振幅の三相交流成分Iuhac,Ivhac,Iwhacの演算方法は、これに限定されず、ローパスフィルタ又はハイパスフィルタを用いてもよい。
Figure JPOXMLDOC01-appb-M000004
 推定位置演算器406は、位置推定用電流振幅の三相交流成分Iuhac,Ivhac,Iwhacが回転子位置θの関数、即ちsin(2θ)又はcos(2θ)の関数となることを利用して、第1の推定回転子位置θ を演算する。具体的に、第1の推定回転子位置θ は、上記(4)式に示される位置推定用電流振幅の三相交流成分Iuhac,Ivhac,Iwhacの何れか1つの信号を逆余弦演算することにより求めることができる。また、第1の推定回転子位置θ は、三相座標上で表現された位置推定用電流振幅の三相交流成分Iuhac,Ivhac,Iwhacを三相二相変換してから逆正接演算することにより求めることができる。また、第1の推定回転子位置θ は、位置推定用電流振幅の三相交流成分Iuhac,Ivhac,Iwhacのそれぞれのゼロクロス点を中心とした電気角60°毎に6つの区間に分け、各区間において位置推定用電流振幅の三相交流成分Iuhac,Ivhac,Iwhacのうちのゼロクロスするものを直線近似することによって演算することも可能である。
 推定位置演算器407は、位置推定用電流振幅の直流成分Ihdcに基づいて、第2の推定回転子位置θ を演算する。図4には、最大トルク(Maximum Torque Per Ampere:MTPA)制御における位置推定用電流振幅の直流成分Ihdcの振る舞いが、電流位相φと、駆動電流ベクトルの大きさ|idqf|とに対する関係で示されている。なお、MTPA制御は、同一のトルクを発生させる駆動電流ベクトルのうちで、駆動電流ベクトルの大きさ|idqf|を最小にする制御である。駆動電流ベクトルの大きさ|idqf|は、三相座標上の回転電機駆動電流iuf,ivf,iwfを回転座標変換し、d軸成分とq軸成分とに分けたときの、d軸成分とq軸成分との二乗和の平方根に相当する。図4では、駆動電流ベクトルの大きさ|idqf|が定格電流の20%から100%電流までの間で、20%の電流刻みで与えられている。また、図4において、電流位相φは、d軸を基準とした進み位相として定義されている。
 実施の形態1において、回転電機1に対するMTPA制御では、電流位相φの制御範囲は、45°から55°を想定する。図4において、駆動電流ベクトルの大きさ|idqf|が定格電流の60%以上であり、且つ、電流位相φが45°から55°の範囲では、位置推定用電流振幅の直流成分Ihdcは、電流位相φに対して単調減少している。実施の形態1における推定位置演算器407では、この特性を利用して、第2の推定回転子位置θ を演算する。
 なお、回転子位置の推定に必要な特性は、位置推定用電流振幅の直流成分Ihdcと駆動電流ベクトルとの一意性であり、位置推定用電流振幅の直流成分Ihdcと電流位相φとの関係は、単調減少に限定されない。
 図5には、図4に示した駆動電流ベクトルの大きさ|idqf|のプロットのうち、駆動電流ベクトルの大きさ|idqf|が100%電流、即ち定格電流である場合の、位置推定用電流振幅の直流成分Ihdcが示されている。また、図5において、電流位相指令値をφとし、電流位相指令値φで駆動した場合の位置推定用電流振幅の直流成分Ihdcを、Ihdc(φ)と表記する。また、回転子位置の位置推定誤差をΔθとし、電流位相指令値φが位置推定誤差Δθだけずれた実際の電流位相φで駆動した場合の位置推定用電流振幅の直流成分IhdcをIhdc(φ)と定義する。このように定義すると、位置推定誤差Δθは、Ihdc(φ)とIhdc(φ)との差に比例する。本明細書では、Ihdc(φ)とIhdc(φ)との差を「Ihdc誤差」と呼ぶ。
 推定位置演算器407を位相同期回路(Phase Locked Loop:PLL)で構成する場合、推定位置演算器407は、Ihdc誤差がゼロとなるようにPLLを動作させることで、第2の推定回転子位置θ を演算する。なお、PLLは、Ihdc誤差がゼロとなるような構成であればよく、比例積分器又は比例積分積分器が例示される。比例積分積分器は、比例積分器の後段に更に積分器を有する構成である。なお、推定位置演算器407では、使用を想定するトルク又は回転電機電流の範囲に対応付けて、電流位相φと位置推定用電流振幅の直流成分Ihdcとの関係が予め記憶されているものとする。
 推定位置演算器408は、図6に示すように、位置推定誤差演算器4080及びPLL4081を備える。前述したように、推定位置演算器408は、回転子1aの突極性により生じる誘起電圧により演算される鎖交磁束に基づいて、第3の推定回転子位置θ を推定する演算器である。
 位置推定誤差演算器4080は、三相座標上の回転電機駆動電圧指令vuf ,vvf ,vwf と、三相座標上の回転電機駆動電流iuf,ivf,iwfとに基づいて、回転子位置の位置推定誤差“-(θ -θ)”を演算する。PLL4081は、回転子位置の位置推定誤差“-(θ -θ)”から第3の推定回転子位置θ を演算する。
 位置推定誤差演算器4080は、三相二相変換器40800,40801、回転座標変換器40802、鎖交磁束インダクタンス交流分演算器40803、鎖交磁束インダクタンス交流分推定器40804及び回転子位置推定誤差演算器40805を備える。
 三相二相変換器40800は、三相座標上の回転電機駆動電圧指令vuf ,vvf ,vwf を、静止二相座標上の回転電機駆動電圧指令vαf ,vβf に変換する。
 三相二相変換器40801は、三相座標上の回転電機駆動電流iuf,ivf,iwfを、静止二相座標上の回転電機駆動電流iαf,iβfに変換する。
 回転座標変換器40802は、推定回転子位置θ を用いて座標変換を行い、静止二相座標上の回転電機駆動電流iαf,iβfを回転二相座標上の回転電機駆動電流idf,iqfに変換する。
 PLL4081は、回転子位置の位置推定誤差“-(θ -θ)”がゼロとなるように、PLL動作することで、第3の推定回転子位置θ を演算する。PLL4081は、位置推定誤差“-(θ -θ)”がゼロとなるように構成されていればよく、比例積分器又は比例積分微分器を用いてもよい。
 次に、鎖交磁束インダクタンス交流分演算器40803の動作について説明する。まず、回転電機1のモデルは、静止二相座標において、下記(5)式から(8)式で表される。
Figure JPOXMLDOC01-appb-M000005
Figure JPOXMLDOC01-appb-M000006
Figure JPOXMLDOC01-appb-M000007
Figure JPOXMLDOC01-appb-M000008
 上記(5)式のvdqはd軸の回転電機電圧vとq軸の回転電機電圧vからなるベクトルであり、idqはd軸回転電機電流iとq軸回転電機電流iからなるベクトルであり、Rは回転電機1の巻線抵抗であり、ωはモデルを表す座標の回転角速度であり、上記(5)式のψdqは鎖交磁束である。上記(6)式のJは、変換行列である。上記(5)式の鎖交磁束ψdqは上記(7)式のように表すことができる。上記(7)式のLdqは、Lsdcと、Lmacと、回転子位置の電気角θとにより、上記(8)式のように行列で表すことができる。上記(8)式のLsdcは、回転子位置によって変化しないインダクタンス直流分であり、Lmacは回転子位置によって変化するインダクタンス交流分である。なお、インダクタンスの変化は一般的に、回転子位置の電気角θに対し、2θの正弦関数又は余弦関数で表される。
 上記(7)式と、上記(8)式とから、鎖交磁束ψdqは、下記(9)式で表される。
Figure JPOXMLDOC01-appb-M000009
 上記(9)式の第1項は、回転子位置によって変化しないインダクタンス直流分Lsdcによる項である。また、上記(9)式の第2項は、回転子位置によって変化するインダクタンス交流分Lmacによる項であり、この項が「鎖交磁束インダクタンス交流分」である。即ち、鎖交磁束インダクタンス交流分は、インダクタンス交流分と回転電機電流とによって生成される鎖交磁束である。
 鎖交磁束インダクタンス交流分演算器40803は、鎖交磁束インダクタンス交流分を演算するため、以下の演算を行う。まず、鎖交磁束インダクタンス交流分演算器40803は、下記(10)式を用いて回転電機の鎖交磁束ψdqを演算する。
Figure JPOXMLDOC01-appb-M000010
 上記(10)式において、vαβ はα軸の回転電機電圧指令vα と、β軸の回転電機電圧指令vβ からなるベクトルである。
 また、上記(10)式の積分は、ラプラス変換におけるs領域において、下記(11)式で表される。
Figure JPOXMLDOC01-appb-M000011
 回転電機1の鎖交磁束ψαβを積分で演算する場合、通常は初期値が不明である。このため、回転電機1の鎖交磁束ψαβの基本波周波数成分に対して十分にカットオフ周波数の低いハイパスフィルタを利用する。ここで、ハイパスフィルタの伝達関数は、カットオフ周波数をωhpfとして、下記(12)式で表される。
Figure JPOXMLDOC01-appb-M000012
 よって、上記(11)式で表される鎖交磁束ψαβを上記(12)式のハイパスフィルタに通すと、フィルタ適用後の鎖交磁束ψ hpfαβは、下記(13)式で計算される。
Figure JPOXMLDOC01-appb-M000013
 また、上記(13)式を変形すると、下記(14)式が得られる。
Figure JPOXMLDOC01-appb-M000014
 更に、鎖交磁束インダクタンス交流分演算器40803は、第3の推定回転子位置θ を用いて、静止二相座標上での鎖交磁束ψ hpfαβを、回転二相座標上での鎖交磁束ψ hpfdqへ座標変換する。回転座標上での鎖交磁束インダクタンス交流分ψ acdqcalcは、上記(9)式に従って、下記(15)式で演算される。
Figure JPOXMLDOC01-appb-M000015
 上記(15)式で演算される鎖交磁束インダクタンス交流分ψ acdqcalcは、以下では「鎖交磁束インダクタンス交流分演算値」と呼ぶ。
 鎖交磁束インダクタンス交流分推定器40804は、下記(16)式に示すように、推定回転子位置θ と回転電機電流idqとを用いて、上記(9)式の第2項である鎖交磁束インダクタンス交流分を推定する。
Figure JPOXMLDOC01-appb-M000016
 ここで、上記(16)式において、回転子位置の推定値θ と真値θがおよそ等しい場合、上記(16)式は下記(17)式に示すように簡略化される。
Figure JPOXMLDOC01-appb-M000017
 上記(17)式のψ acdqは、鎖交磁束インダクタンス交流分推定器40804で演算される鎖交磁束インダクタンス交流分の推定値である。以下、この推定値を「鎖交磁束インダクタンス交流分推定値」と呼ぶ。
 位置推定誤差演算器4080は、鎖交磁束インダクタンス交流分演算値ψ acdqcalcと鎖交磁束インダクタンス交流分推定値ψ acdqとを用いて、回転子位置の位置推定誤差“-(θ -θ)”を演算する。ここで、鎖交磁束インダクタンス交流分演算値ψ acdqcalcと鎖交磁束インダクタンス交流分推定値ψ acdqとの外積は、上記(15)式、即ち上記(9)式の第2項を演算した値と、上記(16)式とを用いて、下記(18)式で表される。
Figure JPOXMLDOC01-appb-M000018
 そして、回転子位置の推定値と真値とがおよそ等しい、即ちθ ≒θとすると、回転子位置の推定誤差は下記(19)式で演算できる。
Figure JPOXMLDOC01-appb-M000019
 以上が、推定位置演算器408による演算処理である。なお、位置推定に用いる回転電機電圧指令及び回転電機電流は、それぞれ回転電機駆動電圧指令及び回転電機駆動電流とする。
 図3に戻り、推定位置切替器409は、第1の推定回転子位置θ 、第2の推定回転子位置θ 及び第3の推定回転子位置θ のうちの何れか1つを選択、即ち切り替えて出力する。或いは、推定位置切替器409は、第1の推定回転子位置θ 、第2の推定回転子位置θ 及び第3の推定回転子位置θ のうちの少なくとも2つの推定回転子位置の情報を選択し、これらを予め設定した割合で加重平均した値を推定回転子位置θとして出力する。このように、推定位置切替器409は、回転子位置の推定情報を選択又は切り替えて出力する。
 以上のように、実施の形態1によれば、位置推定器4は、第1の推定回転子位置θ 、第2の推定回転子位置θ 及び第3の推定回転子位置θ のうちの何れか1つを選択して出力する。或いは、位置推定器4は、第1の推定回転子位置θ 、第2の推定回転子位置θ 及び第3の推定回転子位置θ のうちの少なくとも2つの推定回転子位置を用いて合成した推定回転子位置θを出力する。これにより、所望の推定回転子位置を用いた位置センサレス制御が可能となり、動作範囲のトルク速度領域においての位置センサレス制御が可能となる。
実施の形態2.
 実施の形態1は、第1の推定回転子位置θ 、第2の推定回転子位置θ 及び第3の推定回転子位置θ のうちの何れか1つを選択することで得られる位置情報、又は、少なくとも2つの推定回転子位置を用いて合成することで得られる位置情報を推定回転子位置θとして出力する実施の形態であった。ここで、第1の推定回転子位置θ は、突極方式によって推定される位置情報であり、第2の推定回転子位置θ は、磁気飽和方式によって推定される位置情報であり、第3の推定回転子位置θ は、誘起電圧及び鎖交磁束方式によって推定される位置情報である。
 一方、回転子位置の推定方式には、それぞれの特徴がある。前述の通り、突極方式は、低速域且つ磁気飽和の程度が小さい駆動領域で推定精度が高く、磁気飽和方式は、低速域且つ磁気飽和の程度が大きい駆動領域で推定精度が高く、誘起電圧及び鎖交磁束方式は、高速域で推定精度が高いという特徴がある。実施の形態1では、推定される位置情報と駆動領域との関係性については、特に触れていなかった。そこで、実施の形態2では、回転子位置に加え、回転子1aの回転速度を推定し、推定した速度情報と、磁気飽和情報とに基づいて、推定した回転子位置の情報を切り替える実施の形態について説明する。なお、磁気飽和情報は、回転電機1の磁気飽和の程度と相関のある情報として定義される。
 図7は、実施の形態2に係る回転電機の制御装置100Aの構成例を示す図である。実施の形態2に係る制御装置100Aは、図1に示す実施の形態1に係る制御装置100の構成において、位置推定器4が位置推定器4Aに置換されている。その他の構成については、実施の形態1の構成と同一又は同等であり、同一又は同等の構成部には同一の符号を付して、重複する説明は省略する。
 図8は、図7に示す位置推定器4Aの細部の構成例を示す図である。図8に示す実施の形態2における位置推定器4Aは、図3に示す実施の形態1における位置推定器4の構成において、推定位置演算器406,407,408及び推定位置切替器409が、それぞれ推定位置速度演算器406A,407A,408A及び推定位置切替器409Aに置換されている。推定位置速度演算器406A,407A,408Aのそれぞれを符号無しで区別する場合には、それぞれを「第1の推定位置速度演算器」、「第2の推定位置速度演算器」、及び「第3の推定位置速度演算器」と呼ぶ。
 推定位置速度演算器406Aは、位置推定用電流振幅の三相交流成分Iuhac,Ivhac,Iwhacに基づいて、第1の推定回転子位置θ と、第1の推定速度ω とを演算する。推定位置速度演算器406Aでは、実施の形態1と同様に、突極方式が用いられる。即ち、第1の推定速度ω は、位置推定用電流から検出した回転子1aの突極性に基づいて演算される推定速度である。
 推定位置速度演算器407Aは、位置推定用電流振幅の直流成分Ihdcに基づいて、第2の推定回転子位置θ と、第2の推定速度ω とを演算する。推定位置速度演算器407Aでは、実施の形態1と同様に、磁気飽和方式が用いられる。即ち、第2の推定速度ω は、位置推定用電流振幅の直流成分Ihdcに基づいて演算される推定速度である。
 推定位置速度演算器408Aは、三相座標上の回転電機駆動電流iuf,ivf,iwfと、回転子の突極性により生じる誘起電圧により演算される鎖交磁束と、に基づいて第3の推定回転子位置θ と、第3の推定速度ω とを演算する。推定位置速度演算器408Aでは、実施の形態1と同様に、誘起電圧及び鎖交磁束方式が用いられる。即ち、第3の推定速度ω は、三相座標上の回転電機駆動電流iuf,ivf,iwfと、回転子1aの突極性により生じる誘起電圧により演算される鎖交磁束と、に基づいて演算される推定速度である。
 推定位置切替器409Aには、τ,τ ,τ ,τ ,i ,iff ,iが入力される。τは、前述したトルク指令値であり、iは、駆動電流である。τ は、トルク指令値τに制御系の遅れを考慮したフィルタを通して得られる値である。τ は、回転電機の数学モデルに基づき、駆動電流iと回転電機パラメータとから演算される推定トルクである。τ は、駆動電流iを引数としたルックアップテーブルにより得られる推定トルクである。i は、回転電機駆動電流指令である。iff は、回転電機駆動電流指令i に制御系の遅れを考慮したフィルタを通して得られる値である。駆動電流iでは、駆動電流iの周波数が利用される。推定位置切替器409Aでは、これらの入力情報のうちの少なくとも1つが、磁気飽和情報として使用される。
 図9は、図8に示す推定位置速度演算器406Aの細部の構成例を示す図である。推定位置速度演算器406Aは、図9に示すように、推定位置演算器406及び推定速度演算器4060を備える。図9に示す推定位置演算器406は、図3に示した推定位置演算器406と同等の構成部である。また、推定速度演算器4060は、推定位置演算器406が演算した第1の推定回転子位置θ を擬似微分することにより第1の推定速度ω を演算する。推定位置速度演算器406Aは、演算した第1の推定回転子位置θ 及び第1の推定速度ω を出力する。なお、擬似微分は、入力値に微分演算とフィルタ処理とを施したものである。擬似微分の処理器は、微分器及びローパスフィルタにより実現することができる。
 図10は、図8に示す推定位置速度演算器407Aの細部の構成例を示す図である。推定位置速度演算器407Aは、図10に示すように、ルックアップテーブル(Look Up Table:LUT)4070、LUT4071、減算器4072、除算器4073、比例積分器4074及び積分器4075を備える、2型の演算器として構成される。
 LUT4070には、トルク指令値τと、そのトルク指令値τで駆動した際に得られる位置推定用電流振幅の直流成分IhdcであるIhdc とが記憶されている。推定位置速度演算器407Aは、LUT4070に記憶されたテーブル値Ihdc を参照する。
 減算器4072は、位置推定用電流振幅の直流成分Ihdcとテーブル値Ihdc との偏差ΔIを演算する。
 LUT4071には、トルク指令値τから位置推定誤差Δθを換算するための係数KIhθ2が記憶されている。推定位置速度演算器407Aは、LUT4071に記憶された係数KIhθ2を参照する。
 除算器4073は、偏差ΔIをKIhθ2で除することにより位置推定誤差Δθを得る。比例積分器4074は、より位置推定誤差Δθに基づいて第2の推定速度ω を演算する。積分器4075は、第2の推定速度ω に基づいて第2の推定回転子位置θ を演算する。
 なお、図10は、第2の推定速度ω を先に演算し、第2の推定速度ω に基づいて第2の推定回転子位置θ を演算する構成であるが、これに限定されない。これとは逆に、第2の推定回転子位置θ を先に演算し、第2の推定回転子位置θ を擬似微分することで第2の推定速度ω を演算する構成であってもよい。
 また、図10は、LUT4070及びLUT4071の引数はトルク指令値τとする例であるが、これに限定されない。トルク指令値τに代えて、回転電機駆動電流指令i をLUT4070及びLUT4071の引数としてもよい。
 図11は、図8に示す推定位置速度演算器408Aの細部の構成例を示す図である。推定位置速度演算器408Aは、図11に示すように、位置推定誤差演算器4080及びPLL4081Aを備える。
 図11に示す位置推定誤差演算器4080は、図6に示した位置推定誤差演算器4080と同等の構成部である。PLL4081Aは、位置推定誤差“-(θ^3-θ)”がゼロとなるように演算器が構成される。具体的に、PLL4081Aは、図10に示す推定位置速度演算器407Aと同様に、比例積分器4081Aaと、積分器4081Abとを備えた2型の演算器として構成される。
 PLL4081Aにおいて、比例積分器4081Aaは、位置推定誤差“-(θ -θ)”に基づいて第3の推定速度ω を演算する。積分器4081Abは、第3の推定速度ω に基づいて第3の推定回転子位置θ を演算する。
 なお、図11のPLL4081Aは、第3の推定速度ω を先に演算し、第3の推定速度ω に基づいて第3の推定回転子位置θ を演算する構成であるが、これに限定されない。これとは逆に、第3の推定回転子位置θ を先に演算し、第3の推定回転子位置θ を擬似微分することで第3の推定速度ω を演算するようにPLL4081Aが構成されていてもよい。
 図12は、図8に示す推定位置切替器409Aの細部の構成例を示す図である。推定位置切替器409Aは、図12に示すように、推定速度選択器411A及び推定位置選択器410Aを備える。推定速度選択器411Aは、第1の推定速度ω 、第2の推定速度ω 及び第3の推定速度ω のうちの何れか1つを選択し、選択した推定速度ωを速度情報として推定位置選択器410Aに出力する。推定位置選択器410Aは、推定速度選択器411Aから出力される速度情報と、前述した磁気飽和情報とに基づいて、回転子位置の推定情報を選択又は切り替えて出力する。
 推定位置選択器410Aにおいて、推定回転子位置θの選択又は切替に用いる速度情報は、第1の推定速度ω 、第2の推定速度ω 及び第3の推定速度ω のうちの何れか1つである。即ち、推定速度選択器411Aは、第1の推定速度ω 、第2の推定速度ω 及び第3の推定速度ω のうちの何れか1つの推定速度を選択して推定位置選択器410Aに出力する。
 図13は、図8及び図12に示す推定位置切替器409Aにおける駆動領域の説明に供する図である。図13において、横軸は回転電機の回転速度であり、縦軸は回転電機のトルクである。回転速度は、速度情報の例示である。また、トルクは、前述した磁気飽和情報と相関がある情報である。
 図13において、(1)は突極方式により駆動する駆動領域、(2)は磁気飽和方式により駆動する駆動領域、(3)は誘起電圧及び鎖交磁束方式により駆動する駆動領域である。また、図13における破線は、前述した回転子位置の各推定方式との境界線を意味し、境界判定の閾値となる。この境界線は、磁気飽和情報及び速度情報により規定される。即ち、領域の判定には磁気飽和情報と、速度情報とが用いられ、その領域に規定された回転子位置の推定方式が選択される。
 具体的に説明すると、推定位置切替器409Aは、回転速度が閾値より小さく、且つ、トルクが閾値より小さい場合は、第1の推定回転子位置θ を選択して出力する。また、推定位置切替器409Aは、回転速度が閾値より小さく、且つ、トルクが閾値より大きい場合は、第2の推定回転子位置θ を選択して出力する。そして、推定位置切替器409Aは、回転速度が閾値より大きい場合は、第3の推定回転子位置θ を選択して出力する。
 なお、図13では、切り替えの境界線を直線で示したが、直線である必要はなく、曲線であってもよい。
 以上のように、実施の形態2によれば、回転子の磁気飽和情報及び回転速度に基づいて、適切な回転子位置の推定方式を切り替えることができる。これにより、動作範囲のトルク速度領域において、所望する位置センサレス制御が可能となる。
実施の形態3.
 実施の形態2では、推定位置切替器409Aにおいて、回転子1aの磁気飽和情報及び回転速度に基づいて、回転子位置の推定情報を選択又は切り替える実施の形態であった。一方、回転子位置の推定情報を選択又は切り替える際には、切り替えの前後において、推定情報が不連続に変化する可能性があり、切り替え時のショックが課題となる。そこで、実施の形態3では、切り替え時のショックを小さくできる実施の形態について説明する。
 図14は、実施の形態3に係る回転電機の制御装置100Bの構成例を示す図である。図14に示す実施の形態3に係る制御装置100Bは、図7に示す実施の形態2に係る制御装置100Aの構成において、位置推定器4Aが位置推定器4Bに置換されている。その他の構成については、実施の形態2の構成と同一又は同等であり、同一又は同等の構成部には同一の符号を付して、重複する説明は省略する。
 図15は、図14に示す位置推定器4Bの細部の構成例を示す図である。図15に示す実施の形態3における位置推定器4Bは、図8に示す実施の形態2における位置推定器4Aの構成において、推定位置切替器409Aが、推定位置切替器409Bに置換されている。その他の構成については、図8の構成と同一又は同等であり、同一又は同等の構成部には同一の符号を付して、重複する説明は省略する。
 図16は、図15に示す推定位置切替器409Bの細部の構成例を示す図である。図16に示す実施の形態3における推定位置切替器409Bは、図12に示す実施の形態2における推定位置切替器409Aの構成において、推定位置選択器410Aが推定位置合成器410Bに置換されている。
 推定位置切替器409Bにおいて、推定位置合成器410Bは、磁気飽和情報及び速度情報に基づいて、第1の推定回転子位置θ 、第2の推定回転子位置θ 及び第3の推定回転子位置θ のうちの少なくとも2つの推定回転子位置を用いて合成した位置情報を推定回転子位置θとして出力する。
 推定位置切替器409Bにおいて、推定回転子位置の合成に用いる磁気飽和情報は、実施の形態2で用いる磁気飽和情報と同様である。また、推定回転子位置の切り替えに用いる速度情報は、実施の形態2で用いる速度情報と同様である。
 推定位置合成器410Bは、磁気飽和情報と速度情報とに基づいて、第1の推定回転子位置θ 、第2の推定回転子位置θ 及び第3の推定回転子位置θ を加重平均により合成する。図17は、図15及び図16に示す推定位置切替器409Bにおける駆動領域の説明に供する図である。図17において、横軸は回転電機の回転速度であり、縦軸は回転電機のトルクである。トルクは、実施の形態2と同様に、磁気飽和情報に相関する情報の例示である。
 図17において、(1)は突極方式のみにより駆動する駆動領域、(2)は磁気飽和方式のみにより駆動する駆動領域、(3)は誘起電圧及び鎖交磁束方式のみにより駆動する駆動領域である。また、(1)(2)のように、複数の数字で表記される領域は、それらの数字に対応する推定方式を使用する駆動領域である。
 例えば、(1)(2)は突極方式及び磁気飽和方式により駆動される駆動領域である。この駆動領域では、突極方式により推定された第1の推定回転子位置θ と、磁気飽和方式により推定された第2の推定回転子位置θ とによって、出力すべき推定回転子位置θが合成される。なお、駆動領域の境界は、図13と同様に破線で示している。また、図17では、図13と同様に、駆動領域の境界は直線で示されているが、直線である必要はなく、曲線であってもよい。
 図17において、Wspdは速度情報に基づく重みであり、Wmsは磁気飽和情報に基づく重みである。重みWspd及び重みWmsは、それぞれが速度情報及び磁気飽和情報に応じて、0から1の範囲で値が変化する。具体的に、推定位置合成器410Bは、下記の(20)式及び(21)式を用いて、図17で定義された駆動領域において、推定回転子位置θを合成する。
Figure JPOXMLDOC01-appb-M000020
Figure JPOXMLDOC01-appb-M000021
 上記(21)式におけるθ は、「突極方式と磁気飽和方式との重みWmsによる加重平均値」である。以下、この加重平均値θ を「第4の推定回転子位置」と呼ぶ。
 また、上記(20)式におけるθ は、「第4の推定回転子位置θ と、誘起電圧及び鎖交磁束方式のみによる演算値である第3の推定回転子位置θ との重みWspdによる加重平均値」である。以下、この加重平均により合成したものを「第5の推定回転子位置」と呼ぶ。
 例えば、Wspd=1の場合、(20)式はθ =θ となり、第3の推定回転子位置θ が、推定位置合成器410Bから出力される推定回転子位置θとなる。また、Wspd=0の場合、(20)式はθ =θ となり、(21)式に従って、重みWmsの値に応じて重み付けされたθ 、即ち、突極方式と磁気飽和方式とによる加重平均値である第4の推定回転子位置θ が、推定位置合成器410Bから出力される推定回転子位置θとなる。また、0<Wspd<1の場合、重みWmsの値に応じて、突極方式と、磁気飽和方式と、誘起電圧及び鎖交磁束方式とのうちの少なくとも2つの方式の加重平均によって合成された第5の推定回転子位置θ が推定位置合成器410Bから出力される。なお、(20)式と(21)式の演算の順序を入れ替えてもよい。また、(20)式に(21)式を代入した演算式を用いて第4の推定回転子位置を演算することなく第5の推定回転子位置を演算してもよい。更に、図17で定義された駆動領域において、定義された第1の推定回転子位置と第2の推定回転子位置と第3の推定回転子位置との組み合わせにより推定回転子位置θ^を合成する演算式であれば,推定回転子位置θ^の演算式は(20)式(21)式である必要は無い。
 なお、図14の制御器5における位置推定用電圧演算部5b、及び図15の位置推定器4Bにおける推定位置速度演算器406A,407A,408Aは、常時動作させる必要はなく、図17で規定した駆動領域に応じて停止させてもよい。例えば、図17における数字に対応する推定方式の構成部のみを動作させ、記載されていない数字に対応する推定方式の構成部は動作を停止させることが可能である。また、図17における(3)の領域においては、位置推定用電圧演算部5bの動作を停止させ、その出力を0とすることができる。これにより、演算処理に伴う消費電力の低減が可能になる。
 また、推定位置速度演算器406A,407A,408Aが動作を停止している状態から、動作を開始する際には、動作を開始する位置速度推定器の初期値として、動作を行っている他の位置速度推定器の推定回転子位置、推定速度、積分器出力又は推定位置合成器410Bの出力θ、及び推定速度選択器411Aの出力ωを与えることができることは言うまでもない。
 以上のように、実施の形態3によれば、磁気飽和情報と速度情報とに基づいて、実施の形態2で得られた3つの回転子位置の推定情報を加重平均により合成して出力する。これにより、実施の形態2の効果に加え、回転子位置の推定情報を切り替える際のショックを低減することが可能となる。
実施の形態4.
 実施の形態2及び実施の形態3では、回転子位置の推定方式を切り替える際に用いる磁気飽和情報として、回転電機1に対するトルク指令値τ、トルク指令値τのフィルタ出力値τ 、回転電機の数学モデルにより得られる推定トルクτ 、駆動電流iを引数としたルックアップテーブルにより得られる推定トルクτ 、回転電機1に対する回転電機駆動電流指令i 、回転電機駆動電流指令i のフィルタ出力値iff 、及び駆動電流iの周波数のうちの少なくとも1つを用いる。一方、制御器5の制御において、トルク制御誤差、電流制御誤差、位置推定誤差が生じた場合、実施の形態2及び実施の形態3で用いる磁気飽和情報では、磁気飽和の程度を正確に把握できないことが想定される。そこで、実施の形態4では、回転子位置の推定方式を切り替える際の磁気飽和情報として、前述した位置推定用電流振幅の直流成分Ihdcを使用する。なお、位置推定用電流振幅の直流成分Ihdcは、位置推定用電流振幅Iuh,Ivh,Iwhのうち、回転子位置によって変化しない成分である。即ち、実施の形態4は、回転子位置によって変化しない位置推定用電流振幅の直流成分Ihdcを磁気飽和情報として使用する実施の形態である。
 図18は、実施の形態4に係る回転電機の制御装置100Cの構成例を示す図である。図18に示す実施の形態4に係る制御装置100Cは、図7に示す実施の形態2に係る制御装置100Aの構成において、位置推定器4Aが位置推定器4Cに置換されている。その他の構成については、実施の形態2の構成と同一又は同等であり、同一又は同等の構成部には同一の符号を付して、重複する説明は省略する。
 図19は、図18に示す位置推定器4Cの細部の構成例を示す図である。図18に示す実施の形態4における位置推定器4Cは、図8に示す実施の形態2における位置推定器4Aの構成において、推定位置切替器409Aが、推定位置切替器409Cに置換されている。また、図8に示す磁気飽和情報に代えて、直流成分抽出器405によって抽出される位置推定用電流振幅の直流成分Ihdcが磁気飽和情報として推定位置切替器409Cに入力される構成である。なお、その他の構成については、図8の構成と同一又は同等であり、同一又は同等の構成部には同一の符号を付して、重複する説明は省略する。
 図20は、図19に示す推定位置切替器409Cの細部の構成例を示す図である。図20に示す実施の形態4における推定位置切替器409Cは、図12に示す実施の形態2における推定位置切替器409Aの構成において、推定位置選択器410Aが推定位置合成器410Cに置換されている。また、図12に示す磁気飽和情報に代えて、位置推定用電流振幅の直流成分Ihdcが磁気飽和情報として推定位置合成器410Cに入力されている。
 推定位置合成器410Cは、磁気飽和情報と速度情報とに基づいて、第1の推定回転子位置θ 、第2の推定回転子位置θ 及び第3の推定回転子位置θ を加重平均により合成する。
 図21は、図19及び図20に示す推定位置切替器409Cにおける駆動領域の説明に供する図である。図21において、横軸、縦軸及び駆動領域の定義は実施の形態2のものと同様である。回転電機の回転速度であり、縦軸は回転電機のトルクである。縦軸のトルクは、実施の形態2と同様に、位置推定用電流振幅の直流成分Ihdcに相関する情報の例示である。縦軸を区分するWms(Ihdc)は、位置推定用電流振幅の直流成分Ihdc及び速度情報、即ち磁気飽和情報及び速度情報により規定される重みである。なお、図21では、図13及び図17と同様に、駆動領域の境界は直線で示されているが、直線である必要はなく、曲線であってもよい。
 図21に示すように、重みWspdは、速度情報に応じて、0から1の範囲で値が変化する。また、重みWms(Ihdc)は、磁気飽和情報及び速度情報に応じて、0から1の範囲で値が変化する。推定位置合成器410Cは、下記(22)式及び(23)式を用いて、図21で規定された駆動領域において、推定回転子位置θを合成する。
Figure JPOXMLDOC01-appb-M000022
Figure JPOXMLDOC01-appb-M000023
 上記(23)式は、上記(21)式のWmsをWms(Ihdc)に置換したものである。従って、上記(23)式におけるθ は、Wms(Ihdc)によって重み付けされる点を除き、「突極方式と磁気飽和方式との加重平均値」である「第4の推定回転子位置」を意味している。
 また、上記(22)式は、上記(20)式と同一の式であり、重みWms(Ihdc)によって重み付けされた第4の推定回転子位置θ と、誘起電圧及び鎖交磁束方式のみによる演算値である第3の推定回転子位置θ との重みWspdによる加重平均値である第5の推定回転子位置θ である。なお、(22)式と(23)式の演算の順序を入れ替えてもよい。また、(22)式に(23)式を代入した演算式を用いて第4の推定回転子位置を演算することなく第5の推定回転子位置を演算してもよい。更に、図17で定義された駆動領域において、定義された第1の推定回転子位置と第2の推定回転子位置と第3の推定回転子位置との組み合わせにより推定回転子位置θ^を合成する演算式であれば、推定回転子位置θ^の演算式は(22)式、(23)式である必要は無い。
 なお、実施の形態4では、駆動領域の境界と、重みWms(Ihdc)とは、位置推定用電流振幅の直流成分Ihdcに基づいて規定しているが、これに限定されない。位置推定用電流振幅の直流成分Ihdcに加え、実施の形態2及び実施の形態3で用いた磁気飽和情報を併用して規定してもよい。
 また、実施の形態4では、図21に示すように、加重平均により複数の推定方式を併用する駆動領域を規定しているが、これに限定されない。位置推定用電流振幅の直流成分Ihdc及び速度情報に基づいて、第1の推定回転子位置θ 、第2の推定回転子位置θ 及び第3の推定回転子位置θ のうちの何れか1つが選択されるように、駆動領域を規定してもよい。
 なお、位置推定用電流振幅の直流成分Ihdcは、各推定方式の処理を行う演算器の異常判定、及び、推定方式を切り替える際の可否判定にも利用できることは言うまでもない。なぜなら、推定位置速度演算器407Aは、推定原理上、図4に示すような、駆動電流ベクトルと、位置推定用電流振幅の直流成分Ihdcとの関係を事前に記憶しているためである。即ち、推定位置速度演算器407Aは、駆動条件に応じて得られる位置推定用電流振幅の直流成分Ihdcを把握しているので、記憶している直流成分Ihdcと実際の直流成分Ihdcとを比較することにより、異常判定が可能である。また、駆動電流ベクトルと直流成分Ihdcとの関係から、磁気飽和方式によって駆動可能な駆動領域であるか否かを把握できるので、実際の直流成分Ihdcを用いることにより、推定方式を切り替える際の可否判定が可能となる。
 なお、図18の制御器5における位置推定用電圧演算部5b、及び図19の位置推定器4Cにおける推定位置速度演算器406A,407A,408Aは、常時動作させる必要はなく、図21で規定した駆動領域に応じて停止させてもよい。例えば、図21における数字に対応する推定方式の構成部のみを動作させ、記載されていない数字に対応する推定方式の構成部は動作を停止させることが可能である。また、図21における(3)の領域においては、位置推定用電圧演算部5bの動作を停止させ、その出力を0とすることができる。これにより、演算処理に伴う消費電力の低減が可能になる。
 また、推定位置速度演算器406A,407A,408Aが動作を停止している状態から、動作を開始する際には、動作を開始する位置速度推定器の初期値として、動作を行っている他の位置速度推定器の推定回転子位置、推定速度、積分器出力又は推定位置合成器410Cの出力θ、及び推定速度選択器411Aの出力ωを与えることができることは言うまでもない。
 以上のように、実施の形態4によれば、トルク制御誤差、電流制御誤差、位置推定誤差によらず磁気飽和の程度を正確に把握することができる。このため、実施の形態3の効果に加え、回転子位置の推定精度を高めることが可能になる。
実施の形態5.
 実施の形態2から実施の形態4では、回転子位置の推定情報を切り替えるために用いる速度情報として、第1の推定速度ω 、第2の推定速度ω 及び第3の推定速度ω のうちの何れか1つを用いる実施の形態であった。一方、第1の推定速度ω 、第2の推定速度ω 及び第3の推定速度ω の推定精度は、駆動条件により異なる。例えば、第1の推定速度ω は、低速域且つ磁気飽和の程度が小さい駆動領域では推定精度が高く、第2の推定速度ω は、低速域且つ磁気飽和の程度が大きい駆動領域で推定精度が高く、第3の推定速度ω は、高速域で推定精度が高いという特徴がある。実施の形態2から実施の形態4では、推定される速度情報と駆動領域との関係性については、特に触れていなかった。そこで、実施の形態5では、回転子位置の推定情報を切り替えるために用いる速度情報の精度を高める実施の形態について説明する。
 図22は、実施の形態5に係る回転電機の制御装置100Dの構成例を示す図である。図22に示す実施の形態5に係る制御装置100Dは、図18に示す実施の形態4に係る制御装置100Cの構成において、位置推定器4Cが位置推定器4Dに置換されている。その他の構成については、実施の形態4の構成と同一又は同等であり、同一又は同等の構成部には同一の符号を付して、重複する説明は省略する。
 図23は、図22に示す位置推定器4Dの細部の構成例を示す図である。図23に示す実施の形態5における位置推定器4Dは、図19に示す実施の形態4における位置推定器4Cの構成において、推定位置切替器409Cが、推定位置切替器409Dに置換されている。その他の構成については、図8の構成と同一又は同等であり、同一又は同等の構成部には同一の符号を付して、重複する説明は省略する。
 図24は、図23に示す推定位置切替器409Dの細部の構成例を示す図である。図24に示す実施の形態5における推定位置切替器409Dは、図20に示す実施の形態4における推定位置切替器409Cの構成において、推定位置合成器410Cが推定位置合成器410Dに置換され、推定速度選択器411Aが推定速度合成器411Dに置換されている。また、位置推定用電流振幅の直流成分Ihdcが磁気飽和情報として推定位置合成器410Dに入力されると共に、推定速度合成器411Dの出力が自身にフィードバックされている。
 推定速度合成器411Dは、磁気飽和情報及び速度情報に基づいて、第1の推定速度ω 、第2の推定速度ω 及び第3の推定速度ω のうちの何れか1つを選択することで得られる速度情報、又は、少なくとも2つの推定速度を用いて、加重平均により合成した速度情報を推定速度ωとして出力する。推定速度ωを選択もしくは合成するための速度情報としては、推定速度合成器411Dの出力、即ち自身の出力を利用する。なお、推定速度ωに代えて、第1の推定速度ω 、第2の推定速度ω 及び第3の推定速度ω のうちの何れか1つを用いてもよい。
 また、実施の形態5における磁気飽和情報は、位置推定用電流振幅Iuh,Ivh,Iwhのうち、回転子位置によって変化しない位置推定用電流振幅の直流成分Ihdcとしているが、直流成分Ihdcに加え、もしくは直流成分Ihdcに代えて、実施の形態2及び実施の形態3で説明した磁気飽和情報を使用してもよい。従って、図24に示す推定位置合成器410Dは、図12に示す推定位置選択器410A、図16に示す推定位置合成器410B、及び図20に示す推定位置合成器410Cのうちの何れかに置換することができる。
 次に、推定速度合成器411Dから出力される推定速度ωの演算方法について、図21を参照して説明する。実施の形態5において、図21に示される数値は、各領域において演算される推定速度と見ることができる。例えば、図21において、(1)は第1の推定速度ω を演算する領域、(2)は第2の推定速度ω を演算する領域、(3)は第3の推定速度ω を演算する領域である。また、(1)(2)のように複数の数字で表記される領域は、それらの数字に対応する複数の推定速度を、前述した磁気飽和情報及び速度情報に基づいて合成する領域である。
 例えば、(1)(2)では、第1の推定速度ω と、第2の推定速度ω との加重平均が演算される。また、図中における破線は、何れの推定速度を用いるかを切り替える際の境界であり、この境界は、前述した磁気飽和情報及び速度情報により規定される。なお、この境界は直線で示されているが、直線である必要はなく、曲線であってもよい。
 図21に示すように、重みWspdは、速度情報に応じて、0から1の範囲で値が変化する。また、重みWms(Ihdc)は、磁気飽和情報及び速度情報に応じて、0から1の範囲で値が変化する。推定速度合成器411Dは、下記(24)式及び(25)式を用いて、図21で規定された領域において、推定速度ωを合成する。
Figure JPOXMLDOC01-appb-M000024
Figure JPOXMLDOC01-appb-M000025
 上記(25)式におけるω は、「突極方式と磁気飽和方式との重みWms(Ihdc)による加重平均値」である。以下、この加重平均値ω を「第4の推定速度」と呼ぶ。
 また、上記(24)式におけるω は、「第4の推定速度ω と、誘起電圧及び鎖交磁束方式のみによる演算値である第3の推定速度ω との重みWspdによる加重平均」である。以下、この加重樹平均値により合成したものを「第5の推定速度」と呼ぶ。
 例えば、Wspd=1の場合、(24)式はω =ω となり、「誘起電圧及び鎖交磁束方式のみによる演算値ω 」が、推定速度合成器411Dから出力される推定速度ωとなる。また、Wspd=0の場合、(24)式はω =ω となり、(25)式に従って、重みWms(Ihdc)の値に応じて重み付けされたω 、即ち「突極方式と磁気飽和方式とによる加重平均値ω 」が、推定速度合成器411Dから出力される推定速度ωとなる。また、0<Wspd<1の場合、重みWms(Ihdc)の値に応じて、突極方式と、磁気飽和方式と、誘起電圧及び鎖交磁束方式とのうちの少なくとも2つの方式の加重平均によって合成された推定速度ωが推定速度合成器411Dから出力される。
 また、実施の形態5では、図21に示すように、推定速度を加重平均により合成する領域を規定しているが、これに限定されない。位置推定用電流振幅の直流成分Ihdc及び速度情報に基づいて、第1の推定速度ω 、第2の推定速度ω 及び第3の推定速度ω のうちの何れか1つが選択されるように、領域を規定してもよい。なお、(24)式と(25)式の演算の順序を入れ替えてもよい。また、(24)式に(25)式を代入した演算式を用いて第4の推定速度を演算することなく第5の推定速度を演算してもよい。更に、図21で定義された駆動領域において、定義された第1の推定速度と第2の推定速度と第3の推定速度との組み合わせにより推定速度ω^を合成する演算式であれば、推定速度ω^の演算式は(24)式、(25)式である必要は無い。
 以上のように、実施の形態5によれば、推定位置切替時において、推定速度の急峻な変化を防止し、更に速度情報を正確に把握することができる。加えて、トルク制御誤差、電流制御誤差、位置推定誤差によらず磁気飽和の程度を正確に把握することができる。このため、実施の形態4よりも、回転子位置の推定精度を高めることが可能になる。
 次に、実施の形態1から実施の形態5に係る制御装置100,100A,100B,100C,100Dが備える各機能もしくはその一部を実現するためのハードウェアの構成について説明する。ここで言う各機能とは、電流検出器2、電圧印加器3、位置推定器4,4A,4B,4C,4D及び制御器5が有する機能である。
 図25は、実施の形態1から実施の形態5に係る回転電機の制御装置の各機能もしくはその一部を実現するための第1のハードウェア構成例を示す図である。図26は、実施の形態1から実施の形態5に係る回転電機の制御装置の各機能もしくはその一部を実現するための第2のハードウェア構成例を示す図である。図25には、専用処理回路1000のような専用のハードウェアにより上記の処理回路を実現する例が示される。図26にはプロセッサ1001及び記憶装置1002により上記の処理回路を実現する例が示される。
 図25に示す例では、図1、図7、図14、図18及び図22の電流検出器2、並びに電圧印加器3は、専用のハードウェアを用いてその機能が実現され、位置推定器4,4A,4B,4C,4D、及び制御器5は、専用処理回路1000により実現される。図25に示すように専用の処理回路を利用する場合、専用処理回路1000は、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ASIC、FPGA、又はこれらを組み合わせたものが該当する。上記の各機能のそれぞれを、処理回路で実現してもよいし、まとめて処理回路で実現してもよい。
 図26に示す例では、図1、図7、図14、図18及び図22の電流検出器2、並びに電圧印加器3は、専用のハードウェアを用いてその機能が実現され、位置推定器4,4A,4B,4C,4D、及び制御器5は、記憶装置1002に記録されたプログラムを実行するプロセッサ1001により実現される。なお、複数のプロセッサ1001と複数の記憶装置1002とが連携して、上述した各機能を実現してもよい。
 図26に示すように、プロセッサ1001及び記憶装置1002を利用する場合、上述した各機能は、ソフトウェア、ファームウェア又はこれらの組合せにより実現される。ソフトウェア又はファームウェアはプログラムとして記述され、記憶装置1002に記憶される。プロセッサ1001は、記憶装置1002に記憶されたプログラムを読み出して実行する。また、これらのプログラムは、各機能が実行される手順及び方法をコンピュータに実行させるものであるとも言える。
 記憶装置1002は、RAM(Random Access Memory)、ROM(Read Only Memory)、フラッシュメモリー、EPROM(Erasable Programmable Read Only Memory)、又はEEPROM(登録商標)(Electrically Erasable Programmable Read Only Memory)といった半導体メモリが該当する。半導体メモリは不揮発性メモリでもよいし、揮発性メモリでもよい。また、記憶装置1002は、半導体メモリ以外にも、磁気ディスク、フレキシブルディスク、光ディスク、コンパクトディスク、ミニディスク又はDVD(Digital Versatile Disc)が該当する。
 なお、実施の形態1から実施の形態5では、回転電機1を同期リラクタンスモータとして説明したが、これに限定されない。実施の形態1から実施の形態5に係る制御装置100,100A,100B,100C,100Dは、埋込磁石型同期モータ又は表面磁石型同期モータといった突極性を持つモータにも適応可能である。
 また、実施の形態1から実施の形態5では、位置推定用電流振幅の直流成分から回転子位置を推定しているが、位置推定の高精度化又は低騒音化のために位置推定用電圧指令の大きさを可変にする場合がある。この場合において、位置推定用電流振幅の直流成分の代わりに、上記の直流成分と位置推定用電圧指令との比、例えばインダクタンスの直流成分などを利用できることは自明のことである。
 また、実施の形態1から実施の形態5では、制御装置100,100A,100B,100C,100Dの制御器5がトルクを制御するものとして説明したが、これに限定されない。制御器5は、回転速度を制御する構成とすることもできる。
 また、実施の形態1から実施の形態5では、トルクに対する回転電機1の電流指令が電流実効値、即ち銅損が最小になるように選択されているが、これに限定されない。制御装置100,100A,100B,100C,100Dにおいて、鎖交磁束又は回転電機1の損失が最小になるように設定してよい。
 また、実施の形態1から実施の形態5では、電流検出器2が回転電機1の相電流を検出する構成例を説明したが、これに限定されない。制御装置100,100A,100B,100C,100Dは、相電流を検出することができればよく、電圧印加器3を構成する図示しないインバータに内蔵された電流センサにより、回転電機1の相電流を検出する構成としてもよい。
 なお、以上の実施の形態に示した構成は、本発明の内容の一例を示すものであり、別の公知の技術と組み合わせることも可能であるし、本発明の要旨を逸脱しない範囲で、構成の一部を省略、変更することも可能である。
 1 回転電機、1a 回転子、2 電流検出器、3 電圧印加器、4,4A,4B,4C,4D 位置推定器、5 制御器、5a 駆動電圧指令演算部、5b 位置推定用電圧演算部、6 電流指令演算器、7 電流制御器、8 回転座標逆変換器、9 二相三相変換器、11 駆動電流抽出器、12,40800,40801 三相二相変換器、13,40802 回転座標変換器、14 加算器、41,42 信号処理器、100,100A,100B,100C,100D 制御装置、401 位置推定用電流抽出器、402 駆動電流抽出器、403 位置推定用電流振幅演算器、404 交流成分抽出器、405 直流成分抽出器、406,407,408 推定位置演算器、406A,407A,408A 推定位置速度演算器、409,409A,409B,409C,409D 推定位置切替器、410A 推定位置選択器、410B,410C,410D 推定位置合成器、411A 推定速度選択器、411D 推定速度合成器、1000 専用処理回路、1001 プロセッサ、1002 記憶装置、4060 推定速度演算器、4070,4071 LUT、4072 減算器、4073 除算器、4074,4081Aa 比例積分器、4075,4081Ab 積分器、4080 位置推定誤差演算器、4081,4081A PLL、40803 鎖交磁束インダクタンス交流分演算器、40804 鎖交磁束インダクタンス交流分推定器、40805 回転子位置推定誤差演算器。

Claims (12)

  1.  多相の回転電機の駆動制御を行う回転電機の制御装置であって、
     回転電機に流れる回転電機電流を検出する電流検出手段と、
     前記回転電機電流に基づいて前記回転電機の回転子の回転位置である回転子位置を推定する位置推定手段と、
     前記回転電機電流と、前記回転子位置の推定値である推定回転子位置とに基づいて、前記回転電機を駆動するための回転電機駆動電圧指令と、前記回転子位置を推定するための各相の位置推定用電圧指令とを演算する制御器と、
     を備え、
     前記位置推定手段は、前記位置推定用電圧指令に基づく位置推定用電圧の印加により生じる位置推定用電流の振幅である位置推定用電流振幅の交流成分に基づく第1の推定回転子位置と、前記位置推定用電流振幅の直流成分に基づく第2の推定回転子位置と、前記回転電機電流と、前記回転子の突極性により生じる誘起電圧により演算される鎖交磁束とに基づく第3の推定回転子位置とを演算し、前記第1の推定回転子位置、前記第2の推定回転子位置及び前記第3の推定回転子位置のうちの何れか1つを選択することで得られる推定回転子位置、又は、前記第1の推定回転子位置、前記第2の推定回転子位置及び前記第3の推定回転子位置のうちの少なくとも2つを用いて合成することで得られる推定回転子位置を出力する
     ことを特徴とする回転電機の制御装置。
  2.  前記位置推定手段は、前記回転電機の速度情報が閾値より小さく、且つ、前記回転電機の磁気飽和情報が閾値より小さい場合は、前記第1の推定回転子位置を選択して出力し、前記速度情報が閾値より小さく、且つ、前記磁気飽和情報が閾値より大きい場合は、前記第2の推定回転子位置を選択して出力し、前記速度情報が閾値より大きい場合は、前記第3の推定回転子位置を選択して出力する
     ことを特徴とする請求項1に記載の回転電機の制御装置。
  3.  前記第1の推定回転子位置、前記第2の推定回転子位置及び前記第3の推定回転子位置のうちの何れか1つが選択される場合、選択されない推定回転子位置の演算処理を行う構成部は動作を停止する
     ことを特徴とする請求項2に記載の回転電機の制御装置。
  4.  前記位置推定手段は、前記回転電機の磁気飽和情報及び前記回転電機の速度情報に基づいて前記第1の推定回転子位置と前記第2の推定回転子位置と前記第3の推定回転子位置とを加重平均した推定回転子位置を出力する
     ことを特徴とする請求項1に記載の回転電機の制御装置。
  5.  前記位置推定手段は、前記回転電機の磁気飽和情報に基づいて前記第1の推定回転子位置と前記第2の推定回転子位置とを加重平均した第4の推定回転子位置を演算すると共に、前記回転電機の速度情報に基づいて、前記第4の推定回転子位置と前記第3の推定回転子位置とを加重平均した第5の推定回転子位置を演算し、前記第5の推定回転子位置を出力する
     ことを特徴とする請求項1に記載の回転電機の制御装置。
  6.  前記速度情報は、前記位置推定用電流から検出した回転子の突極性に基づく第1の推定速度と、前記位置推定用電流振幅の直流成分に基づく第2の推定速度と、前記回転電機電流及び前記鎖交磁束に基づく第3の推定速度とのうちの何れか1つである
     ことを特徴とする請求項2から5の何れか1項に記載の回転電機の制御装置。
  7.  前記第1の推定速度、前記第2の推定速度及び前記第3の推定速度のうちの何れか1つが使用される場合、使用されない推定速度の演算処理を行う構成部は動作を停止する
     ことを特徴とする請求項6に記載の回転電機の制御装置。
  8.  前記第1の推定速度、前記第2の推定速度及び前記第3の推定速度の少なくとも2つを用いた加重平均値を前記速度情報として用いる
     ことを特徴とする請求項6に記載の回転電機の制御装置。
  9.  前記位置推定手段は、前記回転電機の磁気飽和情報と前記速度情報に基づいて、前記第1の推定速度と前記第2の推定速度と前記第3の推定速度とを加重平均した推定速度を前記速度情報として用いる
     ことを特徴とする請求項6に記載の回転電機の制御装置。
  10.  前記位置推定手段は、前記回転電機の磁気飽和情報に基づいて前記第1の推定速度と前記第2の推定速度とを加重平均した第4の推定速度を演算すると共に、前記速度情報に基づいて、前記第4の推定速度と前記第3の推定速度とを加重平均した第5の推定速度を演算し、前記第5の推定速度を前記速度情報として用いる
     ことを特徴とする請求項6に記載の回転電機の制御装置。
  11.  前記磁気飽和情報は、前記回転電機に対するトルク指令値、前記トルク指令値のフィルタ出力値、前記回転電機の数学モデルにより得られる推定トルク、駆動電流を引数としたルックアップテーブルにより得られる推定トルク、前記回転電機に対する回転電機駆動電流指令、前記回転電機駆動電流指令のフィルタ出力値、及び前記駆動電流の周波数のうちの少なくとも1つである
     ことを特徴とする請求項2から10の何れか1項に記載の回転電機の制御装置。
  12.  前記磁気飽和情報は、前記位置推定用電流振幅の直流成分である
     ことを特徴とする請求項2から10の何れか1項に記載の回転電機の制御装置。
PCT/JP2019/001516 2019-01-18 2019-01-18 回転電機の制御装置 WO2020148900A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
PCT/JP2019/001516 WO2020148900A1 (ja) 2019-01-18 2019-01-18 回転電機の制御装置
DE112019006688.9T DE112019006688T5 (de) 2019-01-18 2019-01-18 Steuervorrichtung für eine rotierende elektrische maschine
JP2020566081A JP7019077B2 (ja) 2019-01-18 2019-01-18 回転電機の制御装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2019/001516 WO2020148900A1 (ja) 2019-01-18 2019-01-18 回転電機の制御装置

Publications (1)

Publication Number Publication Date
WO2020148900A1 true WO2020148900A1 (ja) 2020-07-23

Family

ID=71613552

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2019/001516 WO2020148900A1 (ja) 2019-01-18 2019-01-18 回転電機の制御装置

Country Status (3)

Country Link
JP (1) JP7019077B2 (ja)
DE (1) DE112019006688T5 (ja)
WO (1) WO2020148900A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022134772A1 (zh) * 2020-12-25 2022-06-30 中车永济电机有限公司 一种永磁辅助同步磁阻电机的控制方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014157628A1 (ja) * 2013-03-28 2014-10-02 アイシン・エィ・ダブリュ株式会社 回転電機制御装置
JP2015149875A (ja) * 2014-02-10 2015-08-20 三菱電機株式会社 電動機制御装置
WO2018179620A1 (ja) * 2017-03-27 2018-10-04 三菱電機株式会社 回転電機の制御装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5145850B2 (ja) 1972-11-08 1976-12-06
JPS51143914A (en) 1975-06-06 1976-12-10 Tokyo Gas Co Ltd Liquid spray proportion controlling means

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014157628A1 (ja) * 2013-03-28 2014-10-02 アイシン・エィ・ダブリュ株式会社 回転電機制御装置
JP2015149875A (ja) * 2014-02-10 2015-08-20 三菱電機株式会社 電動機制御装置
WO2018179620A1 (ja) * 2017-03-27 2018-10-04 三菱電機株式会社 回転電機の制御装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022134772A1 (zh) * 2020-12-25 2022-06-30 中车永济电机有限公司 一种永磁辅助同步磁阻电机的控制方法

Also Published As

Publication number Publication date
JP7019077B2 (ja) 2022-02-14
JPWO2020148900A1 (ja) 2021-09-09
DE112019006688T5 (de) 2021-10-07

Similar Documents

Publication Publication Date Title
JP5324646B2 (ja) 回転電機の制御装置および制御方法
Sayeef et al. Rotor position and speed estimation of a variable structure direct-torque-controlled IPM synchronous motor drive at very low speeds including standstill
JP4665360B2 (ja) 電動機制御装置
JP4425193B2 (ja) モータの位置センサレス制御装置
US7679308B2 (en) Motor control device
JP4674516B2 (ja) 同期モータの磁極位置推定方法
JP5324159B2 (ja) モータ制御装置
JP6324627B2 (ja) 交流回転機の制御装置および電動パワーステアリングの制御装置
JP5774136B2 (ja) 交流回転機の制御装置
US11515821B2 (en) Control device for rotating electrical machine
TWI565219B (zh) 交流旋轉電機之控制裝置
JP4899509B2 (ja) 交流電動機の回転子位相推定装置
WO2020148900A1 (ja) 回転電機の制御装置
JP2010166638A (ja) 回転電機の制御装置
Zhao et al. An extended flux model-based rotor position estimator for sensorless control of interior permanent magnet synchronous machines
Mubarok et al. An Adjustable Wide-Range Speed-Control Method for Sensorless IPMSM Drive Systems
JP7166443B2 (ja) 交流回転機の制御装置及び交流回転機の制御方法
Liu et al. Implementation of a novel high-performance sensorless IPMSM drive for washing machines
JP7199605B1 (ja) 回転機制御装置
JP2010022188A (ja) モータの位置センサレス制御装置
JP6520686B2 (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: 19909799

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2020566081

Country of ref document: JP

Kind code of ref document: A

122 Ep: pct application non-entry in european phase

Ref document number: 19909799

Country of ref document: EP

Kind code of ref document: A1